天然氣水合物沉積物的強(qiáng)度模型

顏榮濤，趙續(xù)月，楊德歡，肖桂元，梁維云，文松松

（桂林理工大學(xué)　　廣西巖土力學(xué)與工程重點(diǎn)實(shí)驗(yàn)室，廣西　　桂林 541004）

天然氣水合物沉積物的強(qiáng)度模型

顏榮濤，趙續(xù)月，楊德歡，肖桂元，梁維云，文松松

（桂林理工大學(xué)廣西巖土力學(xué)與工程重點(diǎn)實(shí)驗(yàn)室，廣西桂林 541004）

建立有效描述含天然氣水合物（以下簡(jiǎn)稱“水合物”）沉積物強(qiáng)度特性的強(qiáng)度模型是水合物開發(fā)前必須解決的關(guān)鍵問題之一。在調(diào)研前人對(duì)含水合物沉積物強(qiáng)度試驗(yàn)結(jié)果的基礎(chǔ)上，分析了水合物存在對(duì)土體強(qiáng)度的兩種影響機(jī)理，分別是水合物作為持力體分擔(dān)土骨架部分應(yīng)力以及土顆粒之間的膠結(jié)作用。在此基礎(chǔ)上，把含水合物沉積物視為復(fù)合材料（土顆粒骨架相和水合物相復(fù)合）進(jìn)行考慮，建立了含水合物沉積物的強(qiáng)度模型，該模型能較好地反映不同水合物賦存模式、不同水合物含量、不同圍壓作用下的含水合物沉積物的強(qiáng)度特性。同時(shí)，該模型還能很好地反映細(xì)粒土中水合物形成后內(nèi)摩擦角的減小。然而，由于前期相關(guān)試驗(yàn)數(shù)據(jù)的缺乏，該模型的一些關(guān)系式還有所欠缺，在后期的工作中需要進(jìn)一步的優(yōu)化。

含水合物沉積物；賦存模式；強(qiáng)度模型

0 引 言

天然氣水合物是由天然氣（主要以甲烷為主）和水在高壓、低溫環(huán)境下形成的一種類似冰狀的化合物（俗稱“可燃冰”），主要存在于永久凍土區(qū)以及海底陸坡等區(qū)域。天然氣水合物儲(chǔ)量豐富，其分解后將釋放出大量的甲烷氣體，學(xué)術(shù)界認(rèn)為它是一種不錯(cuò)的未來替代性能源。然而，天然氣水合物開發(fā)過程可能會(huì)遇到一系列的地質(zhì)安全和環(huán)境污染等問題。例如，天然氣水合物的開采弱化了地層的穩(wěn)定性，有可能導(dǎo)致含天然氣水合物地層大規(guī)模沉降，甚至海底滑坡；再則，如果開采的甲烷氣體來不及捕獲而釋放到大氣中，將會(huì)進(jìn)一步加劇全球變暖等環(huán)境問題［1-5］。

土體強(qiáng)度是控制地層穩(wěn)定性的關(guān)鍵因素之一，深刻認(rèn)識(shí)含天然氣水合物沉積物的強(qiáng)度特性并且建立含天然氣水合物沉積物的強(qiáng)度模型是對(duì)含天然氣水合物地層進(jìn)行穩(wěn)定性評(píng)估的前提條件之一。目前，對(duì)含水合物沉積物的力學(xué)特性研究主要集中在強(qiáng)度的變化規(guī)律及其強(qiáng)度貢獻(xiàn)機(jī)理［6-14］，很少有針對(duì)含天然氣水合物沉積物建立強(qiáng)度模型的研究。Zhang等［15］對(duì)含甲烷水合物和冰粉土進(jìn)行三軸試驗(yàn)后建立了一個(gè)經(jīng)驗(yàn)性的強(qiáng)度模型，該模型能較好地反映水合物對(duì)土體加固的作用，然而該模型缺乏理論基礎(chǔ)，很難推廣到其他工況，并且所需參數(shù)較多；Li等［16］在試驗(yàn)的基礎(chǔ)上建立高圍壓條件下的強(qiáng)度模型，該模型能很好地模擬在高圍壓階段強(qiáng)度隨圍壓降低的特性，但同樣缺乏理論基礎(chǔ)。很明顯，含水合物沉積物強(qiáng)度模型的研究現(xiàn)狀極大地限制了含天然氣水合物地層穩(wěn)定性的分析研究，阻礙水合物商業(yè)開采的進(jìn)程。

