馬文嘉

摘要：天然氣水合物是一種新型非常規(guī)能源。準(zhǔn)確評估資源量，有助于制定后續(xù)開發(fā)利用方式，而選擇儲層參數(shù)與評價(jià)方法對于評估資源量具有決定性作用。概述了儲層評估的方法，包括地球物理法、地球化學(xué)法和勘探取樣，并結(jié)合各類勘探試采配套裝置詳細(xì)介紹了地震數(shù)據(jù)、隨鉆測井法、海洋可控源電磁法、氯離子濃度法、生物地球化學(xué)法和保真取心配套分析技術(shù)。將天然氣水合物的資源量評估方法分成了體積法、面積法和碳平衡法3種，通過對近20年國內(nèi)外天然氣水合物資源量估算結(jié)果的比較，剔除偏差較大的估值，得出了目前全球天然氣水合物資源量的合理估值，數(shù)量級分別為海底陸坡10 15 m3、凍土層10 13 m3。其中中國天然氣水合物資源量的合理估值，數(shù)量級分別為：南海為10 12～10 13 m3，東海沖繩海槽為10 12 m3，青藏高原凍土層為10 12～10 13 m3，漠河為10 12～10 13 m3。

關(guān)鍵詞：天然氣水合物;儲層評估;地球物理;地球化學(xué);勘探試采裝置;資源量估算

中圖分類號：TH455

文獻(xiàn)標(biāo)志碼：A

文章編號：1009-9492f 2022）02-0021-07

0 引言

天然氣水合物（ Natural Cas Hydrate，NCH）是水和甲烷氣體形成的非化學(xué)計(jì)量性籠狀晶體物質(zhì)[1]。它的形成不僅需要時(shí)間和空間，還需要?dú)庠?、水以及低溫和高壓的條件，因而主要在凍土層以下和海底陸坡生成[2]。其在自然界中大量存在且分布廣泛，已在全球的79個(gè)國家超過230個(gè)區(qū)域發(fā)現(xiàn)天然氣水合物[3]。作為一種新型非常規(guī)能源，天然氣水合物有望成為傳統(tǒng)能源的替代品[4]。對于水合物資源量估算，國內(nèi)外做了大量研究，有學(xué)者認(rèn)為天然氣水合物所含天然氣資源量介于3.lx10 12 m3與7.6x10 15 m3之間[5]，也有學(xué)者認(rèn)為介于10 15 m3與10 18 m3之間[6]，均認(rèn)可其巨大資源潛力。但由于估算方法有差異，評估時(shí)的儲層參數(shù)不同，且受到勘探與試采結(jié)果的影響，資源量的估算存在較大差異。因此未來業(yè)內(nèi)仍需加大對儲層性質(zhì)與資源量的研究，注重勘探與試采新裝置的研發(fā)。

本文主要綜述了兩個(gè)方面內(nèi)容：基于地球物理法、地球化學(xué)法，運(yùn)用水合物識別和儲層評估技術(shù)，借助無人遙控潛水器、海洋可控源電磁探測系統(tǒng)、近海底原位多參量地球化學(xué)測量裝置、隨鉆測井裝置、鉆探取樣系統(tǒng)等手段，獲得儲層參數(shù);基于天然氣水合物的資源評估方法體積法、面積法和碳平衡法，結(jié)合水合物儲層參數(shù)對資源量的影響，對近20年內(nèi)世界上其他國家以及中國天然氣水合物資源量的估值進(jìn)行了對比分析。

1 天然氣水合物儲層評估

天然氣水合物的儲層參數(shù)主要包括水合物的飽和度、分布深度、分布面積、離子濃度、孔隙度、滲透率等，而資源量的評估更是受到了水合物飽和度、分布深度、分布面積和孔隙度的影響。通常評價(jià)水合物儲層，可通過地球物理法、地球化學(xué)法，并結(jié)合保真取心技術(shù)。水合物相關(guān)分布位置及深度的信息，主要是通過遙感調(diào)查獲得，而隨著勘探技術(shù)的發(fā)展，無人遙控潛水器技術(shù)與裝備被用于精細(xì)作業(yè)。如圖1所示[7]，我國自主研發(fā)的“海馬”號無人遙控潛水器（ROV），本體配有液動壓力、推動螺旋槳、云臺、視頻和照明、導(dǎo)航定位、供配電、檢測與控制、多功能機(jī)械手和緊急定位裝置。

