韋昌富 顏榮濤 田慧會 周家作 李文濤 馬田田 陳 盼

1.巖土力學與工程國家重點實驗室·中國科學院武漢巖土力學研究所 2.廣西巖土力學與工程重點實驗室·桂林理工大學

0 引言

天然氣水合物（以下簡稱水合物）可存在于海洋大陸架邊緣地層或寒區(qū)永久凍土層中[1]，本文主要討論海底水合物儲層的情況。海底水合物儲層不同于凍土儲層，除了環(huán)境溫度在0 ℃以上以及孔隙中不含冰外，所處水深大、埋深淺，骨架固結(jié)程度弱或沒有固結(jié)，工程地質(zhì)條件復雜。此外，海底水合物儲層不僅具有普通深海沉積物所具有的內(nèi)部結(jié)構(gòu)與工程地質(zhì)特性，而且還因水合物的可相變性，其內(nèi)部結(jié)構(gòu)及工程特性隨著水合物開采過程均在不斷地演變中?？上嘧冃灾冈诤线m溫度、壓力及水化學條件下，沉積物孔隙中固體水合物可以分解成液態(tài)水和天然氣；反之，孔隙中液態(tài)水和天然氣也可生成水合物。正是這種可相變性，使得含水合物沉積物在環(huán)境荷載作用下表現(xiàn)出比普通海底沉積物更為復雜的工程力學行為特征，并伴有顯著的多相、多場、多過程耦合效應。特別地，在水合物開采過程中，水合物分解產(chǎn)生天然氣和液態(tài)水并吸收周圍的熱量，從而引起滲流、擴散、熱傳導等耦合過程；同時，隨著水合物的分解，水合物的膠結(jié)作用逐漸受到破壞，沉積物骨架抗剪強度減弱，壓縮性增大，從而導致沉積物變形或骨架坍塌，最終引發(fā)流泥、流沙、井壁破壞、儲層沉降或失穩(wěn)、大范圍海底滑坡等重大工程或地質(zhì)災害問題[2-3]。因此，有效確定含水合物沉積物的工程力學特性并厘清其隨著水合物分解過程的變化規(guī)律，是實現(xiàn)水合物安全、高效開采的前提保證。

筆者圍繞著水合物開采中所涉及的關(guān)鍵土力學問題，從含水合物沉積物的物理力學特性的關(guān)鍵測試技術(shù)、含水合物沉積物的變形與強度特性、含水合物沉積物的物理力學模型、水合物開采過程的數(shù)值模擬方法等幾個方面，分析和評述國內(nèi)外相關(guān)研究現(xiàn)狀和進展，并在此基礎(chǔ)上探討這些方向的未來發(fā)展趨勢及所面臨的挑戰(zhàn)。從巖土工程的角度來看，水合物開采還涉及許多其他的重要方面，如原位取樣、原位土工測試、原位監(jiān)測、海底穩(wěn)定性分析、深?；A(chǔ)工程（如錨固基礎(chǔ)）等。毫無疑問，這些方面的工作也是不可或缺的，但本文除了在討論未來發(fā)展趨勢時有所涉及外，重點只放在含水合物儲層的工程力學特性方面。本文中，“含水合物沉積物”和“水合物儲層”這對名詞可以互用。在一些文獻中[4]，含水合物沉積物也被稱為“能源土”，但本文將采用國際比較流行的說法，即“含水合物沉積物”（hydratebearing sediment）[2]。

1 物理力學特性的測試方法和技術(shù)

1.1 含水合物沉積物試樣的制備方法

由于水合物穩(wěn)定條件較為苛刻，鉆取和運輸過程中壓力和溫度的變化會導致水合物分解，使得原狀含水合物沉積物試樣的獲取與運輸極為困難，因此，實現(xiàn)人工制備含水合物沉積物試樣顯得尤為重要。在此之前，了解水合物于沉積物中的賦存模式則是人工制備含水合物沉積物試樣的前提條件。

自然界中，天然氣水合物主要存在擴散型、滲漏型及復合型等3種成藏模式[5-6]。在擴散型成藏模式下，水量豐富，天然氣滲漏通量低，稱為富水環(huán)境，氣體只能通過擴散進入水合物穩(wěn)定帶與水結(jié)合形成水合物；而在滲漏型成藏模式下，構(gòu)造裂隙發(fā)育，天然氣滲漏量高，稱為富氣環(huán)境，下部游離氣體沿著構(gòu)造裂隙往上運移至水合物穩(wěn)定區(qū)和當?shù)氐目紫端Y(jié)合形成水合物[7-8]；復合型兼具擴散型和滲漏型兩者的成藏特征。因此，在這3種典型環(huán)境中，水合物在粗粒土中主要以4種模式存在：①孔隙填充模式。水合物顆粒主要懸浮于孔隙水中，與附近土顆粒無實質(zhì)接觸[9]。②持力體模式。這種模式由孔隙填充模式發(fā)展而來，當水合物含量達到一定量后部分水合物起支撐作用并橋接相鄰的顆粒，從而成為承重骨架的一部分。一般地，當水合物飽和度超過25%～40%時，填充型水合物會變成持力型水合物[9-12]。③膠結(jié)模式。水合物優(yōu)先在多孔介質(zhì)顆粒接觸處及附近生成，使土顆粒接觸并黏合在一起，從而膠結(jié)成一個固體骨架[13-14]。④包裹模式。水合物包裹于土顆粒周圍。在細粒沉積物中，水合物一般以分散的結(jié)核狀、透鏡體或脈狀塊體形式存在[15-18]。

為了模擬原位的水合物賦存情況可在實驗室制得含水合物沉積物試樣，目前實驗室制備試樣使用的主要有飽和法、非飽和法、冰種子法、混合法等。飽和法是指在一定壓力條件下制備氣體的超飽和溶液，將其注入沉積物后，通過降低溫度在沉積物孔隙中生成水合物[19]。這種方法可以很好地制得孔隙填充模式和持力體模式試樣，但由于甲烷等天然氣體的溶解度較低，該方法制備的水合物飽和度較低且耗時。為了克服這一困難，有學者嘗試利用添加劑來增大甲烷氣體溶解度[20-21]。非飽和法利用沉積物在非飽和情況下具有較高氣—水接觸面積的特點，使過量氣體與孔隙水在合適溫度、壓力條件下充分反應，從而在沉積物中生成較高飽和度的水合物，該方法已被證實能很好地模擬膠結(jié)模式和包裹模式的試樣[22-24]。冰種子方法是將冰粉與土混合制樣或?qū)⒃嚇宇A先凍結(jié)，在供氣條件下逐漸升溫，孔隙冰融化后與氣體結(jié)合形成水合物[25-26]。該方法可以加快水合物的生成過程，但在制樣過程中存在冰粉壓融現(xiàn)象，如果預先凍結(jié)試樣，則可能因為“凍吸力”引起水分遷移而造成水合物分布不均。上述3種方法均是先制備沉積物，再通過控制溫度和壓力條件在沉積物孔隙中形成水合物。為了進一步提高水合物制備效率及水合物飽和度，Hyodo等[27]通過先制備水合物粉，然后與沉積物按一定比例混合，壓制成含水合物沉積物試樣，這類方法稱為混合制樣法。由于水合物的形成條件較為苛刻，對試驗條件要求較高，應用這種方法的關(guān)鍵在于抑制制備過程中水合物粉的分解，避免水合物含量計算不準確。

