地震屬性融合技術在海拉爾盆地烏東斜坡帶南屯組儲層預測中的應用

朱可丹， 王雅春， 衣啟樊， 王昕堯， 李夢瑤，趙振鐸， 彭瀧毅

(東北石油大學 a.地球科學學院，b.“非常規(guī)油氣成藏與開發(fā)”省部共建國家重點實驗室培育基地 ， 大慶 163318)

地震屬性融合技術在海拉爾盆地烏東斜坡帶南屯組儲層預測中的應用

朱可丹a,b， 王雅春a,b， 衣啟樊a(chǎn),b， 王昕堯a,b， 李夢瑤a,b，趙振鐸a,b， 彭瀧毅a,b

(東北石油大學 a.地球科學學院，b.“非常規(guī)油氣成藏與開發(fā)”省部共建國家重點實驗室培育基地 ， 大慶 163318)

針對海拉爾盆地烏南次凹烏東斜坡帶地區(qū)南屯組砂巖儲層的特征和勘探難點，應用三維地震資料多屬性多元線性回歸融合技術，對南屯組砂巖儲層敏感的地震屬性提取并優(yōu)選，獲得地震屬性與鉆遇砂巖厚度的相關關系，建立砂巖厚度與多地震屬性之間的回歸方程，實現(xiàn)對砂巖儲層厚度高精度地預測，進而明確烏東斜坡帶地區(qū)南二段Ⅰ油組砂巖儲層厚度發(fā)育特征。該方法應用效果良好，為烏東斜坡帶地區(qū)南二段Ⅰ油組隱蔽圈閉識別預測提供借鑒。

烏東斜坡帶； 南屯組； 地震屬性融合； 多元線性回歸； 儲層砂巖預測

0 引言

烏東斜坡帶位于海拉爾盆地烏爾遜凹陷南部地區(qū)，作為一個陸相斷陷湖盆的構造坡折帶，烏東斜坡帶儲層物性好，埋藏淺，距油源較近，部分區(qū)域處于成熟烴源巖范圍，利于油氣的聚集、成藏[1]，是海拉爾盆地的主要勘探地區(qū)之一。下白堊統(tǒng)南屯組的含油層系埋深普遍大于1 500 m，目前鉆遇并獲得工業(yè)油氣流的產(chǎn)層均為有效儲層，是烏東斜坡帶的一個主力油層。對烏東斜坡帶南屯組油藏的研究，目前已經(jīng)從常規(guī)生、儲、蓋要素研究進入到隱蔽圈閉識別預測階段。應用地震屬性和儲層預測方法，可以精細刻畫隱蔽圈閉的平面和空間展布[2]。故儲層的精準預測是隱蔽圈閉識別的前提，探索提高儲層定量預測精度的方法有著具有理論和現(xiàn)實意義。

徐占東[4]針對烏南次凹南屯組儲層，利用均方根振幅和平均反射強度這兩種跟井匹配程度較高的地震屬性開展了定性預測；謝麗華[5]則利用地震相技術分析烏東斜坡帶有利儲集相帶發(fā)育模式，并認為均方根振幅屬性異常能夠較好地反映該區(qū)沉積相內(nèi)部的結構變化。儲層預測的根本目標是確定砂巖分布規(guī)律，而砂巖在不同的沉積體系(如扇三角洲、辮狀河、湖泊等)中，地震屬性響值不盡相同，但是響應規(guī)律和指示特征是統(tǒng)一的，即振幅屬性響應值相對較高的層段，砂巖較為發(fā)育(厚度大或砂地比高)。因此可利用地震屬性響應值，對砂巖厚度分布進行定量評價。

作為數(shù)理統(tǒng)計的一個分支，多元線性回歸分析起源于20世紀初[5-10]，隨后其理論迅速得到完善，在多個領域中廣泛應用，并于本世紀初期應用于地震屬性融合中[11-12]。應用多元線性回歸的數(shù)學方法可以較好地建立各種地震屬性與儲層參數(shù)之間的關系[12]。因此，在前人認識的基礎上，這里將采用地震屬性多元線性回歸融合技術定量預測烏東斜坡帶南屯組砂巖儲層分布特點，為進一步認識烏東斜坡帶地區(qū)南屯組儲層展布特征，及對隱蔽圈閉識別預測打下良好的基礎。

