楊光明，金學良，張憲旭，智 敏，單 蕊

寬頻寬方位處理技術在淮北礦區(qū)全數(shù)字高密度地震勘探中的應用

楊光明1，金學良2，張憲旭1，智 敏1，單 蕊1

(1. 中煤科工集團西安研究院有限公司，陜西 西安 710077；2.淮北礦業(yè)(集團)有限責任公司，安徽 淮北 235000)

與常規(guī)三維地震勘探相比，全數(shù)字高密度地震勘探采用高橫縱比、寬方位觀測系統(tǒng)，使用單點數(shù)字檢波器接收。寬方位地震勘探有利于高陡構造和復雜斷塊成像，但存在各向異性問題，而單點數(shù)字檢波器地震信號頻帶寬、保幅性好、噪聲強。為了充分發(fā)揮全數(shù)字高密度地震資料優(yōu)勢，克服其缺點，必須將寬頻、寬方位處理技術運用到煤炭高密度三維地震數(shù)據(jù)處理中，以提升數(shù)據(jù)處理效果。以淮北礦區(qū)全數(shù)字高密度地震資料為基礎，針對淮北礦區(qū)地質構造復雜特點，以疊前保幅去噪、振幅補償、OVT處理技術以及全方位角偏移成像技術為重點，開展了煤炭全數(shù)字高密度地震資料的寬頻、寬方位處理技術研究。結果表明：保幅去噪在幾乎不損傷有效信號的前提下實現(xiàn)了對噪聲的有效壓制，振幅補償恢復了地震信號高低頻能量的損失，寬方位處理不僅消除了各向異性影響，改善了成像效果，還獲得了豐富的疊前數(shù)據(jù)。寬頻寬方位處理技術是全數(shù)字高密度地震資料處理的必要手段，其處理成果比常規(guī)處理成果頻帶更寬，對復雜構造成像更好，分辨率更高，能夠實現(xiàn)煤田復雜地質條件地震資料精細成像。

數(shù)字檢波器；高密度三維地震；寬頻；寬方位；分頻去噪

淮北煤田煤層傾角大，斷層發(fā)育，底板灰?guī)r含水層嚴重影響煤礦的安全生產(chǎn)，對地震勘探提出了更高的地質任務要求。同時，淮北礦區(qū)地表起伏小，潛水面僅為4～5 m，煤層與圍巖波阻抗差異大，有利于開展三維地震勘探。近20年來淮北礦區(qū)開展了超過500 km2的三維地震勘探，隨著地震勘探儀器的快速發(fā)展，全數(shù)字高密度地震勘探已經(jīng)廣泛應用于油氣勘探，在煤炭系統(tǒng)，也已開展了很多的應用，取得了一定的效果[1]。自2014年以來，淮北礦區(qū)陸續(xù)開展了超過120 km2的全數(shù)字高密度地震勘探[2]，獲得了豐富的煤炭全數(shù)字高密度地震勘探數(shù)據(jù)與驗證資料，為研究煤礦采區(qū)全數(shù)字高密度地震勘探技術奠定了基礎。

全數(shù)字高密度地震勘探與常規(guī)地震勘探主要區(qū)別于采用了數(shù)字檢波器接收，實現(xiàn)了小面元、高覆蓋次數(shù)與寬方位采集[3]。相對于模擬檢波器，數(shù)字檢波器采用MEMS加速度傳感器技術，具有動態(tài)范圍大、信號畸變小、沒有相位延遲的優(yōu)勢[4-5]。但是由于數(shù)字檢波器采用單點接收，對高低頻信號均有良好響應，因此，其對噪聲也更為敏感，導致其單炮記錄較常規(guī)模擬檢波器組合接收單炮記錄更“臟”(噪聲)，對后續(xù)去噪處理提出了更高要求[6]。隨著高密度勘探橫向覆蓋次數(shù)的顯著增加，接收排列橫縱比超過0.5，實現(xiàn)了寬方位采集[7]，有利于改善地下目標的照明度，提高復雜構造的成像效果[8-9]，更容易跨越地表障礙物和地下陰影帶[10]。高密度勘探采用小面元有利于提高縱橫向分辨率，同時高覆蓋次數(shù)也改善了地震剖面的信噪比。

