天然氣水合物實驗模擬裝置和數值模擬器的研究進展

DOI: 10.19911/j.1003-0417.tyn20240613.01　文章編號：1003-0417(2025)01-33-10

天然氣水合物實驗模擬裝置和數值模擬器的研究進展

李 剛 1,2，李小森 1,2，翁一凡 1,2，呂秋楠 1,2*

(1. 中國科學院廣州能源研究所，廣州 510640；2. 中國科學院天然氣水合物重點實驗室，廣州 510640)

摘要：中國的天然氣水合物資源豐富，若實現商業化開採，有助於中國的能源結構調整。天然氣水合物開採的實驗模擬和數值模擬是指導其安全高效開採的重要手段。首先簡要介紹了天然氣水合物的開採方法、國內外野外試驗性開採案例，然后系統總結了天然氣水合物實驗模擬裝置及全球具有代表性的天然氣水合物數值模擬器的開發和應用，最后對未來天然氣水合物實驗模擬和數值模擬的研究方向進行了展望。未來，在實驗模擬裝置方面，應重點研發適用於開採過程中測量和監測熱學、滲流、力學和相變參數的實驗模擬裝置；在數值模擬器方面，應實現天然氣水合物數值模擬過程的精確化，並不斷優化數值模擬器算法，從而建立可服務於中國天然氣水合物產業化開發的高精度數值模擬器。

關鍵詞：天然氣水合物；實驗模擬裝置；數值模擬器；理論模型

中圖分類號：TE52　文獻標誌碼：A

天然氣水合物，俗稱「可燃冰」，是由天然氣和水分子在高壓低温條件下形成的固體化合物，其碳含量為常規化石燃料( 煤炭、石油和天然氣等 ) 中碳含量總和的兩倍 [1]，具有碳儲量大、能量密度高等特點，是一種重要的潛在能源。天然氣水合物大多數賦存在陸地永凍土帶及深海海域，對其進行高效開發和利用有助於緩解能源短缺。因此，越來越多的勘探計劃得以實施。部分天然氣水合物勘探計劃的分佈圖如圖1[2-11]所示。

　　天然氣水合物開採的實驗模擬和數值模擬是指導其安全高效開採的重要手段。基於此，本文首先對現階段常用的幾種天然氣水合物開採方法及全球各種野外試驗性開採(下文簡稱為「試採」)案例進行分析，然后對國內外天然氣水合物實驗模擬裝置和數值模擬器的開發和應用進行介紹，最后對未來天然氣水合物實驗和數值模擬的研究方向進行展望。

1 天然氣水合物開採方法

　　天然氣水合物的開採是一項極具挑戰性且複雜的工程，其中涉及水合物相變、多相滲流、儲層傳熱、力學 4 個相互影響的物理、化學效應 [12]，是通過破壞水合物晶體結構的穩定性，從而促使天然氣水合物分解為天然氣和水。為此，研究人員研究了一系列天然氣水合物的開採方法，包括降壓法、注熱法、注入抑制劑法、CO2 置換法、固態流化法及聯合開採法，這些方法有各自的適用性和優缺點。

　　1) 降壓法。通過抽水實現天然氣水合物儲層壓力的降低，使儲層壓力低於天然氣水合物的相平衡壓力，促使天然氣水合物分解，從而實現對天然氣水合物的開採。此方法具有開採效率高、操作難度低等優點；但其缺點是不適用於滲透率較低的儲層，且其對熱量的需求極大。當顯熱不足時，天然氣水合物的分解速率降低，甚至會出現天然氣水合物二次生成或結冰的現象。此外，利用降壓法開採天然氣水合物通常會伴隨出砂，最終導致開採井口堵塞 [13-14]。

　　2) 注熱法。通過向天然氣水合物儲層注入蒸汽、熱水或其他熱流體，為天然氣水合物的分解提供熱量，從而實現天然氣水合物的開採。該方法操作簡易，對設備的要求低且易於控制；但其缺點是在滲透率低、熱利用率低的環境下，無法實現天然氣水合物的高效開採 [15-16]。

