頻散AVO分析技術(shù)及其在南海深水區(qū)的應(yīng)用*

張益明 牛 聰 韓 利 葉云飛

(中海油研究總院 北京 100028)

頻散AVO分析技術(shù)及其在南海深水區(qū)的應(yīng)用*

張益明 牛 聰 韓 利 葉云飛

(中海油研究總院 北京 100028)

張益明,牛聰,韓利，等.頻散AVO分析技術(shù)及其在南海深水區(qū)的應(yīng)用[J].中國海上油氣，2016,28(6)：21-27.

Zhang Yiming,Niu Cong,Han Li,et al.A method of frequency dispersion AVO analysis and its application in deep water area of South China Sea[J].China Offshore Oil and Gas,2016,28(6):21-27.

深水區(qū)勘探鉆井稀少，由于缺乏測井約束，地震反演等常規(guī)技術(shù)油氣檢測的多解性很強?；诘卣鸩ㄋp的頻散現(xiàn)象，在地震譜分解基礎(chǔ)上，提出應(yīng)用頻率依賴的AVO方程計算地震波的頻散屬性的方法，并設(shè)計理論模型驗證了這種頻散屬性對含氣性識別的高敏感度，進而將AVO屬性與頻散屬性相結(jié)合構(gòu)建了一種烴類指示因子(DvP×DG)。應(yīng)用頻散AVO分析技術(shù)對南海深水區(qū)Y21-1構(gòu)造區(qū)目的層系含氣性進行了鉆前預(yù)測，結(jié)果與實鉆結(jié)果吻合較好，證實了頻散AVO分析技術(shù)可以增強氣層識別的敏感度，降低水層等背景帶來的影響，可以作為烴類檢測指標進行含氣性識別。

地震波衰減；譜分解；AVO方程；頻散屬性；含氣性識別；南海深水區(qū)

在深水等勘探程度較低的區(qū)域，鉆井資料少，如何直接應(yīng)用地震資料進行油氣預(yù)測一直是地球物理學家探索的方向之一。巖石中流體的存在會導致地震波發(fā)生頻散和衰減，并且頻散和衰減與地震波的頻率有關(guān)，這種現(xiàn)象已被巖石物理實驗證實[1-5]。因此，理論上這種由流體的存在而引起的地震波頻散和衰減可以用于地下巖石含流體性預(yù)測。

針對速度頻散與地震反射特征以及AVO之間的關(guān)系，一些學者開展了研究。Chapman等[5]指出在彈性介質(zhì)與頻散介質(zhì)界面處反射與頻率有關(guān)。Ren等[6]討論了垂直入射時彈性介質(zhì)與頻散介質(zhì)界面處的反射系數(shù)振幅與相位隨頻率的變化關(guān)系。Johnston[7]、Sun[8]等證實了速度頻散與地層物性密切相關(guān)，并指出頻散對地層物性的變化比速度更敏感。高剛、吳小羊、張震 等[9-11]開展了速度頻散方面的反演和應(yīng)用研究，探索了利用頻散屬性開展含油氣性分析的條件和可行性。

以上研究表明，地震波在含油氣儲層中的速度頻散和衰減可作為一種流體識別指標，如果地層速度隨頻率變化，那么相應(yīng)的反射系數(shù)也會隨頻率變化，這為應(yīng)用AVO方程計算頻散特性提供了理論依據(jù)。本文基于模型數(shù)據(jù)驗證了頻散AVO分析的有效性，通過計算地震波的頻散屬性進一步分析了氣層響應(yīng)特征與頻率變化的關(guān)系，進而將頻散屬性與AVO屬性相結(jié)合構(gòu)建了烴類指示因子(DvP×DG)，并應(yīng)用該技術(shù)對南海深水Y21-1構(gòu)造區(qū)目的層系含氣性進行了預(yù)測，指出了含氣區(qū)的分布范圍，得到了實鉆結(jié)果的證實。

1 頻散AVO分析原理