在煤炭全數(shù)字高密度地震資料處理方面，對于保護低頻并未得到重視和充分的研究，往往忽略了數(shù)字檢波器接收到的豐富信息。而對于寬方位采集的高密度地震資料，依然采用了常規(guī)窄方位處理技術，未能發(fā)揮寬方位數(shù)據(jù)的優(yōu)勢，難以體現(xiàn)出全數(shù)字高密度地震勘探的優(yōu)點。因此，研究寬頻、寬方位處理技術在煤炭全數(shù)字高密度地震資料處理中的應用很有意義?；诖耍P者以淮北礦區(qū)全數(shù)字高密度地震資料為基礎，對煤炭全數(shù)字高密度地震資料寬頻、寬方位處理技術進行了探討，總結了疊前保幅去噪、振幅補償、OVT處理技術以及全方位角偏移成像技術的適用性，以期形成煤炭高密度三維地震數(shù)據(jù)處理相應的寬頻、寬方位處理技術流程。

1 淮北煤田全數(shù)字高密度地震勘探采集參數(shù)及資料特點

1.1 采集參數(shù)

由于地下地質構造復雜，淮北煤田大部分采區(qū)開展了常規(guī)三維地震勘探，部分重點采區(qū)在常規(guī)勘探的基礎上又開展了全數(shù)字高密度地震勘探?；幢钡V區(qū)常規(guī)地震勘探采集觀測系統(tǒng)大多采用8線8炮(圖1a)、8線10炮的束狀觀測系統(tǒng)，面元大小為10 m×10 m，是典型的窄方位觀測系統(tǒng)。而近5年來開展的全數(shù)字高密度地震勘探，觀測系統(tǒng)由16線10炮160道接收逐步發(fā)展為24線4炮96～160道接收(圖1b)，面元大小為5 m×5 m，甚至在某些勘探區(qū)采用了40線2炮80道接收200次覆蓋的全三維觀測系統(tǒng)(圖1c)。由圖1可知，隨著橫向采樣密度的不斷增加，橫縱比由0.4增大為0.8～1.0。同時高密度勘探觀測系統(tǒng)炮檢距近、中、遠分布均勻，能夠提高速度分析精度，改善偏移成像質量，減小采集腳印影響，小面元不僅能提高縱橫向分辨率，也避免了空間假頻的發(fā)生[11]。

圖1 常規(guī)及高密度觀測系統(tǒng)對比

1.2 地震資料特點

數(shù)字檢波器單點接收單炮較常規(guī)模擬檢波器組合接收單炮干擾更強，對低視速度面波壓制只能通過室內衰減。圖2為某勘探區(qū)相近位置采集的常規(guī)模擬檢波器接收單炮與全數(shù)字高密度勘探單炮記錄，高密度勘探排列長度遠大于常規(guī)勘探，可以明顯看出數(shù)字檢波器接收單炮面波遠強于模擬檢波器單炮，紅框所示目的層反射波也顯示出了低頻復合波的特點，面波與目的層反射波在近道大范圍重疊。圖3b是對圖3a所示全數(shù)字高密度地震典型單炮紅框范圍內的面波進行了頻譜分析，面波視速度低，主頻低于20 Hz，而煤層有效反射波主頻在25~80 Hz，與面波差異較大。除了面波干擾，全數(shù)字高密度地震勘探普遍存在異常振幅干擾，這些干擾能量強，頻帶寬(圖4)，嚴重影響后續(xù)的振幅補償、反褶積、偏移成像等處理。因此，疊前保幅保頻去噪是淮北礦區(qū)全數(shù)字高密度資料寬頻處理的重點。

圖2 模擬檢波器與數(shù)字檢波器單炮記錄

圖3 面波頻譜分析

圖4 淮北礦區(qū)全數(shù)字高密度地震純波疊加剖面

2 寬頻處理關鍵技術