1 地質(zhì)概況

研究區(qū)位于海拉爾盆地烏爾遜凹陷烏南次凹的東部斜坡帶，西起于烏南次凹中心，東至烏爾遜凹陷邊界部位(圖1)，構造上由西向東總體表現(xiàn)為從洼槽到斷裂坡折帶到緩坡帶逐漸抬升的單斜。研究區(qū)南屯組是重要的產(chǎn)油層之一，分為南一段和南二段兩段。截止至2008年，烏南次凹共有探井35口，評價井19口，其中南一段工業(yè)油流井9口，南二段工業(yè)油流井5口，工業(yè)氣流井1口，油藏類型多是巖性、斷層-巖性和斷塊油氣藏[13]。儲層預測精度是制約南二段油藏認識的關鍵因素。故本次研究旨在探索提高南二段砂巖儲層定量預測精度的方法。

圖1 研究區(qū)構造位置圖 (據(jù)大慶油田呼倫貝爾分公司，2014)Fig.1 Structual location of research area

烏南地區(qū)南二段主要發(fā)育扇三角洲—辮狀河—湖泊沉積體系，其頂部Ⅰ油組中，工業(yè)油氣產(chǎn)層主要分布在辮狀河三角洲前緣亞相和扇三角洲前緣亞相[14]；巖性以泥質(zhì)粉砂巖和泥巖為主(圖2)。南二段Ⅰ油組儲層孔隙度主要分布在12%～20%之間，平均約為16.2%，滲透率主要分布區(qū)間為1 mD～50 mD，合計占41.25%，平均約為70 mD。在平面上Ⅰ油組孔隙度具有明顯的分帶性，烏30井區(qū)相對較高，孔隙度主要分布在12%～18%，在烏東斜坡和南部物源體系控制的地區(qū)物性相對較好，油層比較發(fā)育[15]。經(jīng)綜合研究認為烏東斜坡帶南二段Ⅰ油組砂巖儲層性能良好，為中等儲層。

圖2 烏31井單井綜合柱狀圖Fig.2 Cores and log curves characteristics of well Wu31

通過已鉆井南屯組砂巖儲層單層厚度的統(tǒng)計分析，可見南屯組砂巖儲層單層厚度較薄，大部分厚度在1 m～5 m左右，超過61%的單層厚度在1 m～2 m左右(圖3)。而在研究區(qū)勘探開發(fā)過程中，也存在砂巖儲層單層厚度薄、橫向變化快的問題，為該區(qū)儲層預測存在的主要問題。針對這種薄儲層利用常規(guī)波阻抗反演技術分辨率達不到，而單一屬性預測結果又具有較強的多解性[16]。利用地震屬性融合的方法，能夠有效提高此類薄砂巖儲層預測精度[17]。因此，基于南屯組儲層的特點，這里采用地震屬性多元線性回歸融合技術進行儲層厚度預測。

圖3 南屯組砂巖儲層單層厚度統(tǒng)計Fig.3 Thickness distribution charts of sandstone reservoir in Nantun formation

2 地震屬性融合砂巖預測

2.1 地震屬性多元線性回歸融合技術

地震屬性是由地震數(shù)據(jù)衍生出的幾何形態(tài)、運動學、動力學及統(tǒng)計特征的特別量度。地震屬性大多反映某些儲層地質(zhì)條件的變化[18]，但使用單一屬性分析并預測儲層往往會帶來多解性。地震屬性融合技術考慮了不同屬性的物理意義，選取表征不同儲層特征的屬性，將多個屬性經(jīng)過一定的數(shù)學運算融合到一起，以降低儲層預測的多解性，進一步提高儲層預測的精度[19]。

多元線性回歸是數(shù)理統(tǒng)計的一個分支，應用多元線性回歸方法可以建立多種地震屬性與儲層參數(shù)之間的關系。設可預測的隨機變量為y，它受到p個非隨機因素x1，x2，…xp-1，xp，和不可預測的隨機因素ε的影響。多元線性回歸的數(shù)學模型為[20]：

