文昌A凹陷油氣勘探中的頻率域信息應(yīng)用與探索

劉仕友， 李 輝， 楊振建， 閆安菊

(中海石油(中國)有限公司 湛江分公司， 湛江 524057)

?

文昌A凹陷油氣勘探中的頻率域信息應(yīng)用與探索

劉仕友， 李 輝， 楊振建， 閆安菊

(中海石油(中國)有限公司 湛江分公司， 湛江 524057)

地震波能量衰減與巖性、孔隙流體的性質(zhì)有一定的關(guān)聯(lián)，國內(nèi)、外學(xué)者利用高頻成分的頻譜分析對儲層含氣性檢測進(jìn)行了大量的研究：振幅高頻部分的高吸收與巖石孔隙含氣有關(guān)，低頻部分的高吸收則可能預(yù)示著巖石孔隙的含水(油)飽和度。這里總結(jié)了國內(nèi)、外利用瞬時子波吸收衰減技術(shù)的研究成果，闡述了利用頻率域信息進(jìn)行烴類檢測的基本原理。文昌氣田的水層、氣層在地震剖面表現(xiàn)為強(qiáng)振幅的亮點(diǎn)特征，利用頻率域信息對文昌XX構(gòu)造進(jìn)行了目標(biāo)巖性預(yù)測和含氣性預(yù)測，預(yù)測效果經(jīng)鉆井得以驗(yàn)證，取得了較好的效果。

吸收衰減; 瞬時子波; 預(yù)測; 類檢測

0 引言

多年勘探研究和鉆探實(shí)踐表明，南海西部文昌A凹陷生烴潛力較大，具有斷陷湖階段沉積的文昌組、恩平組兩套良好的生油氣源巖，珠海組為勘探的有利目標(biāo)層。但是珠海組地層埋深較大，一般深度為 3 300 m～4 000 m，砂巖成巖性強(qiáng)，孔隙度小，某些砂巖含有灰質(zhì)、有效儲層的波阻抗特征與泥巖疊置，好的儲層一般情況在地震上表現(xiàn)為強(qiáng)振幅特征。水層、氣層、油層的地震響應(yīng)特征(振幅域)基本類似，本區(qū)域利用振幅域技術(shù)(地震反演、AVO、亮點(diǎn)技術(shù)和彈性阻抗反演)進(jìn)行油氣藏評價效果有限，烴類檢測是制約本區(qū)域油氣勘探的一個關(guān)鍵瓶頸。

國內(nèi)、外針對頻率域信息的應(yīng)用研究已經(jīng)做了大量的工作，Mitchell等人[1]提出了一種計算地震信號能量衰減的分析方法(EAA技術(shù))，可以求取信號譜的高頻指數(shù)衰減系數(shù)；Eugene Lichman[2]提出瞬時子波上高頻部分的吸收異常往往預(yù)示著氣藏的存在，而低頻部分的吸收異常則常與巖石孔隙中的流體(油、水)相關(guān)；陳習(xí)峰等[3]對基于瞬時子波吸收分析的油氣檢測方法做了系統(tǒng)性的闡述；王宏語、張益明等[4-5]分別利用瞬時自來波吸收分析技術(shù)對陸上油氣田、海上氣田的應(yīng)用效果、適用范圍進(jìn)行了針對性討論。此外，目前國內(nèi)、外已經(jīng)發(fā)展了多種衰減解釋理論，主要包括以下幾種：①把大地當(dāng)成粘彈性介質(zhì),通過應(yīng)用Kelvin-Voight，Maxwell，SLS等模型模擬大地性質(zhì)來研究大地吸收特性；②把大地當(dāng)作不均勻介質(zhì)來研究，用散射理論來解釋地震波能量的衰減；③把大地當(dāng)作雙相介質(zhì)。相關(guān)的理論包括BIOT理論和噴射流動理論等。

經(jīng)過對文昌A凹的已鉆遇的氣層、水層、泥巖的頻譜特征對比分析發(fā)現(xiàn)，三者的頻率域信息差異明顯，利用頻率域信息是油氣識別的一個有效手段。

1 相關(guān)原理及研究思路

1.1 烴類檢測基本原理

地震波在地下傳播時，振幅強(qiáng)度(A，無量綱)可表示為式(1)。

A=A0e-afr

(1)

