梁 勁 王明君 陸敬安 梁金強(qiáng) 王宏斌 匡增桂

1．廣州海洋地質(zhì)調(diào)查局 2．國土資源部海底礦產(chǎn)資源實驗室 3．中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所

2007年，國土資源部廣州海洋地質(zhì)調(diào)查局在神狐海域進(jìn)行了天然氣水合物（以下簡稱水合物）鉆探取樣，得到水合物實物樣品，取得了系列重大發(fā)現(xiàn)和成果［1］。神狐海域的水合物調(diào)查在未取得測井和鉆井資料之前，主要是通過對地震剖面上似海底反射（BSR）識別和振幅空白帶的分布來判別水合物的存在以及水合物儲層厚度的變化。然而2007年的水合物鉆探巖心和測井資料表明，神狐海域的水合物僅分布在地震剖面BSR上方10～43m的范圍內(nèi)，而其上振幅空白帶相對應(yīng)的沉積層（比較均勻的粉砂質(zhì)泥）并無水合物賦存，這表明依靠BSR和振幅空白不能揭示沉積層內(nèi)水合物的賦存狀態(tài)，不能準(zhǔn)確地圈定水合物的分布面積和儲層厚度。

為此，筆者通過對神狐海域聲波測井速度與地震速度資料的精細(xì)分析，結(jié)合區(qū)域地質(zhì)條件研究陸坡區(qū)各沉積地層的速度特征，根據(jù)沉積層的速度變化來確定水合物的儲層厚度，并對聲波和地震波速度與水合物的關(guān)系進(jìn)行了分析和探討。

1 地質(zhì)概況

研究區(qū)位于南海北部陸坡中部神狐暗沙東南海域珠江口盆地珠二坳陷南翼，水深500～1 500m，水深線走向大體平行于海岸線。該區(qū)海底地形比較復(fù)雜，地形坡度變化較大，從北向南傾斜下降，上陸坡相對較陡，海底坡降平均約為3．2×10－2，下陸坡相對較緩，海底坡降平均約為2．0×10－2。該區(qū)斷裂非常發(fā)育，在陸坡上發(fā)育有海底高原、海底陡坡、海底陡崖、海谷海丘、深海海槽、沖刷槽溝和陸坡臺地等各種特殊構(gòu)造地貌單元。區(qū)內(nèi)發(fā)育的褶皺、滑塌體、增生楔等特殊構(gòu)造易于捕獲甲烷氣體，有利于水合物形成，是水合物發(fā)育的重點(diǎn)區(qū)域；各類扇體、濁積巖以及扇體底部與下部正常沉積相接觸的地方是水合物富集有利區(qū)。扇體通常具有由濱線向海進(jìn)積的特征，從而導(dǎo)致整個地質(zhì)時期中沉積相的橫向遷移，使得富集有機(jī)質(zhì)的較老的深水泥巖（生物成因氣烴源巖）與新沉積未壓實的沉積物（具有較高孔隙度的儲層）呈側(cè)變式接觸，因此生物成因氣不需要作長距離運(yùn)移，在合適的溫壓條件下就可聚集成藏形成水合物。晚中新世以來重力流比較發(fā)育，沉積速率最高達(dá)120 cm／ka，發(fā)育幾千米厚的中、新生代沉積物，累積了大量有機(jī)質(zhì)，為水合物的形成提供了物源［2－3］。

2 含水合物層聲波測井曲線、地震反射特征

2．1 含水合物層聲波測井曲線特征

地球物理測井技術(shù)是了解海底沉積層中水合物與游離氣之間關(guān)系的重要方法，正確識別水合物儲層必須要了解水合物的測井響應(yīng)特征，不同的水合物儲層在測井曲線上有不同程度的反映。含水合物層具有高速異常特征，分析聲波測井曲線是識別水合物的有效手段，可用以確定地震數(shù)據(jù)與水合物和游離氣的關(guān)系［4－5］。

