基于高階累積量時(shí)間延遲估計(jì)的自動(dòng)剩余靜校技術(shù)

安圣培,劉 韜,胡天躍

(1.中國(guó)石油化工股份有限公司石油勘探開(kāi)發(fā)研究院,北京100083;2.北京大學(xué)地球與空間科學(xué)學(xué)院,北京100871)

目前,中國(guó)西部地區(qū)已經(jīng)成為油氣勘探的重點(diǎn)地區(qū)之一,該地區(qū)地表?xiàng)l件相對(duì)復(fù)雜,如山前帶,戈壁和沙漠,近地表非成巖低速介質(zhì)區(qū)往往存在近地表高程、厚度和速度的空間變化,靜校正問(wèn)題突出,從而影響疊加成像效果以及地表風(fēng)化層速度異常的檢測(cè)。

折射類(lèi)方法在解決靜校正問(wèn)題中發(fā)揮了重要作用,如用于基準(zhǔn)面靜校正的折射靜校正方法[1-2]和層析靜校正方法[3-4],以及基于統(tǒng)計(jì)相關(guān)實(shí)現(xiàn)的剩余靜校正方法[5-6]。這些方法的有效性都依賴(lài)于初至拾取的準(zhǔn)確度,但復(fù)雜地表?xiàng)l件往往伴隨著強(qiáng)背景噪聲,降低了初至的信噪比,導(dǎo)致傳統(tǒng)的自動(dòng)初至拾取的準(zhǔn)確度下降,而手動(dòng)初至拾取在面對(duì)大量地震數(shù)據(jù)時(shí)會(huì)產(chǎn)生極大的人力和時(shí)間成本。針對(duì)這一問(wèn)題,有學(xué)者提出使用濾波類(lèi)方法或地震干涉類(lèi)方法提高初至信噪比,以保證自動(dòng)拾取初至的準(zhǔn)確性[7-9]。

在地震干涉類(lèi)方法中,MIKESELL等[10]和ZHANG等[11]使用基于互相關(guān)的地震干涉法,計(jì)算兩點(diǎn)間的時(shí)延,并使用時(shí)延信息進(jìn)行近地表速度建模,再基于速度模型計(jì)算剩余靜校正量。時(shí)延估計(jì)使用了互相關(guān)的方法,但互相關(guān)容易受高斯噪聲的影響,尤其是相干的高斯噪聲[12]。為解決這一問(wèn)題,AN等[8]提出使用高階累積量替代互相關(guān)以提高地震干涉法的抗噪性,并由此發(fā)展了累積量相干積累法用于加強(qiáng)初至波信號(hào)。高斯噪聲的高階累積量為0,可有效提取高斯噪聲中的非高斯信號(hào),更加符合實(shí)際地震數(shù)據(jù)的情況。高階累積量同時(shí)提供信號(hào)的幅度和相位信息,可用于描述信號(hào)系統(tǒng)的非線(xiàn)性特征?；谝陨咸攸c(diǎn),高階累積量已經(jīng)被廣泛應(yīng)用于子波特征分析、噪聲衰減以及初至拾取等地震數(shù)據(jù)處理領(lǐng)域[13-15]。