針對(duì)這個(gè)情況，筆者試圖在前人對(duì)含天然氣水合物沉積物的力學(xué)試驗(yàn)基礎(chǔ)上，分析總結(jié)水合物在沉積物中形成對(duì)其強(qiáng)度的影響機(jī)理，從理論上推導(dǎo)出含水合物沉積物的強(qiáng)度模型，并基于試驗(yàn)結(jié)果建立針對(duì)不同賦存條件、不同水合物含量條件下的強(qiáng)度模型。

1　水合物對(duì)沉積物的強(qiáng)度影響

水合物于土中的形成會(huì)影響沉積物的力學(xué)和水力性質(zhì)［15，17］。一般地，水合物在沉積物中形成，沉積物的強(qiáng)度會(huì)提升，但是提升的幅度以及影響規(guī)律取決于水合物賦存模式、水合物含量等。

根據(jù)大多數(shù)學(xué)者的調(diào)查和研究，自然界中水合物存在3種賦存模式［18］：第1種賦存模式為孔隙懸浮模式，這種模式下水合物顆粒主要懸浮于孔隙水中；第2種賦存模式為持力體模式，水合物顆粒搭連于土顆粒之間成為持力體；第3種賦存模式為膠結(jié)模式，水合物膠結(jié)于土顆粒之間，使兩個(gè)原本相對(duì)獨(dú)立的土顆粒成為一個(gè)整體（圖1）。前兩種賦存模式，可以統(tǒng)一稱之為“填充”模式。在這3種賦存模式之中，第1種填充模式的土顆?？紫吨兴衔锾幱趹腋顟B(tài)，對(duì)于土體強(qiáng)度沒有明顯的提升，但對(duì)于滲透性是存在影響的；而隨著水合物的成長(zhǎng)，當(dāng)水合物的大小可以搭連土顆粒時(shí)，水合物顆粒可以承擔(dān)力之間的傳遞，從而增加了強(qiáng)度，這即為第2種填充模式；當(dāng)然，如果地層地藏條件許可，水合物還可以使土體裂開繼續(xù)生長(zhǎng)，成為透鏡體形式。這種水合物模式下，水合物對(duì)強(qiáng)度的提升只有當(dāng)水合物的含量達(dá)到了一定的量才能體現(xiàn)出來，這個(gè)定值稱之為臨界飽和度。Hyodo等［12-13］根據(jù)試驗(yàn)結(jié)果認(rèn)為臨界飽和度值取25%，而Waite等［17］總結(jié)前人的成果，發(fā)現(xiàn)臨界飽和度處于25%～40%。

第3種賦存模式和前兩種在形成模式上存在本質(zhì)的區(qū)別，這種模式下含量很小的水合物就會(huì)對(duì)土體的強(qiáng)度有很明顯的提升［11，17］，這主要是水合物膠結(jié)于土顆粒之間，產(chǎn)生了膠結(jié)強(qiáng)度。當(dāng)然，在水合物含量增大時(shí)，水合物除了膠結(jié)土顆粒，也會(huì)承擔(dān)相應(yīng)的荷載，與土顆粒共同承擔(dān)剪切應(yīng)力。

圖1 土體中水合物的賦存模式示意圖［17］Fig.1 Occurrence sketches of gas hydrate in sediments

從以上的分析來看，水合物對(duì)于土體的強(qiáng)度增強(qiáng)機(jī)理大致可以總結(jié)為兩類：一是膠結(jié)于土顆粒之間；二是成為傳遞力的載體，充當(dāng)類似于土顆粒的角色，與土顆粒共同承擔(dān)剪切應(yīng)力。為了方便后面敘述，將其分別稱為膠結(jié)貢獻(xiàn)和填充貢獻(xiàn)。對(duì)于低水合物飽和度的情況，這兩種模式可能是分別存在的，但是對(duì)于高水合物飽和度情況，這兩種模式是同時(shí)存在的。這也就可以解釋為什么在不同賦存模式下，對(duì)強(qiáng)度的提升幅度是有區(qū)別的［19-20］。這兩類加固機(jī)制也得到其他學(xué)者的認(rèn)可，并進(jìn)一步加以應(yīng)用。Uchida等［21］分析了兩類增強(qiáng)機(jī)制對(duì)土體的加固影響程度，并且據(jù)此建立了含水合物沉積物的本構(gòu)模型。

