摘 要

摘　要：海底高頻地震儀(HF-OBS)可獲得含天然氣水合物(以下簡稱水合物)層的縱橫波速度、波阻抗等多種有效地震反射信息，對準確估算水合物的資源量具有重要意義。根據我國南海北

摘　要：海底高頻地震儀(HF-OBS)可獲得含天然氣水合物(以下簡稱水合物)層的縱橫波速度、波阻抗等多種有效地震反射信息，對準確估算水合物的資源量具有重要意義。根據我國南海北部海底水合物HF-OBS勘探實例，基于射線追蹤法、雙聚焦法、波動方程的正演模擬法3種方法，對HF-OBS觀測系統進行了優化設計，并利用縱橫波聯合走時反演獲得了該區含水合物層的縱橫波速度關系。結論認為：①采用射線追蹤法模擬目的層射線分布和疊加次數，同時結合波動方程的正演模擬結果，可以有效地進行觀測系統優化設計；②雙聚焦理論的參數模擬能夠進一步驗證觀測系統設計的合理性，針對特定的水合物目標體，為了獲取礦體內部精細的速度結構，HF-OBS的分布形態應采用長方形或正方形的矩陣分布，分布間距以300～500m為宜；③優化后的觀測系統有利于獲得高質量的HF-OBS地震數據，而縱橫波聯合反演獲得的速度信息減小了水合物識別的不確定性。

關鍵詞：中國南海  海底高頻地震儀  觀測系統  優化設計  天然氣水合物  聲波速度  聯合走時反演  分布形態

An optimal design of a high-frequency oceali bottom seismometer(HF-OBS) and its application to the natural gas hydrate exploration in the South China Sea

Abstract：A high-frequency ocean bottom seismometer(HF-OBS)can not only receive multi-component displacement records of the bottom of the sea，but acquire various seismic information of P wave and converted S-wave，which is of great significance to accurately evaluate the hydrate reserves of localities．Therefore，in a case study of natural gas hydrate exploration in the Sou th China Sea．we made an optimal design of the HF-OBS geometry system through three forward modeling simulation methods of ray tracin9，focal beam analysis，and wave equation algorithms．And we adopted a hybrid inversion method of wave-equation travel time inversion(WT)and waveform inversion to obtain the relationship between P and S-wave velocity of hydrate layers in this study area．The following findings were concluded．(1)The ray tracing method simulates the Gamma distribution and stacking folds(superimposition times)of pay zone and integrates the results of wave equation algorithms to effectively improve the optimal design of HF-OBS geometry system．(2)Simulation parameters of the focal beam analysis can further validate the reasonability of the observing system design and as for the specific hydrate target，the matrix distribution of a rectangle or a square with the space between of 300-500m should be employed for the HF-OBS to acquire the very fine inner velocity structure of the hydrate target．(3)The optimized observation system is conducive to obtain HF-OBS seismic data,while velocity obtained by p-and S-wave joint inversion decreases the uncertainty of gas hydrate recognition．

Keywords：South China Sea，HF-OBS，layout design，natural gas hydrate,acoustic velocity

海洋天然氣水合物(以下簡稱水合物)的識別，通常利用地震反射特征——如水合物賦存的四大地震異常特征(即BSR、BZ、速度異常及極性反轉)來確定其分布，而含水合物層物性參數主要通過地震反射信息間接獲取。由于水合物縱波速度比海洋正常沉積物要高，致使水合物沉積物的體積模量增加；同時，沉積物孔隙內的水合物膠結物也改變了沉積物的剪切模量，沉積物內聲波阻抗亦呈“譜白化”現象。因此，隨著勘探程度的深入，需要高精度的速度和高分辨率地震波形資料以獲得含水合物沉積物的空間分布及其物性參數，對水合物的資源量進行估算。由于水合物海底高頻地震儀(HF-OBS)觀測具有震源與檢波器間的距離大、曳航方式獲得低角度反射波、能夠較直觀地觀測BSR反射波波形等優點，因此，利用HF-OBS可獲得含水合物地層的縱波速度(vp)、橫波速度(vs)、速度結構及波阻抗等的變化，對準確估算水合物的資源量具有重要意義^[1-2]。此外，海底高頻地震儀調查獲取的多種資料還可滿足波形反演、走時反演、非線性全波形走時聯合反演等各種研究的需要。目前HF-OBS在我國已經由研制和試驗階段逐漸轉化為成熟應用階段^[3-4]，近年來廣州海洋地質調查局采用HF-OBS探測技術，對水合物進行了詳細勘探和深入研究^[5-7]。

