地震反演在遼河灘海西部儲層預(yù)測中的應(yīng)用

董德勝,郭彥民, 鄒丙方,楊光達,孫 轉(zhuǎn)

(中國石油遼河油田分公司勘探開發(fā)研究院,遼寧盤錦 124010)

?

地震反演在遼河灘海西部儲層預(yù)測中的應(yīng)用

董德勝,郭彥民, 鄒丙方,楊光達,孫 轉(zhuǎn)

(中國石油遼河油田分公司勘探開發(fā)研究院,遼寧盤錦 124010)

遼河灘海西部的勘探目標日趨復(fù)雜，儲層預(yù)測工作越來越重要。闡述了在遼河灘海西部地區(qū)開展地震反演工作的關(guān)鍵技術(shù)環(huán)節(jié)，并對反演效果進行評價分析；利用波阻抗反演結(jié)果定量預(yù)測主要目的層段的砂巖厚度和砂巖含量，并對各層段儲層的發(fā)育程度進行評價，分析儲層的平面分布特征，確定巖性圈閉發(fā)育有利區(qū)，為巖性圈閉識別及評價提供依據(jù)。

遼河灘海；地震反演；λ值；砂巖含量預(yù)測；儲層發(fā)育程度

地震反演是儲層預(yù)測的核心技術(shù)之一，是地質(zhì)綜合評價與地震資料解釋的重要組成部分，在圈閉評價、油藏描述以及井位論證中都扮演著重要角色[1]。在遼河灘海西部地區(qū)，隨著勘探程度的逐步提高，勘探目標日趨復(fù)雜，從以往的構(gòu)造油藏向地層巖性油藏轉(zhuǎn)變，儲層預(yù)測工作也越來越受到重視。研究區(qū)儲集層以薄互層沉積砂體為主，具有單層厚度小、縱橫向變化快的特點，地震預(yù)測多解性強，難度大。本次研究綜合利用地質(zhì)、地震、鉆井等資料，在精細層序地層格架約束下，采用地震波阻抗反演技術(shù)預(yù)測了目的層段(沙一、二段)儲層的縱、橫向發(fā)育特征，取得明顯效果，為地層巖性油藏勘探提供了有力支撐。

1 區(qū)域地質(zhì)背景

遼河灘海地區(qū)地理上位于遼東灣北部, 是遼河盆地陸上向海域的自然延伸，構(gòu)造上位于遼河灘海西部凹陷，包括西部斜坡帶、筆架嶺-嶺南構(gòu)造帶、海南洼陷帶、仙鶴-架東構(gòu)造帶4個構(gòu)造單元，總面積約780 km2。遼河灘海西部凹陷呈箕狀特點，西側(cè)為緩坡帶，東側(cè)為陡坡帶，中間為深陷帶[2]。在東西兩側(cè)高中間低的古地貌背景下，其周邊存在西部凸起、葫蘆島凸起和海月凸起3個主要的物源區(qū)(圖1)，沙一、二期發(fā)育扇三角洲-湖泊相沉積[3-5]。

圖1 遼河灘海西部沙一段頂界構(gòu)造圖

本次研究將遼河灘海西部古近系劃分為4個三級層序，15個準層序組(表1)，其中主要目的層段SQ3(沙一、二段)為一個完整的沉積旋回。低位域為砂礫巖、粗砂巖夾紅色泥巖、灰色泥巖構(gòu)成的正旋回發(fā)育扇三角洲沉積體系，湖盆水體淺，物源供給充足，扇體發(fā)育規(guī)模大；水進域為細砂巖、粉砂巖夾灰色泥巖、油頁巖構(gòu)成的正旋回，仍發(fā)育扇三角洲沉積體系，湖盆水體變深，扇體向物源區(qū)收縮，分布范圍變小。高位域為深灰色泥巖、細砂巖、粉砂巖構(gòu)成的反旋回，也屬于扇三角洲沉積體系，水體再變淺，扇體規(guī)模增大(圖2)。

本次研究細化研究單元，為儲層預(yù)測提供精細的等時層序地層格架。將SQ3細分為5個準層序組，低位域劃分為2個(沙二上、沙二下)，水進域劃分2個(沙一下、沙一中)，高位域劃分1個(沙一上)。鉆探證實，研究區(qū)低位域晚期(沙二上)和水進域早期(沙一下)扇三角洲沉積砂體發(fā)育，儲蓋配置關(guān)系好，為主力出油層段，是本次研究重點目的層段。

圖2 F井-Q井地層對比圖

2 地震品質(zhì)及鉆井巖電特征分析

