摘　要：四川盆地元壩地區(qū)二疊系長興組發(fā)育大型的臺地邊緣生物礁灘氣藏，具有埋藏深、儲層空間展布復(fù)雜多變、一礁或一灘為一個氣藏、氣水分布極其復(fù)雜的特點，給勘探開發(fā)帶來了很大的風(fēng)險。為了精細刻畫元壩氣田生物礁灘氣藏儲層及氣水的展布特征，提高探井成功率，集成創(chuàng)新了儲層氣水識別技術(shù)：采用了疊前道集優(yōu)化技術(shù)提高了道集的信噪比；利用小角度疊前反演技術(shù)解決了超深層彈性參數(shù)的求取；最后從巖石物理參數(shù)分析入手，綜合實驗室?guī)r石物理測試分析及測井、試氣成果統(tǒng)計分析結(jié)果，明確了氣水識別的敏感參數(shù)；在此基礎(chǔ)上進一步利用坐標(biāo)旋轉(zhuǎn)技術(shù)構(gòu)建新的高靈敏度識別因子，從而提高了疊前敏感參數(shù)反演的識別能力和精度。在該地區(qū)l區(qū)塊實際應(yīng)用的結(jié)果表明，該技術(shù)有效刻畫出了元壩氣田氣區(qū)、水區(qū)以及氣水過渡區(qū)，預(yù)測結(jié)果與后期實際鉆井揭示的規(guī)律相一致。

關(guān)鍵詞：四川盆地  元壩氣田  超深層  生物礁灘氣藏  儲集層  敏感參數(shù)  氣水識別  石油地球物理勘探

Gas-water identification technologies for ultra deep reef reservoirs in the Yuanba Gas Field，Sichuan Basin

Abstract：Large platform margin reef gas reservoirs were discovered to be deep，one reef one reservoir，in the Permian Changxing Feixianguan Fm in the Yuanba area，Sichuan Basin．The exploratory drilling reveals that their spatial distributi。n is complex and changeable and their gas-water distribution is extremely complicated，which brings about a high risk in exploration and exploitation there．In order to improve the success ratio of exploratory wells，it is necessary to finely describe the characteristics of reef and beach gas reservoirs and gas water distribution in this study area．Therefore，technologies of gas-water identification were integrated and innovated．First，pre-stack gathers optimization techniques were used to improve the SNR．Then，the small angle stack inversion technology was employed to calculate the elastic paranleters．Finally，based on the analysis of petrophysical parameters from laboratory，logging and well test data，those sensitive gas water parameters were clearly identified．On this basis，the coordinate rotation technology was applied to build a new highly sensitive identification factor，thereby to improve the identificarlon abilitv and accuracy of pre-stack inversion of sensitive parameters．The application result from Block I in this study area indicated that the gas zone，water zone，and gas-water transition zone were precisely describcd，which is in good agreement with the later actual drilling data．

Keywords：Sichuan Basin，Yuanha Gas Field，ultra deep strata，reef reservoir，sensitive parameter，gas water identiffcation，petroleum geophysical prospecting

元壩氣田位于四川盆地東北部，氣藏中部埋深6673m，是目前我國已發(fā)現(xiàn)的埋藏較深的大氣田^[1]。自2007年YB1c1井鉆遇二疊系長興組頂部優(yōu)質(zhì)白云巖儲層以來，部署在該礁灘相上的多口探井均獲百萬立方米以上高產(chǎn)工業(yè)氣流，截至2013年底，累計探明天然氣儲量超過2000×10⁸m³。

在礁灘儲層預(yù)測方面，前期利用“相控三步法”反演技術(shù)實現(xiàn)了礁灘儲層的精細描述，住礁灘儲層氣水識別技術(shù)方面，開展了多種技術(shù)方法的探索，從試驗結(jié)果來看，AV0常規(guī)屬性分析^[2]、電阻率反演^[3]、DR反演^[4]、MEEI反演^[5-6]、瞬時頻譜分析技術(shù)^[7]、低頻伴影^[8]等方法取得了一定的效果，但這些技術(shù)方法還存在多解性和識別陷阱。因此，探索適合于超深層礁灘儲層的疊前氣水識別技術(shù)是實現(xiàn)元壩氣田高效勘探開發(fā)的當(dāng)務(wù)之急。