1　HF-OBS結構及工作原理

1．1　HF-OBS結構和勘探特點

20世紀90年代以前，以拖纜為主體的人工地震探測方法為探測海底地球物理信息的主流方法。但是，在復雜多變的海洋環境中，海水的存在無疑產生了諸多不利影響，例如，微弱的地震信息受到海水層的作用而產生明顯衰減，從而導致信號畸變等。此外，拖纜采集無法獲得有效的橫波信息，從而難以滿足地層速度結構等研究。隨著海洋地震學的誕生與發展，對海洋巖石圈與地殼結構、地球起源、地震起因乃至海洋資源與環境等研究迫切需要一種能夠接收海底地震縱橫波的高保真度、大動態范圍的探測儀器。因此，海底地震儀(Ocean Bottom Seismometer，OBS)應運而生，從而揭開了海洋勘探OBS技術研究的序幕。海底高頻地震探測系統(High Frequency Ocean Bottom SeiSmometer，HF-OBS)，因其具有較寬的頻寬和良好的高頻響應特性，在天然氣水合物勘探方面受到越來越多的關注和重視^[8]。

天然氣水合物海底高頻地震儀(HF-OBS)除了需要具備較寬的勘探頻帶與良好的高頻響應特性之外，還要求儀器與海底之間具有良好的耦合特性。因此，儀器設計主要圍繞以下3個方面開展設計：①高采樣率的模數轉換模塊設計；②姿控檢波器的沒計與性能優化；③海底耦合的研究與HF-OBS的結構設計^[9-10]。其結構由塑料殼玻璃儀器艙、脫鉤機構、沉耦架3部分組成(圖1)。其中，塑料殼玻璃儀器艙是海底地震儀的主體部分，也是其可回收部分，它內部采用單球一體化結構，主要由水聲通訊模塊、檢波器與常平裝置、采集系統的電子電路、地震儀電源、數據提取模塊、無線數傳模塊及電子電路版頻閃等構成。

1．2　HF-OBS工作原理

如圖l-a所示，HF-OBS共有4個地震通道，其中1個通道為水聽器(只能接收縱波，即H分量)，另外3個通道為三分量動圈式速度換能器(位于儀器底部與海底接觸，可接收橫波)，工作頻帶介于5～250Hz，采集器動態范圍超過110dB；地震采樣率不低于1ms；在2ms采樣率下，在水下的連續工作時間可超過60d；采用內置羅盤進行定向，聲學測距設備實現水下定位；工作水深超過4000m，滿足我國南海海域天然氣水合物地震數據采集的需要^[11]。

 

三維拖纜地震與HF-OBS的聯合勘探在近年來水合物地震勘探中經常用到^[12]，其目的是綜合利用縱、橫波信息詳細了解水合物礦體的內部結構。根據勘探目的，將多個HF-OBS以不同的排列方式置于海底目標點以采集地震信號，達到拖纜和HF-OBS聯合記錄縱橫波地震反射信息的目的，其典型的采集方式見圖1-b。為保證地震儀能位于期望位置，通常將儀器隨鋼纜下放至距海底一定深度后，使儀器自由下落至海底。在完成采集后通過指令使地震儀脫鉤上浮^[13]。

2　HF-OBS觀測系統優化設計

