王倩 王鵬 項(xiàng)德貴 馮宇思

中國石油集團(tuán)鉆井工程技術(shù)研究院

頁巖氣儲(chǔ)層中的有機(jī)頁巖在形成過程中具有層理、片理等特征，組成的礦物結(jié)晶顆粒具有不同大小以及不同的組合方式，造成頁巖中具有不同層次的結(jié)構(gòu)構(gòu)造和這些結(jié)構(gòu)構(gòu)造的定向排列，所以有機(jī)頁巖具有明顯的強(qiáng)度和彈性各向異性。如果將頁巖儲(chǔ)層巖石當(dāng)做各向同性體來處理，是顯然不符合實(shí)際情況的。對于存在明顯層理的頁巖，在研究其力學(xué)性質(zhì)時(shí)要考慮其各向異性［1－4］。層理性頁巖力學(xué)性質(zhì)受巖石本體和結(jié)構(gòu)面的力學(xué)性質(zhì)共同控制，通過理論和實(shí)驗(yàn)研究獲得頁巖巖石強(qiáng)度、彈性力學(xué)參數(shù)的各向異性特征，結(jié)合聲波各向異性分析結(jié)果，進(jìn)而完成有機(jī)頁巖力學(xué)性質(zhì)評價(jià)，能實(shí)現(xiàn)頁巖氣水平井的安全鉆井和優(yōu)化水力壓裂過程。

1 層理頁巖彈性參數(shù)變化規(guī)律的理論模型

層理頁巖的各向異性變動(dòng)如圖1所示，垂直構(gòu)造面的壓縮性顯著大于平行構(gòu)造面的壓縮性，所以，這種巖石在層理面的應(yīng)力—應(yīng)變性狀是各向同性的。

圖1 層理頁巖的各向異性變形圖

描述這種性狀的巖石可以用5個(gè)獨(dú)立的彈性常數(shù)。分別用E1和E2表示平行和垂直各向同性面的彈性模量，用ν1和ν2表示平行和垂直各向同性面的泊松比，用G2表示各向同性面內(nèi)的剪切模量。假定在局部坐標(biāo)系（x′，y′，z′）中，z′軸和可壓縮性最大面的法線方向一致，而x′和y′軸落在各向同性面內(nèi)。橫觀各向同性巖石的應(yīng)力—應(yīng)變關(guān)系為：

式中

在實(shí)際工程中，由于地層存在一定的走向和傾角，通常和上述坐標(biāo)系不一致，如圖2所示。需要選用總體坐標(biāo)系（x，y，z）?？傮w坐標(biāo)系（x，y，z）和各向異性定向坐標(biāo)系（x′，y′，z′）的關(guān)系可通過走向角（α）和傾向角（β）建立。

圖2 總體坐標(biāo)系與局部坐標(biāo)系的關(guān)系圖

借助轉(zhuǎn)換矩陣T和T＊，把總體坐標(biāo)系中的應(yīng)力｛σ｝和應(yīng)變｛ε｝分別換算為局部坐標(biāo)系中的應(yīng)力｛σ′｝和應(yīng)變｛ε′｝。

矩陣T的表達(dá)式為：

式中

矩陣T＊和T有如下關(guān)系：

將式（2）和式（3）代入式（1），得：

上式兩邊同乘T＊－1，利用式（5）的關(guān)系，得出：

根據(jù)單軸壓縮邊界條件，若沿z軸方向壓縮，即應(yīng)力σx＝σy＝τxy＝τyz＝τzx＝0，α＝0，而β≠0，由式（7）得到：

由于在β＝45°情況下，ν2E1＝ν1E2，可得與z軸為任意角度β方向的彈性模量Eβ和泊松比νβ：

通過實(shí)驗(yàn)測定層理頁巖垂直層理和平行層理方向的5個(gè)彈性常數(shù)，便可由式（11）和（12）從理論上得出層理頁巖彈性模量E和泊松比v隨層理傾角β的變化關(guān)系。