探究水合物形成對(duì)土體強(qiáng)度的影響機(jī)理必須分析水合物形成對(duì)強(qiáng)度指標(biāo)的影響規(guī)律。根據(jù)目前的試驗(yàn)結(jié)果分析，水合物形成將增強(qiáng)土體的粘聚力強(qiáng)度，只是水合物的賦存形式對(duì)于強(qiáng)度的提升幅度具有一定的差別［15，17，19-20，22］；而水合物的形成對(duì)于內(nèi)摩擦角的影響是和賦存土性有關(guān)。Masui等［19］和顏榮濤等［22］依據(jù)三軸試驗(yàn)得出了水合物在砂土等中形成對(duì)內(nèi)摩擦角影響可以忽略的結(jié)論；而Zhang等［15］通過含水合物粉土的三軸試驗(yàn)得出了其內(nèi)摩擦角隨水合物含量的增大而降低的結(jié)論。這一結(jié)論在Waite等［17］和Yun等［20］的試驗(yàn)中也得到了證實(shí)。圖2為砂土和粉土中形成的水合物對(duì)粘聚力和內(nèi)摩擦角的影響規(guī)律。

圖2 含水合物砂土、粉土的粘聚力以及內(nèi)摩擦角［15，19，22］Fig.2 Cohesions and internal friction angles of hydrate-bearing sands and silts

2　水合物沉積物的強(qiáng)度模型

本模型在不考慮水合物懸浮于孔隙水中對(duì)強(qiáng)度有貢獻(xiàn)的情況下，主要考慮在持力體模式和膠結(jié)模式條件下含天然氣水合物沉積物的強(qiáng)度特性。

2.1前提假設(shè)

目前在對(duì)含水合物沉積物的問題進(jìn)行處理時(shí)，通常的做法是將水合物當(dāng)成土骨架的一部分進(jìn)行處理，這樣處理后含水合物沉積物的固體骨架部分可認(rèn)為是土顆粒和水合物組成的復(fù)合材料。基于這種處理方式，吳二林等［23］建立了含天然氣水合物沉積物的損傷統(tǒng)計(jì)本構(gòu)模型，楊期君等［24］建立含天然氣水合物沉積物的彈塑性損傷本構(gòu)模型，這種方法也應(yīng)用于和含水合物沉積物具有類似力學(xué)特性的凍土問題中［25］。按照以上處理后，存在以下幾種假設(shè)：

①水合物與土顆粒被認(rèn)為是兩種連續(xù)介質(zhì)材料復(fù)合而成的含水合物沉積物，二者具有相同變形。

②水合物與土顆粒共同按照體積分?jǐn)?shù)承擔(dān)上覆荷載。

③在膠結(jié)模式下，考慮膠結(jié)貢獻(xiàn)和填充貢獻(xiàn)；在持力體模式下，只考慮水合物的填充貢獻(xiàn)。

2.2模型建立

根據(jù)2.1節(jié)的前提假設(shè)，把含水合物沉積物假設(shè)為土體和水合物的混合復(fù)合材料，其中一部分水合物膠結(jié)了土顆粒使土顆粒粘聚力提高，一部分水合物成為了類似于土顆粒的持力體，共同承擔(dān)上覆荷載和剪切力。圖3給出了含水合物沉積物強(qiáng)度模型的示意圖。

圖3　含水合物沉積物強(qiáng)度模型示意圖Fig.3 Sketch of strength model on gas hydrate sediments

土體和水合物共同承擔(dān)了上覆正應(yīng)力，同時(shí)也共同承擔(dān)剪應(yīng)力?？擅枋鰹?/p>