水合物勘探區海底地形地貌情況通常較復雜，HF-OBS投放數量有限，為了獲取地層的縱波、轉換橫波等多種有效的地震波信息，對HF-OBS站位的分布形態、分布間距等技術參數選擇非常重要。因此，在海上資料采集前，需要針對水合物探區的海底地質條件，進行海底地震觀測系統模擬。

圖2為南海北部陸坡測線A的地震剖面，剖面上BSR非常清晰、連續，振幅空白帶、與地層斜交等特征明顯。為了進一步弄清水合物礦體的內部速度特征，觀測系統設計以突出地震測線所穿過的水合物礦體作為目標。

 

2．1　射線追蹤法

射線追蹤法是在以往水合物調查的基礎上，針對特定的水合物目標體采用“二維擬三維”的方法建立三維數據體，然后沿特定方向抽取地震剖面，模擬海上采集的作業環境，分析HF-OBS分布區間的射線數量，最終確定HF OBS的分布形態與分布間距，實現觀測系統設計。以地震剖面上水合物為特定目標，建立相應的地質模擬模型，模擬野外采集的作業環境，計算不同測線段內覆蓋次數以及均勻程度，分析HF-OBS分布區間的射線數量，從而最終確定HF-OBS的形態間距與分布形態，實現觀測系統優化設計。

通過目標體的二維地震剖面，按HF-OBS分布間距，通過射線追蹤模擬，分析縱、橫波接收情況，根據覆蓋次數與均勻程度，選取最佳的覆蓋效果(圖3)。圖中藍色表示水合物目標體(剖面)，射線所指的紅色區域表示所對應的反射轉換波射線分布區域。通過對HF-OBS分布間距為100m、200m、300m、400m、500m、600m、700m和1000m等情況下模擬得到的幾組數據的分析比較，分別計算落入不同距離的線段內的疊加次數，并用色標進行表示，然后分析疊加次數的均勻程度和覆蓋效果。當HF-OBS分布間距在500m以下時(300～400m效果最佳)，線段內的疊加次數達超過30次，并且疊加次數的分布比較均勻。最后，確定HF-OBS臺站間距小于500m時效果較好，且分布間距愈小，所使用的HF-OBS臺數愈多，疊加次數愈高，覆蓋效果愈好。

 

因此，針對天然氣水合物目標體，結合軟件模擬和實際投放時海流等因素所造成的投放誤差的影響，最終確定針對特定目標體的HF-OBS最佳觀測系統方案如下：分布間距應該小于500m，300～400m效果最佳；分布方式為長方形或正方形的矩陣分布方式。

2．2　雙聚焦法

雙聚焦法以雙聚焦理論為基礎，通過評價主瓣分辨率、主瓣能量、旁瓣能量及主旁瓣能量比等設計觀測系統的一種方法。通過對射線追蹤法和雙聚焦法兩種方法所設計的觀測系統進行對比驗證，最終合理設計海底地震探測觀測系統^[14]。根據雙聚焦理論^[15]，震源在目標體周圍一定范圍內激發，將檢波器置于目標點。因此，可以用震源陣列聚焦束評估觀測系統中震源對目標點的成像分辨率的影響。具體做法是：①選取不同的炮點間距、測線間距(炮線間距)，進行聚焦模擬；②選取不同的HF-OBS分布間距，進行聚焦模擬；③根據分析結果，確定最佳的觀測系統。

2．2．1評價準則

基于雙聚焦理論研制的軟件評價準則，主瓣分辨率要高(分辨率曲線中最大極值點與左右第一極小值點間距離大者表示分辨率)；主瓣能量大，且聚焦點的位置應為目標點的位置；主旁瓣能量比(主瓣能量與旁瓣能量的比值)大。因此，設計HF-OBS觀測系統時，首先將震源陣列置于目標點位置，HF-OBS當成接收點分布于目標點周圍的海底附近，接收地震信號，當聚焦波束主瓣能量大且波束位置不變，旁瓣能量小且旁瓣數少時，表示HF-OBS的分布間距合理，從而達到HF-OBS觀測系統設計的目的。

2．2．2聚焦模擬

