安振芳，張 進，張建中

1.海底科學與探測技術(shù)教育部重點實驗室，山東 青島 266100 2.中國海洋大學海洋地球科學學院，山東 青島 266100

0 引言

隨著經(jīng)濟和社會的普遍發(fā)展，人類對能源的消耗與日俱增，煤、石油和天然氣等常規(guī)能源儲量逐日遞減，因此迫切需要為我們的后代尋找可替代能源。眾所周知，地球上70%的面積被海水所覆蓋，海底蘊含著豐富的礦物能源，比如錳結(jié)核、熱液礦床與天然氣水合物等。據(jù)有關(guān)資料統(tǒng)計，全球天然氣水合物總量大約相當于已探明傳統(tǒng)化石燃料碳總量的兩倍[1-2]。在標準狀況下，單位體積天然氣水合物可釋放出160～180體積的天然氣[3-4]。天然氣水合物被認為是未來人類賴以生存和發(fā)展的高效清潔能源，受到世界范圍內(nèi)科研人員的普遍關(guān)注。

現(xiàn)今常規(guī)海洋三維地震勘探主要采用一船多纜或多船多纜的方法進行資料采集，在大多數(shù)情況下水平纜能夠滿足成像的基本要求。然而，對于深水崎嶇海底或地下復(fù)雜構(gòu)造，水平纜的成像精度卻不夠理想[5-6]。另外，在采集過程中不可避免地會受到多次波的干擾，且在處理階段無法有效剔除。海洋垂直纜(vertical cable，VC)起源于海洋變偏移距垂直地震剖面(walkaway vertical seismic profiling，WVSP)，它是通過美國海軍用于反潛作戰(zhàn)而發(fā)展起來的一項技術(shù)[7]。海洋VC的底端固定在沉放于海底的錨上，頂端系在漂浮于海面的浮球上，以確保電纜能夠處于垂直狀態(tài)[8]。據(jù)報道[9]，1987年在墨西哥灣首次使用海洋VC進行了地震數(shù)據(jù)采集，用于鹽丘側(cè)翼砂巖儲層的成像。1997年，Moldoveanu等[10]在墨西哥灣兩個小的3D工區(qū)，在水平纜和裝有無線電遙測系統(tǒng)的海洋VC之間進行了一次對比，結(jié)果顯示，海洋VC衰減交混回響的能力更強，信噪比和分辨率更高。實踐[11-15]證明，海洋VC可以識別自由表面多次波和虛反射，還可以對上行波和下行波進行波場分離，將下行波場經(jīng)過處理后疊加到上行波場的能量中可以提高信噪比。此外，由于是多方位角采樣，因此可以真正做到3D成像。與VSP相比，海洋VC具有激發(fā)點均勻、激發(fā)效率高、不受井位限制和成本相對低廉的優(yōu)點。

觀測系統(tǒng)設(shè)計在地震勘探中起著無可替代的作用，最佳觀測系統(tǒng)不僅可以提高資料品質(zhì)，還能夠降低采集成本，因此有必要對觀測系統(tǒng)進行設(shè)計、評價與優(yōu)化。覆蓋次數(shù)和照明分析是評價觀測系統(tǒng)設(shè)計優(yōu)劣的兩個極其重要的指標。本文針對單VC和多VC的特點，分別設(shè)計了相應(yīng)的觀測系統(tǒng)，然后利用覆蓋次數(shù)和照明分析對所設(shè)計的觀測系統(tǒng)進行評價，最后根據(jù)評價結(jié)果優(yōu)選出最佳設(shè)計方案。統(tǒng)計覆蓋次數(shù)時我們采用射線追蹤方法[16]來計算每個面元內(nèi)的反射點數(shù)，覆蓋次數(shù)與反射點數(shù)成正比關(guān)系?；谏渚€理論的照明分析方法[17-19]也同樣是以反射點的疏密為參考值，面元內(nèi)的反射點密集則照明能量強，反之面元內(nèi)的反射點稀疏則照明能量弱。

1 波場隨偏移距的變化特征

