吳丹， 陳文超,3,4， 王曉凱,2,3,4*， 高磊,3， 高靜懷 ,2,3,4

1 西安交通大學(xué)波動(dòng)與信息研究所， 西安　710049　2 西安交通大學(xué)計(jì)算地球物理科學(xué)系， 西安　710049　3 海洋石油勘探國(guó)家工程實(shí)驗(yàn)室， 西安　710049　4 北京數(shù)學(xué)與信息交叉科學(xué)協(xié)同創(chuàng)新中心油氣勘探與開發(fā)平臺(tái)， 西安　710049

?

小波域疊前地震資料衰減參數(shù)定性估計(jì)方法

吳丹1， 陳文超1,3,4， 王曉凱1,2,3,4*， 高磊1,3， 高靜懷1,2,3,4

1 西安交通大學(xué)波動(dòng)與信息研究所， 西安7100492 西安交通大學(xué)計(jì)算地球物理科學(xué)系， 西安7100493 海洋石油勘探國(guó)家工程實(shí)驗(yàn)室， 西安7100494 北京數(shù)學(xué)與信息交叉科學(xué)協(xié)同創(chuàng)新中心油氣勘探與開發(fā)平臺(tái)， 西安710049

摘要疊前地震資料中，高頻分量和低頻分量隨傳播距離的衰減特性不同.本文給出了一種在小波域定性估計(jì)疊前地震資料衰減參數(shù)的方法.該方法利用連續(xù)小波變換提取共反射點(diǎn)道集的高、低頻分量，以低頻分量和高頻分量之差定性反映地震波的衰減.通過累加不同偏移距的衰減，提高了估計(jì)的穩(wěn)定性；采用幅度歸一化方法，降低了信號(hào)幅值對(duì)衰減參數(shù)估計(jì)的影響.將本文提出方法與常用的基于疊后地震資料衰減估計(jì)方法用于某油田的地震資料處理，結(jié)果表明，本文方法得到的衰減估計(jì)結(jié)果能夠更好地反映油氣的空間展布.

關(guān)鍵詞地震衰減; 小波變換; 共反射點(diǎn)道集

1引言

地震波在傳播過程中如果穿過含流體地層，會(huì)引起地震波的衰減.因此提取地震資料的衰減屬性可以指示地層的含流體特性，這對(duì)于直接尋找油氣藏有著重要的意義.

在20世紀(jì)70—90年代，針對(duì)衰減屬性的研究工作主要集中于對(duì)品質(zhì)因子Q的定量計(jì)算.Q值除了能夠指示地層含流體特性外，也可以對(duì)地震記錄進(jìn)行有效的反Q濾波補(bǔ)償(Wang, 2008)，進(jìn)而改善地震剖面分辨率.譜比法是常用的Q值定量估計(jì)方法之一(Bath, 1974)，但該方法估計(jì)的效果受窗長(zhǎng)度的影響很大.此外，受噪聲、波干涉等其它因素的影響，準(zhǔn)確估算Q值的難度很大.地震波在地層傳播時(shí)高、低頻分量存在衰減差異，在地層含流體時(shí)該現(xiàn)象更為明顯.因此，Castagna等人提出了在時(shí)頻域(Mallat and Zhang, 1993)定性估算地震波的衰減屬性.實(shí)際應(yīng)用效果表明，該方法對(duì)地層的含氣性有良好的指示作用(Sun et al., 2002; Castagna et al., 2003).在利用時(shí)頻分析對(duì)薄儲(chǔ)層含氣性進(jìn)行預(yù)測(cè)時(shí)，時(shí)間-頻率分辨率至關(guān)重要.為此，陳文超等人在小波域利用改進(jìn)的最佳匹配子波(MBMSW)分析地震波高、低頻之間衰減差異(高靜懷等，2006；陳文超和高靜懷，2007)，所得到的衰減特性參數(shù)具有較高的時(shí)間分辨率及良好的抗噪性能；陳學(xué)華等利用S變換(Stockwell et al.，1996；高靜懷等，2003)探測(cè)與油氣有關(guān)的低頻陰影(陳學(xué)華等，2009)，有效地刻畫了油氣儲(chǔ)層的巖性邊界和空間展布，減小了油氣儲(chǔ)層檢測(cè)的多解性；Wang Xiaokai等利用自適應(yīng)最優(yōu)核時(shí)頻表示(AOKTFR)定性描述地震衰減特性(Jones and Baraniuk，1995;Wang et al.，2013)，同時(shí)提高了時(shí)間及頻率分辨率.以上研究均是基于疊后地震資料，在實(shí)際地震處理中取得明顯的應(yīng)用效果.