地震反演效果的好壞主要取決于兩個條件，一是地震資料的品質(zhì)，包括信噪比、分辨率及保真度，這一點決定了反演結(jié)果的分辨率和可信度；二是鉆井資料的巖電組合特征，包括聲波時差、密度及縱波阻抗曲線能否很好地反映儲層和圍巖的差異，這決定了波阻抗反演結(jié)果能否有效地區(qū)分儲層與非儲層。

2.1 地震資料品質(zhì)分析

主要以高分辨處理地震資料為基礎(chǔ)，其資料品質(zhì)較以往有了較大改善。在目的層段地震資料主頻達到25 Hz,有效頻帶寬度10～50 Hz。理論上,地震垂向分辨率為調(diào)諧厚度,即H=λ/4,其中λ=V/F[6]。在研究區(qū)沙一、二段儲層層速度約為3 000 m/s,因此,地震資料垂向分辨率大約為30 m，可以識別大于30 m的砂層組。

2.2 目的層巖電組合特征分析

通過對工區(qū)30余口探井巖性資料及電性資料的綜合分析，目的層段巖性組合為砂、泥互層沉積，儲集層以砂巖、礫巖為主，儲集層與上下圍巖的縱波阻抗具有明顯差異，礫巖波阻抗最高，砂巖次之，泥巖最小(圖3)，因此，可以利用波阻抗進行巖性區(qū)分。

3 地震波阻抗反演

3.1 地震反演方法及原理

波阻抗反演處理采用Jason反演軟件中的約束稀疏脈沖反演模塊進行。該方法的主要優(yōu)點是能獲得寬頻帶的反射系數(shù),能較好地解決地震記錄的欠定問題,從而使波阻抗反演得到的數(shù)據(jù)更趨于真實[5]。稀疏脈沖波阻抗反演的主要過程是:首先，通過最大似然反褶積求得一個具有稀疏特性的反射系數(shù)序列；然后，通過最大似然反演導(dǎo)出寬帶波阻抗[7]。

圖3 泥巖、砂巖、礫巖波阻抗直方圖

3.2 子波提取及井震標定分析

地震子波提取及井震標定分析是地震反演工作的基礎(chǔ)[8]。為實現(xiàn)子波在能量和頻率上與實際地震道吻合，研究過程中地震子波提取由實際地震資料出發(fā)，用多道地震記錄進行自相關(guān)統(tǒng)計的方法求取子波[9]。結(jié)合井震標定分析，最終確定研究區(qū)子波為一個零相位的正極性子波，子波主頻25 Hz左右。從井震標定結(jié)果分析，合成地震記錄與原始地震道吻合好，波組對應(yīng)關(guān)系清晰，波峰對應(yīng)砂組的頂界面，波谷對應(yīng)砂組底界面(圖4)。

圖4 C井震綜合標定

3.3 三維初始波阻抗模型建立

測井約束反演需要建立正確的地質(zhì)先驗?zāi)Ｐ停慈S初始波阻抗模型。模型建立的基礎(chǔ)資料是構(gòu)造解釋層位。為了提高模型的精度，筆者把目的層沙一、二段細分5個研究單元，將各準層序組界面引入模型，并對層間產(chǎn)狀及地層接觸關(guān)系進行了準確定義，搭建了一個高精度框架模型。在此基礎(chǔ)上，再利用測井波阻抗曲線進行內(nèi)插外推，產(chǎn)生三維波阻抗模型體(圖5)。

3.4 反演參數(shù)選取

圖5 F井-P井初始阻抗模型剖面

地震波阻抗反演參數(shù)的選擇直接影響反演效果，根據(jù)研究區(qū)地質(zhì)情況及資料特點，通過反復(fù)測試確定了相關(guān)參數(shù)[10]。根據(jù)試驗，地震子波選取由實際資料統(tǒng)計的長度為80 ms的零相位子波；反演處理時窗為沙一段頂界-沙三上段底界；波阻抗橫向變化可拓展范圍為20%；反演迭代次數(shù)為6次；低頻合并分量確定為10 Hz, 殘差權(quán)重因子λ值取18。

3.5 反演效果分析

在對子波提取、層位標定和波阻抗初始模型建立等各項反演關(guān)鍵環(huán)節(jié)進行精細處理之后，選擇合理的控制參數(shù)開展波阻抗反演，并對反演效果進行質(zhì)控分析，以驗證預(yù)測結(jié)果的可靠性。首先，將反演得到的合成地震記錄與原始地震道進行對比發(fā)現(xiàn)，除了在局部斷層發(fā)育的部位存在殘差外，整體對應(yīng)關(guān)系較好，表明反演波阻抗結(jié)果忠實于地震資料，能夠客觀地反映地下巖石波阻抗變化特征。然后，通過井旁反演的波阻抗曲線與原始測井阻抗曲線對比分析發(fā)現(xiàn)，地震有效頻帶內(nèi)兩者的曲線特征基本一致，這表明初始波阻抗模型約束的比較合理，反演結(jié)果客觀地反映了儲層的變化特征。最后，對反演結(jié)果標定分析，波阻抗反演結(jié)果與已知井有較好的對比關(guān)系。