在上述研究的基礎(chǔ)上,本文提出了基于高階累積量時(shí)間延遲估計(jì)的自動(dòng)初至剩余靜校正方法,該方法使用高階累積量和相干疊加的方法,估計(jì)兩點(diǎn)間的初至波時(shí)延,再使用不同時(shí)延信息的組合得到炮檢點(diǎn)的剩余靜校正量。本文方法的優(yōu)勢(shì)在于:①互相關(guān)使用兩道數(shù)據(jù)進(jìn)行時(shí)延估計(jì),只能利用二維地震信息,而高階累積量使用多道數(shù)據(jù),可以利用三維地震數(shù)據(jù)信息增強(qiáng)有效信號(hào);②互相關(guān)是基于噪聲為高斯分布的假設(shè),但實(shí)際地震數(shù)據(jù)很可能不滿(mǎn)足這一假設(shè),導(dǎo)致道間時(shí)延估計(jì)的準(zhǔn)確度不高,而高階累積量并不要求噪聲為高斯分布,更加符合實(shí)際地震數(shù)據(jù)的情況,對(duì)背景噪聲的壓制效果更優(yōu);③再結(jié)合地震干涉法中相干疊加的思路進(jìn)一步提高信號(hào)的有效疊加次數(shù),增強(qiáng)信噪比,保證本文方法獲得穩(wěn)定且可靠的道間時(shí)延估計(jì)結(jié)果;④本文方法采用完全數(shù)據(jù)驅(qū)動(dòng),相比基于初至拾取的剩余靜校正方法,在面對(duì)海量低信噪比地震數(shù)據(jù)時(shí),可以有效降低初至拾取所需的高昂時(shí)間和人力成本。

1 方法原理

本文方法主要包括兩個(gè)部分:①使用高階累積量估計(jì)兩點(diǎn)間的初至波時(shí)延,并使用相干疊加的方法增強(qiáng)信噪比,保證時(shí)延估計(jì)的準(zhǔn)確性;②使用時(shí)延信息的組合得到炮檢點(diǎn)的剩余靜校正量。

1.1 高階累積量時(shí)延估計(jì)

本文方法首先使用高階累積量估計(jì)兩點(diǎn)間的時(shí)延。二階統(tǒng)計(jì)工具,如互相關(guān),常被用于時(shí)延估計(jì),但互相關(guān)對(duì)高斯噪聲敏感,因而人們開(kāi)始使用高階統(tǒng)計(jì)工具進(jìn)行時(shí)延估計(jì),如高階累積量。在實(shí)際問(wèn)題中,加性噪聲往往具有高斯性,因此,使用高階累積量進(jìn)行時(shí)延估計(jì)更加合理,因?yàn)楦咚惯^(guò)程的高階累積量為0[16]。對(duì)零均值的平穩(wěn)隨機(jī)過(guò)程{x0(t),x1(t),…,xk-1(t)},其k階累積量由{x0(t),x1(t+τ1),…,xk-1(t+τk-1)}的聯(lián)合k階累積量表示如下[16]:

Ck(τ1,τ2,…,τk-1)=Cum{x0(t),x1(t+τ1),…,xk-1(t+τk-1)}

(1)

其中,Cum{·}是累積量算子,Ck(·)表示k階累積量函數(shù),τ1,τ2,…,τk-1分別代表x0與x1,x2,…,xk-1之間的時(shí)延。常用的二階、三階和四階累積量函數(shù)的表達(dá)式如下[16]:

C2(τ)=E{x0(t),x1(t+τ)}

(2)

C3(τ1,τ2)=E{x0(t),x1(t+τ1),x2(t+τ2)}

(3)

C4(τ1,τ2,τ3)=E{x0(t),x1(t+τ1),x2(t+τ2),

x3(t+τ3)}-C2(τ1)C2(τ2-τ3)-C2(τ2)·

C2(τ3-τ1)-C2(τ3)C2(τ1-τ2)

(4)

其中,E{·}是期望算子。公式(2)表示協(xié)方差,在零均值的假設(shè)下等價(jià)于互相關(guān)。

基于高階累積量進(jìn)行兩點(diǎn)間初至波時(shí)延估計(jì)的具體物理過(guò)程如下。如圖1所示,對(duì)于炮點(diǎn)S激發(fā),檢波點(diǎn)Rj接收的地震記錄Uj(t),假設(shè)它是由初至波Gj(t)和高斯隨機(jī)噪聲w(t)組成的,其表達(dá)式為:

Uj(t)=Gj(t)+w(t)

(5)