根據(jù)世界范圍內(nèi)已發(fā)現(xiàn)天然氣水合物區(qū)的海水深度、天然氣水合物穩(wěn)定帶厚度、似海底反射(bottom simulating reflection，BSR)埋深及地層速度，并參照我國南海西北陸坡研究區(qū)的地震層序劃分[20-23]，我們建立了含天然氣水合物的七層水平介質(zhì)3D模型(圖1)，模型參數(shù)如表1所示。該模型的海水深度為1 000 m，目標層深度為1 900 m，模型尺寸為100 000 m×100 000 m。VC的坐標為(50 000 m，50 000 m)，接收點間距為6.25 m，共161個接收點。

圖1 含天然氣水合物的七層水平介質(zhì)3D模型Fig.1 Three dimensional model constituted by seven horizontal layers containing natural gas hydrate

層號高程/m縱波速度/(m/s)橫波速度/(m/s)密度/(kg/m3)泊松比品質(zhì)因子10150010000.5001002-10001700100019860.2351003-12002200127020860.2501004-1600300010000.5001005-19002100121020480.2521006-23002700156222000.2481007-28003500202022750.250100

正演模擬方法主要有射線追蹤法[24-25]和波動方程法[26-27]。射線追蹤法具有運行速度快、運算效率高等特點，而且特別適用于模擬自由表面多次波和虛反射，因此采用射線追蹤法合成地震記錄。圖2是表1所示模型在偏移距分別為1 000 m、4 000 m、7 000 m和10 000 m時模擬一次波的合成記錄(暫不考慮臨界角的影響)，其中除直達波以外的其余6個同相軸分別對應(yīng)表1中7個地層的6個界面，地層1與地層2之間為界面1，地層2與地層3之間為界面2，以此類推。當偏移距為1 000 m時(圖2a)，各界面同相軸與地層順序一致，由于第4層是游離氣層，從圖中可以清晰地看出該層頂界面同相軸極性發(fā)生了反轉(zhuǎn)；當偏移距為4 000 m時(圖2b)，界面1、界面2、界面3與直達波同相軸相互交叉、纏繞；當偏移距為7 000 m時(圖2c)，界面4與界面6同相軸相互交叉、合并；當偏移距為10 000 m

a、b、c、d偏移距依次為1 000 m、4 000 m、7 000 m、10 000 m。圖2 VC在不同偏移距的一次波合成記錄Fig.2 Synthetic primary record of vertical cable in the case of different offset

時(圖2d)，雖然各界面同相軸相互之間不交叉，但地層順序發(fā)生了變化，從上到下依次為界面3、界面2、直達波、界面1、界面6、界面4和界面5。由此可見，隨著偏移距的增大，同相軸會出現(xiàn)交叉現(xiàn)象，地層順序也會發(fā)生變化。這是沒有模擬多次波的合成記錄，而實際中多次波非常發(fā)育，波場會更復(fù)雜難辨。因此，在設(shè)計觀測系統(tǒng)時偏移距不宜過大，以上下行波場容易識別為宜。

2 單VC激發(fā)觀測系統(tǒng)

根據(jù)單條VC的特點以及激發(fā)點組合圖案的不同，單VC激發(fā)觀測系統(tǒng)可設(shè)計為“平行形”、“正交形”、“斜交形”、“輪輻形”、“螺旋形”和“圓環(huán)形”，如圖3所示。

由圖3可見，“平行形”、“正交形”、“斜交形”和“輪輻形”激發(fā)觀測系統(tǒng)都屬于“直線型”觀測系統(tǒng)，而“螺旋形”和“圓環(huán)形”激發(fā)觀測系統(tǒng)則屬于“曲線型”觀測系統(tǒng)。在“直線型”觀測系統(tǒng)中，“平行形”、“正交形”和“斜交形”激發(fā)觀測系統(tǒng)又可進一步劃分為一類，“正交形”激發(fā)觀測系統(tǒng)可看成是由“平行形”激發(fā)觀測系統(tǒng)旋轉(zhuǎn)90°后疊加而成，而“斜交形”激發(fā)觀測系統(tǒng)則可看成是由“正交形”激發(fā)觀測系統(tǒng)旋轉(zhuǎn)45°后疊加而成。

從“平行形”到“正交形”再到“斜交形”激發(fā)觀測系統(tǒng)，面元覆蓋次數(shù)近似成倍增加，反射點空白帶逐漸減少，反射點分布越來越均勻，觀測方位角愈來愈豐富。從方位角分布來看，“曲線型”觀測系統(tǒng)比“直線型”觀測系統(tǒng)的方位角富裕。

