高速巖體反射地震探測中的數據采集與處理

李 穩 郭雯雯 侯黎華 酆少英 劉保金

1 中國地震局地球物理勘探中心, 鄭州市文化路75號, 450002 2 中國科學院地質與地球物理研究所, 北京市北土城西路19號, 100029 3 鄭州財經學院信息工程學院, 鄭州市天河路36號, 450049

?

高速巖體反射地震探測中的數據采集與處理

李 穩1,2郭雯雯3侯黎華1酆少英1劉保金1,2

1 中國地震局地球物理勘探中心, 鄭州市文化路75號, 450002 2 中國科學院地質與地球物理研究所, 北京市北土城西路19號, 100029 3 鄭州財經學院信息工程學院, 鄭州市天河路36號, 450049

為了調查內蒙古阿拉善某高放廢物地質處置庫預選場址中目標花崗巖體的空間展布、內部結構和邊界接觸關系，應用人工反射地震探測技術，采集數據并進行處理，獲得了針對大型高速巖體的近地表速度剖面和地震反射波法成像剖面。針對崎嶇地表工區近地表速度建模與層析靜校正、強散射干擾波壓制等問題，引進基于MSFM射線追蹤算法的初至波旅行時層析反演與層析靜校正技術、改進的矢量分解壓噪方法等，并在常規高分辨反射地震數據處理流程的基礎上形成了適合于高速巖體地震探測的精細數據處理流程。

高放廢物地質處置庫；高速巖體；反射地震勘探；初至波層析反演；矢量分解法

對于高放廢物的最終處置，目前普遍接受的方案是把高放廢物埋置在距離地表500～1 000 m的地質體中[1]。深部地質體的完整性與穩定性是選址工作中的一項關鍵因素。為了對高放廢物地質處置備選場址進行準確全面的評價，需要進行構造地質、水文、地球物理、鉆探等多方面的研究，地震勘探是其中一項重要內容。與在沉積區開展的地震勘探工作[2-3]不同，內蒙古阿拉善某高放廢物地質處置庫選址中的地震探測是以大型高速花崗巖體為探測目標，在數據采集和后續資料處理中始終面臨著地表崎嶇、近地表低速層帶橫向變化劇烈和高速巖體頂界面埋深極淺甚至直接出露于地表等諸多不利因素的影響。為了解決這些問題，本次地震勘探在野外數據采集和資料處理階段引進了一些具有針對性的技術。

1 工作區概況與地震測線位置

地震探測工作區位于內蒙古自治區西部，巴丹吉林沙漠東緣，海拔1 300～1 500 m。工作區所在位置見圖1。圖中，紅色線條為結合最新地質研究成果和衛星圖像解譯出的推測斷層，需要進一步的地球物理探測工作和詳細的地質勘查工作來核實。

圖1 推測斷層及地震測線位置分布圖Fig.1 Might exist faults and the location of seismic lines

根據已發表的研究成果[4]，測區內出露的地層主要包括前寒武紀變質表殼巖系、中下侏羅統芨芨溝組(J1-2j)、下白堊統巴音戈壁組(K1b)和第四系。測區總體上可分為西北部的晚古生代巖漿活動帶和東南部的中新生代盆地兩大地質單元(圖2)。中新生代盆地主要為侏羅-白堊紀的砂礫巖沉積；晚古生代巖漿活動帶主要分布著晚古生代的侵入巖，在巖漿巖帶中殘留著前寒武紀變質結晶基底及早古生代巖漿巖，并被印支期巖漿侵入。工作區內巖漿巖在空間上多為北東東向延伸，褶皺構造不發育，不同延伸方向的斷裂構造及韌性剪切帶是區內主要的構造形跡。