按照Shuey、Ostrander、殷八斤 等[12-15]的研究，各向同性介質(zhì)P波反射系數(shù)在入射角小于30°時，AVO截距梯度近似公式為

R(θ)≈P+Gsin2θ

(1)

式(1)中：θ為入射角；P為縱波法線入射的反射系數(shù)，稱為AVO截距；G為AVO梯度。

重新整理式(1)截距項P，得到新的近似表達式

(2)

式(2)中：vP、ρ分別為界面兩側(cè)介質(zhì)的縱波速度和密度平均值；ΔvP、Δρ分別為vP、ρ的差異。

考慮地震波速度隨頻率變化，可將式(2)反射系數(shù)擴展成入射角和頻率的函數(shù),即

(3)

對上述方程的頻率做一階泰勒展開，可以得到

(4)

(5)

RP(θ,f)=RP(θ,f0)+DvPΔf+DGΔfsin2θ

(6)

式(6)中：DvP和DG作為烴類指示參數(shù)(即反演的頻散屬性)，其中DvP表征縱波速度頻散程度，DG相當于AVO梯度，表征縱波頻散程度隨炮檢距的變化梯度。

2 頻散屬性模擬分析

為了驗證頻散屬性進行油氣檢測的可行性，基于Chapman、Wilson、Smith等[5,16-17]的研究設(shè)計一個理想模型進行測試。該模型假設(shè)巖石介質(zhì)處于某種孔隙流體(通常為水)的飽和狀態(tài)，當發(fā)生另一種流體(通常為氣)的替換時可以得到新的具有頻率依賴性質(zhì)的參數(shù)，當介質(zhì)孔隙壓力均衡時具有彈性屬性，且在低頻端(即地震頻帶內(nèi))不發(fā)生頻散。

如圖1、表1所示，該模型上層為不含流體的泥巖，近似于完全彈性，下層為砂巖，左側(cè)為水飽和，右側(cè)為氣飽和，使得模型砂巖具有頻散屬性。假設(shè)該模型中水體積模量為2 000 MPa，氣體積模量為200 MPa。圖2為該理論模型的地震正演響應(yīng)特征，對比圖2a、b可以看到，氣飽和頻散介質(zhì)中反射波振幅隨偏移距增加而增大。

圖1 理論模型

表1 理論模型參數(shù)

對該模型的合成地震記錄進行分頻AVO研究，頻率范圍10～90 Hz，間隔10 Hz。圖3給出了該模型彈性介質(zhì)和頻散介質(zhì)中零偏移距道分頻反演得到的縱波速度變化率(ΔvP/vP)隨頻率的變化情況，對振幅譜進行歸一化處理后可以看出，與彈性介質(zhì)相比，頻散介質(zhì)縱波速度變化率隨頻率變化的特征明顯。由此可見，速度相對變化量隨頻率變化的屬性可用于表征儲層的頻散程度，進而進行烴類檢測。

對該模型水飽和砂巖與氣飽和砂巖分別計算頻散屬性，結(jié)果表明無論是DvP還是DG，氣飽和頻散屬性都明顯大于水飽和頻散屬性(圖4)，這說明這種頻散屬性對含氣性的識別具有更高的敏感度，可作為含氣性的檢測指標。

圖2 理論模型正演地震響應(yīng)

圖3 理論模型零偏移距彈性介質(zhì)和頻散介質(zhì)的速度變化率分析

圖4 理論模型不同條件下的頻散屬性

3 在南海深水區(qū)的應(yīng)用

3.1 區(qū)域概況

研究目標位于南海深水區(qū)Y21-1構(gòu)造區(qū)。前期地質(zhì)研究認為，該構(gòu)造是在基底隆起上發(fā)育的披覆構(gòu)造，主要目的層為珠海組上段深水扇砂體，次要目的層為珠江組上段濁積砂體，構(gòu)造面積由深到淺逐漸變大，幅度逐漸減小。