2 層理頁巖彈性常數(shù)各向異性的實(shí)驗(yàn)研究

為探討層理頁巖的各向異性，按一定的β方向取樣（β角為層理面與巖樣端面間的夾角）。β角從0°～90°，每隔15°取1個(gè)樣，共計(jì)7個(gè)樣（圖3）。巖樣上下端面的不平行度小于0．2mm。對每塊巖樣做單軸壓縮實(shí)驗(yàn)，測定巖石抗壓強(qiáng)度、彈性模量、泊松比，并且在實(shí)驗(yàn)中測量每塊巖樣的彈性波速度Vp、Vs［5－7］。

圖3 巖心取樣示意圖

由聲波實(shí)驗(yàn)結(jié)果根據(jù)式（13）、（14）計(jì)算各巖樣的動(dòng)態(tài)彈性模量和動(dòng)態(tài)泊松比，計(jì)算結(jié)果列于表1。

表1 巖石彈性模量的實(shí)測值Es和計(jì)算值Eβ表

式中Ed為巖石動(dòng)態(tài)彈性模量，MPa；νd為巖石動(dòng)態(tài)泊松比；Vp為縱波速度，m／s；Vs為橫波速度，m／s；ρ為巖石密度，g／cm3。

根據(jù)單軸壓縮實(shí)驗(yàn)所得的巖樣應(yīng)力—應(yīng)變曲線，可求得巖石的靜態(tài)彈性模量Es和靜態(tài)泊松比νs，其值列于表2。

由表2測定的各個(gè)角度巖樣的單軸抗壓強(qiáng)度可以看出，β為60°時(shí)巖石的抗壓強(qiáng)度最小，即此時(shí)巖樣沿層理方向破壞。根據(jù)層理傾角為βmin＝45°＋φs／2時(shí)發(fā)生層理破壞，可以得到層理面的內(nèi)摩擦角φs為30°。還可以看出，如果β小于內(nèi)摩擦角，巖樣不會(huì)發(fā)生沿層理剪切破壞。

表2 巖石力學(xué)參數(shù)實(shí)驗(yàn)和聲波實(shí)驗(yàn)數(shù)據(jù)表

圖4為根據(jù)表2的實(shí)驗(yàn)結(jié)果得出的動(dòng)態(tài)彈性模量、泊松比和靜態(tài)彈性模量、泊松比隨層理傾角的變化圖。

由圖4可見Ed、Es、νd、νs都隨層理傾角（β）的增加呈上升趨勢，靜態(tài)參數(shù)比動(dòng)態(tài)參數(shù)隨層理傾角的變化要陡一些，動(dòng)態(tài)參數(shù)總比靜態(tài)參數(shù)高，動(dòng)態(tài)彈性模量為靜彈性模量的1．35～1．50倍，動(dòng)態(tài)泊松比為靜態(tài)泊松比的1．00～1．15倍。

對橫觀各向同性體來說，5個(gè)獨(dú)立的彈性常數(shù)（E1、E2、ν1、ν2、G2）確定之后，由式（11）可以計(jì)算出任一β角方向上的彈性模量值，表2為彈性模量實(shí)測值和計(jì)算值的結(jié)果數(shù)據(jù)。

從表2可以看出，計(jì)算的彈性模量與實(shí)測值接近，最大誤差為3%。因此，只要對巖石作平行和垂直層理以及45°方向的3組樣品單軸壓縮實(shí)驗(yàn)，就可用計(jì)算值來代替樣品的實(shí)測值，而不需要對每一個(gè)方向（β角）的樣品進(jìn)行單軸壓縮實(shí)驗(yàn)來確定模量，可省去大量樣品加工及實(shí)驗(yàn)工作。結(jié)果表明將層理巖石視為橫觀各向同性介質(zhì)來計(jì)算彈性模量是可行的。

圖4 巖樣動(dòng)、靜態(tài)彈性參數(shù)實(shí)驗(yàn)結(jié)果圖