式中：σ和τ的腳標(biāo)s、h分別代表土體和水合物，ns、n分別代表土體、土體孔隙的體積分?jǐn)?shù)，Sph為貢獻(xiàn)于填充孔隙而作為持力體的水合物的飽和度。

根據(jù)假設(shè)①，水合物和土體具有相同變形，即存在變形協(xié)調(diào)約束。因此，

式中：ε為水合物和土體共同的應(yīng)變，Es和Eh分別為土體和水合物的變形模量。把式（3）和式（4）代入式（1）中，有

土體屬于很典型的摩擦材料，符合摩爾-庫(kù)倫強(qiáng)度準(zhǔn)則，可表示為

式中：cs（Sc

h）為土體相對(duì)于土體部分面積的粘聚力；c（Sch）為土體相對(duì)于整個(gè)土體面積的粘聚力，二者之間存在ns·cs（Sch）=c（Sch）；Sch為貢獻(xiàn)于膠結(jié)強(qiáng)度的水合物的飽和度，當(dāng)水合物賦存模式為填充模式時(shí)，Sch=0。

水合物的剪切性質(zhì)類似于晶體材料，剪切強(qiáng)度受圍壓的影響不大，摩擦角一般取2°～5°［8-9，20］，這里抓住主要因素，忽略次要因素，不考慮純水合物的摩擦強(qiáng)度，可以認(rèn)定水合物的剪切強(qiáng)度τh= 0.5qh，qh為水合物的偏應(yīng)力強(qiáng)度［26］。這里需要說明的是，qh并不是真實(shí)的純水合物的偏應(yīng)力強(qiáng)度。這是因?yàn)樵谕馏w孔隙中形成的水合物微觀、細(xì)觀結(jié)構(gòu)上并不是連續(xù)的一個(gè)整體，這里的處理方式只是從宏觀上等效水合物相，與真實(shí)的純水合物具有一定的差距。

把式（8）和τh代入式（2），可以得出含水合物沉積物的強(qiáng)度模型的基本表達(dá)式：

其中，右邊第1項(xiàng)表示水合物膠結(jié)土顆粒后的粘結(jié)強(qiáng)度；第2項(xiàng)表示水合物填充孔隙后土體摩擦強(qiáng)度；第3項(xiàng)表示孔隙中填充的水合物承擔(dān)的剪切強(qiáng)度。上式從表達(dá)式看較為復(fù)雜，為了使用起來方便，有必要對(duì)其進(jìn)行簡(jiǎn)化：

其中，

很明顯，式（11）表征的是含水合物沉積物的等效摩擦系數(shù)，而式（13）表征的是含水合物沉積物的等效粘聚力。在式（11）中，含水合物沉積物的等效摩擦系數(shù)根據(jù)β值的取值不同來控制水合物含量對(duì)等效摩擦系數(shù)的影響。當(dāng)賦存介質(zhì)為砂土等粗粒土?xí)r，β=0，等效摩擦系數(shù)不受水合物含量的影響；而當(dāng)賦存為粉土或是粘土?xí)r，β≠0，可以根據(jù)水合物對(duì)等效摩擦系數(shù)的影響情況進(jìn)行取值。

由圖4可見，隨著水合物的飽和度增大，等效摩擦角逐漸變小，參數(shù)β控制變小的幅度。這種等效內(nèi)摩擦角變小的情況在一些膠結(jié)土中普遍存在，這一強(qiáng)度模型可以對(duì)此進(jìn)行很好的模擬。

式（13）中，等效粘聚力包括3部分，分別為水合物膠結(jié)作用貢獻(xiàn)的粘聚力、不含水合物沉積物的粘聚力以及水合物作為填充模式成為持力體而貢獻(xiàn)的粘聚力，這3部分很好地反映了含水合物沉積物等效粘聚力的來源機(jī)制。

圖4　等效內(nèi)摩擦角與水合物飽和度之間的關(guān)系（φ=26°）Fig.4 Relation between equivalent internal friction and saturation of hydrate

在膠結(jié)賦存模式下

