含水合物沉積物多相滲流特性研究進展1)

蔡建超 夏宇軒 徐 賽 田海濤

?(中國石油大學油氣資源與探測國家重點實驗室，北京 102249)

?(中國地質大學地球物理與空間信息學院，武漢 430074)

??(中國地質大學資源學院，武漢 430074)

引言

天然氣水合物是氣體分子與水在低溫、高壓條件下形成的類冰狀的非化學計量的籠型結晶水合物，外形如冰雪狀，通常呈現白色，遇火即可燃燒，俗稱“可燃冰”.天然氣水合物全球范圍內主要存在于深水大陸架的淺層沉積物及陸上凍土中[1-2].作為一種資源量豐富的高效清潔能源，天然氣水合物具有潛在的巨大經濟效益，對于世界能源格局、全球氣候環境以及自然災害防治等一系列問題有重要的影響，是近些年來世界范圍內大量水合物研究和開發項目蓬勃發展的驅動力[3-6].全球當前已探明水合物藏主要分布在日本海、印度洋、墨西哥灣、中國南海、白令海峽、韓國海域、特立尼達島和多巴哥島等海域及阿拉斯加、麥肯齊三角洲、西伯利亞和青藏高原凍土區[7-9].20 世紀60 年代末，前蘇聯西伯利亞西北部的麥索亞哈氣田進行商業性開采，氣田上方存在天然氣水合物儲層，后進行降壓法開采[10-11].20 世紀80 年代到90 年代，美國開展對墨西哥灣北部天然氣水合物研究工作[12]；2013 年，日本對南海海槽進行甲烷水合物第一次生產測試[13].中國地質調查局于2007 年4 月在神狐海域首次獲得分散狀天然氣水合物樣品[14-15].2017 年，基于近20年來對南海儲層特征的系統研究，中國地質調查局在南海神狐海域進行首次天然氣水合物試采，歷時60 天，產氣量3.09×105m3，是世界上產氣連續時間最長、累計產氣量最大的水合物試采工程[16].

天然氣水合物對世界能源需求的貢獻取決于開采技術與成本，因此合理的開采方法尤為重要.目前可行的天然氣水合物開采方式主要有降壓開采法、注熱開采法、化學抑制法、二氧化碳置換法以及多種方法聯合開采等[17-23].水合物開發過程是多組分多相流的動態過程，水合物分解過程中會導致儲層力學、熱學及物性等參數發生變化.鉆井開發過程中由于地層水合物發生分解，井壁周圍地層壓力和溫度發生變化，容易引起井壁失穩、井壁坍塌等井眼穩定性問題，造成生產井出砂，影響水合物開采過程[24-27].因此，闡明水合物開采過程中水合物分解及氣水滲流特性有助于生產開發措施的進一步完善.

水合物開采過程中，流體輸運過程決定生產效率.水合物儲層滲流過程伴隨水合物分解引起的氣體釋放及液體生成引起的儲層物性改變.準確刻畫水合物開采過程中的滲流特性成為當前研究熱點之一，國內外許多專家學者對含水合物多孔介質滲流過程進行了理論、實驗及數值模擬等相關研究[28-32].研究含水合物沉積物滲流特性，正確認識開發過程中沉積物內復雜流動規律，對我國天然氣水合物開采具有重要的理論價值和實際意義.本文系統歸納了水合物滲流平行毛細管模型、Kozeny 顆粒模型等理論模型，闡述含水合物沉積物有效滲透率測試、滲流過程中沉積物物性演變以及水合物開采室內模擬實驗研究成果，總結水合物滲流多尺度數值模擬研究及水合物產氣預測模擬最新進展，從理論建模到數值模擬，從室內實驗到大尺度開發，多個角度為水合物儲層滲流特性研究提供思路與展望.

1 基礎滲流理論

含水合物沉積物絕對滲透率及孔隙流體相對滲透率是多相滲流效應的關鍵參數，它們影響沉積物中氣體、水及其中溶解氣與鹽的含量，又反過來控制水合物的合成、分解以及氣體生產的過程.特別地，由于水合物在孔隙空間中的增長而導致的滲透率降低對與這些過程相關的時間尺度具有重要的控制作用[33-34].這種效應通常用歸一化滲透率kr來定量表示，其定義為含水合物沉積物的滲透率與不含水合物的沉積物初始滲透率之比.

含水合物沉積物滲透率的降低，是由沉積物孔隙內水合物的生長改變了孔隙形狀、大小及連通性導致的.水合物呈膜狀附著于孔隙壁面時，對沉積物滲透率影響較小；當水合物生長于孔隙中心時，對沉積物滲透率影響較大，且可能造成孔隙堵塞.目前國內外學者建立了很多含水合物沉積物滲透率與水合物飽和度的理論模型，這些滲透率模型在水合物藏開發模擬中發揮了重要作用.其中最簡化的模型是將含水合物沉積物的孔隙假設為理想的圓管或者球形顆粒堆積形成的粒間孔，然后根據水合物在孔隙中的分布模式來建立的函數關系式，通常主要認為有顆粒包裹型和孔隙填充型(圖1).

圖1 沉積物中水合物分布模式示意圖Fig.1 Schematic diagrams of natural gas hydrate distribution in sediments

1.1 平行毛細管模型

假設沉積物孔隙空間為等徑平行毛細管，當毛管半徑為a，長度為L，毛管數量為n時，流體總流量表示為