1　礁灘儲層特征及氣水識別技術(shù)

元壩氣田長興組生物礁儲層巖石類型主要為溶孔白云巖、白云巖、生屑灰?guī)r等，儲集空間類型包括溶孔、溶洞、晶間孔、裂縫及生物體腔溶孔。而長興組臺緣淺灘儲層巖石類型主要為白云質(zhì)灰?guī)r、含云質(zhì)灰?guī)r，儲集空間類型主要為孔隙和裂縫，呈相對中孔高滲特征。

疊前AVO反演作為AVO應(yīng)用技術(shù)的一個重要組成部分，是合理地提取隱藏在疊前地震信息中的巖性和流體信息的重要途徑之一。該區(qū)預(yù)測的難點主要體現(xiàn)在埋藏深帶來的一系列問題：在如此深的地層AVO效應(yīng)是否明顯?小入射角的疊前反演能否實現(xiàn)?氣水識別的精度能否提高?

針對超深層礁灘儲層預(yù)測的難點，利用方法創(chuàng)新與實踐結(jié)合，形成了一套超深層礁灘儲層氣水識別技術(shù)。即采用疊前道集優(yōu)化技術(shù)提升道集質(zhì)量，并利用小角度疊前反演技術(shù)解決超深層彈性參數(shù)的求取，通過對井中儲層段的流體類型進行巖石物理參數(shù)分析，找到了敏感彈性參數(shù)，明確其響應(yīng)特征，最后利用改進的坐標(biāo)旋轉(zhuǎn)技術(shù)來提高敏感參數(shù)識別精度，最終達到了超深層礁灘儲層氣水識別的目的。

1．1　疊前道集優(yōu)化技術(shù)

高質(zhì)量的疊前道集組合，是地震疊前氣水識別的基礎(chǔ)。疊前道集數(shù)據(jù)優(yōu)化的目的：①凈化道集背景；②拉平道集；③挖掘AVO特征。其要點包括：恢復(fù)和保持CRP道集中各道之間相對振幅關(guān)系，提高疊前道集的信噪比，從而提高道集數(shù)據(jù)的分辨率，尤其是要正確地確定反射界面的位置、剩余時差校正等，最大可能地消除非油氣或巖性因素引起的振幅隨炮檢距的變化，以獲得高信噪比、高分辨率和高保真度的疊前道集數(shù)據(jù)。元壩氣田海相礁灘儲層含氣AVO異常響應(yīng)特征為3類，但是入射角小，同時還受噪音等干擾，其特征不是很明顯(圖1左紅框內(nèi))。經(jīng)過對動校正后的道集進行優(yōu)化處理后，凸現(xiàn)了AVO響應(yīng)特征、提升了信噪比(圖1右紅框內(nèi))，為疊前精細反演提供了相對高質(zhì)量的道集數(shù)據(jù)。

1．2　小角度疊前反演技術(shù)

由于元壩氣田超深層天然氣藏受觀測系統(tǒng)和測線長度的限制，難以獲得大角度信息，僅能獲得小、中角度的信息。從實際資料來看，其入射角一般就27°左右，最大不超過30°。因此，需要研究小角度反演方法來求取彈性參數(shù)。

根據(jù)Mallick等人對合成地震數(shù)據(jù)進行彈性阻抗反演的研究表明：如果彈性阻抗反演時夾雜有少量噪音，則不能恢復(fù)初始速度和密度模型，但這種情況下仍可估算縱、橫波阻抗。Lu Shaoming和McMechan^[9]基于此原理，利用小角度tan²q≈sin²q時得到近似式與經(jīng)驗關(guān)系式：