a.平行形；b.正交形；c.斜交形；d.輪輻形；e.螺旋形；f.圓環(huán)形。圖3 單VC激發(fā)觀測系統(tǒng)示意圖Fig.3 Geometry in the case of single vertical cable

2.1 兩種改善面元覆蓋次數(shù)的效果對比

面元覆蓋次數(shù)由激發(fā)點數(shù)和接收點數(shù)共同貢獻，當單VC接收點數(shù)一定時，可以通過增加激發(fā)點數(shù)的方法以滿足覆蓋次數(shù)的要求。增加激發(fā)點的方法可以分為兩種情況：其中一種情況是當激發(fā)面積一定時，通過減小激發(fā)點間距和激發(fā)線間距，即增大激發(fā)點密度的方法來增加覆蓋次數(shù)；另外一種情況是當激發(fā)點間距和激發(fā)線間距一定時，通過增大激發(fā)面積，即增加激發(fā)線條數(shù)和激發(fā)線長度的方法來增加覆蓋次數(shù)。下面首先看第一種情況下的面元覆蓋次數(shù)分布。

采集表1的模型，第一條激發(fā)線的第一個激發(fā)點坐標為(48 000 m，48 000 m)，最后一條激發(fā)線的最后一個激發(fā)點坐標為(52 000 m，52 000 m)，激發(fā)面積為4 000 m×4 000 m。在激發(fā)面積一定的情況下，設(shè)計不同的激發(fā)點間距與激發(fā)線間距(表2)，觀察不同設(shè)計方案觀測系統(tǒng)的面元覆蓋次數(shù)(圖4)。注意：接收點數(shù)是指每條接收線上的接收點數(shù)，激發(fā)點數(shù)是指每條激發(fā)線上的激發(fā)點數(shù)(下同)。

從圖4可以看出，隨著激發(fā)點間距的縮小，反射點空白帶逐漸減少直至消失。由于海水深度有限，致使接收點數(shù)受到限制，而增加VC條數(shù)意味著高昂成本的增加；因此當激發(fā)面積一定時，只能通過加密激發(fā)點以達到滿次覆蓋的要求。從圖4還可以看出，激發(fā)線間距等于激發(fā)點間距時的覆蓋次數(shù)分布情況(圖4b、d、f)要好于激發(fā)線間距大于激發(fā)點間距時的覆蓋次數(shù)分布情況(圖4a、c、e)。

上面已經(jīng)討論了第一種情況，且其增加覆蓋次數(shù)的效果很好，下面再來分析第二種情況下的覆蓋次數(shù)分布。設(shè)計了4個方案，激發(fā)面積分別為4 000 m×4 000 m，6 000 m×6 000 m，8 000 m×8 000 m，10 000 m×10 000 m。4個方案的面元覆蓋次數(shù)分布情況如圖5所示。

從圖5可以看出，隨著激發(fā)面積的增大，觀測面積隨之增大，面元覆蓋次數(shù)有所增加，但增加的幅度甚微，當激發(fā)面積增大到10 000 m×10 000 m時，目標層仍殘存反射點空白帶。因此當觀測面積一定時，采用增加激發(fā)點密度的方法對面元覆蓋次數(shù)的改善效果要好于采用增加激發(fā)面積的方法對面元覆蓋次數(shù)的改善效果。

由表5說明，土樣1土壤中添加2%骨炭(A)化學修復(fù)劑時土壤中的重金屬鋅、鉛、鉻、銅、砷、鎘含量均有所下降，其中在種有馬鈴薯的土壤區(qū)域主要污染物鋅含量下降幅度最大，下降值為90.1mg/L。在種有油菜和馬鈴薯的土壤區(qū)域主要污染物砷含量下降幅度最大，最大值為1.0mg/L。

表2 不同空間采樣間隔觀測系統(tǒng)設(shè)計方案

a.方案一；b.方案二；c.方案三；d.方案四；e.方案五；f.方案六。圖4 不同空間采樣間隔觀測的面元覆蓋次數(shù)Fig.4 Bin fold of geometry with different spacial sampling interval