式中:A0為初始振幅，無量綱;r為傳播距離，m ;a為地層有效吸收系數(shù)，無量綱;f為地震波頻率Hz。通常地層含有油氣時，地層吸收系數(shù)a會增大。

由公式(1)可以看到，振幅能量同頻率成反比，即如果除去地震波傳播距離的原因，當(dāng)?shù)貙雍杏蜌鈺r，地震波高頻能量的散失高于低頻段能量，呈現(xiàn)“高頻衰減”現(xiàn)象。同時大量實(shí)踐證明，地震波在含油氣介質(zhì)中傳播時，特別是在含氣儲層中，會發(fā)生衰減現(xiàn)象，并且在頻率域上表現(xiàn)為“低頻共振、高頻衰減”的動力學(xué)特征。

對時頻方法和薄互層地震響應(yīng)關(guān)系做過深入研究的學(xué)者們[6-8]，認(rèn)為巖性調(diào)諧(巖性共振)和厚度調(diào)諧并非涇渭分明，兩者常常是共存的。Meckel[9]、Mahrai[10]曾做過地球物理模型實(shí)驗(yàn)和理論計算，認(rèn)為頻率-振幅-巖石粒度之間確有相關(guān)性，即巖性共振的確是存在的。張會星等[11]對頻率-地層粒度關(guān)系研究認(rèn)識認(rèn)為，當(dāng)頻率位于1Hz～0.8fo(fo為地震資料的主頻)中、低頻范圍時，地層或?qū)咏M較厚時，此時的頻率調(diào)諧以巖性調(diào)諧為主，即調(diào)諧頻率越低，反映了地層的粒度越粗。

1.2 地震能量衰減基本理論

衰減屬性分析的主要目的是通過屬性標(biāo)定將定量的地震衰減屬性轉(zhuǎn)化為儲層特征，地震屬性標(biāo)定中最重要的是認(rèn)識和識別能夠反映地質(zhì)意義和物理意義的具有穩(wěn)定統(tǒng)計特征的屬性。

圖1 指定能量比所對應(yīng)的頻率適應(yīng)圖

理論研究表明，與致密的單相地質(zhì)體相比，如圖1所示，當(dāng)?shù)刭|(zhì)體中含流體(油、氣或水)時，會引起地震波能量的衰減;斷層、裂縫等的存在也會引起地震波的散射,造成地震能量的衰減。因此衰減屬性是指示地震波傳播過程中的衰減快慢的物理量，是一個相對的概念。衰減屬性的分析，可以反過來指示這些衰減因素存在的可能性和分布范圍。這里的衰減屬性分析就是要通過計算出的反映地震波衰減快慢的屬性體，來指示油氣存在的可能性和分布范圍。

1.3 研究思路

文昌A凹目前已鉆探井近三十口，鉆遇地層(氣層、水層、泥巖)樣本豐富，大量的樣本統(tǒng)計顯示地震規(guī)律特征有以下特點(diǎn)：①有利儲層在地震上表現(xiàn)為“亮點(diǎn)”，巖性調(diào)諧頻率為9 Hz ～15 Hz，泥巖的巖性調(diào)諧頻率為12 Hz ～18 Hz；②氣層在地震表現(xiàn)為明顯的低頻特征，水層頻率明顯要高于水層；③氣層的高頻段衰減因子遠(yuǎn)遠(yuǎn)高于地層趨勢背景，而且也略高于水層。

鑒于本區(qū)域振幅域技術(shù)的局限性，通過總結(jié)已鉆井規(guī)律認(rèn)識發(fā)現(xiàn)：①頻率域信息與儲層、流體性質(zhì)有著較好的密切關(guān)系，建立了基于頻率域信息的文昌A凹的儲層預(yù)測與烴類檢測技術(shù)序列：②利用巖性調(diào)諧頻率進(jìn)行儲層物性預(yù)測，判斷有利儲層；③在有利儲層的基礎(chǔ)上，利用頻率吸收衰減異常進(jìn)行孔隙流體預(yù)測；④最終根據(jù)有利儲層與孔隙流體預(yù)測結(jié)果優(yōu)選鉆探井位。

2 實(shí)例應(yīng)用

文昌XX構(gòu)造位于文昌A凹陷南部邊緣、珠三南斷裂下降盤，是一個發(fā)育比較完整的背斜構(gòu)造，主要目的層珠海組一段、二段存在明顯的地震屬性異常(圖2)，成藏條件優(yōu)越。為了更好描述，分別把珠海組一段目標(biāo)稱為A目標(biāo)，二段目標(biāo)稱為B目標(biāo)，它們在疊后地震剖面上都表現(xiàn)為強(qiáng)振幅異常。