2007年，我國在神狐海域?qū)嵤┝怂衔镢@探，在SH3、SH2和SH7站位成功鉆獲水合物實物樣品，筆者利用這3個站位的聲波測井曲線來分析含水合物層的特征。

2．1．1 含水合物層聲波測井曲線特征

圖1為神狐海域SH3、SH2和SH7站位先導(dǎo)孔的聲波測井曲線。由圖1可見，SH3站位聲波速度值為1 200～2 050m／s。在100～185m深度段，聲波速度緩慢增大，其中在110～130m深度段聲波速度驟然下降到1 200m／s，分析可能是該段含有氣體的緣故，但并沒有水合物存在；在185～200m深度段，聲波速度隨深度增加而起伏，其中在185～195m深度段聲波速度小幅增高，該段聲波速度平均值為1 899 m／s，并在該段采集到水合物樣品；195～200m深度段聲波速度快速下降到1 450m／s，是水合物之下有游離氣層的緣故。

圖1 神狐海域聲波測井曲線圖

SH2站位聲波測井曲線是一條典型的含水合物層的速度異常曲線（圖1），其速度范圍為1 600～2 250 m／s。在100～195m深度段，聲波速度緩慢升高，范圍在1 600～1 900m／s，是普通沉積物特征；195～215 m深度段，聲波速度隨地層深度的增加而起伏不定，速度值遞增快，之后快速遞減，整段曲線呈典型的三段式含水合物異常特征。其中在195～215m深度段聲波速度明顯增高，該厚度范圍內(nèi)聲波速度平均值為2 105m／s，并在該段成功采集到水合物樣品。215～220m深度段聲波速度快速下降到1 750m／s，是水合物之下存在游離氣層的緣故。

SH7站位聲波測井速度數(shù)據(jù)不完整（圖1）。該站位聲波測井曲線與SH－2站位相似 ，也出現(xiàn)高速、速度倒轉(zhuǎn)等含天然氣水合物沉積層的特征。在100～123 m深度段，此段曲線呈漸增式，聲波速度值為1 500～1 600m／s，此層段無明顯聲波速度變化界面；123～150m深度段聲波速度均值在1 800m／s，整體聲波速度較高；150～175m深度段速度明顯增高，最高達(dá)2 500m／s，平均為2 143m／s，并在該段成功采集到水合物樣品。

2．1．2 天然氣水合物飽和度特征

神狐海域含水合物層的水合物飽和度值是提取巖心沉積物孔隙水現(xiàn)場測試氯離子濃度，并依據(jù)其稀釋程度計算獲得的。水合物的飽和度與沉積物孔隙水的稀釋程度之間存在著一定的相關(guān)性，即自上而下氯離子濃度逐漸減小，這是巖心中水合物分解釋放出來的低鹽度自由水稀釋而形成的。依據(jù)這種方法計算的水合物飽和度，雖然只是某一井口位置垂向的飽和度數(shù)據(jù)，但可以擴(kuò)展到整個研究區(qū)域，并且體現(xiàn)了各種復(fù)雜影響因素，是沉積層水合物飽和度的真實體現(xiàn)［6－8］。

利用巖心中孔隙水測試氯離子的稀釋程度來計算水合物的飽和度首先要求確立水合物分解前原地孔隙水氯離子濃度的背景剖面，控制水合物分解的稀釋程度。假定巖心孔隙水氯離子濃度剖面上小于原地孔隙水氯離子剖面濃度的部分是水合物分解造成的，則水合物飽和度可以用如下經(jīng)驗公式進(jìn)行估算［9－11］：

式中Sh為水合物飽和度；ρh為純水合物密度，取0．9 g／cm3；Clpw為實測的巖心孔隙水中氯離子濃度；Clsw為原地孔隙水中氯離子的濃度。