總體上，雖然目前的實驗室制樣技術(shù)能夠初步滿足含水合物沉積物物理力學試驗要求，但仍有以下幾個關(guān)鍵問題亟待解決：①制備大尺度的孔隙填充模式和持力體模式的含水合物試樣，以供開展物理模型試驗；②克服制樣過程中的重力和吸力的影響[28]，制備物理力學試驗所需的均勻含水合物沉積物試樣；③在水合物形成過程中實現(xiàn)含量的精準控制與測量。

1.2 水合物含量的確定方法

含水合物沉積物的物理力學特性與水合物含量密切相關(guān)[12,23-24,27]，因此準確地確定水合物含量是研究含水合物沉積物物理力學特性的關(guān)鍵問題之一。目前常用的水合物含量的確定方法主要包括[29-31]：壓降法計算、含水量消耗計算、時域反射技術(shù)、核磁共振技術(shù)等。壓降法是指通過水合物生成前后氣壓變化計算出氣體消耗量從而計算得到水合物的生成量。在非飽和法制樣中，客體氣體分子一直保持過量，通過假定沉積物中水分完全參與水合反應生成水合物，即可利用所消耗的水量來計算出水合物飽和度。而時域反射（TDR）技術(shù)則是通過測定水合物生成前后孔隙水介電常數(shù)變化來反算水合物飽和度。近年來，田慧會等[31]利用低場核磁共振技術(shù)確定了天然氣水合物的水合數(shù)，進而求得了水合物含量大小。現(xiàn)有的各種水合物含量的確定方法均存在一定局限性，如測量精度難以確定、不同方法間缺乏一致性等，研發(fā)高精度的沉積物水合物含量測量方法仍是一項亟需開展的任務[32]。

1.3 微細觀組構(gòu)探測技術(shù)

現(xiàn)有的沉積物組構(gòu)的探測技術(shù)主要包括：①聲波技術(shù)[33-35]。即根據(jù)實測的聲波參數(shù)如波速、衰減及其頻散特性反演水合物沉積物的孔隙度、水合物飽和度、彈性模量等參數(shù)，并分析水合物在沉積物中的賦存狀態(tài)。②電阻率法[36-38]?；诔练e物電阻率與復電阻率頻散特性隨著水合物含量及賦存模式的變化規(guī)律，測定沉積物中水合物的含量與分布。③TDR技術(shù)[29,39]。根據(jù)介電常數(shù)與含水量的關(guān)系，推算出沉積物試樣中水合物含量的變化。④計算機輔助成像（CT）技術(shù)[40-42]。利用含水合物沉積物不同組分對X射線衍射的吸收與透過率不同而形成含水合物沉積物的微觀圖像，并由此觀測和表征水合物體積以及剪切過程中微細觀結(jié)構(gòu)變化。⑤核磁共振（NMR）技術(shù)[31,43-44]。利用孔隙水核磁信號與弛豫時間探測含水合物沉積物中水合物含量及賦存狀態(tài)。

在以上5種探測技術(shù)中，前3種探測技術(shù)主要根據(jù)含水合物沉積物的波速、電阻及介電常數(shù)隨沉積物各組分和組分含量的變化關(guān)系來反演其內(nèi)部微細觀組構(gòu)，此類辦法能定性反映體系中水合物相變過程以及微細觀結(jié)構(gòu)變化，但不適用于在微細觀層面進行定量探測和分析。與其他方法相比，CT技術(shù)和NMR技術(shù)是兩種極具潛力的組構(gòu)探測及參數(shù)測定技術(shù)，其最大的優(yōu)勢在于它們是直接、定量、快捷的測試方法。但是，X射線衍射存在輻射效應需要特殊材料進行試驗防護，NMR技術(shù)則必須保證相關(guān)試驗設備對含水合物沉積物無核磁信號干擾。另外，昂貴的價格也限制了CT技術(shù)和NMR技術(shù)在水合物研究領(lǐng)域的應用。

1.4 水力特性的測試方法

含水合物沉積物的水力特性包括滲透性和持水曲線。傳統(tǒng)上，沉積物滲透率的測定方法主要有非穩(wěn)態(tài)法和穩(wěn)態(tài)法。穩(wěn)態(tài)法是在穩(wěn)態(tài)滲透狀態(tài)下通過測定試樣兩端的壓力差及孔隙水的滲流量，再由達西定律求得滲透系數(shù)。在含水合物沉積物滲透系數(shù)的測定試驗中，穩(wěn)態(tài)法比較流行[45-47]。然而，穩(wěn)態(tài)法的測試時間比較長，在測試過程中水合物易發(fā)生分解，從而影響到測試精度。相比之下，非穩(wěn)態(tài)法所需時間較短，只需要測定水流量（或孔隙壓力）隨著時間變化關(guān)系，然后利用數(shù)學模型來反演得到滲透率。目前，非穩(wěn)態(tài)法已開始應用于含水合物沉積物滲透率的測定中，例如：Jaiswal[48]采用JBN非穩(wěn)態(tài)法測定了含水合物沉積物的有效滲透率；張宏源等[49]基于瞬態(tài)壓力脈沖法測定了含水合物粉細砂的滲透率。除了上述兩種方法，一些間接的滲透系數(shù)方法，例如，Kleinberg[50]利用T2分布特征采用Kenyon公式計算了含水合物沉積物中的相對滲透率，這種方法在石油領(lǐng)域測井中廣泛應用；Kneafsey等[51]基于CT技術(shù)監(jiān)測沉積物中不同組分的密度變化，推算出了含水合物砂土的滲透系數(shù)。總之，含水合物沉積物滲透性的測試難點在于水滲流過程會引起水合物的溶解或分解，減小水合物含量。因此，如何有效避免或減小水合物含量的變化是優(yōu)化含水合物沉積物滲透性測試的關(guān)鍵。

