非成巖弱膠結(jié)天然氣水合物沉積物彈性力學(xué)參數(shù)的實驗分析

趙 軍 向薪燃 趙金洲 李海濤 郭 平 張烈輝周守為 魏 納 李清平 龐維新

1. “油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學(xué) 2.中海油研究總院有限責(zé)任公司

0 引言

天然氣水合物（以下簡稱水合物）在低溫高壓的環(huán)境中才能保持穩(wěn)定，當(dāng)儲層環(huán)境發(fā)生變化，就會導(dǎo)致水合物的大量分解，排出大量溫室氣體，甚至還會導(dǎo)致海底大范圍沉降、失穩(wěn)滑坡等工程地質(zhì)災(zāi)害，給人類生活、生產(chǎn)活動帶來巨大的威脅。我國南海90%以上的水合物都屬于非成巖礦體，開采難度極大，所以對其力學(xué)特征的研究，一直都是熱點和重點。目前，由于對水合物的原位取心，技術(shù)尚未成熟，取心成本較高，故實驗室人工合成水合物沉積物來研究物理性質(zhì)是最主要的研究手段，前人也做了不少工作[1-3]。

Winters等[4]利用美國地質(zhì)調(diào)查局的水合物和沉積物模擬實驗裝置（GHASTI低溫高壓水合物三軸儀）測量了3種不同的水合物沉積物的剪切強度；Masui、Hyodo、Miyazaki等[5-7]利用人工合成水合物在不同溫壓條件下進行了三軸壓縮試驗；蔣明鏡[8-9]運用離散元方法，從微觀角度對水合物的力學(xué)特性進行研究，微觀物理模型模擬獲得的三軸壓縮強度與室內(nèi)試驗結(jié)果較為符合；中國科學(xué)院王秀娟[10]利用巖石物理模型并結(jié)合室內(nèi)實驗，對砂目數(shù)相同孔隙飽和度不同、不同水合物類型同一孔隙飽和度以及在不同孔隙飽和度下不同種類和不同粒徑砂的水合物沉積物性質(zhì)規(guī)律進行了研究；但是，利用聲波資料解釋水合物沉積物力學(xué)參數(shù)的研究比較少。利用聲波資料解釋水合物儲層的彈性力學(xué)參數(shù)可解決水合物儲層取心難、成本高的問題，并可以利用測井技術(shù)解釋連續(xù)的井下儲層力學(xué)參數(shù)，故極具研究意義。筆者使用SHW-Ⅲ型水合物聲電力學(xué)測試裝置，利用3種不同粒徑的石英砂，人工合成水合物沉積物，并測得沉積物孔隙中的水合物的飽和度以及不同骨架粒徑、不同泥質(zhì)含量、不同三軸軸壓下的聲波數(shù)據(jù)，再根據(jù)彈性理論，計算彈性力學(xué)參數(shù)，研究不同條件對水合物沉積物彈性力學(xué)參數(shù)的影響規(guī)律。

1 實驗原理

1.1 水合物沉積物的制備

本實驗采用SHW-Ⅲ型水合物聲電力學(xué)測試裝置，在恒溫1 ℃的實驗條件下，用分階段持續(xù)供給天然氣法[11]，先分別制備巖心規(guī)格為?60 mm×20 mm的3種不同骨架粒徑的水合物沉積物；再合成相同規(guī)格、骨架粒徑為0.125～0.600 mm且含有3種不同泥質(zhì)含量的沉積物試件，并對其施加5 MPa的圍壓和10 MPa、15 MPa、20 MPa的軸壓，分別測試各個試件巖電、聲波數(shù)據(jù)，計算出不同沉積物試件的水合物飽和度和動態(tài)彈性參數(shù)，分析不同條件下對動態(tài)彈性參數(shù)的影響。水合物制備步驟如下。

1.1.1 加砂

將處理好的粒徑分別為小于0.125 mm、介于0.125～ 0.600 mm和0.600～1.180 mm的3種石英砂（圖1）砂樣加入反應(yīng)釜，設(shè)置好儀器開關(guān)，連接相應(yīng)的線路。

1.1.2 建立水合物生成的溫度、壓力條件

通過圍壓系統(tǒng)及軸壓系統(tǒng)，設(shè)置壓力環(huán)境，利用循環(huán)制冷系統(tǒng)，降低反應(yīng)釜溫度。