使用炮點(diǎn)S在檢波點(diǎn)Rj,Rj+1,…,Rj+n接收的n+1道地震記錄,計(jì)算其對(duì)應(yīng)的n+1階累積量的公式如下:

Cj;j+1,j+2,…,j+n(τ1,τ2,…,τn)=Cum{Uj(t),

Uj+1(t+τ1),…,Uj+n(t+τn)}=Cum{Gj(t),

Gj+1(t+τ1),…,Gj+n(t+τn)}

(6)

其中,在階數(shù)≥3的高階累積量中,w(t)的響應(yīng)為0。n+1階累積量包含n個(gè)時(shí)延τp(p=1,2,…,n),τp表示檢波點(diǎn)Rj與Rj+p之間的初至波時(shí)延。圖1中的V代表虛源位置,所謂虛源指的是假想的震源,并非真實(shí)存在的震源。S到Rj的旅行時(shí)與S到Rj+1的旅行時(shí)共用了S到V的旅行時(shí),對(duì)應(yīng)圖1中的灰色粗線(xiàn),共用旅行時(shí)在估計(jì)時(shí)延的過(guò)程中被相互抵消,因此,τp只反映了兩個(gè)檢波點(diǎn)之間的時(shí)延,對(duì)應(yīng)圖1中的紅線(xiàn)。高階累積量方法使用了n+1道信號(hào)進(jìn)行時(shí)延估計(jì),而互相關(guān)方法只使用兩道進(jìn)行時(shí)延估計(jì),因此,高階累積量方法的抗噪性更好。n+1階累積量Cj;j+1,j+2,…,j+n(τ1,τ2,…,τn)是n維函數(shù),理論上可以直接求取n個(gè)時(shí)延,但隨著n的增大,存儲(chǔ)量和運(yùn)算量都會(huì)顯著增加。因此,本文方法每次只計(jì)算一個(gè)時(shí)延,實(shí)現(xiàn)方式如下:

圖1 使用高階累積量估計(jì)兩點(diǎn)間時(shí)延示意

(7)

其中,τ是第j+1個(gè)檢波點(diǎn)與第j個(gè)檢波點(diǎn)之間的初至波時(shí)延。公式(7)只以?xún)傻赖卣鹩涗洖檩斎?理論上仍可以得到n階累積量,但由于Gj(t)與自身的時(shí)延為0,因此,不再需要計(jì)算累積量中的除τ以外的其它時(shí)延,大大節(jié)省了運(yùn)算量,并且仍然利用了高階累積量的抗噪性?xún)?yōu)勢(shì)。

1.2 相干疊加

為了保證兩點(diǎn)間時(shí)延估計(jì)的準(zhǔn)確性,本文進(jìn)一步使用相干疊加方法,增強(qiáng)兩點(diǎn)間初至波信噪比,具體過(guò)程如下。如圖2所示,對(duì)炮點(diǎn)S1在檢波點(diǎn)R1和R2的初至波記錄進(jìn)行高階累積量計(jì)算,共用的S1到虛源V的旅行時(shí)在時(shí)延估計(jì)過(guò)程中被抵消,因此,得到了R1和R2間的初至波響應(yīng)(圖2b中紅線(xiàn)所示),圖中實(shí)線(xiàn)表示旅行時(shí)為正,虛線(xiàn)表示旅行時(shí)為負(fù)。使用其它炮點(diǎn)Si(i=2,3,…,N)重復(fù)上述過(guò)程,炮點(diǎn)Si到V的旅行時(shí)仍共用,在時(shí)延估計(jì)時(shí)被抵消,因此,使用不同炮點(diǎn)得到的高階累積量結(jié)果中包含了相同的R1到R2的道間時(shí)差(圖2b)。對(duì)圖2b中的所有結(jié)果進(jìn)行疊加,實(shí)現(xiàn)了初至波響應(yīng)的同相位疊加,而非相干的隨機(jī)噪聲由于具有不同的時(shí)延,不能實(shí)現(xiàn)同相位疊加,疊加能量將小于初至波。至此,實(shí)現(xiàn)了使用相干疊加增強(qiáng)兩點(diǎn)間初至波響應(yīng)的信噪比。