這里需要說明的是，填充賦存模式下，當(dāng)水合物繼續(xù)生長(zhǎng)，處于高飽和度時(shí)，水合物也會(huì)對(duì)土體產(chǎn)生膠結(jié)作用，但是這種膠結(jié)作用的產(chǎn)生條件相當(dāng)苛刻，如必須保證土顆粒的間距以及在高飽和度的條件下，當(dāng)然這種高飽和度在自然條件下 是很 少 見的［26-27］。

2.4強(qiáng)度模型參數(shù)的確定

在強(qiáng)度模型中包括不含水合物沉積物的強(qiáng)度參數(shù)c、φ，二者可以通過不含水合物沉積物的直剪強(qiáng)度試驗(yàn)進(jìn)行確定。

此外，在填充模式和持力體賦存模式下，還需要確定模型參數(shù)χ、β。這兩個(gè)參數(shù)可以通過至少兩個(gè)不同飽和度的含水合物沉積物的直剪強(qiáng)度試驗(yàn)得出等效摩擦角和等效粘聚力后，進(jìn)行擬合得出。

而對(duì)于膠結(jié)賦存模式，需要確定χ、β、α、γ等4個(gè)模型參數(shù)。這些參數(shù)也至少利用4個(gè)不同水合物飽和度的含水合物沉積物，通過與上述類似的方法，在不同豎向正壓力下的直剪強(qiáng)度試驗(yàn)得出。

這里需要強(qiáng)調(diào)說明的是，模型參數(shù)β和χ都是具有一定的物理表征意義的，可以按照相應(yīng)的物理表征公式進(jìn)行計(jì)算得出（式（12））。然而，由于本模型的建立是將宏觀上的水合物材料相假設(shè)為等效的水合物相，因此水合物強(qiáng)度與真實(shí)塊體的水合物強(qiáng)度存在一定差別，并且會(huì)隨賦存土體的不同而發(fā)生改變。因此，為了保證強(qiáng)度模擬的預(yù)測(cè)結(jié)果，這些參數(shù)應(yīng)該通過試驗(yàn)數(shù)據(jù)擬合得出。此外，對(duì)于通過擬合試驗(yàn)數(shù)據(jù)而獲取的強(qiáng)度參數(shù)，應(yīng)在條件滿足的情況下采用更多的數(shù)據(jù)進(jìn)行擬合確定模型參數(shù)，以便使模型獲得更好的精度。

3　強(qiáng)度模型應(yīng)用及分析

以上強(qiáng)度模型是以摩爾-庫(kù)倫強(qiáng)度準(zhǔn)則為基礎(chǔ)推導(dǎo)出來的，把其換算成三軸試驗(yàn)條件下的強(qiáng)度，可以表示為

為了驗(yàn)證以上強(qiáng)度模型的有效性，收集了前人的試驗(yàn)數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證，圖5中給出了Ebinuma等［28］對(duì)于含水合物砂土的強(qiáng)度試驗(yàn)數(shù)據(jù)及其強(qiáng)度模型的模擬。其中，對(duì)于膠結(jié)模式，很小含量的水合物就可以使土體強(qiáng)度有很明顯的提高；而對(duì)于填充模式，當(dāng)水合物含量低于臨界水合物飽和度時(shí)，水合物的形成對(duì)于強(qiáng)度影響很小，可以忽略，當(dāng)水合物飽和度超過臨界飽和度時(shí)，土體強(qiáng)度受水合物含量的增大有明顯的增加。通過對(duì)比發(fā)現(xiàn)該強(qiáng)度模型對(duì)此能有很好的反映。圖6給出了顏榮濤等［22］對(duì)不同水合物飽和度含水合物砂土在1、2和3 MPa圍壓下的剪切強(qiáng)度以及強(qiáng)度模型的模擬曲線，水合物的賦存模式屬于膠結(jié)模式。強(qiáng)度模型的模擬和試驗(yàn)數(shù)據(jù)的對(duì)比表明在同一套參數(shù)下，強(qiáng)度模型也能很好地反映圍壓對(duì)強(qiáng)度的影響。

