戴永壽,張彧豪,張 鵬,張紅倩,王春嫻

(中國(guó)石油大學(xué)(華東)海洋與空間信息學(xué)院,山東青島266580)

隨著油氣勘探目標(biāo)埋深和復(fù)雜程度的加大,地震流體識(shí)別難度也在不斷加大,這對(duì)地震剖面的分辨率及反演解釋的精確度提出了更高要求。由于地震介質(zhì)對(duì)地震波能量的吸收衰減作用,造成子波高頻成分缺失和相位畸變,從而使子波表現(xiàn)出時(shí)變特性。以往的地震子波提取方法往往基于子波時(shí)不變的假設(shè),因此提取出的子波無(wú)法滿(mǎn)足地震勘探所需的精度要求。尚新民等[1]利用時(shí)不變子波和時(shí)變子波分別對(duì)實(shí)際地震資料進(jìn)行處理,結(jié)果利用時(shí)變子波處理的地震數(shù)據(jù)分辨率更高、有效頻帶更寬;李振春等[2]在S域利用時(shí)變子波與地震記錄進(jìn)行反褶積,使地震剖面的同相軸變細(xì)且連續(xù)性變好,淺、中、深層的能量都得到了補(bǔ)償,且深層能量改善明顯;遲喚昭等[3]也指出,相對(duì)前人采用的非時(shí)變子波,采用時(shí)變子波作譜反演處理后的地震剖面分辨率更高。故時(shí)變地震子波的提取是目前高分辨率地震資料處理方法的研究重點(diǎn)之一。目前從非平穩(wěn)地震記錄中提取時(shí)變地震子波的方法主要有兩類(lèi):一類(lèi)是先對(duì)非平穩(wěn)地震記錄進(jìn)行衰減補(bǔ)償或分段處理,使其整體或每段近似成為平穩(wěn)地震記錄,然后采用常規(guī)方法提取地震子波;另一類(lèi)是從非平穩(wěn)地震記錄中提取子波振幅譜與相位譜,再將二者匹配融合從而提取時(shí)變地震子波。本文主要概括了近年主流的時(shí)變地震子波提取方法,對(duì)每種方法的原理、優(yōu)缺點(diǎn)和適用條件進(jìn)行了分析,并對(duì)地震子波提取方法的研究前景進(jìn)行了展望。

1 基于衰減補(bǔ)償和分段處理的子波提取方法

實(shí)際地震資料中的地震子波具有時(shí)變特征,即地震資料具有非平穩(wěn)性[4]。對(duì)非平穩(wěn)地震記錄進(jìn)行衰減補(bǔ)償或分段處理,使其整體或分段成為近似平穩(wěn)地震記錄,然后結(jié)合常規(guī)方法提取地震子波的方法成為目前時(shí)變地震子波提取的研究方向之一。

1.1 衰減補(bǔ)償方法

隨著對(duì)地下介質(zhì)研究的深入,關(guān)于大地對(duì)地震波的吸收機(jī)理研究目前存在多種理論[5],主要包括Biot理論和噴射流動(dòng)理論、散射、電化學(xué)作用和毛細(xì)管作用、巖石顆粒間的摩擦滑移等。地震波在地層中的衰減和頻散一般可以用經(jīng)驗(yàn)公式來(lái)表示,國(guó)內(nèi)外學(xué)者建立了Power Law模型、Kjartansson模型、Kolsky-Futterman模型、Müller模型等多種模型,其中Kolsky-Futterman模型能更好地描述地震波在地下介質(zhì)中的傳播規(guī)律且簡(jiǎn)單方便,在子波的衰減補(bǔ)償中被廣泛應(yīng)用。該模型中Q值基本上不隨頻率變化,可以看作恒Q模型。

目前發(fā)展比較成熟的衰減補(bǔ)償方法主要包括:譜白化方法、反Q濾波方法、時(shí)頻域衰減補(bǔ)償方法和多尺度聯(lián)合分析方法等,本節(jié)對(duì)上述主流方法的優(yōu)缺點(diǎn)進(jìn)行了概述。由于多尺度聯(lián)合分析方法依賴(lài)測(cè)井資料,無(wú)法對(duì)整個(gè)地震記錄進(jìn)行補(bǔ)償,譜白化方法基于反射系數(shù)序列是白噪的假設(shè),在實(shí)際應(yīng)用中具有較大的限制,故在此均不作過(guò)多論述。

