一種地震記錄規(guī)則干擾波的提取方法及裝置的制造方法
【專利摘要】本發(fā)明涉及地震資料處理領(lǐng)域，尤其涉及一種地震記錄規(guī)則干擾波的提取方法及裝置。該方法包括：獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄；對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排齊處理后的地震記錄；計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)；根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，對(duì)所述時(shí)移排齊處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排齊處理后的地震記錄中心道對(duì)應(yīng)的疊加地震道；根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄。本申請(qǐng)實(shí)施例的方法能有效的實(shí)現(xiàn)對(duì)干擾波的提取，為地震資料的解釋奠定了基礎(chǔ)。
【專利說明】
一種地震記錄規(guī)則干擾波的提取方法及裝置
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及地震資料處理領(lǐng)域，尤其涉及一種地震記錄規(guī)則干擾波的提取方法及 裝置。
【背景技術(shù)】
[0002] 在采集野外實(shí)際地震數(shù)據(jù)包含著豐富的信息，除了能夠反映地下構(gòu)造的有效反射 波，通常地震記錄上還包含著各種各樣的干擾波。干擾波可分為兩種，具有明顯傳播規(guī)律的 規(guī)則干擾波，如聲波、面波，多次波等，和沒有明顯傳播規(guī)律性的隨機(jī)干擾波。干擾波的存在 會(huì)降低有效波的質(zhì)量，影響我們對(duì)地下地質(zhì)目標(biāo)的認(rèn)識(shí)，所以，如何在對(duì)地震記錄上的有效 反射波沒有明顯傷害的同時(shí)，將地震記錄中的干擾波去除，從而提高有效信號(hào)的分辨率，始 終是地震資料處理中的一個(gè)關(guān)鍵問題。
[0003] 目前，根據(jù)不同的地震波屬性，有很多種規(guī)則干擾波提取或壓制方法，如中值濾波 法、F-K濾波法、ip變換法、矢量合成法、極化分析法、帶通濾波法等。其中，中值濾波法、F-K 濾波法、T-p變換法主要是利用了規(guī)則干擾波和有效波的視速度差異;矢量合成法和極化分 析法利用了規(guī)則干擾波和有效波的偏振特性的不同；帶通濾波法則是根據(jù)規(guī)則干擾波和有 效波分布的頻率不同進(jìn)行分離。
[0004] 上述規(guī)則干擾波提取方法都是基于規(guī)則干擾波和有效波在視速度、偏振特性和頻 率等屬性方面的明顯差異，當(dāng)規(guī)則干擾波和有效波的上述屬性差異很小或者沒有差異，上 述方法就不能對(duì)干擾波進(jìn)行有效提取或壓制。因此，如何有效的提取地震記錄中的規(guī)則干 擾波仍是目前地震資料處理中亟需解決的問題。

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

