本發(fā)明涉及油氣地球物理技術(shù)領(lǐng)域，更具體地，涉及河道儲層的地震屬性的預(yù)測方法。

背景技術(shù)：

河道砂體是碎屑巖中最重要的一類油氣儲集體。由于河道頻繁改道與遷移,使得多期河道沉積地層中廣泛發(fā)育砂泥混雜相和砂泥交互相。這類砂體不僅相互交叉、疊置或切割，空間分布穩(wěn)定性差，而且儲層厚度薄，因此河道儲層平面分布預(yù)測難度很大。

目前，預(yù)測河道砂巖儲層的地震屬性方法大致可分為兩類，一類主要是利用各類振幅類屬性方法，如利用均方根振幅，最大波峰(或最大波谷)振幅，振幅變化率，以及分頻振幅屬性等來預(yù)測河道的平面展布；另一類是利用邊緣檢測手段刻畫河道的邊界，如利用小波變換的多尺度邊緣檢測、最優(yōu)濾波二階差分邊緣檢測和地震相干分析邊緣檢測等。

前一類預(yù)測方法利用了河道儲層地震響應(yīng)的物理屬性特征，即河道儲層與圍巖之間波阻抗的差異；后一類方法利用了河道砂體展布形態(tài)及其邊界的幾何類波形特征。

顯然，這兩類地震屬性預(yù)測方法都沒有充分利用河道砂體地震響應(yīng)的綜合信息，而只是利用了儲層地震響應(yīng)中某一方面的特征，均存在一定的局限性。因此，這兩類河道儲層地震屬性預(yù)測方法應(yīng)用范圍有限，精度較低。

技術(shù)實(shí)現(xiàn)要素：

為了提高地震屬性預(yù)測河道儲層的精度和適應(yīng)性，本公開提出了一種河道 儲層的地震屬性的預(yù)測方法，該方法以所述河道儲層的屬性參數(shù)數(shù)據(jù)體A來預(yù)測所述河道儲層的地震屬性，其包括如下步驟：S10.估算地震資料的波形相似性數(shù)據(jù)體S；S20.估算并篩選出反映所述河道儲層的敏感地震參數(shù)屬性數(shù)據(jù)體P；S30.根據(jù)公式A＝P/(K·S)計(jì)算所述屬性參數(shù)數(shù)據(jù)體A，其中，K為預(yù)定的加權(quán)系數(shù)。

優(yōu)選地，所述步驟S10中，所述波形相似性數(shù)據(jù)體S的估算方法為：在目標(biāo)層段內(nèi)，垂直于主河道的延伸方向，進(jìn)行地震反射波組的波形相似性分析，從而估算所述波形相似性數(shù)據(jù)體S。其中，所述目標(biāo)層段的確定方法為：根據(jù)鉆井資料進(jìn)行井震標(biāo)定，確定河道沉積在地震剖面上的時間段；在區(qū)域上確定在地震工區(qū)內(nèi)的河道沉積的時間段范圍[Tmin,Tmax]，從而得到所述目標(biāo)層段。所述主河道的延伸方向的確定方法為：根據(jù)沉積相分析，得到所述主河道的延伸方向。

優(yōu)選地，所述步驟S20中，所述敏感地震參數(shù)屬性數(shù)據(jù)體P的估算方法為：根據(jù)鉆井資料的河道厚度、物性或含油氣性與地震振幅、頻率的屬性參數(shù)的相關(guān)性分析，篩選出反映所述河道儲層的敏感地震參數(shù)屬性數(shù)據(jù)體P。

進(jìn)一步地，所述步驟S20中，在估算所述敏感地震參數(shù)屬性數(shù)據(jù)體P之前，還包括：估算所述地震資料的均方根振幅數(shù)據(jù)體P₁、反射能量數(shù)據(jù)體P₂、振幅變化率數(shù)據(jù)體P₃、地震平均振幅數(shù)據(jù)體P₄、以及頻率數(shù)據(jù)體P₅，從P₁、P₂、P₃、P₄、和P₅中篩選出反映所述河道儲層的敏感地震參數(shù)屬性數(shù)據(jù)體P。