2.2 地層傾角對目標層照明的影響

當?shù)貙哟嬖趦A角時對目標層的照明會發(fā)生哪些變化呢？為了得到問題的答案，我們設(shè)計了一個目標層傾角為15°的三層水平介質(zhì)模型(圖6)。目標層深度為2 500 m，目標層方位角為90°，目標層傾斜旋轉(zhuǎn)所繞的定點坐標為(50 000 m，50 000 m)，令目標層傾角在0°、5°、10°、15°之間變化，分析不同傾角時的目標層照明情況，結(jié)果如圖7所示。

由圖7可見，當目標層傾角為0°時(圖7a)，目標層照明強度關(guān)于VC均勻、對稱分布。當目標層存在傾角時(圖7b、c、d)，目標層照明區(qū)域呈梯形，梯形的上底相對于VC位于目標層的下傾方向，梯形的下底相對于VC位于目標層的上傾方向。隨著目標層傾斜角度的增加，相對于VC位于目標層上傾方向上的觀測面積逐漸增大，相對于VC位于目標層下傾方向上的觀測面積逐漸減小，目標層照明強度分布關(guān)于VC不均勻、亦不對稱，照明強度從上傾方向到下傾方向逐級遞增，隨著傾斜角度的增加照明強度逐漸向VC位置聚焦。

a、b、c、d激發(fā)面積依次為4 000 m×4 000 m、6 000 m×6 000 m、8 000 m×8 000 m、10 000 m×10 000 m。圖5 不同激發(fā)面積觀測的面元覆蓋次數(shù)Fig.5 Bin fold of geometry with different shot area

圖6 目標層傾角為15°的三層水平介質(zhì)3D模型Fig.6 Three dimensional model constituted by three horizontal layers when the target inclination is equal to 15°

下面以傾角為15°時的目標層為例，分析觀測系統(tǒng)如何設(shè)計才能補償因存在傾角而導致的目標層照明范圍的損失。圖8a是激發(fā)面積為4 000 m×4 000 m時的目標層照明情況，激發(fā)原點坐標為(48 000.0 m，48 000.0 m)，共321條激發(fā)線，每條激發(fā)線上有321個激發(fā)點，長劃線矩形框表示傾角為0°時的目標層照明范圍。首先保持y方向上的激發(fā)線數(shù)固定不變，激發(fā)原點坐標仍為(48 000.0 m，48 000.0 m)，只增加x方向上的激發(fā)點數(shù)，當每條激發(fā)線上的激發(fā)點數(shù)增加到483個，即激發(fā)面積為6 025 m×4 000 m時，x方向上的照明范圍恰好填充到長劃線矩形框的右邊界(圖8b)。然后再增加y方向上的激發(fā)線數(shù)，激發(fā)原點坐標為(48 000.0 m，47 612.5 m)，當激發(fā)線數(shù)增加到383條，即激發(fā)面積為6 025 m×4 775 m時，y方向上的照明范圍恰好填充到長劃線矩形框的上下邊界(圖8c)。

2.3 背斜與向斜目標層照明分析

圖9a、b分別是背斜和向斜模型，目標層底部高程分別為-2 000 m和-2 500 m，長度與寬度均為2 000 m，高度均為500 m。以VC所在位置為圓心進行圓環(huán)形激發(fā)，最小環(huán)半徑為50 m，最大環(huán)半徑為2 000 m，環(huán)間距為50 m，激發(fā)點間距為25 m，每個環(huán)上的最小激發(fā)點數(shù)為18。

圖10a、b分別為該觀測系統(tǒng)對背斜與向斜目標層的照明情況，可以看出，當目標層為背斜時，只能觀測到背斜的頂部，觀測范圍為以(50 000 m，50 000 m)為圓心、半徑為150 m的圓形區(qū)域，其余一個環(huán)形區(qū)域則為觀測盲區(qū)。當目標層為向斜時，觀測區(qū)域可分為三部分，第一部分是以(50 000 m，50 000 m)為圓心、半徑為62.5 m的圓形區(qū)域，第二部分是以(50 000 m，50 000 m)為圓心、分別以587.5 m和700 m為半徑的環(huán)形區(qū)域，第三部分是以(50 000 m，50 000 m)為圓心、分別以825 m和1 000 m為半徑的環(huán)形區(qū)域，其余兩個環(huán)形區(qū)域則為觀測盲區(qū)。