1.1.1 反Q濾波方法

地震波經(jīng)地層傳播后能量被損耗,高頻部分衰減嚴(yán)重,此時(shí)可將大地看作一個(gè)低通濾波器,此過(guò)程被稱(chēng)為Q濾波。反Q濾波是Q濾波的反過(guò)程,是針對(duì)地層的低通濾波特性,對(duì)地震數(shù)據(jù)進(jìn)行的補(bǔ)償處理。目前反Q濾波方法已經(jīng)較為成熟,但仍存在以下兩方面問(wèn)題:①反Q濾波的準(zhǔn)確度很大程度取決于Q值的求取,而Q值往往難以精確求取;②目前計(jì)算效率高的反Q濾波方法都假設(shè)地下介質(zhì)是一個(gè)常Q或者層狀Q模型,這并不能反映地下介質(zhì)的真實(shí)情況[6]。

1.1.2 時(shí)頻域衰減補(bǔ)償方法

譜白化方法、反Q濾波方法和多尺度聯(lián)合分析方法由于各自的缺陷在非平穩(wěn)地震記錄高頻補(bǔ)償中的應(yīng)用均受到了一定限制,近年來(lái),專(zhuān)家和學(xué)者又提出了許多新的補(bǔ)償方法。其中時(shí)頻分析方法因其可以同時(shí)聯(lián)合時(shí)間域和頻率域信息對(duì)信號(hào)進(jìn)行處理,克服了傅里葉變換時(shí)域和頻域完全分離的缺陷,成為目前地震資料處理的有力工具。但此類(lèi)方法在對(duì)地震記錄做衰減補(bǔ)償時(shí)假設(shè)地震記錄的低頻成分是近似不衰減的,與實(shí)際情況不符[7]。MARGRAVE等[8]在Gabor域進(jìn)行衰減補(bǔ)償和反褶積以估計(jì)地層的反射系數(shù),補(bǔ)償了中、深層的能量衰減并因此拓寬有效頻帶和提高時(shí)間分辨率,但在地層衰減劇烈時(shí),此方法效果不理想。常用的時(shí)頻分析方法還有短時(shí)傅里葉變換、小波變換、S變換、改進(jìn)S變換和Curvelet變換等。

1.1.3 其它方法

WANG等[9]提出了基于貝葉斯反演的地層衰減補(bǔ)償方法,在反射系數(shù)序列稀疏的假設(shè)下,先構(gòu)造衰減矩陣,然后利用貝葉斯理論以實(shí)際地震資料為先驗(yàn)信息求解反演問(wèn)題,實(shí)現(xiàn)了對(duì)非平穩(wěn)地震記錄的補(bǔ)償,無(wú)須精確求取Q值,且能夠較好地壓制隨機(jī)噪聲,具有良好的穩(wěn)定性,但該方法需要從實(shí)際測(cè)井資料中獲取反射系數(shù)的分布特征。BAAN[10]將反Q濾波與時(shí)變Wiener反褶積相結(jié)合進(jìn)行衰減補(bǔ)償,首先利用穩(wěn)定的反Q濾波技術(shù)對(duì)地震記錄的衰減能量進(jìn)行補(bǔ)償,同時(shí)進(jìn)行相位校正,消除其非平穩(wěn)性,然后利用Wiener反褶積來(lái)壓縮地震子波,擴(kuò)展頻帶寬度,提高地震記錄的分辨率[11]。

衰減補(bǔ)償方法在實(shí)際應(yīng)用中有著一定的限制條件,且人為對(duì)非平穩(wěn)地震記錄進(jìn)行高頻衰減補(bǔ)償很容易出現(xiàn)欠補(bǔ)償或過(guò)補(bǔ)償?shù)那闆r,難以保證地震記錄的真實(shí)性。因此,基于衰減補(bǔ)償?shù)臅r(shí)變子波提取方法目前難以提取出高精度的時(shí)變地震子波。

1.2 分段處理方法

在地震資料處理過(guò)程中,處理非平穩(wěn)地震記錄最簡(jiǎn)單常用的方法就是分段處理。分段提取時(shí)變子波的思想是將非平穩(wěn)地震記錄劃分為若干段,每段地震記錄視為近似平穩(wěn),再結(jié)合其它方法提取子波。分段處理方法在子波振幅與相位的提取中均有廣泛應(yīng)用。