當(dāng)K小于0.25時，系數(shù)a＝8.0、b＝0.5；當(dāng)K大于0.25時，a＝3.0、b＝3.0。Lu等人指出，對于所有介于0和0.8之間的K值，式(1)是一個合理的近似式^[9-10]。這樣，在求取縱橫波阻抗時，只需要中小角度信息便可，因而適合小角度入射情況下的參數(shù)求取。根據(jù)元壩地區(qū)測井?dāng)?shù)據(jù)得到經(jīng)驗關(guān)系，發(fā)現(xiàn)該區(qū)K值基本上大于0.25，此時系數(shù)a＝1.41、b＝2.65。

采用該方法進行疊前反演，即可得到縱波阻抗，橫波阻抗，根據(jù)關(guān)系式：

lr＝Ip²-2Is²    　　　　　　(2)

mr＝Is² 　　　　　　　　　 (3)

便可求得該區(qū)敏感彈性參數(shù)。

1．3　敏感參數(shù)優(yōu)選

通過對元壩地區(qū)礁灘相儲層巖心樣品的巖石物理測試與分析，獲得了在地層溫壓等條件下的縱橫波速度、動靜彈性參數(shù)、品質(zhì)因子等大量的物性參數(shù)數(shù)據(jù)，并建立了氣水識別敏感因子(圖2)，可知對該地區(qū)氣水識別比較敏感的因子依次分別是：lr、kr、l、K等，其中lr為最敏感的識別因子。另外從該區(qū)鉆井資料的巖石物理彈性參數(shù)分析來看，lr也是對流體最為敏感(圖3)，圖3選取了1口典型的既產(chǎn)氣又產(chǎn)水，以及氣水同層的井。從圖3上可以看到，當(dāng)儲層含流體后，彈性參數(shù)均產(chǎn)生變化，當(dāng)儲層為氣層時，變化幅度較大，氣水同層次之，當(dāng)儲層為水層時，變化幅度相對前兩者來說最小。通過對各個彈性參數(shù)相對變化率的計算，lr的變化率為最大。因此在該地區(qū)的氣水識別可通過上述的疊前小角度反演方法獲取lr彈性參數(shù)體，進而預(yù)測氣水的展布。

1．4　坐標(biāo)旋轉(zhuǎn)技術(shù)

Goodway認為Lambda(l)和Mu(m)屬性的物理解釋是：l屬性即不可壓縮性，對孔狀流體比較敏感；而m屬性即硬度，對巖石基質(zhì)(礦物成分)比較敏感。因此，對基于參數(shù)l、m和r提出了一種新的方法，稱為Lambda-Mu-Rho方法^[11]。

圖4為元壩地區(qū)長興組礁灘儲層段井上Lambda-Mu-Rho參數(shù)交會圖。從圖4中可以看到：雖然敏感彈性參數(shù)lr能將氣層、水層以及氣水同層基本區(qū)分開，但是識別精度還不夠。因此，為提高lr參數(shù)的識別精度，采用了坐標(biāo)選擇的方法，其基本原理如下。

采用數(shù)學(xué)上的坐標(biāo)系轉(zhuǎn)換方法——以斜直線與橫軸之間的夾角(q)為旋轉(zhuǎn)角度進行坐標(biāo)系旋轉(zhuǎn)，得到新的坐標(biāo)系(旋轉(zhuǎn)后的坐標(biāo)系橫軸與斜直線平行，縱軸與該斜直線垂直)，兩坐標(biāo)系的轉(zhuǎn)換關(guān)系如下所示(x、y為原坐標(biāo)系中的坐標(biāo)；x¢、y¢為相應(yīng)新坐標(biāo)系的坐標(biāo))，即

x¢＝xcosq+ysinq             (4)

y¢＝ycosq+xsinq             (5)

順時針旋轉(zhuǎn)方向時，使用式(4)；逆時針旋轉(zhuǎn)時，使用式(5)。這樣利用新舊坐標(biāo)系之間的數(shù)學(xué)轉(zhuǎn)換關(guān)系，就可以求出新參數(shù)與原參數(shù)的函數(shù)關(guān)系。