圖2是神狐海域SH3、SH2和SH7站位現(xiàn)場測試的水合物飽和度散點(diǎn)圖。SH3站位含水合物層較薄，水合物集中分布在192～196m段，最高飽和度值出現(xiàn)在196m處，飽和度值為25．5%，與聲波曲線高值段相對應(yīng)，含水合物層的最高聲速值為2 051m／s。SH2站位含水合物層大約有30m厚，水合物集中分布在海底190～220m，該站位水合物飽和度相對較高，超過15%高飽和度水合物主要集中在197～220m段，飽和度值為15%～47%，平均為29%。該深度段的聲波速度值相對較高，含水合物層的最高聲速為2 279m／s。SH7站位水合物層大約有22m厚，水合物集中分布在海底155～177m。與SH2站位相似，該站位水合物飽和度較高，水合物的分布主要分為2段，155～164m為高飽和度值段，水合物的飽和度為25%～44%，平均為34%，該深度段的聲波速度值整體較高，最高聲速為2 578m／s；而166～177m深度段的水合物飽和度相對較低，飽和度值為3%～23%，平均為14%，該深度段的聲波速度值有下降的趨勢，最高聲速為2 256m／s。

承插型盤扣式鋼管支架在現(xiàn)澆箱梁中的應(yīng)用…………………………………… 韋新強(qiáng)，吳陽露，陸建培（3-159）

圖2 神狐海域水合物飽和度散點(diǎn)圖

含水合物層速度的大小與水合物的飽和度及地層的孔隙度有關(guān)。研究區(qū)的鉆探結(jié)果表明，含水合物層的巖性主要為粉砂質(zhì)泥，水合物蘊(yùn)藏于未固結(jié)的沉積物中，水合物的飽和度較高，含水合物層的速度與飽和度的對應(yīng)關(guān)系相當(dāng)明顯。由于研究區(qū)海底地質(zhì)構(gòu)造和沉積因素復(fù)雜，海底沉積層存在著較多的不確定因素，導(dǎo)致水合物飽和度隨聲波速度升高而上下波動，但總體趨勢上是隨聲波速度的升高而增高［12－13］。

2．2 含水合物層地震反射特征

2．2．1 地震剖面特征

含水合物層在地震反射剖面上常常會出現(xiàn)一強(qiáng)振幅的連續(xù)反射波，大致與海底反射波平行，被稱為似海底反射（BSR），一般代表天然氣水合物穩(wěn)定域的底界。BSR與地層產(chǎn)狀無關(guān)，當(dāng)?shù)貙赢a(chǎn)狀與海底不一致時，BSR往往與地層斜交［14－15］。另外，由于水合物濃度和厚度的增大，引起含水合物層速度增加，產(chǎn)生上拉構(gòu)造，其下方同時顯示出因低速含氣層引起的速度下拉構(gòu)造，即“眼球狀”的速度振幅異常結(jié)構(gòu)。地震剖面上通常利用BSR和“眼球狀”構(gòu)造來判別水合物層的存在［16］。含水合物層由于孔隙空間被水合物充填膠結(jié)，地層變得“均勻”，地震剖面上通常呈現(xiàn)振幅空白帶或弱振幅帶。研究表明，BSR的子波波形通常呈現(xiàn)負(fù)極性，即極性反轉(zhuǎn)，并具有較大的反射系數(shù)［17］。

圖3是神狐海域過SH2站位的地震反射剖面，從圖3中可看出，沉積物波表現(xiàn)為席狀波狀相，整體呈長條形或伸長狀席狀外形，厚度較為均一，內(nèi)部同相軸呈不對稱波狀反射結(jié)構(gòu)，向上坡方向遷移攀升；在橫切大陸坡的剖面上，呈連續(xù)的波狀或丘狀；其底部有明顯的滑脫面，這是由于底流發(fā)育時具有較強(qiáng)的侵蝕能力，對早期的沉積物進(jìn)行了簸選改造和再沉積，頻繁的侵蝕界面可能反映了底流的脈動性。剖面上TRACE900～TRACE1 000，BSR之上有一明顯的因速度振幅異常而形成的眼球狀結(jié)構(gòu)，是一條比較典型的含水合物層特征的地震剖面。另外，BSR上面有明顯的、連續(xù)性較好的振幅空白帶，極性與海底極性相反。