　　3) 注入抑制劑法。通過向天然氣水合物儲層注入甲醇、乙二醇或鹽水等抑制劑，破壞其相平衡，從而促使天然氣水合物分解，實現天然氣水合物的開採。該方法具有操作簡易、能量消耗小的優點；但其缺點是成本高，對環境也可能會造成污染，且作用效果較慢。注入抑制劑法常用於防止天然氣水合物開採井井口因生成水合物而造成的堵塞，或用於輔助其他的天然氣水合物開採方法 [17-18]。

　　4) CO2 置換法。由於 CO2 水合物更容易生成，因此通常是向天然氣水合物儲層中注入 CO2來置換出水合物籠中的 CH4 氣體分子，從而實現天然氣水合物的開採。此方法不僅可以保證儲層的穩定性，還可以實現 CO2 的封存。然而，該方法受儲層的傳熱和滲流的限制，導致置換速率慢，現階段仍處於實驗模擬研究階段 [19-20]。

　　5) 固態流化法。通過挖掘和機械破碎深水淺層尚未固結的天然氣水合物礦藏，在舉升管道中實現對深水淺層天然氣水合物安全可控的開發。該方法具有天然氣水合物開採產量大、開採效率高等特點；但其缺點是開採過程可能會破壞海牀的穩定性 [21]。

　　6) 聯合開採法。該方法是指基於現場實際情況將上述兩種及兩種以上的開採方法相結合，例如：將降壓法和注熱法相結合的降壓注熱法，以及注熱法和注入抑制劑 ( 鹽水 ) 法相結合的注熱鹽水法，通過充分利用各種方法的優勢，實現天然氣水合物的高效開採 [22-23]。

2 國內外天然氣水合物試採工程

　　全球已經有一些國家分別在陸地和海洋開展了天然氣水合物的試採工程，下文針對部分工程進行詳細介紹。

2.1 中國

　　中國先后在祁連山木里地區和南海海域開展了天然氣水合物試採作業，具體工程包括：

　　1) 祁連山木里地區 ( 凍土區 ) 試採。2011 年9 月，中國地質調查局在祁連山木里地區開展了第 1 次陸地天然氣水合物試採作業。本次試採作業累計時長 101 h；通過利用單直井和降壓注熱法對凍土區天然氣水合物進行開採，累計採氣量達 95 m3 。2016 年 9 月，在該地區開展了第 2 次天然氣水合物試採作業，此次試採作業共持續23 天；通過利用水平井和降壓法進行天然氣水合物開採，累計採氣量達 1078.4 m3 。兩次試採作業的結果表明，水平井可以有效提高採氣量。

　　2) 南海北部荔灣區域試採。2017 年 5 月，中國海洋石油集團有限公司在南海北部荔灣 3 站位利用固態流化法對海洋淺層非成巖天然氣水合物進行了試採 [21]，採氣量達 81 m3 ，驗證了固態流化法開採天然氣水合物的可行性。

　　3) 南海神狐海域試採。2017 年 3 月，中國地質調查局在南海神狐海域第 1 次成功實施了天然氣水合物試採，累計採氣量達 30.9×104 m3 ，日均採氣量為 5151 m3 。本次試采采用垂直井和降壓法，但在試採后期，由於泥質粉砂儲層滲透率低，導致后期採氣量降低 [2]。2019 年 10 月，中國地質調查局在南海神狐海域進行了第 2 次天然氣水合物試採，連續採氣 30 天，累計採氣量達 86.14×104 m3 ，日均採氣量為 2.87×104 m3 。本次試採是國際上首次利用水平井對海域的天然氣水合物進行開採，並配套採用降壓法輔助注熱法和注入抑制劑法。在試採后期，同樣出現了儲層滲透率低，大量產水的現象，可以通過儲層改造來解決該問題 [3]。

2.2 美國

　　美國先后於 2005、2007 年在墨西哥灣海域和阿拉斯加北坡實施了天然氣水合物的鑽探、勘察和取樣工作。2012 年，美國在阿拉斯加北坡採用 CO2 置換法對天然氣水合物進行試採，試採時間持續 6 周，但日均採氣量較低，共計產出天然氣 2.5×104 m3 ，這主要是由 CO2 置換速率低造成的 [4]。