進(jìn)一步地，在所述步驟S10之后、所述步驟S30之前，還包括步驟S15：從估算的所述波形相似性數(shù)據(jù)體S中，在所述河道沉積的時間段范圍[Tmin,Tmax]內(nèi)，確定所述波形相似性數(shù)據(jù)體S的對應(yīng)的數(shù)值范圍[Smin,Smax]，以確定所述加權(quán)系數(shù)K。

進(jìn)一步地，在所述步驟S15之后、所述步驟S30之前，還包括步驟S25：在所述波形相似性數(shù)據(jù)體S上，用所述加權(quán)系數(shù)K，對所述地震相似性數(shù)據(jù)體S進(jìn)行加權(quán)增強(qiáng)。

進(jìn)一步地，在所述步驟S30之后，還包括步驟S40：從所述屬性參數(shù)數(shù)據(jù)體A中提取在所述河道沉積的時間段范圍[Tmin,Tmax]內(nèi)的數(shù)據(jù)，并進(jìn)行顯示。

本發(fā)明的河道儲層的地震屬性的預(yù)測方法，基于地震物理響應(yīng)特征和幾何響應(yīng)特征的混合地震屬性，有效識別預(yù)測出優(yōu)質(zhì)的河道砂巖儲層，極大地降低河道砂巖油氣勘探開發(fā)鉆井的風(fēng)險(xiǎn)性，提升地震采集資料的價值。

附圖說明

通過結(jié)合附圖對本公開示例性實(shí)施方式進(jìn)行更詳細(xì)的描述，本公開的上述以及其它目的、特征和優(yōu)勢將變得更加明顯，其中，在本公開示例性實(shí)施方式中，相同的參考標(biāo)號通常代表相同部件。

圖1示出了根據(jù)本發(fā)明的一個實(shí)施例的河道儲層的地震屬性的預(yù)測方法的流程圖。

圖2示出了根據(jù)本發(fā)明的具體應(yīng)用示例的河道儲層的地震屬性的預(yù)測方法的流程圖。

圖3示出了某工區(qū)地震波形相似性屬性表征的河道儲層分布圖。

圖4示出了某工區(qū)地震敏感屬性表征的河道儲層分布圖。

圖5示出了某工區(qū)利用本發(fā)明的河道儲層的地震屬性的預(yù)測方法表征的河道儲層分布圖。

具體實(shí)施方式

下面將參照附圖更詳細(xì)地描述本公開的優(yōu)選實(shí)施方式。雖然附圖中顯示了本公開的優(yōu)選實(shí)施方式，然而應(yīng)該理解，可以以各種形式實(shí)現(xiàn)本公開而不應(yīng)被這里闡述的實(shí)施方式所限制。相反，提供這些實(shí)施方式是為了使本公開更加透徹和完整，并且能夠?qū)⒈竟_的范圍完整地傳達(dá)給本領(lǐng)域的技術(shù)人員。

通常，河道類砂巖儲層在垂直河道方向具有沉積巖性橫向變化快、連續(xù)性差、非均質(zhì)性強(qiáng)等特點(diǎn)；而在順延河道方向上卻具有砂體連續(xù)性好、巖相相對 穩(wěn)定等特點(diǎn)。此外，隨著河道砂巖的巖性、物性和含油氣性的差異，會與河道間或河道砂巖上下的泥巖沉積形成一定的波阻抗差異，而且不同部位河道沉積砂巖的厚度也有很大的變化。

圖1示出了根據(jù)本發(fā)明的一個實(shí)施例的河道儲層的地震屬性的預(yù)測方法的流程圖。該方法以河道儲層的屬性參數(shù)數(shù)據(jù)體A來預(yù)測河道儲層的地震屬性，其包括如下步驟：