通過對圖4-a進行坐標(biāo)旋轉(zhuǎn)后，氣層、氣水同層、水層分布的范圍區(qū)分更加明顯，構(gòu)建的新敏感因子(lr¢)對氣水的識別能力得到了顯著的提高，而且還可以確定氣水識別的門檻值進行定量預(yù)測，其中含氣層小于100GPa·g／cm³，含水層大于115GPa·g／cm³，氣水同層則介于兩者之間。由此可見，利用旋轉(zhuǎn)后的新敏感因子將大大提高氣水識別的精度。

2　應(yīng)用效果分析

利用該技術(shù)完成了元壩氣田東部區(qū)域的氣水檢測，從識別的結(jié)果來看，與實鉆井的測試情況基本一致，也證實了該方法在該地區(qū)的有效性。圖5為該區(qū)過YBl04、YBll、YBl23、Ybl0cl井連井氣水識別剖面圖。其中YBl04測試獲得氣123×10⁴m³／d，YBll測試獲得氣51×10⁴m³／d，YBl23測試獲得氣5.3×10⁴m³／d，水290m³／d，Ybl0cl測試獲得氣107×10⁴m³／d，從圖5中可以看到，從產(chǎn)氣的YBl04井、YBll井到處于氣水過渡帶的YBl23井，穿過一段灘，間水區(qū)然后一直到產(chǎn)氣的YB10cl井，其氣區(qū)、水區(qū)分布帶及過渡帶清晰，與實鉆情況吻合好。

另外從預(yù)測的元壩氣田l區(qū)塊長興組氣水分布平面圖(圖6)來看，整體上有利含氣區(qū)主要分布在西邊及北邊，東南邊多為產(chǎn)水區(qū)及氣水同層區(qū)，其刻畫的氣水平面展布與實鉆井揭示的規(guī)律基本上保持一致。從圖6中可以看到：含氣區(qū)(紅黃色區(qū)域)與含水區(qū)(藍色區(qū)域)分布清晰。處于黃色區(qū)域(紅圈內(nèi))的井，均獲得了工業(yè)性氣流。尤其處于紅黃色富集區(qū)域的YBl04井以及YB10cl井，獲得了日產(chǎn)天然氣過百萬立方米的產(chǎn)能。藍色區(qū)域則為含水區(qū)(圖中黑圈內(nèi))，處于該區(qū)域的YB9井在長興組鉆遇水層，日產(chǎn)水30m³。綠色區(qū)域為氣水同層區(qū)，該區(qū)域的YBl23井、YBl6井測試結(jié)果既產(chǎn)氣也產(chǎn)水。

在圖6中，識別結(jié)果顯示YB211井與YBl07井也同樣會穩(wěn)獲工業(yè)氣流，這與該井的測井解釋成果吻合，近期測試結(jié)果顯示YB211井產(chǎn)氣5.63×10⁴m³／d，YBl07井為38.35×10⁴m³／d，進一步的說明了預(yù)測結(jié)果的可靠性以及該方法的適應(yīng)性。

3　結(jié)論

該技術(shù)在元壩氣田長興組—飛仙關(guān)組超深層礁灘儲層氣水識別中取得了良好的效果，提高了探井的成功率和經(jīng)濟效益，主要得到以下結(jié)論：

1)對于超深層地區(qū)的疊前道集，盡可能地開展疊前道集預(yù)處理，從而改善道集的質(zhì)量，為后期精細疊前反演打好基礎(chǔ)。

2)小角度疊前反演技術(shù)的應(yīng)用可彌補由于缺失大角度信息而導(dǎo)致的方法適用性問題。

3)坐標(biāo)旋轉(zhuǎn)技術(shù)能夠在巖石物理分析獲得的敏感因子的基礎(chǔ)上，構(gòu)建更為靈敏的流體識別因子，從而可提高疊前反演對氣水的識別能力和精度。

4)該方法技術(shù)在元壩氣田超深層礁灘氣藏儲層氣水的識別中取得了良好的效果，為后續(xù)的勘探開發(fā)提供了重要的支撐，在四川盆地及同類地區(qū)具有一定的推廣應(yīng)用價值。

參考文獻

[1]郭旭升，胡東風(fēng)．川東北礁灘天然氣勘探新進展及關(guān)鍵技術(shù)[J]．天然氣工業(yè)，20ll，3l(10)：6-11．