2.3 加拿大

　　2002、2007 年，加拿大先后在 Mallik 地區的凍土區使用注熱法和降壓法進行了兩次天然氣水合物試採，均出現了產氣速率低、出砂等問題。2008 年，加拿大在該凍土區採用垂直井和降壓法再次進行了天然氣水合物試採，試採共持續 6天，累計採氣量達 13000 m3 ；但本次試採仍是因為出砂原因被迫終止 [5]。

2.4 日本

　　日本於 2013 年 3 月在近海海域進行為期 6天的天然氣水合物試採，為全球首次在近海海域進行的天然氣水合物試採工程，通過使用垂直井和降壓法，得到的總採氣量約為 2×104 m3 ；后期的試採作業因出砂和產水原因被迫終止。2017年 4 月 7 日，日本開展第 2 次近海海域的天然氣水合物試採作業，為期 32 天，累計採氣量約為23.5×104 m3 ，採用兩口生產井交替採氣；后期同樣出現了出砂堵塞井口、大量產水等問題，導致採氣量極低 [6-7]。

2.5 前蘇聯

　　20 世紀 60 年代，前蘇聯西伯利亞北部的麥索亞哈氣田被證實有大量天然氣水合物存在。1972年，前蘇聯通過採用降壓法和注入抑制劑法相結合的方法使該氣田的天然氣水合物儲層開始分解產氣，到 2011 年底，該氣田已累計採集 1.29×1010m3天然氣。由於存在遊離天然氣。因此總採氣量中預計有 41.9% 來自天然氣水合物的分解 [8]。

綜上可知，上述國家開展的天然氣水合物試採工程都存在開採效率低、經濟性差的問題，説明天然氣水合物開採的產業化還有許多關鍵性科學問題和技術問題有待解決。

3 實驗模擬裝置

　　爲了提高開採效率和解決開採技術難題，需要通過實驗模擬天然氣水合物的生成和分解過程，從而解決相應問題。因此，各種中型、大型天然氣水合物實驗模擬裝置得以開發，並結合各種有效的表徵手段，揭示多孔介質中天然氣水合物分解的微觀機理，為理論模型的建立提供基礎參數。

　　目前，全球相對完善和系統化的天然氣水合物模擬裝置主要具有以下特點：1) 耐低温、抗高壓，可進行天然氣水合物的生成和分解實驗；2)可視化程度高，可直接觀察模擬裝置內的天然氣水合物相變過程；3) 測試精度高，可準確測出天然氣水合物生成和分解過程中壓力和温度的變化情況；4) 檢測手段多樣化，包括光、聲、電多種檢測方法。

　　模擬裝置具有上述特點是爲了確保其可以模擬各種環境條件，進而更好地研究多孔介質中天然氣水合物的分解過程，獲得重要的物性參數，為野外調查提供科學依據。全球主要的大型天然氣水合物實驗模擬裝置概況如表 1 所示。

　　中國的全尺寸開採井天然氣水合物三維綜合試驗開採系統是目前全球規模最大的天然氣水合物實驗模擬裝置，其實物圖如圖 2 所示。該裝置包含 245 個測温點、245 個測壓點，可承受 32MPa 的壓力，井筒內配備有攝像頭能夠實時觀察井筒內的情況，其能夠模擬 3000 m 水深條件下天然氣水合物的生成和分解。

　　天然氣水合物生成或分解過程的精準量化直接決定了實驗結果的可靠性，需依據實際的地層温度、壓力及儲層條件，採用不同的布井方式和開採方法進行天然氣水合物的開採實驗研究。通過對多孔介質中天然氣水合物分解進行動態監測，分析產氣、產水規律，研究影響開採效率的關鍵因素，為天然氣水合物數值模擬模型的建立和實際開採提供可靠的實驗數據和參考依據。

4 數值模擬器的開發及應用

4.1 數值模擬器的開發