在以往調查的基礎上，建立水合物目標層與相鄰地層簡化的地球物理模型，將BSR簡化為平行海底的地層。假定炮間距與線間距分別為25m(與野外采集一致)，然后根據不同HF-OBS分布間距聚焦模擬。

2．2．3參數確定

圖4為聚焦模擬效果圖。當HF-OBS分布間距為500m時，出現明顯旁瓣，但模擬震源的位置并未發牛變化(圖4-b)；當分布間距為750m時，模擬震源的位置發生變化，成像效果明顯變差(圖4-c)；當分布間距變小時，模擬成像效果愈來愈好(圖4-a)。

 

在南海北部陸坡進行水合物三維地震與海底地震聯合采集試驗時，由于試驗海域海流變化大、速度快，HF-OBS投放位置精度誤差較大。據試驗統計，如果HF-OBS采用海面投放方式，位置誤差為l00～500m。另外，觀測系統設計時，需要考慮HF-OBS臺站不可能太多。綜合上述因素，雙聚焦理論的模擬結果表明：HF-OBS分布間距應小于500m，并且愈小愈好。

2．3　波動方程的正演模擬法

采用正演模擬方法首先根據研究區資料建立地質模型，然后將震源置于目的層中，根據波動理論正演合成地震記錄，在地表接收地震記錄，再將接收到的地震記錄成像到震源位置。根據正演結果，可以模擬HF-OBS接收信息的地層效果，用于指導野外采集。

以南海北部陸坡地震測線A所對應的二維地震剖面作為波動正演模型的基礎，炮點位于模型中間的海面上，在海上和海底同時接收。圖5為波動方程正演模擬的結果與HF-OBS實際接收資料的效果對比圖。正演模擬結果可以幫助分析野外HF-OBS數據的接收效果，對野外采集具有較好的指導作用。

 

3　效果分析

根據觀測系統設計研究結果，廣州海洋地質調查局于2006年在南海北部陸坡利用HF OBS進行了采集試驗，研究區水深介于1000～1750m。試驗時采用三維地震與海底地震聯合采集方式：把HF-OBS投放到具有明顯BSR顯示的海底位置，以提高采集地震數據的針對性。

試驗一共投放了20臺HF-OBS；針對可能的水合物礦體目標，HF-OBS采用長方形矩陣(5×4陣列)進行布設，分布間距為500m；激發震源為540in³(1in＝25.4mm)的GI槍組合“點震源”系統，炮間距為25m，線間距為25m(圖6)。考慮到HF-OBS與海底的耦合度及靜止于海底后的穩定性，施工時HF-OBS的投放點選擇相對平緩的區域，以盡可能消除海流的影響；而且投放時以鋼纜控制HF-OBS下降，到達一定深度后讓其在水中自由下落至海底，因此，每個HF-OBS在海底的實際位置及其間隔與設計位置會不一致。試驗結束后，實際回收l3臺HF-OBS。

 

圖7為拖纜采集縱波數據與HF-OBS采集轉換橫波數據成像對比。盡管由于HF OBS數量有限并且其實際位置不可控等因素造成縱波的整體成像效果優于橫波，但在沉積物物性發生變化的界面，橫波成像數據表現出更明顯的異常特征。

 

對比可以發現，BSR在橫波成像數據的強振幅特征明顯，也更容易識別。特別是在氣體運移通道位置，縱波成像數據表現為一個模糊反射區域，而橫波成像數據可以確定明顯的反射界面，這也是橫波信息不受流體干擾最好的證明，其成像結果對于水合物形成所需游離氣是如何運移至水合物穩定帶以及整個水合物儲層系統的地質模式有更充分的認識^[16-17]。