圖2 疊后純波地震剖面

2.1 巖性預(yù)測

圖3為過目標(biāo)的單線時頻分析結(jié)果，地震資料主頻約為35 Hz～40 Hz，目標(biāo)A的巖性調(diào)諧頻率約為10 Hz～14 Hz，目標(biāo)B的巖性調(diào)諧頻率為 8 Hz～13 Hz，巖性調(diào)諧頻率偏低，根據(jù)頻率-振幅-巖石粒度之間關(guān)系，可以預(yù)測目標(biāo)A和目標(biāo)B都為有利儲層，而且目標(biāo)B儲層的粒度更粗，其儲層物性可能更好。

圖3 單點(diǎn)時頻分析剖面

2.2 孔隙流體性質(zhì)預(yù)測

根據(jù)瞬時子波吸收分析與孔隙流體之間的關(guān)系，分析高頻信息衰減的變化規(guī)律(圖4)，分別選取了多點(diǎn)進(jìn)行瞬時子波特征分析，地震波的瞬時子波特征發(fā)生了規(guī)律性變化，其中目標(biāo)A所包含的a、b點(diǎn)衰減規(guī)律跟地層衰減趨勢基本一致不存在明顯的衰減異常；目標(biāo)B包含的e點(diǎn)高頻衰減異常明顯，遠(yuǎn)遠(yuǎn)脫離了地層衰減趨勢，這種特性與孔隙含氣的響應(yīng)特征類似。

圖5為過目標(biāo)線的高頻衰減因子計算剖面，依據(jù)1.2所介紹的地震能量衰減基本理論，從圖5中可以明顯看出，目標(biāo)A存在微弱衰減異常，但是特征并不明顯；目標(biāo)B高頻衰減異常強(qiáng)烈。綜合分析，預(yù)測目標(biāo)A是一套水層，目標(biāo)B為含氣儲層概率極高。為了進(jìn)一步分析目標(biāo)B的空間分布，分別在地震體和衰減體上提取了目標(biāo)B的屬性(最大值)，如圖6和圖7所示，對比分析發(fā)現(xiàn)，目標(biāo)體B的異常范圍有所差別，衰減異常范圍明顯小于地震異常范圍。

圖4 利用瞬時子波譜計算的衰減因子

圖5 高頻衰減異常剖面

圖6 目標(biāo)B HFA屬性(最大值)平面展布圖

圖7 目標(biāo)B 地震資料振幅屬性(最大值)平面展布圖

2.3 鉆井結(jié)果

根據(jù)衰減異常范圍重新預(yù)測了含氣范圍，并選擇了新井點(diǎn)位置XX。實(shí)際鉆井與我們預(yù)測基本一致：目標(biāo)A主要是大套砂巖，測井解釋為水層；目標(biāo)B錄井氣測值很高，測井解釋了一百多米的氣層。

3 結(jié)論

1)文昌A凹有利儲層在地震上表現(xiàn)為“亮點(diǎn)”，巖性調(diào)諧頻率為9 Hz ～15 Hz，泥巖的巖性調(diào)諧頻率為12 Hz～18 Hz。

2)文昌A凹?xì)鈱颖憩F(xiàn)為高頻衰減強(qiáng)異常，HFA屬性可以有效地區(qū)分氣層和水層。

3)依據(jù)上述的研究思路和預(yù)測方法，對目標(biāo)A和目標(biāo)B的預(yù)測結(jié)果與實(shí)際鉆井結(jié)果基本一致，反映了基于頻率域信息的文昌A凹的儲層預(yù)測與烴類檢測技術(shù)序列方法是可靠的。

頻率域信息技術(shù)，是文昌A凹儲層預(yù)測和烴類檢測的有效技術(shù)技術(shù)手段。頻率域信息是油氣藏預(yù)測中的一個有效補(bǔ)充，尤其針對亮點(diǎn)型油氣藏，可以有效地降低多解性?；谒矔r子波的吸收衰減分析技術(shù)，通過分別分析主頻兩端頻譜特征變化規(guī)律，可以實(shí)現(xiàn)對孔隙流體性質(zhì)的預(yù)測，但是在地震波的傳播過程中，所有的地層對地震波頻率成分具有衰減吸收作用，孔隙流體衰減是這種衰減背景下的相對異常，不具有絕對值的意義。