顯然,此類(lèi)方法的關(guān)鍵是合理分段的有效實(shí)現(xiàn),目前常用的分段方法有均勻分段、重疊分段、自適應(yīng)分段等方法。MARGRAVE等[12]將地層視為均勻粘彈性介質(zhì),在反射系數(shù)序列白噪的假設(shè)條件下,采用矩形窗平滑均勻分段對(duì)地震記錄作Gabor變換。BAAN[13]將非平穩(wěn)地震記錄分解成若干等長(zhǎng)的重疊片段,在子波常相位的假設(shè)下,在每段內(nèi)基于峰度最大值準(zhǔn)則提取子波相位,再將提取的相位延拓至每一采樣點(diǎn),實(shí)現(xiàn)子波相位的提取和校正。高靜懷等[14]結(jié)合地震記錄的Robinson模型及Goupilaud模型,提出了一種地震記錄的變子波模型,該模型的近似數(shù)學(xué)表達(dá)式為:

(1)

此外,基于分段處理的子波提取方法普遍存在以下兩方面問(wèn)題:一方面是分段數(shù)量的選取,若分段數(shù)太少,無(wú)法充分反映子波的時(shí)變特性;若分段數(shù)太多,則每段地震記錄所含的數(shù)據(jù)量太少,無(wú)法精確提取出地震子波。另一方面是分段方法的應(yīng)用,若分段的方法選取不當(dāng),提取出的子波可能會(huì)掩蓋相鄰地層間子波振幅和頻率成分的變化。

1.3 應(yīng)用于補(bǔ)償或分段處理后地震記錄的子波提取方法

非平穩(wěn)地震記錄經(jīng)過(guò)衰減補(bǔ)償或分段處理后,整體或分段內(nèi)可看作近似平穩(wěn),然后采用自相關(guān)法、同態(tài)法、高階統(tǒng)計(jì)量法和常相位旋轉(zhuǎn)法等傳統(tǒng)的基于時(shí)不變假設(shè)的子波提取方法實(shí)現(xiàn)子波的估計(jì)。

1.3.1 自相關(guān)法

自相關(guān)法假設(shè)反射系數(shù)為白噪聲,則地震記錄的自相關(guān)等于子波自相關(guān),因此利用地震記錄的自相關(guān)即可得到子波的振幅譜。但真實(shí)的反射系數(shù)序列基本無(wú)法滿(mǎn)足白噪的假設(shè),且二階統(tǒng)計(jì)量不包含子波的相位信息,故自相關(guān)法僅能估計(jì)最小相位子波,方法的應(yīng)用受到較大限制。

1.3.2 同態(tài)法

同態(tài)法根據(jù)子波的復(fù)賽譜(對(duì)數(shù)譜)分布在原點(diǎn)附近,而反射系數(shù)序列的復(fù)賽譜遠(yuǎn)離原點(diǎn)的特性,將子波與反射系數(shù)序列進(jìn)行分離,達(dá)到提取地震子波的目的。劉春成等[16]通過(guò)調(diào)節(jié)改進(jìn)廣義S變換參數(shù),獲得衰減地震記錄的時(shí)頻譜,將其變換到復(fù)賽譜時(shí)頻域提取子波振幅譜。

1.3.3 高階統(tǒng)計(jì)量方法

假設(shè)反射系數(shù)序列是一個(gè)獨(dú)立的、同分布的、非高斯的隨機(jī)過(guò)程,則地震數(shù)據(jù)的高階累積量與地震子波的高階累積量?jī)H相差一個(gè)比例系數(shù),通過(guò)估算地震數(shù)據(jù)的高階累積量即可實(shí)現(xiàn)地震子波的提取[17]。戴永壽等[18]先通過(guò)時(shí)變維納濾波法將振幅譜的影響從地震記錄中消除,將處理后僅殘留相位部分的地震記錄視為近似平穩(wěn),然后用高階統(tǒng)計(jì)量的雙譜法重構(gòu)子波相位譜。

1.3.4 常相位旋轉(zhuǎn)方法

