畢麗飛(中國石油大學(北京)地球物理與信息工程學院,北京102249) (中石化勝利油田分公司技術發(fā)展處,山東東營257000)

基于時頻分析疊前地層Q值提取方法研究

畢麗飛(中國石油大學(北京)地球物理與信息工程學院,北京102249) (中石化勝利油田分公司技術發(fā)展處,山東東營257000)

地震波衰減屬性是含氣儲層預測的重要參數(shù)。研究提出了基于時頻分析疊前地層Q值提取方法,該方法基于疊前地震資料中豐富的儲層信息,利用S變換時頻譜具有較高時頻分辨率的特征,提高了地層Q值的估算精度。在模型計算和勝利油田墾71工區(qū)的實際應用表明,該方法能夠較好地預測氣層空間展布,在含氣區(qū)Q值較低,異?，F(xiàn)象明顯,在非含氣區(qū)Q值較高,異?，F(xiàn)象不明顯,尤其是強振幅位置。

疊前道集;地震波衰減;零相位子波;儲層預測

在進行儲層預測,尤其是天然氣藏的預測過程中,地層的吸收特征是一個非常重要的參數(shù),它對巖性的變換具有較高的靈敏度,尤其是含氣時表現(xiàn)出很大的異常值[1,2]。因此,利用地層的吸收特征并結合地震波的速度、波速比等地震屬性參數(shù)和測井信息,可進一步提高儲層預測和氣藏勘探的精度[3]。

地層吸收特征的估算方法有很多,主要有時間域、頻率域、時頻域和基于反演理論的方法[4～10]。目前,基于時頻域和基于反演理論的方法成為研究的熱點。該次研究筆者主要利用S變換時頻分析方法,根據(jù)地震波在一維黏彈性介質中的傳播方程,在時頻域中推導了地層Q值的估算公式,并通過炮檢距歸零處理將其應用到疊前地震資料中。該方法在模型計算和勝利油田墾71工區(qū)中的應用都取得了較好的儲層含氣性預測效果。

1 時頻分析Q值提取方法原理

S變換是由Stockwell等[10]于1996年提出的一種較新的時頻分析方法,它與小波變換、短時窗傅里葉變換都有一定的關系,但又比這兩者具有更多的優(yōu)點。S變換有更高的時頻分辨率,在地層Q值提取中可以充分考慮因地層吸收而引起的振幅、頻率等信息的變化。

S變換公式為:

式中:f表示頻率,s-1;ξ表示時窗函數(shù)的中心點,控制高斯窗函數(shù)在時間軸上的位置,s;t表示傳播時間,s;h(t)表示信號。

式中:Q是品質因子,1。

對式(2)進行S變換,可以得到頻率域表達式[16]:

對式(3)進行逆傅里葉變換,可以得到時間域表達式:

從而得到振幅譜As(t,f)的表達式:

選擇目的層,假定地震波傳播到目的層上、下界面的旅行時分別是t0和t1。對目的層上、下界面的振幅譜作比并取對數(shù),可得到地層Q值的計算公式:

從式 (6)中可以看出,當假定地震子波為一般零相位子波時,對于已選定的目的層,頻譜比對數(shù)與a呈線性關系,進行線性回歸,可得到地層Q值。

上述方法是基于疊后地震資料提出的。在疊前地震資料中,計算的頻譜比斜率是隨著炮檢距的變化而變化的,所以在計算地層Q值之前,需要進行炮檢距歸零處理。假設傳播時間的變化滿足經典的正常時差方程,即:

式中:Δt為正常時差,s;x是炮檢距,m;v是速度,m/s。

由式(5)可知,零炮檢距處的頻譜比斜率為:

式中:p0是炮檢距為零時計算的頻譜比斜率,1;Δt1是地震波在零炮檢距處通過目的層的上、下界面所用的時間差,s;Q0是炮檢距為零時計算的品質因子,1。

因此,在疊前地震資料中,水平方向上的頻譜比斜率近似為:

式中:px是炮檢距為x時計算的頻譜比斜率,1;Qx是炮檢距為x時計算的品質因子,1。將Δt、Δt1代入式(9)可得:

假設在小排列情況下,Qx≈Q0,式(10)近似為:

由式(11)可以看出,在同一位置處,不同炮檢距的px與x2呈線性關系,利用上述關系進行線性回歸計算截距,即為零炮檢距處的斜率值p0,進而根據(jù)式 (5)可求得地層的Q值。在計算的過程中,該方法容易受動校正拉伸的影響,當動校正拉伸比較嚴重時,炮檢距歸零處理中可以考慮適當剔除一部分資料。在具體處理中,時窗的選擇也非常重要,不同的數(shù)據(jù)應選擇不同的時窗,時窗的選擇一般為目的層段時窗的2～3倍。

2 模型計算與實際工區(qū)應用分析

2.1 模型計算

表1是一個簡單的3層模型參數(shù)表,第2層是目的層,地層Q值比較低。圖1是通過表1的參數(shù)利用單程波正演模擬得到的疊前CMP(共中心點)道集。圖2是圖1動校正之后結果。