圖2 相干疊加示意a 對(duì)炮點(diǎn)在檢波點(diǎn)R1和R2的初至波記錄進(jìn)行高階累積量計(jì)算; b 不同炮點(diǎn)得到的高階累積量結(jié)果; c 圖2b的疊加結(jié)果

1.3 時(shí)延組合求取剩余靜校正量

(8)

(9)

(10)

(11)

圖3 使用時(shí)延信息求取剩余靜校正量示意a 使用檢波點(diǎn)左側(cè)炮點(diǎn)的地震記錄進(jìn)行檢波點(diǎn)Rj和Rj+1間的時(shí)延估計(jì); b 使用檢波點(diǎn)右側(cè)炮點(diǎn)的地震記錄進(jìn)行檢波點(diǎn)Rj+1和Rj間的時(shí)延估計(jì); c 使用a和b中得到的時(shí)延信息進(jìn)行相減得到相鄰檢波點(diǎn)間的相對(duì)剩余靜校正量

ΔTj,j+1和ΔTj+1,j都包含了傳播至Rj和Rj+1的兩組上行旅行時(shí)和沿高速層頂傳播的旅行時(shí)。假設(shè)地表高程和近地表低速層厚度引起的長(zhǎng)波長(zhǎng)靜校正量已經(jīng)被消除,則ΔTj,j+1和ΔTj+1,j包含的上、下行旅行時(shí)相同;同時(shí)高速層頂界面的起伏在長(zhǎng)波長(zhǎng)靜校正中被修正,ΔTj,j+1和ΔTj+1,j包含的沿高速層頂?shù)膫鞑ゾ嚯x也近似相等,由于Rj和Rj+1為相鄰檢波點(diǎn),在相對(duì)小的空間范圍內(nèi),可以認(rèn)為高速層頂?shù)乃俣茸兓淮?因此ΔTj,j+1和ΔTj+1,j包含的沿高速層頂傳播的旅行時(shí)也近似相等,ΔTj,j+1和ΔTj+1,j近似相等,ZHANG等[11]給出過(guò)類(lèi)似的證明。因此,公式(10)和公式(11)相減得到:

(12)

1.4 方法流程

本文方法的實(shí)現(xiàn)流程如圖4所示。如果數(shù)據(jù)存在長(zhǎng)波長(zhǎng)靜校正問(wèn)題,需要先進(jìn)行高程靜校正或野外靜校正;針對(duì)數(shù)據(jù)可能存在的異常高頻噪聲,使用高通濾波進(jìn)行噪聲壓制;使用高階累積量和相干疊加方法得到兩點(diǎn)間的初至?xí)r延;最后,使用時(shí)延信息的組合得到炮檢點(diǎn)的剩余靜校正量。

圖4 本文方法流程

2 合成數(shù)據(jù)實(shí)例

2.1 高階累積量時(shí)延估計(jì)實(shí)例