然而，在地震記錄的疊加過程中，由地層吸收引起的疊前記錄波形變化在受到疊加平均效應(yīng)后會(huì)大幅減弱，進(jìn)而影響儲(chǔ)層衰減參數(shù)的預(yù)測(cè)精度.疊前地震記錄是地震信號(hào)最原始的記錄，包含偏移距、旅行時(shí)等豐富信息.為了提高儲(chǔ)層預(yù)測(cè)的精度，一些研究者們選擇疊前記錄進(jìn)行分析處理.Zhang Changjun和Ulrych以疊前共中心點(diǎn)道集(CMP)為數(shù)據(jù)基礎(chǔ)，提出了利用峰值頻率定量估計(jì)Q值的方法(Zhang and Ulrych，2002)，該方法的有效性在模型及實(shí)際地震資料分析中均得到了驗(yàn)證.由于該研究中假設(shè)地震源信號(hào)為Ricker子波且利用信號(hào)的瞬時(shí)頻率估算Q值，其應(yīng)用范圍及方法的穩(wěn)定性均受到了限制.趙靜和高靜懷等利用小波域包絡(luò)峰值處瞬時(shí)頻率法結(jié)合層位信息估計(jì)Q值(趙靜等，2013)，研究表明該方法計(jì)算精度較高，但方法的準(zhǔn)確性依賴于理論子波參數(shù)的匹配.李君君和王志章等通過擬合質(zhì)心頻率-傳播時(shí)間的斜率估算Q值(李君君等，2015)，該方法在近地表Q值估算中取得較好的應(yīng)用效果，但還未能應(yīng)用于更為復(fù)雜的地面地震資料.

在實(shí)際疊前資料處理中，受資料低信噪比的影響，定量估計(jì)Q值的精度與穩(wěn)定性受到了很大的限制.因此本文聚焦于在小波域?qū)ΟB前地震數(shù)據(jù)的衰減屬性進(jìn)行定性描述，利用地震數(shù)據(jù)的高、低頻分量之間的幅度衰減差異來估計(jì)地震衰減：在每個(gè)共反射點(diǎn)(CRP)道集中，利用連續(xù)小波變換提取高、低頻分量的幅度曲線，并利用高、低頻幅度曲線所圍成的面積刻畫地震衰減(即累加道集中每一道信號(hào)的衰減)，以此反映地層的吸收特性.

2小波域疊前地震資料衰減參數(shù)定性估計(jì)

2.1連續(xù)小波變換概述

一個(gè)地震信號(hào)f(t)的連續(xù)小波變換可表示為：

(1)

(2)

連續(xù)小波變換中僅要求母小波滿足容許條件，因此實(shí)際應(yīng)用中可以根據(jù)具體要求選擇小波函數(shù).改進(jìn)的最佳匹配子波(MBMSW)能夠很好地分析振幅快變信號(hào)，具有較高的時(shí)間分辨率，對(duì)于薄儲(chǔ)層是一種較好的時(shí)頻分析工具(高靜懷等，2006；陳文超和高靜懷，2007).因此，本文選擇MBMSW小波作為母小波函數(shù).圖1a與圖1b給出的分別是以Morlet小波及MBMSW小波為母小波對(duì)主頻為45 Hz的Ricker子波作連續(xù)小波變換后的時(shí)-頻分布圖(已經(jīng)將尺度轉(zhuǎn)換為頻率).可以看出采用MBMSW作為母小波所得到的時(shí)-頻圖的時(shí)間分辨率明顯優(yōu)于采用Morlet小波作為母小波所得到的時(shí)-頻分布圖.

圖1　利用連續(xù)小波變換得到的Ricker子波(主頻45 Hz)的時(shí)頻分布圖 (a) 采用Morlet小波作為母小波； (b) 采用MBMSW小波作為母小波.Fig.1　Using continuous wavelet transform to obtain time-frequency distribution of Ricker wavelet (dominant frequency 45 Hz) (a) Adopting Morlet as mother wavelet; (b) Adopting MBMSW as mother wavelet.

2.2衰減參數(shù)定性估計(jì)原理

由于地震波在地層中傳播時(shí)存在衰減與頻散效