式中，μ為流體黏度，?p/L為壓力梯度.考慮孔隙度φ 與單位橫截面積毛管數之間的關系(φ=nπa2),不含水合物的多孔介質的滲透率為

Kleinberg 等[35]在此基礎上，分別推導了水合物生長于孔隙壁面和孔隙中心兩種情況的含水合物沉積物滲透率模型.當假設水合物均勻生長于孔隙壁面時，毛管半徑減小至ar，則含水合物沉積物滲透率為

式中，毛管半徑與水合物飽和度Sh的關系為

故水合物飽和度為Sh的沉積物滲透率模型為

此時含水合物沉積物歸一化滲透率為

當生成的水合物占據在毛細管中心呈圓柱狀分布時，孔隙中的流動空間變為環狀，其中水合物半徑為b，單根毛管中的流量公式為[36]

考慮水合物飽和度Sh與a，b的關系于是該類型含水合物沉積物滲透率為

故歸一化滲透率為

Masuda[37]推廣了水合物生長于毛管壁面的平行毛細管模型，提出含水合物沉積物滲透率與水合物飽和度的指數關系模型(常被稱為Tokyo 模型)

式中指數N是變化的，當考慮水合物在孔喉處生長時，應適當增大N的值.在其數值模擬中指數N定義為10 或15，但具體原因不明確[38].

1.2 Kozeny 顆粒模型

實際含水合物沉積物骨架是由顆粒組成，由于其孔隙空間的不規則性，導致流體實際流動路徑La比定義壓力梯度的直線距離L更長.因此顆粒介質的滲透率預測要比簡單管道介質的滲透率預測要困難得多.Kleinberg 等[35]基于Kozeny-Carman(KC)方程率先提出了Kozeny 顆粒模型，模型中做了如下假設：(1)沉積物的電迂曲度可以代替水力迂曲度；(2)形狀因子不隨水合物飽和度變化；(3)沉積物孔隙可視為毛細管.推導過程基于的KC 方程表示為

式中，υ 為形狀因子，A為孔隙表面積，Vp為孔隙體積，τ 為迂曲度，τ=(La/L)2.

Hearst[39]指出迂曲度與地層因子和孔隙度之間的關系為τ=Fφ.Spangenberg[40]詳細討論了不同水合物生長模式下的電傳導問題，水合物飽和度為Sh的多孔介質地層因子F(Sh)與飽和水的多孔介質地層因子F0間存在以下關系

式中，n為阿爾奇飽和度指數.故含水合物飽和度為Sh時沉積物的歸一化滲透率為

式中，A0與A(Sh)分別為不含水合物和含水合物的沉積物孔隙表面積.

Kleinberg 等[35]基于水合物生成于沉積物顆粒壁面和孔隙中心兩種情況，推導了含水合物沉積物中水的歸一化滲透率模型.若生成的水合物包裹在沉積物顆粒上，假設沉積物孔隙為圓管，初始被水填充的孔隙表面積隨著水合物飽和度的增大而減小.若不含水合物時孔隙半徑為a，水合物飽和度為Sh時孔隙半徑為ar，故沉積物歸一化滲透率為

對于該模型，當00.8 時，指數n變化較大，但此時滲透率已經很小，指數n的變化對其只有很小的影響[40].

若生成的水合物占據孔隙中心，孔隙表面積隨著水合物飽和度的增大而增大，其與水合物飽和度具有如下關系

故含水合物沉積物的歸一化滲透率變為

式中飽和度指數n=0.7Sh+0.3[40].

Dai 和Seol[41]基于數值模擬結果，提出了KC 方程修正模型來描述含水合物沉積物歸一化滲透率與水合物飽和度之間的關系

上述模擬過程中考慮了沉積物孔隙的非均質性以及水合物生長的隨機性，該模型不對水合物生長模式進行假設，綜合了水合物生長引起的孔隙空間和迂曲度的變化，不含其它經驗參數，具有一定的普適性，與實驗結果吻合較好.

Katagiri 等[42]基于KC 方程，分別假設沉積物顆粒為圓柱形立方堆積、球形立方堆積和球形隨機堆積3 種方式，推導了水合物生長于顆粒壁面和孔隙中心時的滲透率模型.對于顆粒為圓柱形立方堆積的沉積物，當水合物生長于壁面時，歸一化滲透率為

當水合物生長于孔隙中心時，含水合物沉積物歸一化滲透率為

考慮到隨飽和度的增大，水合物會與沉積物顆粒發生重疊，因此上式僅適用于時，即Sh<0.628.

對于顆粒為球形立方堆積的沉積物，當水合物生長于壁面時，歸一化滲透率為

當水合物生長于孔隙中心時，含水合物沉積物歸一化滲透率為

對于顆粒為球形隨機堆積的沉積物，當水合物生長于壁面時，歸一化滲透率為

式中，e為孔隙比例，與孔隙度φ 的關系為e=φ/(1?φ).當水合物生長于孔隙中心時，含水合物沉積物歸一化滲透率為

式中，Cp為沉積物顆粒半徑與水合物半徑之比.

1.3 其他模型

美國勞倫斯伯克利國家實驗室開發的EOSHYDR/TOUGH2 模擬器[43]應用非飽和土van Genuchten(vG)模型[44]作為含水合物沉積物歸一化滲透率模型