常相位旋轉(zhuǎn)方法假設(shè)子波相位是與頻率無(wú)關(guān)的常數(shù),利用希爾伯特變換對(duì)地震記錄進(jìn)行相位掃描和校正,當(dāng)評(píng)價(jià)準(zhǔn)則達(dá)到最大值時(shí)即估計(jì)得到子波相位。通常使用的評(píng)價(jià)準(zhǔn)則有峰度最大值準(zhǔn)則[13]與偏斜度最大值準(zhǔn)則[19]。劉俊州等[20]通過(guò)數(shù)值模擬試驗(yàn)說(shuō)明偏斜度準(zhǔn)則函數(shù)在動(dòng)態(tài)范圍、穩(wěn)定性、抗噪效果等多方面都要優(yōu)于峰度準(zhǔn)則函數(shù),但當(dāng)噪聲為一定頻率范圍內(nèi)有色噪聲時(shí),二者在相位表現(xiàn)方面都出現(xiàn)較大的偏差。WANG等[21]提出了一種基于Rényi散度的評(píng)價(jià)準(zhǔn)則,相較于峰度最大值準(zhǔn)則,在低頻時(shí)相位估計(jì)更加準(zhǔn)確。

現(xiàn)有的衰減補(bǔ)償方法很容易出現(xiàn)欠補(bǔ)償或過(guò)補(bǔ)償?shù)那闆r,難以保證地震記錄的真實(shí)性;而分段處理方法又無(wú)法很好的體現(xiàn)地震記錄的非平穩(wěn)性,故目前基于高頻補(bǔ)償和分段處理的子波提取方法的實(shí)際應(yīng)用效果仍有待改進(jìn)。

2 基于振幅相位匹配融合的時(shí)變地震子波提取方法

分別估計(jì)子波振幅譜和相位譜再進(jìn)行匹配融合能夠在無(wú)須衰減補(bǔ)償和分段處理的情況下直接從非平穩(wěn)地震記錄中提取具有時(shí)變特征的地震子波。接下來(lái)將分別總結(jié)子波振幅譜和相位譜提取方法的研究進(jìn)展。

2.1 振幅譜提取方法

對(duì)非平穩(wěn)地震記錄進(jìn)行時(shí)頻變換后再應(yīng)用譜模擬方法擬合振幅譜是目前提取時(shí)變子波振幅譜的有效方法,故本節(jié)將對(duì)現(xiàn)有的時(shí)頻分析和譜模擬方法進(jìn)行總結(jié)。

2.1.1 時(shí)頻分析方法

地震信號(hào)屬于非平穩(wěn)信號(hào),常規(guī)傅里葉變換方法僅能分析信號(hào)的頻域特性,無(wú)法準(zhǔn)確描述頻率分布隨時(shí)間的變化情況。而時(shí)頻分析方法能夠?qū)⒁痪S時(shí)域信號(hào)變換到二維的時(shí)頻平面,全面反映地震信號(hào)的時(shí)頻聯(lián)合特征[22],且無(wú)需任何假設(shè)條件地適用于非平穩(wěn)信號(hào)[23],在非平穩(wěn)地震記錄處理的各個(gè)環(huán)節(jié)中發(fā)揮了重大作用。

經(jīng)典的時(shí)頻分析方法主要有短時(shí)傅里葉變換[24]、小波變換[25]、S變換[26]和廣義S變換[27]等,這些方法通過(guò)構(gòu)造各自的窗函數(shù)以截取信號(hào),并在沿時(shí)間方向移動(dòng)的時(shí)窗內(nèi)對(duì)截取的信號(hào)作傅里葉變換對(duì)信號(hào)進(jìn)行時(shí)頻分析,窗函數(shù)的選擇導(dǎo)致這些方法存在各自的局限性:短時(shí)傅里葉變換的窗函數(shù)固定,無(wú)法兼顧時(shí)間分辨率與頻率分辨率;小波變換是一種時(shí)間-尺度分析方法,尺度因子與頻率的關(guān)系不確定,導(dǎo)致應(yīng)用范圍受到一定限制;S變換窗函數(shù)以固定的趨勢(shì)隨頻率變化,不能根據(jù)具體的應(yīng)用而調(diào)整,缺乏靈活性;而廣義S變換受測(cè)不準(zhǔn)原理的制約,無(wú)法同時(shí)達(dá)到最佳的時(shí)間分辨率和頻率分辨率。為了解決這些問(wèn)題,國(guó)內(nèi)外學(xué)者又不斷研發(fā)了更多時(shí)頻分析方法。HUANG等[28]提出了經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)方法,將一個(gè)具有復(fù)雜頻率成分的多分量信號(hào)自適應(yīng)地按照頻率的高低分解為一系列分量之和,這些分量稱(chēng)為為本征模態(tài)函數(shù)(IMF),即:

(2)