4 儲層定量預(yù)測

為了進一步研究儲層在空間上的分布規(guī)律，精細描述儲層在空間上的分布狀況，利用波阻抗反演結(jié)果，確定砂巖和泥巖阻抗門檻值，分別計算了沙一上、中、下段，沙二上、下段的砂巖厚度和砂巖含量。經(jīng)鉆探證實，研究區(qū)沙一下段是儲層發(fā)育段，也是油氣富集段，因此下面以沙一下段為例對儲層定量預(yù)測效果進行分析。

4.1 砂巖含量及厚度預(yù)測

沙一下段沉積期為SQ3湖侵期早期，研究區(qū)水體加深，可容納空間增大，地層沉積范圍增大，該時期主要發(fā)育扇三角洲沉積體系。從沙一下段砂巖含量預(yù)測結(jié)果分析(圖6)，儲層較為發(fā)育。在研究區(qū)北部，高值區(qū)集中在F井、G井、I井、L井至N井區(qū)，砂巖含量在30%以上，由西向東砂巖含量逐漸降低，在M井、N井區(qū)形成一相對獨立的砂體。結(jié)合地質(zhì)背景分析，推測該套砂體物源主要來自西部凸起，在仙鶴構(gòu)造下傾部位傾沒。該套砂體多井鉆遇，如I井、J井、L井等，發(fā)現(xiàn)了多個巖性油氣藏。在研究區(qū)中部，砂巖含量高值區(qū)集中在B井區(qū)至C井區(qū)附近，結(jié)合古地貌分析，物源主要來自葫蘆島凸起，扇體向北東向延伸。在南部嶺南斜坡區(qū)，砂巖含量高值區(qū)具朵葉狀分布的特征，A井的東側(cè)砂巖含量在30%以上，儲層發(fā)育。研究區(qū)東部仙鶴-架東地區(qū)，砂巖含量普遍偏低，高值區(qū)零星分布，規(guī)模小，推測儲層不發(fā)育。

圖6 沙一下段(SQ33)砂巖含量預(yù)測

根據(jù)沙一下段砂巖厚度預(yù)測結(jié)果(圖7)，研究區(qū)北部砂體較厚的區(qū)域主要分布于L井、M井、N井區(qū)，最大厚度達110 m，具有向東及南北減薄的特征，易形成巖性上傾尖滅等圈閉類型。在筆架嶺斜坡帶東側(cè)，砂體呈朵葉狀分布，最大厚度達70 m，砂體厚度向東側(cè)和南側(cè)減薄，與斜坡區(qū)北東向斷裂相配置，利于形成斷層-巖性圈閉等圈閉類型。

圖7 沙一下段(SQ33)砂巖厚度預(yù)測

4.2 儲層定量預(yù)測效果分析

為評價儲層預(yù)測結(jié)果的可靠性，從已知井出發(fā)，對各層段儲層預(yù)測厚度與實際鉆井揭示厚度進行對比分析。以沙一下段為例，儲層預(yù)測吻合率總體平均值達到87%，取得較好的預(yù)測效果，反映了研究區(qū)儲層的平面展布特征，為儲層綜合評價及巖性圈閉的識別提供了有力支撐。

4.3 儲層發(fā)育程度預(yù)測

如前文所述，通過對各準層序組砂巖含量及砂巖厚度的定量預(yù)測，確定了砂體的橫向展布特征。但是，受地層厚度的影響，有些區(qū)域砂巖含量高，但砂體厚度相對較??；有些區(qū)域砂體厚度比較大，但砂巖含量比較低；因此，上述這些區(qū)域可能都不是尋找?guī)r性圈閉的有利區(qū)域。為此，進一步將各層段砂巖含量與砂巖厚度進行交會分析，對各層段儲層的發(fā)育程度進行綜合評價。根據(jù)層段砂巖厚度與砂巖含量的數(shù)據(jù)分布確定門檻值，將交會結(jié)果劃分成5個區(qū)域(圖8)。I類區(qū)是紅色區(qū)域，對應(yīng)砂巖含量較高(大于50%)，砂巖累計厚度大的區(qū)域(大于30 m)；II類區(qū)砂巖含量中等(大于30%，小于50%)，是砂巖厚度大的區(qū)域；III類區(qū)砂巖含量高(大于50%)，是厚度小的區(qū)域(小于30 m)；IV類區(qū)砂巖含量中等(大于30%，小于50%)，是厚度小的區(qū)域(小于30 m)；V類區(qū)是砂巖含量較小的區(qū)域(小于30%)。

