李洋輝 宋永臣 劉衛(wèi)國 于 鋒

(“海洋能源利用與節(jié)能”教育部重點實驗·大連理工大學(xué)室)

0 前言

天然氣水合物作為一種極具潛質(zhì)的清潔能源，具有儲量大、能量密度高、分布廣等特點，其安全開采已成為石油天然氣工業(yè)新的研究熱點[1，2]。從凍土或海底含天然氣水合物沉積層中采收天然氣會使沉積層強度降低，增加地層的不穩(wěn)定性，可能引起地質(zhì)塌陷，海底滑坡等災(zāi)害。溢出的甲烷氣體具有很強的溫室效應(yīng)，是全球氣候變暖的潛在影響因素[3-5]。在勘探與開采過程中，鉆井作業(yè)會擾動含水合物沉積層的地質(zhì)結(jié)構(gòu)，這可能導(dǎo)致水合物的分解并削弱沉積層的地質(zhì)力學(xué)強度，造成鉆探設(shè)備、海底管道等基礎(chǔ)設(shè)施的破壞及生命財產(chǎn)損失。研究含水合物沉積層的機械力學(xué)強度對天然氣水合物的安全開采具有重要意義。

目前有關(guān)水合物沉積物力學(xué)性質(zhì)的研究很少有報道。Parameswaran等[5]和Cameron等[7]用四氫呋喃水合物代替天然氣水合物研究了應(yīng)變速率對沉積物機械強度的影響。Winters等[8,9]通過聲波測量、三軸壓縮試驗得到水合物沉積物原狀試樣和實驗室人工試樣的應(yīng)力-應(yīng)變曲線和強度等的初步認(rèn)識。Hyodo等[10]和Vanoudheusden等[3]分別對天然氣水合物砂樣和非飽和海洋沉積土進行了一系列的室內(nèi)三軸壓縮試驗，分析了水合物沉積物的力學(xué)性質(zhì)與飽和度、溫度、圍壓、孔隙壓力的關(guān)系。但由于之前的研究一般采用砂土作為水合物沉積物骨架，與凍土區(qū)域或海底含水合物沉積物的粒徑分布及成份組成有很大不同。且四氫呋喃水合物與天然氣水合物相平衡條件差異較大，不適合用來研究溫度對其力學(xué)性質(zhì)的影響。為了更好的模擬含水合物沉積層，評估地層的穩(wěn)定性，本文采用與天然含水合物沉積物參數(shù)更接近的高嶺土和天然氣水合物混合制備試樣。通過在實驗室人工合成天然氣水合物巖芯來模擬含水合物沉積層進行三軸壓縮試驗，分析了溫度、應(yīng)變速率對沉積物強度的影響。

1 實驗設(shè)備及實驗方法

試驗采用TAW-60低溫高壓三軸試驗機(圖1)。本試驗機能夠模擬水合物沉積層原位的壓力和溫度環(huán)境。溫度控制范圍-20℃～25℃，最大圍壓30MPa，最大軸向載荷60KN，試樣尺寸Φ50×100mm。壓縮實驗過程中的圍壓維持和軸向載荷由兩個擴展數(shù)字控制器(EDC，德國DOLI公司)分別自動控制。恒溫槽通過熱交換器控制液壓油的溫度，進而調(diào)節(jié)試樣溫度。通過在壓力室內(nèi)布置熱電偶(精度±0.5℃)和壓力傳感器(精度0.1MPa)分別測量和反饋試驗過程中的試樣溫度和壓力參數(shù)。