地球物理技術(shù)主要是地震反演技術(shù)和測井技術(shù)，地球化學(xué)技術(shù)主要是研究烴類組分比值、分析同位素來判斷天然氣成因，兩種方法相輔相成。地球物理和地球化學(xué)法得出的結(jié)論的準(zhǔn)確性，可通過鉆探取樣驗(yàn)證。在此基礎(chǔ)上，可通過建立數(shù)學(xué)模型對天然氣水合物的資源量做出評估。

1.1 地球物理評估

1.1.1 地震數(shù)據(jù)

地震數(shù)據(jù)可凸顯因天然氣水合物的存在而使得沉積層彈性發(fā)生的改變，用來反演判斷水合物是否存在，其中主要通過海底近似反射（BSR）、振幅空白、振幅隨偏移距的變化和速度異常[8]。Li等[9]研究了廣州海洋地質(zhì)調(diào)查局第二次勘探（GMSG2）的結(jié)果，得到BSR與水合物穩(wěn)定帶（ GHSZ），如圖2所示。BSR分布在具有雙向行程時(shí)間（TWT）的大陸斜坡上，CHSZ厚度根據(jù)BSR和海平面之間的深度差計(jì)算獲得。在BSR發(fā)生的區(qū)域內(nèi)，CHSZ的厚度介于60 -300 m，平均為200 m，并隨著水層的減少而減小[9]。

1.1.2 隨鉆測井法

隨鉆測井是一種原位測量的手段，具有初始壓力和溫度條件保真的特點(diǎn)，是水合物資源勘探過程中必不可少的方式[10]。裝置如圖3所示，相關(guān)技術(shù)目前還得依靠斯倫貝謝、哈里伯頓等國外公司[11]。

趙軍等[12]認(rèn)為在勘探天然氣水合物的過程中，可采用井徑測井、自然電位測井、自然伽馬測井、聲波測井、密度測井、中子孔隙度測井、電阻率測井、鉆井速率測井、鉆井液錄井、電成像測井、核磁共振測井、偶極聲波測井、介電測井、y射線測井和電磁波測井的方法?？捉系萚13]認(rèn)為宜采用電阻率測井、聲波測井、核磁共振測井和y射線測井來評估儲層的飽和度。Wang等[14]研究了CMCS3、CMCS4的測井?dāng)?shù)據(jù)，通過y射線測井和電阻率測井得到了儲層飽和度，如圖4所示。在圖4 （a）中可以看到，y射線測井的API值介于50-80之間，對應(yīng)沉積層為淤泥/淤泥粘土，即細(xì)粒沉積物。而在圖4 （b）中可以發(fā)現(xiàn)在126-238 mbsf的范圍內(nèi)，電阻率測井顯示的值（黑線）高于上覆和下伏剖面（紅線），最大差值為5.8 Ω'm。結(jié)合現(xiàn)場取心發(fā)現(xiàn)的水合物[15]，Wang等[14]認(rèn)定圖4（c）中電阻率發(fā)生增加的位置，是水合物出現(xiàn)的區(qū)域，而電阻增加也是由于水合物的出現(xiàn)，并根據(jù)電阻率估算出水合物的飽和度為15%-45%，平均為30.3%。

1.1.3 海洋可控源電磁法

海洋可控源電磁法（MCSEM）是一種測量地下介質(zhì)電阻率的手段，通過測量周圍空間的電磁場，利用電磁波在不同介質(zhì)中的傳輸速率，反應(yīng)沉積層中電阻率的情況[16]。該探測系統(tǒng)通常包括拖曳式電磁發(fā)射、拖曳式電場接受、海底電磁采集站、甲板監(jiān)控單元以及輔助儀器設(shè)備[17]。為了探測更深更廣的區(qū)域，電磁發(fā)射機(jī)應(yīng)具有更高的功率、更大的電流、更輕便、更高的精度和電能轉(zhuǎn)化效率[18]。

Hsu等[19]在臺灣西南利用海洋可控源電磁法進(jìn)行了研究，得到沿測量軌道的可視電阻率，如圖5所示。由圖可知，C12軌道下面電阻率高達(dá)1 Ω 'm，出現(xiàn)了異常，這部分區(qū)域被認(rèn)為存在水合物，進(jìn)而可以通過電阻率推算出其飽和度[19]。