Delli 和Grozic[46]綜合考慮含水合物沉積物中水合物的兩種生長模式，基于權重推導了綜合的歸一化滲透率模型

Singh 等[30]結合含水合物沉積物中氣水滲流的Navier-Stokes 方程穩態形式與動量守恒方程，提出了含水合物多孔介質的氣水兩相滲透率解析模型.顆粒包裹型含水合物沉積物氣水相滲模型為[30]

孔隙填充型含水合物沉積物氣水相滲模型(kri)PF形式過于復雜，此處不作展示，詳細內容參考文獻[30].

當綜合考慮兩種水合物生長模式時，含水合物沉積物相滲模型為

式(26)～式(28)中，R為孔隙半徑，rg和rw分別為孔隙內供氣體和水流動空間的半徑，α 為氣水黏度比，τ(φ,1-Sh?Swr)為孔隙管道內可流動空間的迂曲度，τ(φ,Sw)與τ(φ,Sg)分別為孔隙內水和氣流動路徑的迂曲度，β 為孔隙形狀校正系數.該模型綜合考慮孔隙形狀、平均孔徑、孔隙度、束縛水飽和度、水合物飽和度等巖石特性，氣、水飽和度和黏度等流體特性以及水合物生長模式.可利用該模型分析孔隙形狀、孔徑、潤濕性和水合物生長模式等對氣水滲流的影響，反過來也可以根據氣水相滲曲線評價巖石物性.

考慮毛管力對氣水兩相滲流的影響，Singh 等[50]推導了含水合物沉積物氣水兩相滲流模型，該模型中引入了4 個經驗參數(氣、水每相各兩個參數)，這4 個參數只需要根據任一給定的水合物飽和度下的實驗參數擬合求取一次，便可在其他水合物飽和度下預測含水合物沉積物氣水滲透率

式中，kri為氣體或水的相對滲透率，kpm為水合物飽和度為Sh時沉積物的絕對滲透率，βi與ηi為經驗參數.

2.隨著人均可支配收入水平的提高，居民更多追求高端的非基本公共服務消費需求。通過實證可以看出人均可支配收入對三者的影響系數分別為0.453、0.245、0.638，且在1%的水平上顯著，說明人均可支配收入對公共服務均等化指數具有正向影響，使得均等化程度降低。

表1 對上述主要模型做了簡要的總結，分析了各模型的優缺點.

2 數值模擬研究

國內外學者已將孔隙網絡、格子玻爾茲曼、分子動力學等數值方法應用于含水合物沉積物滲流模擬研究，不同的模擬方法適用于不同研究尺度及條件.孔隙網絡模擬方法致力于孔隙尺度下構建復雜多孔介質孔隙結構，結合CT 掃描等實驗手段獲取巖心孔隙網絡以模擬多相流體流動，進一步對多孔介質滲流性質進行預測；格子玻爾茲曼方法是針對介觀尺度下廣泛使用的動力學方法，可以描述含水合物多孔介質的內部流場以表征多相流體滲透率變化，研究多孔介質中多相流動規律.分子動力學模擬方法針對微觀尺度下天然氣水合物成核、生長及分解過程，從微觀分子的角度模擬開采過程的產氣效率.除此之外，通過求解控制方程的有限元與離散元等數值方法，也被用于含水合物沉積物滲透率與飽和度等參數間規律關系的研究.

表1 含水合物沉積物歸一化滲透率模型對比Table 1 Comparison of normalized permeability models of hydrate sediments

2.1 孔隙網絡模擬

Liang 等[51]采用孔隙網絡模擬方法研究水合物形成過程中的滲透率變化.模擬過程針對兩個關聯問題展開：一是考慮水合物形成過程中相平衡條件的轉化與孔隙分布參數之間的關系；二是考慮水合物生長習慣對于滲透率變化的影響.Dai 和Seol[41]從孔隙尺度水合物生長模式與介觀尺度非均質性兩方面，對不同含水飽和度下水相滲透率和本征滲透率的比值與水合物飽和度之間存在的規律關系進行研究，利用孔隙網絡模型對水合物分解過程中迂曲度和比表面積變化進行表征，通過修正KC 模型提出孔隙度變化模型.Mahabadi 和Jang[52]采用孔隙網絡模型模擬氣體膨脹過程計算氣水相對滲透率.基于數值模擬得到的氣水相滲結果提出了修正Stone方程的擬合參數.車雯等[53]基于孔隙網絡模型模擬水合物生成于孔隙壁面與中心兩種方式下多孔介質滲流特性變化，結果表明水合物生成于孔隙中心時絕對滲透率小于生成于壁面時的值，水合物飽和度相同時多孔介質孔徑越大滲透率越大，水合物生成于孔隙中心時兩相相對滲透率等滲點低于生成于壁面的情況，當水合物飽和度變化時兩相相對滲透率幾乎不變.Mahabadidi 等[54]采用CT 與孔隙網絡模型相結合的方法，從水合物沉積物CT 圖像中提取孔隙網絡，進一步模擬氣體侵入、水合物分解、氣體膨脹過程以及氣水滲透率變化，達到確定相滲模型中擬合參數取值的目的.Wang 等[55]利用CT 結合孔隙網絡模擬的方法擴展了顆粒粒度和孔隙度對滲透率影響的研究.結果表明顆粒粒徑越大形成的水合物多孔介質孔隙率越大，毛管壓力越小，絕對滲透率越大.在相同的含水飽和度下，孔隙度越大，水相相對滲透率越大，氣相相對滲透率越小.Wang 等[56]為探明水合物沉積層骨架結構特性及內部流體特性等因素對氣水兩相滲流影響規律，將孔隙網絡模型引入到水合物沉積層滲流特性研究，建立了適用于天然氣水合物滲流特性模擬研究的孔隙網絡模型.Ai 等[57]基于多孔介質中孔隙結構的變化特征以及氣水流動的相互制約特性，以最大球算法為基礎建立了適用于氣水兩相流動的孔隙網絡模型，進而計算分析水合物分解過程中氣水的流動以及產氣產水受到的影響.