原狀水合物沉積物巖芯獲取難度較大、成本高，而實驗室人工合成水合物沉積物試樣具有參數(shù)可控性、成本低[11]。研究表明合成樣和原狀樣力學(xué)特性的不同主要是由于初始孔隙比和顆粒級配的不同造成的[12]，所以只要控制合成樣的初始孔隙比和顆粒級配，可以得到與原狀樣接近的實驗結(jié)果。我國青藏高原凍土區(qū)沉積層成分主要為粘土，粒徑較小，且顆粒級配良好[13,14]。本實驗采用烘干后的高嶺土模擬凍土區(qū)含水合物沉積層骨架。所選高嶺土與目標(biāo)沉積物具有相似的粒徑分布及級配，且在一定壓力下成型可以得到與目標(biāo)沉積物接近的初始孔隙比。將冰粉和純甲烷氣體在高壓反應(yīng)釜內(nèi)混合，控制溫度、壓力條件(-5℃, 13MPa)，維持足夠長的時間，生成天然氣水合物。然后將其與高嶺土按一定比例混合，在壓力成型裝置中壓縮成圓柱試樣(直徑為50mm，高度為100mm)，控制溫度-10℃，軸向載荷10MPa。本實驗沉積物試樣水合物飽和度為30%，為水合物礦藏飽和度的平均值。采用高嶺土體積比為40%，試樣密度1.63g/cm3,干密度1.03g/cm3。在試樣表面套厚度為0.5mm的橡皮膜，用以隔絕孔隙壓力和圍壓。剪切試驗開始前，試樣固結(jié)2h，使試樣內(nèi)部溫度分布均勻。試驗采用UU剪切，根據(jù)SL 237-017-1999土工試驗規(guī)程，選取試驗參數(shù)(表1)。

圖1 TAW-60低溫高壓三軸試驗機

表1 天然氣水合物三軸壓縮試驗參數(shù)選取

2 結(jié)果和討論

三軸壓縮試驗過程中，水合物沉積物試樣都表現(xiàn)出明顯的塑性變形現(xiàn)象，圓柱試樣中部凸起呈鼓狀，沒有出現(xiàn)斷面現(xiàn)象(破壞面)。本文主要分析溫度和應(yīng)變速率對沉積物試樣力學(xué)性質(zhì)的影響。

2.1 溫度對天然氣水合物沉積物力學(xué)性質(zhì)的影響

天然氣水合物在自然界處于亞穩(wěn)定狀態(tài)，對溫度變化比較敏感。在一定壓力下，溫度的變化不僅能決定天然氣水合物的生成與分解，還能改變沉積物的膠結(jié)狀態(tài)、含水量、孔隙結(jié)構(gòu)等，從而改變其力學(xué)性質(zhì)。

從沉積物試樣在圍壓5MPa，不同溫度條件下的應(yīng)力應(yīng)變曲線(圖2)?？梢钥吹狡珣?yīng)力隨著軸向應(yīng)變的增大而增大，并逐漸趨于穩(wěn)定，期間并未出現(xiàn)明顯峰值。在壓縮試驗初期，不同溫度條件下的應(yīng)力應(yīng)變曲線基本重合，說明試樣在小應(yīng)變階段的變形特性受溫度影響不大。在壓縮試驗的末段，試樣有明顯的硬化現(xiàn)象。此外，溫度對試樣的破壞強度有顯著的影響，溫度越低，試樣強度越大。

圖2 不同溫度條件下水合物沉積物試樣應(yīng)力應(yīng)變曲線

最大偏應(yīng)力定義為軸向應(yīng)變到達(dá)15%之前偏應(yīng)力的峰值，從試樣最大偏應(yīng)力隨溫度發(fā)生變化的關(guān)系曲線(圖3)，可以看到最大偏應(yīng)力隨著溫度的降低而增大，在一定溫度范圍內(nèi)可用線性表示。隨著溫度的進一步降低，最大偏應(yīng)力的增加趨勢有所減小，在溫度足夠低時，將趨于定值。Hyodo[10]等利用粒徑較大的豐浦砂模擬含水合物沉積層，研究了5℃、-10℃、-30℃條件下的破壞強度，得到了類似的實驗結(jié)果。在0℃以下時，含水合物沉積層強度隨著溫度的降低線性增加，且在實驗過程中，試樣應(yīng)力應(yīng)變曲線未出現(xiàn)明顯峰值。本文通過增加實驗組，進一步分析了溫度對沉積層破壞強度及變形特性的影響。

剛度是材料抵抗變形的能力，一般可用初始彈性模量E0表示材料的剛度。但是由于E0比較難測定，且受實驗誤差的影響較大，試樣在變形過程中的切線模量也不是定值，因此一般用割線模量E50作為試樣的平均剛度。從溫度與沉積物割線模量的關(guān)系曲線(圖4)，可以看到割線模量隨著溫度的降低線性增大。