根據探測目標、地質任務，并經過現場踏勘，布設了2條近平行的地震測線，分別命名為塔木素1線(TMS1)和塔木素2線(TMS2)，在圖1、圖2中以藍色線條標示。兩條地震測線均自南向北探測。其中，TMS1線(40.712 961 4°N、103.131 005 3°E～40.797 889 9°N、103.085 105 8°E)起于新呼都格附近大范圍的三疊紀花崗巖(Tηγ)出露區，穿過巴嘎·額爾崩北部的第四系沖洪積物(Qhpl+eol)覆蓋區，以及海子·阿德爾根韌性剪切帶，總長度10.196 km。TMS2線(40.648 470 1°N、103.258 013 3°E～40.794 809 5°N、103.205 057 1°E)起于測區東南部中新生代盆地內部，穿過陶來阿木韌-脆性變形帶、 扣爾根·呼都格一帶大范圍的晚二疊世花崗巖(P3γδ)出露區、呼和·溫多爾一帶三疊紀花崗巖(Tηγ)出露區、地表基本被第四系上更新統洪積物(Qp3pl)覆蓋的伊和圖克木白堊紀盆地的西南角以及盆地西側的呼和·溫多爾-查干陶勒蓋韌性剪切帶，總長度16.856 km。

根據區域構造研究成果，塔木素2線橫跨測區東南部的中新生代盆地和西北部的晚古生代巖漿活動帶兩大地質單元，沿測線地表施工條件復雜，其探測成果能夠體現目標巖體的南部邊界、風化層或蓋層以下巖體穩定部分的頂界、巖體內部精細結構(尤其是有無斷層、破裂面的存在以及它們的性質)，能夠更好地反映在崎嶇地表工作區開展以高速巖體為探測目標的地震探測工作的特點，因此將其作為本文重點討論的內容。

1.全新統沖洪積物Qhpal 2.全新統洪積+風積物Qhpl+eol 3.上更新統洪積物Qp3pl 4.下白堊統巴音戈壁組K1b 5.中下侏羅統芨芨溝組J1-2j 6.三疊紀中細粒二長花崗巖Tηγ 7.早三疊世中細粒含斑二長花崗巖T1ηγ 8.晚二疊世中粒似斑狀二長花崗巖P3ηγπ 9.晚二疊世中細粒花崗閃長巖P3γδ 10.晚二疊世中粒含斑石英閃長巖P3δο 11.晚二疊世中細粒閃長巖P3δ 12.中二疊世二云母花崗巖P2γ 13.中二疊世似斑狀二長花崗巖P2ηγπ 14.中二疊世片麻狀石英閃長巖P2δο 15.早石炭世片麻狀二長花崗巖C1ηγ 16.扎蓋圖片麻巖Zhgn 17.前寒武片巖巖組Pt3sch 18.前寒武角閃斜長片麻巖巖組Pt3bpg 19.前寒武黑云斜長片麻巖巖組Pt3phg圖2 研究區地質構造綱要圖與地震測線位置Fig.2 Geological structure map of study area and the location of seismic lines

2 觀測系統設計與地震數據采集

2.1 地震儀器與地震觀測系統

本次地震探測使用法國Sercel公司生產的SN408數字地震儀。地震檢波器采用固有頻率為60 Hz的高頻檢波器，以3串1并的方式連接。從提高地震資料信噪比和分辨率的角度出發并兼顧環保和安全問題，地震波激發源選用了M18-612型大噸位可控震源，工作方式為連續變頻掃描，掃描頻率20～160 Hz，掃描長度12 s。

除了儀器因素，合理的地震觀測系統對獲得好的探測結果同樣至關重要。由于是首次對研究區開展地球物理探測工作，施工時沒有可參考的地下速度結構研究成果，因此遵循“密集點距、密集炮距、多道觀測、高頻激發、寬帶接收”的原則設計地震觀測系統，部分采集參數結合現場試驗確定。最終形成的地震觀測方案為：道間距4 m；接收道數300道；炮間距24 m；采樣間隔1 ms；記錄長度3 s；單邊追逐激發或中間激發雙邊不對稱接收；覆蓋次數≥25次；最小和最大偏移距視不同測段的地質條件和干擾波的發育情況現場確定并適時調整(通常最小偏移距為4～16 m，最大偏移距約為1.2 km)。