持水曲線是指沉積物中吸力與液態(tài)水含量的關(guān)系，其受到水合物含量的影響。Mahabadi等[52]在傳統(tǒng)的非飽和土持水曲線測試裝置上增加了控溫系統(tǒng)，測試了含THF水合物砂土在不同水合物飽和度條件下的持水特征曲線；Ghezzehei等[53]對壓力室進行了改造，使試樣處于高壓環(huán)境，并且提高溢出水端口壓力保證氣—水壓力差在高進氣值陶土板可控范圍內(nèi)，同時增加了控溫裝置來控制系統(tǒng)溫度，在此基礎(chǔ)上測得了含水合物沉積物的脫濕和吸濕曲線。由于持水特征曲線測試非常耗時，其中水滲出過程可能導致水合物分解，而施加的氣體壓力則有可能導致水合物進一步形成，從而影響到測試結(jié)果的精度。為了克服這一問題，有學者采用物理仿真模擬來研究水合物賦存對含水合物沉積物持水特征曲線的影響[54-55]。目前有關(guān)含水合物沉積物持水曲線的直接測試工作非常稀缺，主要問題在于如何通過控制溫度、壓力以及其他手段來保證測試過程中水合物不產(chǎn)生二次形成或分解，以保證水合物含量的恒定不變，從而提高測試結(jié)果的精度。

1.5 力學特性的測試技術(shù)

力學特性包括應力應變關(guān)系、壓縮特性、抗剪強度特性等，可通過直剪儀、三軸剪力儀等儀器來測定。盡管直剪儀操作簡便，適合于工程應用，但目前利用直剪儀來測試含水合物沉積物抗剪強度的報道非常少。中國科學院武漢巖土力學研究所是最早開展這方面研究的單位之一，研制了國內(nèi)首臺含水合物土直剪儀（圖1）。該設備將剪切盒安裝于反應釜內(nèi)，通過氣瓶提供生成水合物所需的高壓氣體，用恒溫箱控制反應釜內(nèi)試樣的溫度，最終形成水合物沉積物試樣；通過活塞式加載桿實現(xiàn)法向力和切向力的加載與卸載，用氣瓶輸出壓力提供法向壓力而用恒流泵提供切向推力。Liu等[56-57]采用該設備開展了一系列含CO2水合物沉積物的直剪試驗，揭示了含水合物沉積物的抗剪強度隨水合物飽和度的變化規(guī)律。

三軸剪切試驗是測試含水合物沉積物力學特性最常用的方法，可以測得應力—應變關(guān)系以及變形和強度特性。表1列出了國內(nèi)外各主要研究單位的含水合物沉積物三軸剪切系統(tǒng)的基本信息。各種三軸儀的主要差別在于所采用的試樣制備、水合物含量確定、試樣體變測量等方法上，這些正是含水合物沉積物三軸試驗的難點。在大多數(shù)三軸試驗中，試樣均采用非飽和法制備，通過配置不同初始含水量來控制水合物飽和度。假設試樣中孔隙水完全轉(zhuǎn)化成水合物，Ghiassian等[24]及Hyodo等[66]通過確定試驗結(jié)束后水合物分解所釋放的氣體量反算得到水合物含量大小。目前在三軸試驗中體變測定主要有以下3種方法：①在恒定圍壓時，通過圍壓泵驅(qū)動液體積變化換算獲得[24,65]；②類似于非飽和土試樣體積量測的雙壓力室法[66]；③采用徑向應變傳感器進行量測[67]。為了評價現(xiàn)場所取的含水合物沉積物試樣的力學特性，Yoneda等[68]開發(fā)了一套新型的三軸試驗系統(tǒng)（TACTT），可將原位取得含水合物試樣直接嵌入三軸壓力室進行試驗，并保證整個過程中試樣處于穩(wěn)定狀態(tài)。這套系統(tǒng)采用了圖像處理技術(shù)來監(jiān)測三軸剪切過程中試樣的體積變化。

圖1 含水合物土直剪儀的結(jié)構(gòu)示意圖[56-57]

表1 各種含水合物沉積物三軸試驗系統(tǒng)對比表

目前有關(guān)水合物沉積物的工程力學特性試驗仍然存在許多問題，如測試方法與標準沒有統(tǒng)一、各種試驗數(shù)據(jù)之間缺乏可比性、試樣制備方法單一、試驗結(jié)果不能真實反映原位沉積物的行為特性等。這些均為目前亟待解決的難題，是弄清含水合物沉積物工程力學特性，揭示剛度、強度等物理力學參數(shù)的變化規(guī)律以及識別其關(guān)鍵影響因素所必須解決的關(guān)鍵問題。

2 含水合物沉積物的變形與強度特性

2.1 含水合物沉積物的抗剪強度特性

含水合物沉積物的抗剪強度是指沉積物達到破壞時所能承受的剪應力。Winters等[69]通過試驗證明了水合物沉積物的抗剪強度取決于孔隙中水合物、液體水、天然氣、冰等填充物的相對數(shù)量和空間分布形態(tài)。李洋輝等[70]將水合物顆粒與土顆?；旌现瞥稍嚇硬㈤_展了三軸試驗，發(fā)現(xiàn)水合物含量增加顯著提升了試樣的抗剪強度，且水合物的主要貢獻在于增加黏聚力，而對內(nèi)摩擦角影響較小。Masui等[25]通過開展三軸剪切試驗發(fā)現(xiàn)，冰—砂混合物制得試樣只有在水合物飽和度超過25%時才能提升其抗剪強度，而水—砂混合物制得試樣的抗剪強度隨水合物飽和度的增大而顯著上升。Miyazaki等[71]的試驗結(jié)果表明，含水合物砂的抗剪強度和剛度隨水合物飽和度和有效圍壓的增加而增加，而割線泊松比隨有效圍壓增加而降低。顏榮濤等[23]分別用非飽和法、飽和法制備水合物沉積物試樣并開展三軸試驗，發(fā)現(xiàn)非飽和法制作的水合物試樣抗剪強度隨飽和度增大而顯著提高，而飽和法制作試樣的力學性質(zhì)受水合物含量的影響不大；這種差異性可歸因于水合物賦存模式的不同，非飽和法制備的試樣中水合物以膠結(jié)形式存在，而飽和法產(chǎn)生填充型水合物。Liu等[57]采用直剪儀研究了CO2水合物沉積物的抗剪強度特性，發(fā)現(xiàn)峰值抗剪強度和殘余抗剪強度均隨水合物飽和度和法向應力的增加而提高。綜上所述，水合物含量和賦存模式對含水合物沉積物強度都會產(chǎn)生明顯影響，并且水合物主要貢獻于增加黏聚力，而對內(nèi)摩擦角影響較小[12,23,25,71]。

2.2 含水合物沉積物的剪脹/剪縮特性