1.1.3 加氣反應(yīng)形成水合物

向計量室內(nèi)充入8～10 MPa的甲烷氣體，記錄反應(yīng)釜入口壓力、電阻率、縱橫波的時差和波速等參數(shù)；重復(fù)加氣反應(yīng)步驟，使反應(yīng)釜中水全部生成水合物。

1.2 動態(tài)彈性參數(shù)的計算

試樣在靜態(tài)荷載條件下利用應(yīng)力—應(yīng)變曲線可求取彈性參數(shù)，這些彈性參數(shù)稱為靜態(tài)彈性參數(shù)。根據(jù)彈性波動理論，彈性波在固體介質(zhì)中傳播，它在介質(zhì)中傳播時與介質(zhì)相互作用，使接收波中攜帶了與巖石物理性質(zhì)相關(guān)的各種信息，所以也可利用聲波數(shù)據(jù)求取彈性參數(shù)，而這樣求得的彈性參數(shù)稱為動態(tài)彈性參數(shù)，動態(tài)彈性參數(shù)與靜態(tài)彈性參數(shù)之間存在一定的相關(guān)性，故可利用動彈性參數(shù)求取連續(xù)的靜彈性參數(shù)。

圖1 不同粒徑的松散石英砂樣圖

根據(jù)彈性波動理論，在假定巖石均質(zhì)、各向同性的線彈性體的前提下，利用所測定的縱波速度、橫波速度與試件的密度，便可計算出含水合物試件的動態(tài)泊松比、動態(tài)楊氏模量、動態(tài)剪切模量[12]。即

式中μ表示動態(tài)泊松比，無量綱；vs表示橫波波速，m/s；vp表示縱波波速，m/s； E表示動態(tài)楊氏模量，GPa；ρ表示巖石密度，g/cm3；G表示動態(tài)剪切模量，GPa。

2 實驗結(jié)果與分析

依據(jù)記錄的實驗參數(shù)，結(jié)合理論推導(dǎo)出來的計算公式，即可求得反應(yīng)溫度恒為1 ℃時、水合物在多孔介質(zhì)中的形成過程中3種不同骨架粒徑、3種不同泥質(zhì)含量、3種不同軸壓條件下的水合物沉積物的彈性力學(xué)參數(shù)，如表1～3所示。

2.1 骨架粒徑大小對動態(tài)彈性參數(shù)的影響

根據(jù)實驗，得到如圖2、3所示的關(guān)系。骨架粒徑分別為小于0.125 mm、介于0.125～0.600 mm和大于0.600 mm的水合物沉積物試件，隨著飽和度增加（5%～70%）其泊松比有一定呈線性下降；楊氏模量隨水合物飽和度增大，因為孔隙中水合物飽和度增大，孔隙中水的飽和度減小，所以強度增加。楊氏模量隨著粒徑增大而增大，因為在一定的骨架粒徑范圍內(nèi)，巖土體剪切面主要是切割其骨架部分，所以在一定粒徑范圍內(nèi)，水合物沉積物的強度隨著骨架粒徑增大而增大。

表1 不同骨架粒徑下的水合物沉積物的彈性力學(xué)參數(shù)實驗及計算數(shù)據(jù)表

表2 不同泥質(zhì)含量下的水合物沉積物的彈性力學(xué)參數(shù)實驗及計算數(shù)據(jù)表

表3 不同軸壓下的水合物沉積物的彈性力學(xué)參數(shù)實驗及計算數(shù)據(jù)表

圖2 不同骨架粒徑下泊松比與水合物飽和度的關(guān)系圖

圖3 不同骨架粒徑下楊氏模量與水合物飽和度的關(guān)系圖

2.2 泥質(zhì)含量對動態(tài)彈性參數(shù)的影響

水合物沉積物在相同水合物飽和度條件下，動態(tài)力學(xué)彈性參數(shù)與不同泥質(zhì)含量的關(guān)系如圖4、5所示。從圖4、5可看出，泥質(zhì)含量分別為5%、10%、20%的水合物沉積物樣品隨著泥質(zhì)含量的增加其動態(tài)泊松比有所增大；動態(tài)楊氏模量都隨著泥質(zhì)含量的增加呈冪函數(shù)減小。因為隨著泥質(zhì)含量的增加，泥質(zhì)填充在石英砂骨架中，占據(jù)了一定的孔隙空間，使這些地方無法生成水合物，又因為水合物的強度大于黏土的強度，所以總體上水合物沉積物的楊氏模量隨著泥質(zhì)含量增加而減小。

