一種諧波噪聲的壓制處理方法及裝置的制造方法
【專利摘要】本發(fā)明提供了一種諧波噪聲的壓制處理方法及裝置，其中，該方法包括：首先，利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對諧波樣點(diǎn)進(jìn)行標(biāo)記；再利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理；最后，對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)。在本發(fā)明提供的實施例中，采用兩種諧波噪聲的識別方式對地震高效采集帶來的諧波噪聲進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的信噪比和品質(zhì)。
【專利說明】
一種諧波噪聲的壓制處理方法及裝置
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及地震勘探技術(shù)領(lǐng)域，具體而言，涉及一種諧波噪聲的壓制處理方法及 裝置。
【背景技術(shù)】
[0002] 現(xiàn)階段的巖性地層和復(fù)雜構(gòu)造勘探已經(jīng)對地震勘探提出了薄儲層小斷層識別、精 細(xì)構(gòu)造圈閉目標(biāo)刻畫的高要求，因此，常規(guī)的地震勘探技術(shù)已經(jīng)難以滿足需求。"兩寬一高" 的地震勘探技術(shù)代表著地球物理技術(shù)的發(fā)展方向，它采用可控震源高效采集技術(shù)，具有較 好的經(jīng)濟(jì)可行性，并且，通過大帶道能力、功能強(qiáng)大的地震儀器及寬頻可控震源來實現(xiàn)。
[0003] 目前，高密度三維地震采集技術(shù)普遍采用可控震源滑動掃描方式，雖然大幅度提 高了采集效率，但是，在可控震源地震勘探的過程中，由于可控震源機(jī)械系統(tǒng)和液壓伺服系 統(tǒng)的非線性震動、以及震源與大地耦合等因素的影響，產(chǎn)生了諧波畸變。這種畸變以掃描信 號頻率范圍的倍數(shù)出現(xiàn)，形成了諧波干擾。對于相關(guān)后的地震記錄數(shù)據(jù)而言，諧波干擾將對 本炮其自身產(chǎn)生高次諧波干擾(在初至前)或低次諧波干擾(在初至后的數(shù)據(jù)中）。尤其是當(dāng) 可控震源地震勘探采用滑動掃描方式時，后一炮不等前一炮掃描結(jié)束就已經(jīng)開始激發(fā)，導(dǎo) 致前后兩炮甚至幾炮數(shù)據(jù)部分疊加在一起，后一炮的高次諧波干擾將會產(chǎn)生其自身的諧波 干擾(在初至前)和前一炮(甚至是前幾炮)數(shù)據(jù)的鄰炮干擾，前一炮的低次諧波干擾可能會 產(chǎn)生其自身的諧波干擾(在初至后的數(shù)據(jù)中）和后一炮(甚至是后幾炮)數(shù)據(jù)的鄰炮干擾，從 而形成了自身諧波干擾和鄰炮諧波干擾。以下將自身諧波干擾和鄰炮諧波干擾統(tǒng)稱為諧波 干擾，即諧波噪聲。
[0004] 在可控震源地震勘探的過程中，這種諧波噪聲嚴(yán)重降低了地震資料的信噪比，影 響了地震資料的品質(zhì)。

【發(fā)明內(nèi)容】

[0005] 有鑒于此，本發(fā)明實施例的目的在于提供一種諧波噪聲的壓制處理方法及裝置， 通過有效識別地震道數(shù)據(jù)中的諧波噪聲并對識別出的諧波噪聲進(jìn)行有效壓制，可以提高地 震資料的信噪比和品質(zhì)。
[0006] 第一方面，本發(fā)明實施例提供了一種諧波噪聲的壓制處理方法，該方法包括：
[0007] 利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的多個樣點(diǎn)中 含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對上述諧波樣點(diǎn) 進(jìn)行標(biāo)記；
[0008] 利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣點(diǎn)是否需要 壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理；
[0009] 對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)。
[0010] 結(jié)合第一方面，本發(fā)明實施例提供了第一方面的第一種可能的實施方式，其中，上 述利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的多個樣點(diǎn)中含諧波 噪聲的樣點(diǎn)之前，還包括：
[0011] 設(shè)置諧波噪聲的最小頻庫 并根據(jù)上述諧波噪聲的最 小頻率確定與上述最小頻率對應(yīng)的上述諧波噪聲的最大周期
其中，fmln表示 諧波噪聲的最小頻率，fbegln表示地震道正常子波振幅譜的開始頻率，fend表示地震道正常子 波振幅譜的結(jié)束頻率，fhg表示開始識別的諧波噪聲頻率，nz_表示諧波噪聲的最大周期， 單位為毫秒。
[0012] 結(jié)合第一方面的第一種可能的實施方式，本發(fā)明實施例提供了第一方面的第二種 可能的實施方式，其中，上述利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道 數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)， 并對上述諧波樣點(diǎn)進(jìn)行標(biāo)記，包括：
[0013] 步驟a:在待處理的地震道數(shù)據(jù)中，以切除點(diǎn)作為開始樣點(diǎn)，并在該切除點(diǎn)處創(chuàng)建 時窗長度Snwindows的時窗；
[0014] 步驟b:根據(jù)上述諧波噪聲的最大周期和上述時窗長度確定上述時窗內(nèi)諧波噪聲 的過零點(diǎn)個數(shù)nmin;
[0015] 步驟c:根據(jù)地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)對應(yīng)的振幅是否同號確定上述時窗內(nèi) 地震道數(shù)據(jù)的過零點(diǎn)個數(shù)ntr_;
[0016]步驟d:按照預(yù)設(shè)的時窗移動長度nMve將上述時窗向上述待處理的地震道數(shù)據(jù)的 道結(jié)束點(diǎn)處方向移動；
[0017] 步驟e:判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrace是否大于上述諧波噪聲的過零 點(diǎn)個數(shù)n min，若是，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對確定出的各個含有諧 波噪聲的樣點(diǎn)進(jìn)行標(biāo)識；
[0018] 依次循環(huán)步驟b至步驟e，直到上述時窗的下邊緣點(diǎn)移動至上述地震道數(shù)據(jù)的道結(jié) 束點(diǎn);判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)n tr_是否大于上述諧波噪聲的過零點(diǎn)個數(shù)^^ 若是，則將上述時窗內(nèi)的所有樣點(diǎn)均確定為含有諧波噪聲的樣點(diǎn)，并對確定出的各個含有 諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識。
[0019] 結(jié)合第一方面，本發(fā)明實施例提供了第一方面的第三種可能的實施方式，其中，上 述利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣點(diǎn)是否需要壓制，并 對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理，包括：
[0020] 步驟1:對上述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，以劃分得到多個不同頻率 段的頻帶；
[0021] 步驟2:從劃分得到的多個上述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻帶，并將上述標(biāo)準(zhǔn) 頻帶之后的多個頻帶作為待壓制頻帶；
[0022] 步驟3:從多個上述待壓制頻帶中，選取上述標(biāo)準(zhǔn)頻帶之后的兩個頻帶作為當(dāng)前待 壓制頻帶；
[0023]步驟4:在選取的兩個上述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣點(diǎn)中，分別逐一判斷 兩個上述當(dāng)前待壓制頻帶的各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅 是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若是，則對上述諧波樣點(diǎn)的 地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個上述當(dāng)前待壓制頻帶中與上述標(biāo)準(zhǔn) 頻帶相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶；
[0024] 依次循環(huán)步驟3至步驟4,直到當(dāng)前選取的下一個標(biāo)準(zhǔn)頻帶為最后一個頻帶。
[0025] 結(jié)合第一方面的第三種可能的實施方式，本發(fā)明實施例提供了第一方面的第四種 可能的實施方式，其中，上述在選取的上述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣點(diǎn)中，分別逐 一判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻 帶的半波的平均振幅與插值函數(shù)的乘積，包括：
[0026] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層(lmute<t<ldelaytm)，則判斷各個上述諧波 樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅 與淺層插值函數(shù) ?1(0的乘積；
[0027] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（ldelaytm彡t彡llength)，則判斷各個上述 諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均 振幅與中深層插值函數(shù)P2 (t)的乘積。
[0028] 結(jié)合第一方面的第四種可能的實施方式，本發(fā)明實施例提供了第一方面的第五種 可能的實施方式，其中，上述淺層插值函數(shù)關(guān)
[0029] 上述中深層插值函數(shù)為
[0030] 其中，dtime表示從切除時刻開始壓制的最小門檻值，etime表示距道結(jié)束時刻壓制的 最小門檻，lde3laytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，1^_表示地震道數(shù)據(jù)中的切 除點(diǎn)，1 i^th表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。
[0031] 結(jié)合第一方面的第三種可能的實施方式，本發(fā)明實施例提供了第一方面的第六種 可能的實施方式，其中，對上述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理，包括：
[0032] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t<ldelayt m)時，則根據(jù)公式
對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述當(dāng)前待壓制頻帶k中需要壓制的諧波 樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公式
U)對上述 標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道 數(shù)據(jù)進(jìn)行諧波壓制處理；
[0033] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（ldelaytm彡t彡llength)時，則根據(jù)公式
⑴對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述當(dāng)前待壓制頻帶k中需要壓制的諧波 樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公式
+1(0對上述 標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道 數(shù)據(jù)進(jìn)行諧波壓制處理；
[0034] 其中，t的取值為lmute至1 length，yk-l ( t )表示標(biāo)準(zhǔn)頻帶k_l中樣點(diǎn)t的半波的平均振 幅，y k (t)表示與標(biāo)準(zhǔn)頻帶k -1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，x k (t)表示與標(biāo)準(zhǔn)頻 帶k-1相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yk+1(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個頻帶 k +1中樣點(diǎn)t的半波的平均振幅，x k+1 (t)表示標(biāo)準(zhǔn)頻帶k -1之后的第二個頻帶k+1中樣點(diǎn)t的 原始地震道數(shù)據(jù)，Y表示振幅譜上振幅衰減系數(shù)。
[0035] 結(jié)合第一方面至第一方面的第六種可能的實施方式中的任一種可能的實施方式， 本發(fā)明實施例提供了第一方面的第七種可能的實施方式，其中，對各個頻帶對應(yīng)的壓制后 地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)，包括：
[0036] 根據(jù)公式訴)=￡毛⑴確定壓制后的地震道數(shù)據(jù)，其中，t表示樣點(diǎn)號變量，毛⑴ k 二i 表示頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfnum表示劃分得到的頻帶總數(shù)。
[0037]第二方面，本發(fā)明實施例還提供了一種諧波噪聲的壓制處理裝置，該裝置包括：
[0038] 第一識別模塊，用于利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道 數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)， 并對上述諧波樣點(diǎn)進(jìn)行標(biāo)記；
[0039] 第二識別模塊，用于利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的 諧波樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處 理；
[0040] 數(shù)據(jù)重構(gòu)模塊，用于對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后 的地震道數(shù)據(jù)。
[0041] 結(jié)合第二方面，本發(fā)明實施例提供了第二方面的第一種可能的實施方式，其中，上 述裝置還包括：
[0042] 參數(shù)設(shè)置模塊，用于設(shè)置諧波噪聲的最小頻率 并根 ，. 據(jù)上述諧波噪聲的最小頻率確定與上述最小頻率對應(yīng)的上述諧波噪聲的最大周期
其中，fmin表示諧波噪聲的最小頻率，f be3gin表示地震道正常子波振幅譜的開 始頻率，fend表示地震道正常子波振幅譜的結(jié)束頻率，fharm表示開始識別的諧波噪聲頻率， nz_表示諧波噪聲的最大周期，單位為毫秒。
[0043] 結(jié)合第二方面的第一種可能的實施方式，本發(fā)明實施例提供了第二方面的第二種 可能的實施方式，其中，上述第一識別模塊包括：
[0044] 時窗創(chuàng)建單元，用于在待處理的地震道數(shù)據(jù)中，以切除點(diǎn)作為開始樣點(diǎn)，并在該切 除點(diǎn)處創(chuàng)建時窗長度為nwind?d^時窗；
[0045] 第一確定單元，用于根據(jù)上述諧波噪聲的最大周期和上述時窗長度確定上述時窗 內(nèi)諧波噪聲的過零點(diǎn)個數(shù)nmin;
[0046] 第二確定單元，用于根據(jù)地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)對應(yīng)的振幅是否同號確定 上述時窗內(nèi)地震道數(shù)據(jù)的過零點(diǎn)個數(shù)ntr_;
[0047] 時窗移動單元，用于按照預(yù)設(shè)的時窗移動長度nMve將上述時窗向上述待處理的地 震道數(shù)據(jù)的道結(jié)束點(diǎn)處方向移動；
[0048] 諧波噪聲樣點(diǎn)標(biāo)記單元，用于判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrace是否大于 上述諧波噪聲的過零點(diǎn)個數(shù)n min，若是，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對 確定出的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識;直到上述時窗的下邊緣點(diǎn)移動至上述地震道 數(shù)據(jù)的道結(jié)束點(diǎn);判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)n tr_是否大于上述諧波噪聲的過零 點(diǎn)個數(shù)nmin，若是，則將上述時窗內(nèi)的所有樣點(diǎn)均確定為含有諧波噪聲的樣點(diǎn)，并對確定出 的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識。
[0049] 結(jié)合第二方面，本發(fā)明實施例提供了第二方面的第三種可能的實施方式，其中，上 述弟^?識別t旲塊包括：
[0050] 頻帶劃分單元，用于對上述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，以劃分得到 多個不同頻率段的頻帶；
[0051] 標(biāo)準(zhǔn)頻帶選取單元，用于從劃分得到的多個上述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻 帶，并將上述標(biāo)準(zhǔn)頻帶之后的多個頻帶作為待壓制頻帶；
[0052] 當(dāng)前待壓制頻帶確定單元，用于從多個上述待壓制頻帶中，選取上述標(biāo)準(zhǔn)頻帶之 后的兩個頻帶作為當(dāng)前待壓制頻帶；
[0053] 諧波樣點(diǎn)識別與壓制單元，用于在選取的兩個上述當(dāng)前待壓制頻帶的已標(biāo)記的諧 波樣點(diǎn)中，分別逐一判斷兩個上述當(dāng)前待壓制頻帶的各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制 頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若 是，則對上述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個上述當(dāng)前 待壓制頻帶中與上述標(biāo)準(zhǔn)頻帶相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶;直到當(dāng)前選取的下一個標(biāo) 準(zhǔn)頻帶為最后一個頻帶。
[0054]結(jié)合第二方面的第三種可能的實施方式，本發(fā)明實施例提供了第二方面的第四種 可能的實施方式，其中，上述諧波樣點(diǎn)識別與壓制單元包括：
[0055] 淺層諧波樣點(diǎn)識別子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t <ldelaytm)，則判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于 上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與淺層插值函數(shù) Pl(t)的乘積；
[0056] 中深層諧波樣點(diǎn)識別子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層 (ldelaytmStdlength)，則判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅 是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與中深層插值函數(shù)P2(t)的乘積。
[0057] 結(jié)合第二方面的第四種可能的實施方式，本發(fā)明實施例提供了第二方面的第五種 可能的實施方式，其中，上述淺層插值函數(shù)為
上述中深 層插值函數(shù)為；
[0058] 其中，dtime表示從切除時刻開始壓制的最小門檻值，etime表示距道結(jié)束時刻壓制的 最小門檻，lde3laytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，1^_表示地震道數(shù)據(jù)中的切 除點(diǎn)，1 i^th表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。
[0059]結(jié)合第二方面的第三種可能的實施方式，本發(fā)明實施例提供了第二方面的第六種 可能的實施方式，其中，上述諧波樣點(diǎn)識別與壓制單元包括：
[0060] 淺層諧波樣點(diǎn)壓制子單元，用干當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t Cldelaytm)時，則根據(jù)公式
i.(0對與上述標(biāo)準(zhǔn)頻帶k_ 1相鄰的上述當(dāng)前待壓制 頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；并根據(jù)公式
#)對上述標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的 諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；
[0061] 中深層諧波樣點(diǎn)壓制子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層 (ldelaytmStdlength)時，則根據(jù)公式
.(0.對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述 當(dāng)前待壓制頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根 據(jù)公式
⑴對上述標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要 壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；
[0062] 其中，t的取值為lmute至1 length，yk-l ( t )表示標(biāo)準(zhǔn)頻帶k_l中樣點(diǎn)t的半波的平均振 幅，y k (t)表示與標(biāo)準(zhǔn)頻帶k -1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，X k (t)表示與標(biāo)準(zhǔn)頻 帶k-1相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yk+1(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個頻帶 k +1中樣點(diǎn)t的半波的平均振幅，X k+1 (t)表示標(biāo)準(zhǔn)頻帶k -1之后的第二個頻帶k+1中樣點(diǎn)t的 原始地震道數(shù)據(jù)，Y表示振幅譜上振幅衰減系數(shù)。
[0063] 結(jié)合第二方面至第二方面的第六種可能的實施方式中的任一種可能的實施方式， 本發(fā)明實施例提供了第二方面的第七種可能的實施方式，其中，上述數(shù)據(jù)重構(gòu)模塊包括：
[0064] 數(shù)據(jù)重構(gòu)單元，用于根據(jù)公式所)=￡元(/)確定壓制后的地震道數(shù)據(jù)，其中，t表 示樣點(diǎn)號變量，元(0表示頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfmjm表示劃分得到的頻 帶總數(shù)。
[0065] 在本發(fā)明提供的實施例中，首先，利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待 處理的地震道數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的含諧波噪聲的樣點(diǎn)作為諧 波樣點(diǎn)，并對諧波樣點(diǎn)進(jìn)行標(biāo)記;再利用分頻的半波的平均振幅識別法確定各個頻帶中已 標(biāo)記的諧波樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓 制處理;最后，對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)。 在本發(fā)明提供的實施例中采用兩種諧波噪聲的識別方式對地震高效采集帶來的諧波噪聲 進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的信噪比和 品質(zhì)。
[0066] 為使本發(fā)明的上述目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉較佳實施例，并配合 所附附圖，作詳細(xì)說明如下。
【附圖說明】
[0067] 為了更清楚地說明本發(fā)明實施例的技術(shù)方案，下面將對實施例中所需要使用的附 圖作簡單地介紹，應(yīng)當(dāng)理解，以下附圖僅示出了本發(fā)明的某些實施例，因此不應(yīng)被看作是對 范圍的限定，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù)這 些附圖獲得其他相關(guān)的附圖。
[0068] 圖1示出了本發(fā)明實施例所提供的諧波噪聲的壓制處理方法的流程示意圖；
[0069] 圖2a示出了本發(fā)明實施例所提供的利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別 過程中時窗移動的示意圖；
[0070] 圖2b示出了本發(fā)明實施例所提供的利用分頻的半波的平均振幅識別法識別過程 中分頻后各個頻帶的地震道數(shù)據(jù)的示意圖；
[0071] 圖2c示出了本發(fā)明實施例所提供的利用分頻的半波的平均振幅識別法識別過程 中分頻后各個頻帶的平均振幅的示意圖；
[0072] 圖2d示出了本發(fā)明實施例所提供的利用分頻的半波的平均振幅識別法識別過程 中各個頻帶中相鄰頻帶重疊的示意圖；
[0073] 圖2e示出了本發(fā)明實施例所提供的利用分頻的半波的平均振幅識別法識別過程 中淺層與中深層的差值函數(shù)曲線的示意圖；
[0074] 圖3示出了本發(fā)明實施例所提供的通過實驗得到的原始單炮數(shù)據(jù)、按照相關(guān)技術(shù) 的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)、以及按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓制 后的單炮數(shù)據(jù)的對比示意圖；
[0075] 圖4中示出了本發(fā)明實施例所提供的與圖3中AC區(qū)域的三個單炮數(shù)據(jù)對應(yīng)的頻譜 分析圖；
[0076]圖5中示出了本發(fā)明實施例所提供的與圖3中BCD區(qū)域的三個單炮數(shù)據(jù)對應(yīng)的頻譜 分析圖；
[0077] 圖6中示出了本發(fā)明實施例所提供的諧波壓制前的單炮數(shù)據(jù)的疊加剖面、按照本 發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的疊加剖面、以及疊加后的諧波噪聲剖 面示意圖；
[0078] 圖7示出了本發(fā)明實施例所提供的諧波噪聲的壓制處理裝置的結(jié)構(gòu)示意圖。
【具體實施方式】
[0079] 下面將結(jié)合本發(fā)明實施例中附圖，對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完整 地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例。通常在 此處附圖中描述和示出的本發(fā)明實施例的組件可以以各種不同的配置來布置和設(shè)計。因 此，以下對在附圖中提供的本發(fā)明的實施例的詳細(xì)描述并非旨在限制要求保護(hù)的本發(fā)明的 范圍，而是僅僅表示本發(fā)明的選定實施例。基于本發(fā)明的實施例，本領(lǐng)域技術(shù)人員在沒有做 出創(chuàng)造性勞動的前提下所獲得的所有其他實施例，都屬于本發(fā)明保護(hù)的范圍。
[0080] 考慮到在可控震源地震勘探的過程中，由于可控震源機(jī)械系統(tǒng)和液壓伺服系統(tǒng)的 非線性震動、以及震源與大地耦合等因素的影響，會產(chǎn)生諧波畸變，從而形成諧波噪聲，這 種諧波噪聲會嚴(yán)重降低地震資料的信噪比，影響地震資料的品質(zhì)?；诖?，本發(fā)明實施例提 供了一種諧波噪聲的壓制處理方法，以及一種可以應(yīng)用該方法的諧波噪聲的壓制處理裝 置。
[0081] 諧波噪聲除了具有頻率范圍成倍增長的特性之外，還具有相關(guān)特性、時間特性及 振幅特性。需要對諧波噪聲的特性進(jìn)行分析，搞清諧波噪聲的分布特性、頻率特性及振幅特 性，才能找出識別諧波噪聲并進(jìn)行壓制的方法，從而提高地震資料信噪比。
[0082] 由于大地的吸收衰減作用，地震記錄的主要能量通常集中在12~28Hz范圍內(nèi)，而 采用滑動掃描方式得到的地震道數(shù)據(jù)的諧波噪聲的頻率通常大于30Hz;此外，滑動掃描每 組可控震源的震動要間隔一個滑動時間，而震動間隔時間愈長，諧波噪聲和有效信號的能 量差異就愈大，這就為利用諧波頻率和諧波的半波的平均振幅識別諧波噪聲并對諧波噪聲 進(jìn)行壓制提供了條件?；诖耍鐖D1所示，本發(fā)明實施例提供的一種諧波噪聲的壓制處理 方法，該方法包括步驟S102-S106，具體如下：
[0083] 首先，設(shè)置諧波噪聲的最小頻庫
并根據(jù)上述諧波噪 聲的最小頻率確定與上述最小頻率對應(yīng)的上述諧波噪聲的最大周期
其中， fmin表示諧波噪聲的最小頻率，fbegin表示地震道正常子波振幅譜的開始頻率，fend表示地震 道正常子波振幅譜的結(jié)束頻率，fharm表示開始識別的諧波噪聲頻率，nz_表示諧波噪聲的最 大周期，n z_的單位是毫秒；
[0084]然后，對待處理的地震道數(shù)據(jù)中的多個樣點(diǎn)進(jìn)行諧波噪聲的識別與壓制處理，具 體包括：
[0085]步驟S102:利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的多 個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對上述 諧波樣點(diǎn)進(jìn)行標(biāo)記，其中，待處理的地震道數(shù)據(jù)為WA(t)，樣點(diǎn)t對應(yīng)的諧波噪聲指示標(biāo)志 bf iag (t)，將各個樣點(diǎn)的初始bf iag (t)設(shè)置為0;
[0086] 步驟S104:利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣點(diǎn) 是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理；
[0087] 步驟S106:對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道 數(shù)據(jù)。
[0088] 在本發(fā)明提供的實施例中采用兩種諧波噪聲的識別方式對地震高效采集帶來的 諧波噪聲進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的 f目噪比和品質(zhì)。
[0089] 具體的，上述利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的 多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對上 述諧波樣點(diǎn)進(jìn)行標(biāo)記，包括：
[0090] 步驟a:在待處理的地震道數(shù)據(jù)中，以切除點(diǎn)作為開始樣點(diǎn)，并在該切除點(diǎn)處創(chuàng)建 時窗長度Snwindows的時窗；
[0091] 步驟b:根據(jù)上述諧波噪聲的最大周期和上述時窗長度確定上述時窗內(nèi)諧波噪聲 的過零點(diǎn)個數(shù)nmin，其中
nmin是一個固定值；
[0092] 步驟c:根據(jù)地震道數(shù)據(jù)中是否同號確定上述時窗內(nèi)地震道數(shù)據(jù)的過零點(diǎn)個數(shù) ntr_，具體的，將地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)所對應(yīng)的振幅依次兩兩相乘，即WA(t)*WA(t + 1)，若相乘結(jié)果小于零(相鄰的兩個樣點(diǎn)對應(yīng)的振幅異號，即一正值一負(fù)值），表不該相鄰 兩個樣點(diǎn)之間過零點(diǎn)，則ntr_計數(shù)一次(數(shù)值加一）；若相乘結(jié)果大于零(相鄰的兩個樣點(diǎn)對 應(yīng)的振幅同號，即同為正值或者同為負(fù)值），表示該相鄰兩個樣點(diǎn)之間未過零點(diǎn)，則 值不變;直到t =該時窗中的最后一個樣點(diǎn)，此時ntr_的數(shù)值標(biāo)識該時窗內(nèi)地震道數(shù)據(jù)的過 零點(diǎn)個數(shù)；
[0093]步驟d:按照預(yù)設(shè)的時窗移動長度nMve將上述時窗向上述待處理的地震道數(shù)據(jù)的 道結(jié)束點(diǎn)處方向移動；
[0094] 步驟e:判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrac；e是否大于上述諧波噪聲的過零 點(diǎn)個數(shù)nmin，若是(n tr_大于nmin)，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對確定出 的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識(將移除的個樣點(diǎn)t至t+n m_對應(yīng)的諧波噪聲指示 標(biāo)志bf lag ( t )至bf lag ( t+nmove )設(shè)置為99 );若ntrace小于nmin，則移除的nmove個樣點(diǎn)t至t+nmove對 應(yīng)的諧波噪聲指示標(biāo)志bf lag ( t )至bf lag ( t+nm_ )的數(shù)值不變，仍為0 ;
[0095] 依次循環(huán)步驟C至步驟e，直到上述時窗的下邊緣點(diǎn)移動至上述地震道數(shù)據(jù)的道結(jié) 束點(diǎn);判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntr_是否大于上述諧波噪聲的過零點(diǎn)個數(shù)^^ 若是，則將上述時窗內(nèi)的所有樣點(diǎn)（樣點(diǎn)1 length-nwindows至11 ength )均確定為含有諧波噪聲的 樣點(diǎn)，并對確定出的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識，其中，由于考慮到深層諧波干擾嚴(yán) 重，因而，對于最后一個時窗進(jìn)行單獨(dú)處理，如最后一個時窗內(nèi)ntrace大于nmin，則將該時窗內(nèi) 的所有樣點(diǎn)均標(biāo)記為含有諧波噪聲的樣點(diǎn)，將樣點(diǎn)11(^1^-11?:111〇1。>^至11(^1^對應(yīng)的&￡1 ￡18 (1 length-Ilwindows ) lag ( 1 length )-設(shè)置太?99 〇
[0096] 需要說明的是，考慮到對于每個時窗而言，諧波噪聲的最大周期nze5r。、時窗長度 nwind? s均相同，因而，每個時窗內(nèi)諧波噪聲的過零點(diǎn)個數(shù)nun是一個固定值，即每個時窗中 nmin大小是不變的，因而，依次循環(huán)步驟c至步驟e即可。
[0097]例如，如圖2a所示，待處理的地震道數(shù)據(jù)各個樣點(diǎn)t的取值為lmute至llength，每個時 窗長度為nwind_，時窗移動長度為nm_，上述切除點(diǎn)為lmute(假如未預(yù)先設(shè)置切除點(diǎn)時，則認(rèn) 為從t = 1的樣點(diǎn)作為開始識別的樣點(diǎn)，并從t = 1的樣點(diǎn)處開始創(chuàng)建時窗長度為nwindows的時 窗），舉例說明，若待處理的地震道數(shù)據(jù)的時間序列是7000毫秒，采樣間隔是2毫秒，因而地 震道數(shù)據(jù)中包含有7000/2 = 3500個樣點(diǎn)，即地震道數(shù)據(jù)中道結(jié)束點(diǎn)1 iength = 3500;將樣點(diǎn)1 至樣點(diǎn)lmute-l切除，其主要目的為切除初至波等干擾波，從而提高信噪比，可以將第500個 樣點(diǎn)作為切除點(diǎn)，即地震道數(shù)據(jù)中的切除點(diǎn)l?t e = 500，將作為lmute = 500創(chuàng)建時窗的起始 點(diǎn)，基于單位時間內(nèi)過零點(diǎn)次數(shù)的諧波頻率識別法從第500至3500的樣點(diǎn)中確定出哪些樣 點(diǎn)含有諧波噪聲，將含有諧波噪聲的樣點(diǎn)對應(yīng)的諧波噪聲指示標(biāo)志bn ag(t)的值設(shè)置為99。
[0098] 在本發(fā)明提供的實施例中，若確定出任一樣點(diǎn)t含有諧波噪聲，可以將該樣點(diǎn)t對 應(yīng)的諧波噪聲指示標(biāo)志bflag(t)的值設(shè)置為某一預(yù)設(shè)值，例如，設(shè)置為99,用以指示該樣點(diǎn)t 含有諧波噪聲，bflag(t)未設(shè)置預(yù)設(shè)值之前可以是0。通過不斷的移動時窗，比較時窗內(nèi)地震 道的過零點(diǎn)個數(shù)n trac；e5與諧波噪聲的過零點(diǎn)個數(shù)nmin的大小關(guān)系，進(jìn)而確定待處理的地震道 數(shù)據(jù)中各個樣點(diǎn)中是否含有諧波噪聲，完成基于波形過零點(diǎn)次數(shù)的諧波噪聲識別法的諧波 噪聲識別，即識別完成后，含由諧波噪聲的樣點(diǎn)t對應(yīng)的b flag(t) = 99。
[0099] 具體的，上述利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣 點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理，包括：
[0100] 步驟1:對上述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，以劃分得到多個不同頻率 段的頻帶；
[0101] 其中，根據(jù)公式
角定劃分得到的頻 帶總數(shù)，其中，nf_i表示劃分的頻帶個數(shù)，nre5S(5rve5表示保留頻帶個數(shù)(nre5Sf5rve5 = 2為常數(shù)參 量），f cutoff表示諧波噪聲的截至頻率，f harm表示開始識別的諧波噪聲的頻率，f begin表示地震 道數(shù)據(jù)正常子波振幅譜的開始頻率;計算得到nfnum之后，根據(jù)公式X? Xff)對原始 地震道數(shù)據(jù)進(jìn)行分頻處理，劃分得到nfnum個帶寬相等的頻帶，其中，hk(t)表示濾波因子，X (t)表示原始地震道數(shù)據(jù)，Xk(t)表示濾波后數(shù)據(jù)，即分頻后的數(shù)據(jù)，€)為卷積云算符，t表示 樣點(diǎn)號變量，k表示頻帶變量；
[0102] 其中，仍以若待處理的地震道數(shù)據(jù)的時間序列是7000毫秒，采樣間隔是2毫秒為 例，因而地震道數(shù)據(jù)中包含有7000/2 = 3500個樣點(diǎn)，則分成nfnum個頻帶之后，每一個頻帶也 包含有3500個樣點(diǎn)。
[0103]步驟2:從劃分得到的多個上述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻帶，并將上述標(biāo)準(zhǔn) 頻帶之后的多個頻帶作為待壓制頻帶；
[0104] 步驟3:從多個上述待壓制頻帶中，選取上述標(biāo)準(zhǔn)頻帶之后的兩個頻帶作為當(dāng)前待 壓制頻帶；
[0105] 步驟4:在選取的兩個上述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣點(diǎn)中，分別逐一判斷 兩個上述當(dāng)前待壓制頻帶的各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅 是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若是，則對上述諧波樣點(diǎn)的 地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個上述當(dāng)前待壓制頻帶中與上述標(biāo)準(zhǔn) 頻帶相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶;其中，正常子波振幅譜是將標(biāo)準(zhǔn)能量頻帶內(nèi)以半波 為窗口的平均絕對值能量作為該半波中每一點(diǎn)的標(biāo)準(zhǔn)能量值，根據(jù)公式
確定頻帶k中樣點(diǎn)t的半波的平均振幅，其中，mi和m2分別表示頻帶k的起始樣點(diǎn)號和結(jié)束樣 點(diǎn)號，Xk(j)表示頻帶k中的樣點(diǎn)j對應(yīng)的原始數(shù)據(jù)；
[0106] 依次循環(huán)步驟3至步驟4,直到當(dāng)前選取的下一個標(biāo)準(zhǔn)頻帶為最后一個頻帶。
[0107]例如，如圖2b所示，將待處理的地震道數(shù)據(jù)劃分為6個帶寬相等的頻帶，即nfnum = 6，k的取值為1至6，k=l的低頻頻帶不進(jìn)行基于分頻的諧波平均振幅識別法的諧波噪聲識 別與壓制，如圖2c所示，將k = 2的頻帶作為首次選取的標(biāo)準(zhǔn)頻帶(正常子波振幅譜），對k = 3 至k = 6的頻帶進(jìn)行基于分頻的諧波平均振幅識別法的諧波噪聲識別與壓制;正常子波振幅 譜是以半波為窗口的平均振幅作為該半波中每一點(diǎn)的標(biāo)準(zhǔn)值，為通過比較頻帶間的半波的 平均振幅來進(jìn)行諧波噪聲的識別與壓制提供了基礎(chǔ)，避免了相關(guān)技術(shù)中基于單點(diǎn)能量識別 諧波干擾的缺陷。
[0108]需要說明的是，待壓制頻帶可以是劃分得到的所有頻帶，也可以是劃分得到的頻 率較高的頻帶。優(yōu)選的，考慮到低頻不被破壞，因此，頻率最低的M-1個頻帶不作為待壓制的 頻帶，將k=M-l (如k = 2)的頻帶作為標(biāo)準(zhǔn)頻帶，從k=M(即k = 3)的頻帶開始進(jìn)行壓制，即將 k = M(k = 3)的頻帶至k = nf_(k = 6)的頻帶作為待壓制的頻帶，其中，M為預(yù)設(shè)正整數(shù)。例 如，當(dāng)M = 3時，第1個頻帶不進(jìn)行處理，第2個頻帶作為正常子波振幅譜標(biāo)準(zhǔn)頻帶，而第3個頻 帶作為開始壓制的諧波噪聲頻帶。其中，劃分得到的nfn?個頻帶，按照頻率由低到高的順序 進(jìn)行編號，編號越小的頻帶的頻率越低，即第1個頻帶的頻率最低，第nfnumf頻帶的頻率最 尚。
[0109]具體的，如圖2d所示仍以將待處理的地震道數(shù)據(jù)劃分為6個頻帶為例，將k = 2的頻 帶作為標(biāo)準(zhǔn)頻帶，將k = 3至k = 6的頻帶作為待壓制頻帶，即k = 2的頻帶的后1頻帶(k = 3)、 后2頻帶(k = 4)、后3頻帶(k = 5)、后4頻帶(k = 6);
[0110] 1)第一次識別與壓制:將k = 2的頻帶之后的兩個頻帶作為當(dāng)前待壓制頻帶即后1 頻帶(k = 3)、后2頻帶(k = 4)，在后1頻帶(k = 3)、后2頻帶(k = 4)的已標(biāo)記的諧波樣點(diǎn)中，分 別逐一對k = 3的頻帶和k = 4的頻帶中的已標(biāo)記的諧波樣點(diǎn)進(jìn)行基于分頻的諧波平均振幅 識別法的識別與壓制，即將各個諧波樣點(diǎn)的半波的平均振幅與標(biāo)準(zhǔn)頻帶的半波的平均振幅 與插值函數(shù)的乘積進(jìn)行比較，根據(jù)比較結(jié)果判斷諧波樣點(diǎn)知否需要壓制;然后，將壓制后的 k = 3的頻帶作為下一個標(biāo)準(zhǔn)頻帶；
[0111] 2)第二次識別與壓制:將k = 3的頻帶之后的兩個頻帶作為當(dāng)前待壓制頻帶即后2 頻帶(k = 4)、后3頻帶(k = 5)，在后2頻帶(k = 4)、后3頻帶(k = 5)的已標(biāo)記的諧波樣點(diǎn)中，分 別逐一對k = 4的頻帶和k = 5的頻帶中的已標(biāo)記的諧波樣點(diǎn)進(jìn)行基于分頻的諧波平均振幅 識別法的識別與壓制，即將各個諧波樣點(diǎn)的半波的平均振幅與標(biāo)準(zhǔn)頻帶的半波的平均振幅 與插值函數(shù)的乘積進(jìn)行比較，根據(jù)比較結(jié)果判斷諧波樣點(diǎn)知否需要壓制;然后，將壓制后的 k = 4的頻帶作為下一個標(biāo)準(zhǔn)頻帶(k = 4的頻帶經(jīng)過兩次識別與壓制）；
[0112] 3)第三次識別與壓制:將k = 4的頻帶之后的兩個頻帶作為當(dāng)前待壓制頻帶即后3 頻帶(k = 5)、后4頻帶(k = 6)，在后3頻帶(k = 5)、后4頻帶(k = 6)的已標(biāo)記的諧波樣點(diǎn)中，分 別逐一對k = 5的頻帶和k = 6的頻帶中的已標(biāo)記的諧波樣點(diǎn)進(jìn)行基于分頻的諧波平均振幅 識別法的識別與壓制，即將各個諧波樣點(diǎn)的半波的平均振幅與標(biāo)準(zhǔn)頻帶的半波的平均振幅 與插值函數(shù)的乘積進(jìn)行比較，根據(jù)比較結(jié)果判斷諧波樣點(diǎn)知否需要壓制;然后，將壓制后的 k = 5的頻帶作為下一個標(biāo)準(zhǔn)頻帶(k = 5的頻帶經(jīng)過兩次識別與壓制）；
[0113] 4)第四次識別與壓制:將k = 5的頻帶之后的兩個頻帶作為當(dāng)前待壓制頻帶即后4 頻帶(k = 6)，由于k = 5的頻帶之后只有一個待處理的頻帶，在后4頻帶(k = 6)的已標(biāo)記的諧 波樣點(diǎn)中，分別逐一對k = 6的頻帶中的已標(biāo)記的諧波樣點(diǎn)進(jìn)行基于分頻的諧波平均振幅識 別法的識別與壓制，即將各個諧波樣點(diǎn)的半波的平均振幅與標(biāo)準(zhǔn)頻帶的半波的平均振幅與 插值函數(shù)的乘積進(jìn)行比較，根據(jù)比較結(jié)果判斷諧波樣點(diǎn)知否需要壓制;然后，將壓制后的k =6的頻帶作為下一個標(biāo)準(zhǔn)頻帶(k = 6的頻帶經(jīng)過兩次識別與壓制），所有頻帶均識別與壓 制結(jié)束。
[0114] 綜上可知，當(dāng)將k-1頻帶作為標(biāo)準(zhǔn)頻帶，k頻帶(k-1頻帶的后1頻帶)和k+1頻帶(k-1 頻帶的后2頻帶)分別為當(dāng)前待壓制頻帶：
[0115] 1)第一次識別與壓制:k-1頻帶作為標(biāo)準(zhǔn)頻帶，對k頻帶和k+1頻帶進(jìn)行基于分頻的 諧波平均振幅識別法的識別與壓制；
[0116] 2)第二次識別與壓制:k頻帶作為標(biāo)準(zhǔn)頻帶，對k+1頻帶和k+2頻帶進(jìn)行基于分頻的 諧波平均振幅識別法的識別與壓制；
[0117] 3)第三次識別與壓制:k+1頻帶作為標(biāo)準(zhǔn)頻帶，對k+2頻帶和k+3頻帶進(jìn)行基于分頻 的諧波平均振幅識別法的識別與壓制；
[0118] 4)第四次識別與壓制:若k+3為最后一個頻帶，k+2頻帶作為標(biāo)準(zhǔn)頻帶，只對k+3頻 帶進(jìn)行基于分頻的諧波平均振幅識別法的識別與壓制;其中，k+1頻帶、k+2頻帶、k+3頻帶均 進(jìn)行了兩次進(jìn)行基于分頻的諧波平均振幅識別法的識別與壓制。
[0119] 在本發(fā)明實施例中，采用兩種諧波噪聲的識別方式對地震高效采集帶來的諧波噪 聲進(jìn)行有效識別，即先采用基于單位時間內(nèi)過零點(diǎn)次數(shù)的諧波頻率識別法，再采用基于分 頻的諧波平均振幅識別法，其中，在采用基于分頻的諧波平均振幅識別法進(jìn)行識別時，實際 上對于選取的標(biāo)準(zhǔn)頻帶之后的第二個頻帶起到最后一個頻帶（以劃分6個頻帶，且選第2頻 帶為標(biāo)準(zhǔn)頻帶為例，第4頻帶至第6頻帶)而言，基于分頻的諧波平均振幅識別法進(jìn)行了兩次 識別與壓制，即高頻頻帶需要進(jìn)行兩次識別和壓制，可以對兩個頻帶中重疊部分中的樣點(diǎn) 進(jìn)行有效的壓制，從而取得更好的諧波噪聲壓制效果。
[0120] 其中，考慮到諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層與中深層(包括中層和深層兩部分） 時的諧波噪聲干擾能力存在一定的差異，因而，如圖2e所示，在比較半波的平均振幅時諧波 樣點(diǎn)位于地震道數(shù)據(jù)的淺層或者中深層引入不同的差值函數(shù)，提高了對諧波噪聲進(jìn)行識別 與壓制的準(zhǔn)確度，基于此，上述在選取的上述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣點(diǎn)中，分別 逐一判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn) 頻帶的半波的平均振幅與插值函數(shù)的乘積，包括：
[0121] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層(lmute<t<ldelaytm)，則判斷各個上述諧波 樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅 與淺層插值函數(shù) ?1(0的乘積；
[0122] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（ldelaytm彡t彡llength)，則判斷各個上述 諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均 振幅與中深層插值函數(shù)P2 (t)的乘積。
[0123] 在本發(fā)明提供的實施例中，考慮到淺層和中深層諧波噪聲的能量有差異，淺層弱， 中深層強(qiáng)，因此，為了有效壓制諧波噪聲，而不傷及淺層和中深層的有效信號，本發(fā)明實施 例的方法中采用不同的差值函數(shù)來識別淺層、中深層，相較于相關(guān)技術(shù)的方法采用淺層和 中深層的識別參數(shù)不變的方法，本發(fā)明實施例的方法取得的效果更優(yōu)。
[0124] 具體的，上述淺層插值逑
[0125] 上述中深層插值函數(shù)為：
[0126] 其中，dtime5表示從切除時刻開始壓制的最小門檻值，表示距道結(jié)束時刻壓制的 最小門檻，l de3laytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，1^_表示地震道數(shù)據(jù)中的切 除點(diǎn)，1 i^th表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。
[0127] 其中，對上述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理，包括：
[0128] 當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t<ldelaytm)時，則根據(jù)公式
#)對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述當(dāng)前待壓制頻帶k中需要壓制的諧波樣 點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公式
對上述標(biāo) 準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù) 據(jù)進(jìn)行諧波壓制處理；
[0129]當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（lde5laytm<t彡llength)時，則根據(jù)公式
(0對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述當(dāng)前待壓制頻帶k中需要壓制的諧波 樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公式
k#)對上述 標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道 數(shù)據(jù)進(jìn)行諧波壓制處理；
[01 30]其中，t的取值為lmute至1 length，yk-l ( t )表示標(biāo)準(zhǔn)頻帶k-1中樣點(diǎn)t的半波的平均振 幅，y k (t)表示與標(biāo)準(zhǔn)頻帶k -1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，X k (t)表示與標(biāo)準(zhǔn)頻 帶k-1相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yk+1(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻 帶k+1中樣點(diǎn)t的半波的平均振幅，xk+i (t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻帶k+1中樣點(diǎn) t的原始地震道數(shù)據(jù)，Y表示振幅譜上振幅衰減系數(shù)。
[0131] 具體的，將k-1頻帶作為首次選取的標(biāo)準(zhǔn)頻帶，當(dāng)頻帶k中的諧波樣點(diǎn)位于地震道 數(shù)據(jù)的淺層（ImuteCtCldelaytm)時，判斷 yk(t)是否大于 yk-l(t) Xpi(t)，若 yk(t)>yk-l(t) X Pl(t)，則根據(jù)公式
對與標(biāo)準(zhǔn)頻帶k-1后1頻帶k中需要壓制的諧波樣點(diǎn)t對 應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；以及，當(dāng)頻帶k+1中的諧波樣點(diǎn)位于地震道數(shù)據(jù)的 淺層（lmute < t < ldelaytm)時，判斷yk+1 ( t )是否大于yk-1 ( t ) X pi ( t )，若yk+1 ( t ) > yk-1 ( t ) X pi (t)，則根據(jù)公式
+1(〇對與標(biāo)準(zhǔn)頻帶k-1后2頻帶k+1中需要壓制的諧波樣點(diǎn) t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；
[0132] 當(dāng)頻帶k中的諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（ldelaytm<t<llength)時，判斷y k (t)是否大于yk-i(t) Xp2(t)，若yk(t)>yk-i(t) Xp2(t)，則根據(jù)公式
對 與標(biāo)準(zhǔn)頻帶k-1后1頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處 理；以及，當(dāng)頻帶k+1中的諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（1 delaytm彡t彡1 length )時，判斷 yk + i(t)是否大于yk-i(t) Xp2(t)，若 yk + i(t)>yk -i(t) Xp2(t)，則根據(jù)公式
1(0對與標(biāo)準(zhǔn)頻帶k-1后2頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地 震道數(shù)據(jù)進(jìn)行諧波壓制處理。
[0133] 具體的，對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù) 據(jù)，包括：
[0134] 根據(jù)公式雄)=￡毛(0確定壓制后的地震道數(shù)據(jù)，其中，t表示樣點(diǎn)號變量乂 W k=i 表示頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfnum表示劃分得到的頻帶總數(shù)。
[0135] 在本發(fā)明提供的實施例中采用兩種諧波噪聲的識別方式對地震高效采集帶來的 諧波噪聲進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的 信噪比和品質(zhì);進(jìn)一步的，在采用基于分頻的諧波平均振幅識別法進(jìn)行諧波噪聲識別時，對 高頻頻帶中的諧波樣點(diǎn)進(jìn)行了兩次基于分頻的諧波平均振幅識別法的識別和壓制，因而可 以對兩個頻帶的重疊部分中的樣點(diǎn)進(jìn)行有效的壓制，從而取得更好的諧波噪聲壓制效果； 更進(jìn)一步的，考慮到諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層與中深層時的諧波噪聲干擾能力存在 一定的差異，因而，在比較半波的平均振幅時諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層或者中深層 引入不同的差值函數(shù)，提高了對諧波噪聲進(jìn)行識別與壓制的準(zhǔn)確度。
[0136] 通過實驗得到的實驗結(jié)果表明，相較于相關(guān)技術(shù)的方法，本發(fā)明實施例的方法能 夠降低地震資料的信噪比，取得更好的地震資料品質(zhì)。圖3中示出了原始單炮數(shù)據(jù)、按照相 關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)、以及按照本發(fā)明實施例的方法進(jìn)行諧波噪 聲壓制后的單炮數(shù)據(jù)的對比示意圖，圖4中示出了圖3中AC區(qū)域?qū)?yīng)的原始單炮數(shù)據(jù)的頻譜 分析、按照相關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜分析、以及按照本發(fā)明 實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜分析的對比示意圖。圖5中示出了圖3 中BCD區(qū)域?qū)?yīng)的原始單炮數(shù)據(jù)的頻譜分析、按照相關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的 單炮數(shù)據(jù)的頻譜分析、以及按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻 譜分析的對比示意圖。從圖3、圖4和圖5中可以明顯看出，與相關(guān)技術(shù)的方法相比，按照本發(fā) 明實施例的方法進(jìn)行諧波噪聲壓制后的數(shù)據(jù)更加干凈。圖4所示的頻譜分析圖表明，由于諧 波噪聲能量較強(qiáng)，30Hz-90Hz頻率之間的諧波高頻噪聲比較突出；在低頻lHz-5Hz中，按照本 發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜的有效成分要高于原始單炮數(shù) 據(jù)的頻譜以及按照相關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜；在5Hz-30Hz 中，按照相關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜與按照本發(fā)明實施例的方 法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜的有效成分較接近;在30Hz-90Hz中主要是諧波 噪聲，按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)中的諧波噪聲得到了有效 的壓制。圖5所示的頻譜分析圖表明，由于諧波能量較強(qiáng)，同樣30Hz-90Hz頻率之間諧波高頻 噪聲比較突出。按照相關(guān)技術(shù)的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜與按照本發(fā)明 實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的頻譜的低頻有效成分得到加強(qiáng)，同時 10Hz-30Hz有效成分也得到加強(qiáng)。但在lHz-30Hz內(nèi)，本發(fā)明實施例的方法的保幅程度高于相 關(guān)技術(shù)的方法;在30Hz-90Hz諧波噪聲按照本發(fā)明實施例的方法得到了有效壓制，從圖4和 圖5中可以看出，與相關(guān)技術(shù)的方法相比，按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的 單炮數(shù)據(jù)中的諧波噪聲得到了有效的壓制，并且按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓 制后的單炮數(shù)據(jù)的頻譜的光滑型更高，同時，單頻干擾壓制要更好，整個頻譜包絡(luò)線也更加 理想。圖6中示出了諧波壓制前的單炮數(shù)據(jù)的疊加剖面、按照本發(fā)明實施例的方法進(jìn)行諧波 噪聲壓制后的單炮數(shù)據(jù)的疊加剖面、以及疊加后的諧波噪聲剖面示意圖，從圖中可以看出， 壓制前與壓制后相比，按照本發(fā)明實施例的方法進(jìn)行諧波噪聲壓制后的單炮數(shù)據(jù)的疊加剖 面，其諧波噪聲得到了有效的壓制，信噪比提高。從主要勘探目的層成像效果改善明顯，取 得了令人滿意的效果。
[0137] 本發(fā)明實施例還提供了一種諧波噪聲的壓制處理裝置，如圖7所示，該裝置包括：
[0138] 第一識別模塊702,用于利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地 震道數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的上述含諧波噪聲的樣點(diǎn)作為諧波樣 點(diǎn)，并對上述諧波樣點(diǎn)進(jìn)行標(biāo)記；
[0139] 第二識別模塊704,用于利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo) 記的諧波樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制 處理；
[0140] 數(shù)據(jù)重構(gòu)模塊706,用于對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制 后的地震道數(shù)據(jù)。
[0141] 在本發(fā)明提供的實施例中采用兩種諧波噪聲的識別方式對地震高效采集帶來的 諧波噪聲進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的 f目噪比和品質(zhì)。
[0142] 進(jìn)一步的，上述裝置還包括：
[0143] 參數(shù)設(shè)置模塊，用于設(shè)置諧波噪聲的最小頻率 并根 據(jù)上述諧波噪聲的最小頻率確定與上述最小頻率對應(yīng)的上述諧波噪聲的最大周期
其中，fmin表示諧波噪聲的最小頻率，f be3gin表示地震道正常子波振幅譜的開 始頻率，fend表示地震道正常子波振幅譜的結(jié)束頻率，fharm表示開始識別的諧波噪聲頻率， nz_表示諧波噪聲的最大周期，單位為毫秒。
[0144] 進(jìn)一步的，上述第一識別模塊702包括：
[0145] 時窗創(chuàng)建單元，用于在待處理的地震道數(shù)據(jù)中，以切除點(diǎn)作為開始樣點(diǎn)，并在該切 除點(diǎn)處創(chuàng)建時窗長度為nwind?d^時窗；
[0146] 第一確定單元，用于根據(jù)上述諧波噪聲的最大周期和上述時窗長度確定上述時窗 內(nèi)諧波噪聲的過零點(diǎn)個數(shù)nmin;
[0147] 第二確定單元，用于根據(jù)地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)對應(yīng)的振幅是否同號確定 上述時窗內(nèi)地震道數(shù)據(jù)的過零點(diǎn)個數(shù)ntr_;
[0148] 時窗移動單元，用于按照預(yù)設(shè)的時窗移動長度nMve將上述時窗向上述待處理的地 震道數(shù)據(jù)的道結(jié)束點(diǎn)處方向移動；
[0149] 諧波噪聲樣點(diǎn)標(biāo)記單元，用于判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntr_是否大于 上述諧波噪聲的過零點(diǎn)個數(shù)n min，若是，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對 確定出的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識;直到上述時窗的下邊緣點(diǎn)移動至上述地震道 數(shù)據(jù)的道結(jié)束點(diǎn);判斷上述時窗內(nèi)地震道的過零點(diǎn)個數(shù)n tr_是否大于上述諧波噪聲的過零 點(diǎn)個數(shù)nmin，若是，則將上述時窗內(nèi)的所有樣點(diǎn)均確定為含有諧波噪聲的樣點(diǎn)，并對確定出 的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識。
[0150] 進(jìn)一步的，上述第二識別模塊704包括：
[0151]頻帶劃分單元，用于對上述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，以劃分得到 多個不同頻率段的頻帶；
[0152]標(biāo)準(zhǔn)頻帶選取單元，用于從劃分得到的多個上述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻 帶，并將上述標(biāo)準(zhǔn)頻帶之后的多個頻帶作為待壓制頻帶；
[0153] 當(dāng)前待壓制頻帶確定單元，用于從多個上述待壓制頻帶中，選取上述標(biāo)準(zhǔn)頻帶之 后的兩個頻帶作為當(dāng)前待壓制頻帶；
[0154] 諧波樣點(diǎn)識別與壓制單元，用于在選取的兩個上述當(dāng)前待壓制頻帶的已標(biāo)記的諧 波樣點(diǎn)中，分別逐一判斷兩個上述當(dāng)前待壓制頻帶的各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制 頻帶的半波的平均振幅是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若 是，則對上述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個上述當(dāng)前 待壓制頻帶中與上述標(biāo)準(zhǔn)頻帶相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶;直到當(dāng)前選取的下一個標(biāo) 準(zhǔn)頻帶為最后一個頻帶。
[0155] 進(jìn)一步的，上述諧波樣點(diǎn)識別與壓制單元包括：
[0156] 淺層諧波樣點(diǎn)識別子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層(lmute<t <ldelaytm)，則判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅是否大于 上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與淺層插值函數(shù) Pl(t)的乘積；
[0157] 中深層諧波樣點(diǎn)識別子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層 (ldelaytmStdlength)，則判斷各個上述諧波樣點(diǎn)的上述當(dāng)前待壓制頻帶的半波的平均振幅 是否大于上述標(biāo)準(zhǔn)頻帶的半波的平均振幅與中深層插值函數(shù)P2(t)的乘積。
[0158] 進(jìn)一步的，上述淺層插值函數(shù)為
上述中深層插值 函數(shù)為
[0159] 其中，dtime5表示從切除時刻開始壓制的最小門檻值，eti me5表示距道結(jié)束時刻壓制的 最小門檻，lde3laytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，1^_表示地震道數(shù)據(jù)中的切 除點(diǎn)，1 i^th表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。
[0160] 進(jìn)一步的，上述諧波樣點(diǎn)識別與壓制單元包括：
[0161]淺層諧波樣點(diǎn)壓制子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層(lmute5<t Cldelaytm)時，則根據(jù)公另
i.W對與上述標(biāo)準(zhǔn)頻帶k_ 1相鄰的上述當(dāng)前待壓制 頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；并根據(jù)公式
對上述標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要壓制的 諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；
[0162]中深層諧波樣點(diǎn)壓制子單元，用于當(dāng)上述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層 (ldelaytmStdlength)時，則根據(jù)公式
對與上述標(biāo)準(zhǔn)頻帶k_l相鄰的上述 當(dāng)前待壓制頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根 據(jù)公式
_對上述標(biāo)準(zhǔn)頻帶k_l之后的第二個上述待壓制頻帶k+1中需要 壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；
[01 63]其中，t的取值為lmute至1 length，yk-l ( t )表示標(biāo)準(zhǔn)頻帶k-1中樣點(diǎn)t的半波的平均振 幅，y k (t)表示與標(biāo)準(zhǔn)頻帶k -1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，x k (t)表示與標(biāo)準(zhǔn)頻 帶k-1相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yk+1(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻 帶k+1中樣點(diǎn)t的半波的平均振幅，xk+i (t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻帶k+1中樣點(diǎn) t的原始地震道數(shù)據(jù)，Y表示振幅譜上振幅衰減系數(shù)。
[0164] 進(jìn)一步的，上述數(shù)據(jù)重構(gòu)模塊706包括：
[0165] 數(shù)據(jù)重構(gòu)單元，用于根據(jù)公式⑴確定壓制后的地震道數(shù)據(jù)，其中，t表 k=l 示樣點(diǎn)號變量，.YA (0表示頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfmjm表示劃分得到的頻 帶總數(shù)。
[0166] 在本發(fā)明提供的實施例中采用兩種諧波噪聲的識別方式對地震高效采集帶來的 諧波噪聲進(jìn)行有效識別，再對識別出的諧波噪聲進(jìn)行了有效壓制，從而提高了地震資料的 信噪比和品質(zhì);進(jìn)一步的，在采用基于分頻的諧波平均振幅識別法進(jìn)行諧波噪聲識別時，對 高頻頻帶中的諧波樣點(diǎn)進(jìn)行了兩次基于分頻的諧波平均振幅識別法的識別和壓制，因而可 以對兩個頻帶的重疊部分中的樣點(diǎn)進(jìn)行有效的壓制，從而取得更好的諧波噪聲壓制效果； 更進(jìn)一步的，考慮到諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層與中深層時的諧波噪聲干擾能力存在 一定的差異，因而，在比較半波的平均振幅時諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層或者中深層 引入不同的差值函數(shù)，提高了對諧波噪聲進(jìn)行識別與壓制的準(zhǔn)確度。
[0167] 本發(fā)明的諧波噪聲的壓制處理裝置中包括的各個模塊可以為設(shè)備上的特定硬件 或者安裝于設(shè)備上的軟件或固件等。所屬領(lǐng)域的技術(shù)人員可以清楚地了解到，為描述的方 便和簡潔，前述描述的模塊和單元的具體工作過程，均可以參考上述方法實施例中的對應(yīng) 過程。在本申請所提供的幾個實施例中，應(yīng)該理解到，所揭露的諧波噪聲的壓制處理方法及 裝置，可以通過其它的方式實現(xiàn)。以上所描述的諧波噪聲的壓制處理裝置的實施例僅僅是 示意性的，例如，上述單元的劃分，僅僅為一種邏輯功能劃分，實際實現(xiàn)時可以有另外的劃 分方式，又例如，多個單元或組件可以結(jié)合或者可以集成到另一個系統(tǒng)，或一些特征可以忽 略，或不執(zhí)行。另一點(diǎn)，所顯示或討論的相互之間的耦合或直接耦合或通信連接可以是通過 一些通信接口，裝置或單元的間接耦合或通信連接，可以是電性，機(jī)械或其它的形式。
[0168] 所述作為分離部件說明的單元可以是或者也可以不是物理上分開的，作為單元顯 示的部件可以是或者也可以不是物理單元，即可以位于一個地方，或者也可以分布到多個 網(wǎng)絡(luò)單元上?？梢愿鶕?jù)實際的需要選擇其中的部分或者全部單元來實現(xiàn)本實施例方案的目 的。另外，在本發(fā)明各個實施例中的各功能單元可以集成在一個處理單元中，也可以是各個 單元單獨(dú)物理存在，也可以兩個或兩個以上單元集成在一個單元中。
[0169] 所述功能如果以軟件功能單元的形式實現(xiàn)并作為獨(dú)立的產(chǎn)品銷售或使用時，可以 存儲在一個計算機(jī)可讀取存儲介質(zhì)中?；谶@樣的理解，本發(fā)明的技術(shù)方案本質(zhì)上或者說 對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分或者該技術(shù)方案的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計 算機(jī)軟件產(chǎn)品存儲在一個存儲介質(zhì)中，包括若干指令用以使得一臺計算機(jī)設(shè)備(可以是個 人計算機(jī)，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本發(fā)明各個實施例所述方法的全部或部分步驟。 而前述的存儲介質(zhì)包括:U盤、移動硬盤、只讀存儲器(R0M，Read-0nly Memory)、隨機(jī)存取存 儲器(RAM，Random Access Memory)、磁碟或者光盤等各種可以存儲程序代碼的介質(zhì)。
[0170]以上所述，僅為本發(fā)明的【具體實施方式】，但本發(fā)明的保護(hù)范圍并不局限于此，任何 熟悉本技術(shù)領(lǐng)域的技術(shù)人員在本發(fā)明揭露的技術(shù)范圍內(nèi)，可輕易想到變化或替換，都應(yīng)涵 蓋在本發(fā)明的保護(hù)范圍之內(nèi)。因此，本發(fā)明的保護(hù)范圍應(yīng)所述以權(quán)利要求的保護(hù)范圍為準(zhǔn)。
【主權(quán)項】
1. 一種諧波噪聲的壓制處理方法，其特征在于，所述方法包括： 利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù)的多個樣點(diǎn)中含諧 波噪聲的樣點(diǎn)，將識別出的所述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對所述諧波樣點(diǎn)進(jìn)行 標(biāo)記； 利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波樣點(diǎn)是否需要壓制， 并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理； 對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述利用波形過零點(diǎn)次數(shù)的諧波噪聲識別 法識別待處理的地震道數(shù)據(jù)的《個樣點(diǎn)中含諧波噪聲的樣點(diǎn)之前，還包括： 設(shè)置諧波噪聲的最小頻率牛根據(jù)所述諧波噪聲的最小頻 率確定與所述最小頻率對應(yīng)的所述諧波噪聲的最大周巧庚中，fmin表示諧波, 噪聲的最小頻率，打egin表示地震道正常子波振幅譜的開始頻率，fend表示地震道正常子波振 幅譜的結(jié)束頻率，fW"表示開始識別的諧波噪聲頻率，rw。表示諧波噪聲的最大周期，單位 為毫秒。3. 根據(jù)權(quán)利要求2所述的方法，其特征在于，所述利用波形過零點(diǎn)次數(shù)的諧波噪聲識別 法識別待處理的地震道數(shù)據(jù)的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的所述含諧波噪聲 的樣點(diǎn)作為諧波樣點(diǎn)，并對所述諧波樣點(diǎn)進(jìn)行標(biāo)記，包括： 步驟a:在待處理的地震道數(shù)據(jù)中，W切除點(diǎn)作為開始樣點(diǎn)，并在所述切除點(diǎn)處創(chuàng)建時 窗長度為riwindDws的時窗； 步驟b:根據(jù)所述諧波噪聲的最大周期和所述時窗長度確定所述時窗內(nèi)諧波噪聲的過 零點(diǎn)個數(shù)Ilmin; 步驟C:根據(jù)地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)對應(yīng)的振幅是否同號確定所述時窗內(nèi)地震 道數(shù)據(jù)的過零點(diǎn)個數(shù)ntrace ; 步驟d:按照預(yù)設(shè)的時窗移動長度nmDve將所述時窗向所述待處理的地震道數(shù)據(jù)的道結(jié)束 點(diǎn)處方向移動； 步驟e:判斷所述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntraw是否大于所述諧波噪聲的過零點(diǎn)個 數(shù)nmin，若是，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對確定出的各個含有諧波噪 聲的樣點(diǎn)進(jìn)行標(biāo)識； 依次循環(huán)步驟C至步驟e，直到所述時窗的下邊緣點(diǎn)移動至所述地震道數(shù)據(jù)的道結(jié)束 點(diǎn)；判斷所述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrace是否大于所述諧波噪聲的過零點(diǎn)個數(shù)nmin，若 是，則將所述時窗內(nèi)的所有樣點(diǎn)均確定為含有諧波噪聲的樣點(diǎn)，并對確定出的各個含有諧 波噪聲的樣點(diǎn)進(jìn)行標(biāo)識。4. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述利用分頻的半波的平均振幅識別法確 定各個頻帶中已標(biāo)記的諧波樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波 樣點(diǎn)進(jìn)行諧波壓制處理，包括： 步驟1 :對所述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，W劃分得到多個不同頻率段的 頻帶； 步驟2:從劃分得到的多個所述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻帶，并將所述標(biāo)準(zhǔn)頻帶 之后的多個頻帶作為待壓制頻帶； 步驟3:從多個所述待壓制頻帶中，選取所述標(biāo)準(zhǔn)頻帶之后的兩個頻帶作為當(dāng)前待壓制 頻帶； 步驟4:在選取的兩個所述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣點(diǎn)中，分別逐一判斷兩個 所述當(dāng)前待壓制頻帶的各個所述諧波樣點(diǎn)的所述當(dāng)前待壓制頻帶的半波的平均振幅是否 大于所述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若是，則對所述諧波樣點(diǎn)的地震 道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個所述當(dāng)前待壓制頻帶中與所述標(biāo)準(zhǔn)頻帶 相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶； 依次循環(huán)步驟3至步驟4,直到當(dāng)前選取的下一個標(biāo)準(zhǔn)頻帶為最后一個頻帶。5. 根據(jù)權(quán)利要求4所述的方法，其特征在于，所述在選取的所述當(dāng)前待壓制頻帶的已標(biāo) 記的諧波樣點(diǎn)中，分別逐一判斷各個所述諧波樣點(diǎn)的所述當(dāng)前待壓制頻帶的半波的平均振 幅是否大于所述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積，包括： 當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層(lmute<t<ldelaytm)，則判斷各個所述諧波樣點(diǎn) 的所述當(dāng)前待壓制頻帶的半波的平均振幅是否大于所述標(biāo)準(zhǔn)頻帶的半波的平均振幅與淺 層插值函數(shù)Pl ( t )的乘積； 當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層(ldelayt"《t《llength)，則判斷各個所述諧波 樣點(diǎn)的所述當(dāng)前待壓制頻帶的半波的平均振幅是否大于所述標(biāo)準(zhǔn)頻帶的半波的平均振幅 與中深層插值函數(shù)P2 ( t )的乘積。6. 根據(jù)權(quán)利要求5所述的方法，其特征在于，所述淺層插值函數(shù)為 所述中深層插值函數(shù)為其中，dtime表示從切除時刻開始壓制的最小口檻值，etime表示距道結(jié)束時刻壓制的最小 口檻，Idelaytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，Uute表示地震道數(shù)據(jù)中的切除 點(diǎn)，Ilength表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。 7 .根據(jù)權(quán)利要求4所述的方法，其特征在于，對所述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪 聲的壓制處理，包括： 當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t< Idelaytm)時，則根據(jù)公式f與所述標(biāo)準(zhǔn)頻帶k-1相鄰的所述當(dāng)前待壓制頻帶k中需要壓制的諧波樣 點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公支4片)對所述標(biāo) 準(zhǔn)頻帶k-1之后的第二個所述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù) 據(jù)進(jìn)行諧波壓制處理； 當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層Ilength)時，則根據(jù)公式j(luò)與所述標(biāo)準(zhǔn)頻帶k-1相鄰的所述當(dāng)前待壓制頻帶k中需要壓制的諧波 樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理;并根據(jù)公支>]所述 標(biāo)準(zhǔn)頻帶k-1之后的第二個所述待壓制頻帶k+1中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道 數(shù)據(jù)進(jìn)行諧波壓制處理； 其中，t的取值為Imute至1 length，yk-1 ( t )表示標(biāo)準(zhǔn)頻帶k-1中樣點(diǎn)t的半波的平均振幅，yk (t)表示與標(biāo)準(zhǔn)頻帶k-1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，Xk(t)表示與標(biāo)準(zhǔn)頻帶k-1 相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yw(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個頻帶k+1中 樣點(diǎn)t的半波的平均振幅，Xk+l(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個頻帶k+1中樣點(diǎn)t的原始地 震道數(shù)據(jù)，丫表示振幅譜上振幅衰減系數(shù)。8. 根據(jù)權(quán)利要求1-7任一項所述的方法，其特征在于，對各個頻帶對應(yīng)的壓制后地震道 數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地震道數(shù)據(jù)，包括： 根據(jù)公JI定壓制后的地震道數(shù)據(jù)，其中，t表示樣點(diǎn)號變量，％似表示 頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfn?表示劃分得到的頻帶總數(shù)。9. 一種諧波噪聲的壓制處理裝置，其特征在于，所述裝置包括： 第一識別模塊，用于利用波形過零點(diǎn)次數(shù)的諧波噪聲識別法識別待處理的地震道數(shù)據(jù) 的多個樣點(diǎn)中含諧波噪聲的樣點(diǎn)，將識別出的所述含諧波噪聲的樣點(diǎn)作為諧波樣點(diǎn)，并對 所述諧波樣點(diǎn)進(jìn)行標(biāo)記； 第二識別模塊，用于利用分頻的半波的平均振幅識別法確定各個頻帶中已標(biāo)記的諧波 樣點(diǎn)是否需要壓制，并對確定出的各個頻帶中需要壓制的諧波樣點(diǎn)進(jìn)行諧波壓制處理； 數(shù)據(jù)重構(gòu)模塊，用于對各個頻帶對應(yīng)的壓制后地震道數(shù)據(jù)進(jìn)行重構(gòu)，得到壓制后的地 震道數(shù)據(jù)。10. 根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述裝晉還包巧： 參數(shù)設(shè)置模塊，用于設(shè)置諧波噪聲的最小頻聾 并根據(jù)所述諧 ? 波噪聲的最小頻率確定與所述最小頻率對應(yīng)的所述諧波噪聲的最大周期岸 中，fmin表示諧波噪聲的最小頻率，fbegin表示地震道正常子波振幅譜的開始頻率，fend表示地 震道正常子波振幅譜的結(jié)束頻率，fharm表示開始識別的諧波噪聲頻率，nzer。表示諧波噪聲的 最大周期，單位為毫秒。 11 .根據(jù)權(quán)利要求10所述的裝置，其特征在于，所述第一識別模塊包括： 時窗創(chuàng)建單元，用于在待處理的地震道數(shù)據(jù)中，W切除點(diǎn)作為開始樣點(diǎn)，并在所述切除 點(diǎn)處創(chuàng)建時窗長度為nwindDws的時窗； 第一確定單元，用于根據(jù)所述諧波噪聲的最大周期和所述時窗長度確定所述時窗內(nèi)諧 波噪聲的過零點(diǎn)個數(shù)nmin; 第二確定單元，用于根據(jù)地震道數(shù)據(jù)中相鄰的兩個樣點(diǎn)對應(yīng)的振幅是否同號確定所述 時窗內(nèi)地震道數(shù)據(jù)的過零點(diǎn)個數(shù)ntrace ; 時窗移動單元，用于按照預(yù)設(shè)的時窗移動長度nmDve將所述時窗向所述待處理的地震道 數(shù)據(jù)的道結(jié)束點(diǎn)處方向移動； 諧波噪聲樣點(diǎn)標(biāo)記單元，用于判斷所述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrace是否大于所述 諧波噪聲的過零點(diǎn)個數(shù)nmin，若是，則將移出的樣點(diǎn)確定為含有諧波噪聲的樣點(diǎn)，并對確定 出的各個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識;直到所述時窗的下邊緣點(diǎn)移動至所述地震道數(shù)據(jù) 的道結(jié)束點(diǎn);判斷所述時窗內(nèi)地震道的過零點(diǎn)個數(shù)ntrase是否大于所述諧波噪聲的過零點(diǎn)個 數(shù)nmin，若是，則將所述時窗內(nèi)的所有樣點(diǎn)均確定為含有諧波噪聲的樣點(diǎn)，并對確定出的各 個含有諧波噪聲的樣點(diǎn)進(jìn)行標(biāo)識。 12 .根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述第二識別模塊包括： 頻帶劃分單元，用于對所述待處理的地震道數(shù)據(jù)進(jìn)行頻帶劃分處理，W劃分得到多個 不同頻率段的頻帶； 標(biāo)準(zhǔn)頻帶選取單元，用于從劃分得到的多個所述頻帶中選取一個頻帶作為標(biāo)準(zhǔn)頻帶， 并將所述標(biāo)準(zhǔn)頻帶之后的多個頻帶作為待壓制頻帶； 當(dāng)前待壓制頻帶確定單元，用于從多個所述待壓制頻帶中，選取所述標(biāo)準(zhǔn)頻帶之后的 兩個頻帶作為當(dāng)前待壓制頻帶； 諧波樣點(diǎn)識別與壓制單元，用于在選取的兩個所述當(dāng)前待壓制頻帶的已標(biāo)記的諧波樣 點(diǎn)中，分別逐一判斷兩個所述當(dāng)前待壓制頻帶的各個所述諧波樣點(diǎn)的所述當(dāng)前待壓制頻帶 的半波的平均振幅是否大于所述標(biāo)準(zhǔn)頻帶的半波的平均振幅與插值函數(shù)的乘積;若是，貝U 對所述諧波樣點(diǎn)的地震道數(shù)據(jù)進(jìn)行諧波噪聲的壓制處理;將壓制后的兩個所述當(dāng)前待壓制 頻帶中與所述標(biāo)準(zhǔn)頻帶相鄰的頻帶作為下一個標(biāo)準(zhǔn)頻帶;直到當(dāng)前選取的下一個標(biāo)準(zhǔn)頻帶 為最后一個頻帶。13.根據(jù)權(quán)利要求12所述的裝置，其特征在于，所述諧波樣點(diǎn)識別與壓制單元包括： 淺層諧波樣點(diǎn)識別子單元，用于當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t< IdelaytJ，則判斷各個所述諧波樣點(diǎn)的所述當(dāng)前待壓制頻帶的半波的平均振幅是否大于所 述標(biāo)準(zhǔn)頻帶的半波的平均振幅與淺層插值函數(shù)Pl(t)的乘積；中深層諧波樣點(diǎn)識別子單元，用于當(dāng)所述諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層(ldelayt"《 t《llength)，則判斷各個所述諧波樣點(diǎn)的所述當(dāng)前待壓制頻帶的半波的平均振幅是否大于 所述標(biāo)準(zhǔn)頻帶的半波的平均振幅與中深層插值函數(shù)P2(t)的乘積。 14 .根據(jù)權(quán)利要求1 3所述的裝置，其特征在于，所述淺層插值函數(shù)文 所述中深層插值函數(shù)為其中，dtime表示從切除時刻開始壓制的最小口檻值，etime表示距道結(jié)束時刻壓制的最小 口檻，Idelaytm表示地震道數(shù)據(jù)的淺層與中深層間的臨界點(diǎn)，Imute表示地震道數(shù)據(jù)中的切除 點(diǎn)，Ilength表示地震道數(shù)據(jù)中道結(jié)束點(diǎn)。15.根據(jù)權(quán)利要求12所述的裝置，其特征在于，所述諧波樣點(diǎn)識別與壓制單元包括：淺層諧波樣點(diǎn)壓帶ISK ^ 述諧波樣點(diǎn)位于地震道數(shù)據(jù)的淺層（lmute<t< ldelaytm)時，則根據(jù)公^^ 巧與所述標(biāo)準(zhǔn)頻帶k-1相鄰的所述當(dāng)前待壓制頻 帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；并根據(jù)公式>]所述標(biāo)準(zhǔn)頻帶k-1之后的第二個所述待壓制頻帶k+1中需要壓制的 諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理； 中深層諧波樣點(diǎn)壓制子單元，巧干巧所沐諧波樣點(diǎn)位于地震道數(shù)據(jù)的中深層（ldelaytm《 *《116。加)時，則根據(jù)公^5與所述標(biāo)準(zhǔn)頻帶k-1相鄰的所述當(dāng)前待壓 制頻帶k中需要壓制的諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理；并根據(jù)公式 如計所述標(biāo)準(zhǔn)頻帶k-1之后的第二個所述待壓制頻帶k+1中需要壓制的 諧波樣點(diǎn)t對應(yīng)的原始地震道數(shù)據(jù)進(jìn)行諧波壓制處理； 其中，t的取值為Imute至1 length，yk-1 ( t )表示標(biāo)準(zhǔn)頻帶k-1中樣點(diǎn)t的半波的平均振幅，yk (t)表示與標(biāo)準(zhǔn)頻帶k-1相鄰的頻帶k中樣點(diǎn)t的半波的平均振幅，Xk(t)表示與標(biāo)準(zhǔn)頻帶k-1 相鄰的頻帶k中樣點(diǎn)t的原始地震道數(shù)據(jù)，yw(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻帶k+1 中樣點(diǎn)t的半波的平均振幅，Xk+l(t)表示標(biāo)準(zhǔn)頻帶k-1之后的第二個的頻帶k+1中樣點(diǎn)t的原 始地震道數(shù)據(jù)，丫表示振幅譜上振幅衰減系數(shù)。16.根據(jù)權(quán)利要求9-15任一項所述的裝置，其特征在于，所述數(shù)據(jù)重構(gòu)模塊包括： 數(shù)據(jù)重構(gòu)單元，用于根據(jù)公^!I定壓制后的地震道數(shù)據(jù)，其中，t表示樣 點(diǎn)號變量，表示頻帶k中樣點(diǎn)t對應(yīng)的壓制后地震道數(shù)據(jù)，nfn?表示劃分得到的頻帶總 數(shù)。
【文檔編號】G01V1/28GK105911590SQ201610299277
【公開日】2016年8月31日
【申請日】2016年5月9日
【發(fā)明人】劉宏杰
【申請人】劉宏杰