研究結(jié)果表明，水合物的形成會影響含水合物沉積物的剪脹/剪縮特性。Miyazaki等[71]用人工合成水合物砂試樣進行三軸試驗，在軸向荷載逐漸增加，試樣逐漸從剪縮變?yōu)榧裘?，偏應力也在到達最高點后迅速下降呈現(xiàn)應變軟化現(xiàn)象。Hyodo等[66]的試驗也證實了水合物沉積物的剪脹性隨著水合物飽和度的升高而明顯增大。Masui等[72]利用水合物沉積物的鉆井巖心進行三軸試驗，發(fā)現(xiàn)與室內(nèi)合成的水合物砂相比，剪脹性弱于人工合成水合物試樣。從微觀層面，Pinkert等[73]認為在剪切作用下水合物顆粒破碎損傷引起了試樣應變軟化的出現(xiàn)，破碎后的水合物顆粒與砂顆粒一起運動，相互翻轉(zhuǎn)跨越導致試驗出現(xiàn)剪脹變形。顯然，水合物生成使含水合物沉積物試樣從開始階段的應變硬化向應變軟化轉(zhuǎn)變，同時伴隨著剪脹現(xiàn)象。隨著水合物飽和度的增加，應變軟化和體積剪脹性會更加明顯[66,74-77]。另外，增大的圍壓會限制顆粒之間的相互翻轉(zhuǎn)跨越運動，使沉積物的剪脹性變?nèi)?，甚至出現(xiàn)剪縮現(xiàn)象[78-79]。

水合物沉積物的剪脹機理可用圖2來進行說明。剪切面兩側(cè)土顆粒受到法向力FN和剪切力FT，作用在每個顆粒上。當顆粒松散排列時切向變形使其排列更密集，從而使土體收縮并顯示出應變硬化行為（圖2a→b）；相反，當顆粒處于密實狀態(tài)時，在剪切過程中則可能發(fā)生體積膨脹即剪脹（圖2b→a）。一旦水合物形成，孔隙空間就會被水合物部分填充，從而使土體的固體骨架變得更加致密，因此剪脹性隨水合物飽和度增大得以提高（圖2-c）。如圖2-d所示，取剪切面上局部土顆粒和水合物顆粒對其進行受力分析，F(xiàn)N和FT的合力與土顆粒和水合物顆粒接觸面構(gòu)成夾角α。在較低圍壓時FN較小，夾角α較小，土顆粒和水合物顆粒發(fā)生滑動跨越導致土體膨脹增大而發(fā)生剪脹，同時由于膠結(jié)結(jié)構(gòu)喪失導致應變軟化變形特征（圖2-e）；當圍壓增大時，F(xiàn)N增大導致夾角α增大，當α大于土顆粒與水合物顆粒之間的摩擦角時，土顆粒與水合物顆粒組成的系統(tǒng)將產(chǎn)生自鎖效應而無法沿接觸面滑動，水合物顆粒被不斷增大的FT壓碎，從而形成更緊密的結(jié)構(gòu)表現(xiàn)為剪縮（圖2-f），在此過程中水合物沉積物應變軟化特征減弱甚至轉(zhuǎn)變?yōu)閼冇不?/p>

2.3 含水合物沉積物的小應變剛度

小應變剛度指的是應變?yōu)?0－5～10－8情況下試樣的變形模量或剪切模量，聲波探測是研究試樣小應變剛度的有效方法。Lee等[81]在砂土、粉土和黏土中形成THF水合物，通過測試試樣的P波和S波研究了小應變剛度特性。其試驗結(jié)果表明，含水合物沉積物的P波和S波受水合物飽和度、有效應力和沉積物比表面積的影響顯著。由于水合物的膠結(jié)和填充作用，含水合物沉積物的小應變體積模量和剪切模量會隨著水合物增大呈現(xiàn)明顯的上升趨勢[82-84]。Priest等[84]采用“富氣法”和“富水法”分別制備含水合物砂土試樣，發(fā)現(xiàn)兩種試樣顯示明顯不同的波速響應。對于“富氣法”形成試樣，水合物主要起到膠結(jié)作用，少量的水合物就可以明顯提高含水合物沉積物的剛度；而“富水法”形成試樣，水合物主要填充于孔隙中，只有當水合物達到一定含量時，水合物賦存模式從填充型轉(zhuǎn)變?yōu)槌至w型，含水合物沉積物的剛度才能得以提升。這一結(jié)論在Choi等[34,82]所進行的實驗室中也得到了驗證。可見，水合物含量及賦存模式對沉積物的小應變（10－5～10－8）剛度具有顯著的影響[9-11,85]。

圖2 水合物沉積物剪切變形機理示意圖

基于水合物含量及賦存模式對含水合物沉積物試樣的小應變剛度的影響機制，研究人員針對于不同水合物賦存模式建立了預測含水合物沉積物波速的理論模型，通過測定剪切波速、壓縮波速與水合物飽和度關(guān)系來綜合判別沉積物中水合物的賦存模式[86-87]。這種水合物賦存模式的識別方法目前已得到一定程度的應用[34,88-90]。

綜合上述分析可知，水合物含量及其賦存模式對沉積物強度和變形特性的影響已經(jīng)成為當前重點研究的問題。盡管現(xiàn)有的研究成果已初步揭示了水合物含量和賦存模式對沉積物物理力學特性的影響規(guī)律，但囿于制樣方法與實驗技術(shù)，目前對含水合物沉積物物理力學特性的研究仍然存在以下幾個問題：①對在水合物含量變化下填充型向持力體型轉(zhuǎn)變過程中沉積物的物理力學變化規(guī)律缺乏深入系統(tǒng)的認識；②對水合物分解過程中因水合物局部分解而導致的含水合物沉積物抗剪強度的衰減及壓縮性的增加規(guī)律幾乎沒有涉及；③流泥流沙情況下含水合物沉積物骨架的物理力學特性衰化規(guī)律仍然沒有厘清。另外，目前有關(guān)沉積物物理力學特性的試驗數(shù)據(jù)均在特定試驗條件下獲得，缺乏一致性、可比性和系統(tǒng)性，需要針對典型水合物儲層介質(zhì)的物理力學特性開展系統(tǒng)深入的試驗研究工作。

3 含水合物沉積物的物理力學模型

3.1 沉積物體系中水合物相平衡模型

沉積物中水合物相平衡模型描述了水合物在平衡時溫度、壓力與其他狀態(tài)變量（如水合物含量、液態(tài)水含量、含鹽量等）之間的關(guān)系，是含水合物沉積物最基礎(chǔ)的本構(gòu)關(guān)系之一。最早的水合物相平衡模型是由 van der Waals 等[91]所建立的，主要預測非沉積物中純水合物的相平衡壓力和溫度條件。在沉積物中，由于沉積物顆粒表面物理吸附和孔隙毛細管效應的影響，孔隙水活度降低影響了相平衡條件，因此相平衡模型應該考慮沉積物孔隙大小及分布特征的影響。Clennell等[92-93]將Gibbs-Thomson方程引入vdWP模型，從而考慮了孔徑引起的毛細管效應，建立了沉積物中水合物相平衡模型。該模型只適合于描述具有單一或窄孔徑分布的沉積物中水合物相平衡條件。為突破這一局限性，顏榮濤等[94]引入了孔徑分布函數(shù)，并根據(jù)水合物分解的孔徑大小順序，確定了孔徑大小分布與水合物飽和度的關(guān)系，從而得到了p—T—Sh三維相平衡關(guān)系。最近，利用沉積物孔徑分布與持水特性的內(nèi)在關(guān)聯(lián)性，顏榮濤等[95]通過引入持水曲線，建立了沉積物中水合物的宏觀相平衡模型。