式中:x(t)為目標(biāo)信號(hào);ci為第i個(gè)IMF分量;rn為剩余分量。而后,WU等[29]再次提出了一種更為穩(wěn)定的集成經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)方法,該方法通過(guò)多次向信號(hào)中添加固定百分比的高斯白噪聲,再對(duì)經(jīng)過(guò)EMD處理后的數(shù)據(jù)進(jìn)行多次平均處理,在一定程度上解決了模態(tài)混疊效應(yīng),但是仍存在不能保證完全重構(gòu)信號(hào)以及分解出的IMF不唯一等缺點(diǎn)。YEH等[30]通過(guò)在EMD分解信號(hào)過(guò)程中的每一階段添加特定的白噪聲,同時(shí)計(jì)算一個(gè)唯一殘差來(lái)獲得每個(gè)符合定義的IMF,將EMD發(fā)展為互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解(CEEMD),解決了EMD中的模態(tài)混疊問(wèn)題,并實(shí)現(xiàn)了原始信號(hào)精確重構(gòu)。DAUBECHIES等[31]提出了應(yīng)用于時(shí)頻域信號(hào)的同步壓縮變換,該方法通過(guò)提取時(shí)頻譜的瞬時(shí)頻率信息,對(duì)時(shí)頻平面的能量在頻率方向進(jìn)行重新分配,極大地提高了頻率方向的清晰度。WANG等[32]開(kāi)發(fā)了廣義β小波族(GBWs)作為連續(xù)小波變換的母小波,通過(guò)改變控制小波形狀的兩個(gè)參數(shù)(α和β),實(shí)現(xiàn)自適應(yīng)地權(quán)衡時(shí)間分辨率與頻率分辨率。姚振岸等[33]首先對(duì)非平穩(wěn)地震道作基追蹤譜分解,將地震信號(hào)分解到預(yù)先定義的波形字典的各個(gè)基函數(shù)上,然后求得基函數(shù)表示系數(shù),經(jīng)尺度到頻率的映射,得到高分辨率譜分解結(jié)果,再逐個(gè)時(shí)間點(diǎn)提取地震子波,該方法能夠同時(shí)獲得較高的時(shí)間和頻率分辨率。時(shí)頻分析方法將時(shí)間域的信號(hào)映射到時(shí)頻域,在二維時(shí)頻譜中進(jìn)行相關(guān)處理,能夠較為準(zhǔn)確地同時(shí)定位時(shí)間和頻率信息,是時(shí)變子波提取的有力工具。

2.1.2 譜模擬方法

譜模擬方法的主要思想是對(duì)地震子波振幅譜進(jìn)行參數(shù)化建模,從地震記錄的振幅譜中擬合出子波的振幅譜[34]。該方法最早由ROSA等[35]提出,假設(shè)子波振幅譜是類(lèi)似于雷克子波的光滑單峰曲線,在頻率域建立數(shù)學(xué)模型:

(3)

式中:k為常數(shù);N為階數(shù);an為關(guān)于f的多項(xiàng)式的系數(shù)。一般地,0

尚新民等[1]采用改進(jìn)S變換對(duì)地震記錄進(jìn)行時(shí)頻分析,然后利用(3)式對(duì)不同時(shí)刻地震記錄的振幅譜進(jìn)行最小二乘擬合實(shí)現(xiàn)時(shí)變子波振幅譜的提取。譜模擬方法能夠突破反射系數(shù)序列白噪的假設(shè),原理簡(jiǎn)單、效果明顯,在子波振幅譜的提取中得到廣泛應(yīng)用。李振春等[36]通過(guò)改變譜模擬參數(shù)擬合出一簇地震子波振幅譜,然后計(jì)算各擬合子波的自相關(guān)與從地震記錄中估算出的子波自相關(guān)之間的誤差,誤差最小時(shí)即確定為最優(yōu)的子波振幅譜,提高了譜模擬方法的適應(yīng)性和準(zhǔn)確性。MAMASANI等[37]通過(guò)可分非線性最小二乘算法找出多個(gè)不同峰值頻率的最優(yōu)Ricker子波以重建子波振幅譜,避免了多項(xiàng)式譜模擬中參數(shù)的選擇,求解效率較高。唐博文等[38]建立了二次譜模擬方法的基本原理,假設(shè)地震子波的振幅譜是光滑的,且與反射系數(shù)的二次譜可分,則可設(shè)計(jì)一個(gè)低通濾波器作用于地震記錄的二次譜,濾除高頻成分,重構(gòu)低頻成分,實(shí)現(xiàn)子波振幅譜的估計(jì)。