　　為便於評估不同天然氣水合物開採方法的可行性及天然氣水合物的開採潛力，國內外開發了一系列天然氣水合物數值模擬器(或模擬軟件)。數值模擬器可以根據實際天然氣水合物的賦存條件，對不同開採方式下天然氣水合物的開採過程進行模擬。當前的天然氣水合物開採數值模擬主要採用兩種技術路徑：1) 基於傳統石油和天然氣開採的相關數學模型，例如：在傳統黑油 (blackoil) 模型和成熟商業軟件 COMSOL 等的基礎上，將液相的石油直接轉變為固相的天然氣水合物。2) 專門針對天然氣水合物的成藏和開採，以及CO2 水合物地質封存進行數值模擬，並考慮了水合物生成或分解過程的放熱或吸熱化學反應過程，以及多孔介質中多相滲透和傳熱的複雜過程。國際上主要的天然氣水合物數值模擬器 ( 或模擬軟件 ) 如表 2 所示。

　　Wilder 等 [44] 和 White 等 [34] 分別於 2008 年和 2020 年開展了天然氣水合物數值模擬器的研究對比，結果顯示：從國際主流數值模擬器的建模方法、主要參數選取及求解結果的準確性來看，目前國際上尚缺乏公認的、能夠精確有效分析含水合物沉積物這一複雜體系的實用工具。尤其是求解涉及複雜的傳熱 ( 導熱、滲流換熱和擴散換熱 )、多相滲流 ( 氣液兩相的相對滲透率、含水合物多孔介質的絕對滲透率 )，以及相變化學反應 ( 水合物生成和分解反應 ) 等的強耦合問題時，在物理和化學過程描述和迭代收斂等方面還存在較大缺陷。主要體現在温度、壓力、各相飽和度分佈等求解結果不合理，收斂性欠佳，模型參數的選取隨水合物變化範圍較大等，因此，有待於進一步深入研究。

4.2 數值模擬器的應用

　　通過天然氣水合物數值模擬器對實驗結果的復現，可以驗證天然氣水合物物理模型和物性參數的合理性和適用性，這些模型和參數也將用於野外天然氣水合物的評估模擬。

　　Li 等 [45] 基於實驗數據提出了 1 種新的動力學模型，將其植入到「TOUGH+HYDRATE」模擬器中，通過數值模擬對比了 3 套不同尺度實驗裝置中天然氣水合物分解過程的產氣、產水、壓力及温度的變化規律。Li 等 [40-41] 使用全隱式水合物模擬器分別模擬了天然氣水合物三維綜合試驗開採系統和立方體天然氣水合物模擬裝置(CHS) 中的天然氣水合物生成和分解過程，並得到了一系列能夠用於后續研究的物性參數。

　　Moridis[46] 使用「TOUGH+HYDRATE」模擬器模擬研究了阿拉斯加北坡的天然氣水合物儲層的產氣潛力。模擬結果表明：在井底壓力恆定的情況下，水平井的產氣量與垂直井相比幾乎增加了兩個數量級，但產量仍然很低，而沉積層的初始温度可能是決定產氣潛力的重要因素。Konno等 [47] 基於第 1 次在日本南海海槽東部進行的海域天然氣水合物試採作業的數據進行了數值模擬研究，模擬結果表明：有效滲透率相對較高的含高飽和度天然氣水合物的砂岩層更適合採用降壓法開採，具有該特點的天然氣水合物儲層也是商業開採的首選。Yamamoto 等 [48] 基於日本南海海槽第 2 次 (2017 年 ) 海域天然氣水合物試採的温度和壓力數據，使用 MH21-HYDRES 數值模擬器對試採過程進行模擬研究，發現儲層內天然氣水合物的不均勻分佈是造成產氣量模擬結果和實際結果存在較大差異的主要原因，儲層內水的流動給天然氣水合物開採過程的分析增加了複雜性。Yu 等 [49-53] 根據日本南海海槽的地質勘探資料及海域天然氣水合物試採現場數據，利用「TOUGH+HYDRATE」模擬器預測了長達 5 年開採時間下天然氣水合物開採過程中的產氣行為，分別研究了採用雙垂直井和降壓法開採，以及水平井和降壓法開採，發現通過提高天然氣水合物儲層滲透率和採用降壓注熱法等方式可提高天然氣的產出率。