2.2 格子玻爾茲曼模擬

喻西崇等[58]根據沉積物中水合物分解過程中流體運移和孔隙介質的特點，提出將格子玻爾茲曼方法(lattice Boltzmann method,LBM)用于含天然氣水合物沉積物中多相滲流規律研究，使用LBM 對多孔介質中水合物生成、分解過程中飽和度變化對滲透率特性的影響進行模擬，結果表明多孔介質中單相流動的流場分布與孔隙直徑及滲透率有關，水合物生成使得多孔介質滲透率大大降低.Kang 等[59]提出考慮毛細效應的含水合物沉積物滲透率模型，針對顆粒包裹型和孔隙填充型兩類水合物賦存形式開展LBM 流動模擬，研究水合物形成過程中水合物飽和度與滲透率之間的關系.Hou 等[60]認為孔隙網絡模型應用于流動模擬時，對多孔介質的簡化會丟失水合物分布形態的信息，LBM 在流動模擬研究方面是一種有效的手段，并利用LBM模擬研究水合物對沉積物滲流特性的影響，討論了礦物顆粒與水合物在孔隙尺度的分布形式對滲透率變化的影響.Chen 等[61]在氣體過量條件下合成氙氣水合物，利用CT 掃描圖像獲取時間與空間上孔隙尺度水合物分布形態，結合LBM 方法計算氣相相對滲透率，獲取水合物飽和度與滲透率之間的規律關系.Zhang 等[62]建立了一種基于LBM 的多孔介質中反應輸運與非等溫多物理化學過程的孔隙尺度框架，較好地描述了不同水合物生長模式下甲烷水合物分解相態變化和孔隙結構演化.結果表明，孔隙尺度甲烷水合物分解過程數值模擬有助于加深對連續介質模型滲透率-飽和度變化關系的理解.

2.3 分子動力學模擬

分子動力學方法在分子水平模擬水合物合成分解過程中應用廣泛.Fang 等[63]利用分子動力學方法研究多孔介質沉積物中水合物分解動力學過程.Yasuoka 等[64]對一定數目水分子與甲烷分子進行分子動力學模擬觀察水合物分解過程，以揭示水合物分解過程中的結構類型.Hirai 等[65]使用分子動力學方法研究二氧化碳水合物形成過程.萬麗華等[66]通過分子動力學模擬在乙二醇溶液作用下甲烷水合物分解的微觀過程，從原子水平上研究水合物分解的微觀過程.Moon 等[67]提出甲烷水合物在合成條件下的分子動力學模擬，模擬過程考慮甲烷與水接觸面，避免以往研究中不合理的高濃度甲烷設置和傳質限制.Kondori 等[68]綜述了水合物不同開采方式的關鍵理論與不足之處，介紹了分子動力學模擬理論在各開采方法中的應用，進一步闡述了分子動力學方法在正確理解水合物分子結構與分解過程中的實際意義.

除上述方法外，Kossel 等[69]在甲烷水合物合成實驗中結合核磁共振成像技術獲得天然氣水合物三維分布圖，然后把圖像輸入到有限元算法程序(COMSOL Multiphysics)測試對比5 種滲透率方程，得到水合物飽和度與滲透率之間關系.Katagiri 等[42]基于KC 方程提出滲透率模型，使用離散元方法構建不同水合物生長模式下立方堆積結構，結合計算流體力學模擬對模型假設進行驗證.

3 多相滲流實驗

實驗室內進行含水合物沉積物滲流實驗難度較大，目前實驗方法主要有利用特殊裝置對含水合物沉積物保壓巖心進行滲流實驗，測定含水合物巖心滲透率，或者在含人工合成水合物的沉積物中進行單相或氣水兩相滲流實驗，研究含水合物沉積物滲流特性，但十分耗費時間，而且需要精確控制壓力和溫度[70-71].沉積物中水合物飽和度的確定是實驗中的另一難點，目前主要是利用CT掃描成像技術或根據水合物分解釋放的氣體量計算[72].不同介質中水合物飽和度測試結果表明，水合物飽和度隨著沉積物顆粒粒度的增大而增大，黏土中水合物飽和度為0～10%，粉砂黏土中水合物飽和度約為45%，富含砂的沉積物中水合物飽和度為65%～85%，而在礫石堆積的沉積物中，飽和度可達到100%[73-74].目前國內外進行的含水合物沉積物滲流實驗主要包括有效滲透率測試、滲流過程中沉積物物性演變表征以及水合物開采室內模擬等方面.

3.1 有效滲透率測試