大量的實驗研究結(jié)果顯示，凍土中的液態(tài)水在凍結(jié)過程中并沒有全部轉(zhuǎn)化為固態(tài)冰，始終保持一定數(shù)量的未凍水[15]。且由于壓融現(xiàn)象的存在，在高壓下固結(jié)時，顆粒間的冰粒會出現(xiàn)部分融化，導(dǎo)致冰點以下仍可能存在液態(tài)水。這種現(xiàn)象同樣在水合物沉積物試樣中存在。研究表明未凍水含量與負(fù)溫絕對值呈冪函數(shù)關(guān)系，未凍水含量隨著溫度的降低而減小。當(dāng)溫度低于-15℃時，未凍水含量隨溫度的減小速率減慢，并逐漸趨于定值[16]。液態(tài)水在顆粒間能起到潤滑的作用，降低剪切過程中的摩擦阻力，進而使沉積物試樣強度降低。固態(tài)冰能夠增加顆粒間的膠結(jié)作用力，且固態(tài)冰作為試樣骨架能夠承受外部載荷，固態(tài)冰的存在，能夠增加試樣強度。溫度對水合物沉積物強度的影響主要是由于溫度變化改變了試樣中的未凍水含量及固態(tài)冰含量。當(dāng)溫度降低時，含水量降低，降低了水在土顆粒間的潤滑作用，增加了摩擦阻力。同時，含水量的減少增加了冰與土顆粒的膠結(jié)面積，增加了試樣顆粒與冰粒間的膠結(jié)作用力。且冰的結(jié)構(gòu)也隨著溫度的降低向更緊密堅硬的有序結(jié)構(gòu)發(fā)展，增大了試樣抵抗變形的能力。正是由于未凍水含量與溫度的冪函數(shù)關(guān)系，造成了圖4曲線斜率隨著溫度的降低而逐漸減小的現(xiàn)象，在-15℃時斜率變化尤其明顯。

圖3 溫度對水合物沉積物強度的影響

圖4 溫度對水合物沉積物割線模量的影響

剛度受顆粒間相互作用、顆粒接觸面積以及摩擦阻力的影響。溫度的降低增加了土顆粒、冰、水合物之間的膠結(jié)作用力，增加了顆粒間的膠結(jié)面積以及摩擦阻力，因此割線模量E50隨著溫度的降低而增大。

2.2 應(yīng)變速率對天然氣水合物沉積物力學(xué)性質(zhì)的影響

應(yīng)變速率會影響材料的屈服強度和失效模式(圖5、圖6)，試樣的強度隨應(yīng)變速率的增加有較明顯的增加趨勢，而在軸向應(yīng)變達(dá)到0.5%之前，試樣的應(yīng)力應(yīng)變曲線幾乎重合。另外，應(yīng)變速率可明顯的改變材料的塑性性能。當(dāng)應(yīng)變速率為1.0%/min和0.5%/min時，偏應(yīng)力隨著應(yīng)變的增大而增大，期間沒有峰值出現(xiàn)，呈現(xiàn)明顯的硬化趨勢。當(dāng)應(yīng)變速率為0.1%/min時，偏應(yīng)力隨著應(yīng)變先增大后減小，并逐漸趨于定值，表現(xiàn)出軟化現(xiàn)象。這種現(xiàn)象與Baker[17]、馬巍等[18]在冰和凍土方面的實驗研究結(jié)果一致。Parameswaran等[6]和Cameron等[7]用四氫呋喃水合物代替天然氣水合物研究應(yīng)變速率對沉積物破壞強度的影響，也得到類似的實驗結(jié)果。但是含四氫呋喃水合物沉積物的應(yīng)力應(yīng)變曲線呈現(xiàn)明顯的峰值，存在軟化和屈服現(xiàn)象，這可能是四氫呋喃水合物和天然氣水合物自身物理性質(zhì)的不同造成的。

圖5 應(yīng)變速率對合物沉積物應(yīng)力應(yīng)變曲線的影響

圖6 應(yīng)變速率對水合物沉積物強度的影響

Grady[19]基于巖石動態(tài)力學(xué)特性的實驗研究，認(rèn)為材料內(nèi)部存在的裂紋擴展和聚合是材料破壞的主要原因。應(yīng)變速率較低時，只有那些能在低應(yīng)力水平下被激活的裂紋發(fā)生擴展，這些裂紋的擴展和聚合使得材料在應(yīng)力水平達(dá)到能使其他裂紋擴展之前已經(jīng)發(fā)生破壞，因此材料具有較低的強度。而在高應(yīng)變速率下，在低應(yīng)力水平下被激活的裂紋聚合之前，應(yīng)力已經(jīng)達(dá)到一個較高的水平，這時需要更多的裂紋參與擴展，消耗外力功，從而導(dǎo)致材料強度的增加。

