張繁昌, 張汛汛, 翁 斌, 謝秀川, 鄭 煬

(1.中國石油大學地球科學與技術學院,山東青島 266580; 2.中海油研究總院,北京 100028;3.華北油田,河北任丘 062552; 4.中海石油(中國)有限公司天津分公司,天津 300452)

基于匹配追蹤的疊前道集吸收補償和頻散校正

張繁昌1, 張汛汛1, 翁 斌2, 謝秀川3, 鄭 煬4

(1.中國石油大學地球科學與技術學院,山東青島 266580; 2.中海油研究總院,北京 100028;3.華北油田,河北任丘 062552; 4.中海石油(中國)有限公司天津分公司,天津 300452)

基于匹配追蹤算法,對地震道集吸收和頻散進行補償。首先,運用匹配追蹤算法對疊前道集作匹配分解,以零偏移距匹配子波建立隨子波中心時間和中心頻率乘積值變化的振幅衰減曲線,進而得到不同匹配子波衰減振幅的增益系數(shù);然后,基于射線路徑獲取零偏移距吸收衰減與非零偏移距吸收衰減間的關系,運用零偏移距振幅增益系數(shù)與道集角度計算各個角度地震道的匹配子波振幅增益系數(shù);最后,用振幅增益系數(shù)對不同匹配子波衰減振幅分別進行吸收補償,并重構地震信號;對于頻散影響,通過匹配子波求取品質因子進行頻散校正。數(shù)據(jù)資料測試發(fā)現(xiàn),深部吸收衰減能量得到有效補償,同時消除了吸收對振幅隨偏移距變化特征的影響。

地震勘探; 疊前道集; 吸收衰減; 匹配追蹤; 振幅增益系數(shù); 頻散校正

地震波在地下地層中傳播時,由于介質的黏彈性性質,地震波會出現(xiàn)明顯的吸收衰減現(xiàn)象,表現(xiàn)為振幅衰減和波形畸變,降低地震資料的分辨率。自1962年Futterman首次提出吸收衰減理論以及補償方法[1]以來,在此基礎上發(fā)展了許多補償算法,主要有頻率域[2-3]、時間域[4]和時頻域[5]的Q值求取方法,但這些方法求取的Q值精度不佳,為提高疊前疊后衰減補償?shù)木?非線性最優(yōu)化反演[6]、正則化反演[7-9]、反Q濾波[10-13]、Gabor域振幅平滑[14]等方法被發(fā)展起來。對疊前地震資料來說,吸收衰減使振幅的橫向反射特征發(fā)生改變,很容易影響AVA的正確分析,針對這一問題,許多基于NMO拉伸校正[15]、常Q模型[16]、射線路徑[17-20]的補償方法發(fā)展了起來。雖然這些方法都能對疊前衰減進行補償,但是都有一個共同的問題,即在補償過程中高頻成分振幅會被過度放大,產(chǎn)生高頻噪音,降低地震資料信噪比。匹配追蹤算法[21]是一種高精度的地震信號分解與重構算法,但計算效率一直制約著該算法的應用,為提高算法計算效率,對其時頻原子庫進行了改進[22-24],得到了較好的效果。筆者在前人的研究基礎上,結合匹配追蹤算法構建隨匹配子波中心時間和中心頻率乘積值變化的振幅衰減曲線,對疊前道集進行振幅補償和頻率校正,將二維信息時間和頻率降低為一維信息時間和頻率的乘積值,以增強算法的穩(wěn)健性。

1 方法原理

1.1 吸收衰減機制

在黏彈性介質中,波數(shù)是一個滿足因果關系的與頻率相關的復函數(shù)。Wang[25]通過研究發(fā)現(xiàn),為增加頻散校正的準確度,復波數(shù)可表示為

(1)

式中,γ和v分別表示與頻率有關的衰減因子和相速度;f為頻率。在Kolsky[26]衰減模型中,假設衰減是一個關于頻率的線性函數(shù),衰減因子表示為

(2)

相速度表示為

(3)

式中,Q為品質因子;f0為參考頻率;v0為參考頻率處的參考相速度。根據(jù)Bickel和Natarajan[27]的理論,公式(3)可以簡寫為

(4)

對公式(3)作一階級數(shù)近似可得到慢度P(f):

(5)

將公式(2)、(5)代入式(1)中,可得復波數(shù)為

(6)

在頻率域內(nèi)求解單程波動方程的平面波,可得衰減信號X(f)和原信號X0(f)有如下關系:

X(f)=X0(f)exp[ik(f)ξ] .

(7)

(8)