含水合物沉積物的有效滲透率是對水合物的開發具有重要意義的關鍵參數，其主要與沉積物本身的孔隙結構、水合物飽和度以及水合物的生長模式等因素有關，含水合物沉積物的有效滲透率測試一直是國內外學者研究的熱點問題之一，對此已有大量相關實驗研究.

Yoneda 等[75]對來自印度海域的保壓取芯樣品在PCATS 設備中完成密度、孔隙度、水合物飽和度及CT 掃描等測試后，進行了原位含水合物沉積物樣品初始滲透率測定.實驗發現，受沉積物粒徑及水合物飽和度影響，含水合物沉積物滲透率分布在0.01 毫達西到幾十毫達西之間.Minagawa 等[76]在砂中人工合成了水合物并進行單相水滲流實驗獲得了滲透率數據，同時利用核磁共振技術分析了含水合物沉積物的孔隙分布，并根據前人模型計算了理論滲透率，獲得的理論滲透率值與實測滲透率值保持一致.Konno 等[77]在設計的鋁合金夾持器中進行了含水合物沉積物滲流實驗，分別進行了注氣和注水滲流實驗，滲流過程中測量樣品兩端壓差及流量計算各相滲透率，并根據質量守恒定律計算了沉積物中各相飽和度，同時利用CT 掃描技術檢測了注水滲流過程中驅替前緣的移動過程.對于注氣滲流實驗，最終沉積物中水合物、水和氣體的占比分別為69%，5%和26%，根據滲流壓差及氣體流量計算得到含水合物沉積物的滲透率為0.417 mD；注水滲流實驗最終沉積物中水合物、水和氣體的占比分別為76%，21%和0，計算得到含水合物沉積物的滲透率為0.423 mD.Delli 和Grozic[70]在粒徑為20～30 目的砂中人工合成了二氧化碳水合物.他們先使樣品飽和蒸餾水，然后使用二氧化碳驅替樣品至流出預定的水量，從而使砂具有既定的初始含水飽和度.水合物合成后，使用氦氣驅出樣品中殘余的二氧化碳，而后使用蒸餾水測量含水合物沉積物的滲透率.實驗測試得到的滲透率數據與已存在滲透率模型對比結果表明，隨著水合物飽和度的增大，水合物在沉積物孔隙中的生長模式趨向于由顆粒包裹型向孔隙填充型轉變.Marinakis 等[78]對取自阿姆斯特丹海底火山泥沉積物進行蒸餾水洗鹽、烘干以及重塑，使用蒸餾水和甲烷、乙烷與丙烷的混合氣體合成水合物，對含水合物沉積物進行恒壓水驅滲流實驗，實驗結果量化了水合物飽和度對含水合物沉積物滲透率的影響，發現水合物分解前后沉積物的滲透率保持在同一數量級，即使對于極高的水合物飽和度，滲透率也只是增加三倍.Konno 等[79]對日本南海海槽東部的含水合物與不含水合物的沉積物樣品進行滲流實驗，以人工配置的鹽水為測試流體，采用恒定流速法測量樣品滲透率.劉樂樂等[80]和張宏源等[81]運用瞬態壓力脈沖法對含水合物沉積物進行滲透性測試，實驗獲取的沉積物滲透率與數值模擬結果一致.Johnson 等[71]在取自美國阿拉斯加北部陸坡埃爾伯山脈的松散沉積物中利用甲烷和鹵水人工合成了甲烷水合物，在含水合物沉積物中分別進行了鹵水單相滲流實驗以及非穩態法氣水兩相滲流實驗.單相滲流實驗表明沉積物的滲透性隨著水合物飽和度的增大而明顯降低.在水合物飽和度為12.5%時，氣水兩相滲流實驗發現含水合物沉積物中的兩相流動為活塞式驅替.沉積物末端有氣體流出時，水的流量很小，沉積物中的可動水已基本被驅出.由于含水合物沉積物滲透率很低，為防止測試中生成額外的水合物，該氣水兩相滲流中測試氣體為氮氣.水合物合成過程中甲烷是充足的，不再向沉積物中補充水，當滿足合成的條件時水合物將會一直合成.然而實際上，甲烷氣體難以全部進入沉積物孔隙中與水形成水合物，同時水合物的形成也會限制氣體流動，因此水合物的合成會逐漸減慢并最終停止.

目前為止，國內外對于某一水合物飽和度下沉積物的單相滲透率的測試實驗發展迅速但尚未成熟，在沉積物中精確控制人工合成水合物飽和度難度較大，測試過程中是否發生水合物的分解與二次合成以及對測試滲透率的影響尚不明晰.此外，含水合物沉積物中氣水兩相滲流實驗難以開展，獲取氣水兩相滲流曲線困難.

3.2 滲流過程中沉積物物性演變研究

雖然海底大多具有適合水合物生成的低溫高壓條件，但實際水合物一般存在于靠近大陸邊緣的沉積物中[82].該種沉積物由于埋深較淺，多為未膠結成巖、富含黏土的泥質粉砂，其結構容易在水合物開采過程中發生改變，從而造成沉積物物性演變，影響水合物的開采效果.