盡管橫波的整體成像效果尚無法與縱波相比，但其走時信息可用于縱橫波速度聯合反演以獲得HF-OBS所在位置較可靠的縱波和橫波速度，進而進行定量分析并確定水合物存在與否。圖8為HF-OBShl7站位的速度反演結果。縱橫波速度在BSR上部明顯增大，縱波速度值介于1788～2742m／s，增幅為420～1090m／s；橫波速度介于400～600m／s，增幅為200～260m／s，反映水合物礦體的存在。淺部的縱波速度急劇下降而橫波幾乎沒有變化，表明此處存在淺部游離氣或游離氣運移通道。從圖8可以看出，HF-OBS反演速度的垂向分辨率比傳統的速度譜分辨率有了很大的提高，充分揭示了天然氣水合物的縱、橫波速度結構，并且速度結構異常與BSR及水合物礦體的分布層位有很好的對應關系^[18]，從而增強了對水合物礦體的分辨能力。

 

4　結論

1)采用射線追蹤法模擬目的層射線分布和疊加次數，同時結合波動方程正演模擬顯示剖面的合成效果，可以進行有效的觀測系統的優化設計。雙聚焦理論的參數模擬能夠進一步驗證觀測系統設計的合理性；針對特定的天然氣水合物目標體并結合實際地質情況，為了獲取礦體內部精細的速度結構，HF-OBS的分布形態應采用長方形或正方形的矩陣分布方式，分布間距以300～500m為宜。

2)通過HF-OBS探測獲取了含水合物地層的縱波、轉換橫波等多種有效的地震波信息，野外采集數據質量較高。

3)HF-OBS獲得的橫波信息是對常規地震資料的有效補充，可以減小常規地震資料解釋的不確定性。

 

參考文獻

[1]KORENAGA J，HOLBROOK W S，SINGH S C，et al．Natural gas hydrates on the Southeast U．S．margin：Constraints from full waveform and travel time inversions of wide angle seismic data[J]．Journal of Geophysical Research，l997，102(B7)：15345-15365．

[2]GREVEMEYER I，ROSENBERGER A，VILLINGER H．Natural gas hydrates on the continental slope off Pakistan：Constrains from seismic techniques[J]．Geophysical Journal International，2000，140(2)：295-310．

[3]吳志強，文麗，煮思友，等．海域天然氣水合物的地震研究進展[J]．地球物理學進展，2007，22(1)：218-227．

WU Zhiqiang，WEN Li，TONG Siyou，ct al．Advance in seismic researches on natural gas hydrate in ocean[J]．Progress in Geophysics，2007，22(1)：218-227．

[4]張光學，黃永樣，陳邦彥，等．海域天然氣水合物地震學[M]．北京：海洋出版社，2003．

ZHANG Guangxue，HUANG Yongyang，CHEN Bangyan，et al．Marine gas hydrate seismology[M]．Beijing：Ocean Press，2003．

[5]郝小柱，伍忠良，曾憲軍．海底地震儀及其發展現狀[J]．中國科技成果，2008，9(17)：33-36．

HAO Xiaozhu，WU Zhongliang，ZENG Xianjun．HF-OBS and its developing situation[J]．China Science and Technology Achievements，2008，9(17)：33-36．

[6]王偉巍．海底地震勘探最新方法與技術發展[J]．中國科技成果，2010，11(6)：7-11．

WANG Weiwei．Marine seismic prospecting method and the latest technology development[J]．China Science and Technology Achievements，2010，11(6)：7-11．

[7]欒錫武，趙克斌，孫冬勝，等．海域天然氣水合物勘測的地球物理方法[J]．地球物理學進展，2008，23(1)：210-219．

LUAN Xiwu，ZHAO Kebin，SUN Dongshen9，et al．Geophysical methods for marine gas hydrates exploration[J]．Progress in Geophysics，2008，23(1)：210-219．

[8]FORSYTH D，DETRICK B．Ocean mantle dynamics science plan[R]．Woods Hole：Woods Hole；ceanographic Institution，2000．

[9]郭建，劉光鼎．無纜存儲式數字地震儀的現狀及展望[J]．地球物理學進展，2009，24(5)：1540-1549．

GUO Jian，LIU Guangding．Current situation and expectation of cable less seismic acquisition system[J]．Progress in Geophysics，2009，24(5)：1540-1549．