如圖3所示，本次地震勘探沿測線地表起伏明顯、高程變化劇烈，屬于崎嶇地表的范疇。當遇到障礙測段時，為了保證對地下反射信息的連續追蹤，需要根據障礙類型、規模、現場激發接收條件等采取合理的改變觀測系統設計(簡稱變觀設計)。在本次野外施工過程中，根據現場障礙情況并結合中國地震局地球物理勘探中心自主研發的淺層地震勘探過障礙變觀設計軟件進行變觀處理，最終形成的TMS2線地震觀測系統綜合平面圖見圖4。

圖3 TMS2線沿測線地表高程變化曲線Fig.3 Surface elevations along TMS2 line

圖4 TMS2線地震觀測系統綜合平面圖Fig.4 Seismic geometry of TMS2 line

圖5為TMS2線上3幅典型的野外單炮地震記錄(TWT為雙程走時)。可以看出，當地震排列位于中新生代盆地上方時(圖5(a))，可接收到大量的反射信息，地震記錄中反射波同向軸豐富；而當地震排列位于高速巖體上方時，如圖5(b)和圖5(c)所示，地震記錄上的反射波信號稀疏而且微弱。這是因為對于較為致密的高速巖體單元，雖然其表面波阻抗較大，但其內部的波阻抗差卻很小，因此在巖體內部能夠形成的反射信息不足。

圖5 典型原始野外單炮記錄Fig.5 Typical original single-shot seismic records

2.2 地震數據處理

對于本次地震探測來說，沿測線地表高程落差大，地貌橫向變化顯著，老地層、巖體頂界面埋深淺甚至直接出露于地表，并且大部分巖體出露部分風化嚴重。因此復雜地表工區近地表速度建模與層析靜校正問題、強散射干擾波壓制問題是本次地震數據處理工作中的重點和難點。我們在常規高分辨反射波法地震數據處理流程的基礎上，引進了基于MSFM(multi-stencils fast marching，多模板快速步進)射線追蹤算法的初至波旅行時層析反演技術，建立了準確、精細的近地表速度模型，并利用所獲得的速度模型成果完成了后續對于陸上復雜地表工區地震勘探至關重要的層析靜校正工作；同時，利用改進的矢量分解壓噪方法在疊前濾除地震散射干擾波。

2.2.1 崎嶇地表工區近地表速度建模與層析靜校正

地震波走時層析成像是解決速度建模問題的重要技術，高效、高精度且能適應復雜模型的射線追蹤算法[5-6]是該技術的關鍵。有別于傳統的試射(Shooting)法和彎曲(Bending)法，MSFM射線追蹤算法是一種將波傳播描述為波前面擴展，通過求解用程函方程描述的波前面演化，再從接收點沿地震波走時梯度方向追蹤得到接收點和震源點之間路徑的射線追蹤算法[7]。利用該算法進行射線追蹤突破了層狀假設，使能夠利用的地震初至波包括直達波、回轉波、折射波或多個層折射波的組合。

圖6 塔木素2線近地表速度層析反演結果Fig.6 TMS2 line near-surface velocity tomography result

圖6為利用該技術得到的TMS2線近地表速度層析反演結果。這一高精度、高分辨的速度建模結果可以直接應用于后續的層析靜校正處理(圖7)，還可以作為進行地震解釋工作時重要的參考資料。

圖7 層析靜校正處理效果Fig.7 Tomographic static correction processing effect

2.2.2 利用改進的矢量分解壓噪方法濾除近地表散射波

地震散射波是當散射障礙體的線度與入射波波長相比很小或相近時，由于干涉效應產生的與散射障礙體幾何形狀的規則性有關的相干或不相干波列[8]。當表層介質小尺度非均勻性強時，易產生強烈的地震散射波干擾，并且往往以眾多散射波互相疊加的復雜形態出現。本次地震探測工作區內不同尺度的非均勻地質體共生，老地層、巖體頂界埋深淺且橫向變化大，沿地震測線戈壁、陡坎以及風化嚴重的巖體露頭隨處可見，因此強散射干擾波在本次獲得的原始地震記錄上普遍存在，嚴重影響著地震數據的品質。根據前人的研究成果，散射波的速度與面波相當，但散射波的頻段位于有效波的頻段范圍內[9]，因此對散射波的處理不能采用簡單的頻率濾波手段。在本次數據處理過程中，利用散射干擾波與有效反射信號傳播方向不同的特點，選用了改進的矢量分解法保幅去噪技術[10-11]在疊前濾除散射噪聲。