Konno 等[79]測試了含水合物沉積物的有效滲透率及不含水合物沉積物的絕對滲透率，發現含水合物沉積物的滲透率與有效孔隙度具有較強的相關性.滲流過程中對沉積物加載有效應力并改用測試流體為蒸餾水，在10 MPa 應力下砂沉積物的滲透率降低為原來的一半，粉砂質和黏土質沉積物的滲透率受蒸餾水的影響明顯.分析認為水合物降壓開采過程中地層有效應力改變，使未固結的松散沉積物結構變得致密，同時水合物分解產生的水會改變沉積物孔隙中地層水的鹽度，可能會導致沉積物中蒙脫石等黏土礦物膨脹，堵塞沉積物孔隙，從而降低沉積物滲透率.Jarrar 等[83]在包含高嶺石粉末的石英砂中進行氣水滲流實驗，模擬水合物開采過程沉積物中的滲流行為，利用同步輻射計算機斷層掃描成像手段觀察到樣品顆粒運移現象和裂縫的形成.Lu 等[84]對中國南海神狐海域泥質粉砂沉積物進行了注水滲流實驗，實驗中隨著樣品兩端滲流壓差的增大，滲透率逐漸減小，且該過程是不可逆的.實驗中利用CT 掃描觀察到滲流過程中樣品中孔隙體積蠕變的現象，滲流過程不同壓差下沉積物的有效滲透率與孔隙度存在明顯的冪律關系.值得注意的是，滲流壓力梯度小于每米3 MPa 時，沉積物孔隙變形較小；壓力梯度大于每米3 MPa 時，沉積物孔隙空間迅速壓縮，滲透率迅速下降.Yoneda 等[75]對取自印度海域的含水合物沉積物樣品進行了滲透率測試，對比前人研究，認為該地區水合物主要生長于沉積物顆粒表面.水合物的分解會導致沉積物滲透率變大，在樣品上加載應力時，滲透率逐漸減小，最終可能會小于含水合物時的初始滲透率.

針對水合物沉積物物性演變特性，國內外學者開展了大量實驗，研究了沉積物結構及物性演化現象，分析了對氣水滲流的影響，對解釋實際生產中出現的問題提供了參考，也為進一步改進水合物開采技術，提高開采效率提供了新思路.

3.3 水合物開采室內模擬

Konno 等[85]通過沉積物中水合物的合成與分解實驗，認為沉積物性質直接影響著水合物降壓分解速度及產氣量.李淑霞等[86]通過多孔介質中水合物降壓分解實驗，對影響降壓開采過程的因素進行了分析.李淑霞等[87]采用自制的一維填砂管水合物開采模擬實驗裝置，在石英砂中人工合成了不同飽和度的甲烷水合物并進行降壓開采模擬實驗.實驗結果表明整個開采過程可分為自由氣產出、水合物降壓分解產氣和最后產出已分解氣3 個階段.初始水合物飽和度越大，降壓分解過程中溫度下降幅度越大.實驗范圍內中等飽和度的水合物藏降壓分解產氣速率較大，開采效果最好.Lee 等[88]進行水合物降壓分解實驗，設置不同的降壓幅度，對多孔介質中水合物進行降壓分解實驗，分析了降壓幅度對水合物分解過程的影響.認為水合物的分解受分解壓差控制，壓差越大，水合物分解速率越慢.此外，降壓級次也是影響產氣速率和產氣量的重要因素，越靠近相平衡壓力時產氣速率越快.實驗結果表明實際開采過程中可通過多級降壓使水合物分解壓力處于相平衡壓力附近，從而加快產氣速率.考慮到水合物分解是吸熱過程，Konno 等[89]通過水合物降壓分解實驗，發現降壓開采過程中還應適當補充熱量，否則可能導致水合物分解產生的水迅速在沉積物孔隙中結冰，堵塞氣體滲流通道，從而影響開采效率.Sakamoto 等[90]通過實驗向含水合物多孔介質中注熱水模擬水合物開采，根據分解過程中的產氣情況將水合物注熱開采模擬實驗過程分為熱水驅替沉積物自由氣、溫度上升水合物開始分解、溫度達到平衡水合物穩定分解、分解結束等4 個階段.李淑霞等[91]在自制的一維填砂管模型中進行了注熱鹽水開采水合物實驗，根據實驗結果分析認為影響能量效率的因素由大到小依次為：注熱水溫度、水合物飽和度、初始溫度、注熱水時間、注熱水速度，提出在注熱開采時需合理優化注熱水溫度、時間和速度，從而實現水合物的高效開發.Wang 等[92]在實驗中模擬了降壓、注熱以及兩者相結合的水合物開采方式，認為降壓注熱相結合是水合物開采的最佳方法.

4 產氣模擬研究

目前國內外使用的主流水合物數值模擬軟件工具包括美國勞倫斯伯克利國家實驗室開發的TOUGH+HYDRATE 可用于模擬天然氣水合物多組分多相態流體流動，具有較強優勢，應用廣泛；日本甲烷水合物資源聯盟開發的MH21-HYDRES 可以模擬熱激發、降壓法以及各種方法相互結合的開發過程，適用于4 種組分、5 種相態；加拿大CMG公司(Computer Modelling Group LTD)開發的CMGSTAR 可用于多組分流體流動模擬；美國西北太平洋國家實驗室開發的STOMP-HYD 可用于模擬地下多相流體的運動和變化規律.

