基于Ricker類地震子波匹配追蹤

1、基于Ricker類地震子波的匹配追蹤摘要Ricker類地震子波是由Ricker子波演變?nèi)〉玫牧阆辔蝗齾?shù)理論子波，通過適當(dāng)?shù)暮喕?，取得二參?shù)C子波，參量獨立地操縱波形特點。以兩個波形參數(shù)（主旁瓣極值比和主旁瓣寬度比）來考察子波，發(fā)覺Ricker類地震子波比常規(guī)的Ricker子波、寬帶Ricker子波波形更為豐碩?；贑子波（Ricker類子波）的匹配追蹤算法，以該類子波構(gòu)建原子庫，成立包絡(luò)瞬時頻率與子波參量的定量關(guān)系，通過三步法實現(xiàn)地震信號的分解。模型試算結(jié)果說明了該算法的有效性，時頻譜具有較高的時頻聚焦性，C子波比Ricker子波更適用于時變記錄的匹配?；贑子波庫對實際資料進行匹配追蹤，分

2、頻結(jié)果取得了較好的說明成效。關(guān)鍵詞Ricker類子波主旁瓣極值比主旁瓣寬度比時變記錄匹配追蹤ABSTRACTRicker-likeseismicwaveletisakindofzerophasewaveletwith3parametersdevelopedfromRickerwavelet.Throughpropersimplification,obtainanewkindofwaveletnamedCwavelet,whichhas2parameterscontrollingwaveformindependently.Studyingwaveletbymagnituderatioandwid

3、thratioofmainlobeandsidelobe,wefindCwaveletswaveformisricherthanRickerwaveletandwide-bandRickerwavelet,andCwaveletismoreadapttotime-varyingsignal.ComparedwithSTFTandST,thematchingpursuitalgorithmbasedonCwavelethasabettertime-frequencyresolution.Spectraldecompositionofseismicdataachievedgoodresults.K

4、eywords:Ricker-likewavelet,mainlobeandsidelobemagnituderatio,mainlobeandsidelobewidthratio,time-varyingseismogram,matchingpursuit1引言最近幾年來，信號時頻分解技術(shù)在地震資料處置和說明中取得了普遍的應(yīng)用，如高分辨率處置1、姓類檢測2向等。常規(guī)的時頻分解（如窗口傅里葉變換和小波變換等），時頻窗受不確信準(zhǔn)那么約束，時頻分辨率難以同時達到最正確。為了克服常規(guī)方式的缺點，Mallat等4（1993）和Qian等（1994）5提出了匹配追蹤分解算法（Matchingpursui

5、t,簡稱MP。MPW法基于最正確匹配原那么，通過不斷的迭代，從原子庫當(dāng)選取與殘余信號結(jié)構(gòu)最正確匹配的原子，將原始信號分解為一系列原子。由于在迭代尋優(yōu)進程中，時窗長度均由信號局部特點決定，因此這些最優(yōu)原子代表了信號的局部特點，能較好地描述時變地震信號的時頻散布特性。經(jīng)典的MP算法是基于Gabor函數(shù)構(gòu)建超完備子波庫，為了適應(yīng)信號不同的結(jié)構(gòu)特點，LiuJ67等前后提出了基于Ricker子波（2004）和Morlet子波（2005）的匹配追蹤算法。對研究薄互層儲層特點而言，Morlet小波分析薄層的成效欠佳；宋新武等8（2020）以為Ricker子波波形簡單，延遲時刻短，收斂較快，基于Ricker子

6、波匹配追蹤算法分解技術(shù)在薄互層分析中具有較高的時頻分辨率，能較好地刻畫地質(zhì)體形態(tài)。原子庫的合理選取，關(guān)于信號的匹配追蹤分解成效具有重要作用。Ricker類地震子波是由Ricker子波演變?nèi)〉玫牧阆辔焕碚撟硬?，具有三個操縱參量，旁瓣能量的收斂速度可變，延遲時刻可控，子波庫波形豐碩，與地震子波的匹配較為靈活90筆者在研究子波波形隨三操縱參量轉(zhuǎn)變關(guān)系的基礎(chǔ)上，對該子波作適當(dāng)?shù)暮喕?，成立參量與波形特點之間的半定量關(guān)系。最后，以簡化的子波構(gòu)建原子庫，通過三步法原理10實現(xiàn)信號的匹配追蹤分解。2大體原理2.1Ricker類地震子波Ricker類地震子波頻域解析式kWfAfmexpff0（1）g1 TOC