　　Xiao 等 [54] 以 2020 年中國南海神狐海域第 2次天然氣水合物試採數據為基礎，研究了滲透率對天然氣水合物儲層生產行為的影響，提出了滲透率修正模型，首次在野外尺度上分析了天然氣水合物分解對儲層滲透率的影響。分析結果顯示：天然氣水合物的分解可以提高儲層滲透率，從而促進流體流動，提高產氣效率。Zhao 等 [55] 使用「TOUGH+HYDRATE」模擬器在中國瓊東南海域的地質勘探數據基礎上，對天然氣水合物與淺層氣共存體系進行了長期的開採模擬，對比了多氣合採和單一開採的產氣效率。研究結果表明：遊離氣的存在可以顯著提高天然氣水合物儲層的產氣能力。Li 等 [42] 使用全隱式水合物模擬器對中國南海神狐海域 SHSC2-6 區域進行了 8 年內的天然氣開採潛力評估，發現天然氣水合物分解過程的吸熱效應遏制了降壓法開採的持續進行，提出降壓注熱法是具有潛力的海域天然氣水合物開採方法。

　　利用天然氣水合物數值模擬器復現實驗結果，獲得可靠的物理模型和物性參數，並將其用於野外天然氣水合物的產氣能力評估及獲得高效的開採方法。

　　在未來的數值模擬器研究中，一方面，應當基於實驗數據獲取更加可靠的物理模型和物性參數，建立可靠的多相滲流模型及相變 - 傳熱 - 滲流 - 力學強耦合理論模型，實現數值模擬過程的精確化；另一方面，應當實現數值模擬器的優化集成，通過並行計算、中央處理器 (CPU) 與圖形處理器 (GPU) 結合、多處理器協同計算等方式，在保障計算精度的條件下實現超大網格矩陣的求解，並實現三相四場耦合模型的快速求解。

5 展望

　　現階段，中國在天然氣水合物野外試採和實驗室尺度下實驗模擬裝置的研發方面都處於領先水平，但針對數值模擬器的研究起步較晚。未來在實驗模擬裝置方面，應重點研發適用於開採過程中測量和監測熱學、滲流、力學和相變參數的實驗模擬裝置，通過相變 - 傳熱 - 滲流 - 力學多場耦合一體化監測裝備的建立，實現熱導率、熱容、相變熱、應力、滲透率等多項物性參數的長期一體化監測。在數值模擬器方面，未來應當將前期忽略或尚未考慮的影響天然氣水合物開採的因素引入到物理模型中進行修正，並基於實驗數據獲取更為精準的物理模型和物性參數，建立多相滲流模型及相變 - 傳熱 - 滲流 - 力學強耦合理論模型，實現天然氣水合物開採數值模擬過程的精確化；突破常規數值模擬器在均質、小尺度模型網格方面的限制，實現天然氣水合物的精細化描述；實現儲層、井筒、管道的全過程集成化模擬；優化數值模擬器算法，從而建立可服務於中國天然氣水合物產業化開發的高精度數值模擬器。

6 結論

　　本文首先對現階段常用的幾種天然氣水合物開採方法及全球各種野外試採案例進行了分析，然后對國內外天然氣水合物實驗模擬裝置和數值模擬器的開發及應用進行了介紹，最后對未來天然氣水合物實驗模擬和數值模擬的研究方向進行了展望。未來，在實驗模擬裝置方面，應重點研發適用於開採過程中測量和監測熱學、滲流、力學和相變參數的實驗模擬裝置；在數值模擬器方面，應實現天然氣水合物開採數值模擬過程的精確化，並不斷優化數值模擬器算法，從而建立可服務於中國天然氣水合物產業化開發的高精度數值模擬器。

[ 參考文獻 ]

[1] KVENVOLDEN K A. Gas hydrates—geologicalperspective and global change[J]. Reviews of geophysics，1993，31(2)：173-187.

[2] YE J L，QIN X W，QIU H J，et al. Data report：molecular and isotopic compositions of the extracted gasfrom China’s first offshore natural gas hydrate productiontest in South China Sea[J]. Energies，2018，11(10)：2793.