圖9 背斜目標層(a)與向斜目標層(b)3D模型Fig.9 Three dimensional model constituted by anticline (a) or syncline (b)

a、b. 最大環(huán)半徑為2 000 m；c、d. 最大環(huán)半徑為3 000 m；e、f. 最大環(huán)半徑為4 000 m。圖10 背斜(a、c、e)與向斜(b、d、f)目標層照明Fig.10 Illumination of geometry when the target is anticline or syncline

當最大環(huán)半徑增大到3 000 m時(圖10c、d)，對背斜的照明半徑增加到187.5 m，對向斜的照明范圍沒有變化，增加的只是向斜邊界外的區(qū)域；當最大環(huán)半徑增大到4 000 m時(圖10e、f)，對背斜的照明半徑增加到212.5 m，觀測到了背斜邊界，對向斜的照明范圍還是沒有變化；當最大環(huán)半徑增大到5 000 m時，對背斜的照明半徑不再增加，對向斜的照明范圍仍舊沒有變化。因此，通過增大最大環(huán)半徑來增加背斜和向斜照明范圍的方法奏效甚微，即使將VC的位置任意移動也無法對背斜和向斜進行全面觀測，這時就需要在以VC坐標為圓心的圓環(huán)上布設(shè)多條VC以補償照明損失。

3 多VC接收觀測系統(tǒng)

地質(zhì)構(gòu)造變化劇烈的復(fù)雜勘探區(qū)，單條VC接收無法得到地下沉積地層的有效反射。根據(jù)VC的數(shù)量及其相互之間的位置關(guān)系，我們設(shè)計了“n階方陣”接收觀測系統(tǒng)(圖11)、“正N邊形”接收觀測系統(tǒng)(圖12)和“圓環(huán)形”接收觀測系統(tǒng)(圖13)。

所謂“n階方陣”接收觀測系統(tǒng)，就是在等間距的n行n列的交點位置布設(shè)n2條VC，每相鄰兩行兩列交點上的4條VC構(gòu)成一個接收單元，共有(n-1)2個接收單元，所有這些接收單元組合在一起我們稱之為接收區(qū)域，接收區(qū)域的面積稱為接收面積。注意：這里的接收面積指的是由外圍邊界上2n+2(n-2)條VC所圍成的自由表面面積，而并非我們通常所認為的這n2條VC實際所觀測的地下反射界面面積(觀測面積)。

采用表1的模型，觀測系統(tǒng)設(shè)計方案如表3所示，4條VC的坐標見表4。4種設(shè)計方案的接收面積分別為1 000 m×1 000 m、2 000 m×2 000 m、3 000 m×3 000 m和4 000 m×4 000 m，每種設(shè)計方案的激發(fā)面積等于接收面積，選擇“正交形”激發(fā)觀測系統(tǒng)。注意：激發(fā)線數(shù)是指y方向上的激發(fā)線數(shù)。

a. n=2；b. n=3；c. n=4。圖11 “n階方陣”接收觀測系統(tǒng)Fig.11 Geometry of n order square matrix in the case of multiple vertical cable

a. N=3；b. N=4；c. N=5；d. N=6。圖12 “正N邊形”接收觀測系統(tǒng)Fig.12 Geometry of regular polygon in the case of multiple vertical cable

a.同一圓上的接收線間距相等，不同圓上的方位角相同、接收線條數(shù)相等、接收線間距不相等；b.同一圓上的方位角不相同、接收線間距相等，不同圓上的方位角不相同、接收線間距不相等；c.同一圓上的接收線間距相等，不同圓上的接收線條數(shù)不相等、外圓比內(nèi)圓上接收線條數(shù)多；d.相鄰兩圓外圓方位角數(shù)是內(nèi)圓的兩倍，接收線條數(shù)亦是如此。圖13 “圓環(huán)形”接收觀測系統(tǒng)Fig.13 Geometry of rings in the case of multiple vertical cable

設(shè)計方案接收線數(shù)接收線間距/m接收點數(shù)接收點間距/m激發(fā)線數(shù)激發(fā)線間距/m激發(fā)點數(shù)激發(fā)點間距/m方案一410001596.2562008112.5方案二420001596.251120016112.5方案三430001596.251620024112.5方案四440001596.252120032112.5