式中,t為地震波傳播時間。定義衰減算子如下:

(9)

式中,實部為振幅衰減項,虛部為速度頻散項。

1.2 零偏移距數(shù)據(jù)的吸收補償

對于零偏移距地震道,其吸收衰減只與深度有關。當只考慮振幅衰減時,零偏移距衰減算子為

(10)

式中,t0為零偏移距反射時間。由式(10)不難發(fā)現(xiàn),衰減算子是關于時間和頻率的函數(shù),且與兩者的乘積有關,因此可以尋找衰減振幅與時間和頻率乘積值間的關系進行吸收補償。

匹配追蹤是一種自適應信號分解算法,可以對信號進行匹配分解。對零偏移距地震信號匹配分解后,其可表示為一系列匹配子波的線性疊加:

(11)

式中,Ai、t0i、φi和f0i分別表示不同匹配子波的振幅、中心時間、相位和中心頻率。將各個匹配子波中心時間和中心頻率乘積值t0if0i和相應匹配子波振幅投影到坐標平面內(nèi);把t0if0i分為等間隔的條段并求取各個條段內(nèi)的子波振幅平均值,即得到隨t0if0i變化的振幅衰減曲線a(t0if0i);然后由振幅衰減曲線計算振幅增益曲線:

(12)

式中,α為白噪系數(shù);Amax為振幅衰減曲線的最大值。通過式(12)即可求得各個匹配子波衰減振幅的增益系數(shù),對衰減振幅進行補償:

(13)

1.3 非零偏移距數(shù)據(jù)的吸收補償

對于非零偏地震道,其振幅衰減算子可以表示為

(14)

式中,tf為地震波非零偏移距反射時間。分析可知,對于同一反射層,零偏移距反射時間t0與非零偏移距反射時間tf有如下關系:

(15)

式中,θ為地震波傳播方向與垂直方向的夾角,(°)。則式(14)可寫為

(16)

由式(16)可知:對于同一反射層,θ越大,振幅衰減越嚴重;并且非零偏移距振幅衰減算子可由零偏移距振幅衰減算子與夾角θ表示,即其衰減同時與深度和偏移距有關。

結合零偏移距吸收補償方法,首先計算與深度有關的匹配子波振幅增益系數(shù),然后結合夾角θ得到非零偏移距各個匹配子波振幅增益系數(shù),可表示為

(17)

式(17)同時補償了深度衰減和偏移距衰減。當僅補償與偏移距有關衰減,即非零偏移距與零偏移距的吸收差異時,增益系數(shù)可寫為

(18)

然后,用各個匹配子波振幅增益系數(shù)對衰減振幅進行補償,最后,用振幅補償后的各個匹配子波重構地震道。

1.4 頻散校正

對于速度頻散,通過匹配追蹤算法求取地層品質因子[28]:

(19)

式中,χ=t0if0i為匹配子波中心時間和中心頻率的乘積值;a(χ)=a(t0if0i)為隨χ變化的振幅衰減曲線;a0為未衰減振幅。這樣就可以通過擬合斜率得到地層品質因子值,然后運用公式(9)的虛部進行頻散校正。

2 方法測試與應用效果分析

2.1 理論模型測試

為測試本文吸收補償方法的有效性,建立衰減道集模型對方法進行測試。圖1為建立該衰減道集模型所需參數(shù),包括密度、橫波速度、縱波速度以及品質因子模型。根據(jù)佐普里茲方程計算得到反射系數(shù),并用60 Hz雷克子波合成疊前道集;在此過程中不考慮透射損失和幾何擴散,振幅衰減和速度頻散用式(9)給出,其中速度頻散的參考頻率為200 Hz。

圖1 模型參數(shù)Fig.1 Model parameters

圖2(a)為沒有地層吸收的疊前道集,圖2(b)為包含地層吸收的疊前道集。通過對比可以發(fā)現(xiàn),由于受地層吸收的影響,縱向上,隨著傳播時間的增加,地震波振幅衰減越嚴重,分辨率降低;橫向上,地震反射特征發(fā)生明顯變化,改變了地震記錄的同相性,不能反映真實的振幅隨入射角變化特征。

圖3(a)是對圖2(b)通過匹配追蹤求取品質因子然后進行頻散校正后的結果,由于頻散造成的同相軸下拉現(xiàn)象得到了校正,改善了疊前道集的同相性;圖3(b)是對圖3(a)的結果進一步作與偏移距相關吸收補償后的結果,由于消除了與偏移距相關的吸收差異,地震記錄的橫向一致性得到明顯的提高,AVA反射特征得到恢復;圖3(c)是圖3(b)作與深度有關吸收補償后的結果,可以看出深部衰減能量得到有效補償,并且與沒有地層吸收及速度頻散效應的圖2(a)完全一致。