1.2 地球化學(xué)評估

地球化學(xué)法在水合物勘探領(lǐng)域的最主要貢獻(xiàn)在于查找氣源以及天然氣水合物成藏模式的判定，例如同位素標(biāo)記。雖然暫未存在一種像BSR或者M(jìn)CSEM這樣能直接確定天然氣水合物存在的技術(shù)[20]，但是可用來預(yù)測儲層水合物的飽和度并確定水合物所在的位置，而生物地球化學(xué)方法更是通過碳在自然界中的傳遞反應(yīng)，可直接用于資源量的預(yù)測。以我國自主研發(fā)的近海底原位多參量地球化學(xué)測量拖體為例，可原位探測獲取水合物靶區(qū)地球化學(xué)方面的基礎(chǔ)數(shù)據(jù)和資料，裝置如圖6所示[21]]。

1.2.1 氯離子濃度法

氯離子濃度法是利用水合物所在孔隙與海水中的氯離子濃度具有差異，因?yàn)樘烊粴馑衔镌谛纬蓵r(shí)會使原位孔隙中的氯化物脫除，并隨之?dāng)U散開，巖心取樣期間隨著水合物的分解，氯離子濃度會進(jìn)一步降低。Yuan等[22]利用氯離子濃度和水合物密度，對海洋鉆探計(jì)劃（ODP）的水合物飽和度進(jìn)行了預(yù)測。Liu等[23]利用聲波和電阻在南海北部神狐海域測定的水合物飽和度變化，與利用氯離子濃度測定得到的水合物飽和度規(guī)律呈現(xiàn)出一致性。Chen等[24]基于徑向基函數(shù)插值的數(shù)學(xué)描述法，以南海神狐SH2的41個(gè)取心樣品以及SH7的21個(gè)取心樣品作為數(shù)據(jù)，使用氯離子濃度相關(guān)的阿爾奇公式計(jì)算水合物飽和度，得到的飽和度與利用聲波或者電阻計(jì)算的結(jié)果近似。

1.2.2 生物地球化學(xué)法

生物地球化學(xué)方法是利用地球化學(xué)以及熱力學(xué)的方法確定CHSZ，利用地球化學(xué)的方法尋找天然氣水合物的氣源，根據(jù)海底顆粒有機(jī)碳（POC）以及海底微生物降解和甲烷生成的速率，考慮孔隙流體和甲烷氣體進(jìn)入CHSZ并發(fā)生沉降來計(jì)算天然氣水合物的資源且[25]。Wallmann等[25]建立了地球化學(xué)傳遞反應(yīng)模型，假設(shè)沉降速率為POC的1 wt%，利用該模型對全球海域天然氣水合物的資源量進(jìn)行了評估。Burwicz等[26]研究了ODP的鉆探樣品以及模擬了千萬年CHSZ中沉降與溫壓變化，發(fā)現(xiàn)甲烷厭氧氧化（AOM）指示的地球化學(xué)測量結(jié)果能夠支持單一原位甲烷氣源模擬的結(jié)果，模擬生成的天然氣水合物的4-7 voL.%是由甲烷熱流貢獻(xiàn)，其余為顆粒有機(jī)碳沉降貢獻(xiàn)。

1.3 鉆探取樣

付強(qiáng)等[27]根據(jù)適用性，將海域水合物開采技術(shù)設(shè)備分成了通用關(guān)鍵技術(shù)設(shè)備和專用技術(shù)設(shè)備，前者包括海底采礦車、疏松淺表層雙梯度鉆井技術(shù)、高強(qiáng)度輕質(zhì)柔性管、舉升泵和井下多相分離器，后者包括防砂技術(shù)設(shè)備、預(yù)造斜導(dǎo)向鉆進(jìn)技術(shù)、水合物破碎回收分離一體化鉆采技術(shù)、水合開采實(shí)驗(yàn)?zāi)M平臺和“三氣合采”技術(shù)裝備。鉆探取樣是一種直接勘探研究的技術(shù)，所用的取心及配套分析設(shè)備包括配有保壓取樣器的高壓取心設(shè)備、保壓巖心分析與轉(zhuǎn)移系統(tǒng)和保壓巖心表征工具。根據(jù)鉆探區(qū)域的不同，凍土層所需設(shè)備與傳統(tǒng)油氣開采相同，而海底陸坡則需使用中深水半潛式鉆井平臺或者鉆井船[28]。在南海神狐海域，2017年試采時(shí)使用的“藍(lán)鯨I號”和2020年試采時(shí)使用的“藍(lán)鯨II號”均屬于半潛式鉆井平臺。鉆探取樣可分為保溫保壓和非保溫保壓，要獲得保真巖心樣品，就必須采用保溫保壓取心裝置，在切割、分離、轉(zhuǎn)移過程中維持溫度壓力不變，從而分析樣品的表觀體積、密度、總重量和含氣含水量來確定水合物的飽和度。