Moridis[93]提出TOUGH2 多組分多相流、熱流和地下輸運模擬軟件的新模塊EOSHYDR，通過求解質量與熱量守恒耦合方程來模擬常規甲烷水合物儲層條件下非等溫甲烷釋放、相態和流動特征.Wilder 等[94]比較了5 種儲層模擬軟件在非等溫多流體平衡態問題、封閉域水合物分解、一維水合物分解、一維徑向水合物分解和二維徑向水合物分解5 類問題上的表現.Zhang 等[95]提出一種基于粗粒度并行計算的水合物模擬的區域分解方法，該方法將模擬區域細化，模型域由離散子域組成，適用多處理器同時模擬.這種區域分解方法補充到TOUGH+HYDRATE 代碼中，速度與拓展性方面有所提升.Moridis[96]發展并更新了水合物儲層系統模擬代碼TOUGH+HYDRATE，通過求解耦合質量與熱量平衡方程，TOUGH+HYDRATE 可以模擬天然甲烷水合物儲層條件下復雜地質介質非等溫氣體釋放、相行為、流體與熱流動等過程.

4.1 TOUGH+HYDRATE

Reagan 和Moridis[97]利用并行TOUGH+HYDRATE 代碼模擬了西斯匹次卑爾根島邊緣條件下淺層水合物分解過程，以驗證觀測到的氣體釋放來自水合物的假設.模擬結果與假設一致，也與觀測結果具有顯著的一致性.李剛等[98]以2007 年南海北部神狐海域SH2、SH3 和SH7 站位的鉆探測井資料為基礎，利用TOUGH+HYDRATE 建立實際水合物藏分層地質模型，通過設計不同開采井的方式，進行了單井降壓和降壓注熱開采模擬.李小森等[99]利用TOUGH+HYDRATE軟件模擬深海條件下天然氣水合物藏單井降壓的開采過程與甲烷氣體擴散過程.Reagan 等[100]利用并行TOUGH+HYDRATE 代碼評估大尺度深海水合物儲層的生產潛力以制定有效的生產策略.模擬過程涵蓋38 km2的三維系統，考察水平井、單井和多井的產能，并探索儲層性質的變化.胡立堂等[101]針對南海神狐海域天然氣水合物的成藏條件，利用并行TOUGH+HYDRATE 軟件結合水平井開采方法對天然氣水合物注熱和注熱降壓開采進行數值模擬，模擬結果表明注熱降壓開采模式比單純注熱模式效果有較大改善，溫度對于提高甲烷產量效果不明顯.Sun 等[102]基于日本南海海槽東部第一海上生產井位可用地質數據，利用TOUGH+HYDRATE 模擬軟件建立與生產測試方案匹配的模擬模型，包括水合物形成與分解的平衡與動力學模型，考慮4 相(氣、液體、冰和水合物)與4 種質量組分(水、甲烷氣、水合物和水溶性抑制劑)，研究長時段生產表現以及射孔段對產氣的影響.驗證結果發現，模擬產氣量與實測產氣量基本一致，但產水率存在偏差.Sun 等[103]基于2015 年中國第三次天然氣水合物鉆探GMGS3-W19 井位地質資料，結合TOUGH+HYDRATE 與商業代碼FLAC3D 對南海該井位儲層進行儲層產量評價及穩定性模擬研究.Yu 等[104]基于2017 年日本南海海槽甲烷水合物生產資料信息，利用TOUGH+HYDRATE 模擬軟件從多層儲層模型中對甲烷水合物生產進行短期模擬及模型驗證，對一定生產周期的長期產氣表現進行預測，研究生產間隔對水合物分解產氣過程的影響.

4.2 MH21-HYDRES

Masuda等[105]利用MH21-HYDRES 方法針對甲烷水合物井使用甲醇吞吐法模擬產氣量，并與降壓方法比較，分析甲烷水合物井吞吐增產和連續產氣過程中水合物的分解過程，模擬結果表明增產措施使井筒附近水合物分解速度顯著增大.Kurihara 等[106]基于日本南海海槽東部甲烷水合物儲層測井資料解釋與三維地震數據分析結果，構建三維儲層模型利用MH21-HYDRES 數值模擬程序進行生產預測.根據數值模擬結果討論不同尺度儲層非均質性、滲透率變化以及井位對生產的影響.Konno 等[107]為研究實際含水合物沉積物滲透率與降壓法水合物分解過程，利用MH21-HYDRES 模擬水合物分解過程中累積產氣、產水和壓力變化，并結合南海海槽東部砂質濁積巖沉積物中水合物分解實驗，證實降壓法的優越性以及數值模擬方法對分析水合物分解行為的有效性，可用于儲層產氣量的評價.

4.3 STOMP-HYD

White 和McGrail[108]證實了實驗室尺度上直接交換技術對地質介質產氣與二氧化碳封存的可行性.采用STOMP-HYD 數值模擬的方法，研究了地質沉積物中二氧化碳和甲烷的降壓、熱刺激和直接分子交換聯合作用下天然氣水合物的開采.

4.4 CMG STAR

Zatsepina 等[109]利用CMG STAR 數值模擬方法研究III型水合物氣藏產氣量，結果表明水合物分解能力不僅取決于傳熱速率，還取決于地層允許流體流動的能力，并探索了流體流動與傳熱之間的相互作用，討論了III 型水合物儲層數值模擬面臨的挑戰.Myshakin 等[110]對墨西哥灣沃克嶺水合物儲層利用TOUGH+HYDRATE 與CMG STAR 分別進行了降壓法產氣預測的二維和三維數值模擬.李淑霞等[111]基于中國南海神狐海域水合物藏地質資料，采用CMG STAR 模擬軟件研究降壓開采過程分解前緣的形狀和移動規律，并對影響分解前緣移動的敏感性因素進行分析.