圖8為某散射噪聲干擾嚴重的典型單炮地震記錄去噪前、去噪后以及濾出的噪聲圖像。可以看出，經過處理，有效反射波同相軸顯得更加清晰和連續，地震記錄信噪比得到顯著提高。在圖8(c)展示的濾出的噪聲圖像上可以直觀地看到，濾出的噪聲散射特征明顯，并且在噪聲剖面上基本觀察不到有效波同相軸的影像。

圖8 改進的矢量分解壓噪方法濾除地震散射波Fig.8 De-noising scattered wave using the improved vector resolution noise removal approach

在前文基礎上，本次反射地震探測最終形成如下的常規數據處理流程(圖9(a))和精細數據處理流程(圖9(b))。兩者之間的主要區別為：在常規處理流程中，靜校正工作通過高程靜校正+折射靜校正的方法完成，對于近地表散射波干擾未進行針對性處理。在精細處理流程中，通過高精度速度建模與層析靜校正更為合理、有效地解決了陸上復雜工區地震資料的靜校正問題，并對地震散射干擾波進行了重點壓制。特殊處理方法所得的中間處理結果對于改善后續速度分析、疊加、偏移成像等諸多處理環節的效果均有有益影響。

3 數據處理結果對比分析

圖10為對TMS2線地震數據進行常規處理(圖10(a))和精細處理(圖10(b))后得到的反射波時間偏移剖面。圖11是根據圖10(b)以及圖6所示的近地表速度層析反演結果作出的地震解釋剖面圖。可以看出，在常規處理結果中，由于靜校正問題和近地表散射波問題的存在，直接影響了對風化層和蓋層以下巖體穩定部分頂界面的追蹤。剖面淺部來自巖體頂界的反射波組同相軸連續性不佳；對于某些同相軸錯斷的位置(如5 500、9 000 m剖面樁號的下方)，難以判斷其是否預示著小型裂隙的存在。在剖面南端反映中新生代盆地內部結構的剖面部分，反射波同相軸的連續性、清晰度及其所揭示的地層產狀的可靠性等均有待進一步加強。在進行精細處理所得的反射波成像剖面上，由1～2個強能量同相軸組成的來自巖體穩定部分頂界的反射波組相對連續。通過對比可以發現，該反射波組在橫向上的起伏變化形態、連續性特征與圖11(a)所示的高速體頂面特征相一致，并且與圖2所揭示的不同巖性出露區之間的分界存在良好的對應關系。在剖面樁號約4 240 m以南反映中新生代盆地內部結構的剖面部分，圖10(b)也顯示了更好的成像效果。

圖9 數據處理流程圖Fig.9 Data processing flow

圖10 塔木素2線數據處理結果對比Fig.10 The contrast of data processing results of TMS2 line

利用本次地震勘探成果進行地震地質解釋(見圖11)，首先在剖面南端可以可靠地識別出AU1、AU2、和AU3三組不整合構造關系。其中，AU3為分隔了目標巖體與其南側中新生代盆地的目標巖體邊界。在AU1、AU2兩側，上覆與下伏的反射波組之間存在明顯的交角，反映出上覆地層與下伏地層的產狀不一致并以一定的角度相交，這種現象正是角度不整合構造關系在地震剖面上的典型反映。同時，在近地表速度剖面(圖11(a))上，與AU1、AU2、AU3所對應的位置存在梯度特征明顯的速度分界。根據反射波場特征，可在地震剖面上識別出4個斷點，分別以FP1～FP4標示，其位置如圖中紅色虛線所示。此外，還可識別出一組來自巖體深部的弧狀反射波組，以R1標示。其反映了巖體在成巖或之后的過程中受到過擠壓應力環境的作用，該剖面特征與對該區的地質構造背景認識相吻合。

