頁(yè)巖水平井井眼軌跡方位與層理面產(chǎn)狀的關(guān)系

夏宏泉 ，史亞紅 ，王瀚瑋 ，趙昊 ，朱偉厚

（1.西南石油大學(xué)油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610500；2.中國(guó)石油西南油氣田公司蜀南氣礦，四川 瀘州 646000；3.中國(guó)石化中原油田分公司工程技術(shù)管理部，河南 濮陽(yáng) 457001）

0 引言

世界頁(yè)巖氣資源量約為456×1012m3，相當(dāng)于煤層氣（256×1012m3）與致密砂巖氣（209×1012m3）的總和。我國(guó)頁(yè)巖氣資源量約為26×1012m3，占全球頁(yè)巖氣資源量的5.7%［1-3］。水平井是頁(yè)巖氣開(kāi)發(fā)的核心技術(shù)之一，已廣泛應(yīng)用于頁(yè)巖氣開(kāi)發(fā)。頁(yè)巖儲(chǔ)層具有層理和天然裂隙，易膨脹破碎，水平段鉆井易發(fā)生井漏、垮塌等。針對(duì)頁(yè)巖儲(chǔ)層特征，開(kāi)展頁(yè)巖氣水平井井壁穩(wěn)定機(jī)理研究，避免工程事故發(fā)生，對(duì)現(xiàn)場(chǎng)具有重要指導(dǎo)意義［4］。國(guó)內(nèi)外學(xué)者在研究井壁穩(wěn)定問(wèn)題方面，已經(jīng)取得了顯著成就。例如國(guó)外的Okland等［5-6］通過(guò)實(shí)驗(yàn)得出頁(yè)巖的井壁穩(wěn)定與層理面的強(qiáng)度各向異性有很大關(guān)系，國(guó)內(nèi)的金衍等［7-8］采用單一弱面強(qiáng)度準(zhǔn)則研究了弱面地層直井與斜井井壁的坍塌壓力計(jì)算模型，劉向君等［9］研究了巖石弱面傾角和傾斜方位角改變對(duì)直井、斜井井壁穩(wěn)定性的影響，但是對(duì)于頁(yè)巖水平井井壁穩(wěn)定性問(wèn)題的研究不夠完善，難以保證頁(yè)巖氣的經(jīng)濟(jì)高效開(kāi)發(fā)。

本文以長(zhǎng)寧頁(yè)巖氣示范區(qū)為研究對(duì)象，基于頁(yè)巖巖石強(qiáng)度單一弱面剪切破壞準(zhǔn)則和頁(yè)巖水平井井壁應(yīng)力分布模型，建立了計(jì)算不同層理面產(chǎn)狀和不同井眼方位的頁(yè)巖水平井坍塌壓力模型。在此基礎(chǔ)上，定量分析層理面產(chǎn)狀對(duì)水平井井壁坍塌壓力的影響，旨在尋找頁(yè)巖層理面產(chǎn)狀與水平井井眼軌跡方位的關(guān)系規(guī)律，用以指導(dǎo)工區(qū)頁(yè)巖氣水平井的鉆井設(shè)計(jì)和施工。

1 井壁剪切破壞準(zhǔn)則和應(yīng)力分布模型

1.1 頁(yè)巖巖石強(qiáng)度單一弱面剪切破壞準(zhǔn)則

頁(yè)巖強(qiáng)度在很大程度上受層理弱面的影響［10］。國(guó)外對(duì)于層理面巖石的研究已經(jīng)比較成熟，為了更好地描述巖石的各向異性，許多學(xué)者提出了各向異性破壞準(zhǔn) 則［11-13］。尤 其 是 J.C.Jaeger［14］對(duì) 摩 爾-庫(kù) 倫（Mohr-Coulomb）破壞準(zhǔn)則進(jìn)行推廣而提出了單一弱面準(zhǔn)則，對(duì)研究巖石的各向異性產(chǎn)生了重大影響。結(jié)合該準(zhǔn)則，經(jīng)過(guò)推導(dǎo)計(jì)算，得到頁(yè)巖巖石強(qiáng)度單一弱面剪切破壞準(zhǔn)則，如圖1所示。

圖1 含單一弱面頁(yè)巖巖石剪切破壞理論分析示意

巖石本體的剪切破壞模型為

巖石層理面的剪切破壞模型為