應(yīng)，因此地震波的幅度與相位都會(huì)發(fā)生相應(yīng)的變化.傳播相同距離后，地震波的高頻分量比低頻分量衰減的更快；另外地震波的高、低頻分量的衰減差異會(huì)隨著儲(chǔ)層品質(zhì)因子Q的減小而增大(Aki and Richards，1980).因而，可以利用地震波的高、低頻分量的衰減差異估計(jì)地層吸收特性.

圖2　合成CRP道集的地層模型(五角星對(duì)應(yīng)炮點(diǎn)位置， 倒三角對(duì)應(yīng)接收點(diǎn)的位置)Fig.2　The model for generating CRP gather (star represents the position of shot, and the inverse triangular represents the corresponding receiver).

(3)

(4)

式中：Q表示品質(zhì)因子，v是地層速度，d(n)是第n個(gè)激發(fā)點(diǎn)到第n個(gè)檢波器的傳播距離.如果用l表示檢波器間距，h表示層厚，則d(n)可表示為

(5)

觀察圖3中的兩個(gè)CRP道集，可以看出信號(hào)的能量衰減隨著偏移距的增加而增加.另外，由于波在傳播過程中高頻分量衰減比低頻分量衰減的快，隨著偏移距的增大，導(dǎo)致信號(hào)主頻向低頻方向移動(dòng)，在時(shí)間軸上表現(xiàn)為子波展寬.比較圖3a與圖3b可以看出，信號(hào)衰減以及子波展寬現(xiàn)象在低Q值的CRP道集(圖3a)中尤為明顯.這也表明，在低Q值CRP道集中高低頻之間的衰減差異大于高Q值CRP道集中的高低頻分量之間的衰減差異.因此通過估算高低頻分量之間的衰減差異可以定性地指示出地層Q值變化.為了便于分析CRP道集中來自同一反射層的反射波隨偏移距的變化，后續(xù)分析均采用時(shí)差校正處理后CRP道集(如圖4a和圖4b所示).

圖3　不同Q值下合成的CRP道集 (a) Q=30; (b) Q=200.Fig.3　Synthetic CRP gathers with different Q factors

以上研究表明，對(duì)地層吸收特性的定性估計(jì)依賴于對(duì)CRP道集中信號(hào)高、低頻分量衰減差異的估算.在選擇高頻fH和低頻fL時(shí)，使得兩者滿足條件

(6)

圖4　不同Q值下的合成CRP道集經(jīng)過時(shí)差校正后的結(jié)果 (a) Q=30; (b) Q=200.Fig.4　CRP gathers with different Q factors (after time correction)

(7)

(8)

此時(shí)，所求AP(n)除地層傳播速度外，只與Q值及傳播距離有關(guān).

由于CRP道集中有N道信號(hào)，可以估計(jì)出N個(gè)衰減參數(shù).將道集中N個(gè)衰減參數(shù)進(jìn)行累加作為這個(gè)CRP點(diǎn)的衰減參數(shù)，記為AP，表達(dá)式為

(9)

衰減參數(shù)AP的估算，考慮了地震信號(hào)的衰減隨傳播距離影響，同時(shí)，對(duì)于實(shí)際信號(hào)，降低了異常道及隨機(jī)噪聲的影響，增強(qiáng)了方法的穩(wěn)定性.由于采用了多道平均的方法來提高抗噪性能，上式中的AP與Q并沒有簡(jiǎn)單明確的直接對(duì)應(yīng)關(guān)系，但可用來定性刻畫地層吸收作用.

在CRP道集中，遠(yuǎn)偏移距地震道在高吸收區(qū)域的傳播路程大于近偏移距在高吸收區(qū)域的傳播路程，因此遠(yuǎn)偏移距地震道上的高低頻差異與近偏移距相比較大.然而，與吸收較弱的區(qū)域相比，吸收較強(qiáng)的地層區(qū)域無論在近偏移距還是遠(yuǎn)偏移距地震道上高低頻差異均比較明顯.在實(shí)際處理中假設(shè)地層的吸收效應(yīng)在橫向變化緩慢，另外本文的主要目標(biāo)并不是嚴(yán)格估計(jì)地層Q值，而是定量估計(jì)地層的吸收.在如上基礎(chǔ)上，對(duì)一定偏移距范圍內(nèi)的高低頻進(jìn)行平均，高低頻之間的差值能夠定性反映地層吸收的相對(duì)強(qiáng)弱.