表1 模型參數(shù)表

圖1 疊前地震記錄

圖2 動校正后的疊前地震記錄

圖3是利用圖2的數(shù)據(jù)選擇第2層為目的層,得到的頻譜比對數(shù)與a的變化關系?？梢钥闯?頻譜比對數(shù)隨著a的變化呈線性關系,通過線性回歸可得到頻譜比斜率。

圖4是目的層頻譜比斜率隨偏移距平方的變化,可以看出,頻譜比斜率隨著偏移距平方的變化呈線性關系,通過線性回歸求取截距,即為零偏移距的頻譜比斜率,通過計算可得到地層Q值。

圖5是地層Q值的求取結果,可以看出,在目的層的位置,Q值比較低,為29.7;在非目的層的位置,Q值比較高,為497.1,這與表1中的模型參數(shù)值基本一致,說明該方法的正確性。

圖3 頻譜比對數(shù)與a的關系圖

圖4 頻譜比斜率與偏移距平方的關系圖

圖5 地層Q值曲線

2.2 實際工區(qū)應用

勝利油田墾71工區(qū)位于沾化凹陷中部,是墾西油田的主力開發(fā)區(qū)。其主要氣層為新近系館陶組二段3+4亞段(Ng32+4),氣層主要發(fā)育在工區(qū)的中部,含氣范圍呈近東西向條帶狀展布。對該區(qū)進行地震波衰減研究,有助于該區(qū)的含氣儲層預測。針對墾71工區(qū)過井的2條線提取地層Q值,進行含氣儲層的預測。

圖6 疊后地震剖面A

圖6是墾71工區(qū)的某一疊后地震剖面A,k71-36井經過該條線,氣層的位置大約在時間1.15s處,道號是86。圖7是利用本文方法提取的地層Q值剖面,可以看出,在氣藏的位置 (圖中方框處)Q值比較低,異?，F(xiàn)象比較明顯;而在非含氣的位置,異常現(xiàn)象得到了一定程度的壓制。

圖8是墾71工區(qū)的某一疊后地震剖面B,k71井經過該條線,氣層的位置大約在時間1.15s處,道號是126。圖9是利用本文方法提取的地層Q值剖面,可以看出,在氣藏的位置(圖中方框處)Q值比較低,異?，F(xiàn)象比較明顯,表現(xiàn)出強烈的吸收特征。

實際應用結果表明,該方法能夠準確地預測氣藏的位置,為研究區(qū)的天然氣藏的預測提供了有力的技術支持。

3 結論和認識

基于時頻分析的疊前地層Q值提取方法一方面利用了S變換時頻分析地層吸收引起的地震波振幅、頻率等信息的變化,另一方面利用了疊前地震資料信息豐富的特點,提高了地層Q值提取的精確度。在模型計算和勝利油田墾71工區(qū)的實際應用表明,該方法能夠較好地預測氣層空間展布,在含氣區(qū)Q值較低,異?，F(xiàn)象明顯;在非含氣區(qū)Q值較高,異?，F(xiàn)象不明顯,尤其是強振幅位置,取得了較好的應用效果。

圖7 基于疊前地震數(shù)據(jù)提取的地層Q值剖面A

圖8 疊后地震剖面B

圖9 基于疊前地震數(shù)據(jù)提取的地層Q值剖面B

[1]劉祖沅.巖石中地震波衰減的實驗研究[J].地震科學研究,1983,(6):34～40.

[2]Quan Youli,Harris J M.Seismic attenuation tomography using the frequency shift method[J].Geophysics,1997,63(3):895～905.

[3]Zhang C J,Ulrychz T J.Estimation of quality factors from CMP record[J].Geophysics,2002, 67(5):1542～1547.

[4]李宏兵,趙文智,曹宏,等.小波尺度域含氣儲層地震波衰減特征[J].地球物理學報,2004,47 (5):892～898.

[5]Sun Shengjie.Determination of seismic attenuation from surface and downhole measurements[D]. Oklahoma:University of Oklahoma,2000.

[6]Dasgupta R,Clarkz R A.Estimation of Q from surface seismic reflection data[J].Geophysics,1998,63(12):2120～2128.

[7]Hackert C L,Parra J O.Improving Q estimation from seismic reflection data using well-log-based localized spectral correction[J]. Geophysics,2004,69(6):1521～1529.

[8]Gao Jinghuai,Yang Senlin.Quality factors estimation using wavelet’s envelope peak instantaneous frequency[J].Geophysics, 2009,74:2457～2461.

[9]Stockwell R G,Mansinha L,Lowe R P.Localization of the complex spectrum:The S transform[J].IEEE Transactions on Signal Processing,1996,44(4):998～1001.

[10]鄒文.S-變換時頻分析技術及其在地震勘探中的應用研究[D].北京:中國地質大學,2005.

[編輯] 龔丹

P631.44

A

1000-9752(2014)12-0070-04

2014-11-25

畢麗飛(1976-),男,1999年大慶石油學院畢業(yè),高級工程師,博士生,現(xiàn)主要從事地震資料處理方法研究及科研項目管理等工作。