本節(jié)通過(guò)合成數(shù)據(jù)實(shí)例,對(duì)比了高階累積量和互相關(guān)方法應(yīng)用于時(shí)延估計(jì)的抗噪性。圖5a和圖5b展示了僅包含一組折射波信號(hào)的兩個(gè)單炮記錄,折射波視速度為2500m/s,檢波點(diǎn)間距為15m,兩個(gè)炮點(diǎn)相距60m,雷克子波主頻為60Hz。圖5c和圖5d分別為對(duì)圖5a和圖5b加入相干高斯噪聲的單炮記錄,其信噪比為-9dB。圖5a和圖5b的互相關(guān)結(jié)果如圖5e所示。圖5f、圖5g和圖5h分別為對(duì)含噪數(shù)據(jù)(圖5c和圖5d)使用互相關(guān)、三階累積量和四階累積量得到的結(jié)果。由于互相關(guān)對(duì)相干噪聲敏感,圖5f的互相關(guān)結(jié)果中除了包含折射波的響應(yīng),還混入了相干噪聲的響應(yīng),信噪比為3dB,圖5g和圖5h中的相干噪聲被有效壓制,殘余噪聲也不具有相同的時(shí)延,因此無(wú)法實(shí)現(xiàn)相干疊加,其信噪比分別達(dá)到了9dB和15dB。圖5i、圖5j、圖5k和圖5l分別為對(duì)圖5e、圖5f、圖5g和圖5h的結(jié)果沿檢波點(diǎn)疊加得到的疊加單道。圖5e使用不同檢波點(diǎn)得到的折射波響應(yīng)具有相同的時(shí)延,因此在疊加過(guò)程中實(shí)現(xiàn)了相干疊加,疊加結(jié)果(圖5i)中的峰值位置與理論時(shí)延(紅色虛線(xiàn)標(biāo)示)一致。圖5j中相干噪聲響應(yīng)的峰值大于折射波響應(yīng)的峰值,因此難以準(zhǔn)確得到折射波響應(yīng)的時(shí)延估計(jì)。圖5k和圖5l中包含的折射波響應(yīng)的峰值位置與理論時(shí)延一致。該實(shí)例證明了三階和四階的高階累積量方法壓制相干白噪聲的效果優(yōu)于互相關(guān)方法,同時(shí),高階累積量的階數(shù)越高,得到的折射波響應(yīng)的信噪比也越高。

圖5 合成數(shù)據(jù)的測(cè)試結(jié)果a,b 僅包含一組折射波信號(hào)的單炮記錄; c,d 加入相干高斯噪聲的單炮記錄; e 圖5a和圖5b的互相關(guān)結(jié)果; f 圖5c和圖5d的互相關(guān)結(jié)果; g 圖5c和圖5d的三階累積量結(jié)果; h 圖5c和圖5d的四階累積量結(jié)果; i 對(duì)圖5e疊加得到的疊加單道; j 對(duì)圖5f疊加得到的疊加單道; k 對(duì)圖5g 疊加得到的疊加單道; l 對(duì)圖5h疊加得到的疊加單道

2.2 求取剩余靜校正量實(shí)例

采用聲波有限差分正演得到合成數(shù)據(jù),所用模型參數(shù)見(jiàn)表1,炮距為40m,道距為10m,炮點(diǎn)數(shù)100,檢波點(diǎn)數(shù)400。剩余靜校正量由隨機(jī)序列產(chǎn)生,滿(mǎn)足均值為0和標(biāo)準(zhǔn)差為8ms的正太分布,正演得到的單炮記錄如圖6a所示,由于剩余靜校正量的存在,圖6a中同相軸的連續(xù)性較差。將采用本文方法得到的剩余靜校正量應(yīng)用到單炮記錄中,得到的結(jié)果如圖6b 所示。由圖6b可見(jiàn),單炮記錄同相軸的連續(xù)性顯著提高。圖7給出了對(duì)不含噪合成數(shù)據(jù)采用本文方法得到的剩余靜校正量。其中,灰線(xiàn)表示真實(shí)靜校正量,紅線(xiàn)表示計(jì)算得到的靜校正量。由圖7可見(jiàn),本文方法計(jì)算得到的剩余靜校正量與真實(shí)值完全吻合。

圖6 對(duì)不含噪合成數(shù)據(jù)(a)應(yīng)用本文方法進(jìn)行剩余靜校正處理后的結(jié)果(b)

圖7 對(duì)不含噪合成數(shù)據(jù)應(yīng)用本文方法得到的剩余靜校正量

表1 模型參數(shù)