圖2 疊前道集對比Fig.2 Pre-stack seismic gathers comparison

圖4是圖2(b)模型數(shù)據(jù)中480 ms處反射子波與偏移距相關吸收補償前后對比。圖4(a)是補償前后反射振幅隨偏移距的變化對比,可以看出:由于吸收的影響,補償前反射振幅隨偏移距增大而減小,改變了橫向反射特征,產(chǎn)生錯誤的AVA信息,補償后反射振幅隨偏移距增大而增大,恢復了正確的橫向反射特征。圖4(b)是補償前后反射主頻隨偏移距的變化對比,可以看出:由于頻散的影響,反射波主頻隨偏移距增大逐漸降低,削弱了疊前道集橫向一致性,補償后不同偏移距地震反射波主頻沒有發(fā)生變化,這可以增強疊前道集的聚焦性。

圖3 圖2(b)數(shù)據(jù)吸收補償?shù)慕Y果Fig.3 Results of absorption compensation for Fig.2(b)

圖4 圖2(b)中480 ms處反射子波與偏移距相關吸收補償前后對比Fig.4 Comparison of reflection wavelet at 480 ms before and after offset-related absorption compensation for Fig.2(b)

2.2 實際資料處理

圖5(a)為某氣井附近的疊前道集,由于地層吸收的影響,深部能量衰減嚴重,分辨率較低;橫向上反射振幅特征并沒有隨偏移距的增大而增大,不能代表真實的AVA特征。圖5(b)是進行頻散校正和與深度有關吸收補償后的道集,頻散引起的同相軸相位畸變問題得到解決,且深部衰減能量得到有效補償,分辨率得到提高;圖5(c)是對圖5(b)結果進一步作與偏移距有關吸收補償后的結果,從紅框中對比可以看出,補償后反射強度隨偏移距的增大而增大,顯現(xiàn)出典型的含氣儲層AVA反射特征。

圖5 疊前道集經(jīng)吸收補償和頻散校正結果Fig.5 Results of a pre-stack gather after absorption compensation and dispersion correction

圖6是與深度有關吸收衰減不同補償方法的結果對比,圖6(c)是反Q濾波方法同時進行振幅補償和頻散校正的結果?？梢钥闯?由于頻散校正是無條件穩(wěn)定的,而振幅補償是不穩(wěn)定的,所以該方法在補償深部信號的同時,高頻噪音也被相應放大,嚴重降低了補償后地震數(shù)據(jù)的信噪比;而本文方法的補償結果如圖6(b)所示,在補償衰減能量的同時,有效抑制了高頻噪音的放大,使補償后的地震數(shù)據(jù)具有更高的分辨率和信噪比。

圖6 與深度有關的吸收補償結果對比Fig.6 Comparison of different depth-related absorption compensation methods

3 結束語

基于匹配追蹤算法,提出一種疊前資料的吸收補償和頻散校正方法。首先,對疊前道集進行匹配分解,用零偏移距匹配子波求取隨匹配子波中心時間和中心頻率乘積變化的匹配子波振幅衰減曲線,以此獲取匹配子波振幅增益系數(shù),對與深度有關吸收進行補償;然后通過零偏移距地震數(shù)據(jù)匹配子波振幅增益系數(shù)和非零偏地震道集角度計算非零偏移距匹配子波振幅增益系數(shù),對與偏移距有關的吸收衰減進行補償。本文中方法通過降維處理,將二維信息時間和頻率降低為一維信息時間和頻率的乘積,增加了算法的穩(wěn)健性。對于速度頻散,同樣通過匹配子波求取地層品質因子,從而進行頻散校正。模型數(shù)據(jù)和實際資料的處理結果表明,本文中方法有效補償了深部地層的衰減能量,極大地改善了地震數(shù)據(jù)的AVA反射特征,有利于提取正確的AVA信息,為疊前反演提供更可靠的道集數(shù)據(jù)。

[1] FUTTERMAN W I. Dispersive body waves[J]. Journal of Geophysical Research, 1962,7(13):5279-5291.

[2] BATH M. Spectral analysis in geophysics[M]. New York: Elsevier, 1974.

[3] ZHANG CH J, ULRYCHZ T J. Estimation of quality factors from CMP record[J]. Geophysics, 2002,67(5):1542-1547.

[4] JANNSEN D, VOSS J, THEILEN F. Comparison of methods to determineQin shallow marine sediments from vertical reflection seismograms[J]. Geophysical Prospecting,1985,23:479-497.