地震數(shù)據(jù)低頻信號在數(shù)據(jù)處理和解釋中有重要意義，低頻信號能夠提高分辨率與成像精度、改善反演質量。同時低頻成分穿透力強，比高頻穩(wěn)定、傳播距離遠，在巖性識別，薄煤層成像和全波形反演中起到至關重要的作用，能夠有效提升煤層下部奧灰等地層成像效果。寬頻處理工作的重點是拓寬頻帶，尤其是保護低頻，寬頻處理更應該重視提高倍頻程，通常要求倍頻程能夠達到3～5，而增強低頻成分是確保頻帶寬度和提高倍頻程數(shù)的最有效做法[12]?；幢钡V區(qū)全數(shù)字高密度勘探使用數(shù)字檢波器保證了低頻信號的高質量接收，因此，必須通過寬頻處理盡可能挖掘地震資料的解釋潛力，提供更豐富可靠的地震解釋成果，保證煤礦的安全生產(chǎn)。寬頻處理的重點是采用對低頻及高頻有效信號損傷小的保幅去噪手段，并綜合應用振幅補償、反Q濾波、反褶積等方法補償高低頻信號損失、拓寬頻帶。

2.1 疊前保幅綜合去噪

由圖3及圖4分析可知，淮北礦區(qū)全數(shù)字高密度地震數(shù)據(jù)處理疊前去噪的重點是壓制面波及異常振幅干擾，該區(qū)地震資料面波普遍具有低頻特征，而異常振幅干擾頻帶寬，能量強。在常規(guī)處理中常使用帶通濾波壓制面波及高頻的異常振幅干擾，必然損傷了地震資料的高低頻有效信號，因此，在寬頻處理中是不可取的。

由于面波與有效反射信號之間在分布范圍、頻率及視速度等方面有較大差異，因此，對于淮北礦區(qū)全數(shù)字高密度地震資料，可以采用自適應面波衰減和十字交叉排列去噪技術壓制面波。這2種方法都能很好壓制面波，且對有效信號損傷小。由于寬方位觀測系統(tǒng)橫向采樣更充足，因此，更有利于十字排列域壓制面波。對于障礙物少的勘探區(qū)，推薦使用十字交叉排列法去面波。

和有效反射波相比異常振幅干擾能量通常遠大于有效反射能量，大部分異常振幅噪聲在空間和時間上具有隨機性特征，其能量強，頻帶寬，且與有效波在空間和時間上重疊，難以區(qū)分?；诋惓Ｕ穹肼曇陨咸攸c，在選取去噪方法時，采用“多道識別，分頻單道處理”的思想對疊前的異常強能量進行壓制，這種分頻處理方法具有較好的保真性。

2.2 振幅補償

能夠引起地震波衰減的不確定因素有很多[13]，但主要是球面擴散及地層吸收，此外就是由于激發(fā)接收條件的差異造成的炮間、道間的振幅不均衡?；幢钡V區(qū)全數(shù)字高密度地震資料振幅補償重點是球面擴散補償、地表一致性振幅補償及反Q濾波。

球面擴散和地層吸收是引起地震波能量衰減的最主要因素，會導致地震波能量從淺到深，從近到遠不斷減弱。球面擴散補償以速度和時間為參數(shù)計算補償參數(shù)，對地震道各采樣點的振幅進行加權補償。

地表一致性振幅補償在給定的時窗內計算全部地震道的能量，采用統(tǒng)計方法求出各炮集、各接收道集、各炮檢距道集上的補償系數(shù)，然后將不同道集補償系數(shù)同時應用于各數(shù)據(jù)道完成補償處理。地表一致性振幅補償在多域進行多道統(tǒng)計，保幅效果好，能較好解決空間上由于激發(fā)和接收條件的差異引起的能量差異，有效補償由于炮檢距變化所引起的能量衰減?；幢钡V區(qū)全數(shù)字地震資料存在較強異常振幅及面波干擾，會嚴重影響振幅統(tǒng)計結果，因此，地表一致性振幅補償應在去噪及球面擴散補償后進行。

反Q濾波是補償大地吸收作用的一種有效方法，可以補償?shù)卣鸩ǖ恼穹p和頻率損失，還可改善地震記錄的相位特性，從而改善同相軸的連續(xù)性，增強弱反射波的能量，進而提高地震資料的信噪比和分辨率。反Q濾波方法可分為級數(shù)展開作近似高頻補償?shù)姆碤濾波、基于波場延拓的反Q濾波和其他類型的反Q濾波等三大類[14]。目前，Q補償最大的難點在于求取準確的Q值，雖然很多處理軟件中提供了各種Q值估算工具，但是Q補償依然要進行大量的測試，才有可能取得良好的補償效果。

2.3 拓頻處理