應(yīng)用頻散AVO分析方法對Y21-1構(gòu)造區(qū)珠江組和珠海組儲層分布及儲層含流體性質(zhì)進行預(yù)測，以期降低勘探風險，為勘探?jīng)Q策提供依據(jù)。目前該構(gòu)造區(qū)勘探程度較低，區(qū)內(nèi)無鉆井資料，無法基于地震反演開展相應(yīng)的油氣識別研究，但是研究區(qū)內(nèi)具有2011年采集的3D地震資料，地震品質(zhì)較高(過目標體實際地震疊加剖面如圖5所示，其中強振幅異常處即為珠江組—珠海組深海重力流沉積體)，這為利用疊前地震數(shù)據(jù)開展頻散AVO分析工作提供了必要的數(shù)據(jù)基礎(chǔ)。

圖5 Y21-1構(gòu)造區(qū)地震剖面

3.2 實現(xiàn)流程

針對Y21-1構(gòu)造區(qū)地質(zhì)特點與資料條件，結(jié)合現(xiàn)有的頻散屬性提取方法的研究成果[18-25]，本文采用如下提取流程：

1) 對疊前地震記錄進行球面擴散補償、速度分析、構(gòu)建速度模型、動校正等預(yù)處理；

2) 用廣義S變換對預(yù)處理得到的地震記錄進行頻譜分解，獲得一系列不同頻率的等頻剖面[24-25]；

3) 提取子波，以子波主頻為參考頻率，求取譜振幅均衡系數(shù)，并對其余等頻率剖面進行譜均衡處理(這是因為譜分解以后的數(shù)據(jù)中有子波的影響，需要進行譜均衡來消除“子波疊印”對地層反射系數(shù)的影響)；

4) 選取地震子波主頻為參考，進行分頻AVO反演，得到頻散屬性。

3.3 含氣性預(yù)測效果

應(yīng)用頻率依賴的AVO方程計算地震波的頻散屬性，并與P×G屬性結(jié)合，開展了Y21-1構(gòu)造區(qū)目標儲層的含氣性鉆前預(yù)測研究。由于剪切模量對填充流體性質(zhì)的變化不敏感，因此只考慮P波頻散。

圖6為該構(gòu)造區(qū)常規(guī)AVO的P和G交會分析得到的含氣異常剖面，可以看出，該構(gòu)造頂部目標一氣層異常明顯，而下部目標二氣層只出現(xiàn)在構(gòu)造邊緣部位，主要原因是上覆厚氣層的屏蔽作用引起下部氣層的地震反射系數(shù)隨角度變化而發(fā)生了變化。分析表明，該構(gòu)造目標一氣層表現(xiàn)為III類AVO異常，目標二氣層表現(xiàn)為弱IV類AVO異常(圖7)，所以常規(guī)AVO屬性分析技術(shù)基本無法應(yīng)用。

圖6 Y21-1構(gòu)造區(qū)P和G交會剖面

圖7 Y21-1構(gòu)造區(qū)地震道集和AVO曲線

圖8是通過頻散AVO流體識別技術(shù)流程得到的該構(gòu)造區(qū)頻散屬性剖面，可以看出，在目標一和目標二處有2套較強的頻散屬性異常，分別對應(yīng)珠江組和珠海組目的層。由式(1)和式(6)可知，頻散介質(zhì)的AVO屬性包括P、G、DvP和DG，因此可以通過將P×G與DvP×DG交會來開展研究，圖9為只考慮氣層頂面情況下P×G與DvP×DG的交會結(jié)果，其中紅色區(qū)域表示P×G值較大，頻散程度DvP×DG值也較大，相當于常規(guī)AVO的III類異常；黃色區(qū)域表示P×G值較小，但頻散程度DvP×DG較大，相當于常規(guī)AVO的非III類異常；藍色區(qū)域表示背景，其中左上部藍色區(qū)域表示氣層底界面的響應(yīng)，下半部藍色區(qū)域表示泥巖和含水砂巖的響應(yīng)。 將圖9中的3個區(qū)域反投影到剖面，結(jié)果表明：對于目標一氣層，常規(guī)P×G就能夠識別；而對于目標二氣層，只有頻散屬性才能很好地識別，而且與圖6相比，構(gòu)造頂部目標二氣層特征更為完整(圖10)。由此可見，AVO屬性與頻散屬性的結(jié)合應(yīng)用可以增強氣層識別的敏感度，降低水層等背景帶來的干擾。