引入巖石沿層理破壞面法向與井壁最大主應(yīng)力的夾角φw，對(duì)式（2）變形整理，可得另一種形式的巖石層理面剪切破壞模型：

其中

式中：σ1，σ3分別為最大、最小主應(yīng)力，MPa；c0，cw分別為本體、 層理面的黏聚力，MPa；φ0，φw分別為本體、層理面的內(nèi)摩擦角，（°）。

巖石沿層理面產(chǎn)生滑移破壞的條件為φ1＜φw＜φ2，此時(shí)用式（3）計(jì)算頁(yè)巖的坍塌壓力；頁(yè)巖本體發(fā)生破壞的條件為 0＜φw＜φ1和 φ2＜φw＜90°，此時(shí)用式（1）計(jì)算頁(yè)巖的坍塌壓力。 其中，φ1和 φ2的值可由式（4）和式（5）確定［15］。

1.2 頁(yè)巖水平井井壁應(yīng)力分布模型

水平井井壁穩(wěn)定性問(wèn)題為定向井井壁穩(wěn)定性問(wèn)題的特例［16］。研究頁(yè)巖定向井壁穩(wěn)定性問(wèn)題的關(guān)鍵是建立井壁應(yīng)力分布模型。令σv為上覆地應(yīng)力，σH和σh分別為水平最大主應(yīng)力和水平最小主應(yīng)力。為方便起見(jiàn)，建立地應(yīng)力坐標(biāo)系（1，2，3）和井眼坐標(biāo)系（x，y，z），選取坐標(biāo)系（1，2，3）分別與主地應(yīng)力 σH，σh，σv方向一致。oz軸對(duì)應(yīng)于井軸，ox和oy位于與井軸垂直的平面之中。其中：定義ψ為井斜角（z軸與地層上覆層壓力σv方向的夾角），β為井筒在水平面內(nèi)的投影與最大水平主應(yīng)力σH的夾角。任意的井斜角ψ和β角就可以確定出井筒的方位。如圖2所示，可以建立井眼坐標(biāo)系（x，y，z）和地應(yīng)力坐標(biāo)系（1，2，3）之間的轉(zhuǎn)化關(guān)系［17］。令ψ=90°，得到頁(yè)巖水平井井壁圍巖應(yīng)力分布模型［18］：

圖2 地應(yīng)力坐標(biāo)系和井眼坐標(biāo)系間的轉(zhuǎn)換關(guān)系

式中：σr，σθ，σz，τθz分別為與水平井對(duì)應(yīng)的徑向應(yīng)力、周向應(yīng)力、軸向應(yīng)力和剪應(yīng)力，MPa；pm為鉆井液柱壓力，MPa；pp為地層孔隙壓力，MPa；θ為井周角，（°）；ν為泊松比；α為有效應(yīng)力系數(shù)；φ為地層孔隙度；δ為滲透系數(shù)（當(dāng)井壁不可滲透時(shí)，δ=0；當(dāng)井壁滲透時(shí)，δ=1）；ζ為滲流效應(yīng)系數(shù)。

根據(jù)應(yīng)力狀態(tài)分析，井壁處3個(gè)主應(yīng)力分別為

將式（6）代入式（7）中，并將其從大到小排序，可得到三向應(yīng)力狀態(tài)下的井壁主應(yīng)力分別為σ1，σ2，σ3，即求得頁(yè)巖水平井井壁最大、最小主應(yīng)力計(jì)算模型。

2 水平井井壁坍塌壓力求解

在大地坐標(biāo)系中，斜井的井斜角為ψ、方位角為β1，弱面走向?yàn)門(mén)R，弱面傾角為DIP，則弱面法線方向矢量n為

井壁最大主應(yīng)力σ1的方向矢量N在斜井直角坐標(biāo)系中可表示為

斜井直角坐標(biāo)系中井壁最大主應(yīng)力σ1的方向矢量N在大地坐標(biāo)系中可表示為

其中：b1=cos β1cosψsinθ-sinβ1cosθ+cos β1sinψcos γ；b2=sinβ1cosψsinθ+cosβ1cosθ+sinβ1sinψcosγ；b3=-sinψ·sinθ+cosψcosγ。

而井壁最大主應(yīng)力σ1的作用面與井軸oz間夾角γ 的計(jì)算公式為［19］

最后得到井壁最大主應(yīng)力與弱面法向的夾角為