圖5　Q值為30和200的CRP道集的高頻幅度曲線與低頻幅度曲線.紅色實(shí)線與紅色虛線表示Q值為200的CRP道集的低頻幅度曲線與高頻幅度曲線，藍(lán)色實(shí)線與藍(lán)色虛線表示Q值為30的CRP道集 的低頻幅度曲線與高頻幅度曲線Fig.5　High-frequency and low-frequency amplitude curves of CRP gathers with Q=30 and Q=200. Red solid line and red dotted line indicate low-frequency amplitude curve and high-frequency amplitude curve in CRP gather with Q=200, and blue solid line and blue dotted line indicate low-frequency amplitude curve and high-frequency amplitude curve in CRP gather with Q=30.

單層模型中其它參數(shù)不變，取多組Q值對(duì)該疊前道集進(jìn)行實(shí)驗(yàn)，按照公式(9)求得衰減參數(shù)AP，如圖6藍(lán)色星號(hào)“*”所示.從圖中可看出衰減參數(shù)AP值隨著Q增大而減小.所以本文定義的衰減參數(shù)AP可以指示Q值的變化.圖6中紅色星號(hào)“*”為利用疊后地震記錄估計(jì)衰減參數(shù)估計(jì)結(jié)果(疊后數(shù)據(jù)高低頻之差)，通過比較可以發(fā)現(xiàn)，疊前道集所獲得的衰減參數(shù)AP隨著Q值變化范圍更大，這表明在實(shí)際地震資料處理中，利用疊前地震數(shù)據(jù)進(jìn)行衰減估計(jì)可以獲得更高的精度.

圖6　疊前道集與疊后數(shù)據(jù)對(duì)應(yīng)不同Q值的 衰減參數(shù)曲線Fig.6　Attenuation parameter curves of synthetic CRPs and stacked data with different Q-factors

2.3小波域衰減參數(shù)的定性估計(jì)

上節(jié)利用單層模型合成數(shù)據(jù)檢驗(yàn)了在CRP道集中進(jìn)行地震衰減定性估計(jì)的有效性.考慮到實(shí)際疊前地震資料的復(fù)雜性(存在有多個(gè)反射子波)及時(shí)變特性，利用傅里葉變換無法獲得信號(hào)的時(shí)間局部化特性，因此必須選擇合適的時(shí)頻分析工具來獲得時(shí)頻分布以提取高、低頻分量.本文利用連續(xù)小波變換為時(shí)頻分析工具并以前文提到MBMSW作為母小波，對(duì)地震信號(hào)進(jìn)行分析以獲得時(shí)頻分布.

(10)式中：n表示CRP道集中的道號(hào).由于b為時(shí)移因子，與時(shí)間具有相同的單位和相同的含義.因此為了描述簡(jiǎn)單，用變量t表示b，則s(t,n)的時(shí)間尺度分布為SW(a,t,n).利用尺度與頻率的對(duì)應(yīng)關(guān)系(高靜懷等，2006)，將尺度a變換為頻率坐標(biāo)，即得到時(shí)間頻率分布為

(11)

σ0為MBMSW的調(diào)制頻率(單位：弧度/秒).從SW(f,t,n)中提取高頻分量SW(fH,t,n)與低頻分量SW(fL,t,n)，并將小波變換引入到衰減參數(shù)定性估計(jì)公式(9)中，得

(12)

式中：SW(fH,t)為小波域震源信號(hào)的高頻分量幅度.在實(shí)際地震資料處理中用CRP道集中N道信號(hào)高頻分量幅度的平均代替源信號(hào)的高頻分量幅度以提高估計(jì)的穩(wěn)定性，則可得到

(13)

利用公式(13)，逐CRP道集進(jìn)行計(jì)算，可以得到對(duì)應(yīng)數(shù)據(jù)的衰減參數(shù)AP的剖面.該剖面橫軸為CRP號(hào)，縱軸為時(shí)間，剖面中的亮點(diǎn)區(qū)域(大AP值)反映地震波的強(qiáng)衰減區(qū)域，可以用來指示地層的吸收異常區(qū)(低Q值)，進(jìn)一步可用來指示儲(chǔ)層的含油氣性.另外，我們對(duì)圖2模型生成的合成數(shù)據(jù)(Q=30)利用不同的小波來進(jìn)行處理(Morlet小波與MBMSW小波)以分析不同小波的時(shí)間分辨率.利用Morlet小波和MBMSW小波對(duì)合成數(shù)據(jù)做衰減參數(shù)估計(jì)，獲得實(shí)驗(yàn)結(jié)果如圖7所示.從圖中可以看出MBMSW小波處理結(jié)果的時(shí)間分辨率明顯高于Morlet小波處理結(jié)果,且具有較高的AP值.