以上的模型應(yīng)用都是針對(duì)賦存介質(zhì)為砂土等粗粒土的情況，含水合物沉積物的內(nèi)摩擦角不受水合物的形成和含量的影響。在圖7中，利用本文提出的強(qiáng)度模型模擬了文獻(xiàn)［15］中含甲烷水合物粉土的摩擦系數(shù)（內(nèi)摩擦角的正切值）和粘聚力隨水合物含量的變化［15］。通過與試驗(yàn)數(shù)據(jù)的對(duì)比，發(fā)現(xiàn)模型能很好地反映隨水合物飽和度的增加，含水合物粉土的摩擦系數(shù)降低和粘聚力上升的趨勢(shì)。但是模擬結(jié)果與試驗(yàn)結(jié)果還是存在一定的差距，這是因?yàn)樵谖墨I(xiàn)［15］的試驗(yàn)中，含水合物粉土孔隙里除了水合物外，還存在冰的膠結(jié)，從而影響了結(jié)果。

圖5　強(qiáng)度模擬　?。═.Ebinuma，2005）［28］Fig.5 Stength modelling

圖6　強(qiáng)度模擬　?。仒s濤等，2012）［22］Fig.6 Stength modelling

圖7　含水合物粉土摩擦系數(shù)及粘聚力的模擬（Zhang et al，2012）［15］Fig.7 Simulations of friction coefficients and cohesions of hydrate-bearing silts

從以上強(qiáng)度模型模擬的3組試驗(yàn)強(qiáng)度案例來看，該強(qiáng)度模型能很好地反映水合物賦存形式、水合物含量、圍壓對(duì)含水合物沉積物強(qiáng)度的影響。同時(shí)還能有效的考慮當(dāng)賦存介質(zhì)為細(xì)粒土?xí)r，含水合物沉積物內(nèi)摩擦角變小的特征。然而，有些情況下模型模擬的強(qiáng)度值與試驗(yàn)值還是存在較大的差距。這一問題主要有兩方面的原因：（1）該模型理想化水合物為非摩擦材料，內(nèi)摩擦角為零；（2）由于試驗(yàn)數(shù)據(jù)的欠缺，等效粘聚力的經(jīng)驗(yàn)公式的設(shè)定有所欠缺。

4　結(jié) 論

本文在調(diào)研前人對(duì)含水合物沉積物強(qiáng)度特性的研究基礎(chǔ)上，分析了不同賦存模式下水合物存在對(duì)含水合物沉積物的強(qiáng)度加固機(jī)理，得出水合物對(duì)賦存土體強(qiáng)度影響機(jī)理主要有兩種：第1種是水合物作為持力體，分擔(dān)了土骨架的一部分應(yīng)力，相當(dāng)于在土體孔隙中添加水合物顆粒來提高土體的密實(shí)度，從而提高其抗剪強(qiáng)度；第2種是土顆粒之間的膠結(jié)作用，這種模式能在很大程度上提高土體的粘聚力。

基于水合物賦存對(duì)含水合物沉積物的強(qiáng)度影響機(jī)理，把含水合物沉積物視為復(fù)合材料進(jìn)行處理，理論上推導(dǎo)出了含水合物沉積物的強(qiáng)度模型。該模型把含水合物沉積物分為土顆粒骨架相和水合物相，二者具有相同的變形。同時(shí)，土顆粒骨架為典型的摩擦材料，符合傳統(tǒng)的摩爾-庫(kù)倫強(qiáng)度準(zhǔn)則，而水合物相理想為晶體類材料，忽略內(nèi)摩擦角影響。以此為基礎(chǔ)推導(dǎo)出的強(qiáng)度模型能較好的模擬含水合物沉積物強(qiáng)度的一些特性，如不同賦存模式、不同水合物含量、不同圍壓對(duì)含水合物沉積物的強(qiáng)度的影響趨勢(shì)，以及細(xì)粒土中水合物的形成將導(dǎo)致內(nèi)摩擦角變小的特性等。

該模型是以試驗(yàn)揭示的影響機(jī)理為基礎(chǔ)，從理論上推導(dǎo)出的強(qiáng)度模型，因此具有很好的擴(kuò)展適用性。但由于模型中采用了一些簡(jiǎn)化措施以及強(qiáng)度試驗(yàn)數(shù)據(jù)的缺乏，一些關(guān)系式還需要進(jìn)行優(yōu)化，將在今后的工作中進(jìn)一步研究完善。