井壁最大、最小主應(yīng)力σ1，σ3都是鉆井液柱壓力pm的非線性函數(shù)，將σ1，σ3按條件代入巖石剪切破壞模型式（1）或式（3），則可求出含弱面的地層水平井井壁坍塌壓力的計(jì)算模型（方程）。通過(guò)解此非線性方程可求出維持弱面穩(wěn)定的鉆井液柱壓力安全下限值，此值就是所求的坍塌壓力值pb。利用軟件進(jìn)行數(shù)值模擬，可以得到不同傾角和傾向含層理面的頁(yè)巖井壁失穩(wěn)風(fēng)險(xiǎn)分布的演變過(guò)程。這為進(jìn)一步研究頁(yè)巖水平井段坍塌風(fēng)險(xiǎn)隨層理面產(chǎn)狀變化的規(guī)律提供了依據(jù)。

3 層理面產(chǎn)狀對(duì)水平段方位選取的影響

以H1井為研究對(duì)象，根據(jù)建立的頁(yè)巖水平井井壁剪切破壞模型和井壁應(yīng)力分布模型，編制相應(yīng)的程序。利用測(cè)井資料計(jì)算其2 750～3 050 m的地應(yīng)力和三壓力剖面，由剖面圖可知其井深2 800 m巖石力學(xué)和地應(yīng)力數(shù)據(jù)（見(jiàn)表1）。

表1 H1井2 800 m巖石力學(xué)參數(shù)和地應(yīng)力統(tǒng)計(jì)

利用表1數(shù)據(jù)，模擬層理面傾角和傾向?qū)?yè)巖水平井井眼軌跡方位的影響程度，得到了層理面產(chǎn)狀對(duì)頁(yè)巖井壁坍塌失穩(wěn)的影響規(guī)律，如圖3和圖4所示。圖3 為頁(yè)巖層理面傾向?yàn)?30°， 傾角分別為 0°，15°，45°，75°時(shí)頁(yè)巖水平井井壁坍塌壓力當(dāng)量密度的變化規(guī)律。圖4為頁(yè)巖層理面傾角為60°，傾向分別為0°，45°，90°，135°，180°，225°時(shí)頁(yè)巖水平井井壁坍塌壓力當(dāng)量密度的變化規(guī)律（從藍(lán)色到紫色坍塌壓力依次升高）。

由此，可以得出坍塌壓力隨層理面與井眼軌跡方位的變化關(guān)系規(guī)律：

1）由圖3可知，當(dāng)層理面傾角小于45°時(shí)，井周360°應(yīng)力云圖呈不對(duì)稱(chēng)分布，紅色區(qū)域明顯大于藍(lán)色區(qū)域，井壁巖石坍塌掉塊嚴(yán)重，井壁巖石發(fā)生層理面的剪切滑移；當(dāng)層理面傾角大于45°時(shí)，井周360°應(yīng)力云圖基本呈對(duì)稱(chēng)分布，藍(lán)色區(qū)域逐漸開(kāi)始變大，井壁巖石變得越加穩(wěn)定，井壁巖石發(fā)生本體破壞。

圖3 傾向30°、傾角不同時(shí)頁(yè)巖井壁坍塌失穩(wěn)規(guī)律

圖4 傾角60°、傾向不同時(shí)頁(yè)巖井壁坍塌失穩(wěn)規(guī)律

隨層理面傾角的逐漸增大，頁(yè)巖水平井井壁坍塌壓力當(dāng)量密度范圍逐漸變寬。當(dāng)井壁巖石發(fā)生層理面破壞時(shí)（DIP＜45°），坍塌壓力當(dāng)量密度在 1.00～1.37 g/cm3變化，變化范圍較窄；當(dāng)井壁巖石發(fā)生本體破壞時(shí)（DIP＞45°）， 坍塌壓力當(dāng)量密度變化范圍為0.80～1.45 g/cm3，變化范圍變寬。

2）井壁巖石發(fā)生本體破壞時(shí)，在特定的方位鉆進(jìn)時(shí)井壁坍塌壓力明顯低于其他方位（見(jiàn)圖4，以DIP=60°為例分析）。例如，當(dāng)傾向?yàn)?5°時(shí)，沿傾向相反的方向225°鉆進(jìn)較安全；當(dāng)傾向?yàn)?35°時(shí)，沿傾向相反的方向315°鉆進(jìn)較安全。因此，當(dāng)層理面傾角大于45°時(shí)，頁(yè)巖水平井井壁巖石發(fā)生本體破壞，此時(shí)沿層理面傾向的相反方向鉆進(jìn)，可降低井壁坍塌風(fēng)險(xiǎn)。