表4 “2階方陣”接收觀測系統(tǒng)設(shè)計方案中的VC坐標

注：VC-01、VC-02、VC-03、VC-04分別代表4條不同VC的編號。

圖14為4種設(shè)計方案的面元覆蓋次數(shù)，可以看出，當激發(fā)面積等于接收面積時，觀測面積內(nèi)的面元覆蓋次數(shù)分布不均勻，4個角點附近的“矩形”區(qū)域面元覆蓋次數(shù)較高，而“十字形”中心區(qū)域的面元覆蓋次數(shù)較低。此外，隨著激發(fā)面積和接收面積的增大，4個角點區(qū)域的面積逐漸增大，中心區(qū)域的面積也會隨之增大?？梢?，當激發(fā)面積和接收面積相等時，通過同時增大激發(fā)面積和接收面積來提高中心區(qū)域面元覆蓋次數(shù)的方法行不通。

下面以接收面積為1 000 m×1 000 m的“2階方陣”接收觀測系統(tǒng)為例，保持接收面積固定不變，通過增大激發(fā)面積來觀察中心區(qū)域內(nèi)面元覆蓋次數(shù)的變化情況(圖15)。

圖15中4幅圖的激發(fā)面積分別為1 400 m×1 400 m、2 200 m×2 200 m、3 000 m×3 000 m、3 800 m×3 800 m。由圖15可見，當激發(fā)面積大于接收面積時，接收區(qū)域內(nèi)的面元覆蓋次數(shù)高于接收區(qū)域外的面元覆蓋次數(shù)，隨著激發(fā)面積的增加，觀測面積逐漸增加，接收區(qū)域內(nèi)的面元覆蓋次數(shù)逐漸升高。當激發(fā)面積增大到3 800 m×3 800 m時，反射點空白帶消失殆盡，接收區(qū)域內(nèi)的面元覆蓋次數(shù)較為均勻。

a、b、c、d接收面積依次為1 000 m×1 000 m、2 000 m×2 000 m、3 000 m×3 000 m、4 000 m×4 000 m。圖14 “2階方陣”接收觀測系統(tǒng)的面元覆蓋次數(shù)(激發(fā)面積等于接收面積)Fig.14 Bin fold of geometry of two order square matrix when shot area is equal to receiver area

a、b、c、d激發(fā)面積依次為1 400 m×1 400 m、2 200 m×2 200 m、3 000 m×3 000 m、3 800 m×3 800 m。圖15 “2階方陣”接收觀測系統(tǒng)的面元覆蓋次數(shù)(激發(fā)面積大于接收面積)Fig.15 Bin fold of geometry of two order square matrix when shot area is greater than receiver area

4 結(jié)論

1)隨著偏移距的增大，同相軸會出現(xiàn)交叉現(xiàn)象，地層順序也會發(fā)生變化。在偏移距增大過程中，有時幾個相鄰?fù)噍S之間的時差會很小，個別鄰近同相軸出現(xiàn)合并現(xiàn)象。因此在設(shè)計觀測系統(tǒng)時偏移距不宜過大，以上下行波場容易識別為宜。

2)采用增加激發(fā)點密度的方法對面元覆蓋次數(shù)的改善效果要好于采用增加激發(fā)面積的方法對面元覆蓋次數(shù)的改善效果，而且當激發(fā)線間距等于激發(fā)點間距時，每個面元內(nèi)的覆蓋次數(shù)分布均勻性最好。

3)當目標層存在傾角時，目標層照明區(qū)域呈梯形，隨著傾斜角度的增加，照明強度逐漸向VC位置聚焦。通過在構(gòu)造走向上增加激發(fā)線條數(shù)，同時在下傾方向上增加激發(fā)線長度可以補償照明損失。

4)當目標層為背斜或向斜時，通過增大最大環(huán)半徑來增加背斜和向斜照明范圍的方法奏效甚微，即使將VC的位置任意移動也無法對背斜和向斜進行全面觀測，這時就需要在以VC坐標為圓心的圓環(huán)上布設(shè)多條VC以補償照明損失。