y=β0+β1x1+β2x2+…+βp-1xp-1+βpxp+ε

ε～N(0,σ2)

(1)

其中:β0、β1、…、βp為回歸系數(shù)。

設砂巖厚度y對p種地震屬性x1、x2、…、xp的p元線性回歸方程為：

(2)

其中:b0、b1、…、bp為β0、β1、…、βp的最小二乘估計值。即b0、b1、…、bp應使砂巖厚度觀測值y與砂巖厚度回歸估計值的偏差平方和最小。

(3)

Q為b0、b1、…、bp的p+1元函數(shù)。經(jīng)整理可得到關于偏回歸系數(shù)b0、b1、…、bp的正規(guī)方程組，進而解正規(guī)方程組即可得到偏回歸系數(shù)b0、b1、…、bp的解，于是解出p元線性回歸方程(2)[21]。

2.2 地震屬性的提取與優(yōu)選

地震屬性多元線性回歸融合技術受敏感屬性的約束作用較大，地震屬性的優(yōu)選對預測結果的準確性起到很大的影響[22]。精細標定烏東斜坡帶南二段Ⅰ油組油氣層的頂?shù)捉缑?，并對其頂?shù)走M行精細解釋。在圖4中，烏154-68井南二段Ⅰ油組頂部砂巖儲層不發(fā)育，在地震剖面上表現(xiàn)為弱反射、弱振幅；烏31井南二段Ⅰ油組頂部砂巖儲層發(fā)育，在地震剖面上表現(xiàn)為強反射、強振幅。

圖4 烏31井—烏154-68井南屯組地震反射剖面Fig.4 Seismic reflection section of Nantun formation of well Wu31～Wu154-68

由此可見，剖面上地震屬性與砂巖儲層發(fā)育程度存在一定的正相關性。選擇合理的時窗，提取南二段Ⅰ油組頂?shù)字g三維地震數(shù)據(jù)體的振幅、頻率、相位、能量和相關等7類20余種地震屬性，對這些地震屬性進行優(yōu)選，其中物理意義相同的屬性保留相關性高的，再把相關系數(shù)小于0.5的屬性剔除掉，這樣就得到了物理意義不同、能夠反映儲層本質(zhì)特征的并且與已鉆井南二段Ⅰ油組砂巖儲層厚度相關性高的屬性。經(jīng)過分析，最終優(yōu)選出5個敏感的地震屬性(圖5)。

2.3 地震多屬性融合砂巖厚度預測

通過對烏東斜坡帶南屯組砂巖儲層發(fā)育井的巖心、鉆井及錄井資料展開分析研究，統(tǒng)計了82口井在南二段Ⅰ油組解釋時窗內(nèi)砂巖發(fā)育厚度。將其中56口井點處地震屬性標準化并代入理論模型，計算出回歸系數(shù)，并將回歸系數(shù)帶入理論模型，余下26口井作為后驗井，得到砂巖厚度預測模型如下：

圖5 優(yōu)選出的敏感屬性及預測厚度平面分布圖Fig.5 Sensitive attributes against sandstone reservoir and thickness map as prediction result(a)平均反射強度;(b)平均能量;(c)均方根振幅;(d)總絕對振幅;(e)弧長

ym= 135.65-0.525×Arms-0.251×

Astr+0.151×Aabs+0.427×Aene-

0.062×Aarc

(4)

式中：ym為預測砂巖厚度；Arms為均方根振幅；Astr為平均反射強度；Aabs為總絕對振幅；Aene為平均能量；Aarc為弧長。

計算砂巖預測模型(4)預測井點處砂巖厚度與井點處南二段Ⅰ油組砂巖厚度之間的誤差ε。把ε代回理論模型(1)得到了校正誤差后的砂巖厚度預測模型(5)。

ym= 135.65-0.525×Arms-0.251×

Astr+0.151×Aabs+0.427×Aene-

0.062×Aarc+ε

(5)

把烏東斜坡帶標準化好的地震屬性帶入校正后的模型(5)，即把預測砂巖厚度與地震屬性的關系在整個研究區(qū)推開，得到烏東斜坡帶地區(qū)南二段Ⅰ油組砂巖厚度平面分布(圖6)。