［1］錢伯章，朱建芳.天然氣水合物：巨大的潛在能源 ［J］.天然氣與石油，2008，26（4）：47-52.

［2］江懷友，喬衛(wèi)杰，鐘太賢，等.世界天然氣水合物資源勘探開發(fā)現(xiàn)狀與展望［J］.中外能源，2008，13（6）：19-25.

［3］吳時(shí)國(guó)，陳珊珊，王志君，等.大陸邊緣深水區(qū)海底滑坡及其不穩(wěn)定性風(fēng)險(xiǎn)評(píng)估［J］.現(xiàn)代地質(zhì)，2008，22（3）：430-437.

［4］Davie M K，Buffett B A.Sources of methane for marine gas hydrate：Inferences from a comparison of observations and numerical models［J］.Earth and Planetary Science Letters，2003，206（1）：51-63.

［5］Milkov A V，Dickens G R，Claypool G E，et al.Co-existence of gas hydrate，free gas，and brine within the regional gas hydrate stability zone at Hydrate Ridge（Oregon margin）：Evidence from prolonged degassing of a pressurized core［J］. Earth and Planetary Science Letters，2004，222（3）：829-843.

［6］Yu F，Song Y C，Liu W G，et al.Analyses of stress strain behavior and constitutive model of artificial methane hydrate［J］.Journal of Petroleum Science and Engineering，2011，77（2）：183-188.

［7］Miyazaki K，Tenma N，Aoki K，et al.Effects of confining pressure on mechanical properties of artificial methane-hydratebearing sediment in triaxial compression test［J］.International Journal of Offshore and Polar Engineering，2011，21（2）：148-154.

［8］Miyazaki K，Masui A，Sakamoto Y，et al.Triaxial compressive properties of artificial methane-hydrate-bearing sediment［J］.Journal of Geophysical Research，2011，116：B06102，doi：10.1029/2010 JB008049.

［9］Song Y C，Yu F，Li Y H，et al.Mechanical property of artificial methane hydrate under triaxial compression［J］.Journal of Natural Gas Chemistry，2010，19（3）：246-250.

［10］Lee J Y，Santamarina J C，Ruppel C.Mechanical and electromagnetic properties of northern Gulf of Mexico sediments with and without THF hydrates［J］.Marine and Petroleum Geology，2008，25（9）：884-895.

［11］Hyodo M，Nakata Y，Yoshimoto N，et al.Shear strength of methane hydrate bearing sand and its deformation during dissociation of methane hydrate［C］//Proceedings of the 4th International Symposium on Deformation Characteristics of Geomaterials.2008：549-556.

［12］Hyodo M，Nakata Y，Yoshimoto N，et al.Shear behavior of methane hydrate-bearing sand［C］//The Seventeenth International Offshore and Polar Engineering Conference.International Society of Offshore and Polar Engineers，2007：1326 -1333.

［13］Hyodo M，Nakata Y，Yoshimoto N，et al.Basic research on the mechanical behavior of methane hydrate-sediments mixture［J］.Soils and Foundations，2005，45（1）：75-85.

［14］Durham W B，Kirby S H，Stern L A，et al.The strength and rheology of methane clathrate hydrate［J］.Journal of GeophysicalResearch，2003， 108（B4），2182， doi：10.1029/2002JB 001872.

［15］Zhang X H，Lu X B，Zhang L M，et al.Experimental study on mechanical properties of methane-hydrate-bearing sediments［J］.Acta Mechanica Sinica，2012，28（5）：1356-1366.

［16］Li Y H，Song Y C，Liu W G，et al.A new strength criterion and constitutive model of gas hydrate-bearing sediments under high confining pressures［J］.Journal of Petroleum Science and Engineering，2013，109：45-50.

［17］Waite W F，Santamarina J C，Cortes D D，et al.Physical properties of hydrate-bearing sediments［J］.Reviews of Geophysics，2009，47（4）：1-38.

［18］Dvorkin J，Helgerud M B，Waite W F，et al.Introduction to physical properties and elasticity models［M］//Natural Gas Hydrate.Springer Netherlands，2003：245-260.