5)當激發(fā)面積和接收面積相等時，通過同時增大激發(fā)面積和接收面積來提高中心區(qū)域面元覆蓋次數(shù)的方法行不通；而當激發(fā)面積大于接收面積時，接收區(qū)域內(nèi)的面元覆蓋次數(shù)則會高于接收區(qū)域外的面元覆蓋次數(shù)。

[1] 王淑紅, 宋海斌, 顏文. 全球與區(qū)域天然氣水合物中天然氣資源量估算[J]. 地球物理學進展, 2005, 20(4): 1145-1154.

Wang Shuhong, Song Haibin, Yan Wen. The Global and Regional Estimation of Gas Resource Quantity in Gas Hydrates[J]. Progress in Geophysics, 2005, 20(4): 1145-1154.

[2] 談順佳, 于常青, 聶逢君, 等. 青海南部天然氣水合物二維地震探測技術(shù)方法試驗研究[J]. 地球物理學進展, 2015, 30(2): 817-828.

Tan Shunjia, Yu Changqing, Nie Fengjun, et al. The Experimental Study of 2D Seismic Exploration Technology for Gas Hydrate in the Southern of Qinghai[J]. Progress in Geophysics, 2015, 30(2): 817-828.

[3] 孫嘉鑫, 寧伏龍, 鄭明明, 等. 室內(nèi)沉積物中天然氣水合物形成數(shù)值模擬研究[J]. 天然氣地球科學, 2015, 26(11): 2172-2184.

Sun Jiaxin, Ning Fulong, Zheng Mingming, et al. Numerical Simulation on Natural Gas Hydrate Formation Within Porous Media Using Constant Volume Method[J]. Natural Gas Geoscience, 2015, 26(11): 2172-2184.

[4] 唐志遠, 胡云亭, 郭清正, 等. 天然氣水合物勘探開發(fā)新技術(shù)進展[J]. 地球物理學進展, 2015, 30(2): 805-816.

Tang Zhiyuan, Hu Yunting, Guo Qingzheng, et al. Advances in Natural Gas Hydrate Exploration and Development Technology[J]. Progress in Geophysics, 2015, 30(2): 805-816.

[5] 李鵬. 海上全方位觀測系統(tǒng)的采集設(shè)計方法技術(shù)研究[D]. 武漢: 中國地質(zhì)大學(武漢), 2014.

Li Peng. Study on the Geometries Design of Marine Full-Zimuth Acquisition[D]. Wuhan: China University of Geosciences, 2014.

[6] 呂曉春. 面向目標照明與成像的海上多方位及寬方位觀測系統(tǒng)的評價與優(yōu)化[D]. 武漢: 中國地質(zhì)大學(武漢), 2014.

Lü Xiaochun. Evaluation and Optimization of Marine Wide-Azimuth and Multi-Azimuth Geometry Acquisition for Target-Oriented Illumination and Imaging[D]. Wuhan: China University of Geosciences, 2014.

[7] Krail P M. Vertical Cable as a Subsalt Imaging Tool[J]. The Leading Edge, 1994, 13(8): 885-887.

[8] Krail P M. Measurements of Cable Shape in a Cross Current[J]. Geophysics, 1994, 59(7): 1156.

[9] Rodriguez-Suarez C, Stewart R R. Survey Design for Vertical Cable Seismic Acquisition[C]//61st EAGE Conference and Exhibition International. Houten: EAGE, 1999.

[10] Moldoveanu N, Henman R, Vlasin J, et al. Bottom Referenced Vertical Hydrophone Arrays-Towed Streamers-a Comparison Study[C]//59th EAGE Conference & Exhibition. Houten: EAGE, 1997: B019.

[11]Ikelle L T. Attenuating Primaries and Free-Surface Multiples of Vertical Cable(VC) Data While Preserving Receiver Ghosts of Primaries[J]. Geophysics, 2001, 66(3): 953-963.

[12] Ikelle L T, Amundsen L, Yoo S. An Optimization of the Inverse Scattering Multiple Attenuation Method for OBS and VC Data[J]. Geophysics, 2002, 67(4): 1293-1303.

[13]Wilson R J. Potential Impacts of Vertical Cable Seismic: Modeling, Resolution and Multiple Attenuation[D]. Texas: Texas A & M University, 2002.