圖6 南二段Ⅰ油組砂巖厚度預測圖Fig.6 Sandstone thickness prediction

3 預測結果分析

對于南二段Ⅰ油組預測結果的合理性和準確性，本次研究從井上厚度符合率以及平面地質(zhì)認識兩個角度進行了分析。把研究區(qū)內(nèi)參與地震屬性融合的56口井南二段Ⅰ油組井上砂巖厚度與屬性融合預測砂巖厚度進行交會(圖7)，相關系數(shù)達到0.88。對比26口后驗井預測厚度與井上厚度，發(fā)現(xiàn)二者基本一致(圖8)。

在平面上南二段Ⅰ油組砂巖主要分布于研究區(qū)東部緩坡構造帶上，整體上砂體有呈北東—南西展布的趨勢。砂巖厚度預測結果顯示，研究區(qū)北部砂巖普遍發(fā)育，其中東北部存在一個面積較大的扇狀砂巖發(fā)育區(qū)，其應受到東北部物源控制；東南側局部砂巖呈近條帶狀形態(tài)展布；西北部砂巖薄而連片。對比該區(qū)南二段沉積相圖(圖9)可以發(fā)現(xiàn)：整體輪廓吻合度較高，富砂區(qū)與南屯組發(fā)育的扇三角洲沉積形態(tài)和范圍一致；而研究區(qū)中西部砂巖厚度較小，為主要的湖相發(fā)育區(qū)；西北部砂巖呈三角洲前緣席狀砂發(fā)育特點。主要物源方向與沉積上的認識契合，說明了地震屬性多元線性回歸融合技術，在該區(qū)南屯組砂巖儲層預測中具有良好的適用性。

圖7 儲層預測結果分析Fig.7 Analyses of reservoir prediction result

圖8 后驗井預測情況分析Fig.8 Analyses of check wells prediction result

圖9 研究區(qū)南二段沉積相圖(據(jù)鄧海，有改動)Fig.9 Sediment facies of N2 member in research area

4 結論

1)平面上烏東斜坡帶南二段Ⅰ油組砂巖主要分布于東部緩坡構造帶上，整體上砂體有呈北東—南西展布的趨勢。北部砂巖普遍發(fā)育，其中東北部存在一個面積較大的扇狀砂巖發(fā)育區(qū)，其應受到東北部物源控制；東南側局部砂巖呈近條帶狀形態(tài)展布；西北部砂巖薄而連片。砂巖分布特征與沉積上的認識一致。

2)通過地震屬性多元線性回歸融合技術，對烏東斜坡帶南二段砂巖儲層開展儲層厚度預測，優(yōu)選出平均反射強度、平均能量、均方根振幅、總絕對振幅和弧長5個南二段砂巖敏感的地震屬性。其預測結果較好，能夠清楚地反映烏東斜坡帶南屯組砂巖儲層厚度變化特征，有效提高了預測精度。該技術在研究區(qū)南屯組砂巖儲層預測中有良好的適用性，值得推廣使用。

[1] 董萬百, 吳河勇, 齊玉林, 等. 海拉爾盆地坡折帶成藏類型與勘探突破方向[J]. 石油學報, 2006, 27(S1): 7-13.DONGWB,WUHY,QIYL,etal.reservoir-fomingtypesandexplorationtargetsinslopebreakbeltsofhailaerbasin[J].ActaPetroleiSinica, 2006，27(s1): 7-13.(InChinese)

[2] 張斌. 烏南次凹隱蔽圈閉識別方法研究[J]. 油氣藏評價與開發(fā), 2013, 2(6): 1-5.ZHANGB.ResearchontheidentificationmethodsofsubtletrapsinsouthWuerxunsub-sag[J].Reservoirevaluationanddevelopment, 2013,2(6): 1-5. (InChinese)

[3] 謝利華, 吳朝東, 孔凡忠, 等. 烏爾遜凹陷烏東斜坡南屯組沉積相及有利儲集相帶[J]. 大慶石油地質(zhì)與開發(fā), 2011, 30(3): 1-7.XIELH,WUCD,KONGFZ,etal.Sedimentaryfaciesandfavorablereservoirfaciesbeltsof;lowerZuunbayanformationinwudongslopeinwuerxunsag[J].PetroleumGeologyandOilfieldDevelopmentinDaqing, 2011,30(3):1-7.(InChinese)

