于 杰

疊前去噪技術(shù)在煤礦采區(qū)全數(shù)字高密度三維地震中的應(yīng)用

于 杰

(中國(guó)礦業(yè)大學(xué) 資源與地球科學(xué)學(xué)院，江蘇 徐州 221116)

高密度三維地震采用數(shù)字檢波器單點(diǎn)接收，數(shù)字檢波器靈敏度高，對(duì)信號(hào)采取寬頻接收，使得更多弱背景噪聲和有效信號(hào)一起被記錄下來(lái)。因此，全數(shù)字高密度三維地震資料單炮信噪比偏低，需要在疊前采用針對(duì)性去噪技術(shù)以提高處理效果。以淮北礦區(qū)全數(shù)字高密度三維地震數(shù)據(jù)為例，討論了高密度空間采樣的充分性以及高密度三維地震數(shù)據(jù)中面波、異常振幅與多次波等主要干擾波的特點(diǎn)，結(jié)合高密度三維地震數(shù)據(jù)寬方位、寬頻帶、高密度等特點(diǎn)，針對(duì)性地壓制面波、異常振幅與多次波。在縱測(cè)線和非縱測(cè)線方向上對(duì)面波進(jìn)行十字排列分析，根據(jù)面波空間分布規(guī)律設(shè)計(jì)錐形濾波器進(jìn)行濾除；通過(guò)計(jì)算多個(gè)地震道的包絡(luò)并設(shè)定閾值自動(dòng)識(shí)別后進(jìn)行單道壓制以處理異常振幅。應(yīng)用拋物Radon變換把CMP道集變換到-域剔除多次波。原始數(shù)據(jù)經(jīng)過(guò)處理后，有效信號(hào)得以保留，信噪比得到較大提升，取得了良好的效果，對(duì)今后的全數(shù)字高密度三維地震數(shù)據(jù)處理工作有借鑒意義。

高密度三維地震；干擾波；十字排列；自動(dòng)識(shí)別；Radon變換

在地震勘探中，由于采集過(guò)程中的環(huán)境條件、施工因素和儀器設(shè)備等影響，導(dǎo)致地震記錄不同程度地存在各種噪聲干擾，常使有效信號(hào)能量被強(qiáng)能量的干擾波湮沒(méi)[1]，原始記錄信噪比較低。近年來(lái)，全數(shù)字高密度三維地震在煤田地震勘探中得到推廣[2]。全數(shù)字高密度三維地震采用數(shù)字檢波器單點(diǎn)無(wú)組合接收，能夠最大限度地保護(hù)信號(hào)頻率成分，但地震記錄中的噪聲也顯著增強(qiáng)，信噪比降低，而高信噪比是實(shí)現(xiàn)高分辨率和高保真度的基礎(chǔ)[3]，因此，壓制噪聲、提高信噪比是全數(shù)字高密度地震資料處理中的重要環(huán)節(jié)。

淮北礦區(qū)全數(shù)字高密度地震資料中干擾波主要有面波、異常振幅干擾及多次波[4]。由于數(shù)字檢波器單點(diǎn)接收對(duì)低頻響應(yīng)良好，使得地震數(shù)據(jù)中面波普遍較強(qiáng)，而面波頻率低，一般為20 Hz以下，與反射波低頻信號(hào)重疊。當(dāng)前應(yīng)用最為廣泛的壓制面波的方法是二維F-K濾波方法。王金龍等[5]將三維錐形濾波與二維濾波在頻率–太波數(shù)域和時(shí)間–太空間域進(jìn)行對(duì)比，認(rèn)為三維錐形濾波的優(yōu)勢(shì)在于對(duì)遠(yuǎn)端排列的處理能力比二維濾波更強(qiáng)，三維濾波克服了二維濾波僅能壓制有限區(qū)域內(nèi)噪聲的缺點(diǎn)，避免了二維濾波產(chǎn)生的蚯蚓化現(xiàn)象，去噪能力與二維濾波相比有大幅度提升；李文花[6]在黃土塬地區(qū)應(yīng)用錐形濾波方法壓制面波，取得了良好的效果。但是，三維濾波要求數(shù)據(jù)能夠抽成十字排列，在野外施工中，由于地表障礙物的存在會(huì)使得炮點(diǎn)或檢波點(diǎn)偏離原設(shè)計(jì)位置，金丹提出三維錐形濾波法在野外施工情況下炮點(diǎn)和檢波點(diǎn)偏離設(shè)計(jì)位置時(shí)進(jìn)行觀測(cè)系統(tǒng)規(guī)則化的方法[7]，但對(duì)于變觀較多的區(qū)域該方法應(yīng)用效果不佳。