對(duì)圖6a數(shù)據(jù)加入相干高斯噪聲,得到信噪比為-7dB的含噪數(shù)據(jù)(圖8a)。將采用本文方法得到的剩余靜校正量應(yīng)用到含噪單炮記錄中,得到如圖8b所示的結(jié)果。由圖8b可見(jiàn),反射波同相軸連續(xù)性明顯提高。圖8c和圖8d分別是圖8a和圖8b方框標(biāo)示區(qū)域的局部放大圖。圖9展示了對(duì)含噪合成數(shù)據(jù)使用互相關(guān)和本文方法得到的剩余靜校正量。由圖9可以看出,本文方法得到的剩余靜校正量與真實(shí)值基本一致,互相關(guān)方法受相干噪聲影響,其結(jié)果偏離真實(shí)值較多。該實(shí)例驗(yàn)證了本文方法求取剩余靜校正量的有效性,同時(shí),相比互相關(guān)方法,本文方法不易受相干噪聲影響,抗噪性更好。

圖8 對(duì)含噪合成數(shù)據(jù)應(yīng)用本文方法進(jìn)行剩余靜校正處理后的結(jié)果a 未做剩余靜校正的單炮記錄; b 采用本文方法進(jìn)行剩余靜校正處理后的單炮記錄; c 圖8a方框區(qū)域局部放大結(jié)果; d 圖8b方框區(qū)域局部放大結(jié)果

圖9 對(duì)含噪合成數(shù)據(jù)使用互相關(guān)和本文方法得到的剩余靜校正量

3 實(shí)際資料應(yīng)用

選取中國(guó)西部某地區(qū)的實(shí)際地震資料驗(yàn)證本文方法的有效性。該工區(qū)地表存在較強(qiáng)的近地表速度異常,剩余靜校正問(wèn)題突出。對(duì)比了本文方法和某商業(yè)軟件中的基于初至拾取的初至波剩余靜校正模塊的結(jié)果,該模塊采用的是廣義互換的方法。基于初至拾取的剩余靜校正方法要求拾取全部道集的初至到時(shí),但在遠(yuǎn)偏移距處,受該工區(qū)復(fù)雜近地表影響,初至信噪比極低,不得不進(jìn)行耗時(shí)耗力的手動(dòng)初至拾取,而本文方法可以避免初至拾取的工作量。圖10a和圖10b對(duì)比了本文方法和基于手動(dòng)初至拾取方法得到的檢波點(diǎn)和炮點(diǎn)的剩余靜校正量,結(jié)果基本吻合,僅在檢波點(diǎn)號(hào)為900附近位置存在微弱差異。將剩余靜校正量應(yīng)用到道集數(shù)據(jù),選取剩余靜校正問(wèn)題比較突出的原始單炮(圖11a)進(jìn)行比較,該單炮初至波同相軸明顯不連續(xù),圖11b和圖11c分別是應(yīng)用本文方法和基于初至拾取方法的結(jié)果,兩種方法都明顯改善了初至波同相軸的連續(xù)性,圖11d至圖11f分別是圖11a至圖11c黑色方框處的局部放大圖,反射波同相軸的連續(xù)性也明顯改善。圖12對(duì)比了本文方法和基于手動(dòng)初至拾取方法得到的疊加剖面。圖12a為未進(jìn)行剩余靜校正的疊加剖面。由圖12a可見(jiàn),該疊加剖面同相軸連續(xù)性較差。采用本文方法得到的疊加剖面(圖12b)同相軸的連續(xù)性顯著改善,并且與基于手動(dòng)初至拾取的剩余靜校正方法相比,本文方法避免了手動(dòng)初至拾取的繁重工作量。

圖10 本文方法和基于手動(dòng)初至拾取方法應(yīng)用于實(shí)際數(shù)據(jù)得到的剩余靜校正量對(duì)比a 檢波點(diǎn)剩余靜校正量對(duì)比; b 炮點(diǎn)剩余靜校正量對(duì)比