[3]葉建良，秦緒文，謝文衞，等 . 中國南海天然氣水合物第二次試採主要進展 [J]. 中國地質，2020，47(3)：557-568.

[4] BOSWELL R，SCHODERBEK D，COLLETT T S，et al. The Inik Sikumi field experiment，Alaska NorthSlope：design，operations，and implications for CO2-CH4exchange in gas hydrate reservoirs[J]. Energy & fuels，2017，31(1)：140-153.

[5] ANDERSON B，COLLETT T，LEWIS R，et al. Usingopen hole and cased-hole resistivity logs to monitor gashydrate dissociation during a thermal test in the Mallik5L-38 research well，Mackenzie Delta，Canada[J].Petrophysics，2008，49(3)：285-294.

[6] YAMAMOTO K，TERAO Y，FUJII T，et al. Operationaloverview of the first offshore production test of methanehydrates in the Eastern Nankai Trough[C]//OffshoreTechnology Conference 2014，May 5-8，2014，Houston，Texas，USA. [S.l.：s.n.]，2014.

[7] YAMAMOTO K，WANG X X，TAMAKI M，et al.The second offshore production of methane hydrate inthe Nankai trough and gas production behavior froma heterogeneous methane hydrate reservoir[J]. RSCadvances，2019，9(45)：25987-26013.

[8] MAKOGON Y F，OMELCHENKO R Y. Commercial gasproduction from Messoyakha deposit in hydrate conditions[J].Journal of natural gas science and engineering，2013，11：1-6.

[9] A S K M，CAMERLENGHI A，KOPF A，e t a l .Scientific ocean drilling to assess submarine geohazardsalong European margins[C]//2008 American GeophysicalUnion(AGU) Fall Meeting，December 15-19，2008，SanFrancisco, United States. [S.l.：s.n.]，2008.

[10] 祝有海，張永勤，文懷軍，等 . 青海祁連山凍土區發現天然氣水合物 [J]. 地質學報，2009，83(11)：1762-1771.

[11] 申鵬飛 . 含天然氣水合物沉積物滲流實驗及增產方法研究 [D]. 重慶：重慶大學，2021.

[12] MORIDIS G，MORIDIS G J，KOWALSKY M B，et al. TOUGH+HYDRATE v1.0 user’s manual：a codefor the simulation of system behavior in hydrate-bearinggeologic media[R]. Berkeley：Lawrence Berkeley NationalLaboratory，2008.

[13] LI G，MORIDIS G J，ZHANG K N，et al. Evaluation ofgas production potential from marine gas hydrate depositsin Shenhu area of South China Sea[J]. Energy & fuels，2010，24(11)：6018-6033.

[14] ZHANG Y，LI X S，CHEN Z Y，et al. Methane hydratedissociation by depressurization in sediments with differenthydrate saturations in cubic hydrate simulator[J]. Energyprocedia，2014，61：990-994.

[15] WANG Y，FENG J C，LI X S. Experimental investigationinto methane hydrate dissociation by thermal stimulationwith dual vertical well[J]. Energy procedia，2017，105：4738-4744.

[16] LI G，MORIDIS G J，ZHANG K N，et al. The use ofhuff and puff method in a single horizontal well in gasproduction from marine gas hydrate deposits in the Shenhuarea of South China Sea[J]. Journal of petroleum scienceand engineering，2011，77(1)：49-68.

[17] PALERMO T，SLOAN D. Artificial and natural inhibitionof hydrates[M]//SLOAN D，KOH C，SUM A K，et al. Natural gas hydrates in flow assurance. Houston：Gulf professional publishing，2011.

[18] FARHADIAN A，NAEIJI P，VARFOLOMEEV M A，etal. Reconsideration of the micellization theory：promotionor inhibition of gas hydrate formation for gas storage andflow assurance applications[J]. Chemical engineeringjournal，2022，427：131852.

[19] XU C G，CAI J，YU Y S，et al. Effect of pressure onmethane recovery from natural gas hydrates by methanecarbon dioxide replacement[J]. Applied energy，2018，217：527-536.