近年來，Wei[96]基于多孔材料的連續(xù)介質(zhì)理論建立了能考慮骨架與孔隙溶液之間復雜物理化學作用的多相巖土介質(zhì)化學—力學耦合理論，并提出了能考慮滲透、毛細管、吸附等效應的孔隙溶液各組分化學勢的一般數(shù)學表達式。在此基礎(chǔ)上，Zhou等[97]建立了能夠考慮復雜物理化學效應的水合物相平衡方程，該方程給出了平衡溫度偏移量與毛細管吸力和鹽溶液濃度的關(guān)系，揭示了沉積物持水特性與相平衡條件的內(nèi)在聯(lián)系。結(jié)合所提出的相平衡方程和沉積物的持水特征曲線，進而提出了水合特征曲線（SHCC）的概念，即含水合物沉積物中液態(tài)水含量與平衡溫度偏移量的關(guān)系曲線。在p—T—w三維空間中（w為液態(tài)水含量），SHCC描述了沉積物中水合物的相邊界。如圖3所示，對于任意指定的水合物沉積物，在平衡時含水合物沉積物中的溫度、壓力和液態(tài)水之間存在唯一關(guān)系，這與實驗觀測結(jié)果是完全一致的。該模型克服了傳統(tǒng)水合物相平衡模型不能預測任意平衡的溫度和壓力下沉積物中液態(tài)水含量問題，不僅形式簡單（沒有涉及微觀參數(shù)），而且還能綜合考慮鹽分、毛細管以及物理化學吸附的影響。

圖3 在p—T—w空間中含CO2水合物沉積物的相平衡面及其在各平面上的投影[97]

3.2 含水合物沉積物水力模型

含水合物沉積物的水力模型包括滲透率函數(shù)模型和持水曲線模型，其中，滲透率函數(shù)不是嚴格意義上的本構(gòu)關(guān)系，而只是一條輔助的經(jīng)驗關(guān)系式。在處理土體滲透性時，通常的做法是根據(jù)孔隙分布特征把孔隙簡化成不同孔徑大小的平行的毛細管，再利用流體力學的原理來確定其滲透性[98]?；谶@種思想，考慮到孔隙包裹型和孔隙填充型對平行毛細管特征的影響，Kleinberg[50]分別建立了能夠考慮孔隙包裹型與填充型水合物賦存模式的水合物沉積物滲透率模型。然而，這類模型沒有考慮土體中滲流路徑的曲折性，也沒有考慮水合物形成后對孔隙聯(lián)通性的改變，因此該模型的預測結(jié)果與實際情況有一些差距。在自然界中，根據(jù)所處環(huán)境條件與水合物含量，包裹型和填充型水合物可同時出現(xiàn)于海底沉積物中，并且同時對滲透率存在影響，只是程度不同而已。針對于這種情況，Delli等[99]提出了一個混合式模型，該模型可以考慮包裹型和填充型兩種水合物模式共存的情況。此外，對于水合物分布不均勻的情況，Dai等[100]修正了Kozeny-Carman模型，較好地模擬了滲透率隨水合物飽和度的變化。

含水合物沉積物的持水特性控制了氣—水滲流過程，并影響著沉積物的其他工程力學特性。在有關(guān)水合物開采的邊界值問題分析中[101-102]，一般均采用非飽和土的van Genuchten模型來描述含水合物沉積物的持水特性[103]。然而，由于水合物的生成/分解，含水合物沉積物的孔隙大小及分布都在不斷變化中，因此沉積物持水特性也隨之發(fā)生變化。Gupta等[104]通過引入修正系數(shù)來考慮水合物飽和度和孔隙率對毛細管吸力的影響。該模型主要是從模仿孔隙中微生物生長對骨架產(chǎn)生的影響修正而來，并沒有真正考慮水合物生成對土體結(jié)構(gòu)的影響。Yan等[105]分析了水合物形成對含水合物土毛細管吸力的影響機理，考慮了水合物填充的效應，建立了含水合物沉積物的土水特征曲線模型。雖然該模型很好地描述了含水合物沉積物的持水特征關(guān)系，特別是水合物對進氣值和殘余飽和度的影響，但是其有效性仍需要進一步驗證。

總體上，現(xiàn)有的多數(shù)含水合物沉積物的水力模型（滲透率函數(shù)和持水曲線模型）都是基于高度簡化的含水合物沉積物結(jié)構(gòu)，對現(xiàn)有的水力模型修正來建立的，所以其模擬精度及適用性受到限制。因此，后期發(fā)展的水力模型應該能夠考慮水合物含量及賦存模式對孔隙水/氣遷移過程的影響，在大量試驗結(jié)果基礎(chǔ)上，建立能統(tǒng)一考慮水合物含量及賦存模式影響的含水合物沉積物的滲透率函數(shù)和持水曲線模型。

3.3 應力應變本構(gòu)關(guān)系模型

含水合物沉積物的應力應變本構(gòu)關(guān)系大體上可分為非線性彈性模型和彈塑性模型兩大類。前者主要包括Duncan-Chang類[106-109]和彈性損傷模型[110-113]等。Miyazaki等[106]建立了含水合物沉積物強度和剛度與飽和度的經(jīng)驗關(guān)系式，修正Duncan-Chang模型來模擬了含水合物沉積物的應力應變關(guān)系。Yu等[107]把含水合物沉積物在受荷過程中的破壞分成結(jié)構(gòu)性快速破壞和結(jié)構(gòu)性完全破壞（即屈服）兩個階段，并且在初始切線模量和極限偏應力中考慮了溫度的影響，建立了包含溫度、圍壓、應變速率等不同參數(shù)的修正Duncan-Chang本構(gòu)模型。這類模型的優(yōu)勢在于形式簡單，參數(shù)易于測定，易于在數(shù)值模型中進行求解。然而，該類模型雖然考慮了水合物飽和度對試樣強度和剛度的增強效應，但不能模擬應變軟化及剪脹行為。下一步工作中，可以重點考慮對Duncan-Chang本構(gòu)模型進行修正，使其能反映含水合物沉積物的應變軟化或剪脹行為。

吳二林等[110-111]較早引入彈性損傷概念來模擬含水合物沉積物的本構(gòu)行為，后期經(jīng)過了李彥龍、祝效華、顏榮濤等[112-114]研究者的發(fā)展。該模型的主要思想是將沉積物視為由水合物和土體骨架組成的復合材料，利用混合律計算出沉積物的等效彈性模量，將統(tǒng)計損傷、概率等概念考慮進去，使其能反映沉積物受力時水合物對微結(jié)構(gòu)和應力應變的影響，從而建立含水合物沉積物的彈性損傷模型。此類模型較好地反映了剪切過程中含水合物沉積物中的損傷過程，從而有效模擬了含水合物沉積物的應變軟化行為，但該模型不能很難描述剪切過程中試樣體積的變化。