[5] 王小杰,印興耀,吳國忱.基于疊前地震數(shù)據(jù)的地層Q值估計[J]. 石油地球物理勘探, 2011,46(3):423-428. WANG Xiaojie, YIN Xingyao, WU Guochen. Estimation of stratigraphic quality facrors on pre-stack seismic data[J]. Oil Geophysical Prospecting, 2011,46(3):423-428.

[6] 郭振波,田坤,李振春,等.利用非線性最優(yōu)化方法構建近似常Q黏彈模型[J]. 中國石油大學學報(自然科學版), 2014,38(2):52-58. GUO Zhenbo, TIAN Kun, LI Zhenchun, et al. Constructing nearly constant-Qviscoelastic model using the nonlinear optimization method[J]. Journal of China University of Petroleum(Edition of Natural Science), 2014,38(2):52-58.

[7] 王守東.基于反演的衰減補償方法[J]. 應用地球物理, 2011,8(2):150-157. WANG Shoudong. Attenuation compensation method based on inversion[J]. Applied Geophysics, 2011,8(2):150-157.

[8] 尹洪東,楊懷章,薛亞茹,等.基于混合正則化的最小二乘三維電阻率反演成像[J]. 中國石油大學學報(自然科學版), 2015,39(5):72-81. YIN Hongdong, YANG Huaizhang, XUE Yaru, et al. Least square method based on hybrid regularization for 3D resistivity inversion imaging[J]. Journal of China University of Petroleum(Edition of Natural Science), 2015,39(5):72-81.

[9] 高陽,肖立志,謝慶明. 基于BG理論確定Phillips方法的正則化參數(shù)并評價反演結果[J]. 中國石油大學學報(自然科學版), 2014,38(3):50-56. GAO Yang, XIAO Lizhi, XIE Qingming. Evaluating regular parameter of Phillips method and assessing inversions based on BG theory[J]. Journal of China University of Petroleum(Edition of Natural Science), 2014,38(3):50-56.

[10] WANG Y. A stable and efficient approach to inverseQfiltering[J]. Geophysics, 2002,67:657-663.

[11] WANG Y. InverseQfilter for seismic resolution enhancement[J]. Geophysics, 2006,71:V51-V60.

[12] 李國發(fā),張小明,彭更新,等.與炮檢距有關的地層吸收對AVO分析的影響及其補償方法[J]. 石油地球物理勘探, 2014,49(1):89-94. LI Guofa, ZHANG Xiaoming, PENG Gengxin, et al. Influence of offset-related absorption on AVO analysis and its compensation[J]. Oil Geophysical Prospecting, 2014,49(1):89-94.

[13] LI G, LIU Y, ZHENG H, et al. Absorption decomposition and compensation via a two-step scheme[J]. Geophysics, 2015,80(6):V145-V155.

[14] MARGRAVE G F. Gabor deconvolution: estimating reflectivity by nonstationary deconvolution of seismic data[J]. Geophysics, 2011,76:W15-W30.

[15] LAZARATOS S, FINN C. Deterministic spectral balancing for high fidelity AVO[C]// Society of Exploration Geophysicists. Expanded Abstracts of 74th SEG Exposition. Tulsa: Society of Exploration Geophysicists, 2004:219-223.

[16] CARLOS A M, MARGRAVE G F. Phase correction in Gabor deconvolution[C]// Society of Exploration Geophysicists. Expanded Abstracts of 75th SEG Exposition, Tulsa: Society of Exploration Geophysicists, 2005:2173-2177.

[17] LIU Y, WEI X. IntervalQinversion from CMP gathers: part 1 absorption equation[C]//Society of Exploration Geophysicists. Expanded Abstracts of 75th SEG Exposition, Tulsa: Society of Exploration Geophysicists, 2005:1717-1720.

[18] LIU Y, WEI X. IntervalQinversion from CMP gathers: part 2 inversion method[C]// Society of Exploration Geophysicists. Expanded Abstracts of 75th SEG Exposition. Tulsa: Society of Exploration Geophysicists, 2005: 1721-1724.

[19] LI G, CAO M, ZHOU H. Effects of near-surface absorption on the reflection characteristics of continental interbedded strata: the Dagang oilfield as an example[J]. Acta Geologica Sinica, 2010,84:1306-1314.

[20] DUTTA G, SCHUSTER G T. Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation[J]. Geophysics, 2014,79(6):S251-S262.

[21] MALLAT S,ZHANG Z. Matching pursuit with time-frequency dictionaries[J]. IEEE Transactions on Signal Processing, 1993,41(12):3397-3415.