2 天然氣水合物資源量評估

學(xué)者們對天然氣水合物資源量評估的方法有很多種[29]，大致可以分為3種：體積法、面積法和碳平衡法。體積法[30]如下所示：

y= f（Z×φ×H×G×E）dA

（1）式中：V為標(biāo)況下天然氣水合物的體積;A為天然氣水合物沉積層的面積;Z為水合物穩(wěn)定帶的厚度;咖為沉積層中的孔隙度;H為孔隙中天然氣水合物的飽和度;G為每立方米水合物釋放甲烷氣體的標(biāo)況體積;E為天然氣水合物的聚集率。

面積法[31]如下式所示：

V=（A×R）shelf+（A×R）slope+（A×R）a byssalplane

（2）式中：y為標(biāo)況下天然氣水合物的體積;A為天然氣水合物沉積層的面積;R為6標(biāo)況下天然氣水合物的資源密度。

碳平衡法則是利用有機(jī)碳的沉降以及熱流的影響來模擬天然氣水合物的生成，它被發(fā)展成一個(gè)非常復(fù)雜的模型[25]，先通過溶解無機(jī)碳濃度、溶解甲烷濃度和甲烷抑制有機(jī)物降解常數(shù)得到顆粒有機(jī)碳沉降速率，再通過熱力學(xué)性質(zhì)確定海洋中的CHSZ，最后將沉積層表面的顆粒有機(jī)碳濃度、沉降速率與CHSZ相結(jié)合，得到海域天然氣水合物的資源量。

2.1 全球天然氣水合物資源量

隨著研究深入，發(fā)現(xiàn)學(xué)者們早期的研究可能存在很大局限性，表1所示為近20年全球海底陸坡、凍土層以及中國以外部分著名地區(qū)天然氣水合物資源量，這里折算成標(biāo)況下甲烷氣體的體積。在確定CHSZ上，3種評價(jià)方法沒有太大差異，均是利用水合物在海洋中的熱力學(xué)穩(wěn)定性來判斷可能生成的區(qū)域，并根據(jù)地球物理、地球化學(xué)等勘探數(shù)據(jù)進(jìn)行進(jìn)一步的篩選。早期的面積法和體積法，采用單一的儲層參數(shù)進(jìn)行預(yù)測，而碳平衡法的出現(xiàn)，開始以物質(zhì)平衡的角度、天然氣水合物氣源以及成藏原因?qū)訁?shù)進(jìn)行生成。天然氣水合物的資源量在全球海底陸坡為0.82x10 13 -774.7x10 15 m3，全球凍土層為1.87x10 13 m3。其中克里希納戈達(dá)瓦里（KG）盆地、韓國郁陵盆地、羅斯海、日本南開海槽弧前熊野盆地、東海烏倫盆地、加利福尼亞州邊境地區(qū)和亞得里亞海的天然氣資源數(shù)量在101i m3，南極陸緣、墨西哥北部灣和智利巴塔哥尼亞海岸的天然氣資源量在1013 m3。

Klauda等[36]采用了體積法的公式來模擬碳沉積的總量，采用的碳沉積速率較大因而估算總量偏大。Wood等[37]采用永凍圈中沉積物的平均甲烷濃度，使得估算總量偏大。Burwicz等[39]的估值差異是因?yàn)檫x取了不同的沉積模型，低估值對應(yīng)全新世沉積速率，高估值對應(yīng)第四紀(jì)沉積速率和較高的甲烷通量。剔除這些估算偏差較大的，剩下的估值出現(xiàn)了較好的一致性，因而目前的研究結(jié)果表明，海底陸坡天然氣水合物資源量在10is m3，高出凍土層中天然氣資源量兩個(gè)數(shù)量級。