7、o 1-5 h z _2m112其中，m、g和k為待定參量，均取正數(shù)，Ak212m1kgk為能量歸一化系數(shù)。一樣地，關(guān)于確信的k值，m值越大，高頻能量相對越強；關(guān)于確信的m值,k值越大，高頻能量相對越弱，同時頻譜的有效頻帶變窄，致使子波信號的旁瓣轉(zhuǎn)變復(fù)雜。文獻9建議k值不宜過大。事實上，參量k和m都是操縱波譜高低頻相對強弱的參數(shù)，一樣情形下，相同的增（減）量，m對頻譜形態(tài)的改變較為平緩。卜面的討論中，k取固定值2,將式（1）子波簡化并作如下改寫:2c2ff7(2)Wk2fAfmexp-g2AfmexpmAf2exp2fm其中，待定參量cm,fmgJm,Cf為二參量Ricker類地震子波,22簡

8、稱C子波。通過簡化和改寫，將子波譜表述成Ricker譜（主頻為fm）的c次方（c1對應(yīng)Ricker子波）。fm為波譜主頻，操縱波譜的等比例伸縮；c操縱波譜主頻雙側(cè)能量的衰減速度。下面考察參量c與子波波形的關(guān)系。第一，針對持續(xù)相位數(shù)小于1的對稱波形（見圖1）,概念主旁瓣極值比PR和主旁瓣寬度比WR兩個參量，用于描述子波的波形特點。主旁瓣極值比（PR）：旁瓣極值與主瓣極值之比的絕對值；主旁瓣寬度比（WR）：旁瓣寬度（極小值點間距）與主瓣寬度（零值點問距）之比。圖2為C子波波形特點值（紅色為極值比PR線，藍色為寬度比WR線）隨參量c的轉(zhuǎn)變曲線。隨著c的增加，主旁瓣極值比增加，旁瓣能量越強；主旁瓣寬度

9、比增加，即在相同主瓣寬度的情形下，旁瓣能量達到極值需要的時刻變長。因此，對C子波有如下結(jié)論：參數(shù)fm操縱子波的橫向伸縮，fm越大，子波越短;參數(shù)c操縱子波的波形形態(tài)，c越大，子波的主旁瓣極值比、寬度比越大，持續(xù)相位數(shù)增加。當(dāng)c0時，C子波趨近于脈沖信號，PR0,WR1,當(dāng)c時，C子波趨近于單頻信號，PR1,WR2。W和5電i就極值比PR上參辮電度比HE圖1對稱C子波示意圖2主旁瓣寬度比（藍）和極值比（紅）曲線考察Ricker子波和寬帶Ricker子波，那么有：Ricker子波子波的PR值恒為2e32,WR包為43;寬帶Ricker子波的PR值可變，WR恒為。3。對照三類子波發(fā)覺,C子波（Ric

10、ker類地震子波）的PR和WR值都可變，波形比Ricker子波、寬帶Ricker子波更豐碩，操縱參量意義明確，選取較為方便，更適用于實際時變地震子波的匹配?；贑子波的匹配追蹤匹配追蹤是基于原子庫掃描的信號自適應(yīng)分解。記D為進行信號分解的超完備子波庫，概念為Dw:,w為單個原子，知足歸一化條件|w|1。記待分解信號為st,經(jīng)N次迭代分解如下：N1stanwntRNs（3）n,nn0其中，w,n為第n次迭代取得的匹配子波，an為w,n對應(yīng)的振幅，RNs為迭代N次后的殘差。匹配追蹤算法中子波庫的波形與原信號結(jié)構(gòu)越相似，分解的成效越好，那個地址選擇C子波來構(gòu)建匹配子波庫。匹配子波w,n由參量ntn,

11、fm,n,Cn,n來操縱，其中tn、n為子波的中心時刻和相位，%n、g為波譜操縱參量。匹配追蹤算法是基于不斷迭代的貪婪算法，常規(guī)算法的每一次迭代都要從子波庫中尋覓最正確子波。目前，通常采納“三步法”原理，提高分解速度。第一步，確信子波參量初始值0t0,f；,n，c0，:。對輸入信號RnS進行Hilbert轉(zhuǎn)變，提取“三瞬”屬性，將最大瞬時振幅對應(yīng)的時刻作為匹配子波的初始延遲時t：，瞬時相位作為子波的初始相位0，瞬時頻率作為子波的初始瞬時頻率fcn。按文獻11，包相位子波包絡(luò)峰上取得的瞬時頻率fcn即振幅譜的質(zhì)心頻率，那么有0c,nfC f df 0C f df o(4)初始值cn取為1,f；n