除面波以外，淮北礦區(qū)地表障礙物極多，村莊、工廠、公路、河道等障礙物極易產(chǎn)生強(qiáng)異常振幅噪聲。與有效反射波相比，異常能量往往大于有效反射能量，大部分異常振幅噪聲在空間和時(shí)間上具有隨機(jī)性特征，其能量強(qiáng)，頻帶寬[8]，且與有效波在空間和時(shí)間上重疊而難以區(qū)分。強(qiáng)異常噪聲的存在導(dǎo)致無(wú)法準(zhǔn)確計(jì)算振幅補(bǔ)償因子，反褶積處理無(wú)法獲取準(zhǔn)確子波、偏移及DMO處理畫(huà)弧現(xiàn)象等，嚴(yán)重影響處理效果。對(duì)于異常振幅噪聲，劉洋等提出一種改進(jìn)的時(shí)變窗口中值濾波方法，通過(guò)選取合理的閾值使中值濾波去能夠在消除隨機(jī)噪聲時(shí)采用較大濾波窗口進(jìn)行處理，在處理有效信號(hào)時(shí)采用小濾波窗口進(jìn)行處理既能有效消除隨機(jī)噪聲，又能最大限度保護(hù)有用信號(hào)[9]；刁塑等[10]在常規(guī)曲波變換方法基礎(chǔ)上提出了局部自適應(yīng)曲波變換去噪方法，在不同尺度下設(shè)置不同閾值，壓制每個(gè)尺度內(nèi)的噪聲；陳志德等[11]提出了在不同頻段的數(shù)據(jù)子集上利用噪聲在不同頻段與有效反射信號(hào)的振幅差異進(jìn)行去噪的噪聲分頻壓制技術(shù)。

在衰減多次波的諸多方法中，Radon變換以理論簡(jiǎn)單、計(jì)算效率高、壓制效果突出等優(yōu)點(diǎn)，一直是衰減多次波的常用方法。熊登等[12]提出結(jié)合時(shí)間域高分辨率和頻率域高計(jì)算效率計(jì)算的優(yōu)點(diǎn)進(jìn)行混合域計(jì)算的方法，減小了采集孔徑效應(yīng)引起的近偏移距假象，能夠更好地聚焦一次波，壓制多次波。

與傳統(tǒng)三維地震相比，高密度三維地震數(shù)據(jù)具有全數(shù)字、寬方位、高密度、多偏移距和海量數(shù)據(jù)等特點(diǎn)，地震數(shù)據(jù)的炮檢距與方位角信息更加豐富，在傳統(tǒng)去噪方法的基礎(chǔ)上，開(kāi)發(fā)新的去噪技術(shù)是十分必要的[13-14]?；诖?，筆者總結(jié)了適用于全數(shù)字高密度三維地震數(shù)據(jù)的壓制面波、異常振幅與多次波的去噪方法，并將其應(yīng)用于淮北礦區(qū)的地震資料中，取得了較好的效果。

1 工區(qū)概況

淮北礦區(qū)位于淮北平原，地勢(shì)平坦，地表村莊、工廠建筑物密集、面積大、居民多，村莊周邊分布大量養(yǎng)殖場(chǎng)等障礙物，大型河流眾多，小型溝渠縱橫發(fā)達(dá)，公路鐵路交錯(cuò)，能通行汽車的土公路較多，局部地段對(duì)地震波的激發(fā)和接收會(huì)產(chǎn)生負(fù)面影響。礦區(qū)內(nèi)新生界地層厚度較大(200~600 m)，對(duì)地震波的高頻成分吸收衰減嚴(yán)重，地震資料的分辨率偏低，多次波干擾嚴(yán)重；深層反射的信噪比低，反射系數(shù)低、球面擴(kuò)散與大地濾波導(dǎo)致高頻衰減嚴(yán)重。典型單炮記錄如圖1所示，可見(jiàn)工區(qū)內(nèi)面波普遍發(fā)育，并且異常振幅明顯。

針對(duì)淮北礦區(qū)地震地質(zhì)條件，在分析全數(shù)字高密度三維地震原始資料的基礎(chǔ)上，認(rèn)為淮北礦區(qū)全數(shù)字高密度地震數(shù)據(jù)含有的主要干擾波是面波、異常振幅以及多次波干擾。