2.2 中國天然氣水合物資源量

在中國，天然氣水合物主要分布在南海海域、東海沖繩海槽、青藏高原凍土層和東北漠河凍土層，表2所示為近20年來的中國天然氣水合物資源量。與全球海域中儲層參數(shù)取值相比，姚伯初[ss】在儲層參數(shù)上取值較大，導(dǎo)致估值較大。唐勇等[60]采用體積法時(shí)沒有考慮聚集率，導(dǎo)致估值較大。Chen等[61]和Zhao等[63]只是做了粗略估計(jì)，導(dǎo)致估值有偏差。目前公認(rèn)相對合理的資源量，是在對估值偏差較大區(qū)域的數(shù)值進(jìn)行調(diào)整以后得出的，天然氣水合物在南海數(shù)量級在10 12～10 13 m3，在東海沖繩海槽的數(shù)量級在1012 m3，在青藏高原凍土層的數(shù)量級10 12～10 13 m3，在漠河的數(shù)量級在10 12～10 13 m3。

3 結(jié)束語

在天然氣水合物儲層評價(jià)方法中，地球物理法適用于判斷是否存在水合物，包括BSR和MCSEM技術(shù)，地球化學(xué)法適用于其氣源和成藏模式的判斷。其中地震數(shù)據(jù)，隨鉆測井法，海洋可控源電磁法，氯離子濃度法和生物地球化學(xué)法均可確定天然氣水合物存在的區(qū)域并計(jì)算其飽和度，但單獨(dú)使用時(shí)都會出現(xiàn)判斷錯(cuò)誤，因而通常聯(lián)合使用，并結(jié)合勘探與采樣數(shù)據(jù)進(jìn)行校對。天然氣水合物資源評估的3種方法，可采用相同的方法來確定CHSZ，但由于碳平衡法能夠用于解釋天然氣水合物形成的原因，似乎是評估手段未來發(fā)展的趨勢。通過比較近20年國內(nèi)外天然氣水合物資源量的估算結(jié)果，得出了目前全球天然氣水合物資源量的合理估值，數(shù)量級分別為：海底陸坡為10 15 m3、凍土層為10 13 m3。其中中國天然氣水合物資源量的合理估值，數(shù)量級分別為：南海為10 12～10 13 m3、東海沖繩海槽為10 12 m3、青藏高原凍土層為10 12～10 13 m3、漠河為10 12～10 13 m3。

儲層參數(shù)因?yàn)榈赜虿煌?、氣源以及成藏模式各異出現(xiàn)很大的差別，因而在資源評估時(shí)，需要細(xì)分每個(gè)水合物礦藏區(qū)域，努力探索更加準(zhǔn)確的資源量，為進(jìn)一步制定經(jīng)濟(jì)、合理的開采方案打下堅(jiān)實(shí)的基礎(chǔ)?？梢灶A(yù)期，將來水合物勘探會不斷加強(qiáng)，試采工藝流程也會不斷完善和提高。業(yè)內(nèi)對于勘探、試采裝置不斷增加投入，加快研發(fā)與創(chuàng)新，各種更加精密的探測設(shè)備以及鉆探裝置會不斷涌現(xiàn)，必將會為水合物的勘探、試采提供更加有力的保障。同時(shí)隨著國際交流合作進(jìn)一步加深，各國之間互聯(lián)互通、數(shù)據(jù)共享，人工智能等先進(jìn)技術(shù)手段不斷應(yīng)用到實(shí)踐中去，水合物勘探與試采會迎來一個(gè)嶄新的局面，實(shí)現(xiàn)跨越式發(fā)展。

參考文獻(xiàn)

[1]Sloan E D.Fundamental principles and applications of naturalgas hydrates[J]. Nature， 2003. 426（6964）： 353-359.

[2] Kvenvolden K A.Gas hydrates-geological perspective and glob-al change[J]. Reviews of geophysics. 1993， 31（2）： 173-187.

[3] Makogon Y F.Naturalgas hydrates -A promising source of ener-gy[J]. Joumal of natural gas science and engineering， 2010， 2（1）：49-59.

[4] Kerr R A.Gas hydrate resource： smaller but sooner[J]. Science，2004. 303（5660）：946-947.

[5] Collett T S.Energy resource potential of natural gas hydrates[J].AAPG Bulletin. 2002， 86（11）：1971-1992.