需要注意的是，在實(shí)際資料處理中，地震道集經(jīng)過動(dòng)校正之后在遠(yuǎn)偏移距會(huì)產(chǎn)生明顯的拉伸現(xiàn)象，影響本文方法對(duì)衰減參數(shù)的估計(jì).因此，在對(duì)實(shí)際疊前地震資料處理時(shí)，應(yīng)切除道集中拉伸畸變較為嚴(yán)重的部分，這樣可以有效減少動(dòng)校正的影響且改善衰減參數(shù)估計(jì)結(jié)果.另外，實(shí)際疊前數(shù)據(jù)受噪聲等因素影響，信號(hào)的幅度變化劇烈，不易直接估算幅度的高、低頻差異.所以，在實(shí)際數(shù)據(jù)分析中，可對(duì)高、低頻幅度曲線進(jìn)行平滑處理.

3實(shí)際資料算例

圖8a為某油田一條疊后測(cè)線，該測(cè)線共901個(gè)CRP道集，整個(gè)疊前數(shù)據(jù)共87911道信號(hào).有兩口探井經(jīng)過(Well 1與Well 2)該測(cè)線，分別位于CRP700與CRP800附近，其中Well 1在儲(chǔ)層附近獲得工業(yè)氣流，而Well2在儲(chǔ)層附近未見工業(yè)氣流.圖8b為圖8a中第700、702及704個(gè)CRP道集(已經(jīng)過時(shí)差校正)，可以看出相較于疊后剖面疊前記錄顯示出更豐富的信息.

圖8　實(shí)際地震記錄 (a) 疊后地震剖面； (b) 疊前地震剖面(CRP道集).Fig.8　Real seismic records (a) Stacked seismic section; (b) Pre-stack seismic section (CRP gathers).

圖7　不同小波估計(jì)的衰減參數(shù)分辨率比較 (Morlet小波和MBWSW小波)Fig.7　Comparison of time-resolution of different wavelets

將本文中衰減參數(shù)估計(jì)方法應(yīng)用于該疊前數(shù)據(jù)，得到衰減參數(shù)剖面，如圖9所示.圖中亮點(diǎn)區(qū)域(紅色點(diǎn))為強(qiáng)衰減區(qū)域，說明該處吸收異常，表明該段儲(chǔ)層含氣.圖中紅色框中指示了儲(chǔ)層位置，兩條紅色曲線分別為Well 1與Well 2的電阻率曲線.可以看出Well1的電阻率曲線高值區(qū)域與衰減參數(shù)剖面亮點(diǎn)區(qū)域位置吻合，而Well2的電阻率曲線沒有明顯高值，這也同樣與衰減參數(shù)剖面吻合.

圖9　衰減參數(shù)(AP)剖面以及Well 1與Well 2的電阻率曲線Fig.9　Estimated attenuation parameter (AP) section and resistivity curves of Well 1 and Well 2

分別分析Well 1與Well 2對(duì)應(yīng)的CRP700與CRP800疊前道集，分別如圖10a與10b所示.圖中紅色框?yàn)閮?chǔ)層對(duì)應(yīng)的時(shí)間，可以看出，隨偏移距的增加，CRP700中信號(hào)比CRP800中衰減更明顯.對(duì)兩個(gè)CRP道集進(jìn)行小波分析，提取出相應(yīng)的高頻幅度曲線與低頻幅度曲線，如圖11a與11b所示.CRP700中高低頻幅度曲線之間的差異更大，這與該位置獲得工業(yè)氣流非常一致.

圖10　疊前地震道集 (a) CRP700; (b) CRP800.Fig.10　Pre-stacked gathers

圖11　CRP道集高頻幅度曲線與低頻幅度曲線 (a) CRP 700; (b) CRP 800.Fig.11　Amplitude curves of high-frequency and low frequency