[1] MITCHELL J.T.,DERZHI N.,LIDMAN E.,et al. Energy absorption analysis A case study. SEG 66-th Annual Meeting Expanded abstracts , 1996：1785-1788.

[2] EUGENE LICHMAN. Unified approach to gas and fluid detection on instantaneous seismic wavelets[C].SEG Technical Program Expanded abstracts ,2003: 1699-1702.

[3] 陳習(xí)峰,毛樹禮.基于瞬時子波吸收分析的油氣檢測方法[J].小型油氣藏,2007,12(2)：29-31. CHEN X F, MAO SH L. A method of detect oil and gas based on the instantaneous wavelet [J].Small Hydrocarbon Reservoirs, 2007,12 (2): 29 -31.(In Chinese)

[4] 王宏語,趙全國. 瞬時子波吸收分析技術(shù)在復(fù)雜地區(qū)的應(yīng)用—以吐哈盆地鄯勒地區(qū)為例[J].天然氣地球科學(xué), 2007, 18(2): 289-292. WANG H Y, LI Q G. Application of absorption analysis technique for instantaneous wavelet in complex area：a case in shanle area in Tuha basin[J]. Natural Gas Geoscience, 2007, 18(2): 289-292.(In Chinese)

[5] 張益明，李緒宣，符力耘.瞬時子波吸收分析技術(shù)的應(yīng)用――以中國近海天然氣檢測為例[J].石油物探，2009，48(4)：368-376. ZHANG Y M, LI X X, FU L Y. Application of instantaneous wavelet absorption analysis technology [J]. Geophysical Prospecting for Petroleum 2009, 48 (4): 368-376.(In Chinese)

[6] 印興耀, 張奎, 張廣智.聯(lián)合時頻分布及其屬性的應(yīng)用[J].石油地球物理勘探, 2003,38(5):522-526. YIN X Y, ZHANG K, ZHANG G ZH.Application of joint time-frequency distribution and its attribution[J].OGP, 2003,38(5):522-526.(In Chinese)

[7] 夏竹, 劉超穎, 魏文博, 等. 地震道時間域頻率屬性特征和地層層序劃分[J.石油地球物理勘探, 2005,40(5):550-560. XIA ZH,LIU CH Y,WEI W B,et al. Frequency attributes characters of seismic traces in time domain and division of stratigraphic sequence[J].OGP, 2005,40(5):550-560.(In Chinese).

[8] MECKE L D , NATH A K. Geologic considerations for stratigraphic modeling and interpretation. In:CE Paytioned, Seismic Stratigraphy-Applications to Hydrocarbon Exploration[C]. AAPG Memoir 26, 1977:417-438

[9] MAHRADI. Physical modeling studies of thin beds[D]. USA:University of Houston, 1983.

Application of seismic frequency domain information in theoil and gas exploration in sag A of Wenchang

LIU Shi-you, LI Hui, YAN Zhen-jian, YAN An-ju

(CNOOC Ltd.-Zhanjiang, Zhanjiang 524057, China)

Seismic attenuation is connected to the lithology and the nature of pore fluid. Great deals of researches have been done by scholars at home and abroad on detecting reservoir gas saturation by using the spectral analysis of high frequency component. The researches have showed that high absorption of the high-frequency amplitude is often relevant to the rock porosity gas content, and high absorption of the low-frequency amplitude may signal the water (oil) saturation of pores of the rock. This paper summarizes the results on using instantaneous wavelet attenuation technology and expounds the principle of hydrocarbon detection by using frequency domain information. The water-bearing and gas-bearing reservoirs in Wenchang gas field display have similar seismic responses with strong amplitude bright spot. This paper predicts the lithology and physical property of the XX target in Wenchang by using frequency domain information, and predictive results has been verified in the actual drilling.

absorption attenuation; instantaneous wavelet; lithologic prediction; hydrocarbon detection

2015-04-02 改回日期：2015-07-06

國家科技重大專項課題(2011ZX05023-004)

劉仕友(1982- )，男，碩士，主要從事儲層預(yù)測及烴類檢測工作研究，E-mail：liushiyou@139.com。

1001-1749(2015)04-0508-04

P 631.4

A

10.3969/j.issn.1001-1749.2015.04.15