在振幅補償后，需要在疊前應用反褶積技術壓縮子波，在疊后使用譜白化、藍色濾波等技術進一步拓寬頻帶，拓頻往往會造成高頻噪聲能量增強，可以采用俞氏濾波進行拓頻后的改善處理[15]。

反褶積通常應用于疊前資料，也可廣泛用于疊后資料。反褶積不僅能夠壓縮子波，它也能從記錄上消除大部分的短周期多次波，通過反褶積處理，能夠提高縱向分辨率。目前應用比較廣泛的有脈沖反褶積、地表一致性反褶積和預測反褶積。地表一致性反褶積能夠消除由于近地表條件的變化對地震子波波形的影響。脈沖反褶積應用于野外資料通常難以取得預期效果[16]，在大部分情況下尖脈沖反褶積提起來的主要是高頻的噪聲能量，因此，在常規(guī)處理中，預測反褶積應用更為廣泛。

淮北礦區(qū)全數(shù)字高密度地震勘探的目標不僅包含高主頻的新生界各地層、中高主頻的煤層反射波，也要對低頻為主的煤層下伏太灰、奧灰頂成像，因此在淮北礦區(qū)全數(shù)字高密度地震資料處理中，往往會采用不同預測步長的時變預測反褶積進行疊前反褶積處理。圖5是淮北礦區(qū)全數(shù)字高密度地震資料在應用了地表一致性預測反褶積前后單炮(圖5a)、自相關(圖5b)以及頻譜對比圖(圖5c)，地表一致性反褶積消除了子波的橫向差異，預測反褶積壓制了交混回響并提高了分辨率，反褶積后頻譜(藍色線)明顯拓寬，反射波波組特征明顯改善。

圖5 地表一致性預測反褶積效果對比

3 寬方位處理技術

如前所述，寬方位采集相對窄方位有很多優(yōu)勢，但是如果采用常規(guī)窄方位的地震資料處理技術，勢必浪費了寬方位采集帶來的很多有效信息。寬方位地震資料的速度、旅行時差會隨著觀測方位角的變化而有一定差異，存在各向異性問題[17]。2002年A. Cordsen認為由于地下介質的非均質性和方位各向異性等原因，導致地震波在地下傳播時在不同的方向上速度不同，僅采用過去的單一速度難以解決不同方位的偏移成像問題，因而要考慮基于方位的速度分析[18]。國內石油地震勘探最早使用分方位處理技術進行寬方位資料的處理，取得了較好的效果[10,19]。分方位處理首先按照一定的角度范圍將地震數(shù)據(jù)劃分為分方位道集，扇區(qū)劃分要保證各個方位上的覆蓋次數(shù)均勻分布，速度分析工作也需要在分方位道集上分別進行，以獲得各個方位的最佳速度。偏移后的分方位道集可繼續(xù)做資料解釋，如通過方位切片和地震屬性分析等來判斷地層變化或者裂縫的發(fā)育情況[20]。分方位技術減小了速度的方位各向異性對處理效果的影響，但各自偏移后的分方位數(shù)據(jù)體無法重構為一個統(tǒng)一的成果，只能在各個分方位數(shù)據(jù)體上進行綜合性解釋，其本質還是窄方位處理。近年來，已經(jīng)發(fā)展出了OVT、全方位角偏移成像等寬方位處理技術，這些技術能更好的對寬方位地震數(shù)據(jù)進行處理。

3.1 OVT處理技術

段文勝等[21]在2013年首次將炮檢距向量片(Offset Vector Tiles，簡稱OVT)技術應用到國內塔里木油田，提高了碳酸鹽巖縫洞體成像和裂縫預測精度。OVT道集是對十字排列道集按照炮線距和檢波線距等距離劃分得到的子道集，十字排列道集是以相同炮線和檢波線抽取的道集，因此，其個數(shù)與炮線和檢波線交點的數(shù)目一致；對每個十字排列道集矩形按照炮線距和檢波線距等距離劃分，得到的小矩形就是一個OVT，其大小由炮線距和檢波線距決定，每個十字排列中OVT的個數(shù)等于覆蓋次數(shù)。將每個十字排列道集中相同位置的OVT抽取出來，就形成了一個覆蓋全區(qū)的單次覆蓋OVT道集，這些單次覆蓋的OVT道集具有相近的炮檢距和方位角。