圖11 塔木素2線地震解釋剖面Fig.11 TMS2 line interpretation profile

4 結 語

以大型高速巖體為目標開展反射地震探測工作時，需要根據特殊的地質任務、探測目標以及通常情況下較為不利的施工環境，引進有針對性的數據采集、處理方法。在本次針對內蒙古阿拉善某大型目標花崗巖體的地震探測中，重點解決了野外靈活變觀數據采集、復雜地表工區高精度近地表速度建模與層析靜校正、強散射干擾波壓制等問題，提高了最終地震成像剖面的質量。

本次地震探測實際采用了反射波法地震勘探與初至波旅行時層析反演相結合的工作方法。可以發現，對于高速巖體內部的小型斷層、裂隙，其在反射波成像剖面上的展布形態與在淺層速度剖面上的形態特征之間具有良好的對應關系。可以認為，對于巖體頂界面埋深淺甚至直接出露于地表的大型高速巖體進行地震勘探時，采用反射波法界面成像和初至波速度成像相結合的方法，可以獲得更加準確、可靠的探測與解釋成果。

[1] 王駒,徐國慶,金遠新. 論高放廢物地質處置庫圍巖[J]. 世界核地質科學, 2006, 23(4):222-231(Wang Ju, Xu Guoqing, Jin Yuanxin. On the Host Rock for the Geological Repositories of High Level Radioactive Waste[J]. World Nuclear Geoscience, 2006, 23(4):222-231)

[2] 劉保金,柴熾章,酆少英,等. 第四紀沉積區斷層及其上斷點探測的地震方法技術——以銀川隱伏活動斷層為例[J]. 地球物理學報, 2008, 51(5):1 475-1 483(Liu Baojin, Chai Chizhang, Feng Shaoying, et al. Seismic Exploration Method for Buried Fault and Its Up-Breakpoint in Quaternary Sediment Area-An Example of Yinchuan Buried Active Fault[J]. Chinese J Geophys, 2008, 51(5):1 475-1 483)[3] Feng S Y, Zhang P Z, Liu B J, et al. Deep Crustal Deformation of the Longmenshan, Eastern Margin of the Tibetan Plateau, from Seismic Reflection and Finite Element Modeling[J]. Journal of Geophysical Research: Solid Earth, 2016, 121(2):767-787

[4] 劉治博,張維杰. 內蒙古阿拉善右旗杭嘎勒中二疊世石英閃長巖地球化學特征和LA-ICP-MS鋯石U-Pb定年[J]. 地質學報, 2014, 88(2):198-207(Liu Zhibo, Zhang Weijie. Geochemical Characteristics and LA-ICP-MS Zircon U-Pb Dating of the Middle Permian Quartz Diorite in Hanggale, Alax Right Banner, Inner Mongolia[J]. Acta Geologica Sinica, 2014, 88(2):198-207)

[5] 鄒長橋,段永紅,唐秋惠,等. Ronge-Kutta法在計算復雜構造區域地震波旅行時中的應用[J]. 大地測量與地球動力學, 2015, 35(5):894-897(Zou Changqiao, Duan Yonghong, Tang Qiuhui, et al. The Applications of Ronge-Kutta Method to Seismic Wave Traveltimes in Complex Structural Model[J]. Journal of Geodesy and Geodynamics, 2015, 35(5):894-897)

[6] 孫章慶,孫建國,韓復興．針對復雜地形的三種地震波走時算法及對比[J]. 地球物理學報, 2012, 55(2):560-568(Sun Zhangqing, Sun Jianguo, Han Fuxing. The Comparison of Three Schemes for Computing Seismic Wave Travel Times in Complex Topographical Conditions[J]. Chinese J Geophys, 2012, 55(2):560-568)