我們將本文方法與常用的基于疊后資料衰減估計(jì)方法進(jìn)行了對(duì)比.圖12給出了基于AOKTFR的疊后衰減特性剖面(Wang et al.，2013).此方法在儲(chǔ)層高含氣區(qū)域也顯示相對(duì)高值，但與本文方法結(jié)果(圖9)相比，高含氣區(qū)域異常變化不明顯.另外，圖12中箭頭所示位置的明顯異常區(qū)域?qū)?yīng)于疊后數(shù)據(jù)的高振幅區(qū)域，與地層衰減無關(guān).本文方法包括了高頻分量幅度歸一化，因此能夠有效地減弱了信號(hào)幅度對(duì)衰減參數(shù)的影響.

圖12　基于疊后資料衰減估計(jì)方法得到的衰減參數(shù)剖面Fig.12　Attenuation section obtained by attenuation estimating method for stacked data

4結(jié)論

本文提出的基于疊前地震數(shù)據(jù)的衰減屬性定性估計(jì)方法是有效的.通過典型的合成數(shù)據(jù)及實(shí)際地震資料分析可以看出，本文方法定義的衰減參數(shù)AP用于定性分析地層吸收特性效果良好.相比于以往的基于疊后地震數(shù)據(jù)衰減參數(shù)估計(jì)方法，本文方法通過累加CRP道集中每道信號(hào)的衰減作為一個(gè)CRP點(diǎn)的衰減.由于本文方法考慮了地震信號(hào)衰減隨波傳播距離的變化，因此提高了預(yù)測(cè)精度和信噪比.另外，通過采用平均高頻分量幅度歸一化的方法，有效地降低了信號(hào)幅值對(duì)吸收參數(shù)估計(jì)精度的影響.

致謝感謝審稿專家和編輯部的大力支持.

References

Aki K, Richards P G. 1980. Quantitative Seismology. San Francisco: W. H. Freeman and Company. Bath M. 1974. Spectral Analysis in Geophysics: Developments in Solid Earth Geophysics. Amsterdam: Elsevier.

Castagna J P, Sun S J, Siegfried R W. 2003. Instantaneous spectral analysis: Detection of low-frequency shadows associated with hydrocarbons. The Leading Edge, 22(2): 120-127.

Chen, W C, Gao J H. 2007. Characteristic analysis of seismic attenuation using MBMSW wavelets. Chinese J. Geophys. (in Chinese), 50(3): 837-843.

Chen X H, He Z H, Huang D J, et al. 2009. Low frequency shadow detection of gas reservoirs in time-frequency domain. Chinese J. Geophys. (in Chinese), 52(1): 215-221.

Gao J H, Chen W C, Li Y M, et al. 2003. Generalized S transform and seismic response analysis of thin interbeds. Chinese J. Geophys. (in Chinese), 46(4): 526-532.

Gao J H, Wan T, Chen W C, et al. 2006. Three parameter wavelet and its applications to seismic data processing. Chinese J. Geophys. (in Chinese), 49(6): 1802-1812.

Jones D L, Baraniuk R G. 1995. An adaptive optimal-kernel time-frequency representation. IEEE Transactions on Signal Processing, 43(10): 2361-2371. Li J J, Wang Z Z, Zhang Z H, et al. 2015. Quality factor Qestimation based on linear fitting of centroid frequency curve. Oil Geophysical Prospecting (in Chinese), 50(2): 254-259.

Mallat S G, Zhang Z F. 1993. Matching pursuits with time-frequency dictionaries. IEEE Transactions on Signal Processing, 41(12): 3397-3415. Stockwell R G, Mansinha L, Lowe R P. 1996. Localization of the complex spectrum: the S transform. IEEE Transactions on Signal Processing, 44(4): 998-1001.

Sun S J, Castagna J P, Seigfried R W. 2002. Examples of wavelet transform time-frequency analysis in direct hydrocarbon detection. ∥ 2002 SEG Annual Meeting. Salt Lake City, Utah: SEG, 457-460.Wang X K, Gao J H, Chen W C, et al. 2013. Seismic attenuation qualitative characterizing method based on adaptive optimal-kernel time-frequency representation. Journal of Applied Geophysics, 89: 125-133.

Wang Y H. 2008. Inverse-Q filtered migration. Geophysics, 73(1): S1-S6.

Zhang C J, Ulrych T J. 2002. Estimation of quality factors from CMP records. Geophysics, 67(5): 1542-1547.

Zhao J, Gao J H, Wang D X, et al. 2013. Estimation of quality factor Q from pre-stack CMP records. Chinese J. Geophys. (in Chinese), 56(7): 2413-2428, doi: 10.6038/cjg20130727.