Pinkert及其合作者主要發(fā)展了含水合物沉積物的Mohr-Coulomb類本構(gòu)模型[73,115-116]，該模型很好地模擬了含水合物沉積物的應變軟化效應，采用的Rowe應力剪脹理論對含水合物砂土剪脹行為有較好的描述。盡管如此，Mohr-Coulomb類模型還是沒有充分考慮體積屈服效應以及剪切過程中水合物膠結(jié)作用衰減或喪失所帶來的影響。

為了更全面地描述含水合物沉積物的本構(gòu)行為，一些學者嘗試在臨界狀態(tài)彈塑性理論框架內(nèi)建立水合物沉積物的本構(gòu)模型。目前比較常用的方式主要是通過改變屈服面大小來反映水合物的影響，結(jié)合硬化規(guī)律、屈服函數(shù)、流動法則和加卸載準則及彈性行為模擬，給出了一套內(nèi)容完整、物理意義明確、結(jié)構(gòu)清晰的理論體系。例如，Klar等[102]所建立的理想彈塑性模型、Sultan等[117]和Uchida等[118]提出的修正劍橋模型、Zhang等[119]提出的修正亞塑性模型等、Lin等[120]基于SMP準則的彈塑性本構(gòu)模型等。在目前發(fā)表的臨界狀態(tài)彈塑性模型中，主要關(guān)注含水合物沉積物的以下幾個典型特征：①強度和剛度隨水合物含量的增大；②應變軟化效應會隨著水合物含量的增加而增強；③隨著水合物含量的增大，試樣的剪脹效應會明顯增大[117-118,120-123]。這些模型主要通過修正原劍橋模型中的屈服函數(shù)和剪脹函數(shù)來考慮水合物含量對含水合物沉積物力學特性的影響規(guī)律。這些屈服函數(shù)和剪脹函數(shù)的建立主要依賴于理論假設，缺少足夠的試驗數(shù)據(jù)支持。此外，一些其他影響因素也在臨界狀態(tài)模型中進行了考慮。例如，Shen等[124]分析了水合物對臨界狀態(tài)以及剪脹函數(shù)的作用，建立了依賴密實度的含水合物砂土的臨界狀態(tài)模型；Yan等[121]重新定義了含水合物土的有效應力，并且在屈服函數(shù)、流動法則和硬化規(guī)律中嵌入了水合物膠結(jié)和毛細管吸力的影響，提出了有效考慮水合物飽和度和賦存模式對力學特性、模擬分解誘發(fā)的土體變形的含水合物土的本構(gòu)模型；鄒遠晶等[125]考慮重點分析了受荷情況下含水合物沉積物膠結(jié)結(jié)構(gòu)的破壞情況，基于擾動狀態(tài)概念建立了描述膠結(jié)結(jié)構(gòu)破壞的全過程的含水合物土的應力—應變模型。

近年來，含水合物沉積物的應力應變本構(gòu)關(guān)系模型發(fā)展迅速，新的模型變得越來越完備，除了能夠考慮一般土體的主要力學行為特征外，還能模擬水合物含量及賦存模式的影響，但代價是模型變得越來越復雜，包含了許多物理意義不明確的且難以獲取的模型參數(shù)，從而導致了這些新模型很難應用于解決實際工程問題。

4 水合物開采土力學多場耦合過程的數(shù)值模擬方法

數(shù)值模擬方法為水合物儲層的開采潛力評價及開采場地安全性預測提供了一條重要途徑[126-127]。然而，針對水合物開采的數(shù)值模擬需要求解復雜的非線性多過程耦合問題。由于存在強耦合作用，常常會出現(xiàn)物理上的奇異和數(shù)值計算上的不穩(wěn)定。另外，水合物儲層是一種多相孔隙介質(zhì)，其物理和化學性質(zhì)較強地依賴于水合物含量、溫度、壓力等狀態(tài)變量，強烈的非線性特征使得求解的困難程度大為增加。因此，發(fā)展針對水合物開采模擬的穩(wěn)定高效數(shù)值計算方法目前依然是重要的科學難題。

針對于開采潛力評價的數(shù)值分析（針對于產(chǎn)能預測），目前國際上已研發(fā)出一系列較為成熟的水合物開采數(shù)值模擬器，如美國勞倫斯伯克利國家實驗室的TOUGH+HYDRATE[128]、日本石油工程公司和東京大學合作開發(fā)的MH-21HYDRES[129]和加拿大CMG公司的CMG-STARS[130]等。這些模擬器主要是基于孔隙介質(zhì)的質(zhì)量和能量守恒方程，考慮了相變、傳熱、液氣運移和組分擴散過程，忽略了土體變形/破壞對其他物理場的影響。筆者主要針對開采過程中含水合物儲層的穩(wěn)定性和井壁穩(wěn)定性的研究進行論述。寧伏龍[3]研究了水合物地層鉆探過程中出現(xiàn)的井壁穩(wěn)定性問題，指出實際鉆井過程中井壁圍巖地層的孔壓和應力分布情況是流體流動和固體變形共同作用的結(jié)果。程遠方等[131]基于有限元法，采用流固耦合模型對降壓開采過程中儲層的應力狀態(tài)以及物性參數(shù)的動態(tài)變化規(guī)律進行了研究。然而，該模型忽略了溫度場的變化，并采用了解耦方法對滲流場和變形場進行了單向順次求解。Rutqvist等[132-133]采用單向耦合的方法將TOUGH+HYDRATE和FLAC結(jié)合起來，采用前者計算溫度、孔壓、水合物含量等，模擬水合物儲層的傳熱、滲流以及水合物分解過程，采用后者計算應力應變，模擬分析儲層的變形和穩(wěn)定性問題。兩者之間通過專門設計的通信模塊來實現(xiàn)數(shù)據(jù)的單向傳遞。與此思路相同，Qiu等[134]通過單向耦合MH21-HYDRES和3D FE地質(zhì)力學模擬器分析評價了日本Nankai海槽水合物藏降壓開采過程中的井筒完整性。