12、按式（4）求?。窗j(luò)日ft時頻率fcn與譜峰值頻率f；n之間進行換算）。第二步，優(yōu)化子波參量。對匹配子波的每一個操縱參量在初始參量周圍掃描,如參量t,其掃描范圍為t0tn，t0么，取得匹配子波的最優(yōu)參量now t arg n第三步，估算最優(yōu)子波的幅值。按式tn o在掃描范圍內(nèi)，依照式(5)準(zhǔn)那maxARn(5)WnWn(6)取得最優(yōu)子波的幅值an ,Rns-wn2w每一次迭代，都將生成匹配子波anW,n和殘差信號Rn1s0分解完成后，即可用提取的匹配子波近似表示原始信號，N1stanW,nt(7)n0時頻表征常規(guī)的匹配追蹤時頻表征通過各匹配子波的Wigner-Ville散布疊加實現(xiàn)，但只能給

13、出振幅散布，無法取得相位信息。那個地址，采納式(8)所示的時頻表征方法12,對每一個匹配子波的復(fù)譜求和，獲取信號的時頻譜。N1At,fanW,nfenvw,nt(8)n0其中，W,nf是匹配子波w,nt的頻譜，envw,nt是w,nt的瞬時包絡(luò)。圖3是單個C子波(主頻幻取50,c取)的時頻譜。該時頻表征法具有較高的時頻分辨率，能同時獲取振幅和相位等信息的時頻散布，而且幸免了窗口截斷效應(yīng)。圖3 C原子時頻譜專工一寺3模型試算設(shè)計兩個合成信號進行基于C子波庫的匹配追蹤算法測試，時刻采樣率均為1ms信號1由多個C原子組成，第1、2層是兩個不同參量的C原子，第3層由兩個C原子合成，第4、5、6層為薄反

14、射層，層數(shù)慢慢增加，信號組成方式見圖5左，第一道為合成記錄，為其余四道的疊加。圖4為信號1基于C子波的匹配追蹤重構(gòu)結(jié)果，左為合成記錄，中為重構(gòu)記錄，右為重構(gòu)相對誤差百分比，維持在2%：內(nèi)，對信號的損傷較?。ㄖ貥?gòu)精度與算法終止條件有關(guān)，精度越高，運行時刻越長）。圖5右為信號1的匹配追蹤時頻譜圖，。圖6左為信號1的短時傅里葉變換（STFD譜，右為S變換譜。對照三類時頻分解譜圖發(fā)覺，在600ms處，S變換和STFT譜的時刻分辨率不夠，匹配追蹤那么較好地識別出兩個反射；在700m麗800ms處，匹配追蹤譜的時頻分辨能力明顯高于STFT和S變換譜。因此，匹配追蹤時頻譜具有較高的時頻聚焦性。2003004

15、00500600700800900匹配結(jié)果相對誤差Signal圖4信號1重構(gòu)圖（左：原始信號；中：重構(gòu)信號；右：相對誤差為a n9 SilO o O noTime-frequency Spectrum50100150Frequency/HzJI200圖5信號1的分解結(jié)果（左：合成記錄組成示意；右：基于C子波的匹配追蹤時頻譜）圖6時頻譜（左：短時Fourier變換；右：S變換）信號2由Ricker類理論子波與隨機反射系數(shù)褶積生成。為符合實際情形，由淺到深理論子波的主頻和帶寬慢慢減小，且主頻與帶寬的比值隨機改變。圖7從左到右依次編為a、b、c、d,其中，a為隨機反射系數(shù)，b為合成記錄，從上往下信號

16、的分辨率降低，c為匹配追蹤結(jié)果，d為相對誤差百分比。為作對照，對信號2同時進行基于Ricker子波的匹配追蹤（算法的終止條件不變），時頻譜圖見圖8。其中，左為信號，中為基于C子波的匹配分解，右為基于Ricker子波的匹配分解。對照發(fā)覺，兩個分解結(jié)果存在必然的不同：C子波匹配追蹤譜主頻隨時刻呈近似下降的趨勢，符合模型設(shè)計，Ricker子波匹配譜在600ms后主頻呈遞增的趨勢，與實際不符；C子波匹配譜的強能量團與信號的峰值散布符合度較高，如信號2在680ms-740ms之間有兩個強峰值，C子波匹配譜上有兩個強能量團，而Ricker子波匹配譜上那么有三個，在800ms周圍存在強峰值，C子波匹配譜上存

