胥良君，李錄明，黨錄瑞，鐘 峙，歐陽(yáng)煬

(1.成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，成都 610059;2.西南油氣田分公司 重慶氣礦，重慶 400021)

頻譜分解在生物礁儲(chǔ)層油氣檢測(cè)中的應(yīng)用

胥良君1，2，李錄明1，黨錄瑞2，鐘 峙1，歐陽(yáng)煬1

(1.成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，成都 610059;2.西南油氣田分公司 重慶氣礦，重慶 400021)

生物礁是特殊的沉積體，其生長(zhǎng)形態(tài)、發(fā)育規(guī)模通常因古地理位置及古氣候環(huán)境的差異而不盡相同，由于其與圍巖存在波阻抗差異，在地震反射剖面上常形成上隆、下凹、內(nèi)部雜亂反射等特征，解釋人員一般也通過(guò)該特征來(lái)尋找生物礁。但由于振幅異常往往受多種因素的影響，因此降低了生物礁氣藏勘探的成功率。這里利用短時(shí)傅里葉變換、連續(xù)小波變換對(duì)生物礁異常體進(jìn)行頻譜分解，通過(guò)頻降分析，與實(shí)鉆情況吻合，驗(yàn)證了該方法的有效性。

生物礁;頻譜分解;短時(shí)傅里葉變換;連續(xù)小波變換

0 前言

生物礁是特殊的地質(zhì)沉積體，其主要發(fā)育于低緯溫暖的淺海環(huán)境，按沉積環(huán)境、形態(tài)、與海岸關(guān)系等標(biāo)準(zhǔn)可有多種分類(lèi)。生物礁是良好的油氣儲(chǔ)集場(chǎng)所，具有孔隙度高，儲(chǔ)量豐度大，油氣產(chǎn)能高等特點(diǎn)。世界上日產(chǎn)萬(wàn)噸級(jí)的油氣井，半數(shù)為生物礁儲(chǔ)層，美國(guó)找到與巨大的馬尼拉礁有關(guān)的油田就有三十七個(gè)。加拿大西部找到與生物礁有關(guān)的油氣田總數(shù)達(dá)到三百多個(gè)，占加拿大油氣儲(chǔ)量的60%以上［1］。

礁體的生長(zhǎng)形態(tài)及發(fā)育規(guī)模，通常因古地理位置及古氣候環(huán)境的差異而不盡相同。因此，其非均質(zhì)性較強(qiáng)，巖性橫向變化快。由于與圍巖存在的波阻抗差異，所以在地震反射剖面上常形成上隆、下凹、減弱、內(nèi)部雜亂反射等特征［2～4］。解釋人員通常對(duì)該特征進(jìn)行分類(lèi)總結(jié)，以形成生物礁地震響應(yīng)模式，用于尋找生物礁。但振幅異常往往受多種因素的影響，且在川東地區(qū)，礁灰?guī)r若未經(jīng)白云巖化作用改造，亦不能形成良好的儲(chǔ)層，但其地震響應(yīng)特征仍表現(xiàn)為異常，因此降低了生物礁氣藏勘探的成功率。

根據(jù)地震波的傳播理論可知，地震波在地下巖層傳播時(shí)，由于巖層的非完全彈性，使地震波的能量轉(zhuǎn)化為熱能，造成振幅、頻率衰減，高頻被吸收。當(dāng)巖石中含有油氣時(shí)，衰減尤為顯著［5～8］，作者利用短時(shí)傅里葉變換及連續(xù)小波變換，對(duì)生物礁異常體進(jìn)行頻譜分解，通過(guò)頻降來(lái)判斷地質(zhì)體的含油氣性，實(shí)際資料的處理結(jié)果與實(shí)鉆情況吻合，驗(yàn)證了該方法的有效性。

1 方法原理

1.1 短時(shí)傅里葉變換

短時(shí)傅里葉變換(STFT)，是一種利用固定時(shí)窗進(jìn)行傅里葉變換的時(shí)頻分析方法［9］。它假設(shè)在時(shí)窗內(nèi)，信號(hào)是平穩(wěn)的，對(duì)時(shí)窗內(nèi)信號(hào)進(jìn)行傅氏變換后，滑動(dòng)時(shí)窗再進(jìn)行下一段信號(hào)的計(jì)算，直至完成整段信號(hào)的時(shí)頻分析。

連續(xù)信號(hào)x(t)∈L2(R)的短時(shí)傅里葉變換定義為:

由式(1)可知，h(x)窗口函數(shù)是有時(shí)頻局域性的時(shí)窗函數(shù)，使得STFT具有局部特性，h(x)函數(shù)的形狀及長(zhǎng)度，直接影響變換的分辨率。常用的時(shí)窗函數(shù)有 Rectangular窗、Bartlett窗、Kaiser窗、Gaussian窗、Hanning窗、Blackman窗等。

窗口函數(shù)滿(mǎn)足如下條件:

x(t)的反變換公式為:

短時(shí)傅里葉變換是經(jīng)典的線(xiàn)性時(shí)頻分析方法，其變換結(jié)果表征的頻譜特征比小波變換更容易被理解。但由于其時(shí)窗固定，時(shí)頻網(wǎng)格不能隨信號(hào)自適應(yīng)改變，因此具有單一的頻率分辨率。

1.2 連續(xù)小波變換

連續(xù)小波變換(CWT)是由法國(guó)科學(xué)家Morlet于1980年在進(jìn)行地震資料分析時(shí)引入，并在此基礎(chǔ)上由Morlet和物理學(xué)家Grossma等人建立了完整的連續(xù)小波變換的幾何體系［10］。與STFT比較，CWT是一種多分辨分析方法，具有平移伸縮不變性。它繼承了STFT的局部化思想，能在時(shí)間、頻率域同時(shí)對(duì)信號(hào)進(jìn)行局部化分析，并針對(duì)高頻信號(hào)、低頻信號(hào)自動(dòng)采用窄、寬時(shí)窗進(jìn)行分析。

連續(xù)信號(hào)x(t)∈L2(R)的連續(xù)小波變換定義為:

由式(4)可知，Ψa、b(t)是母小波經(jīng)過(guò)伸縮移動(dòng)后的結(jié)果，當(dāng)a、b不斷改變，便能產(chǎn)生函數(shù)族{Ψ，(t)}，人們稱(chēng)之為小波基。其中，每個(gè)Ψa，b(t)被稱(chēng)為一個(gè)小波。經(jīng)典類(lèi)小波又可以分為Haar小波、Morlet小波、Mexicanhat小波和Gaussian小波等。Ψ(t)滿(mǎn)足式(5)。

其反變換定義為:

2 實(shí)際應(yīng)用效果

四川盆地從二疊紀(jì)開(kāi)始，地殼全面下沉廣泛海侵，早二疊世末，受東吳運(yùn)動(dòng)影響，盆地露出水面遭受剝蝕，晚二疊世晚期，海侵規(guī)模擴(kuò)大，以長(zhǎng)興組為代表的淺海碳酸鹽巖廣泛分布。WL構(gòu)造位于四川盆地川東南中隆高陡構(gòu)造區(qū)，上二疊統(tǒng)沉積期位于開(kāi)江～梁平海槽西南側(cè)開(kāi)闊碳酸鹽臺(tái)地之上。長(zhǎng)興組地層巖性以灰色、褐灰色泥粉晶、細(xì)粉晶灰?guī)r為主，夾薄層硅質(zhì)灰?guī)r、泥灰?guī)r、燧石結(jié)核灰?guī)r，局部地區(qū)發(fā)育臺(tái)內(nèi)點(diǎn)礁，以溶孔白云巖、生屑灰?guī)r為主。

對(duì)WL構(gòu)造區(qū)塊內(nèi)生物礁氣井S15井、W118－1井、生物礁鉆探失利井W124井及生物礁專(zhuān)層建議井J6井過(guò)井地震剖面，進(jìn)行STFT、CWT頻譜分解形成分頻剖面，然后以井位所在CDP產(chǎn)生頻率道集。

下頁(yè)圖1(a)～圖1(d)分別為S15井、W118－1井、W124井、J6井的過(guò)井局部時(shí)間剖面，矩形框內(nèi)為地震異常響應(yīng)。不難看出，其不同程度的存在同相軸減弱、上隆、內(nèi)部反射雜亂等特征，因此，極易將其判定為生物礁引起的地震響應(yīng)。但從實(shí)鉆情況來(lái)看，S15井、W118－1井均獲氣，而W124井為干層。

后面圖2(a)～圖2(d)分別為針對(duì)S15井、W118－1井、W124井及J6井過(guò)井剖面采用 STFT方法進(jìn)行頻譜分解，并以井位所在CDP產(chǎn)生的頻率道集。從圖2上不難看出，在地質(zhì)目標(biāo)的下方，S15井、W118－1井及 J6井存在明顯頻降，而W124井則沒(méi)有此特征。