在獲得OVT道集之后，可以在OVT域實現(xiàn)插值、疊前偏移等處理，相對于常規(guī)共炮點域、共炮檢距域有很多優(yōu)勢。常規(guī)窄方位地震資料由于各向異性影響小，在共炮檢距域進行插值問題不大，但是寬方位采集數(shù)據(jù)可能會存在嚴重各向異性問題，在處理中必須考慮各向異性問題，而OVT道集具有相近方位角，消除了各向異性影響，因此，OVT域的插值更適應寬方位的地震資料規(guī)則化處理。應用OVT道集進行疊前時間/深度偏移，獲得的CRP道集保留了炮檢距和方位角信息，可用于各向異性分析、疊前反演和裂縫預測。OVT域偏移后的道集可抽選成“蝸?！钡兰M行各項異性分析和校正，能夠有效消除全數(shù)字高密度地震資料寬方位采集帶來的各向異性問題，改善成像效果[22]。

3.2 基于地下局部角度域的全方位角偏移技術

全方位角成像技術是基于地下局部角度域的成像方法，是針對寬方位資料處理的一種新思路，OVT道集的方位角是建立在地表觀測系統(tǒng)之上的，而全方位角度域輸出的方位角、傾角和反射角均為地下反射點的角度，精度更高。全方位角成像技術是將地面采集的地震數(shù)據(jù)映射到地下局部角度域系統(tǒng)中進行偏移成像，每個地下成像點處的射線都包含開角、開方位角、傾角和方位角坐標分量，這4個量形成了一個4D局部角度域系統(tǒng)，可以對成像點處的射線進行準確的描述[23]。通過全方位角偏移成像處理，可以產(chǎn)生傾角道集和反射角道集。傾角道集通過加權疊加能夠改善斷面等非連續(xù)目標的成像效果，而反射角道集可以用于層析速度反演，有利于深度域速度建模和疊前反演。

4 應用實例分析

青東礦位于淮北礦區(qū)北部，區(qū)內地下構造復雜，構造主體表現(xiàn)為一走向NW，局部略有轉折，向北傾斜的單斜。該區(qū)全數(shù)字高密度地震勘探觀測系統(tǒng)采用16線8炮160道接收64次覆蓋線束狀觀測系統(tǒng)，面元大小5 m×5 m，線距80 m。

4.1 處理流程

該礦區(qū)地表起伏較小，低降速帶橫向有變化，存在一定的靜校正問題。全數(shù)字高密度地震資料目的層段信噪比較低，面波和異常振幅是主要干擾波。由于激發(fā)因素、接收條件影響，炮/道間能量有一定差異，受煤層等上覆地層屏蔽影響，煤層下部奧灰頂、太灰等目的層反射波能量弱。根據(jù)地震資料特點，結合地質任務，確定寬頻、寬方位處理流程如圖6所示，寬頻處理的重點是疊前保幅去噪和拓頻處理，而寬方位處理采用了全方位角偏移成像技術。

圖6 處理流程

4.2 寬頻處理

寬頻處理重點是疊前保頻噪聲壓制、振幅補償及反褶積處理。圖7是疊前應用自適應面波壓制及分頻去噪技術前后純波疊加剖面對比，面波及高頻強能量異常噪聲得到顯著壓制，去噪過程未使用帶通濾波，有效保護了低頻信號。振幅補償采用球面擴散補償、地表一致性振幅補償?shù)燃夹g，消除了縱橫向能量差異。

圖7 疊前保幅去噪前后剖面對比

疊前反褶積處理重點對比了脈沖反褶積、單道預測反褶積、多道預測反褶積，處理效果如圖8所示。由圖8可知：反褶積前剖面存在明顯短周期多次，剖面波組特征差，頻率低；而脈沖反褶積剖面雖然主頻提高，但部分目的層同相軸連續(xù)性變差；多道預測反褶積在400 ms處目的層反射波依然存在明顯復合波特征，其分辨率較低，不利于后續(xù)解釋；而單道預測反褶積分辨率更高、同相軸連續(xù)性以及波組特征更好。

圖8 反褶積方法對比

為了兼顧淺部及深部目的層成像，對預測步長(8、12、16、20、24、28 ms)進行了測試(圖9)，淺部使用12 ms預測步長時600 ms反射波分辨率最高、連續(xù)性好，而深部使用20 ms步長時700 ms反射波成像效果好。

4.3 全方位角偏移