3 結(jié)論

(1)溫度對水合物沉積物試樣的破壞強度有顯著影響，溫度越低，試樣強度越大，在一定溫度范圍內(nèi)可用線性表示。隨著溫度的進一步降低，破壞強度的增加趨勢有所減小，在溫度足夠低時，將趨于定值。

(2)水合物沉積物試樣的破壞強度隨著應(yīng)變速率的增大而增大。另外，應(yīng)變速率可明顯的改變材料的塑性性能。當(dāng)應(yīng)變速率為0.5%/min和1.0%/min時，試樣的應(yīng)力應(yīng)變曲線有明顯的硬化現(xiàn)象，而當(dāng)應(yīng)變速率為0.1%/min時，存在軟化現(xiàn)象。

(3)割線模量E50隨著溫度的降低線性增加。

1 Sangyong Lee, Gerald D. Holder. Methane Hydrates Potential as a Future Energy Source [J]. Fuel Processing Technology, 2001, 71(1-3):181-186.

2 李棟梁, 樊栓獅. 天然氣水合物資源開采方法研究[J]. 化工學(xué)報, 2003, 54(增刊): 108-112.

3 Vanoudheusden E, Sultan N, Cochonat P. Mechanical behaviour of unsaturated marine sediments: experimental and theoretical approaches [J]. Marine Geology, 2004, 213(1-4): 323-342.

4 李洋輝, 宋永臣, 劉衛(wèi)國. 天然氣水合物三軸壓縮試驗研究進展[J]. 天然氣勘探與開發(fā). 2010, 33(2): 51-55.

5 Glasby G P. Potential impact on climate of the exploitation of methane hydrate deposits offshore [J]. Marine and Petroleum Geology, 2003, 20(2003): 163-175.

6 Parameswaran V R, Paradis M, Handa Y P. Strength of frozen sand containing tetrahydrofuran hydrate [J]. Canadian Geotechnical Journal, 1989, 26(3): 479-483.

7 Cameron I, Handa Y P, Baker T H. Compressive strength and creep behavior of hydrate-consolidated sand [J]. Canadian Geotechnical Journal, 1990, 27(2): 255-258.

8 Winters W J, Waite W F, Mason D H, et al. Methane gas hydrate effect on sediment acoustic and strength properties [J]. Journal of Petroleum Science and Engineering, 2007, 56(1-3): 127-135.

9 Waite W F, Kneafsey T J, Winters W J, et al. Physical property changes in hydrate-bearing sediment due to depressurization and subsequent repressurization [J]. Journal of Geophysical Research-Solid Earth, 2008, 113(B7): B07102.

10 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.

11 王淑云, 魯曉兵. 水合物沉積物力學(xué)性質(zhì)的研究現(xiàn)狀[J]. 力學(xué)進展, 2009, 39(2): 176-188.

12 Akira Masui?, Kuniyuki Miyazaki, Hironori Haneda, et al. Mechanical characteristics of natural and artificial gas hydrate bearing sediments [J]. Proceedings of the 6th International Conference on Gas Hydrates (ICGH 2008), Canada, 2008.

13 馬巍, 朱元林, 馬文婷, 等. 凍結(jié)粘性土的變形分析[J].冰川凍土, 2000, 22(1): 43-47.

14 胡偉, 齊吉琳, 馬巍. 凍土前期固結(jié)壓力的試驗研究[J]. 巖土力學(xué), 2009, 30(5): 1313-1316.

15 Bouyocous G J. The freezing point method as a new means of measuring the concent ration of the soil solution directly in the soil [J]. Michigan Agricultural Experiment Station Technical Bulletin, 1916, 24: 1-44.

16 徐敩祖, 奧利奮特 J L, 泰斯 A R. 土水勢、未凍水含量和溫度[J]. 冰川凍土, 1985, 7(1): 1-11.

17 Baker T H W. Strain rate effect on the compressive strength of frozen sand [J]. Engineering Geology, 1979, 13(1-4): 223-231.

18 馬巍, 吳紫汪, 常小曉, 等. 高圍壓下凍結(jié)砂土的強度特性[J]. 冰川凍土, 1996, 18(3): 268-272.

19 Grady D E. The mechanics of fracture under high-rate stress loading [J]. Mechanics of Geomaterials, 1985, 129-155.