為了提高計算精度，Kim等[135]進一步發(fā)展了TOUGH+HYDRATE和FLAC的交錯耦合模型，模擬水合物分解過程中的流體滲流和巖土體變形之間的耦合作用。與單向耦合法相比，交錯耦合法（即雙向耦合法）將信息由單向傳遞變?yōu)殡p向傳遞，使計算過程能夠符合不同物理場之間是相互影響的這一基本事實。為了克服FLAC作為商業(yè)軟件在源代碼和內(nèi)存共享等方面的局限，Queiruga等[136]開發(fā)了有限元程序Millstone以替換FLAC來計算儲層的力學行為，并提出了新的不同網(wǎng)格間的數(shù)據(jù)傳遞方法，有效提高了計算效率。需要指出的是，單向耦合法和交錯耦合法同屬于解耦算法，雖然在計算方法上相對容易實現(xiàn)，可以利用已有的較為成熟的熱傳導、流體滲流和固體變形等單場求解程序，但未能將各個物理場放在同一地位進行考慮，理論上不能稱為嚴格的耦合算法。更重要的是，在處理像水合物開采這種瞬態(tài)非線性問題時，要統(tǒng)一協(xié)調(diào)不同軟件之間的時域積分和非線性迭代，解耦合法或多或少存在一些“相延遲”問題，通常不得不采用較小的時間步長，很難把握計算過程的收斂性，甚至會出現(xiàn)無法模擬某些重要的物理現(xiàn)象的情況，例如Biot流固耦合問題中經(jīng)典的Mandel-Cryer效應。為此，不少學者開始嘗試將全耦合方法引入到水合物開采的多場耦合模擬中。

Salehabadi等[137]利用力學商業(yè)軟件ABAQUS結(jié)合Heriot-Watt大學開發(fā)的水合物分解模擬軟件HWHYD，研究了水合物分解條件下套管周圍應力場的變化，分析了水合物分解對井眼及套管的穩(wěn)定性的影響。該研究中對流固耦合問題采用了全耦合方法，但對熱流耦合和熱固耦合問題只進行了單向耦合求解，即只考慮傳熱對滲流和變形的作用，而忽略了滲流和變形對溫度場的影響。Klar等[138]推導了氣液兩相滲流方程和動量平衡方程的顯式有限差分格式，采用了理想彈塑性模型來描述固體骨架的變形，最終通過對FLAC軟件的二次開發(fā)實現(xiàn)了解耦求解。不久前，Klar等[138]繼續(xù)補充了熱傳導方程，進一步考慮了溫度場的變化對開采的影響。Uchida等[139]首次將顆粒運移問題耦合到水合物開采過程中，使得該程序可模擬儲層出砂過程。模型中考慮了沉積物變形、顆粒分離、顆粒運移、顆粒沉降等過程，通過數(shù)值分析評估了出砂相關(guān)參數(shù)對固體骨架體變和水合物含量的敏感性。最近，Zhou等[140]在該程序中嵌入了修正的臨界狀態(tài)模型來模擬水合物儲層的力學行為，并對日本Nankai海槽東部海域的水合物開采過程進行了分析。由于采用了顯示算法，該數(shù)值模型雖然考慮了眾多參數(shù)的非線性特征，但并未引入迭代過程，而是直接采用上一時間步的結(jié)果對當前時間步的參數(shù)進行更新，其計算結(jié)果的準確性及其對時間步長的依賴性仍值得進一步研究。

除此之外，Kimoto等[101,141]、Raz[142]、吳二林等[143]也自主發(fā)展了水合物儲層開采模擬的多場耦合數(shù)學模型，并基于實測數(shù)據(jù)開展了地層力學參數(shù)演化及沉降變形預測的研究。與之前類似，上述研究并不完善，仍有部分計算過程采用了解耦策略，忽略或淡化了某些耦合作用，這與實際發(fā)生的強耦合物理過程是不相符合的。此外，當采用解耦合方法時，往往是直接調(diào)用現(xiàn)有的單場軟件。不同商業(yè)軟件對輸入和輸出數(shù)據(jù)的開放程度有很大的不確定性。受到可選用的單場軟件的限制，不同物理場的耦合求解經(jīng)常不得不采用不同的數(shù)值方法，如滲流場采用有限體積法、位移場采用有限元法、溫度場采用有限差分法。如何在不同類型的計算軟件之間協(xié)調(diào)網(wǎng)格剖分、統(tǒng)一邊界條件、設計數(shù)據(jù)傳遞渠道，對于解耦方法來說目前仍然困難重重。

5 發(fā)展趨勢與未來的挑戰(zhàn)

5.1 物理力學特性測試方法和技術(shù)

目前發(fā)展出的多功能探測技術(shù)主要以定性為主，圍繞沉積物組構(gòu)的精細定量探測手段相對缺乏。CT、聲波和核磁共振技術(shù)可以對土組構(gòu)進行無損、定量的探測，但由于在信號解析上缺乏成熟的理論模型，這些技術(shù)在應用中仍存在困難。隨著CT、聲波技術(shù)與核磁共振技術(shù)的不斷發(fā)展，沉積物組構(gòu)探測會逐漸走向定量化，并將與其他探測技術(shù)以及沉積物工程力學特性測試相結(jié)合，這是一種發(fā)展趨勢。在沉積物工程力學特性測試方面，試樣制備—組構(gòu)探測—力學測試一體化、多種測試功能集成化、面向原位取樣的測試方法、模擬原位開采工況的大型模型試驗等將成為未來的重要研究方向。另外，由于含水合物沉積物的結(jié)構(gòu)易損性，沉積物組構(gòu)探測與工程特性的原位和船載測試技術(shù)會得到進一步發(fā)展，除了傳統(tǒng)的海洋原位土工試驗、地球物理、聲吶等技術(shù)會進一步完善外，新的測試技術(shù)（如水下光纖、機器人等）將越來越多地被應用。同時，室內(nèi)與原位測試會更加緊密結(jié)合，相互取長補短。

5.2 含水合物沉積物的力學特性研究

在含水合物沉積物的力學特性方面，目前的試驗數(shù)據(jù)比較零散且各類數(shù)據(jù)相互獨立，缺乏針對典型水合物儲層介質(zhì)的系統(tǒng)性試驗數(shù)據(jù)。這些問題對全面、準確地認識含水合物沉積物的力學特性及其變化規(guī)律是非常不利的。解決這一問題需要針對一些典型水合物儲層，系統(tǒng)開展相關(guān)沉積物的水力特性、力學特性以及相變特性的研究，建立完備的沉積物儲層物理力學特性數(shù)據(jù)庫。另外，隨著水合物分解/生成，水合物含量及其賦存模式在不斷變化中，使得水合物儲層介質(zhì)具有明顯的動態(tài)結(jié)構(gòu)性，因此，厘清水合物儲層介質(zhì)的組構(gòu)特性及其關(guān)鍵影響因素、揭示組構(gòu)演化對沉積物力學特性的影響機理和規(guī)律也將成為未來有關(guān)水合物研究的重點方向。

5.3 含水合物沉積物的物理力學模型