圖11 本文方法和基于手動(dòng)初至拾取方法進(jìn)行剩余靜校正的單炮記錄效果對(duì)比a 未進(jìn)行剩余靜校正的單炮記錄; b 采用本文方法進(jìn)行剩余靜校正后的單炮記錄; c 基于手動(dòng)初至拾取方法進(jìn)行剩余靜校正后的單炮記錄; d 圖11a黑色方框區(qū)域的局部放大結(jié)果; e 圖11b黑色方框區(qū)域的局部放大結(jié)果; f 圖11c黑色方框區(qū)域的局部放大結(jié)果

圖12 本文方法和基于手動(dòng)初至拾取方法進(jìn)行剩余靜校正后的疊加剖面效果對(duì)比a 未進(jìn)行剩余靜校正的疊加剖面; b 采用本文方法進(jìn)行剩余靜校正后的疊加剖面; c 基于手動(dòng)初至拾取方法進(jìn)行剩余靜校正后的疊加剖面

4 結(jié)論與認(rèn)識(shí)

本文發(fā)展了基于高階累積量時(shí)間延遲估計(jì)的自動(dòng)剩余靜校正方法,采用高階累積量得到兩點(diǎn)間的初至波時(shí)延;使用相干疊加提高信噪比,保證時(shí)延估計(jì)的準(zhǔn)確性;再使用時(shí)延信息的組合得到炮檢點(diǎn)的剩余靜校正量。與互相關(guān)類(lèi)方法相比,使用高階累積量進(jìn)行時(shí)延估計(jì)可有效利用三維地震數(shù)據(jù)信息增強(qiáng)有效信號(hào),高階累積量不再要求噪聲為高斯分布,與實(shí)際地震資料的情況更加吻合;采用相干疊加方法提高信號(hào)有效疊加次數(shù)和信噪比,保證本文方法獲得穩(wěn)定且可靠的道間時(shí)延估計(jì)結(jié)果;本文方法采用完全數(shù)據(jù)驅(qū)動(dòng),與基于初至拾取的剩余靜校正方法相比,可有效避免初至拾取所需的高昂時(shí)間和人力成本。模型試驗(yàn)和實(shí)際資料的應(yīng)用結(jié)果表明,本文方法可以明顯改善疊加剖面中反射波同相軸的連續(xù)性,為后續(xù)地震數(shù)據(jù)處理提供高品質(zhì)資料。

本文方法涉及到高階累積量階數(shù)的選取,選取的階數(shù)越高,噪聲壓制效果越好,但運(yùn)算量也同步增加,本文選取了常用的四階累積量,在具體應(yīng)用于實(shí)際地震資料時(shí),可事先通過(guò)參數(shù)測(cè)試確定適合的累積量階數(shù),實(shí)現(xiàn)高效穩(wěn)定的剩余靜校正量求取。另外,本文方法要求觀測(cè)系統(tǒng)規(guī)則化,以實(shí)現(xiàn)兩點(diǎn)間初至波響應(yīng)的相干疊加,針對(duì)不規(guī)則觀測(cè)系統(tǒng),可事先使用空間插值技術(shù)重構(gòu)得到規(guī)則化道集。本文方法實(shí)現(xiàn)了自動(dòng)剩余靜校正,未來(lái)可結(jié)合人工智能類(lèi)方法實(shí)現(xiàn)流程參數(shù)的自動(dòng)優(yōu)化選取,進(jìn)一步提高自動(dòng)化程度以減少人工成本,更加高效地應(yīng)用于實(shí)際工區(qū)的地震資料處理。本文方法可進(jìn)一步與深度學(xué)習(xí)結(jié)合,改進(jìn)傳統(tǒng)的基于特征值的時(shí)延估計(jì)方法,提高道間時(shí)延估計(jì)的準(zhǔn)確率和對(duì)更加復(fù)雜噪聲環(huán)境的適用性,并同時(shí)降低特征提取和篩選的人工成本,提高計(jì)算效率。