[6] Reagan M T，Moridis G J，Collett T S，et al.Toward ProductionFrom Gas Hydrates： Current Status. Assessment of Resources.and Simulation-Based Evaluation of Technology and Potential[J]. Spe Reservoir Evaluation&Engineering， 2009， 12（5）：745- 771。

[7]陳宗恒，田烈余，胡波，等，”海馬”號ROV在天然氣水合物勘查中的應(yīng)用[J].海洋技術(shù)學(xué)報(bào)，2018，37（2）：6.

[8]景鵬飛，胡高偉，卜慶濤，等.天然氣水合物地球物理勘探技術(shù)的應(yīng)用及發(fā)展[J].地球物理學(xué)進(jìn)展，2019，34（5）：19.

[9] Li L，Liu H，Zhang X，et al.BSRs， estimated heat flow. hydrate-relatedgas volume and their implications for methane seepageand gas hydrate in the Dongsha region， northem South China Sea[J]. Marine& Petroleum Geology， 2015（67）：785-794.

[10]高興軍，于興河，李勝利，等，地球物理測井在天然氣水合物勘探中的應(yīng)用[J].地球科學(xué)進(jìn)展，2003，18（2）： 305-311.

[11]王麗忱，李男.國內(nèi)外天然氣水合物測井方法應(yīng)用現(xiàn)狀及啟示[J].中外能源，2015（4）：7.

[12]趙軍，武延亮，周燦燦，等.天然氣水合物的測井評價(jià)方法綜述[J].測井技術(shù)，2016 （4）： 392-398.

[13]孫建孟，羅紅，焦滔，等.天然氣水合物儲層參數(shù)測井評價(jià)綜述[J].地球物理學(xué)進(jìn)展，2018， 33（2）： 715-723.

[14] Wang J， Wu S，Kong X， et al. Geophysical characterization of afine-grained gas hydrate reservoir in the Shenhu area. northernSouth China Sea： Integration of seismic data and downhole logs[J]. Marine and Petroleum Geology， 2018（92）： 895-903.

[15] Yang S，Liang J， Lei Y. et al. GMGS4 gas hydrate drilling expe-dition in the South China Sea[J]. Fire in the Ice. Methane Hy-drate Newsletter. National Energy Technology Laboratory. USDepartment of Energy， 2017， 17（1）： 7-11.

[16]陳杰，胡高偉，卜慶濤.BSR與CSEM識別天然氣水合物的優(yōu)缺點(diǎn)對比[J].新能源進(jìn)展，2018， 6（5）：112-124.

[17]趙鐵虎，齊君，胡剛，等，天然氣水合物勘查相關(guān)技術(shù)集成研發(fā)進(jìn)展[J].海洋地質(zhì)前沿，2015， 31（6）：7.

[18]楊蜀冀，胡楊承鈺，宋紅喜，等.海洋可控源電磁探測系統(tǒng)概述[J].電氣時(shí)代，2018（4）：3.

[19] Hsu S K，Chiang C W， Evans R L，et al. Marine controlledsource electromagnetic method used for the gas hydrate investi-gation in the offshore area of SW Taiwan[J]. Journal of AsianEarth Sciences. 2014（92）： 224-232.

[20]張富貴，孫忠軍，楊志斌，等，中國凍土區(qū)天然氣水合物的地球化學(xué)調(diào)查[J].礦物巖石地球化學(xué)通報(bào)，2019， 38（6）：11.

[21]單晨晨，國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目”近海底高精度水合物探測技木”完成首次深海聯(lián)合測試[EB/OL]. （2019-06-07）. ttps：//www.cgs.gov.c n/gzdt/z sdw/201906/t20190626_485 25 5.html.

[22] Yuan T，Hyndman R D，Spence G D，et al. Seismic velocity in-crease and deep-sea gas hydrate concentration above a bottom-simulating reflector on the northern Cascadia continental slope[J]. Journal of Geophysical Research： Solid Earth. 1996. 101（B6）： 13655-13671.

[23] Liu J， Zhang J， Ma F，et al. Estimation of seismic velocities andgas hydrate concentrations：A case study from the Shenhu area.northem South China Sea[J]. Marine and Petroleum Geology，2017（88）： 225-234.

[24] Chen Y，Dunn K J，Liu X， et al. New method for estimating gashydrate saturation in the Shenhu area[J]. Geophysics. 2014， 79（5）： IM11-IM22.