［19］Masui A，Haneda H，Ogata Y，et al.Effects of methane hydrate formation on shear strength of synthetic methane hydrate sediments［C］//Proceedings of the 15th International Off-shore and Polar Engineering Conference.Seoul：The International Society of Offshore and Polar Engineers，2005：364 -369.

［20］Yun T S，Santamarina J C，Ruppel C.Mechanical properties of sand，silt，and clay containing tetrahydrofuran hydrate［J］.JournalofGeophysicalResearch，2007，112：B04106，doi：10.1029/2006JB 004484.

［21］Uchida S，Soga K，Yamamoto K.Critical state soil constitutive model for methane hydrate soil［J］.Journal of Geophysical Research，2012，117：B03209，doi：10.1029/2011JB 008661.

［22］顏榮濤，韋昌富，魏厚振，等.水合物形成對(duì)含水合物砂土強(qiáng)度影響 ［J］.巖土工程學(xué)報(bào)，2012，34（7）：1234-1240.

［23］吳二林，韋昌富，魏厚振，等.含天然氣水合物沉積物損傷統(tǒng)計(jì)本構(gòu)模型 ［J］.巖土力學(xué)，2013，34（1）：60-65.

［24］楊期君，趙春風(fēng).含氣水合物沉積物彈塑性損傷本構(gòu)模型探討［J］.巖土力學(xué)，2014，35（4）：991-997.

［25］寧建國(guó)，王慧，朱志武，等.基于細(xì)觀力學(xué)方法的凍土本構(gòu)模型研究［J］.北京理工大學(xué)學(xué)報(bào)，2005，25（10）：847-851.

［26］Santamarina J C，Ruppel C.The impact of hydrate saturation on the mechanical，electrical，and thermal properties of hydrate-bearing sand，silts，and clay［J］.Geophysical Characterization of Gas Hydrates，Geophys.Dev.Ser，2010，14：373-384.

［27］蔣明鏡，賀潔，周雅萍.考慮水合物膠結(jié)厚度的深海能源土粒間膠結(jié)模型研究 ［J］.巖土力學(xué)，2014，35（5）：1231-1240.

［28］Ebinuma T，Kamata Y，Minagawa H，et al.Mechanical properties of sandy sediment containing methane hydrate［C］//Fifth International Conference on Gas Hydrates. 2005：958-961.

Strength model for methane hydrate-bearing soil

YAN Rong-tao，ZHAO Xu-yue，YANG De-huan，XIAO Gui-yuan，LIANG Wei-yun，WEN Song-song
（Guangxi Key Laboratory of Geomechanics and Geotechnical Engineering，Guilin University of Technology，Guilin 541004，China）

The strength model of methane hydrate-bearing soil is one of the vital problems to be solved in hydrate exploiting.In this paper，based on previous tests of methane hydrate-bearing soils on strength，two reinforcing mechanisms are analyzed.The reinforcing mechanisms include hydrate transfer into bearing body which can load part stress from soil frame and cementing action between soil particles.Based on the above analyses，methane hydrate-bearing soil is considered to be composite material，which consists of soil particle phase and hydrate phase，then a strength model is established.The strength model can predict the strength values in condition of various occurrence models，various saturations and various confining pressures for methane hydratebearing soils satisfactorily.Meanwhile，the dropping characteristics of friction angles of fine-grained soil with hydrate are taken into account as increasing saturation of hydrate.However，in short test data on strength，some formulas are not good enough and should be optimized in the future works.

hydrate-bearing soil；occurrence model；strength model

TU411

A

1674-9057（2016）03-0514-07

10.3969/j.issn.1674-9057.2016.03.015

2015-01-21

國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目（51309055；11372078；11562007）；廣西自然科學(xué)基金項(xiàng)目 （2014GXNSFBA118236）

顏榮濤（1984—），男，博士，副教授，研究方向：天然氣水合物開采過程中關(guān)鍵巖土工程問題，yrt301@163.com。

引文格式：顏榮濤，趙續(xù)月，楊德歡，等.天然氣水合物沉積物的強(qiáng)度模型［J］.桂林理工大學(xué)學(xué)報(bào)，2016，36（3）：514-520.