[22] LIU J, WU Y, HAN D, et al. Time-frequency decomposition based on Ricker wavelet[C]// Society of Exploration Geophysicists. Expanded Abstracts of 74th SEG Exposition, Tulsa: Society of Exploration Geophysicists, 2004:1937-1940.

[23] 張繁昌,李傳輝,印興耀.基于動態(tài)匹配子波庫的地震數(shù)據(jù)快速匹配追蹤[J]. 石油地球物理勘探, 2010,45(5):667-673. ZHANG Fanchang, LI Chuanhui, YIN Xingyao. Seismic data fast matching pursuit based on dynamic matching wavelet library[J]. Oil Geophysical Prospecting, 2010,45(5):667-673.

[24] 張繁昌,李傳輝.地震信號復數(shù)域高效匹配追蹤分解[J].石油地球物理勘探,2013,48(2):171-175. ZHANG Fanchang, LI Chuanhui. Complex domain efficient matching pursuit decomposition of seismic signals[J]. Oil Geophysical Prospecting, 2013,48(2):171-175.

[25] WANG Y. Inverse-Qfiltered migration[J]. Geophysics, 2008,73:S1-S6.

[26] KOLSKY H. The propagation of stress pulses in viscoelastic solids[J]. Philosophical Magazine, 1956,1(8):693-710.

[27] BICKEL S H, NATARAJAN R R. Plane-waveQdeconvolution[J]. Geophysics, 1985,77(3):1426-1439.

[28] 張繁昌,張汛汛,張立強,等.基于自適應子波分解的品質因子Q提取方法[J].石油物探,2016,55(1):41-48. ZHANG Fanchang, ZHANG Xunxun, ZHANG Liqiang, et al. Extraction method for quality factorQbased on adaptive wavelet decomposition[J]. Geophysical Prospecting for Petroleum, 2016,55(1):41-48.

(編輯 修榮榮)

Pre-stack seismic gather absorption compensation and dispersion correction based on matching pursuit

ZHANG Fanchang1, ZHANG Xunxun1, WENG Bin2, XIE Xiuchuan3, ZHENG Yang4

(1.SchoolofGeosciencesinChinaUniversityofPetroleum,Qingdao266580,China;2.CNOOCResearchCenter,Beijing100028,China;3.HuabeiOilField,Renqiu062552,China;4.TianjinBranchofCNOOC(China)Limited,Tianjin300452,China)

Pre-stack seismic gathers are firstly decomposed into a variety of matching wavelets by utilizing the matching pursuit algorithm. The amplitude attenuation curve, which changes with the products of wavelets central times and central frequencies, can be established by using zero-offset matching wavelets. Based on this, the gain coefficients that belong to different matching wavelets attenuation amplitude can be obtained. Then, according to the absorption differences between zero-offset event and non-zero offset event, the zero-offset amplitude gain coefficients and the gather angels are combined to calculate non-zero offset events amplitude gain coefficients. Lastly, different matching wavelets attenuation amplitudes can be compensated by utilizing amplitude gain coefficients accordingly to reconstruct the seismic signal. To correct the influence of dispersion,Qvalues are estimated by matching pursuit. The comparison between the synthetic seismic signal and the real data shows that this new algorithm is efficient in both compensating the deep part seismic attention energy, and eliminating the influence of absorption in offset direction.

seismic prospecting;pre-stack gather; absorption attenuation; matching pursuit; amplitude gain coefficient; dispersion correction

2015-05-24

國家科技重大專項(2016ZX05006004)；中海油重點科研項目(YXKY-2016-ZY-09)；中國石油大學(華東)研究生創(chuàng)新工程項目(YCX2015006)

張繁昌(1972-)，男，教授，博士，研究方向為地球物理反演與儲層預測方法。E-mail：zhangfch@upc.edu.cn。

張汛汛(1991-)，男，碩士研究生，研究方向為地震波衰減補償方法。E-mail：zhang1417xun@163.com。

1673-5005(2017)04-0071-07

10.3969/j.issn.1673-5005.2017.04.009

P 631

A

張繁昌,張汛汛,翁斌,等.基于匹配追蹤的疊前道集吸收補償和頻散校正[J].中國石油大學學報(自然科學版),2017,41(4):71-77.

ZHANG Fanchang, ZHANG Xunxun, WENG Bin, et al. Pre-stack seismic gather absorption compensation and dispersion correction based on matching pursuit[J]. Journal of China University of Petroleum (Edition of Natural Science), 2017,41(4):71-77.