[25] Wallmann K，Pinero E，Burwicz E，et al. The Global Inventoryof Methane Hydrate in Marine Sediments：A Theoretical Ap-proach[J]. Energies， 2012， 5（7）：2449-2498.

[26] urwicz E，Rupke L.Thermal State of the Blake Ridge Gas Hy-drate Stability Zone （GHSZ）-Insights on Gas Hydrate Dynam-ics from a New Multi-Phase Numerical Model[J]. Energies，2019， 12（17）： 3403.

[27]付強(qiáng)，王國榮，周守為，等.海洋天然氣水合物開采技術(shù)與裝備發(fā)展研究[J].中國工程科學(xué)，2020， 22（6）：8.

[28]張輝，金偉晨，趙羿羽，等，天然氣水合物開發(fā)裝備技術(shù)體系研究與展望[C]//北京造船工程學(xué)會2018-2019年學(xué)術(shù)論文集， 2020.

[29]叢曉榮，吳能友，蘇明，等.天然氣水合物資源量估算研究進(jìn)展及展望[J].新能源進(jìn)展，2014， 2（6）： 462-470.

[30] Trofimuk A A，Cherskiy N V，Tsarev V P.Accumulation of natu-ral gases in hydrate formation zones of the world's oceans[J].DokladyAkademiinauk SSSR. 1973.

[31] Trofimuk A A，Cherskiy N V，Tsarev V P.The biogenic meth-ane resources in the oceans[J]. DokladyAkademiinauk SSSR.1975（225）：936-939.

[32] Soloviev V A.Global estimation of gas content in submarine gashydrate accumulations[J]. Geologiya I Geofizika. 2002， 43（7）：648-661.

[33] Milkov A V，Claypool G E，Lee Y J，et al.In situ methane con-centrations at Hydrate Ridge， offshore Oregon： New constraintson the global gas hydrate inventory from an active margin[J]. Ge-ology， 2003， 31（10）：804-6.

[34] Milkov A V.Global estimates of hydrate-bound gas in marinesediments： How much is really out there？[J]. Earth-Science Re-views， 2004. 66（3）：183-197.

[35] Buffett B，Archer D.Global inventory of methane clathrate： sen-sitivity to changes in the deep ocean[J]. Earth and Planetary Sci-ence Letters. 2004. 227（3-4）：1 85-199.

[36] Klauda J B，Sandler S I.Global distribution of methane hydratein ocean sediment[J]. Energy& Fuels. 2005， 19（2）： 459-470.

[37] Wood W T，Jung W Y.Modeling the extent of Earth's marinemethane hydrate cryosphere[J]. 2008.

[38] Archer D，Buffett B，Brovkin V.Ocean methane hydrates as aslow tipping point in the global carbon cycle[J]. Proceedings ofthe National Academy of Sciences， 2009， 106（49）： 20596-20601.

[39] Burwicz E B，Rupke L H， Wallmann K.Estimation of the globalamount of submarine gas hydrates formed via microbial meth-ane formation based on numerical reaction-transport modelingand a novel parameterization of Holocene sedimentation[J]. Geo-chimica et Cosmochimica Acta. 2011. 75（16）： 4562-4576.

[40] Boswell R，Collett T S.Current perspectives on gas hydrate re-sources[J]. Energy&environmental science. 2011（4）： 1206-1215.

[41] Pinero E，Marquardt M.Hensen C.et al.Estimation of the globalinventory of methane hydrates in manne sediments using trans-fer functions [J]. Biogeosciences， 2013. 10（2）： 959-975.

[42] Ruppel C. Permafrost-associated gas hydrate： is it really ap-proximately lo-/o of the global system？[J]. Journal of Chemical &Engineering Data， 2015， 60（2）： 429-436.

[43] Kretschmer K. Biastoch A， RupkeL， et al. Modeling the fate ofmethane hydrates under global warming [J]. Global Biogeochem-ical Cycles， 2015， 29（5）：610-625.

[44] Shankar U. Riedel M. Gas hydrate saturation in the Krishna -Godavari basin from P-wave velocity and electrical resistivitylogs[J]. Marine & Petroleum Geology， 2011， 28（10）：1768-1778.