2 疊前面波衰減方法

高密度三維地震采用數(shù)字檢波器單點(diǎn)接收，單炮記錄的面波干擾能量強(qiáng)、視速度低、頻率低、分布范圍大，并與目的層反射波大范圍重疊。工區(qū)內(nèi)面波視速度普遍低于600 m/s，能量主要集中在20 Hz以下，且7 Hz能量最強(qiáng)，超過(guò)10 Hz后面波能量迅速降低，20 Hz以上面波能量趨于零。這與高密度三維地震保護(hù)弱信號(hào)、保護(hù)寬頻的初衷相違背；F-K濾波能夠有效壓制線性干擾，但是面波在單炮記錄上不完全呈線性特征，尤其是在遠(yuǎn)離炮點(diǎn)的排列上往往呈現(xiàn)出雙曲線特征，利用去除線性噪聲的F-K濾波難以達(dá)到預(yù)期的效果，而且在去除面波的同時(shí)在一定程度上損傷了有效波。目前，在淮北礦區(qū)全數(shù)字高密度地震數(shù)據(jù)處理中，最常用的面波衰減方法有自適應(yīng)面波衰減法和十字交叉排列去噪技術(shù)2種。

圖1 研究區(qū)內(nèi)為原始單炮記錄

2.1 自適應(yīng)面波衰減

自適應(yīng)面波衰減技術(shù)利用時(shí)頻分析的方法，根據(jù)面波和反射波在頻率分布特征、空間分布范圍、能量等方面的差異，先檢測(cè)出面波在時(shí)間和空間上的分布范圍，以確定面波能量的頻率分布特征，并根據(jù)這種特征對(duì)其進(jìn)行加權(quán)壓制。該方法只壓制面波干擾，對(duì)有效反射損傷小。

2.2 十字交叉排列去噪

十字交叉排列去噪技術(shù)是壓制面波干擾的新技術(shù)，其核心思想是打破傳統(tǒng)二維線的約束，將一條炮線與一條檢波線組合為一個(gè)十字排列，通過(guò)對(duì)十字排列道集中面波的分布特征來(lái)設(shè)計(jì)相應(yīng)的濾波器，達(dá)到從三維上壓制面波的目的。處理時(shí)將炮線和接收線號(hào)組合產(chǎn)生新的道頭字，按照該道頭字對(duì)數(shù)據(jù)進(jìn)行抽道，形成由橫向的紅色炮線與縱向的檢波線共同組成的一個(gè)十字排列道集(圖2)，兩條線的交點(diǎn)可視為十字交叉排列壓制面波處理中的一個(gè)面元。

在十字排列道集中，由于抽取的炮集之間的能量有差異，為了保證后續(xù)工作的展開(kāi)，首先，要對(duì)十字排列道集內(nèi)各道之間的能量進(jìn)行均衡處理[15]；然后，分別沿時(shí)間、縱向測(cè)線和橫向測(cè)線3個(gè)方向?qū)Φ卣鹩涗涍M(jìn)行三維傅里葉變換，變換后面波噪聲在三維頻率－波數(shù)域上呈現(xiàn)出圓錐狀(圖3)，因此，十字排列去噪方法也叫三維錐形濾波法。

圖2 十字排列道集

圖3 面波空間分布

地震記錄的褶積關(guān)系在頻率–波數(shù)域中變?yōu)椋?/p>

加入高通傾角濾波器，得到高通傾角濾波方程[17]：

高通濾波方程在頻率–波數(shù)域中是一個(gè)圓錐體。根據(jù)面波和有效波的視速度不同，體現(xiàn)在三維空間中就是圓錐的圓心角不同，適當(dāng)調(diào)整圓錐的圓心角，可以將面波與聲波干擾一同壓制并保留有效波[18]。在此基礎(chǔ)上，進(jìn)行三維反傅里葉變換，完成十字排列去噪處理。

對(duì)于勘探區(qū)障礙物少，變觀少的區(qū)域，推薦使用十字交叉排列法去面波；而對(duì)于變觀多的工區(qū)，常用的方法是采用三維數(shù)據(jù)規(guī)則化，即將不在十字排列上的炮點(diǎn)和檢波點(diǎn)進(jìn)行移位后再去面波，但是該方法會(huì)產(chǎn)生誤差。圖4為淮北礦區(qū)全數(shù)字高密度地震資料面波衰減前后單炮對(duì)比，可見(jiàn)自適應(yīng)去噪方法壓制面波效果良好，經(jīng)過(guò)處理后，單炮記錄上的掃帚狀面波頻散消失，有效波突出且同相軸連續(xù)性好，深部信號(hào)也得以保留，對(duì)有效反射波損傷小。