[7] 盧回憶,劉伊克,常旭. 基于MSFM的復雜近地表模型走時計算[J]. 地球物理學報, 2013, 56(9):3 100-3 108(Lu Huiyi, Liu Yike,Chang Xu. MSFM-Based Travel-Times Calculation in Complex Near-Surface Model[J]. Chinese J Geophys, 2013, 56(9):3 100-3 108)

[8] 李燦蘋,劉學偉,李敏鋒. 非均勻地質體散射波特征初探[J]. 石油物探, 2006, 45(2):134-140(Li Canping, Liu Xuewei, Li Minfeng, et al. Preliminary Study on Seattering Wave Characteristics of Heterogeneous Geologic Bodies[J]. Geophysical Prospecting for Petroleum, 2006, 45(2): 134-140)

[9] 李戰業,尹軍杰,王赟. 地震散射波模擬成像在金屬礦勘探中的應用[J]. 地質與勘探, 2009, 45(2):80-84(Li Zhanye, Yin Junjie, Wang Yun. The Application of the Numerical Simulation and Imaging of Seismic Scattered Wave in Exploration of Metal Mineral Deposits[J]. Geology and Exploration, 2009, 45(2):80-84)

[10]王宏偉. 矢量分解法壓噪[J]. 石油地球物理勘探, 1989, 24(1):16-29(Wang Hongwei. Noise Suppression Using Vector Resolution Method[J]. Oil Geophysical Prospecting, 1989, 24(1):16-29)

[11]何銀娟,李穩,劉保金等. 改進的矢量分解壓噪方法[J]. 石油地球物理勘探, 2015, 50(2):243-253(He Yinjuan, Li Wen, Liu Baojin, et al. Improved Vector Resolution Noise Reduction Method[J]. Oil Geophysical Prospecting, 2015, 50(2):243-253)

Reflection Seismic Data Acquisition and Processing Technology Used in Detecting the High-Velocity Rock Body

LIWen1,2GUOWenwen3HOULihua1FENGShaoying1LIUBaojin1,2

1 Geophysical Exploration Center, CEA, 75 Wenhua Road, Zhengzhou 450002, China 2 Institute of Geology and Geophysics, CAS, 19 Beituchengxi Road, Beijing 100029, China 3 College of Information Engineering,Zhengzhou Institute of Finance and Economics, 36 Tianhe Road, Zhengzhou 450049, China

In order to investigate the spatial distribution, internal structure, and boundary contact relations of target granite rock at Inner Mongolia Alxa pre-selected high-level radioactive waste geological)

high-level radioactive waste geological repository; high-velocity rockc body; reflection seismic exploration; first-arrival wave tomography inversion; vector resolution noise removal

National Natural Science Foundation of China, No. 40874045,41374100.

LIU Baojin, researcher, PhD supervisor, majors in exploration geophysics and seismology,E-mail: LBJ001@126.com.

2016-04-27

項目來源：國家自然科學基金(40874045,41374100)。

李穩,博士生,主要從事地震勘探數據處理方法和地震資料解釋研究，E-mail: liwen_work@163.com。

劉保金,研究員,博士生導師,主要從事勘探地球物理和地震學研究，E-mail: LBJ001@126.com。

10.14075/j.jgg.2016.12.001

1671-5942(2016)012-1035-06

P315；P631

A

repository site, seismic detecting work is carried out. Near-surface velocity tomography inversion profiles and seismic reflection profiles are obtainedby acquiring and processing the seismic data.This seismic exploration work aimsat technical problems such as ruggedness surface work area near-surface velocity model-building, tomographic statics correction, and strong seismic scattered wave pressing. To this end, we adoptakind of MSFM(multi-stencils fast marching) based first-arrival travel-time tomography inversion technology, and an improved vector resolution noise removal approach.Meanwhile, based on routine high resolution reflection seismic exploration data processing flow, a set of customized fine data processing flow is set up.

About the first author:LI Wen, PhD candidate, majors in seismic data processing and seismic interpretation,E-mail: liwen_work@163.com.