圖4 不同水合物飽和度下泊松比與泥質(zhì)含量的關(guān)系圖

圖5 不同水合物飽和度下楊氏模量與泥質(zhì)含量的關(guān)系圖

利用在相同飽和度下不同泥質(zhì)含量與泊松比、楊氏模量回歸分析可得到他們之間的關(guān)系分別為：

式中Ssh表示泥質(zhì)含量；a、b、c、d分別表示與水合物飽和度有關(guān)的系數(shù)。a、b、c、d取值如表4所示。

2.3 軸壓對動態(tài)彈性參數(shù)的影響

根據(jù)實驗在5 MPa的圍壓下，施加10 MPa、15 MPa、20 MPa的軸壓所得數(shù)據(jù)（圖6、7）表明，天然氣水合沉積物在不同水合物飽和度下，施加軸壓對泊松比的影響并不大；軸壓大小對泊松比影響不大，與Masui研究發(fā)現(xiàn)水合物沉積物的泊松比與水合物飽和度關(guān)系的結(jié)論一致；而楊氏模量在相同水合物飽和度下和軸壓呈一定對數(shù)關(guān)系，隨著軸壓增大而增大；水合物沉積物的強度主要由水合物強度、巖石顆粒間摩擦力和基質(zhì)顆粒間的相互作用共同決定，增加軸壓抑制了壓實過程中的裂隙發(fā)育，使其顆粒間的摩擦力及咬合力增加，使得水合物沉積物顆粒之間的膠結(jié)能力增強，從而增加了試樣的強度。

表4 a、b、c、d在不同水合物飽和度情況下的取值表

圖6 不同水合物飽和度下泊松比與軸壓的關(guān)系圖

圖7 不同水合物飽和度下楊氏模量比與軸壓的關(guān)系圖

利用相同水合物飽和度在不同軸壓情況下的動態(tài)楊氏模量回歸分析，可得到水合物沉積物在相同水合物飽和度下動態(tài)楊氏模量與軸壓的關(guān)系如下：

式中σ表示軸壓，MPa；e、f分別表示與水合物飽和度有關(guān)的系數(shù)。e、f取值如表5所示。

表5 e、f在不同水合物飽和度情況下的取值表

3 結(jié)論

本次研究采用SHW-Ⅲ型水合物聲電力學(xué)測試裝置，合成了不同飽和度下3種不同骨架粒徑、3種不同泥質(zhì)含量、3種不同三軸軸壓下的水合物沉積物試件，并測得了各個試件的水合物飽和度、巖電和聲波數(shù)據(jù)；再根據(jù)彈性波動理論計算出不同試件的泊松比和動態(tài)彈性模量，并進行了分析。數(shù)據(jù)分析顯示：在線性坐標(biāo)系下，一定水合物飽和度范圍內(nèi)水合物沉積物中的水合物飽和度和沉積物試件的彈性參數(shù)呈一定的正比關(guān)系，隨著水合物飽和度的增大，沉積物試件的泊松比減小，動態(tài)楊氏模量增大；水合物飽和度在5%～70%時，隨著試件的骨架粒徑增大（0.125～1.180 mm），沉積物試件的泊松比降低，動態(tài)楊氏模量增大；隨著泥質(zhì)含量的增大（5%～20%），沉積物試件的動態(tài)楊氏模量呈現(xiàn)冪函數(shù)減??；在5 MPa的圍壓條件下，隨著軸壓的增大水合物沉積物的動態(tài)楊氏模量大幅增大，而泊松比與軸壓的關(guān)系不大。根據(jù)實驗數(shù)據(jù)，建立動態(tài)泊松比與泥質(zhì)含量、動態(tài)楊氏模量與泥質(zhì)含量及軸壓的關(guān)系模型。

利用聲波測井方法分析水合物的動態(tài)彈性力學(xué)參數(shù)，可解決水合物取心難、造價高、無法獲取連續(xù)數(shù)據(jù)的問題。本次實驗所得到的結(jié)果符合客觀的巖石物理實驗現(xiàn)象及規(guī)律，但仍需要進一步研究水合物沉積物的動態(tài)彈性參數(shù)與靜態(tài)彈性參數(shù)之間的關(guān)系。