全方位角疊前深度偏移核心是速度模型的精確求取。如圖10a所示，首先在疊前時間偏移數(shù)據(jù)體上進行構造解釋，拾取9個主要反射層位建立地質模型，在層析反演時作為射線追蹤點。圖10b是構建的地質模型，該模型縱向層位密度合理，橫向變化與區(qū)域構造特征吻合，縱橫向測線解釋方案閉合，層位解釋滿足后續(xù)處理要求。

其次，以構造模型為約束條件，通過網(wǎng)格層析成像技術進行速度模型迭代，直到CRP道集拉平，剩余延遲為零。最后，結合鉆井數(shù)據(jù)，應用基于模型的層析成像技術，獲得較準確的各向異性深度域速度體。

圖9 預測步長測試

最終，選擇合理偏移參數(shù)應用全方位角疊前深度偏移，獲得了全方位角深度域偏移數(shù)據(jù)體。通過與原處理剖面對比(圖11)可見，資料頻帶由原來的40～78 Hz拓寬到25～85 Hz(圖11c)，尤其是20～40 Hz低頻成分得到保留和加強，采用針對寬方位資料的全方位角疊前深度偏移處理技術，改善了復雜斷塊資料的成像效果，斷層反映更為清晰，更有利于構造解釋，寬頻寬方位處理成果400 ms振幅切片(圖12)分辨率更高。

5 結論

a. 寬頻寬方位處理技術是淮北礦區(qū)全數(shù)字高密度地震資料處理的關鍵技術，通過寬頻、寬方位處理，能夠提高地震數(shù)據(jù)的縱橫向分辨率和復雜構造解釋可靠程度，改善煤層下伏太灰、奧灰頂成像效果，更有利于開展后續(xù)反演、屬性分析。

b. 寬頻處理的重點在于疊前疊后的保幅保頻處理，尤其是在疊前去噪、振幅補償以及反褶積處理環(huán)節(jié)。疊前去噪必須采用不損傷有效信號的去噪手段，對于低頻面波宜采用自適應面波壓制或十字排列去噪方法進行去除，對于高頻異常振幅干擾，可采用多道分頻分時去噪技術壓制。疊前地表一致性預測反褶積是寬頻處理重要的拓頻手段，應根據(jù)資料特點確定適當反褶積參數(shù)。為了保證不同深度目的層的分辨率和信噪比，可以采用時變步長預測反褶積。疊后還可靈活使用譜白化、藍色濾波、Q補償?shù)燃夹g，進一步拓寬頻帶。

c.寬方位采集數(shù)據(jù)包含豐富的地質信息，但是也帶來了方位各向異性問題，OVT和全方位角處理技術，可以較好地解決方位各向異性問題，能夠提高復雜構造成像效果和縱橫向分辨率，是煤炭全數(shù)字高密度地震資料寬方位處理的關鍵技術。

圖11 常規(guī)與寬頻寬方位處理剖面

圖12 常規(guī)與寬頻寬方位處理水平振幅切片(400 ms)

[1] 程建遠，聶愛蘭，張鵬. 煤炭物探技術的主要進展及發(fā)展趨勢[J]. 煤田地質與勘探，2016，44(6)：136–141. CHENG Jianyuan，NIE Ailan，ZHANG Peng. Outstanding progress and development trend of coal geophysics[J]. Coal Geology & Exploration，2016，44(6)：136–141.

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

[3] 趙鐠，武喜尊. 高密度采集技術在西部煤炭資源勘探中的應用[J]. 中國煤炭地質，2008，20(6)：11–14. ZHAO Pu，WU Xizun. High density acquisition technology and it’s application in Western China coal resource exploration[J]. Coal Geology of China，2008，20(6)：11–14.

[4] 張宏樂，曹金棟. 用于地球物理勘探的MEMS加速度傳感器[J]. 石油儀器，2002，16(4)：1–4.ZHANG Hongle，CAO Jindong. MEMS acceleration sensor for geophysical prospecting[J]. Petroleum Instruments，2002，16(2)：1–4.

[5] 王喜雙，董世泰，王梅生. 全數(shù)字地震勘探技術應用效果及展望[J]. 中國石油勘探，2007，12(6)：32–36. WANG Xishuang，DONG Shitai，WANG Meisheng. Application and prospect of full digital seismic exploration technology[J]. China Petroleum Exploration，2007，12(6)：32–36.