圖4 面波衰減前后單炮對(duì)比

3 異常振幅噪聲的分頻自適應(yīng)壓制

由于工區(qū)內(nèi)村莊工廠多，公路密布，由人為或機(jī)械振動(dòng)產(chǎn)生的干擾波在地震剖面上呈現(xiàn)為隨機(jī)分布的異常振幅，大部分異常振幅噪聲在空間和時(shí)間上具有隨機(jī)性特征，其能量強(qiáng)，頻帶寬，異常振幅能量是有效波能量的10倍以上。由于激發(fā)條件不一，頻率分布廣泛，且與有效波在空間和時(shí)間上重疊，具有非地表一致性，為了壓制異常振幅干擾并盡可能保護(hù)有效波不受損傷，李雪英等[19]提出了連續(xù)硬閾值函數(shù)與自適應(yīng)閾值相結(jié)合的地震資料高頻噪聲壓制方法來(lái)提高信噪比的保真度；周懷來(lái)在一種改進(jìn)的軟閾值法的基礎(chǔ)上，提出了一種改進(jìn)的閾值函數(shù)和自適應(yīng)閾值選取相結(jié)合的去噪方法[20]，都取得了較好的效果。

計(jì)算中使用樣點(diǎn)數(shù)據(jù)包絡(luò)的橫向加權(quán)中值作為干擾識(shí)別參數(shù)，將橫向加權(quán)中值[18]定義為：

圖5為分頻去除異常振幅干擾前后對(duì)比圖，圖5a為去噪前的剖面，局部有貫穿整個(gè)時(shí)間道的高能量異常振幅，也有連續(xù)多道的傾斜異常振幅，噪聲貫穿了有效波同相軸并壓制了深部的地震信號(hào)。經(jīng)過(guò)處理后(圖5b)的地震記錄上反射波同相軸連續(xù)且清晰，深部無(wú)噪聲干擾，去噪效果良好。

圖5 分頻去噪前后對(duì)比

4 多次波處理

當(dāng)?shù)叵掠卸鄠€(gè)反射系數(shù)較大的反射界面時(shí)，地震波在反射界面之間多次反射，形成多次反射波。多次波是淮北礦區(qū)地震資料處理中常見(jiàn)的干擾波，多次波的出現(xiàn)給處理和解釋帶來(lái)很多麻煩，多次波通常與一次反射波相互疊加，不利于有效波的識(shí)別與追蹤，對(duì)構(gòu)造解釋影響很大。更為突出的矛盾是在巖性解釋，如做AVO、亮點(diǎn)、頻譜比等特殊方法處理時(shí)，由于多次波存在造成解釋可信度降低?；幢钡V區(qū)地震資料中，由新生界底產(chǎn)生的水平多次波較嚴(yán)重地影響了速度分析精度，干擾了露頭、傾斜地層的成像，因此，壓制多次波是非常必要的。

4.1 多次波壓制原理

常規(guī)壓制多次波的方法主要有3種：① 利用一次波與多次波之間正常時(shí)差的差別。由于多次波是由淺層的全程或?qū)娱g的多次反射所形成，其傳播速度為淺層的速度，低于深層一次波的速度。因此，可以利用速度差異進(jìn)行多次波壓制，比如在速度分析時(shí)避開(kāi)低速能量團(tuán)或使用線性去噪方法壓制低速多次波。② 減去法。通過(guò)在原始資料中減去由準(zhǔn)確模型所得多次波，從而取得良好壓制效果。減去法要求多次波具有準(zhǔn)確的時(shí)間、振幅、相位以及周期性等，因此，實(shí)際應(yīng)用起來(lái)難以取得良好效果。一般而言，海水全程多次波、鳴震等多次波用這一類壓制方法較易見(jiàn)效。③ 利用多次波的重復(fù)性和統(tǒng)計(jì)特性采用預(yù)測(cè)反褶積等方法壓制短周期多次波。預(yù)測(cè)反褶積方法要求多次波具有穩(wěn)定的周期性，如果預(yù)測(cè)步長(zhǎng)選擇不合適，很容易損傷一次反射，預(yù)測(cè)反褶積適于壓制短周期如海水鳴震等類型的多次波，該類方法更多地在疊后應(yīng)用。