[14] Sun C W, Stratton J, Anderson J, et al. Separation of Up-Going and Down-Going Wave Fields of Vertical Cable Data[J]. Chinese Journal of Oceanology and Limnology, 2005, 23(3): 259-268.

[15] Leach P E. Strathspey Vertical-Cable Seismic Survey: A North-Sea First[C]//Petroleum Geology of Northwest Europe: Proceedings of the 5th Conference. London: Geological Society, 1999: 1235-1242.

[16] Zhang J Z, Huang Y Q, Song L P, et al. Fast and Accurate 3-D Ray Tracing Using Bilinear Traveltime Interpolation and the Wave Front Group Marching[J]. Geophysical Journal International, 2011, 184(3): 1327-1340.

[17]Muerdter D, Ratcliff D. Understanding Subsalt Illumination Through Ray-Trace Modeling, Part 1: Simple 2-D Salt Models[J]. The Leading Edge, 2001, 20(6): 578-594.

[18] Muerdter D, Kelly M, Ratcliff D. Understanding Subsalt Illumination Through Ray-Trace Modeling, Part 2: Dipping Salt Bodies, Salt Peaks, and Nonreciprocity of Subsalt Amplitude Response[J]. The Leading Edge, 2001, 20(7): 688-697.

[19]Muerdter D, Ratcliff D. Understanding Subsalt Illumination Through Ray-Trace Modeling, Part 3: Salt Ridges and Furrows, and the Impact of Acquisition Orientation[J]. The Leading Edge, 2001, 20(8): 803-816.

[20] 楊木壯, 王明君, 呂萬軍. 南海西北陸坡天然氣水合物成礦條件研究[M]. 北京: 氣象出版社, 2008.

Yang Muzhuang, Wang Mingjun, Lü Wanjun.The Research of Metallogenic Conditions of Natural Gas Hydrate in the Northwest Slope of the South China Sea [M]. Beijing: Meteorological Press, 2008.

[21] 肖鋼, 白玉湖, 董錦. 天然氣水合物綜論[M]. 北京: 高等教育出版社, 2012.

Xiao Gang, Bai Yuhu, Dong Jin.Review of Natural Gas Hydrate[M]. Beijing: Higher Education Press, 2012.

[22] Kvenvolden K A. Gas Hydrate-Geological Perspective and Global Change[J]. Reviews of Geophysics, 1993, 31(2): 173-187.

[23] Max M D, Lowrie A. Oceanic Methane Hydrates: A “Frontier” Gas Resource[J]. Journal of Petroleum Geology, 1996, 19(1): 41-56.

[24] 鄧飛, 劉超穎. 三維射線快速追蹤及高斯射線束正演[J]. 石油地球物理勘探, 2009, 44(2): 158-165.

Deng Fei, Liu Chaoying. 3-D Rapid Ray-Tracing and Gaussian Ray-Beam Forward Simulation[J]. Oil Geophysical Prospecting, 2009, 44(2): 158-165.

[25] 宋維琪, 楊曉東. 基于射線追蹤的微地震多波場正演模擬[J]. 地球物理學進展, 2012, 27(4): 1501-1508.

Song Weiqi, Yang Xiaodong. The Multi-Wave Field Forward Simulation of Micro-Seismic Based on Ray Tracing[J]. Progress in Geophysics, 2012, 27(4): 1501-1508.

[26] 郭宏偉, 王尚旭, 孫文博. 基于非均質(zhì)體的波動方程有限元正演模擬[J]. 石油物探, 2012, 51(4): 319-326.

Guo Hongwei, Wang Shangxu, Sun Wenbo.Wave Equation Modeling by Finite-Element Method for Heterogeneous Body[J]. Geophysical Prospecting for Petroleum, 2012, 51(4): 319-326.

[27] 李慶洋, 李振春, 黃建平, 等. 基于貼體全交錯網(wǎng)格的起伏地表正演模擬影響因素[J]. 吉林大學學報(地球科學版), 2016, 46(3): 920-929.

Li Qingyang, Li Zhenchun, Huang Jianping, et al.Factor Analysis of Seismic Modeling with Topography Based on a Fully Staggered Body-Fitted Grids[J]. Journal of Jilin University(Earth Science Edition), 2016, 46(3): 920-929.