附中文參考文獻(xiàn)

陳文超, 高靜懷. 2007. 基于改進(jìn)的最佳匹配地震子波的地震資料衰減特性分析. 地球物理學(xué)報(bào), 50(3): 837-843.

陳學(xué)華, 賀振華, 黃德濟(jì)等. 2009. 時(shí)頻域油氣儲(chǔ)層低頻陰影檢測(cè). 地球物理學(xué)報(bào), 51(1): 215-221.

高靜懷, 陳文超, 李幼銘等. 2003. 廣義S變換與薄互層地震響應(yīng)分析. 地球物理學(xué)報(bào), 46(4): 526-532.

高靜懷, 萬濤, 陳文超等. 2006. 三參數(shù)小波及其在地震資料分析中的應(yīng)用. 地球物理學(xué)報(bào), 49(6): 1802-1812.

李君君, 王志章, 張枝煥等. 2015. 質(zhì)心頻率線性擬合法估算品質(zhì)因子Q. 石油地球物理勘探, 50(2): 254-259.

趙靜, 高靜懷, 王大興等. 2013. 利用疊前CMP資料估計(jì)介質(zhì)品質(zhì)因子. 地球物理學(xué)報(bào), 56(7): 2413-2428, doi: 10.6038/cjg20130727.

(本文編輯劉少華)

基金項(xiàng)目國(guó)家自然科學(xué)基金 (41390454、41274125、41504092)資助.

作者簡(jiǎn)介吳丹，女，漢族，碩士，主要從事地震信號(hào)處理方法研究. *通訊作者王曉凱，男，漢族，博士，主要從事地震信號(hào)處理和地震儲(chǔ)層解釋方法研究. E-mail: xkwang@mail.xjtu.edu.cn

doi:10.6038/cjg20160730 中圖分類號(hào)P631

收稿日期2015-08-18，2016-03-03收修定稿

Qualitative estimation of seismic attenuation parameters based on pre-stack data in the wavelet domain

WU Dan1, CHEN Wen-Chao1,3,4, WANG Xiao-Kai1,2,4*, GAO Lei1,3, GAO Jing-Huai1,2,3,4

1Instituteofwaveandinformation,Xi′anJiaotongUniversity,Xi′an710049,China2Departmentofcomputationalgeophysicalscience,Xi′anJiaotongUniversity,Xi′an710049,China3Nationalengineeringlaboratoryforoffshoreoilexploration,Xi′an710049,China4BeijingCenterforMathematicsandInformationInterdisciplinaryScience,Xi′an710049,China

AbstractIn seismic wave propagation, the high-frequency component and low-frequency component have different attenuating properties with propagating distance, especially in strong attenuating areas. Many time-frequency analyzing tools are used to extract high-frequency component and low-frequency component of stacked seismic data to characterize attenuation quantitatively or qualitatively. However, the stacking operation on pre-stacked seismic data would affect characterizing precision seriously. In order to improve precision and robustness, we propose a qualitative attenuation estimating method for pre-stack seismic data with the help of continuous wavelet transform. This method uses continuous wavelet transform to extract the high-frequency component and low-frequency component from the common reflection point gather data, and adopts the modified best matching seismic wavelet as mother wavelet. The difference between low-frequency component and high-frequency component is used to reflect seismic attenuation qualitatively. The summation of attenuation along offset is adopted to enhance the robustness of estimation, while the amplitude normalization is adopted to reduce the affection of amplitude.

The proposed method is applied to a field data which contains 901 CRP gathers. Two wells are located in this area. The results show that the proposed method obtains a seismic attenuation parameter section which accords with drilling results very well. The qualitative seismic attenuation estimating method based on stacked data is also applied to a field data. Compared with the method based on stacked data, the proposed method has better anti-noise ability and is less likely to be affected by amplitude.

KeywordsSeismic attenuation； Continuous wavelet transform; Common reflection gather

吳丹， 陳文超， 王曉凱等. 2016. 小波域疊前地震資料衰減參數(shù)定性估計(jì)方法.地球物理學(xué)報(bào)，59(7):2674-2683,doi:10.6038/cjg20160730.

Wu D, Chen W C, Wang X K, et al. 2016. Qualitative estimation of seismic attenuation parameters based on pre-stack data in the wavelet domain. Chinese J. Geophys. (in Chinese),59(7):2674-2683,doi:10.6038/cjg20160730.