[6] 金丹，程建遠，張憲旭. 高密度全數(shù)字地震勘探技術在煤田中的應用及效果分析[J]. 煤炭技術，2015，34(8)：89–92.JIN Dan，CHENG Jianyuan，ZHANG Xianxu. Application and effect analysis of single-point high-density digital seismic exploration in coal seismic exploration[J]. Coal Technology，2015，34(8)：89–92.

[7] 張曉江. 寬、窄方位角三維地震勘探采集方法研究與應用[D]. 北京：中國石油大學，2007. ZHANG Xiaojiang. The research and application of wide and narrow azimuth 3D survey[D]. Beijing：China University of Petroleum，2007.

[8] 劉依謀，印興耀，張三元，等．寬方位地震勘探技術新進展[J]. 石油地球物理勘探，2014，49(3)：596–610. LIU Yimou，YIN Xingyao，ZHANG Sanyuan，et al. Recent advances in wide-azimuth seismic exploration[J]. Oil Geophysical Prospecting，2014，49(3)：596–610.

[9] 王華忠．“兩寬一高”油氣地震勘探中的關鍵問題分析[J]. 石油物探，2019，58(3)：313–324. WANG Huazhong. Key problem analysis in seismic exploration based on wide-azimuth，high-density，and broadband seismic data[J]. Geophysical Prospecting for Petroleum，2019，58(3)：313–324.

[10] 夏常亮，張紅軍，郝建波，等. K區(qū)塊寬方位角地震資料處理[J]. 石油地球物理勘探，2010，45(增刊1)：74–79. XIA Changliang，ZHANG Hongjun，HAO Jianbo，et al. Wide azimuth seismic data processing for K block[J]. Oil Geophysical Prospecting，2010，45(Sup.1)：74–79.

[11] 劉勝，程增慶，代鳳紅，等. 煤炭高密度三維地震勘探方法及應用[J]. 煤炭工程，2009(3)：67–68. LIU Sheng，CHENG Zengqing，DAI Fenghong，et al. Method and application of high density 3D seismic coal exploration[J]. Coal Engineering，2009(3)：67–68.

[12] 李慶忠. 走向精確勘探的道路[M]. 北京：石油工業(yè)出版社，1993. LI Qingzhong. The way to obtain a better resolution in seismic prospecting[M]. Beijing：Petroleum Industry Press，1993.

[13] 李振春，王清振. 地震波衰減機理及能量補償研究綜述[J]. 地球物理學進展，2007，22(4)：1147–1152. LI Zhenchun，WANG Qingzhen. A review of research on mechanism of seismic attenuation and energy compensation[J]. Progress in Geophysics，2007，22(4)：1147–1152.

[14] 張瑾. 地震波能量補償反Q濾波方法研究[D]. 長春：吉林大學，2013.ZHANG Jin. Study on the inverse Q filtering method for seismic wave energy compensation[D]. Changchun：Jilin University，2013.

[15] 蔡希玲. 俞氏子波在地震數(shù)據(jù)處理中的應用研究[J]. 石油地球物理勘探，2000，35(4)：497–507. CAI Xiling. Application of Yu wavelet to seismic data processing[J]. Oil Geophysical Prospecting，2000，35(4)：497–507.

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