淮北礦區(qū)地震資料多次波與一次波之間速度差異較大，應(yīng)用高精度Radon變換壓制多次波技術(shù)可有效去除多次波干擾。該方法通過(guò)拋物Radon變換把動(dòng)校正后的CMP道集變換到–域。–域的地震數(shù)據(jù)經(jīng)過(guò)Radon正變換到到–域以后，聲波、直達(dá)波、淺層折射波等線性干擾波變成了一個(gè)能量團(tuán)，反射波變成了橢圓。對(duì)于經(jīng)過(guò)一次波動(dòng)校正后的CMP道集來(lái)說(shuō)，Radon正變換后，一次波能量分布在零道附近，而多次波能量呈彎曲條帶狀分布在遠(yuǎn)離零道的地方。這樣可以在–域中把多次波從地震資料中分離出來(lái)，將分離出的多次波進(jìn)行Radon反變換后再將它從原始地震資料中減去，能夠達(dá)到消除多次波目的。

4.2 應(yīng)用實(shí)例

由淮北工區(qū)多次波濾除前后疊加剖面(圖6)進(jìn)行對(duì)比可以發(fā)現(xiàn)：處理前的疊加剖面(圖6a)紅色箭頭處有一條不清晰的水平同相軸，是上部強(qiáng)反射界面的多次波。在處理后的剖面上(圖6b)多次波消失，上下地層的角度不整合接觸關(guān)系更加清晰，下部?jī)A斜地層的同相軸也變得連續(xù)。

5 結(jié)論

a. 根據(jù)面波在頻率–波數(shù)域中的分布特點(diǎn)，首先通過(guò)十字交叉排列分析，確定錐形濾波器的濾波參數(shù)，再在全區(qū)或某一區(qū)域應(yīng)用自適應(yīng)面波衰減技術(shù)，以壓制面波。

b. 對(duì)于異常振幅進(jìn)行多道識(shí)別、計(jì)算振幅的包絡(luò)并設(shè)定閾值、分頻單道處理，對(duì)振幅高于閾值的異常振幅進(jìn)行單道壓制。

c. 通過(guò)Radon變換在–域中將有效波與多次波進(jìn)行分離，根據(jù)有效波與多次波能量分布的差異來(lái)壓制多次波。

圖6 多次波濾除前后疊加剖面

d. 疊前去噪處理作為煤礦采區(qū)高密度處理必不可少的基本流程，需要根據(jù)采集原始資料特點(diǎn)，選擇合理方法與處理模塊組合，以提高資料質(zhì)量，在淮北工區(qū)得到了成功應(yīng)用，并取得了較好成果。

[1] 張成玉. 基于多域的去噪技術(shù)方法研究[D]. 北京：中國(guó)石油大學(xué)，2009. ZHANG Chengyu. Technical research based on multi-domain noise elimination[D]. Beijing：China University of Petroleum，2009.

[2] 王琦. 全數(shù)字高密度三維地震勘探技術(shù)在淮北礦區(qū)的應(yīng)用[J]. 煤田地質(zhì)與勘探，2018，46(增刊1)：41–45. WANG Qi. Application of all digital high-density 3D seismic ex-ploration technology in Huaibei Mining Area[J]. Coal Geology & Exploration，2018，46(Sup.1)：41–45.

[3] 馮曉輝，陳怡斌. 組合去噪在東嶺地區(qū)資料處理中的應(yīng)用[J]. 天然氣技術(shù)與經(jīng)濟(jì)，2012，6(1)：33–36. FENG Xiaohui，CHEN Yibin. Application of combined denoising to seismic-data processing in Dongling area[J]. Natural Gas Technology and Economy，2012，6(1)：33–36.

[4] 趙立明. 淮北礦區(qū)高密度三維地震勘探巖性解釋技術(shù)研究[D]. 徐州：中國(guó)礦業(yè)大學(xué)，2015. ZHAO Liming. Study on lithology Interpretation technology of high-density 3D seismic data in Huaibei Mining Area[D]. Xuzhou：China University of Mining and Technology，2015.

[5] 王金龍，胡治權(quán). 三維錐形濾波方法研究及應(yīng)用[J]. 石油地球物理勘探，2012，47(5)：705–711. WANG Jinlong，HU Zhiquan. Three dimensional cone filtering[J]. Oil Geophysical Prospecting，2012，47(5)：705–711.