[45] Lee G H. BoYY. Yoo D G， et al. Estimation of the gas-hydrateresource volume in a small area of the Ulleung Basin， East Seausing seismic inversion and multi-attribute transform tech- niques [J]. Marine and petroleum geology， 2013（47）： 291-302.

[46]王力峰，鄧希光 ，沙志彬，等 .南極陸緣熱流分布與天然氣水合物資源量研究[J].極地研究， 2013， 25（3）：241-248.

[47]王威 ，高金耀 ，沈中延，等 .羅斯海天然氣水合物成藏條件及資源量評估[J]. 海洋學(xué)研究， 2015， 33（1）： 16-24.

[48] Taladay K， Boston B. Moore G F. Gas-in-place estimate for po-tential gas hydrate concentrated zone in the Kumano Basin.Nankai Trough Forearc， Japan[J]. Energies， 2017（10）： 1552.

[49] Burwicz E B. Reichel T. Wallmann K. et al. 3-D basin-scale re-construction of natural gas hydrate system of the Green Canyon，Gulf of Mexico[J]. Geochemistry， Geophysics， Geosystems，2017. 18（5）：1959-1985.

[50] Majumdar U. Cook A E. The Volume of Gas Hydrate-BoundGas in the Northern Gulf of Mexico[J]. Ceochemistry， Ceophys-ics. Geosystems， 2018， 19（11）： 4313-4328.

[51] Villar-Munoz L. Bento J P. Klaeschen D. et al. A first estima-tion of gas hydrates offshore Patagonia （Chile）[J]. Marine and Pe-troleum Geology， 201 8（96）： 232-239.

[52] Bo Y Y， Lee G H， Kang N K. et al. Deterministic estimation ofgas-hydrate resource volume in a small area of the Ulleung Ba-sin. East Sea （Japan Sea） from rock physics modeling and pre-stack inversion[J]. Marine and Petroleum Geology， 2018（92）：597-608.

[53] Kannberg P K， Constable S. Characterization and quantificationof gas hydrates in the California Borderlands[J]. GeophysicalResearch Letters， 2020， 47（6）： 1-8.

[54] Obhodas J，Tinivella U.Giustiniani M，et al.Past and presentpotential of the Adriatic deep sea sediments to pIDduce methanehydrates[J].Joumal of Soils and Sediments，2020，20（6）：2724—2732.

[55]姚伯初，南海的天然氣水合物礦藏[J].熱帶海洋學(xué)報(bào)，2001，20（21：20—28.

[56]Guan J，Liang D，Wu N，et al.The methane hydrate formationand the resource estimate resulting fbm free gas migration inseeping seanoor hydrate stability zone[J].Journal of Asian EarthSciences，2009，36（4—5）：277—288.

[57]Tnlng N N.The gas hydrate potential in the South China Sea[J].Joumal of petroleum science and engineering，2012（88）：41—47.

[58]Wa“g J，Wu S，Yao Y.Quantifyi“g gas hydrate from microbialmethane in the South China Sea[J].Joumal of Asian earth sci—ences，2018（168）：48—56.

[59]Li L，Lei X，Zhang X，et al.Heat flow derived from BSR and itsimplications for gas hydIate stability zone in Shenhu Area ofnorthem South China Sea[J].Marine Geophysical Research，2012，33（1）：77—87.

[60]唐勇，方銀霞，高金耀，等，沖繩海槽天然氣水合物穩(wěn)定帶特征及資源量評價(jià)[J]，海洋地質(zhì)與第四紀(jì)地質(zhì)，2005，25（4）：79—84.

[6l]CHEN D F，WANG M C，xia B.Formation condition and distri—bution prediction of gas hydrate in Qinghai—Tibet Plateau per—mafrost[J].Chinese Journal of Geophysics，2005，48（1）：179一187.

[62]祝有海，趙省民，盧振權(quán)，等，中國凍土區(qū)天然氣水合物的找礦選區(qū)及其資源潛力[J].天然氣工業(yè)，2011，31（1）：13一l9.

[63]Zhao X，Deng J，Li J，et al.Gas hydrate formation and its accu—mulation potential in Mohe permafrost，China[J].Marine and pe.troleum geology，2012，35（1）：166一l75.

[64]劉杰，楊睿，鄔黛黛，等.基于生烴思路的微生物成因水合物資源量估算——以瓊東南盆地西南深水區(qū)為例[J].天然氣地球科學(xué)，2019，30（4）：539—548.