圖8 Y21-1構(gòu)造區(qū)頻散屬性剖面

圖9 Y21-1構(gòu)造區(qū)P×G與DvP×DG交會分析

圖10 Y21-1構(gòu)造區(qū)P×G與DvP×DG交會剖面

之后，Y21-1井鉆探揭示的2套目的層均為深水扇砂體沉積體。其中，上部目的層為深水扇遠端濁積砂體，以粉砂—細砂為主，泥質(zhì)含量偏高，測井解釋孔隙度在19%～25%之間，在該層段縱向識別出3套高氣測異常段，綜合測井解釋含氣層段累計厚度60 m；下部目的層為深水扇砂體，以中砂—粗砂為主，部分含礫，分選中等，為高孔隙優(yōu)質(zhì)儲層(平均孔隙度23%)，氣層異常較高，綜合測井解釋氣層累計厚度40 m。這表明，該井2套異常位置都鉆遇到了儲層物性較好的氣層，證實了基于頻散屬性分析的AVO含氣性預(yù)測結(jié)果。

4 結(jié)論

1) 考慮地震波傳播過程的速度頻散因素，提出了應(yīng)用頻散依賴的AVO方程計算地震波的頻散屬性的方法，并設(shè)計理論模型進行了測試，所得到的縱波速度頻散程度和頻散梯度屬性對氣層識別效果明顯。

2) 南海深水區(qū)Y21-1構(gòu)造區(qū)應(yīng)用表明，常規(guī)AVO屬性分析與頻散屬性的結(jié)合，既考慮了地震振幅隨角度的變化，又考慮了地震振幅的頻散影響，含氣性預(yù)測結(jié)果更為可靠，并且得到了實鉆結(jié)果的證實。

3) 考慮到影響速度頻散因素的非唯一性，基于頻散AVO屬性的含氣性預(yù)測結(jié)果可作為定性分析依據(jù)，同時應(yīng)用中要重視頻散屬性與地質(zhì)、測井和反演等多物探資料的綜合分析，以提高含氣性預(yù)測的可靠性。

[1] 王海洋,孫贊東,CHAPMAN M. 巖石中波傳播速度頻散與衰減[J].石油學報,2012,33(2):332-342. Wang Haiyang,Sun Zandong,CHAPMAN M.Velocity dispersion and attenuation of seismic wave propagation in rocks[J].Acta Petroleum Sinica,2012,33(2):332-342.

[2] 鄧繼新,史謌,俞軍.流體飽和巖石超聲速度頻散的特性研究[J].北京大學學報(自然科學版),2003，39(6):835-843. Deng Jixin,Shi Ge,Yu Jun.The research of the characteristics of ultrasonic velocity dispersion in fluid saturated rock[J].Acta Scientiarum Naturalium Universitaties Pekinensis,2003,39(6):835-843.

[3] BATZLE M,HOFMANN R,HAN D H,et al.Fluids and frequency dependent seismic velocity of rocks[J].The Leading Edge,2001,20(2):168-171.

[4] 李勝軍,高建虎,雍學善,等.小尺度體反射系數(shù)近似公式研究[J].巖性油氣藏,2014,26(1):96-99. Li Shengjun,Gao Jianhu,Yong Xueshan,et al.Approximate formula of reflection coefficient for small-scale body[J].Lithologic Reservoirs,2014,26(1):96-99.

[5] CHAPMAN M,LIU Enru,LI Xiangyang.The influence of fluid-sensitive dispersion and attenuation on AVO analysis[J].Geophysical Journal International,2006,167:89-105.