[4] 徐占東. 烏爾遜凹陷烏南地區(qū)南屯組儲層預測[D]. 大慶石油學院, 2010.XUZD.StudyonreservoirpredictionofNantungroupinsouthernareaofwuerxunhollow[D].DaqingPetroleumInstitute, 2010. (InChinese)

[5]FISHERRA.Thefoundationsoftheoreticalstatistics[J].PhilosophicalTransactionsoftheRoyalSocietyaMathematicalPhysical&EngineeringSciences, 1922, 222(1): 309-368.

[6]FISHERRA.Theoryofstatisticalestimation[J].MathematicalProceedingsoftheCambridgePhilosophicalSociety.CambridgeUniversityPress,1925,22(05):700-725.

[7]HOTELLINGH.Ageneralizedttestandmeasureofmultivariatedispersion[J].ProceedingsoftheSecondBerkeleySymposiuminMathematicsStatisticsandProbability,UniversityofCaliforniaLosAngeles&Berkeley, 1951, (2): 23-41.

[8]HOTELLINGH.Therelationsofthenewermultivariatestatisticalmethodstofactoranalysis[J].BritishJournalofStatisticalPsychology, 1957, 10(2): 69-79.

[9]ROYSN.Onaheuristicmethodoftestconstructionanditsuseinmultivariateanalysis[J].AnnalsofMathematicalStatistics, 1953, 24(2):220-238.

[10]ROYSN,BARGMANNRE.Testsofmultipleindependenceandtheassociatedconfidencebounds[J].AnnalsofMathematicalStatistics, 1958, 29(2): 491-503.

[11]季玉新, 歐欽. 優(yōu)選地震屬性預測儲層參數(shù)方法及應用研究[J]. 石油地球物理勘探, 2003, 38(1): 5726-5726.JIYX,OUQ.Researchonmethodandapplicationofoptimizingseismicattributetopredictreservoirparameter[J].OilGeophysicalProspecting, 2003, 38(1): 5726-5726. (InChinese)

[12]謝雄舉, 季玉新. 優(yōu)選地震屬性預測儲層參數(shù)的方法及其應用[J]. 石油物探, 2004 (S1): 127-131.XIEXJ,JIYX.Methodandapplicationofoptimizingseismicattributetopredictreservoirparameter[J].GeophysicalProspectingforPetroleum, 2004 (S1): 127-131. (InChinese)

[13]王雅春, 張曉玲, 袁紅旗, 等. 烏東斜坡帶南一段油富集的主控因素分析[J]. 山東科技大學學報 (自然科學版), 2012(3): 1-10.WANGYC,ZHANGXL,YUNAHQ,etal.Dominativefactorsofoilenrichmentinthefirstsectorofnantungroupinwudongslopezone[J].JournalofShandongUniversityofScienceandTechnology(NaturalScience), 2012(3): 1-10. (InChinese)

[14]鮑園, 羅群, 韋重韜, 等. 海拉爾盆地烏南次凹南屯組有效儲集層形成條件[J]. 新疆石油地質(zhì), 2011, 32(5): 461-463.BAOY,LUOQ,WEICD,etal.Formingconditionsofeffectivereservoirofnantunformationinwunansub-sagofHailaerbasin[J].XinjiangPetroleumGeology, 2011, 32(5): 461-463. (InChinese)

[15]鄧海. 海拉爾盆地烏南次凹儲層綜合評價研究[J]. 內(nèi)蒙古石油化工, 2013, 39(16): 128-130.DENGH.Studyonreservoircomprehensiveevaluationofwunansaginhailaerbasin[J].InnerMongoliaPetrochemicalIndustry, 2013, 39(16): 128-130. (InChinese)

[16]秦月霜, 陳顯森. 用優(yōu)選后的地震屬性參數(shù)進行儲層預測[J]. 大慶石油地質(zhì)與開發(fā), 2000, 19(6): 44-45.QINYS,CHENXS.Performingareservoirpredictionwithoptimizedseismicattributeparameters[J].PetroleumGeologyandOilfieldDevelopmentinDaqing, 2000,19(6): 44-45.(InChinese)