圖3 通過SH2站位的地震反射剖面圖

2．2．2 速度譜特征

速度是識別水合物的重要物性數(shù)據(jù)，速度異常是水合物地震響應(yīng)特征之一，是判斷沉積地層中是否含水合物的重要條件［18－22］。在神狐海域速度譜特征中，與水合物有關(guān)的速度異常一般可分為兩種類型：①水合物之下的地層含有游離氣型，其特點(diǎn)是含水合物層的速度高于上、下地層的速度，水合物之下游離氣地層的速度低于水合物之上的地層速度，且層速度下降較大；②另一種是水合物之下的地層不含游離氣型，其特點(diǎn)是與水合物位置相對應(yīng)地層的速度高于上、下地層的速度，水合物之下的地層速度高于水合物之上的地層速度，且層速度下降較小。

圖4為神狐海域含水合物層的速度譜。速度譜中紅線代表拾取的疊加速度。圖4左邊為含游離氣型譜特征，在1．90～1．98s段，疊加速度增加明顯，速度值從1 520m／s增加到1 550m／s，至1．98～2．02s段，疊加速度突然下降到1 500m／s，下降幅度較大，即水合物下地層比水合物上地層疊加速度要低。圖4右邊為不含游離氣型譜特征，在1．43～1．46s的水合物段疊加速度有一明顯的增加過程，到1．46～1．48s段，疊加速度有所下降，但下降幅度較小，即水合物下地層比水合物上地層疊加速度要高，無明顯的游離氣存在特征。

圖4 神狐海域含水合物層速度譜特征圖

2．2．3 層速度剖面特征

圖5 神狐海域130測線層速度剖面圖

圖5為神狐海域130測線用Dix公式計算的層速度剖面。在距海底大約250ms上方附近（地震剖面上解釋為BSR位置），明顯顯示出一段近似平行于海底的相對高速層，高速層上、下都是明顯的低速層。高速層的層速度在1 800～2 500m／s，其上面低速層的層速度為1 500～1 700m／s，而下面的低速層的層速度在1 100～1 500m／s，有明顯的游離氣存在的低速特征，可以認(rèn)為該測線的相對高速層速度是水合物層的反映。

2．2．4 速度反演剖面特征

對地震速度剖面進(jìn)行綜合分析，是水合物地震資料解釋的關(guān)鍵技術(shù)。高精度的速度分析可幫助判定水合物儲層及其厚度和游離氣的分布情況。通常利用測井、鉆井及地質(zhì)資料，可以較準(zhǔn)確地獲得井孔及其周圍小范圍內(nèi)某一層段的巖性及含水合物信息，但卻難以描述井間以及井孔周圍的儲層、水合物情況。而高分辨率地震資料在一定程度上可以反映儲層物性信息，在橫向上具有鉆井資料無可比擬的連續(xù)性。因此通過高精度的反演對地震資料進(jìn)行巖性識別和物性預(yù)測，是目前行之有效的手段之一［25－27］。

圖6是神狐海域HS622A測線速度反演剖面。距海底大約200ms的沉積層中有一與海底平行高速帶，厚度約為25m，速度值在2 050m／s左右，是沉積層中有水合物存在的體現(xiàn)。高速層之下有一速度值約為1 800m／s的低速帶，是水合物層之下存在游離氣層的緣故。與層速度剖面相比，有井約束速度反演剖面具有更高的分辨率，能更準(zhǔn)確地判定水合物的富集層位，礦體厚度和游離氣的分布情況，對水合物的速度精細(xì)解釋起到關(guān)鍵的作用。

圖6 神狐海域HS622A測線有井約束速度反演剖面圖