[6] REN H,GOLOSHUBIN G,HILTERMAN F J.Poroelastic analysis of amplitude-versus-frequency variations[J].Geophysics,2009,74(6):41-48.

[7] JOHNSTON D H.The attenuation of seismic waves in dry and saturated rocks[D].Cambridge:Massachusetts Institute of Technology,2005.

[8] SUN L.Attenuation and velocity dispersion in the exploration seismic frequency band[D].Toronto:University of Toronto,2009.

[9] 高剛,李玉海,桂志先,等.基于廣義S變換頻散AVO屬性提取方法研究[J].巖性油氣藏,2015,27(4):84-90. Gao Gang,Li Yuhai,Gui Zhixian,et al.Abstraction of frequency-dependent AVO attributes based on generalized S transform[J].Lithologic Reservoirs,2015,27(4):84-90.

[10] 吳小羊.基于頻散分析技術(shù)的頻散AVO反演研究[D].武漢：中國地質(zhì)大學,2010. Wu Xiaoyang.Frequency dependent AVO inversion using spectral decomposition techniques[D].Wuhan:China University of Geosciences,2010.

[11] 張震,印興耀,郝前勇.AVO反演的頻變流體識別方法[J].地球物理學報,2014,57(12):4171-4184. Zhang Zhen,Yin Xingyao,Hao Qianyong.Frequency-dependent fluid identification method based on AVO inversion[J].Chinese Journal of Geophysics,2014,57(12):4171-4184.

[12] SHUEY R T.A simplification of the Zoeppritz equations[J].Geophysics,1985,50(4)：609-614.

[13] OSTRANDER W J.Plane-wave reflection coefficients for gas sands at non-normal angles of incidence[J].Geophysics,1984,49(10):1637-1648.

[14] 殷八斤,曾灝,楊在延.AVO 技術(shù)的理論與實踐[M].北京:石油工業(yè)出版社,1995:218-224. Yin Bajin,Zeng Hao,Yang Zaiyan.Theory and practice of AVO[M].Beijing:Petroleum Industry Press,1995:218-224.

[15] 陳兆明.疊前AVO特征分析方法在惠州地區(qū)油檢測中的應(yīng)用[J].中國海上油氣,2015,27(5):43-48.DOI:10.11935/j.issn.1673-1506.2015.05.007. Chen Zhaoming．Application of pre-stack AVO characteristics analysis on oil detection in Huizhou area[J].China Offshore Oil and Gas,2015,27(5):43-48.DOI:10.11935/j.issn.1673-1506.2015.05.007.

[16] SMITH G,GIDLOW P.Weighted stacking for rock property estimation and detection of gas[J].Geophysical Prospecting,1987,35(9):993-1014.

[17] WILSON A.Theory and methods of frequency-dependent AVO inversion[M].England:University of Edinburgh,2010.

[18] 趙萬金,楊午陽,張巧鳳,等.一種頻率域AVO油氣檢測方法[J].石油地球物理勘探,2012,47(3):436-441. Zhao Wanjin,Yang Wuyang,Zhang Qiaofeng,et al.A frequency AVO method for hydrocarbon detection[J].Oil Geophysical Prospecting,2012,47(3):436-441.

[19] 李緒宣,王建花,張金淼,等.南海深水區(qū)地震資料采集設(shè)計和處理關(guān)鍵技術(shù)及其野外試驗效果[J].中國海上油氣,2013,25(6):8-14． Li Xuxuan,Wang Jianhua,Zhang Jinmiao,et al.Some seismic acquisition designs and key processing techniques and their testing effects in the deep water areas,South China Sea[J].China Offshore Oil and Gas,2013,25(6):8-14．

[20] MULLER T M,LAMBERT G,GUREVICH B.Dynamic permeability of porous rocks and its seismic signatures[J].Geophysics,2007,72(5):149-158.