[20] WA N G X H，SUN Y F，WANG Y F，e t a l . G a sproduction from hydrates by CH4-CO2/H2 replacement[J].Applied energy，2017，188：305-314.

[21] 周守為，陳偉，李清平，等 . 深水淺層非成巖天然氣水合物固態流化試採技術研究及進展 [J]. 中國海上油氣，2017，29(4)：1-8.

[22] S O N G Y C，CHENG C X，ZHAO J F，e t a l .Evaluation of gas production from methane hydrates usingdepressurization，thermal stimulation and combinedmethods[J]. Applied energy，2015，145：265-277.

[23] FENG J C，WANG Y，LI X S. Hydrate dissociationinduced by depressurization in conjunction with warm brinestimulation in cubic hydrate simulator with silica sand[J].Applied energy，2016，174：181-191.

[24] LI G，LI B，LI X S，et al. Experimental and numericalstudies on gas production from methane hydrate inporous media by depressurization in pilot-scale hydratesimulator[J]. Energy & fuels，2012，26(10)：6300-6310.

[25] SCHICKS J M，SPANGENBERG E，GIESE R，et al.New approaches for the production of hydrocarbons fromhydrate bearing sediments[J]. Energies，2011，4(1)：151-172.

[26] LI N，ZHANG J，XIA M J，et al. Gas productionfrom heterogeneous hydrate-bearing sediments bydepressurization in a large-scale simulator[J]. Energy，2021，234：121183.

[27] LI G，LI X S，WANG Y，et al. Production behaviorof methane hydrate in porous media using huff and puffmethod in a novel three-dimensional simulator[J]. Energy，2011，36(5)：3170-3178.

[28] KONNO Y，JIN Y，SHINJOU K，et al. Experimentalevaluation of the gas recovery factor of methane hydrate insandy sediment[J]. RSC advances，2014，4(93)：51666-51675.

[29] LOH M，TOO J L，FALSER S，et al. Gas productionfrom methane hydrates in a dual wellbore system[J]. Energy& fuels，2015，29(1)：35-42.

[30] MEREY S，SINAYUC C. Experimental set-up design forgas production from the Black Sea gas hydrate reservoirs[J].Journal of natural gas science and engineering，2016，33：162-185.

[31] MORIDIS G J，KOWALSKY M B，PRUESS K.TOUGH+HYDRATE v1.1 user’s manual：a code for thesimulation of system behavior in hydrate bearing geologicmedia[R]. Berkeley：Lawrence Berkeley NationalLaboratory，2009.

[32] MORIDIS G J，KOWALSKY M B，PRUESS K.Depressurization-induced gas production from class 1hydrate deposits[J]. SPE reservoir evaluation & engineering，2007，10(5)：458-481.

[33] WHITE M，OOSTROM M. STOMP subsurface transportover multiple phases，version 4.0，user’s guide[R].Washington：United States Department of Energy，2006.

[34] WHITE M D，KNEAFSEY T J，SEOL Y，et al. Aninternational code comparison study on coupled thermal，hydrologic and geomechanical processes of natural gashydrate-bearing sediments[J]. Marine and petroleumgeology，2020，120：104566.

[35] NARITA H. Introduction of MH21 (research consortiumfor methane hydrate resources in Japan) and current topicsin production method & modeling of methane hydrate[C]//The Proceedings of the Fifth(2003) ISOPE Ocean MiningSymposium，September 15-19，2003，Tsukuba，Japan. [S.l.：s.n.]，2003.

[36] K O N N O Y，MASUDA Y，AKAMINE K，e t a l .Sustainable gas production from methane hydrate reservoirsby the cyclic depressurization method[J]. Energy conversionand management，2016，108：439-445.

[37] K I M A R，KIM J T，CHO G C，e t a l . M e t h a n eproduction from marine gas hydrate deposits in Korea：thermal-hydraulic-mechanical simulation on productionwellbore stability[J]. Journal of geophysical research：solidearth，2018，123(11)：9555-9569.