[6] 李文花. 三維十字交叉排列去噪技術(shù)的應(yīng)用實(shí)踐[J]. 工程地球物理學(xué)報(bào)，2020，17(4)：414–420. LI Wenhua. Application of 3D cross array denoising technique[J]. Chinese Journal of Engineering Geophysics，2020，17(4)：414–420.

[7] 金丹. 煤炭全數(shù)字高密度三維地震勘探關(guān)鍵技術(shù)研究[D]. 北京：煤炭科學(xué)研究總院，2015. JIN Dan. Research on key technology of the digital high-density 3D seismic exploration in coal mine[D]. Beijing：China Coal Research Institute，2015.

[8] 陳美年. 地震資料隨機(jī)干擾壓制技術(shù)研究[D]. 青島：中國(guó)海洋大學(xué)，2010. CHEN Meinian. Research on random interference suppression tec-hnology of seismic data[D]. Qingdao：Ocean University of China，2010.

[9] 劉洋，劉財(cái)，王典，等. 時(shí)變中值濾波技術(shù)在地震隨機(jī)噪聲衰減中的應(yīng)用[J]. 石油地球物理勘探，2008，43(3)：327–332.LIU Yang，LIU Cai，WANG Dian，et al. Application of time-variant median filtering technique to attenuation of seismic random noises[J]. Oil Geophysical Prospecting，2008，43(3)：327–332.

[10] 刁塑，張恒琪，王冬年. 基于曲波變換在地震數(shù)據(jù)隨機(jī)噪聲壓制的應(yīng)用[J]. 能源研究與管理，2020(1)：87–90. DIAO Su，ZHANG Hengqi，WANG Dongnian. Application of curvelet transform in suppressing random noise of seismic data[J]. Energy Research and Management，2020(1)：87–90.

[11] 陳志德，鄭錫娟，關(guān)昕，等. 長(zhǎng)垣油田地震異常振幅噪聲分頻壓制技術(shù)[J]. 大慶石油地質(zhì)與開(kāi)發(fā)，2010，29(5)：164–167. CHEN Zhide，ZHENG Xijuan，GUAN Xin，et al. Frequency division suppression technology of seismic abnormal amplitude noise in Changyuan oil field[J]. Daqing Petroleum Geology and Development，2010，29(5)：164–167.

[12] 熊登，趙偉，張劍鋒. 混合域高分辨率拋物Radon變換及在衰減多次波中的應(yīng)用[J]. 地球物理學(xué)報(bào)，2009，52(4)：1068–1077. XIONG Deng，ZHAO Wei，ZHANG Jianfeng. Hybrid-demain high resolution parabolic transform and its application to demultiple[J]. Chinese Journal of Geophysics，2009，52(4)：1068–1077.

[13] 寧宏曉，唐東磊，皮紅梅，等. 國(guó)內(nèi)陸上“兩寬一高”地震勘探技術(shù)及發(fā)展[J]. 石油物探，2019，58(5)：645–653. NING Hongxiao，TANG Donglei，PI Hongmei，et al. The technology and development of “WBH” seismic exploration in land，China[J]. Geophysical Prospecting for Petroleum，2019，58(5)：645–653.

[14] 馬在田. 地震偏移剖面的假頻、頻散和橫向分辨力[J]. 石油地球物理勘探，1982，4(4)：16–29. MA Zaitian. False frequency，dispersion and lateral resolution of seismic migration profile[J]. Petroleum Geophysical Exploration，1982，4(4)：16–29.

[15] FRANKLINVILLE A，MEEK R. A procedure for optimally removing localized coherent noise[J]. Geophysics，1995，60(1)：191–203.

[16] 渥·伊爾馬滋. 地震資料分析[M]．劉懷山，王克斌，童思友，等，譯. 北京：石油工業(yè)出版社，2006：120． ILMAZ W. Seismic data analysis[M]. LIU Huaishan，WANG Kebin，TONG Siyou，et al. Translate. Beijing：Petroleum Industry Press，2006：120.

[17] 薛超. 十字排列去噪方法研究與應(yīng)用[J]. 內(nèi)蒙古石油化工，2019，45(7)：11–16. XUE Chao.The study of Cross-spreading de-noise method and its application[J]. Inner Mongolia Petrochemical Industry，2019，45(7)：11–16.