[21] 郝前勇,張世鑫,張峰,等.基于頻變AVO反演的頻散屬性估算方法及其應(yīng)用[J].石油地球物理勘探,2013,48(2):255-261. Hao Qianyong,Zhang Shixin,Zhang Feng,et al.Dispersion attributes estimation based on frequency dependent AVO inversion and its application in hydrocarbon detection[J].Oil Geophysical Prospecting,2013,48(2):255-261.

[22] 莊祖垠,陳繼宗,王征,等.深水地震資料特性及相關(guān)處理技術(shù)探析[J].中國海上油氣,2011,23(1):26-32． Zhuang Zuyin,Chen Jizong,Wang Zheng,et al.A discussion on characteristics of deep-water seismic data and their processing techniques[J].China Offshore Oil and Gas,2011,23(1):26-32．

[23] 周竹生,楊鑫.基于譜分解技術(shù)的分頻AVO反演[J].中南大學學報(自然科學版),2014,45(4):1197-1205. Zhou Zhusheng,Yang Xin.Frequency-dependent AVO inversion using spectral decomposition technique[J].Journal of Central South University (Science and Technology),2014,45(4):1197-1205.

[24] 孫萬元,張會星，孫楊.地震波衰減和頻散屬性的提取及其在油氣檢測中的應(yīng)用[J].中國海洋大學學報,2013,43(10):82-87. Sun Wanyuan,Zhang Huixing,Sun Yang.Apply S-transform to extract the attenuation and dispersion attributes of the seismic wave to detect oil and gas[J].Periodical of Ocean University of China,2013,43(10):82-87.

[25] 孫雷鳴,萬歡,陳輝,等.基于廣義S變換地震高分辨率處理方法的改進及在流花11-1油田的應(yīng)用[J].中國海上油氣,2011,23(4):234-237． Sun Leiming,Wan Huan,Chen Hui,et al.An improved method of seismic high-resolution processing based on generalized S transform and its application in LH11-1 oilfield[J].China Offshore Oil and Gas,2011,23(4):234-237．

(編輯：馮 娜)

A method of frequency dispersion AVO analysis and its application in deep water area of South China Sea

Zhang Yiming Niu Cong Han Li Ye Yunfei

(CNOOCResearchInstitute,Beijing100028,China)

Conventional technologies such as seismic inversion can generate multiple solutions for hydrocarbon detection because of few drilled wells in deep water area and lack of logging constrained. According to the seismic wave attenuation and velocity dispersion, the frequency-dependent AVO equation is used to calculate the dispersive properties of seismic wave based on seismic spectrum decomposition. By designing theoretical model, the high sensitivity of dispersion properties to gas detection is verified, and hydrocarbon indicator factor (DvP×DG) is constructed by combining frequency dispersion analysis with AVO attributes. The method is applied to forecast the target formation gas-bearing possibility of Y21-1 structure area in deep water area of China South Sea before drilling, and the predicted results are in good correspondence with actual drilling results, which imply that the frequency-dispersion AVO can increase the sensitivity of the gas reservoir, reduce the effects of water background, and can be used to detect gas as a hydrocarbon detection indicator.

seismic wave attenuation; spectrum decomposition; AVO equation; frequency dispersion attribution; hydrocarbon identification; deep water area of South China Sea

1673-1506(2016)06-0021-07

10.11935/j.issn.1673-1506.2016.06.004

張益明，男，高級工程師，2011年畢業(yè)于中國科學院地質(zhì)與地球物理研究所油儲地球物理專業(yè)，獲博士學位，主要從事地震反演和儲層預(yù)測方面的研究。地址：北京市朝陽區(qū)太陽宮南街6號院A座(郵編：100028)。E-mail:zhangym1@cnooc.com.cn。

P631.4

A

2016-03-28 改回日期：2016-04-29

*“十二五”國家科技重大專項“少井或無井條件下小波域儲層流體識別方法研究(編號:2011ZX05025-001-07)”部分研究成果。