目前譜模擬方法是子波振幅譜提取最常用的方法,但此類(lèi)方法要求子波振幅譜滿(mǎn)足類(lèi)似于雷克子波的單峰、光滑的曲線的假設(shè),限制了子波振幅譜的形態(tài),且多項(xiàng)式中的參數(shù)N難以選擇;基于優(yōu)勢(shì)Ricker分量的譜模擬方法在子波振幅譜形態(tài)不滿(mǎn)足假設(shè)條件時(shí),估計(jì)誤差較大;二次譜模擬方法構(gòu)造的低通濾波器的截止頻率等參數(shù)只能依靠人工設(shè)定。

2.2 相位譜提取方法

為了解決分段方法受限于子波分段平穩(wěn)假設(shè)的問(wèn)題,局部相似度和純相位濾波器尋優(yōu)等方法常被應(yīng)用于子波相位估計(jì)中。

2.2.1 局部相似度方法

局部相似度方法的主要思想是先將地震道進(jìn)行相位旋轉(zhuǎn),計(jì)算旋轉(zhuǎn)后的地震道與標(biāo)準(zhǔn)道的相似度,相似度達(dá)到最大時(shí)的旋轉(zhuǎn)相位即為最優(yōu)相位。BAAN等[39]和FOMEL等[40]提出將相位估計(jì)看成最小二乘反演問(wèn)題,將地震道的包絡(luò)看作標(biāo)準(zhǔn)道,不斷改變地震數(shù)據(jù)的相位,利用相位與地震記錄包絡(luò)的局部相似度作為判別準(zhǔn)則來(lái)估計(jì)子波相位。DAI等[41]對(duì)這一類(lèi)方法做了進(jìn)一步的改進(jìn),首先利用高階累積量的雙譜估計(jì)相位范圍,然后在預(yù)估范圍內(nèi)進(jìn)行常相位旋轉(zhuǎn),并利用局部相似性?xún)?yōu)化法準(zhǔn)確提取時(shí)變子波相位譜,解決了單獨(dú)使用雙譜估計(jì)相位不準(zhǔn)確以及局部相似性?xún)?yōu)化方法計(jì)算量過(guò)大的問(wèn)題。該方法以常相位旋轉(zhuǎn)法為基礎(chǔ),須假設(shè)子波為常相位或不同時(shí)刻子波相位的變化量為常數(shù),存在一定的局限性。

2.2.2 純相位濾波器尋優(yōu)方法

純相位濾波器尋優(yōu)方法的思想是在子波振幅譜已知的條件下,通過(guò)改變子波Z變換零極點(diǎn)的分布情況,構(gòu)造由一系列相同振幅譜不同相位譜組成的子波庫(kù),再通過(guò)優(yōu)化算法從中選擇使得準(zhǔn)則函數(shù)達(dá)到極值的子波。張亞南等[42]采用單位化自回歸滑動(dòng)平均(ARMA)模型振幅譜來(lái)構(gòu)造純相位濾波器,并在最大方差模準(zhǔn)則約束下,采用改進(jìn)的粒子群算法對(duì)子波相位殘余進(jìn)行非線性尋優(yōu)。WANG等[43]提出了基于蟻群算法的子波相位估計(jì)方法,首先對(duì)子波振幅譜進(jìn)行希爾伯特變換估計(jì)最小相位子波,然后將最小相位子波Z域的零點(diǎn)關(guān)于單位圓進(jìn)行對(duì)稱(chēng)變換,建立候選子波庫(kù),最后利用蟻群算法實(shí)現(xiàn)子波相位的準(zhǔn)確尋優(yōu)。該方法在有色噪聲的影響下具有較好的穩(wěn)定性,但蟻群、粒子群算法等在確定尋優(yōu)初值方面還有待進(jìn)一步研究。

2.2.3 其它相位估計(jì)方法