運(yùn)用多種特殊地震成像技術(shù)進(jìn)行綜合分析，是水合物地震資料解釋的有效手段。與層速度剖面相比，反演速度剖面的分辨率有了很大的提高，能夠精細(xì)刻畫地層內(nèi)部速度變化。有井約束速度反演采用聲波測井速度、伽馬和電阻率等測井曲線來做合成記錄，既克服了反演過程中可能出現(xiàn)的多解性，又補(bǔ)償了地震記錄中缺失的高頻與低頻成分，剖面分辨率比常規(guī)速度剖面和其他層速度剖面有很大的提高，能完全滿足高分辨地震剖面的分辨率（采樣率為1ms）需要。

3 討論

3．1 神狐海域BSR上方振幅空白帶與水合物的關(guān)系

當(dāng)沉積層中含有水合物時，由于水合物地震波速較高，導(dǎo)致該地層與下伏地層間的反射系數(shù)增大，地震剖面上就會出現(xiàn)相應(yīng)的強(qiáng)反射界面，含水合物層由于孔隙空間被水合物充填膠結(jié)，地層變得“均勻”，地震剖面上通常呈現(xiàn)振幅空白帶或弱振幅帶。一般情況下BSR之上的振幅空白帶是水合物存在的地震響應(yīng)，是水合探測物的重要地震證據(jù)。同時，由于振幅減弱程度與水合物含量直接相關(guān)，因此可以利用地震反射振幅信息，間接地估計水合物的含量及儲量大?。?8］。然而2007年水合物鉆探的巖心和測井資料表明，神狐海區(qū)的地震剖面上弱振幅帶和振幅空白帶相對應(yīng)的沉積層為比較均勻的粉砂質(zhì)泥，水合物儲層只有11～25m厚，大大低于預(yù)測的水合物儲層厚度。這表明振幅空白帶與其他地質(zhì)因素也有一定的關(guān)系，如沉積環(huán)境、沉積物巖性等，神狐海域水合物的沉積層只是位于BSR之上較淺的深度范圍內(nèi)，這可能是因為形成水合物的甲烷氣體很可能來源于原位微生物成因甲烷，分散狀地賦存于沉積物深部［29］。所以要判斷地震剖面上出現(xiàn)的振幅空白帶是否與水合物有關(guān)，必須結(jié)合地層速度特征進(jìn)行綜合分析。

3．2 神狐海域地層速度與水合物飽和度的關(guān)系

神狐海域水合物賦存于海底100～200m的沉積層中。純水合物的聲波速度在3 600m／s左右［30］，理論上水合物的飽和度與地層速度呈線性的正比例關(guān)系，水合物的飽和度越高，地層速度越高。但由于海底地層的構(gòu)造和沉積因素比較復(fù)雜，水合物飽和度的理論計算，是假設(shè)沉積物均勻變化，忽略了復(fù)雜地質(zhì)因素對水合物的影響，是一種理想化的計算結(jié)果，誤差的存在是不可避免的。而現(xiàn)場測試的水合物飽和度，是利用巖心中孔隙水測試氯離子的淡化程度來計算，是沉積地層中水合物飽和度的真實體現(xiàn)。在海底復(fù)雜的地質(zhì)因素作用下，水合物的飽和度與聲波速度并不是單一的比例關(guān)系，而是飽和度隨聲波速度的升高而上下波動，波動幅度在20%以內(nèi)，但總體趨勢上隨聲波速度的升高而升高，并集中分布在理論曲線附近［31］。這表明水合物飽和度是影響地層速度的重要因素。

4 結(jié)論

1）地震反射剖面上，由于水合物飽和度、厚度增大，引起含水合物層速度增大產(chǎn)生上拉構(gòu)造，其下方同時顯示出因低速含氣層引起的速度下拉構(gòu)造，即“眼球狀”的速度振幅異常結(jié)構(gòu)。

2）含水合物層速度的大小與沉積物孔隙度和水合物飽和度密切相關(guān)，水合物飽和度隨聲波速度升高而上下波動，總體趨勢上隨聲波速度的升高而增高。