17、在強能量顯示，Ricker子波匹配譜上能量不明顯。d（a：反射系數(shù)；b：合成記錄；c：匹配結(jié)果；d：相對誤差）O90o o O 5 6 7基于匚子速基于Ricker子波50100160200。501001S0200Freauenev/HzFreauenev/Hz圖8匹配追蹤時頻圖（左：信號：中：基于C子波；右：基于Ricker子波）對照結(jié)果說明，匹配追蹤的分解成效依托于原子庫與信號大體結(jié)構(gòu)的匹配度。C原子庫的參量調(diào)控靈活，波形較Ricker子波庫豐碩，對實際地震信號大體結(jié)構(gòu)（時變地震子波）的適應(yīng)性更強。因此，關(guān)于子波特點值（主旁瓣極值比和寬度比）時變的地震信號，基于C子波的匹配成效優(yōu)于基于Ri

18、cker子波的匹配成效。4實際資料處置圖9為陸上某探區(qū)的過井測線剖面。該區(qū)地下構(gòu)造簡單，地層平緩，斷裂較少，淺表氣遍及全區(qū)，第四系沉積以砂泥巖互層為主，成巖性較差，特殊的地質(zhì)條件造就了該區(qū)地震資料的特殊性。經(jīng)勘探發(fā)覺該井存在工業(yè)氣流，地質(zhì)綜合研究以為該地域氣藏為鼻狀構(gòu)造背景上的巖性氣藏。依照以往體會，含油氣地域一樣會有明顯地震異樣。圖9顯示，在井周圍同相軸的持續(xù)性變差，有下拉趨勢，主頻變低，在主測線和聯(lián)絡(luò)線上地震相外形呈塔狀。這種地震異樣盡管存在，但并非專門明顯。針對該地域的氣藏特點，應(yīng)用基于C子波的匹配追蹤算法，并對資料進行分頻處置，為該地域開展巖性氣藏研究提供依據(jù)。圖9某工區(qū)過井測線剖面圖

19、10是對過井測線進行分頻處置取得的單頻調(diào)諧體剖面。圖中紅線為氣層所處深度，觀看井周圍能量，在12Hz左右達到最大，以后隨著頻率的增加而降低，這是因為氣層的調(diào)諧頻率在12Hz左右。在紅線與測井相交位置周圍，右邊的能量在12Hz以后慢慢減小，而左側(cè)在25Hz剖面上仍有較強的能量，這是由于左側(cè)含氣層的阻礙。6Hz8Hz10Hz12Hz14Hz16Hz18Hz25Hz30Hz圖10單頻調(diào)諧體剖面5小結(jié)Ricker類子波是一類新型的理論子波，波譜主頻和高低頻的衰減速度由三參量聯(lián)合操縱。依如實際需要，將該類子波簡化成二參量C子波，波譜主頻、帶寬由二參量獨立操縱。針對子波波形，提出以主旁瓣極值比PR和主旁瓣

20、寬度比WR來定量表征，并論述了C子波波形比Ricker子波、寬帶Ricker子波更豐碩，更能匹配時變地震子波。依照C子波成立原子庫，用于信號的匹配追蹤算法。模型試算說明，與短時Fourier變換和S變換相較，匹配追蹤取得的時頻譜圖具有較好的時頻聚焦性；關(guān)于時變地震記錄，C子波與信號大體結(jié)構(gòu)的匹配度較Ricker子波高，匹配追蹤成效更好。最后，將基于C子波的匹配追蹤應(yīng)用于實際資料的分頻處置，取得了較好的說明成效。參考文獻熊曉軍，賀振華等廣義S變換在地震高分辨處置中的應(yīng)用勘探地球物理進展，2006，29(6)：415418CastagnaJP，SunS，SiegfriedRWInstantaneo

21、usspectralanalysis：Detectionoflow-frequencyshadowsassociatedwithhydrocarbonsJTheLeadingEdge，2003，22(2)：120127SinhaS，RouthPS，etalSpectraldecompositionofseismicdatawithcontinuous-wavelettransformsJGeophysics，2005，70(6)：1925MallatS，ZhangZMatchingpursuitswithtime-frequencydictionariesIEEETransactionsonSignalProcessing，1993，41(12)：33973415QianS，Che