[0005] 本申請(qǐng)實(shí)施例提供了一種地震記錄規(guī)則干擾波的提取方法及裝置，以有效的提取 地震記錄中的規(guī)則干擾波。
[0006] 為達(dá)到上述目的，一方面，本申請(qǐng)實(shí)施例提供了一種地震記錄規(guī)則干擾波的提取 方法，所述方法包括：
[0007] 獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄；
[0008] 對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排齊處理后的地震記 錄；
[0009] 計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和 能量權(quán)系數(shù)；
[0010] 根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，對(duì)所述時(shí)移排齊 處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排齊處理后的地震記錄中心 道對(duì)應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干擾波；
[0011] 根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄;其中，所述剩余地 震記錄為去除規(guī)則干擾波的地震記錄。
[0012] 進(jìn)一步地，在根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄之后， 所述方法還包括：
[0013] 逐步減小所述滑動(dòng)窗的大小，對(duì)所述剩余地震記錄進(jìn)行規(guī)則干擾波去除處理。
[0014] 進(jìn)一步地，采用以下計(jì)算公式獲得所述時(shí)移排齊處理后地震記錄中心道對(duì)應(yīng)的疊 加地震道：
[0015] ^(j) = ^ii,(t)W(t)PV(t) k=l.
[0016] 其中，s(t)表示疊加地震道；表示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)； 表示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù);uk(t)表示時(shí)移排齊處理后的地震記錄上 第k道地震道在時(shí)間樣點(diǎn)序號(hào)為t時(shí)的振幅值;k表示地震道的序號(hào);Mi表示滑動(dòng)窗內(nèi)地震記 錄的道數(shù);t表示時(shí)間樣點(diǎn)序號(hào)，N為每道地震記錄的時(shí)間采樣點(diǎn)數(shù)。
[0017] 進(jìn)一步地，根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄，包括：
[0018] 對(duì)所述疊加地震道進(jìn)行反時(shí)移校正處理，獲得反時(shí)移校正處理后的疊加地震道；
[0019] 根據(jù)所述疊前地震記錄和所述反時(shí)移校正處理后的疊加地震道，獲取振幅匹配系 數(shù)；
[0020] 根據(jù)所述振幅匹配系數(shù)對(duì)所示疊加地震道進(jìn)行振幅匹配處理，利用所述疊前地震 記錄減去振幅匹配處理后的疊加地震道，獲得剩余地震記錄。
[0021] 進(jìn)一步地，采用以下公式計(jì)算振幅匹配系數(shù)：
[0022] 乞 V-、. I i=k-
[0023] 其中，a」為第j道的振幅匹配系數(shù)，Ujl為第j道疊前地震記錄第i個(gè)采樣點(diǎn)的振幅 值，Sl為第j道地震記錄對(duì)應(yīng)的疊加地震道第i個(gè)采樣點(diǎn)的振幅值，i b和h為振幅計(jì)算時(shí)窗的 起止樣點(diǎn)序號(hào)。
[0024]進(jìn)一步地，計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似 權(quán)系數(shù)，包括：
[0025]計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的相關(guān)系數(shù)；
[0026]根據(jù)所述地震道的相關(guān)系數(shù)計(jì)算所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)；
[0027]對(duì)所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理 后的相似權(quán)系數(shù)計(jì)算所述地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)。
[0028] 進(jìn)一步地，采用以下公式計(jì)算所述地震道的相關(guān)系數(shù)：
[0029] , = /^,^!< (0' Uj(/), k, j = 1,? ? ?,
[0030] 其中，k和j是滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是滑動(dòng)窗內(nèi)地震記錄道數(shù)，i是滑動(dòng)窗內(nèi) 的時(shí)間樣點(diǎn)序號(hào)，ib和分別是當(dāng)前滑動(dòng)窗第一個(gè)和最后一個(gè)時(shí)間樣點(diǎn)序號(hào)，u k(i)為第k道 第i個(gè)時(shí)間樣點(diǎn)的振幅值，^⑴為第j道第i個(gè)時(shí)間樣點(diǎn)的振幅值，Cu為第k道和第j道地震 記錄的相關(guān)系數(shù)。
[0031] 進(jìn)一步地，采用以下公式計(jì)算所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)：
[0033]其中，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)，Ckj為第k道和 第j道的互相關(guān)系數(shù)，Cw和Ckk分別為第j道和第k道的自相關(guān)系數(shù)，i。為中心采樣點(diǎn)的時(shí)間序 號(hào)，wib為第k道在i。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)。
[0034]進(jìn)一步地，采用以下公式對(duì)所述地震道中心樣點(diǎn)的相似權(quán)系數(shù)進(jìn)行歸一化處理：
[0036]其中，為第k道在i。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)，為第j道在i。時(shí)間樣點(diǎn) 處的相似權(quán)系數(shù)，為第k道在i。時(shí)間樣點(diǎn)處的歸一化后的相似權(quán)系數(shù)，i。為當(dāng)前滑動(dòng) 窗中心采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。 [0037]進(jìn)一步地，計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量 權(quán)系數(shù)，包括：
[0038] 根據(jù)所述時(shí)移排齊處理后的地震記錄獲取能量濾波參數(shù)；
[0039] 根據(jù)所述能量濾波參數(shù)獲得能量濾波算子；
[0040] 根據(jù)所述能量濾波算子計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道中心采樣 點(diǎn)的能量權(quán)系數(shù)；
[0041] 對(duì)所述地震道中心采樣點(diǎn)的能量權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理 后的能量權(quán)系數(shù)計(jì)算所地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)。
[0042] 進(jìn)一步地，根據(jù)所述時(shí)移排齊處理后的地震記錄獲取能量濾波參數(shù)，包括：
[0043]計(jì)算所述時(shí)移排齊處理后的地震記錄中每一地震道的能量；
[0044] 根據(jù)預(yù)設(shè)的能量比例閾值，將所述時(shí)移排齊處理后的地震記錄中能量小于閾值的 地震道去除；
[0045] 計(jì)算所述時(shí)移排齊處理后的地震記錄中剩余地震道的平均能量值，并根據(jù)所述每 一地震道的能量與所述平均能量值計(jì)算所述時(shí)移排齊處理后的地震記錄的能量標(biāo)準(zhǔn)偏差；
[0046] 根據(jù)所述能量標(biāo)準(zhǔn)偏差獲取能量濾波參數(shù)。
[0047] 進(jìn)一步地，采用以下公式計(jì)算所述能量濾波參數(shù)：
[0048] \ " ^ e e$^ef"+l.5rae e4^em + 3-<tc
[0049] 其中，ei、e2、e3和e4表示能量濾波參數(shù)，em表示平均能量值，〇e表示能量標(biāo)準(zhǔn)偏差。
[0050] 進(jìn)一步地，采用以下公式計(jì)算所述能量濾波算子： 1 < ec < e3 、 0,5-().5*cos(/T(e. -e,)/(e. -e,)) e： < e(. < e.
[0051 ] f{e,-cn..a；)^< ^ " ! ' 1 ~ 0.5-0,^-cos(/T(ec ~e4)/(e. -e_J) e. < ec < e4 , 〇 ec <ei,ec >e4
[0052] 其中，ei、e2、e3和e4表示能量濾波參數(shù)，em表示平均能量值，e c表示地震道能量值。
[0053]進(jìn)一步地，采用以下公式對(duì)所述地震道中心樣點(diǎn)的能量權(quán)系數(shù)進(jìn)行歸一化處理：
[0055] 其中，為第k道在i。時(shí)間樣點(diǎn)處的能量權(quán)系數(shù)，Wje)為第j道在i。時(shí)間樣點(diǎn) 處的能量權(quán)系數(shù)，〇}ie)為第k道在i。時(shí)間樣點(diǎn)處的歸一化后的能量權(quán)系數(shù)，i。為當(dāng)前滑動(dòng) 窗中心采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。
[0056] 另一方面，本申請(qǐng)實(shí)施例還提供了一種地震記錄規(guī)則干擾波的提取裝置，所述裝 置包括：
[0057] 獲取單元，用于獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄；
[0058]時(shí)移排齊處理單元，用于對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí) 移排齊處理后的地震記錄；
[0059] 計(jì)算單元，用于計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng) 的相似權(quán)系數(shù)和能量權(quán)系數(shù)；
[0060] 加權(quán)疊加處理單元，用于根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量 權(quán)系數(shù)，對(duì)所述時(shí)移排齊處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排 齊處理后的地震記錄中心道對(duì)應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干擾 波；
[0061] 剩余地震記錄獲取單元，用于根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩 余地震記錄;其中，所述剩余地震記錄為去除規(guī)則干擾波的地震記錄。
[0062] 本申請(qǐng)實(shí)施例提供了一種地震記錄規(guī)則干擾波的提取方法及裝置，該方法采用多 道相關(guān)的方法計(jì)算獲得滑動(dòng)窗內(nèi)地震道每個(gè)樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，然后再 利用多相關(guān)疊加技術(shù)提取規(guī)則干擾波，從而可以獲得去除規(guī)則干擾波的有效的地震記錄。 現(xiàn)有技術(shù)是利用干擾波和有效波在屬性方面的明顯差異實(shí)現(xiàn)干擾波的提取，本申請(qǐng)實(shí)施例 的方法利用了干擾波自身的波形和能量的相關(guān)性，即使在規(guī)則干擾波和有效波的屬性差異 性很小或者沒有差異時(shí)，也能實(shí)現(xiàn)對(duì)規(guī)則干擾波的提取。本申請(qǐng)實(shí)施例的方法能有效的實(shí) 現(xiàn)了規(guī)則干擾波的提取，為地震資料的解釋奠定了基礎(chǔ)。
【附圖說明】
[0063] 為了更清楚地說明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)性的前提 下，還可以根據(jù)這些附圖獲得其他的附圖。
[0064] 圖1是本申請(qǐng)實(shí)施例的地震記錄規(guī)則干擾波的提取方法的示意圖；
[0065] 圖2是本申請(qǐng)實(shí)施例的地震記錄規(guī)則干擾波的提取裝置的示意圖；
[0066] 圖3是本申請(qǐng)一實(shí)施例的合成疊前記錄的共炮點(diǎn)道集圖；
[0067] 圖4是本申請(qǐng)一實(shí)施例的干擾波時(shí)移排齊后的地震記錄圖；
[0068]圖5是本申請(qǐng)一實(shí)施例的時(shí)移排齊后的地震記錄做鑲邊處理后的圖；
[0069] 圖6是本申請(qǐng)一實(shí)施例的提取的干擾波的圖；
[0070] 圖7是本申請(qǐng)一實(shí)施例的提取的干擾波做反時(shí)移校正后的圖；
[0071] 圖8是本申請(qǐng)一實(shí)施例的去除規(guī)則干擾波后的剩余地震記錄的圖；
[0072] 圖9是本申請(qǐng)一實(shí)施例的野外實(shí)際地震記錄的圖；
[0073] 圖10是本申請(qǐng)一實(shí)施例的實(shí)際地震記錄提取的線性干擾波場(chǎng)的圖；
[0074] 圖11是本申請(qǐng)一實(shí)施例的實(shí)際地震記錄去除線性干擾波場(chǎng)后的剩余地震記錄的 圖。
【具體實(shí)施方式】
[0075] 為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí) 施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施 例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通 技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù) 的范圍。
[0076]下面結(jié)合附圖，對(duì)本申請(qǐng)實(shí)施例的【具體實(shí)施方式】作進(jìn)一步的詳細(xì)說明。
[0077] 參考圖1，本申請(qǐng)實(shí)施例提供了一種地震記錄規(guī)則干擾波的提取方法，該方法包 括：
[0078] S1、獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄。
[0079] 本申請(qǐng)實(shí)施例中，首先預(yù)設(shè)滑動(dòng)窗的大小，即確定滑動(dòng)窗內(nèi)地震道的道數(shù)，然后通 過射線追蹤正演模擬的方法得到合成的疊前地震記錄u r(t)，其中r = l，…，當(dāng)前滑 動(dòng)窗內(nèi)地震記錄道數(shù)，t為地震記錄時(shí)間序號(hào)，t=l, ...，N，N為地震記錄的時(shí)間采樣點(diǎn)數(shù)。 如圖3所示為正演模擬的合成疊前地震記錄的共炮點(diǎn)道集，其中地震記錄道數(shù)為240道，時(shí) 間采樣點(diǎn)數(shù)為2000個(gè)。在本申請(qǐng)實(shí)施例中，疊前地震記錄可以為野外的實(shí)際地震記錄，也可 以為模擬合成的地震記錄，本申請(qǐng)實(shí)施例中并不限定。
[0080] S2、對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排齊處理后的地震 記錄。
[0081] 本申請(qǐng)實(shí)施例中，對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排齊 處理后的地震記錄，包括：
[0082]確定所述疊前地震記錄的干擾波的時(shí)距曲線；
[0083]根據(jù)所述干擾波的時(shí)距曲線對(duì)所述疊前地震記錄進(jìn)行處理，獲得時(shí)移排齊處理后 的地震記錄。
[0084]本申請(qǐng)實(shí)施例中，所述干擾波的時(shí)距曲線可以根據(jù)不同類型的干擾波的特點(diǎn)，采 用如下幾種方式確定：
[0085] (1)如果干擾波出現(xiàn)的時(shí)間與空間位置無關(guān)，則干擾波的時(shí)距曲線T (r)可由時(shí)間 常量T。來確定，即T(r)=Tc，r為地震道的序號(hào)，T。為干擾波在地震記錄上出現(xiàn)的時(shí)間常量； (2)如果干擾波為線性波場(chǎng)（如直達(dá)波），則干擾波的時(shí)距曲線可由線性函數(shù)
來確定，其中X(r)為序號(hào)為r的地震道對(duì)應(yīng)的炮檢距，f為干擾波的視速 度，To為零炮檢距對(duì)應(yīng)的干擾波的出現(xiàn)時(shí)間；（3)如果干擾波為雙曲型波場(chǎng)(如反射波），則 干擾波的時(shí)距曲線可由雙曲函數(shù)
來確定，其中X(r)為序號(hào)為r的地震 道對(duì)應(yīng)的炮檢距，V為干擾波的傳播速度，To為零炮檢距對(duì)應(yīng)的干擾波在地震記錄上出現(xiàn)的 時(shí)間；（4)如果干擾波沒有上述特征，則可以通過拾取每道地震記錄上的干擾波同相軸來確 定干擾波的時(shí)距曲線T(r)，其中r為地震道的序號(hào)。
[0086]根據(jù)確定的干擾波的時(shí)距曲線，對(duì)每道地震記錄ur(t)做時(shí)移，每道地震記錄的時(shí) 移量AT(r) = T(r)-Td，Td為時(shí)間常量，時(shí)移排齊后的干擾波同相軸被校正到該時(shí)間上。在本 申請(qǐng)一實(shí)施例匯總，將Td設(shè)為150,則時(shí)移排齊后的干擾波同相軸被校正到該時(shí)間上，如圖4 為干擾波時(shí)移排齊后的地震記錄。
[0087]在本申請(qǐng)另一實(shí)施例中，為減小邊界效應(yīng)，在步驟S2后，還包括對(duì)時(shí)移排齊處理后 的地震記錄進(jìn)行鑲邊處理。如果采用的滑動(dòng)窗道數(shù)為施，則鑲邊道數(shù)為(Mi-D/2道，也就是 在地震的第一道前面和最后一道后面，分別鑲嵌(MrD/2道。地震記錄的前(MrD/2道關(guān)于 第1道鏡像放置形成前邊的鑲邊地震道;地震記錄的后(Mi-D/2道關(guān)于最后一道鏡像放置 形成后邊的鑲邊地震道，圖5為對(duì)時(shí)移排齊后的地震記錄做鑲邊處理的結(jié)果。
[0088] S3、計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系 數(shù)和能量權(quán)系數(shù)。
[0089] 計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)，包 括：
[0090] 計(jì)算所述時(shí)移排齊處理后的地震記錄計(jì)算中地震道的相關(guān)系數(shù)；
[0091] 根據(jù)所述地震道的相關(guān)系數(shù)計(jì)算所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)；
[0092] 對(duì)所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理 后的相似權(quán)系數(shù)計(jì)算所述地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)。
[0093] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算所述地震道的相關(guān)系數(shù)：
[0094] CLi： = ^uk(/) ? UjO'), k, j - 1,- ? A/,
[0095]其中，k和j是滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是滑動(dòng)窗內(nèi)地震記錄道數(shù)，i是滑動(dòng)窗內(nèi) 的時(shí)間樣點(diǎn)序號(hào)，ib和分別是當(dāng)前滑動(dòng)窗第一個(gè)和最后一個(gè)時(shí)間樣點(diǎn)序號(hào)，uk(i)為第k道 第i個(gè)時(shí)間樣點(diǎn)的振幅值，^⑴為第j道第i個(gè)時(shí)間樣點(diǎn)的振幅值，Cu為第k道和第j道地震 記錄的相關(guān)系數(shù)。
[0096]本申請(qǐng)實(shí)施例中，在獲得地震道的相關(guān)系數(shù)后，可以采用以下公式計(jì)算所述地震 道中心采樣點(diǎn)的相似權(quán)系數(shù)：
[0098] 其中，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)，Ckj為第k道和 第j道的互相關(guān)系數(shù)，Cw和Ckk分別為第j道和第k道的自相關(guān)系數(shù)，i。為當(dāng)前滑動(dòng)窗中心采樣 點(diǎn)的時(shí)間序號(hào)，wjl ；)為第k道在i。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)。
[0099] 本申請(qǐng)實(shí)施例中，可以采用以下公式對(duì)所述地震道中心采樣點(diǎn)的相似權(quán)系數(shù)進(jìn)行 歸一化處理：
[0101] 其中，wjb為第k道在i。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)，1為第j道在i。時(shí)間樣點(diǎn) 處的相似權(quán)系數(shù)，為第k道在i。時(shí)間樣點(diǎn)處的歸一化后的相似權(quán)系數(shù)，i。為當(dāng)前滑動(dòng) 窗中心采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。
[0102] 本申請(qǐng)實(shí)施例中，在時(shí)間方向上，滑動(dòng)窗按照給定的時(shí)間步長滑動(dòng)，逐個(gè)計(jì)算每一 個(gè)滑動(dòng)窗內(nèi)每一道時(shí)間中心樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)，在兩個(gè)時(shí)窗的時(shí)間中心樣點(diǎn)之間，利 用線性插值計(jì)算每個(gè)樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)，其中，t為時(shí)間樣點(diǎn)序號(hào)，t = l,…，N，N 為每道地震記錄的時(shí)間采樣點(diǎn)數(shù)，k是當(dāng)前滑動(dòng)窗內(nèi)地震道的序號(hào)，k=l，…％屬是當(dāng)前滑 動(dòng)窗內(nèi)地震道記錄道數(shù)。
[0103] 本申請(qǐng)實(shí)施例中，計(jì)算所述時(shí)移排齊處理后的地震記錄計(jì)算中地震道每個(gè)采樣點(diǎn) 對(duì)應(yīng)的能量權(quán)系數(shù)，包括：
[0104] 根據(jù)所述時(shí)移排齊處理后的地震記錄獲取能量濾波參數(shù)；
[0105] 根據(jù)所述能量濾波參數(shù)獲得能量濾波算子；
[0106] 根據(jù)所述能量濾波算子計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道中心樣點(diǎn) 的能量權(quán)系數(shù)；
[0107] 對(duì)所述地震道中心采樣點(diǎn)的能量權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理 后的能量權(quán)系數(shù)計(jì)算所述地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)。
[0108] 本申請(qǐng)實(shí)施例中，根據(jù)所述時(shí)移排齊處理后的地震記錄獲取能量濾波參數(shù)，包括：
[0109] 計(jì)算所述時(shí)移排齊處理后的地震記錄中每一地震道的能量；
[0110] 根據(jù)預(yù)設(shè)的能量比例閾值，將所述時(shí)移排齊處理后的地震記錄中能量小于閾值的 地震道去除；
[0111] 計(jì)算所述時(shí)移排齊處理后的地震記錄中剩余地震道的平均能量值，并根據(jù)所述每 一地震道的能量與所述平均能量值計(jì)算所述時(shí)移排齊處理后的地震記錄的能量標(biāo)準(zhǔn)偏差；
[0112] 根據(jù)所述能量標(biāo)準(zhǔn)偏差獲取能量濾波參數(shù)。
[0113] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算所述滑動(dòng)窗的能量標(biāo)準(zhǔn)偏差％，
[0115] 其中，< 表示地震道的能量，1 = 1，~沁，111是當(dāng)前滑動(dòng)窗內(nèi)地震道記錄道數(shù)，如表 示平均能量值。
[0116] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算所述能量濾波參數(shù)： ^1 - _ 3 CT 「 n -1-5 <j
[0117] - T c " ei^em+\.5l-(je e4=em + 3-cre
[0118] 其中，ei、e2、e3和e4表示能量濾波參數(shù)，em表示平均能量值，〇e表示能量標(biāo)準(zhǔn)偏差。
[0119] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算所述能量濾波算子： 1 e2<ec <e3
[0120] ^ ' 0.5- ()_5- cos(;r(e.」e4)/(，,- e4)) e-; < ec < e_t 、 0 ec <el5ec >e4
[0121] 其中，f (ec^em^e)表示能量濾波算子，ei、e2、e3和e4表示能量濾波參數(shù)，e。表示當(dāng)前 地震道的能量，~表示能量標(biāo)準(zhǔn)偏差。
[0122] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算所述滑動(dòng)窗內(nèi)地震道中心樣點(diǎn)的能量權(quán) 系數(shù)：
[0123] V^:.)=/(d〇；.) k = ~-(k.)
[0124] 其中，為第k道在ic時(shí)間樣點(diǎn)處的能量權(quán)系數(shù)，f(ec;em，〇e)表示能量濾波算 子，i。為當(dāng)前滑動(dòng)窗中心樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地 震記錄道數(shù)。
[0125] 本申請(qǐng)實(shí)施例中，可以采用以下公式對(duì)所述地震道中心采樣點(diǎn)的能量權(quán)系數(shù)進(jìn)行 歸一化處理：
[0127] 其中，Wefe)為第k道在i。時(shí)間樣點(diǎn)處的能量權(quán)系數(shù)，W^}為第j道在i。時(shí)間樣點(diǎn) 處的能量權(quán)系數(shù)，^為第k道在i。時(shí)間樣點(diǎn)處的歸一化后的能量權(quán)系數(shù)，i。為當(dāng)前滑動(dòng) 窗中心采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。
[0128] 本申請(qǐng)實(shí)施例中，在時(shí)間方向上，滑動(dòng)窗按照給定的時(shí)間步長滑動(dòng)，逐個(gè)計(jì)算每一 個(gè)滑動(dòng)窗內(nèi)每一道時(shí)間中心樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)，在兩個(gè)時(shí)窗的時(shí)間中心樣點(diǎn)之間，利 用線性插值方法計(jì)算每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)其中，t為時(shí)間樣點(diǎn)序號(hào)，t = 1，…，N，N為每道地震記錄的時(shí)間采樣點(diǎn)數(shù)，k是當(dāng)前滑動(dòng)窗內(nèi)地震道的序號(hào)， 是當(dāng)前滑動(dòng)窗內(nèi)地震道總道數(shù)。
[0129] S4、根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，對(duì)所述時(shí)移 排齊處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排齊處理后的地震記錄 中心道對(duì)應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干擾波。
[0130] 本申請(qǐng)實(shí)施例中，可以采用以下計(jì)算公式獲得所述時(shí)移排齊處理后地震記錄中心 道對(duì)應(yīng)的疊加地震道： ik) ik)
[0131] s(r) = (，)價(jià):(r)%,(，）t=-\，…N k=l
[0132] 其中，s(t)表示疊加地震道；表示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)； 表示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù);Uk(t)表示時(shí)移排齊處理后的地震記錄上 第k道地震道在時(shí)間樣點(diǎn)序號(hào)為t時(shí)的振幅值;k表示地震道的序號(hào);Mi表示滑動(dòng)窗內(nèi)地震記 錄的道數(shù);t表示時(shí)間樣點(diǎn)序號(hào)，N為每道地震記錄的時(shí)間采樣點(diǎn)數(shù)。
[0133] 根據(jù)地震波理論可知，地震記錄來自同一激發(fā)源的規(guī)則干擾波，除了在視速度、偏 振特性和頻率等屬性上具有相關(guān)性外，不同道上的干擾波的波形和能量也具有相關(guān)性。本 申請(qǐng)實(shí)施例根據(jù)規(guī)則干擾波的這一特點(diǎn)，采用多道相關(guān)的方法計(jì)算出每道的波形相似權(quán)系 數(shù)和能量權(quán)系數(shù)，然后應(yīng)用多相關(guān)加權(quán)疊加技術(shù)從而可以有效的提取規(guī)則干擾波。如圖6所 示為本申請(qǐng)一實(shí)施例得到的疊加地震記錄，即需要提取的干擾波。
[0134] S5、根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄;其中，所述剩 余地震記錄為去除規(guī)則干擾波的地震記錄。
[0135] 本申請(qǐng)實(shí)施例中，根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄， 包括：
[0136] 對(duì)所述疊加地震道進(jìn)行反時(shí)移校正處理，獲得反時(shí)移校正處理后的疊加地震道；
[0137] 根據(jù)所述疊前地震記錄和所述反時(shí)移校正處理后的疊加地震道，獲取振幅匹配系 數(shù)；
[0138] 根據(jù)所述振幅匹配系數(shù)對(duì)所述疊加地震道進(jìn)行振幅匹配處理，利用所述疊前地震 記錄減去振幅匹配處理后的疊加地震道，獲得剩余地震記錄。
[0139] 本申請(qǐng)實(shí)施例中，可以根據(jù)步驟S2中確定的干擾波的時(shí)距曲線，對(duì)所示疊加地震 道進(jìn)行反時(shí)移校正，獲得反時(shí)移校正處理后的疊加地震道。其中，時(shí)移量為A T(r)=Td_T (r)，本申請(qǐng)實(shí)施例中，Td為150,如圖7所示為提取的干擾波做反時(shí)移校正后的結(jié)果。
[0140] 本申請(qǐng)實(shí)施例中，可以采用以下公式計(jì)算振幅匹配系數(shù)：
[0141] ?；/ /乞'.， i =k / i-=ib
[0142] 其中，a」為第j道的振幅匹配系數(shù)，Ujl為第j道疊前地震記錄第i個(gè)采樣點(diǎn)的振幅 值，Sl為第j道地震記錄對(duì)應(yīng)的疊加地震道第i個(gè)采樣點(diǎn)的振幅值，i b和h為振幅計(jì)算時(shí)窗的 起止樣點(diǎn)序號(hào)。
[0143] 本申請(qǐng)實(shí)施例中，根據(jù)所述振幅匹配系數(shù)對(duì)所述疊加地震道進(jìn)行振幅匹配處理， 然后用疊前地震記錄減去振幅匹配處理后的疊加地震道，得到剩余地震記錄，可以采用以 下公式計(jì)算剩余地震記錄：
[0144] rj(t) =uj(t)-a ? s(t),tb^t^te
[0145] 其中，a表示振幅匹配系數(shù)，s(t)表示疊加地震道，Uj(t)表示疊前地震記錄，〇(t) 表示剩余地震記錄，j是滑動(dòng)窗內(nèi)地震道的序號(hào)，t b和te為地震記錄相減應(yīng)用時(shí)窗的起止時(shí) 間。這樣，就實(shí)現(xiàn)了將規(guī)則干擾波從疊前地震記錄中分離出來的目的。如圖8所述為將規(guī)則 干擾波去除后的剩余地震記錄。
[0146] 在本申請(qǐng)一實(shí)施例中，在步驟S5之后，還包括：
[0147] 逐步減小所述滑動(dòng)窗的大小，對(duì)所述剩余地震記錄進(jìn)行規(guī)則干擾波去除處理。
[0148] 本申請(qǐng)實(shí)施例中干擾波的提取是在預(yù)設(shè)的滑動(dòng)窗內(nèi)進(jìn)行的，在開始進(jìn)行干擾波提 取時(shí)可以采用較大的滑動(dòng)窗，從而將疊前地震記錄上的干擾波的主要能量提取出來，然后 再逐步減小滑動(dòng)窗的大小，將剩余地震記錄中的干擾波的殘余能量逐步地提取出來，從而 得到質(zhì)量更好的有效波。如圖9所示為野外實(shí)際地震記錄，圖10所示為本申請(qǐng)實(shí)施例的方法 提取的線性干擾波場(chǎng)，如圖11所示為本申請(qǐng)實(shí)施例的方法得到的去除線性干擾波后的剩余 地震記錄，從圖中可以看出，本申請(qǐng)實(shí)施例可以有效的提取干擾波，從而為地震資料解釋奠 定基礎(chǔ)。
[0149] 另一方面，參考圖2,本申請(qǐng)實(shí)施例還提供了一種地震記錄規(guī)則干擾波的提取裝 置，所述裝置包括：
[0150] 獲取單元21，用于獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄；
[0151]時(shí)移排齊處理單元22,用于對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取 時(shí)移排齊處理后的地震記錄；
[0152]計(jì)算單元23,用于計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì) 應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)；
[0153]加權(quán)疊加處理單元24,用于根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能 量權(quán)系數(shù)，對(duì)所述時(shí)移排齊處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移 排齊處理后的地震記錄中心道對(duì)應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干 擾波；
[0154] 剩余地震記錄獲取單元25,用于根據(jù)所述疊前地震記錄和所述疊加地震道，獲取 剩余地震記錄;其中，所述剩余地震記錄為去除規(guī)則干擾波的地震記錄。
[0155] 本實(shí)施例的裝置的各組成部分分別用于實(shí)現(xiàn)前述實(shí)施例的方法的各步驟，由于在 方法實(shí)施例中，已經(jīng)對(duì)各步驟進(jìn)行了詳細(xì)說明，在此不再贅述。
[0156] 本申請(qǐng)實(shí)施例提供了一種地震記錄規(guī)則干擾波的提取方法及裝置，該方法采用多 道相關(guān)的方法計(jì)算獲得滑動(dòng)窗內(nèi)地震道每個(gè)樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，然后再 利用多相關(guān)疊加技術(shù)及振幅平滑相減技術(shù)提取規(guī)則干擾波，從而獲得去除規(guī)則干擾波的有 效的地震記錄?，F(xiàn)有技術(shù)是利用規(guī)則干擾波和有效波在屬性方面的差異性實(shí)現(xiàn)對(duì)規(guī)則干擾 波的提取，當(dāng)規(guī)則干擾波和有效波的屬性差異很小或沒有差異時(shí)，現(xiàn)有技術(shù)的方法就不能 對(duì)干擾波進(jìn)行有效提取或壓制;而本申請(qǐng)實(shí)施例的方法是利用了干擾波自身的波形和能量 的相關(guān)性，因此即使規(guī)則干擾波和有效波的屬性差異很小或沒有差異，也可以利用干擾波 自身波形和能量的相關(guān)性提取出規(guī)則干擾波。本申請(qǐng)實(shí)施例的方法能有效的實(shí)現(xiàn)對(duì)干擾波 的提取，從而為地震資料的解釋奠定了基礎(chǔ)。
[0157] 在一個(gè)或多個(gè)示例性的設(shè)計(jì)中，本申請(qǐng)實(shí)施例所描述的上述功能可以在硬件、軟 件、固件或這三者的任意組合來實(shí)現(xiàn)。如果在軟件中實(shí)現(xiàn)，這些功能可以存儲(chǔ)與電腦可讀的 媒介上，或以一個(gè)或多個(gè)指令或代碼形式傳輸于電腦可讀的媒介上。電腦可讀媒介包括電 腦存儲(chǔ)媒介和便于使得讓電腦程序從一個(gè)地方轉(zhuǎn)移到其它地方的通信媒介。存儲(chǔ)媒介可以 是任何通用或特殊電腦可以接入訪問的可用媒體。例如，這樣的電腦可讀媒體可以包括但 不限于1^11、1?(觀、￡￡?1?(通、^-1?(通或其它光盤存儲(chǔ)、磁盤存儲(chǔ)或其它磁性存儲(chǔ)裝置，或其它任 何可以用于承載或存儲(chǔ)以指令或數(shù)據(jù)結(jié)構(gòu)和其它可被通用或特殊電腦、或通用或特殊處理 器讀取形式的程序代碼的媒介。
[0158] 以上所述的具體實(shí)施例，對(duì)本申請(qǐng)的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳 細(xì)說明，所應(yīng)理解的是，以上所述僅為本申請(qǐng)實(shí)施例的具體實(shí)施例而已，并不用于限定本申 請(qǐng)的保護(hù)范圍，凡在本申請(qǐng)的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包 含在本申請(qǐng)的保護(hù)范圍之內(nèi)。
【主權(quán)項(xiàng)】
1. 一種地震記錄規(guī)則干擾波的提取方法，其特征在于，所述方法包括： 獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄； 對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排齊處理后的地震記錄； 計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能 量權(quán)系數(shù)； 根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系數(shù)，對(duì)所述時(shí)移排齊處理 后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排齊處理后的地震記錄中屯、道對(duì) 應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干擾波； 根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地震記錄;其中，所述剩余地震記 錄為去除規(guī)則干擾波的地震記錄。2. 如權(quán)利要求1所述的方法，其特征在于，在根據(jù)所述疊前地震記錄和所述疊加地震 道，獲取剩余地震記錄之后，所述方法還包括： 逐步減小所述滑動(dòng)窗的大小，對(duì)所述剩余地震記錄進(jìn)行規(guī)則干擾波提取處理。3. 如權(quán)利要求1所述的方法，其特征在于，采用W下計(jì)算公式獲得所述時(shí)移排齊處理后 地震記錄中屯、道對(duì)應(yīng)的疊加地震道：其中，s(t)表示疊加地震道;表示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù);表 示地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù);Uk(t)表示時(shí)移排齊處理后的地震記錄上第k道地 震道在時(shí)間樣點(diǎn)序號(hào)為t時(shí)的振幅值;k表示地震道的序號(hào);Mi表示滑動(dòng)窗內(nèi)地震記錄的道 數(shù);t表示時(shí)間樣點(diǎn)序號(hào)，N為每道地震記錄的時(shí)間采樣點(diǎn)數(shù)。4. 如權(quán)利要求1所述的方法，其特征在于，根據(jù)所述疊前地震記錄和所述疊加地震道， 獲取剩余地震記錄，包括： 對(duì)所述疊加地震道進(jìn)行反時(shí)移校正處理，獲得反時(shí)移校正處理后的疊加地震道； 根據(jù)所述疊前地震記錄和所述反時(shí)移校正處理后的疊加地震道，獲取振幅匹配系數(shù)； 根據(jù)所述振幅匹配系數(shù)對(duì)所述疊加地震道進(jìn)行振幅匹配處理，利用所述疊前地震記錄 減去振幅匹配處理后的疊加地震道，獲得剩余地震記錄。5. 如權(quán)利要求4所述的方法，S賠化#干.要田W下/A式計(jì)算振幅匹配系數(shù)：其中，W為第j道的振幅匹配系數(shù)，UW為第j道疊前地震記錄第i個(gè)采樣點(diǎn)的振幅值，Si為 第j道地震記錄對(duì)應(yīng)的疊加地震道第i個(gè)采樣點(diǎn)的振幅值，ib和ie為振幅計(jì)算時(shí)窗的起止樣 點(diǎn)序號(hào)。6. 如權(quán)利要求1所述的方法，其特征在于，計(jì)算所述時(shí)移排齊處理后的地震記錄中地震 道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)，包括： 計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的相關(guān)系數(shù)； 根據(jù)所述地震道的相關(guān)系數(shù)計(jì)算所述地震道中屯、采樣點(diǎn)的相似權(quán)系數(shù)； 對(duì)所述地震道中屯、采樣點(diǎn)的相似權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理后的 相似權(quán)系數(shù)計(jì)算所述地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)。7. 如權(quán)利要求6所述的方法，其特征在于，采用W下公式計(jì)算所述地震道的相關(guān)系數(shù)：其中，k和j是滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是滑動(dòng)窗內(nèi)地震記錄道數(shù)，i是滑動(dòng)窗內(nèi)的時(shí)間 樣點(diǎn)序號(hào)，ib和ie分別是當(dāng)前滑動(dòng)窗第一個(gè)和最后一個(gè)時(shí)間樣點(diǎn)序號(hào)，Uk(i)為第k道第i個(gè) 時(shí)間樣點(diǎn)的振幅值，11如）為第j道第i個(gè)時(shí)間樣點(diǎn)的振幅值，Cm為第k道和第j道地震記錄 的相關(guān)系數(shù)。8. 如權(quán)利要求6所述的方法，其特征在于，采用W下公式計(jì)算所述地震道中屯、采樣點(diǎn)的 相似權(quán)系數(shù)：其中，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)，CkJ為第k道和第j 道的互相關(guān)系數(shù)，Cw和Ckk分別為第j道和第k道的自相關(guān)系數(shù)，ic為當(dāng)前滑動(dòng)窗中屯、采樣點(diǎn) 的時(shí)間序號(hào)，^\^]^為第4道在1。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)。9. 如權(quán)利要求6所述的方法，其特征在于，采用W下公式對(duì)所述地震道中屯、采樣點(diǎn)的相 似權(quán)系數(shù)進(jìn)行歸一化處理：其中，W:策：)為第k道在1。時(shí)間樣點(diǎn)處的相似權(quán)系數(shù)，\滿。；> 為第j道在1。時(shí)間樣點(diǎn)處的 相似權(quán)系數(shù)，為第k道在i。時(shí)間樣點(diǎn)處的歸一化后的相似權(quán)系數(shù)，i。為當(dāng)前滑動(dòng)窗中 屯、采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。10. 如權(quán)利要求1所述的方法，其特征在于，計(jì)算所述時(shí)移排齊處理后的地震記錄中地 震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)，包括： 根據(jù)所述時(shí)移排齊處理后的地震記錄獲取能量濾波參數(shù)； 根據(jù)所述能量濾波參數(shù)獲得能量濾波算子； 根據(jù)所述能量濾波算子計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道中屯、采樣點(diǎn)的 能量權(quán)系數(shù)； 對(duì)所述地震道中屯、采樣點(diǎn)的能量權(quán)系數(shù)進(jìn)行歸一化處理，并根據(jù)所述歸一化處理后的 能量權(quán)系數(shù)計(jì)算所述地震道每個(gè)采樣點(diǎn)對(duì)應(yīng)的能量權(quán)系數(shù)。11. 如權(quán)利要求10所述的方法，其特征在于，根據(jù)所述時(shí)移排齊處理后的地震記錄獲取 能量濾波參數(shù)，包括： 計(jì)算所述時(shí)移排齊處理后的地震記錄中每一地震道的能量； 根據(jù)預(yù)設(shè)的能量比例闊值，將所述時(shí)移排齊處理后的地震記錄中能量小于闊值的地震 道去除； 計(jì)算所述時(shí)移排齊處理后的地震記錄中剩余地震道的平均能量值，并根據(jù)所述每一地 震道的能量與所述平均能量值計(jì)算所述時(shí)移排齊處理后的地震記錄的能量標(biāo)準(zhǔn)偏差； 根據(jù)所述能量標(biāo)準(zhǔn)偏差獲取能量濾波參數(shù)。12.如權(quán)利要求11所述的方法，其據(jù)佈在干，采巧Pi下公式計(jì)算所述能量濾波參數(shù)： 其中，ei、62、63和64表示能量濾波參數(shù)，em表示平均能量值，Oe表示能量標(biāo)準(zhǔn)偏差。 13 .如權(quán)利要求10所述的方法，其特征在于，采用W下公式計(jì)算所述能量濾波算子：其中，ei、62、63和64表示能量濾波參數(shù)，Gm表示平均能量值，Gc表示地震道能量值。14. 如權(quán)利要求10所述的方法，其特征在于，采用W下公式對(duì)所述地震道中屯、采樣點(diǎn)的 能量權(quán)系數(shù)進(jìn)行歸一化處理：其中，為第k道在1。時(shí)間樣點(diǎn)處的能量權(quán)系數(shù)，為第j道在1。時(shí)間樣點(diǎn)處的 能量權(quán)系數(shù)，\\5；^；)為第4道在1。時(shí)間樣點(diǎn)處的歸一化后的能量權(quán)系數(shù)，1。為當(dāng)前滑動(dòng)窗中 屯、采樣點(diǎn)的時(shí)間序號(hào)，k為滑動(dòng)窗內(nèi)地震道的序號(hào)，Mi是當(dāng)前滑動(dòng)窗內(nèi)地震記錄道數(shù)。15. -種地震記錄規(guī)則干擾波的提取裝置，其特征在于，所述裝置包括： 獲取單元，用于獲取預(yù)設(shè)滑動(dòng)窗內(nèi)的疊前地震記錄； 時(shí)移排齊處理單元，用于對(duì)所述疊前地震記錄進(jìn)行干擾波時(shí)移排齊處理，獲取時(shí)移排 齊處理后的地震記錄； 計(jì)算單元，用于計(jì)算所述時(shí)移排齊處理后的地震記錄中地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相 似權(quán)系數(shù)和能量權(quán)系數(shù)； 加權(quán)疊加處理單元，用于根據(jù)所述地震道的每個(gè)采樣點(diǎn)對(duì)應(yīng)的相似權(quán)系數(shù)和能量權(quán)系 數(shù)，對(duì)所述時(shí)移排齊處理后的地震記錄中地震道進(jìn)行加權(quán)疊加處理，獲得所述時(shí)移排齊處 理后的地震記錄中屯、道對(duì)應(yīng)的疊加地震道;其中，所述疊加地震道為提取的規(guī)則干擾波； 剩余地震記錄獲取單元，用于根據(jù)所述疊前地震記錄和所述疊加地震道，獲取剩余地 震記錄;其中，所述剩余地震記錄為去除規(guī)則干擾波的地震記錄。
【文檔編號(hào)】G01V1/28GK105911585SQ201610522410
【公開日】2016年8月31日
【申請(qǐng)日】2016年7月5日
【發(fā)明人】張文波, 李合群, 張文棟
【申請(qǐng)人】中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司, 中油油氣勘探軟件國家工程研究中心有限公司