3）利用高精度的速度反演對地震資料進(jìn)行巖性識別和物性預(yù)測，可得到具有更高分辨率的速度剖面，能更準(zhǔn)確地判定水合物儲層及其厚度和游離氣的分布情況，對水合物的速度精細(xì)解釋起到關(guān)鍵作用。

4）水合物層內(nèi)部高速層呈平行于海底的帶狀分布，底部速度最高，由底部往上速度逐漸降低。若水合物層下伏地層有明顯的低速特征，則水合物層下含有游離氣體。

［1］沙志彬，王宏斌，楊木壯，等．天然氣水合物成礦帶的識別技術(shù)［J］．現(xiàn)代地質(zhì)，2008，22（3）：447－451．SHA Zhibin，WANG Hongbin，YANG Muzhuang，et al．Study on recognizing technology of gas hydrates zone［J］．Geoscience，2008，22（3）：447－451．

［2］王宏斌，梁勁，龔躍華，等．基于天然氣水合物地震數(shù)據(jù)計算南海北部陸坡海底熱流［J］．現(xiàn)代地質(zhì)，2005，19（1）：67－73．WANG Hongbin，LIANG Jin，GONG Yuehua，et al．Estimation of the heat flow in the northern of the South China Sea based on the seismic data of gas hydrate［J］．Geoscience，2005，19（1）：67－73．

［3］梁勁，王宏斌，沙志彬．剩余層速度分析在南海天然氣水合物解釋中的指示意義［J］．南海地質(zhì)研究，2008：68－77．LIANG Jin，WANG Hongbin，SHA Zhibin．The signification of remaining interval velocity analysis in interpretation of gas hydrates in South China Sea［J］．Geological Research South China Sea，2008：68－77．

［4］陸敬安，楊勝雄，吳能友，等．南海神狐海域天然氣水合物地球物理測井評價［J］．現(xiàn)代地質(zhì)，2008，22（3）：447－451．LU Jing＇an，YANG Shengxiong，WU Nengyou，et al．Well logging evaluation of gas hydrates in Shenhu area，South China Sea［J］．Geoscience，2008，22（3）：447－451．

［5］陸敬安，梁金強(qiáng)．分散型天然氣水合物對儲層的物性影響特征分析［J］．南海地質(zhì)研究，2008：54－57．LU Jing＇an，LIANG Jinqiang．The impact analysis of disseminated gas hydrates on reservoir properties［J］．Geological Research South of China Sea，2008：54－57．

［6］胡高偉，業(yè)渝光，張劍，等．沉積物中天然氣水合物微觀分布模式及其聲學(xué)響應(yīng)特征［J］．天然氣工業(yè)，2010，30（3）：120－124．HU Gaowei，YE Yuguang，ZHANG Jian，et al．Micro－models of gas hydrate and their impact on the acoustic properties of the host sediments［J］．Natural Gas Industry，2010，30（3）：120－124．

［7］蘇丕波，雷懷彥，梁金強(qiáng)，等．神狐海域氣源特征及對天然氣水合物成藏的指示意義［J］．天然氣工業(yè)，2010，30（10）：103－108．SU Pibo，LEI Huaiyan，LIANG Jinqiang，et al．Characteristics of gas source in the waters of Shenhu and their significance to gas hydrate accumulation［J］．Natural Gas Industry，2010，30 （10）：103－108．

［8］王祝文，李周波，劉菁華．天然氣水合物評價的測井響應(yīng)特征［J］．物探與化探，2003，27（1）：13－17．WANG Zhuwen，LI Zhoubo，LIU Jinghua．Logging response characteristics in gas hydrates evaluation［J］．Geophysical ＆ Geochemical Exploration，2003，27（1）：13－17．

［9］甘華陽，王家生，陳建文，等．海底天然氣水合物的飽和度預(yù)測技術(shù)［J］．海相油氣地質(zhì)，2004，9（1）：111－115．GAN Huayang，WANG Jiasheng，CHEN Jianwen，et al．Methods of predicting saturation of natural gas hydrate within submarine sediments［J］．Marine Origin Petroleum Geology，2004，9（1）：111－115．