見(jiàn)后面，圖3(a)～圖3(d)分別為針對(duì)S15井、W118－1井、W124井及J6井采用 CWT方法進(jìn)行頻譜分解，并以井位所在CDP產(chǎn)生的頻率道集。從圖3上依然可以看出，在地質(zhì)目標(biāo)的下方，S15井、W118－1井及J6井存在明顯頻降，而W124井沒(méi)有此特征。并且，由于CWT方法的多分辨特性，因此比STFT方法具有更高的分辨率。

圖4(見(jiàn)后面)為利用20Hz與70Hz沿層切片求差產(chǎn)生的平面預(yù)測(cè)圖，黑色區(qū)域?yàn)樯锝赴l(fā)育有利區(qū)。從圖4中可以看到，S15井、W118－1井及J6井均位于有利區(qū)，而W124則在不利區(qū)，同時(shí)，工區(qū)左下方邊界處亦表現(xiàn)為有利區(qū)，但由于其屬于工區(qū)邊界信噪比相對(duì)較低，因此未在該區(qū)定井。

圖1 過(guò)井地震剖面Fig． 1 Cross － well seismic profile

3 結(jié)論

通過(guò)對(duì)過(guò)井剖面進(jìn)行頻譜分解分析可以看出，STFT、CWT方法均可以反映儲(chǔ)層，由于含油氣后對(duì)高頻成份的吸收作用，驗(yàn)證了該方法用于儲(chǔ)層含油氣性檢測(cè)的可行性。同時(shí)，由于CWT方法具有多分辨特性，因此比STFT方法具有更高的分辨率，而STFT則具有較好的宏觀性。在進(jìn)行實(shí)際計(jì)算時(shí)，對(duì)頻率分析范圍、步長(zhǎng)、小波類(lèi)型等參數(shù)的選擇，都會(huì)對(duì)頻譜特征產(chǎn)生影響。因此需進(jìn)行參數(shù)試驗(yàn)，以便合理選取。在對(duì)待預(yù)測(cè)結(jié)果時(shí)，應(yīng)該采取綜合分析的原則，以提高鉆探成功率，降低風(fēng)險(xiǎn)。

作者在本文運(yùn)用該方法對(duì)建議J6井進(jìn)行了分析，從結(jié)果來(lái)看，其與同區(qū)其它獲氣井具有相似的頻降特征，并且在預(yù)測(cè)平面圖上位于有利區(qū)，為該井實(shí)施提供了進(jìn)一步有利依據(jù)。

［1］ 衛(wèi)平生，劉全新，張景廉，等.再論生物礁與大油氣田的關(guān)系［J］.石油學(xué)報(bào)，2006，27(2):38

［2］ 劉殊，郭旭升，馬宗晉，等.礁灘相地震響應(yīng)特征和油氣勘探遠(yuǎn)景［J］.石油物探，2006，45(5):452.

［3］ 汪澤成，趙文智，張林，等.四川盆地構(gòu)造層序與天然氣勘探［M］.北京:地質(zhì)出版社，2002.

［4］ 王興志，張帆，馬青，等.四川盆地東部晚二疊世－早三疊世飛仙關(guān)期礁、灘特征與海平面變化［J］.沉積學(xué)報(bào)，2002，20(2):249.

［5］ 董寧，楊立強(qiáng).基于小波變換的吸收衰減技術(shù)在塔河油田儲(chǔ)層預(yù)測(cè)中的應(yīng)用研究［J］.地球物理學(xué)進(jìn)展，2008，23(2):533.

［6］ CASTAGNA J P，SUN S J，SIEGFRIED R W． Instantaneousspectral analysis: detection of low － frequencyshadows associated with hydrocarbons［J］． The LeadingEdge，2003，22( 2) : 120．

［7］ 李子順.地震波衰減規(guī)律及其恢復(fù)方法［J］.地球物理學(xué)進(jìn)展，2007，22(5):1545.

［8］ 馬昭君，劉洋.地震波衰減反演研究綜述［J］.地球物理學(xué)進(jìn)展，2005，20(4):1074.

［9］ JIANG Y Q，HEY G． New digital algorithm for measuringelectric frequency based on adaptive short time fouriertransform［J］． journal of Electronic in Measurementand Instrument，2006，20( 2) : 10．

［10］李宏兵，趙文智，曹宏，等.小波尺度域含氣儲(chǔ)層地震波衰減特征［J］.地球物理學(xué)報(bào)，2004，47(5):892.

圖4 生物礁預(yù)測(cè)平面圖Fig． 4 Prediction map of reef

TE121.3+4

A

1001—1749(2011)05—0544—05

博士學(xué)科專(zhuān)項(xiàng)基金(20070616004)

2011－02－28 改回日期:2011－06－23

胥良君(1981－)，男，博士，主要從事地質(zhì)勘探及地球物理信號(hào)與信息處理研究。