[17] 寧宏曉，唐東磊，皮紅梅，等. 國內陸上“兩寬一高”地震勘探技術及發(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，2019，58(5)：645–653.

[18] CORDSEN A. Narrow-versus wide-azimuth land 3D seismic surveys[J]. Leading Edge，2002，21(8)：764–770.

[19] 張保慶，周輝，左黃金，等. 寬方位地震資料處理技術及應用效果[J]. 石油地球物理勘探，2011，46(3)：396–400. ZHANG Baoqing，ZHOU Hui，ZUO Huangjin，et al. Wide azimuth data processing techniques and their applications[J]. Oil Geophysical Prospecting，2011，46(3)：396–400.

[20] 張公社，馬國光，宋玉龍，等. 利用全三維縱波資料進行裂縫檢測[J]. 石油地球物理勘探，2004，39(1)：41–44. ZHANG Gongshe，MA Guoguang，SONG Yulong，et al. Fracture detection by using full 3-D P-wave data[J]. Oil Geophysical Prospecting，2004，39(1)：41–44.

[21] 段文勝，李飛，王彥春，等. 面向寬方位地震處理的炮檢距向量片技術[J]. 石油地球物理勘探，2013，48(2)：206–213. DUAN Wensheng，LI Fei，WANG Yanchun，et al. Offset vector tile for wide-azimuth seismic processing[J]. Oil Geophysical Prospecting，2013，48(2)：206–213.

[22] 馬淵明，楊雪霖，張麗梅，等. GeoEast-OVT寬方位處理技術[J]. 石油工業(yè)計算機應用，2018，26(2/3/4)：31–34. MA Yuanming，YANG Xuelin，ZHANG Limei，et al. GeoEast-OVT wide azimuth processing technology[J]. Computer Applications of Petroleum，2018，26(2/3/4)：31–34.

[23] 龔翟韶. 全方位角成像技術在躍滿三維中的應用研究[D]. 北京：中國石油大學(北京)，2017. GONG Zhaishao. The Applied research of Full-azimuth imaging in Yueman 3D[D]. Beijing：China University of Petroleum(Beijing)，2017.

Application of broadband and wide azimuth processing technology in full digital high density seismic exploration in Huaibei mining area

YANG Guangming1, JIN Xueliang2, ZHANG Xianxu1, ZHI Min1, SHAN Rui1

(1. Xi’an Research Institute Co. Ltd., China Coal Technology and Engineering Group Corp.,Xi’an 710077, China; 2. Huaibei Mining(Group)Co.,Ltd.,Huaibei 235000, China)

Compared with conventional 3D seismic exploration, full-digital high density seismic exploration adopts high horizontal and vertical ratio, wide azimuth observation system and single-point digital detectors .Wide azimuth seismic exploration is beneficial to high and steep structure and complex fault block imaging, but there are anisotropy problems, while single-point digital geophone seismic signal has wide frequency band, good amplitude preservation and strong noise. In order to give full play to the advantages of the full digital high-density seismic data and overcome its disadvantages, broadband and wide azimuth processing technology must be applied to high density 3D seismic data processing in coal mining districts to improve the data processing effect. Based on the full digital high density seismic data of Huaibei mining area, this paper studies the broadband and wide-azimuth processing technology of the full digital high-density seismic data of coal mining districts, focusing on pre-stack amplitude preservation denoising, amplitude compensation, OVT processing technology and omni-directional angular migration imaging technology according to the complex geological structure characteristics of Huaibei mining area. The actual digital high density seismic data processing in Huaibei mining area shows that the image denoising in almost no damage under the premise of effective signals to realize the effective suppression of noise, restore the seismic amplitude compensation signal of high frequency energy loss, wide azimuth processing not only eliminates the anisotropy effect to improve the imaging effect, but also obtains rich prestack data. Broadband azimuth processing technology is a necessary means for processing full digital high density seismic data. Its processing results are wider in frequency band than conventional processing results, better in imaging complex structures, higher in resolution, and able to achieve fine imaging of seismic data of complex geological conditions in coal mining districts.

digital detector; high density 3D seismic exploration; broadband; wide azimuth; frequency division denoising

請聽作者語音介紹創(chuàng)新技術成果等信息，歡迎與作者進行交流

P631

A

10.3969/j.issn.1001-1986.2020.06.008

1001-1986(2020)06-0055-09

2020-10-25；

2020-11-08

淮北礦業(yè)科研計劃項目(2019-45，2019-46)；國家重點研發(fā)計劃課題(2018YFC0807804)

Research Project of Huaibei Mining(Group) Co.，Ltd. (2019-45，2019-46)；National Key R&D Program of China (2018YFC0807804)

楊光明，1979年生，甘肅天水人，高級工程師，從事地震數(shù)據(jù)處理工作. E-mail：yangguangming@cctegxian.com

楊光明，金學良，張憲旭，等. 寬頻寬方位處理技術在淮北礦區(qū)全數(shù)字高密度地震勘探中的應用[J]. 煤田地質與勘探，2020，48(6)：55–63.

YANG Guangming，JIN Xueliang，ZHANG Xianxu，et al. Application of broadband and wide azimuth processing technology in full digital high density seismic exploration in Huaibei mining area[J]. Coal Geology & Exploration，2020，48(6)：55–63.

(責任編輯 聶愛蘭)