[10]曾寶宏．多源多纜深海三維地震勘探技術[M]//中國海洋石油“三新、三化”優秀論文集．廣州：中國海洋石油南海東部公司，1997．

ZENG Baohong．Multiple source and cable 3D deep seismic exploration teChnology[M]//CNOOC“Sanxin，Sanhua” Excellent Essays．Guangzhou：China Offshore Oil Nanhai East Corporation，1997．

[11]羅文造，韋成龍，王立明，等．海上地震勘探主要采集參數的選取與驗證[J]．熱帶海洋學報，2009，28(4)：93-101．

LUO Wenzao，WEI Chenglong，WANG Liming，et al．Selection and validation of the main acquisition parameters in marine seismic exploration[J]．Journal of Tropical Oceanography，2009，28(4)：93-101．

[12]CORDSEN A，GALBRAITH M．Narrow-versus wide-azimuth land 3D seismic surveys[J]．The Leading Edge，2002，21(8)：764-770．

[13]張光學，徐華寧，劉學偉，等．三維地震與HF-OBS聯合勘探揭示的神狐海域含水合物地層聲波速度特征[J]．地球物理學報，2014，57(4)：1169-1176．

ZHANG Guangxue，XU Huaning，LIU Xuewei，et al．The acoustic velocity characteristics of sediment with gas hydrate recovered from integrated exploration of 3D seismic and HF-OBS data in Shenhu Area[J]．Chinese Journal of Geophysics，2014，57(4)：1169-1176．

[14]曾憲軍，伍忠良，郝小柱．天然氣水合物礦體的三維地震與海底高頻地震聯合采集技術一海底地震儀觀測系統研究[J]．海洋技術，2010，29(1)：124-130．

ZENG Xianjun，WU Zhongliang，HAO Xiaozhu．Three dimensional seismic and OBS joint acquisition technology for gas hydrate investigation[J]．Ocean Technology，2010，29(1)：124-130．

[15]VOLKER A W F．Analysis and optimization of 3D seismic acquisition geometries by focal beams[C]//2000 SEG Annual Meeting，6-11 August 2000，Calgary，Alberta，Canada．Tulsa：American Society of Exploration Geophysicists，2000．

[16]何家雄，顏文，祝有海，等．南海北部邊緣盆地生物氣／亞生物氣資源與天然氣水合物成礦成藏[J]．天然氣工業，2013，33(6)：121-134．

HE Jiaxiong，YAN Wen，ZHU Youhai，et al．Bio-genetic and sub biogenetic gas resource potential and genetic types of natural gas hydrates in the northern marginal basins of South China Sea[J]．Natural Gas Industry，2013，33(6)：121-134．

[17]梁永興，曾濺輝，楊智峰，等．南海北部神狐海域輸導體系特征及其對天然氣水合物成藏的影響[J]．地球科學與環境學報，2013，3 5(4)：30-38．

LIANG Yongxing，ZENG Jianhui，YANG Zhifeng，et al．Characteristics of passway systems and their influence on gas hydrate accumulation in Shenhu Area of Northern South China Sea[J]．Journal of Earth Sciences and Environment，2013，35(4)：30-38．

[18]梁勁，王明君，陸敬安，等．南海北部神狐海域含天然氣水合物沉積層的速度特征[J]．天然氣工業，2013，33(7)：29-35．

LIANG Jin，WANG Mingjun，LU Jing^，an，et al．Characteristics of sonic and seismic velocities of gas hydrate bearing sediments in the Shenhu Area，Northern South China Sea[J]．Natural Gas Industry，2013，33(7)：29-35．

 

 

 

本文作者：沙志彬  鄭濤  張光學  劉學偉  伍忠良  梁金強  蘇丕波  王靜麗

作者單位：中國地質大學(武漢)資源學院

  國土資源部海底礦產資源重點實驗室·廣州海洋地質調查局

  中國地質大學(北京)地球物理與信息技術學院