天然氣水合物儲層與資源量評價現狀與展望

馬文嘉

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

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

中圖分類號：TH455

文獻標志碼：A

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

0 引言

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

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

1 天然氣水合物儲層評估

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

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

1.1 地球物理評估

1.1.1 地震數據

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

1.1.2 隨鉆測井法

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

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

1.1.3 海洋可控源電磁法

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

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

1.2 地球化學評估

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

1.2.1 氯離子濃度法

氯離子濃度法是利用水合物所在孔隙與海水中的氯離子濃度具有差異，因為天然氣水合物在形成時會使原位孔隙中的氯化物脫除，并隨之擴散開，巖心取樣期間隨著水合物的分解，氯離子濃度會進一步降低。Yuan等[22]利用氯離子濃度和水合物密度，對海洋鉆探計劃（ODP）的水合物飽和度進行了預測。Liu等[23]利用聲波和電阻在南海北部神狐海域測定的水合物飽和度變化，與利用氯離子濃度測定得到的水合物飽和度規律呈現出一致性。Chen等[24]基于徑向基函數插值的數學描述法，以南海神狐SH2的41個取心樣品以及SH7的21個取心樣品作為數據，使用氯離子濃度相關的阿爾奇公式計算水合物飽和度，得到的飽和度與利用聲波或者電阻計算的結果近似。

1.2.2 生物地球化學法

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

1.3 鉆探取樣

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3 結束語

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

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

參考文獻

[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在天然氣水合物勘查中的應用[J].海洋技術學報，2018，37（2）：6.

[8]景鵬飛，胡高偉，卜慶濤，等.天然氣水合物地球物理勘探技術的應用及發展[J].地球物理學進展，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]高興軍，于興河，李勝利，等，地球物理測井在天然氣水合物勘探中的應用[J].地球科學進展，2003，18（2）： 305-311.

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

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

[13]孫建孟，羅紅，焦滔，等.天然氣水合物儲層參數測井評價綜述[J].地球物理學進展，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識別天然氣水合物的優缺點對比[J].新能源進展，2018， 6（5）：112-124.

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

[18]楊蜀冀，胡楊承鈺，宋紅喜，等.海洋可控源電磁探測系統概述[J].電氣時代，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]張富貴，孫忠軍，楊志斌，等，中國凍土區天然氣水合物的地球化學調查[J].礦物巖石地球化學通報，2019， 38（6）：11.

[21]單晨晨，國家重點研發計劃項目”近海底高精度水合物探測技木”完成首次深海聯合測試[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]付強，王國榮，周守為，等.海洋天然氣水合物開采技術與裝備發展研究[J].中國工程科學，2020， 22（6）：8.

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

[29]叢曉榮，吳能友，蘇明，等.天然氣水合物資源量估算研究進展及展望[J].新能源進展，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]. 海洋學研究， 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].熱帶海洋學報，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]唐勇，方銀霞，高金耀，等，沖繩海槽天然氣水合物穩定帶特征及資源量評價[J]，海洋地質與第四紀地質，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]祝有海，趙省民，盧振權，等，中國凍土區天然氣水合物的找礦選區及其資源潛力[J].天然氣工業，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]劉杰，楊睿，鄔黛黛，等.基于生烴思路的微生物成因水合物資源量估算——以瓊東南盆地西南深水區為例[J].天然氣地球科學，2019，30（4）：539—548.