從沙一下段儲層發(fā)育程度預(yù)測圖分析(圖9)，主要以I類、II類和V類區(qū)為主。從已知井標定來看，預(yù)測結(jié)果與井對應(yīng)關(guān)系較好，反映了儲層的橫向變化情況。I類區(qū)主要集中西部斜坡帶D井區(qū)、E井區(qū)，在A井區(qū)、H井區(qū)、J井區(qū)及M井區(qū)局部分布；II類區(qū)分布范圍較廣，I井區(qū)、N井區(qū)及B井區(qū)均處于該區(qū)域。I類區(qū)面積約100 km2，II類區(qū)面積約223 km2。鉆探證實，該層段I類和II類區(qū)均是尋找?guī)r性、構(gòu)造-巖性圈閉的有利區(qū)域，具有較大勘探潛力。

圖8 沙一下段(SQ33)砂巖含量與砂巖厚度交會圖

圖9 沙一下段(水進域準層序組I)儲層發(fā)育程度預(yù)測

5 結(jié)論

(1)本次研究綜合利用了地質(zhì)、地震、鉆井等數(shù)據(jù)，根據(jù)研究區(qū)資料及儲層特點，采用Jason反演軟件中的約束稀疏脈沖反演模塊進行了層序約束下的波阻抗反演，取得明顯效果。

(2)在準層序組格架約束下建立波阻抗初始模型，小層約束模型較以往的大層約束模型精度有較大的改善，提高了預(yù)測結(jié)果的可靠性。

(3)基于波阻抗反演結(jié)果，通過分析井點處波阻抗與巖性的對應(yīng)關(guān)系，建立巖性量版，進行巖性解釋，定量預(yù)測了目的層的砂巖厚度和砂巖含量，落實了儲層的發(fā)育特征。

(4)通過砂巖含量與砂巖厚度交會分析，對目的層段儲層的發(fā)育程度進行評價，定量劃分出5個區(qū)域，結(jié)合鉆探情況預(yù)測I類和II類區(qū)是巖性圈閉及構(gòu)造-巖性圈閉發(fā)育有利區(qū)。

[1] 魏冬，王宏語，杜羽書,等.地震反演技術(shù)在儲層預(yù)測中的應(yīng)用——以遼河灘海東部凹陷東營組為例[J].油氣藏評價與開發(fā)，2013，3(1)：9-10.

[2] 周宗明,于炳松,張聰,等.遼河灘海地區(qū)西部凹陷古近系巖性和地層圈閉形成條件與發(fā)育模式[J].油氣地質(zhì)與采收率,2011,18(1):11-14.

[3] 劉媛媛,于炳松,朱金富,等.遼河灘海西部凹陷古近系碎屑巖儲層成巖作用及其對儲層物性的影響[J].現(xiàn)代地質(zhì),2009,23(4):731-738.

[4] 田克勤,于志海,馮明,等.渤海灣盆地下第三系深層油氣地質(zhì)與勘探[M].北京:石油工業(yè)出版社,2000：89-120.

[5] 孫洪斌,謝文彥.遼河灘海油氣勘探與開發(fā)[M].北京:地質(zhì)出版社,2008：89-100.

[6] 姜書林,梁浩,毛世權(quán),等.地震反演技術(shù)在三塘湖盆地油氣勘探中的應(yīng)用[J].吐哈油氣，2007，12(4)：308-309.

[7] 穆星.稀疏脈沖波阻抗反演參數(shù)對反演效果的影響研究[J].工程地球物理學報，2005，2(2)：104-108.

[8] 安鴻偉,李正文,李仁甫, 等.稀疏脈沖波阻抗反演在XY 油田開發(fā)中的應(yīng)用[J].石油物探,2002, 41 (1): 56-60.

[9] 梁光河.地震子波提取方法研究[J].石油物探,1998,37(1):31-39.

[10] 趙華，江潔，鄭建軍.Jason地震反演技術(shù)應(yīng)用中的幾個重要環(huán)節(jié)[J].油氣地球物理,2003，1(1):41-42.

編輯：趙川喜

1673-8217(2016)06-0040-06

2016-06-02

董德勝，工程師，1983年生，2006年畢業(yè)于東北石油大學電子信息工程專業(yè),目前從事地震解釋及勘探部署工作。

P631.445.9

A