S10.估算地震資料的波形相似性數(shù)據(jù)體S；

S20.估算并篩選出反映河道儲層的敏感地震參數(shù)屬性數(shù)據(jù)體P；

S30.根據(jù)公式A＝P/(K·S)計(jì)算屬性參數(shù)數(shù)據(jù)體A，其中，K為預(yù)定的加權(quán)系數(shù)。

本實(shí)施例利用河道砂巖沉積的前述特點(diǎn)，先估算出地震資料的波形相似性數(shù)據(jù)體S，提取出河道砂巖大致的沉積分布；再采用估算加權(quán)系數(shù)K，突出河道響應(yīng)；然后估算出多種地震屬性數(shù)據(jù)體(P₁，P₂，P₃，P₄，P₅)，并根據(jù)鉆井資料從(P₁，P_2、，P₃，P₄，P₅)中篩選出反映河道儲層物性或含油氣性的敏感地震參數(shù)屬性數(shù)據(jù)體P；最后把這兩種屬性復(fù)合在一起，計(jì)算出河道儲層屬性參數(shù)數(shù)據(jù)體A，達(dá)到實(shí)現(xiàn)了沉積相控下河道儲層預(yù)測的目的。

為便于理解本發(fā)明實(shí)施例的方案及其效果，以下給出一個具體應(yīng)用示例。本領(lǐng)域技術(shù)人員應(yīng)理解，該示例僅為了便于理解本發(fā)明，其任何具體細(xì)節(jié)并非意在以任何方式限制本發(fā)明。

圖2示出了根據(jù)本發(fā)明的具體應(yīng)用示例的河道儲層的地震屬性的預(yù)測方法的流程圖。在該實(shí)施例中，該方法包括如下步驟：

S10.根據(jù)鉆井資料(聲波、密度測井等數(shù)據(jù))和地震資料(主要是地震疊偏數(shù)據(jù))，確定目標(biāo)層時間段。然后在目標(biāo)層段的時間窗口范圍[Tmin,Tmax]內(nèi)，即河道儲層在地震資料中出現(xiàn)的最小時間Tmin和最大時間Tmax范圍內(nèi)，垂直于主河道的延伸方向，進(jìn)行地震反射波組的波形相似性分析，估算出地震資料的波形相似性數(shù)據(jù)體S。

河道儲層目標(biāo)層段的確定方法為：根據(jù)鉆井資料進(jìn)行井震標(biāo)定，確定河道沉積在地震剖面上的時間段；在區(qū)域上確定河道沉積在地震工區(qū)內(nèi)的時間段，即目標(biāo)層段的時間范圍。主河道的延伸方向的確定方法為：根據(jù)沉積相分析，得到所述主河道的延伸方向。

S15和S25.從地震資料中估算河道沉積的地震反射波形的相似性數(shù)數(shù)據(jù)體S，確定河道沉積對應(yīng)的時間段范圍[Tmin,Tmax]內(nèi)相似性數(shù)據(jù)體S的數(shù)值范圍[Smin,Smax]，即S數(shù)據(jù)體的最大值Smax和最小值Smin；在河道沉積對應(yīng)的時間段范圍[Tmin,Tmax]內(nèi)，根據(jù)[Smin,Smax]對波形相似性數(shù)據(jù)體S進(jìn)行加權(quán)增強(qiáng)，以突出河道的沉積響應(yīng)。其加權(quán)系數(shù)為預(yù)定的加權(quán)系數(shù)K，即得到K·S。

S20.估算地震資料的均方根振幅數(shù)據(jù)體P₁、反射能量數(shù)據(jù)體P₂，振幅變化率數(shù)據(jù)體P₃、地震平均振幅數(shù)據(jù)體P₄、以及頻率數(shù)據(jù)體P₅；根據(jù)河道厚度、物性或含油氣性與地震振幅、頻率屬性的相關(guān)性分析，從P₁，P₂，P₃，P₄，和P₅多種地震屬性數(shù)據(jù)體中篩選出反映所述河道儲層的敏感地震參數(shù)屬性數(shù)據(jù)體P；