[18] HALE D，CLAERBOUTJ F. Butterworth dip filters[J]. Geophysics，2012，48(8)：1033–1038.

[19] 李雪英，張晶，孔祥琦，等. 基于離散小波變換的地震資料自適應(yīng)高頻噪聲壓制[J]. 物探與化探，2013，37(1)：165–170. LI Xueying，ZHANG Jing，KONG Xiangqi，et al. High frequency seismic noise adaptive suppression based on discrete wavelet transform[J]. Geophysical &Geochemical Exploration，2013，37(1)：165–170.

[20] 周懷來(lái)，李錄明，羅省賢，等. 基于小波變換的改進(jìn)閾值函數(shù)自適應(yīng)去噪方法[J]. 物探化探計(jì)算技術(shù)，2006，28(2)：173–177. ZHOU Huailai，LI Luming，LUO Shengxian，et al. A method of improved threshold function denoising based on wavelet analysis[J]. Computing techniques for geophysical and geochemical exploration，2006，28(2)：173–177.

[21] 蔡希玲. 聲波和強(qiáng)能量干擾的分頻自適應(yīng)檢測(cè)與壓制方法[J]. 石油地球物理勘探，1999，34(4)：373–380. CAI Xiling. An effective method to suppress acoustic wave and high energy noise frequency-divisionally and adaptively[J]. Oil Geophysical Prospecting，1999，34(4)：373–380.

[22] 秦婕，李輝峰，王宏偉，等. 疊前去噪技術(shù)在鄂爾多斯黃土塬區(qū)地震資料的應(yīng)用[J]. 物探化探計(jì)算技術(shù)，2015，37(5)：644–650.QIN Jie，LI Huifeng，WANG Hongwei，et al. The application ofpre-stack noise attenuation technology for seismic data of Soil Yuan area in Ordos basin[J]. Computing Techniques for Geophysical and Geochemical Exploration，2015，37(5)：644–650.

Application of pre-stack denoising technique in full digital high density 3D seismic technique in coal mining districts

YU Jie

(School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, China)

High density 3D seismic technique uses digital detectors to receive data. The high sensitivity of digital detectors makes weaker background noises and the effective signals are recorded by adopting broadband acceptance. Therefore, in the full digital high density seismic exploration, the data collected by digital detectors have higher noise level and lower SNR. Taking the full digital high density 3D seismic data of Huaibei mining area as an example, the paper explains the characteristics of the three main noises in the high density 3D seismic data, including surface wave, abnormal amplitude and multiple wave, and combines the characteristics of the high density 3D seismic data including wide azimuth, high density and multi-offset to makes cross-alignment analysis of surface wave in the direction of in-line and off-line, and designs conical filter according to the characteristics of surface wave to filter it, sets a threshold to identify abnormal amplitude and performs single channel suppression after calculating the envelope of multiple seismic channels, and applies the Radon transform to eliminate multiple waves from the CMP channel set into the domain. After processing the original data, the effective signals are preserved and the SNR is greatly improved.

high density 3D seismic exploration; interference wave; cross array; automatic identification; Radon transformation

請(qǐng)聽(tīng)作者語(yǔ)音介紹創(chuàng)新技術(shù)成果等信息，歡迎與作者進(jìn)行交流

P631

A

10.3969/j.issn.1001-1986.2020.06.007

1001-1986(2020)06-0048-07

2020-10-25；

2020-11-18

國(guó)家重點(diǎn)研發(fā)計(jì)劃課題(2018YFC0807803)；中國(guó)礦業(yè)大學(xué)雙一流建設(shè)自主創(chuàng)新專項(xiàng)項(xiàng)目(2018ZZCX05)

National Key R&D Program of China(2018YFC0807803)；Independent Innovation Project of Double First Class Construction of China University of Mining and Technology(2018ZZCX05)

于杰，1997年生，男，山東日照人，碩士研究生，研究方向?yàn)榈厍蛭锢韺W(xué). E-mail：juxiandiyi@163.com

于杰. 疊前去噪技術(shù)在煤礦采區(qū)全數(shù)字高密度三維地震中的應(yīng)用[J]. 煤田地質(zhì)與勘探，2020，48(6)：48–54.

Application of pre-stack denoising technique in full digital high density 3D seismic technique in coal mining districts[J]. Coal Geology & Exploration，2020，48(6)：48–54.

(責(zé)任編輯 聶愛(ài)蘭)