[38] ZATSEPINA O Y，POOLADI-DARVISH M. Storage ofCO2 as hydrate in depleted gas reservoirs[J]. SPE reservoirevaluation & engineering，2012，15(1)：98-108.

[39] J I N G R，LEI H W，XU T F，e t a l . S i m u l a t e dgeomechanical responses to marine methane hydraterecovery using horizontal wells in the Shenhu area，SouthChina Sea[J]. Marine and petroleum geology，2018，92：424-436.

[40] LI G，LI X S，LV Q N，et al. Full implicit simulatorof hydrate (FISH) and analysis on hydrate dissociation inporous media in the cubic hydrate simulator[J]. Energy，2023，280：128191.

[41] LI G，Lyu Q N，LI X S，et al. Analysis of methanehydrate dissociation experiment in a pilot-scale hydratesimulator by a full implicit simulator of hydrate[J].Industrial & engineering chemistry research，2023，62(19)：7704-7720.

[42] LI G，LI X S，LV Q N，et al. Long-term analysis of the2020 gas hydrate production test in the South China Seaby full implicit simulator of hydrate[J]. Energy & fuels，2023，37(5)：3785-3798.

[43] LI G，ZHAO Z C，LI X S，et al. Control factors for gasrecovery in the second methane hydrate field test in theNankai trough[J]. Energy & fuels，2023，37(18)：13985-13996.

[44] WILDER J W，MORIDIS G J，WILSON S J，et al.An international effort to compare gas hydrate reservoirsimulators[C]//Proceedings of the 6th InternationalConference on Gas Hydrates，July 6-10，2008，Vancouver，British Columbia，Canada. [S.l.：s.n.]，2008.

[45] LI B，LI X S，LI G. Kinetic studies of methane hydrateformation in porous media based on experiments in a pilotscale hydrate simulator and a new model[J]. Chemicalengineering science，2014，105：220-230.

[46] MORIDIS G J. Gas production from a cold，stratigraphicallybounded hydrate deposit at the Mount Elbert site，NorthSlope，Alaska[J]. Marine and petroleum geology，2009，28(2)：517-534.

[47] KONNO Y，FUJII T，SATO A，et al. Key findings of theworld’s first offshore methane hydrate production test off thecoast of Japan：toward future commercial production[J].Energy & fuels，2017，31(3)：2607-2616.

[48] YAMAMOTO K，WANG X X，TAMAKI M，et al.The second offshore production of methane hydrate inthe Nankai trough and gas production behavior froma heterogeneous methane hydrate reservoir[J]. RSCadvances，2019，9(45)：25987-26013.

[49] YU T，GUAN G Q，ABUDULA A，et al. Heatassisted production strategy for oceanic methane hydratedevelopment in the Nankai trough，Japan[J]. Journal ofpetroleum science and engineering，2019，174：649-662.

[50] YU T，GUAN G Q，ABUDULA A，et al. Application ofhorizontal wells to the oceanic methane hydrate productionin the Nankai trough，Japan[J]. Journal of natural gasscience and engineering，2019，62：113-131.

[51] YU T，GUAN G Q，ABUDULA A，et al. Enhanced gasrecovery from methane hydrate reservoir in the Nankaitrough，Japan[J]. Energy procedia，2019，158：5213-5218.

[52] YU T，GUAN G Q，ABUDULA A，et al. Gas recoveryenhancement from methane hydrate reservoir in the Nankaitrough using vertical wells[J]. Energy，2019，166：834-844.

[53] Y U T，GUAN G Q，A B U D U L A A . P r o d u c t i o nperformance and numerical investigation of the 2017offshore methane hydrate production test in the Nankaitrough of Japan[J]. Applied energy，2019，251：113338.

[54] XIAO C W，LI X S，LI G，et al. Numerical analysis ofproduction behaviors and permeability characteristics onthe second gas hydrate production test in the South ChinaSea[J]. Energy & fuels，2022，36(18)：10960-10974.

[55] ZHAO Q，LI X S，CHEN Z Y，et al. Numericalinvestigation of production characteristics and interlayerinterference during co-production of natural gas hydrate andshallow gas reservoir[J]. Applied energy，2024，354：122219.

（轉自：太陽能雜誌）