4 應(yīng)用實(shí)例分析

同理，以H2井龍馬溪組為研究對(duì)象，其地層弱面φw為 25°、cw為 6 MPa，本體 φ0為 38°、c0為 18 MPa。 該區(qū)地層最小水平主應(yīng)力方位為NE25°，最大水平主應(yīng)力方位為NE115°，層理面傾向?yàn)镾E25°，層理面傾角為60°，且工區(qū)地應(yīng)力狀態(tài)為正常地應(yīng)力類(lèi)型。利用測(cè)井資料并結(jié)合測(cè)井精細(xì)解釋?zhuān)傻玫? 750～4 150 m地應(yīng)力和三壓力剖面。由剖面圖可知其井深為3 850 m和3 950 m的巖石力學(xué)和地應(yīng)力數(shù)據(jù)（見(jiàn)表2）。

表2 H2井龍馬溪組巖石力學(xué)參數(shù)和地應(yīng)力統(tǒng)計(jì)

以表2數(shù)據(jù)為依據(jù)，結(jié)合頁(yè)巖水平井井壁坍塌壓力求解公式，使用軟件模擬井深為3 850 m和3 950 m，井眼軌跡分別沿最小水平主應(yīng)力方位（NE25°）、最大主應(yīng)力方位（NE115°）和層理面相反方位（NW25°）的安全鉆井液密度窗口，模擬結(jié)果如圖5所示。由圖5可知，井深為3 850 m時(shí)，沿層理面傾向相反方向（NW25°）鉆進(jìn)時(shí)的鉆井液安全密度窗口（1.0～2.5 g/cm3），明顯寬于沿最小水平主應(yīng)力和最大水平主應(yīng)力鉆進(jìn)時(shí)鉆井液安全密度窗口（2.0～2.5 g/cm3和 1.4～2.5 g/cm3）。井深為3 950 m時(shí)，規(guī)律相似。因此，通過(guò)分析，將H2井水平段方位角選為NW25°有利于井壁穩(wěn)定。H2井水平段的實(shí)際鉆進(jìn)方位為NE8°，此方位鉆進(jìn)井壁坍塌掉塊嚴(yán)重，建議后期根據(jù)本文的方法優(yōu)化鉆井方位為NW25°。

H2井水平段頁(yè)巖地層實(shí)際鉆井過(guò)程中，鉆至3 750～3 850 m時(shí)，未考慮層理面的影響，采用密度為1.49 g/cm3的鉆井液鉆進(jìn)，井壁發(fā)生了明顯垮塌。為了阻止井壁繼續(xù)坍塌，考慮到層理面對(duì)頁(yè)巖水平井鉆井施工的影響，在3 850 m將鉆井液密度提高到1.99 g/cm3后，井壁垮塌掉塊有明顯減輕。因此，相比實(shí)際鉆井液密度，采用優(yōu)化后的地層鉆井液當(dāng)量密度，更有利于安全鉆井。

圖5 井深3 850 m和3 950 m不同方位坍塌壓力對(duì)比

5 結(jié)論

1）基于頁(yè)巖巖石強(qiáng)度單一弱面剪切破壞準(zhǔn)則和頁(yè)巖水平井井壁應(yīng)力分布模型，綜合考慮巖石本體強(qiáng)度、頁(yè)巖層理面產(chǎn)狀等因素，建立層理性頁(yè)巖水平井井壁地層的坍塌壓力計(jì)算模型。該模型可為頁(yè)巖水平井井眼軌跡方位和鉆井液密度的優(yōu)選提供參考依據(jù)。

2）層理面產(chǎn)狀對(duì)頁(yè)巖水平井井壁穩(wěn)定有顯著影響。當(dāng)頁(yè)巖層理面傾向一定時(shí)，隨層理面傾角的逐漸增大，頁(yè)巖水平井井壁坍塌壓力當(dāng)量密度范圍逐漸變寬。當(dāng)層理面傾角小于45°時(shí)，井壁巖石發(fā)生層理面的剪切破壞；當(dāng)層理面傾角大于45°時(shí)，井壁巖石發(fā)生本體破壞，且當(dāng)井壁巖石發(fā)生本體破壞時(shí)，頁(yè)巖水平井沿層理面傾向的相反方向鉆進(jìn)有利于井壁穩(wěn)定。