GUO Xusheng，Hu Dongfeng．Newest progress and key techniques in gas exploration of reef-bank reservoirs in the northeastern Sichuan Basin[J]．Natural gas industry，2011，31(10)：6-11．

[2]郭旭升，凡睿．AVO技術(shù)在普光氣田鮞灘儲層預(yù)測中的應(yīng)用[J]．石油與天然氣地質(zhì)，2007，28(2)：198-202．

GUO Xusheng，FAN Rui．Application of AVO techniques to reservoirs prediction of oolitic beach facies in Puguang Gas Field[J]．Oil&Gas Geology，2007，28(2)：198-202．

[3]黃先平，張健，謝芳，等．鮞灘儲層預(yù)測技術(shù)、流體識別方法在川東北地區(qū)飛仙關(guān)天然氣勘探中的應(yīng)用[J]，中國石油勘探，2003，8(1)：41-51．

HUANG Xianping，ZHANG Jian，XIE Fang，et al．Application of prediction technology about oolitic beach reservoir and fluid identification method in exploration of Tif natural gas in the Northeast of Sichuan Basin[J]．China Petroleum Exploration，2003，8(1)：41-51．

[4]孫巧玲．川東北地區(qū)飛仙關(guān)組鮞灘氣藏氣水識別[J]．天然氣勘探與開發(fā)，2005，28(6)：20-26．

SUN Qiaoling．Gas-water discernnlent of Feixianguan Formarion oolitic beach gas reservoir in the northeast of Sichuan Basin[J]．Natural Gas Exploration and Development．2005，28(6)：20-26．

[5]彭真明，李亞林，巫盛洪，等．碳酸鹽巖儲層多角度彈性阻抗識別方法[J]．地球物理學(xué)報，2008，51(3)：881-885．

PENG Zhenming，LI Yalin，WU Shenghong，et al．Discriminating gas and water using multi-angle extended elastic impedance inversion in carbonate reservoirs[J]．Chinese Journal of Geophysics，2008，5l(3)：881-885．

[6]胡偉光，李發(fā)貴，楊鴻飛．疊前彈性波阻抗反演在四川FL地區(qū)礁灘型儲層預(yù)測中的應(yīng)用[J]．海相油氣地質(zhì)，2010，15(4)：62-67．

HU Weiguang，LI Fagui，YANG Hongfei．Application of prediction of pre-stack elastic impedance inversion in reefshoal reservoir at FL Area，Sichuan Basin[J]．Marine Origin Petroleum Geology，2010，15(4)：62-67．

[7]徐天吉，程冰潔．瞬時頻譜分析技術(shù)與氣水識別方法研究[J]．石油天然氣學(xué)報，2008，30(4)：73-78．

XU Tianji，CHENG Bingjie．Research of instantaneous spectral analysis technology and application in distinguishing water from gas reservoir[J]．Journal of Oil and Gas Technology，2008，30(4)：73-78．

[8]HE Zhenhua，XIONG Xiaojun，BIAN Li¢en．Numerica lsimulation of seismic low frequency shadows and its application[J]．Applied Geophysics，2008，5(4)：301-306．

[9]LU Shaoming，MCMECHAN G A．Elastic impedance in version of multichannel seismic data from unconsolldated sediments containing gas hydrate and free gas[J]．Geophysics，2004，69(1)：164-179．

[10]李愛山，印興耀，陸娜，等．兩個角度彈性阻抗反演在中深層含氣儲層預(yù)測中的應(yīng)用[J]．石油地球物理勘探，2009，44(1)：87-92．

LI Aishan，YIN Xingyao，LU Na，et al．Application of elastic impedance inversion with two angle stack gathers to predict gas bearing reservoir of mid-deep layer[J]．Oil Geophysical Prospecting，2009，44(1)：87-92．

[11]CHOPRA S，PRUDEN D．Multiattribute seismic analysis on AVO-derived parameters-a case study[J]．The Leading Edge，2003，22(10)：998．

本文作者：尹正武

作者單位：中國石化勘探南方分公司