VOSOUGHI等[44]將混合相位子波視為最小相位子波與全通子波的褶積,首先從地震記錄中求出最小相位子波,再使用累積量匹配法從譜白化后的地震數(shù)據(jù)中得到全通子波,最后將二者褶積得到混合相位子波。此方法能夠增加數(shù)據(jù)帶寬,且算法簡(jiǎn)單、計(jì)算速度快,但對(duì)數(shù)據(jù)帶寬與中心頻率之比有嚴(yán)格要求。ZHANG等[45]根據(jù)粘彈性介質(zhì)中地震波的傳播規(guī)律,以地層品質(zhì)因子Q為中間橋梁,推導(dǎo)建立了子波振幅譜和相位譜的關(guān)系,從而利用所估計(jì)的子波振幅譜準(zhǔn)確估計(jì)時(shí)變子波相位譜。

目前基于振幅相位匹配融合的時(shí)變地震子波提取方法通常存在限制子波振幅譜形態(tài),子波相位局限于分段平穩(wěn)或常相位假設(shè)的問(wèn)題,但此類(lèi)方法相較于基于衰減補(bǔ)償?shù)淖硬ㄌ崛》椒?提取的子波更為準(zhǔn)確,是今后研究的重點(diǎn)方向之一。

3 總結(jié)與展望

本文分析了主流時(shí)變子波提取方法的研究進(jìn)展,詳細(xì)總結(jié)了現(xiàn)有方法的優(yōu)缺點(diǎn)及適用條件,從而得出結(jié)論:高頻衰減補(bǔ)償與分段處理方法嚴(yán)重依賴(lài)于地層衰減的估計(jì)精度,分段數(shù)目與分段方法的選擇等因素,因而應(yīng)用范圍受限。譜模擬等時(shí)變子波振幅譜提取方法限制了子波振幅譜的形態(tài),當(dāng)子波振幅不滿(mǎn)足假設(shè)條件時(shí),提取效果不佳?，F(xiàn)有時(shí)變子波相位譜提取方法大多局限于子波常相位或分段平穩(wěn)假設(shè),難以應(yīng)用于復(fù)雜地層衰減、強(qiáng)噪聲干擾的地區(qū)。如何改進(jìn)現(xiàn)有方法從而有效提高子波估計(jì)精度是目前亟待解決的問(wèn)題。同時(shí),隨著地震勘探技術(shù)的不斷發(fā)展,研究目標(biāo)已經(jīng)逐漸轉(zhuǎn)向深層復(fù)雜構(gòu)造的準(zhǔn)確成像,對(duì)地震子波提取以及高分辨率地震資料處理的方法研究也提出了更高的要求。

綜上所述,目前時(shí)變子波提取方法的實(shí)現(xiàn)普遍建立在地下介質(zhì)線性衰減、噪聲高斯分布等假設(shè)條件下,但實(shí)際情況往往更加復(fù)雜,非線性衰減與非高斯噪聲干擾等會(huì)對(duì)地震子波的準(zhǔn)確提取產(chǎn)生嚴(yán)重的影響。同時(shí),因深層相對(duì)于中淺層具有信噪比更低、傳播路徑更復(fù)雜、衰減更嚴(yán)重等特點(diǎn),現(xiàn)有子波提取方法難以直接應(yīng)用于深層地震資料處理。此外,深層復(fù)雜構(gòu)造的研究對(duì)勘探精度有著更高的要求,利用疊后地震數(shù)據(jù)提取出的子波無(wú)法完全滿(mǎn)足全波形反演和偏移成像等處理環(huán)節(jié)的實(shí)際需求,而利用包含更多信息的疊前地震數(shù)據(jù)提取子波能夠更可靠地揭示地下儲(chǔ)層的分布情況。因此今后的研究工作將主要從以下4方面展開(kāi)。

1) 將現(xiàn)有時(shí)變子波提取方法擴(kuò)展應(yīng)用至非線性衰減、非高斯噪聲干擾條件下,研究品質(zhì)因子Q隨頻率變化時(shí)的子波提取方法。

2) 進(jìn)一步關(guān)注智能信息處理技術(shù)的前沿發(fā)展,研究適用范圍更廣的時(shí)變地震子波提取方法。

3) 針對(duì)部分地區(qū)深層油氣勘探面臨的信號(hào)弱、分辨能力不足等問(wèn)題,開(kāi)展強(qiáng)衰減、低信噪比條件下的子波提取方法研究,并利用提取的子波對(duì)地震數(shù)據(jù)進(jìn)行高低頻補(bǔ)償,拓寬頻帶,提高地震資料的分辨率。

4) 開(kāi)展疊前地震數(shù)據(jù)中子波提取方法的研究,并應(yīng)用于全波形反演和疊前偏移,有效提高速度建模與成像精度。