3）將該方法應(yīng)用于長(zhǎng)寧龍馬溪頁(yè)巖地層，發(fā)現(xiàn)其預(yù)測(cè)的理想鉆井液密度和實(shí)際維持井壁穩(wěn)定的鉆井液密度基本一致。新方法克服了傳統(tǒng)未考慮層理面產(chǎn)狀對(duì)井壁穩(wěn)定的影響，可為頁(yè)巖水平井高效、安全鉆進(jìn)與壓裂改造提供可靠力學(xué)參考。

［1］鄒才能，董大忠，王社教，等.中國(guó)頁(yè)巖氣形成機(jī)理，地質(zhì)特征與資源前景［J］.石油勘探與開(kāi)發(fā)，2010，37（6）：641-653.

［2］閆存章，黃玉珍，葛春梅，等.頁(yè)巖氣是潛力巨大的非常規(guī)天然氣資源［J］.天然氣工業(yè)，2009，29（5）：1-6.

［3］王中華.國(guó)內(nèi)頁(yè)巖氣開(kāi)采技術(shù)進(jìn)展［J］.中外能源，2013，18（2）：23-32.

［4］陳安明，張輝，宋占偉.頁(yè)巖氣水平井鉆完井關(guān)鍵技術(shù)分析［J］.石油天然氣學(xué)報(bào)，2012，34（11）：98-103.

［5］OKLAND D，COOK J M.Bedding-related borehole instability in highangle wells［R］.SPE 47285，1998.

［6］FILDANI A，HANSON A D，CHEN Z，et al.Geochemical characteristics of oil and source rocks and implications for petroleum systems，Talara Basin，northwest Peru［J］.AAPG Bulletin，2005，89（11）：1519-1545.

［7］金衍，陳勉，陳治喜，等.弱面地層的直井井壁穩(wěn)定力學(xué)模型［J］.鉆采工藝，1999，22（3）：13-14.

［8］金衍，陳勉，柳貢慧，等.弱面地層斜井井壁穩(wěn)定性分析［J］.石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1999，23（4）：33-35.

［9］劉向君，陳一健，肖勇.巖石軟弱面產(chǎn)狀對(duì)井壁穩(wěn)定性的影響［J］.西南石油學(xué)院學(xué)報(bào)（自然科學(xué)版），2001，23（6）：12-13.

［10］閆傳梁，鄧金根，蔚寶華，等.頁(yè)巖氣儲(chǔ)層井壁坍塌壓力研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2013，32（8）：1595-1602.

［11］MCLAMORE R，GRAY K E.The mechanical behavior of anisotropic sedimentaryrocks［J］.Journalof EngineeringforIndustry，1967，89（1）：62-73.

［12］HILL R.The mathematical theory of plasticity［M］.New York：Oxford University Press，1998：122-187.

［13］MCLAMORE R，GRAY K E.The mechanical behavior of anisotropic sedimentaryrocks［J］.JournalofEngineeringfor Industry，1967，89（1）：62-73.

［14］JAEGER J C.Shear failure of anistropic rocks［J］.Geological Magazine，1960，97（1）：65-72.

［15］蔡美峰，何滿(mǎn)潮，劉東燕.巖石力學(xué)與工程［M］.北京：科學(xué)出版社，2009：104-108.

［16］趙凱，袁俊亮，鄧金根，等.層理產(chǎn)狀對(duì)頁(yè)巖氣水平井井壁穩(wěn)定性的影響［J］.科學(xué)技術(shù)與工程，2013，13（3）：580-583.

［17］金衍，陳勉，柳貢慧，等.弱面地層斜井井壁穩(wěn)定性分析［J］.石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1999，23（4）：33-35.

［18］金衍，齊自立，陳勉，等.水平井試油過(guò)程裂縫性?xún)?chǔ)層失穩(wěn)機(jī)理［J］.石油學(xué)報(bào)，2011，32（2）：295-298.

［19］陳勉，金衍，張廣清.石油工程巖石力學(xué)［M］.北京：科學(xué)出版社，2008：99-118.