由表1的數(shù)據(jù)通過式（12）計(jì)算泊松比隨角β的變化規(guī)律，泊松比的實(shí)測值和計(jì)算值結(jié)果數(shù)據(jù)見表3。

表3 巖石泊松比的實(shí)測值vs和計(jì)算值vβ表

從表3可以看出，泊松比的實(shí)測值與計(jì)算值也比較接近，誤差在合理的范圍內(nèi)，所以將層理巖石視為橫觀各向同性介質(zhì)來計(jì)算泊松比是可行的。

定義頁巖的各向異性指數(shù)e為：

式中A1為平行于各向同性面上的彈性模量或泊松比；A2為垂直于各向同性面上的彈性模量或泊松比。

利用實(shí)驗(yàn)數(shù)據(jù)，可以得到巖石動(dòng)態(tài)彈性模量和泊松比的各向異性程度分別為0．088和0．089，巖石靜態(tài)彈性模量和泊松比的各向異性程度分別為0．180和0．155?？梢钥闯鰩r石各向異性在靜態(tài)時(shí)更為突出，而動(dòng)態(tài)時(shí)巖石各向異性有所減弱。巖石各向異性指數(shù)可評價(jià)巖石的各向異性程度，各向異性指數(shù)越大，表明巖石的各向異性程度越高。

3 結(jié)論

通過理論分析和實(shí)驗(yàn)研究頁巖力學(xué)性質(zhì)的各向異性，主要得出以下結(jié)論：

1）層理頁巖平行于層理面的彈性模量和泊松比大于垂直于層理面的彈性模量和泊松比。

2）層理頁巖的彈性常數(shù)Ed、Es、vd、vs都隨層理與水平方向夾角的減小呈上升趨勢，動(dòng)態(tài)彈性參數(shù)總比靜態(tài)彈性參數(shù)高，巖石各向異性在靜態(tài)時(shí)更明顯。

3）根據(jù)不同層理傾角頁巖巖樣的單軸抗壓強(qiáng)度結(jié)果可以看出，層理傾角為60°時(shí)巖樣的抗壓強(qiáng)度最小，巖樣沿層理方向破壞，若巖樣層理傾角小于層理面內(nèi)摩擦角，巖樣不會(huì)發(fā)生沿層理剪切破壞。

4）不同層理傾角頁巖巖樣的彈性模量、泊松比實(shí)驗(yàn)結(jié)果和將層理頁巖作為橫觀各向同性介質(zhì)理論計(jì)算得出的彈性模量、泊松比基本一致，說明將層理頁巖視為橫觀各向同性介質(zhì)是合理的。

5）根據(jù)巖石的各向異性指數(shù)可評價(jià)頁巖的各向異性程度，各向異性指數(shù)越大，表明巖石的各向異性程度越高。

［1］尹宏錦．地層各向異性—巖石可鉆性研究［J］．石油鉆采工藝，1989，11（1）：15－22．

［2］趙平勞，姚增．層狀巖石動(dòng)靜態(tài)力學(xué)參數(shù)相關(guān)性的各向異性［J］．蘭州大學(xué)學(xué)報(bào)：自然科學(xué)版，1993，29（4）：225－229．

［3］劉斌，席道瑛，葛寧潔，等．不同圍壓下巖石中泊松比的各向異性［J］．地球物理學(xué)報(bào)，2002，45（6）：880－890．

［4］何沛田，黃志鵬．層狀巖石的強(qiáng)度和變形特性研究［J］．巖土力學(xué)，2003，24（增刊1）：1－5．

［5］楊圣奇，張學(xué)民，蘇承東．巖塊聲學(xué)特性的實(shí)驗(yàn)研究［J］．遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2003，22（6）：772－775．

［6］鄭淑芬．地層速度各向異性及其應(yīng)用［J］．測井技術(shù)，2000，24（3）：223－226．

［7］王文冰，季榮生，喬艷紅，等．波形特征與層理巖石各向異性關(guān)系分析［J］．潔凈煤技術(shù)，2010，16（6）：83－85．