［10］王祝文，李周波，劉菁華．天然氣水合物的測井識別和評價［J］．海洋地質(zhì)和第四紀(jì)地質(zhì)，2003，23（2）：97－102．WANG Zhuwen，LI Zhoubo，LIU Jinghua．Logging identification and evaluation methods for gas hydrates［J］．Marine Geology ＆ Quaternary Geology，2003，23（2）：97－102．

［11］王秀娟，吳時國，郭旋，等．南海陸坡天然氣水合物飽和度估計［J］．海洋地質(zhì)和第四紀(jì)地質(zhì)，2005，25（3）：89－95．WANG Xiujuan，WU Shiguo，GUO Xuan，et al．Estimation of gas hydrates saturation in the continental slope，South China Sea［J］．Marine Geology ＆ Quaternary Geology，2005，25（3）：89－95．

［12］王秀娟，吳時國，劉學(xué)偉．天然氣水合物和游離氣估算的影響因素［J］．地球物理學(xué)報，2006，49（2）：505－511．WANG Xiujuan，WU Shiguo，LIU Xuewei．Factors affecting the estimation of gas hydrate and free gas saturation［J］．Chinese Journal of Geophysics，2006，49（2）：505－511．

［13］梁勁，王明君，王宏斌，等．南海神狐海域天然氣水合物聲波測井速度與飽和度關(guān)系分析［J］．現(xiàn)代地質(zhì)，2009，23（2）：217－223．LIANG Jin，WANG Mingjun，WANG Hongbin，et al．Relationship between the sonic logging velocity and saturation for gas hydrate in Shenhu area，northern slope of South China Sea［J］．Geoscience，2009，23（2）：217－223．

［14］梁勁，王宏斌，趙志超．射線追蹤法在南海天然氣水合物速度分析中的應(yīng)用［J］．物探化探計算技術(shù)，2007，29（6）：486－491．LIANG Jin，WANG Hongbin，ZHAO Zhichao．The application of ray－tracing method in velocity analysis of gas hydrate of the South China Sea［J］．Computing Techniques for Geophysical and Geochemical Exploration，2007，29（6）：486－491．

［15］劉學(xué)偉，李敏峰，張聿文，等．天然氣水合物地震響應(yīng)研究——中國南海 HD152測線應(yīng)用實例［J］．現(xiàn)代地質(zhì)，2005，19（1）：33－38．LIU Xuewei，LI Minfeng，ZHANG Yuwen，et al．Studies of seismic characteristics about gas hydrate：A case study of line HD152in the South China Sea［J］．Geoscience，2005，19（1）：33－38．

［16］徐華寧，梁蓓雯，張光學(xué)，等．南海北部陸坡天然氣水合物地震識別研究［J］．天然氣工業(yè)，2006，26（9）：49－51．XU Huaning，LIANG Beiwen，ZHANG Guangxue，et al．Seismic identification of gas hydrates on Northern slope of South China Sea［J］．Natural Gas Industry，2006，26（9）：49－51．

［17］SHIPEY T H，HOUSTON，BULLER R T，et al．Seismic evidence for wide－spread possible gas hydrates horizons on continental slopes and rises［J］．AAPG Bulletin，1979，63（12）：2204－2213．

［18］吳大奎，張本全，戴勇，等．地震資料氣水識別方法及其應(yīng)用［J］．天然氣工業(yè)，2011，31（12）：54－58．WU Dakui，ZHANG Benquan，DAI Yong，et al．Seismic identification methods for gas and water layers and their application［J］．Natural Gas Industry，2011，31（12）：54－58．

［19］ECKER C，DVORKIN J，NUR A．Sediments with gas hydrates：Internal structure from seismic AVO［J］．Geophys，1998，63（5）：1659－1669．