開采過程中水合物分解會導致儲層介質(zhì)會發(fā)生變形沉降，分析水合物儲層沉降對開采系統(tǒng)（特別是海底管線）所造成的影響，需要進一步發(fā)展能夠描述水合物開采過程的儲層介質(zhì)多場、多過程耦合理論，建立能有效模擬在壓力、溫度、化學等復雜荷載作用下含水合物沉積物力學行為的實用本構(gòu)模型。由于含水合物沉積物具有顯著的相變性和動態(tài)結(jié)構(gòu)性（其結(jié)構(gòu)隨著水合物分解、流泥流沙過程而不斷發(fā)生變化），需要對水合物沉積物組構(gòu)進行定量表征，并在構(gòu)建強度準則和本構(gòu)模型時充分考慮組構(gòu)演化的影響。盡管開采場地的沉降問題需要考慮，但相對而言，儲層破壞（如海底滑坡）、井壁失穩(wěn)等破壞問題更加突出。針對這些問題，需要發(fā)展相應的計算模型和分析方法，包括：井壁穩(wěn)定性分析、海底穩(wěn)定性分析、錨固基礎(chǔ)抗拔承載力計算等，因此，含水合物沉積物的強度特性及其影響因素和變化規(guī)律將是未來需要重點關(guān)注的問題。

5.4 水合物開采土力學多場耦合數(shù)值模擬

水合物開采中產(chǎn)能預測以及開采場地力學安全性分析涉及到復雜的非線性問題求解。隨著水合物儲層特征刻畫的精細化、沉積物多過程耦合理論的進一步發(fā)展以及數(shù)值方法和計算技術(shù)的不斷進步，這一方向的發(fā)展趨勢具有以下兩種特征：①從分析尺度上看，逐漸實現(xiàn)從儲層介質(zhì)的微細結(jié)構(gòu)演變到宏觀場地尺度的變形破壞的跨尺度模擬分析；②從數(shù)值求解方法看，正在從傳統(tǒng)的單向耦合、交錯耦合等解耦算法轉(zhuǎn)到全耦合算法上，逐漸實現(xiàn)水合物開采邊值問題數(shù)值解法的根本性轉(zhuǎn)變。

5.5 含水合物沉積物力學研究面臨的挑戰(zhàn)

經(jīng)過近20多年來的發(fā)展，有關(guān)水合物開采的理論、方法和技術(shù)均取得了顯著的進展，但從巖土力學的角度看，目前仍然面臨諸多挑戰(zhàn)。

5.5.1 大尺度含水合物沉積物試樣的人工制備

開展室內(nèi)試驗和物理模型試驗需要在實驗室中制備大尺度均勻的含水合物沉積物試樣。由于含水合物沉積物物理力學特性強烈依賴于沉積物的結(jié)構(gòu)，因此所制備試樣必須能夠真實反映原位含水合物沉積物的結(jié)構(gòu)特征，特別是孔隙水合物的賦存模式應與原位一致。迄今為止，這仍然是沒有完全解決的難題。

5.5.2 含水合物沉積物微細觀組構(gòu)的精細探測與定量表征

含水合物沉積物具有顯著的結(jié)構(gòu)性，而且水合物分解/生成使得沉積物的微細觀結(jié)構(gòu)處于動態(tài)演化中。因此，厘清含水合物沉積物的物理力學特性的變化規(guī)律必須能夠?qū)ζ湮⒓氂^結(jié)構(gòu)進行精細探測與定量表征，這正是目前有關(guān)含水合物沉積物力學研究中面臨的挑戰(zhàn)之一。

5.5.3 組構(gòu)變化對含水合物沉積物力學特性的影響機理與規(guī)律

水合物分解會導致其膠結(jié)作用破壞，從而引起沉積物的強度和變形特性的變化。因此，闡明水合物含量、水合物賦存模式等組構(gòu)因素變化對沉積物宏觀力學特性的影響機理與規(guī)律，是有效評價水合物儲層力學穩(wěn)定性的重要前提。由于缺乏相應的試驗數(shù)據(jù)，這一關(guān)鍵問題目前尚未解決。

5.5.4 開采擾動下流沙發(fā)生條件及其對儲層穩(wěn)定性的影響機制

實現(xiàn)水合物高效開采需要對開采擾動下流（出）沙過程進行有效控制，因此需要闡明流沙過程的發(fā)生機理和觸發(fā)條件。另外，流沙過程也會導致儲層或開采井壁破壞，從而影響到開采過程。目前對水合物開采過程的流沙問題認識有限，研究手段單一（主要基于離散元法），沒有建立有效的流沙判別準則，缺乏考慮流沙影響的儲層穩(wěn)定性分析方法。

5.5.5 開采擾動下多相多組分含水合物儲層的多過程耦合問題

含水合物儲層是一種多相多組分孔隙介質(zhì)，在開采擾動下，會發(fā)生熱傳導、氣—液流動、物質(zhì)擴散、細粒流失（流泥流沙）、骨架變形、水合物生成/分解等耦合的多物理化學過程。雖然傳統(tǒng)理論很好模擬熱傳導—滲流/擴散—骨架變形等耦合過程，但對存在相變等化學過程情況仍存在明顯的局限性。另外，現(xiàn)有的多孔介質(zhì)連續(xù)介質(zhì)力學模型沒有充分考慮流沙問題的影響。

5.5.6 水合物分解/生成條件下含水合物沉積物的本構(gòu)響應

在水合物分解或生成條件下，含水合物沉積物的內(nèi)部組構(gòu)可以在膠結(jié)模式—填充模式—持力體模式之間相互轉(zhuǎn)化。目前囿于試驗技術(shù)的限制，對這類組構(gòu)不斷變化的土體的本構(gòu)響應特征缺乏了解。已有的本構(gòu)模型無法準確模擬組構(gòu)變化的影響，而近年來發(fā)展起來的能考慮賦存模式影響的沉積物本構(gòu)關(guān)系模型不僅結(jié)構(gòu)復雜，而且涉及到模式判別問題，這限制了新模型的應用。因此，發(fā)展能有效模擬水合物分解/生成條件下含水合物沉積物本構(gòu)關(guān)系模型仍然是目前亟需解決的難題。

5.5.7 有關(guān)水合物開采過程數(shù)學模型的適定性問題

目前模擬水合物開采過程的數(shù)學模型都是通過對原本只適用于描述單一過程的數(shù)學模型的簡單組裝而得到，如熱傳導采用的是經(jīng)典的傅里葉法則，滲流過程采用的是達西定律，擴散過程則是菲克斯第一法則等。為了有效模擬水合物開采過程，必須首先回答該類數(shù)學模型解的存在性、穩(wěn)定性和唯一性問題。

5.5.8 高效穩(wěn)定的水合物開采過程的數(shù)值模型的全耦合解法

有關(guān)水合物開采過程模擬所采用的解耦或半解耦解法與實際物理過程嚴重不符，不僅精度低而且難以控制。因此，發(fā)展全耦合計算方法將是一種理想選擇。但是，采用全耦合計算方法除了需要從頭開始進行軟件設計開發(fā)，還需要發(fā)展高效穩(wěn)定的全耦合數(shù)值算法，這是一項亟待解決的重大科學問題。