[17]趙繼龍, 王俊鵬, 張先龍, 等. 基于地震多屬性融合與地質(zhì)統(tǒng)計學反演的薄層砂巖預測[J]. 地質(zhì)科技情報, 2014(05): 93-99.ZHAOJL,WANGJP,etal.Thinsandstonepredictionbasedonmulti-seismicattributeinformationfusionandgeostatisticalinversion[J].GeologicalScienceandTechnologicalInformation, 2014(05): 93-99. (InChinese)

[18]于建國, 姜秀清. 地震屬性優(yōu)化在儲層預測中的應用[J]. 石油與天然氣地質(zhì), 2003,24(3): 291-295.YUJG,JIANGXQ.ApplicationofSeismicAttributeOptimizationinReservoirPrediction[J].Oil&GasGeology, 2003,24(3): 291-295. (InChinese)

[19]李婷婷, 王釗, 馬世忠, 等. 地震屬性融合方法綜述[J]. 地球物理學進展, 2015(01): 378-385.LITT,WANGZ,MASZ,etal.summaryofseismicattributesfusionmethod[J].ProgressinGeophysics, 2015(01): 378-385. (InChinese)

[20]林彬. 多元線性回歸分析及其應用[J]. 中國科技信息, 2010(09): 60-61.LINB.Multiplelinearregressionanalysisanditsapplication[J].ChinaScienceandTechnologyInformation, 2010(9): 60-61. (InChinese)

[21]李操. 基于地震屬性的儲層預測方法研究及初步應用[D]. 大慶: 大慶石油學院, 2006.LIC,Reservoirparameterpredictedmethodresearchandprimaryapplicationbasingonseismicattribute[D].Daqing:DaqingPetroleumInstitute, 2006. (InChinese)

[22]張娟. 基于多元線性回歸分析的薄儲層預測技術在勝利探區(qū)的研究與應用[J]. 工程地球物理學報, 2013, 10(1): 91-94.ZHANGJ.ResearchandApplicationofthinreservoirpredictiontechniquesbasedonmultiplelinearregressioninShengliexplorationarea[J].ChineseJournalofEngineeringGeophysics, 2013, 10(1): 91-94. (InChinese)

Application of multi-seismic attribute information fusion in sandstone reservoir prediction in nantun formation in the eastern Wu slope in hailaer basin

ZHU Kedana,b, WANG Yachuna,b, YI Qifana,b, WANG Xinyaoa,b,LI Mengyaoa,b, ZHAO Zhenduoa,b, PENG Longyia,b

(a.College of Earth Science b.The state key laboratory base of unconventional oil and gas accumulation and exploitation, Northeast Petroleum University, Daqing 163318,China)

The characteristic of sandstone reservoir in Nantun formation in the eastern Wu slope in Hailaer basin is analyzed in this paper. It points out the exploration difficulties. Against the features and difficulties, the authors apply multiple linear regression in multi-seismic attribute information fusion to make an accurate thickness prediction of sandstone reservoir in N2I Bed. The general process of this method is as follows. Firstly, extract seismic attributes from 3D-seismic data, optimize sensitive attributes against sandstone reservoir among them. Then analyses the relationship between seismic attributes and thickness of sandstone reservoir. After that, establish the regression equation between seismic attributes and thickness of sandstone reservoir. Base on this regression equation, mapping the plane thickness distribution of sandstone reservoir in N2I bed. This method has good application results. With this result, the thickness feature of sandstone reservoir in N2I bed can be clearly recognized, which lay a solid foundation for the prediction of subtle traps in N2I bed in the eastern Wu slope in Hailaer basin.

eastern Wu slope, Nantun formation; multi-seismic attribute information fusion; multiple linear regression; sandstone reservoir prediction

2016-02-20 改回日期：2016-04-25

國家自然科學基金(41440017)；黑龍江省自然科學基金(D201406)

朱可丹(1991-)，男，碩士，研究方向為油氣田開發(fā)地質(zhì)，E-mail：redcloudszkd@163.com。

1001-1749(2017)01-0109-07

P 631.4

A

10.3969/j.issn.1001-1749.2017.01.16