［20］張聿文，劉學(xué)偉，金玉潔．含天然氣水合物地層的速度和衰減研究［J］．石油地球物理勘探，2004，39（2）：205－214．ZHANG Yuwen，LIU Xuewei，JIN Yujie．Study of velocity and attenuation for gas－bearing hydrate formation［J］．Oil Geophysical Prospecting，2004，39（2）：205－214．

［21］祝有海，趙省民，盧振權(quán)．中國凍土區(qū)天然氣水合物的找礦選區(qū)及其資源潛力［J］．天然氣工業(yè)，2011，31（1）：13－19．ZHU Youhai，ZHAO Shengmin，LU Zhenquan．Resource potential and reservoir distribution of natural gas hydrate in permafrost areas of China［J］．Natural Gas Industry，2011，31（1）：13－19．

［22］傅寧，林青，劉英麗，等．白云凹陷天然氣水合物氣源條件再認(rèn)識［J］．中國石油勘探，2011，16（4）：48－53．FU Ning，LIN Qing，LIU Yingli，et al．A rethinking upon gas source condition of gas hydrate in Baiyun Sag［J］．China Petroleum Exploration，2011，16（4）：48－53．

［23］DVORKIN J，NUR A，YIN H．Effective properties of cemented granular material［J］．Mechanics of Materials，1994，18（1）：351－366．

［24］宋海斌．天然氣水合物的地球物理研究［M］．北京：海洋出版社，2003．SONG Haibin．Geophysical researches on gas hydrates［M］．Beijing：Ocean Press，2003．

［25］梁勁，王宏斌，梁金強(qiáng)．Jason反演技術(shù)在天然氣水合物速度分析中的應(yīng)用［J］．南海地質(zhì)研究，2006：114－120．LIANG Jin，WANG Hongbin，LIANG Jinqiang．The Application of Jason inversion technology in velocity analysis of gas hydrate［J］．Geological Research South China Sea，2006：114－120．

［26］楊睿，吳能友，雷新華，等．波阻抗反演在南海北部神狐海域天然氣水合物勘探中的應(yīng)用［J］．現(xiàn)代地質(zhì)，2010，24（3）：495－500．YANG Rui，WU Nengyou，LEI Xinhua，et al．Impedance Inversion and its application in gas hydrate exploration in Shenhu area，Northern South China Sea［J］．Geoscience，2010，24（3）：495－500．

［27］肖昆，鄒長春，邱禮泉，等．漠河凍土區(qū)天然氣水合物科學(xué)鉆探MK－2孔地層巖性的測井識別［J］．天然氣工業(yè)，2013，33（5）：46－50．XIAO Kun，ZOU Changchun，QIU Liquan，et al．Welllog lithology identification in well MK－2for scientific drilling and exploration of gas hydrate in Mohe permafrost，China［J］．Natural Gas Industry，2013，33（5）：46－50．

［28］LEE M W，HUTCHINSON D R，AGENA W F，et al．Seismic character of gas hydrates on the southeastern U．S．continental margin［J］．Marine Geophysical Researches，1994，16（3）：163－184．

［29］吳能友，張海啟，楊勝雄，等．南海神狐海域天然氣水合物成藏系統(tǒng)初探［J］．天然氣工業(yè)，2007，27（9）：1－6．WU Nengyou，ZHANG Haiqi，YANG Shengxiong，et al．Preliminary discussion on natural gas hydrate（NGH）reservoir system of Shenhu area，northern slope of South China Sea［J］．Natural Gas Industry，2007，27（9）：1－6．

［30］MATHEWS M．Logging Characteristics of Methane Hydrate［J］．The Log Analyst，1986，27（3）：26－63．

［31］梁勁，王明君，陸敬安，等．南海神狐海域含水合物地層測井響應(yīng)特征［J］．現(xiàn)代地質(zhì)，2010，24（3）：506－514．LIANG Jin，WANG Mingjun，LU Jing′an，et al．Logging response characteristics of gas hydrae formation in Shenhu area of the South China Sea［J］．Geoscience，2010，24（3）：506－514．