S30.根據(jù)公式A＝P/(K·S)來計(jì)算所述屬性參數(shù)數(shù)據(jù)體A，其中，K為前述步驟S15和S25中預(yù)定的加權(quán)系數(shù)。

S40.從所述屬性參數(shù)數(shù)據(jù)體A中提取所述河道沉積對應(yīng)的時間段范圍[Tmin,Tmax]內(nèi)的數(shù)據(jù)，并進(jìn)行顯示，以便進(jìn)行分析。

為了便于理解本發(fā)明實(shí)施例的技術(shù)效果，圖3示出了某工區(qū)地震波形相似性屬性表征的河道儲層分布圖；圖4示出了某工區(qū)地震敏感屬性表征的河道儲層分布圖；圖5示出了某工區(qū)利用本發(fā)明的河道儲層的地震屬性的預(yù)測方法表征的河道儲層分布圖。對比這三幅圖可知，應(yīng)用本發(fā)明的河道儲層的地震屬性的預(yù)測方法得到的屬性參數(shù)數(shù)據(jù)體A的儲層分布顯示(圖5所示)，比圖3和圖4中示出的現(xiàn)有技術(shù)中的用單一屬性的參數(shù)來表征河道儲層的方案，其分布圖中的邊界更為清晰，整個預(yù)測圖像的河道儲層識別的精度更高。

本發(fā)明的河道儲層的地震屬性的預(yù)測方法，融合了有關(guān)河道形態(tài)的地震反射幾何學(xué)方面的特征，又突出了與儲層巖性、物性和含油氣性變化相關(guān)的物理 特性，真正意義上實(shí)現(xiàn)了沉積相控下的儲層預(yù)測。應(yīng)用該方法得到的新的復(fù)合屬性參數(shù)與單一屬性參數(shù)相比，表征河道儲層細(xì)節(jié)分布的能力明顯增強(qiáng)，刻畫出的河道分布圖的精度大幅提高，降低了后續(xù)油氣勘探開發(fā)井位鉆探的風(fēng)險(xiǎn)。

附圖中的流程圖顯示了根據(jù)本公開的多個實(shí)施例的方法的可能實(shí)現(xiàn)的功能和操作。在這點(diǎn)上，流程圖中的每個方框可以代表一個模塊、程序段或指令的一部分，所述模塊、程序段或指令的一部分包含一個或多個用于實(shí)現(xiàn)規(guī)定的邏輯功能的可執(zhí)行指令。在有些作為替換的實(shí)現(xiàn)中，方框中所標(biāo)注的功能也可以以不同于附圖中所標(biāo)注的順序發(fā)生。例如，兩個連續(xù)的方框?qū)嶋H上可以基本并行地執(zhí)行，它們有時也可以按相反的順序執(zhí)行，這依所涉及的功能而定。也要注意的是，流程圖中的每個方框、以及流程圖中的方框的組合，可以用執(zhí)行規(guī)定的功能或動作的專用的基于硬件的系統(tǒng)來實(shí)現(xiàn)，或者可以用專用硬件與計(jì)算機(jī)指令的組合來實(shí)現(xiàn)。

以上已經(jīng)描述了本公開的各實(shí)施例，上述說明是示例性的，并非窮盡性的，并且也不限于所披露的各實(shí)施例。在不偏離所說明的各實(shí)施例的范圍和精神的情況下，對于本技術(shù)領(lǐng)域的普通技術(shù)人員來說許多修改和變更都是顯而易見的。本文中所用術(shù)語的選擇，旨在最好地解釋各實(shí)施例的原理、實(shí)際應(yīng)用或?qū)κ袌鲋械募夹g(shù)的技術(shù)改進(jìn)，或者使本技術(shù)領(lǐng)域的其它普通技術(shù)人員能理解本文披露的各實(shí)施例。