4.5 其他模擬方法

沈海超[112]針對天然氣水合物藏降壓開采相關問題開展流固耦合數值模擬研究，建立了天然氣水合物藏降壓開采相變滲流模型，根據流固耦合滲流理論，建立了考慮相態變化的滲透率各向異性天然氣水合物藏降壓開采氣水兩相非等溫流固耦合數學模型.對流固耦合模型控制方程進行離散，基于有限元平臺(Finite Element Program Generator)，開發了天然氣水合物藏降壓開采氣水兩相非等溫流固耦合數值模擬軟件.程遠方等[113]基于多孔介質中水合物分解動力學、傳熱及氣水兩相流理論，同時考慮水合物分解引起的滲透率及有效孔隙度變化，建立水合物、氣、水三相的天然氣水合物降壓分解模型，并進行有限元程序開發及驗證.于峰等[114]在考慮水合物二次生成及滲透率變化的基礎上，建立了實驗室尺度下水合物降壓開采數學模型，利用該模型對水合物降壓開采過程中水合物二次生成進行有限差分模擬，分析評價水合物二次生成及滲透率變化對水合物分解產氣的影響.Liu 等[115]開發了GrapeFloater 數值模擬器，實現了熱對流-擴散、兩相流動、水合物相變、沉積物顆粒變形的多場耦合，研究了水合物降壓開采過程中沉積物變形特征以及氣體生產情況.

雖然天然氣水合物試采工作已初步進行，但數值模擬方法仍然是目前的主要研究手段，針對開采方法的不同，模擬軟件展示出不同的適用性，對于降壓開采、注熱開采方法來說，TOUGH+HYDRATE 適用性較好；STOMP-HYD 對于二氧化碳置換法模擬應用廣泛.另外，模擬過程與實際現場開采對比，考慮的影響因素會有減少或簡化，需要增加模擬過程的影響因素敏感性分析，使模擬結果更加具有指導意義.

5 結論與展望

天然氣水合物的開采涉及傳熱、水合物相變、多相滲流、地層變形等物理化學效應，其中沉積物內多相滲流行為影響其對流傳熱效率及孔隙壓力消散速度，對水合物開采過程的重要性不言而喻.本文著重介紹了含水合物沉積物基礎滲流理論、滲流數值模擬方法、多相滲流實驗以及產氣模擬研究等方面的內容，以期為含水合物沉積物多相滲流特性研究以及水合物開采工程提供一定的參考和借鑒.對于含水合物沉積物滲流特性的研究，未來可側重于以下幾個方面：

(1)含水合物沉積物多相滲流理論的深入研究在今后仍然是一個重要方向.目前絕大多數含水合物滲流理論的研究依然借鑒沿用達西定律，眾多已存在滲透率模型大多包含擬合參數，或在沉積物孔隙結構以及水合物生長模式等方面作理想化假設，在實際應用中局限性較大.為進一步豐富和完善含水合物沉積物多相滲流理論，還需從含水合物沉積物孔隙結構特征出發，深入分析水合物對滲流行為的影響，建立更具有普適性的滲流模型.

(2)含水合物沉積物滲透率的測定能為水合物開采數值模擬以及實際開采工程提供基礎資料與指導.目前已發展了諸多含水合物沉積物滲透率測試實驗方法，但實驗所用沉積物樣品取自實際水合物沉積物的較少，獲取的實驗結果對于水合物實際開采的指導意義有限.此外，如何縮短含水合物沉積物滲透率測試時間、減少測試過程對水合物的擾動以及開展多尺度氣水兩相滲流實驗也是未來需要研究的重要方向.

(3)天然氣水合物多賦存于海底松散未固結沉積物中，受水合物開采過程中多種因素影響，沉積物極易發生結構變化從而導致物性演變，制約水合物的高效開采.但該研究目前大多還是處于定性描述階段，且對沉積物物性演化機理的認識尚不完善，在未來仍是一個極具挑戰性的課題.

(4)天然氣水合物開采具有成本高、難度大、周期長等特點，開展室內開采實驗模擬及數值模擬相對簡單，研究不同開采方式過程中沉積物滲流特性以及氣體產量，是研究天然氣水合物開采技術的基礎工作，對實際開采技術與方案的設計具有重要參考價值.將來對于室內實驗和數值模擬方面還應更加綜合考慮各方面因素，以求更貼近實際生產情況，增大模擬結果的可信度.

(5)現有各國的天然氣水合物試采或開發多使用的是直井開采，而多應用于低滲透儲層的水平井壓裂技術可有效提高生產井與水合物儲層的接觸面積，提高產氣效率.同時，水平井壓裂技術應用于海域天然氣水合物儲層開采還可有效減少平臺數目，增大儲層的動用程度.但水平井壓裂技術的應用不可避免的會帶來更加復雜的滲流問題，如天然氣水合物儲層壓裂過程中裂縫起裂擴展的力學及滲流力學機理，降壓開采過程孔隙-裂縫雙重介質系統中多尺度多相流體流動規律.因此，將水平井壓裂技術應用于天然氣水合物開采的研究是十分具有現實意義的前瞻性課題.