各向異性度對水平井井壁穩定性的影響

肖欽萍，賀云婷，何雅文

（1.西安石油大學機械工程學院，陜西西安 710065；2.西安石油大學石油工程學院，陜西西安 710065）

在鉆井過程中由于原位應力場，往往存在力學穩定性問題，因此需要進行井壁穩定力學研究[1]。在以往的工程設計中，往往粗略的將巖體近似為各向同性體，這對于淺部地層來說是合理的；但是隨著鉆井深度的增加，深度井井壁圍巖表現出強烈的各向異性特性從而導致井壁失穩，此時采用簡單的各向同性應力方程已經不能滿足安全鉆井的要求[2，3]。在各向異性地層井壁應力分布研究方面，Lekhnitskii、Amadei、Aadnoy 和Ong[4-7]建立了考慮各向異性介質的井壁應力計算模型；衡帥等[8]研究了巖石的力學特性、強度特性和破裂模式的各向異性，給出了橫向各向同性體的五個材料參數；王漢等[9]建立了各向異性巖石的井壁應力計算模型，分析了巖石力學參數各向異性對井壁應力分布的影響；王倩等[10]通過實驗和理論計算的方法驗證了各向異性度可評價巖石的各向異性程度。但是這些研究都沒有進一步分析各向異性度對井壁穩定性的影響，鑒于此，本文基于層理地層巖石力學特性，用Matlab軟件獲得了各向異性度對水平井壁應力分布的規律，然后通過考慮拉伸失效準則和剪切失效準則來研究各向異性度對水平井壁穩定性的影響。

1 各向異性度

各向異性度k定義為層理面的楊氏模量E與垂直于層理面的楊氏模量E'之比，即k＝E/E'，即k=1為各向同性。各向異性度越大，各向異性的程度越高，不同的材料其各向異性度也不一樣。

2 各向異性度對水平井壁應力分布的影響

當井眼軸線平行于最大主應力方向，作用在水平井壁上兩個相互垂直的σv和σh，min之間的值差異較大時，各向異性度對水平井壁應力分布的影響較為突出[11]，所以本文在建立的橫向各向異性巖石水平井壁應力分布模型的基礎上，只計算此種情況下水平井壁應力分布情況。取最大水平地應力σh，max=45 MPa，最小水平地應力 σh，min=35 MPa，上覆地應力 σv=70 MPa，αb=45°，βb=90°，αs=ατ=0°，βs=45°，βτ=30°，彈性模量 E=31.75 GPa，泊松比v=0.28，鉆井液液柱壓力pw=45 MPa。

由于進行井壁穩定性分析時，只用到最小主應力，所以這里只分析水平井壁上的情況。最小主應力分布用圓柱坐標系畫出（見圖1）。由于圖1關于中心對稱，所以分析圖形時只討論井周角θ在0°～90°這部分圓弧。從圖1可以觀察到：最小主應力受各向異性度的影響較小。只有當θ在0°～20°時，對于不同的各向異性度最小主應力的值才存在差異，各向異性度越大，其最小主應力的值越小，θ=0°時最小主應力的值最小；θ≥20°時，對于不同的各向異性度最小主應力的值近乎一樣。

考慮了不同各向異性度水平井壁處的切應力分布圖（見圖2）。切應力分布圖關于中心對稱，所以分析圖形時也只討論井周角θ為0°～90°這部分圓弧。從圖2可以觀察到：當θ角接近0°和 90°時，各向異性度越大其切應力的值越小；當θ=0°時切應力總是最小的，θ=90°時切應力總是最大的；每四分之一個圓弧上切應力曲線有兩個交點，第一個交點大約在θ=20°處相交，相交之后各向異性度越大其切應力的值越大；第二個交點大約在θ=60°處相交，相交之后各向異性度越大其切應力的值越小；當k≥4時，切應力曲線近乎一致。

3 水平井壁穩定性分析

分析了水平井壁處的主應力分布和切應力分布，還需要進一步分析水平井壁在該應力狀態下是否穩定。當作用在井壁上的應力狀態超過地層強度時，認為井壁開始失效，在進行水平井壁穩定性分析時通常只關注拉伸失效和剪切失效[12]。

圖1 井壁處的最小主應力分布

圖2 井壁處的切應力分布

拉伸失效準則選用最大主應力強度理論來分析水平井壁穩定性，只要一個主應力達到巖石的抗壓強度或者巖石的抗拉強度，巖石便失效。井壁的拉伸破壞通常是由最小主應力來定義的，所以本文用最小主應力與巖石的抗拉強度來做比較。而常用的剪切失效準則有Mohr-Coulomb準則、Tresca's準則、Lade準則和Drücker-Prager準則，Mohr-Coulomb 準則是使用時間最長和最常用的準則，本文選用Mohr-Coulomb準則分析水平井壁穩定性。

3.1 拉伸失效準則

當井壁處的最小主應力大于或等于巖石的拉伸強度時，巖石開始失效，即當σ3≥-σt時，發生失效。

當θ=0°時最小主應力達到最小值，所以只分析θ=0°時不同各向異性度的最小主應力曲線。不同各向異性度的拉伸失效準則曲線（見圖3），水平黑線在縱軸上的值為巖石的拉伸強度，如果最小主應力在黑線對應值的下方，表明水平井壁拉伸失效。從圖3可以觀察到：當k≤2時，對于任何方位角α，水平井壁都是穩定的；當k=3時，水平井壁在α為0°～20°時發生拉伸失效，α為20°～90°時穩定；當k=4時，水平井壁在α為0°～30°時發生拉伸失效，α 為 30°～90°時穩定；當 k=5 時，水平井壁在 α 為 0°～40°時發生拉伸失效，α 為 40°～90°時穩定；當k≥3時，隨著各向異性度的增加，水平井壁發生拉伸失效的角度在逐漸增大。

3.2 剪切失效準則

剪切失效準則為[13]：

其中Φ為內摩擦角，當F≤0時發生失效。

當θ=90°時切應力達到最大值，所以只分析θ=90°時不同各向異性度的切應力曲線，不同各向異性度的Mohr-Coulomb失效準則曲線（見圖4），根據Mohr-Coulomb失效準則，F≤0時發生剪切失效。從圖4可以觀察到：對于k=1的情況，大約在α≤35°時水平井壁發生剪切失效，α≥35°時水平井壁穩定，所以用k=1假設獲得的井壁穩定性預測僅適用于α≥35°的情況；當k≥2時井壁不發生剪切失效，并且隨著各向異性度的增加，曲線遠離失效線（F=0）。

圖3 不同各向異性度的拉伸失效準則曲線

圖4 巖石不同各向異性度的Mohr-Coulomb失效準則

4 結論

（1）最小主應力在 θ為 0°～20°時，各向異性度越小，其應力分布越接近k=1的應力分布，隨著各向異性度的增加，其應力分布與k=1的應力分布之間的差距增大；θ≥20°時不受各向異性度的影響，其應力分布與k=1的應力分布一致。所以在實驗中確定各向異性度時務必謹慎，以提高穩定性預測的準確性。

（2）當各向異性度 k≤2，α≥35°時水平井壁既不發生拉伸失效也不發生剪切失效，這時巖石可以做各向同性假設；隨著各向異性度k的增加，井壁的拉伸失效趨勢增加，而剪切失效趨勢減小。

參考文獻：

[1]陳勉，金衍，張廣清.石油工程巖石力學[M].北京：科學出版社，2008.

[2]Lajtai E.Z.Strength of discontinuous rocks in shear[J].Geotechnique，1969，19（2）：218－233.

[3]陳新，楊強，何滿潮，等.考慮深部巖體各向異性強度的井壁穩定性分析[J].巖石力學與工程學報，2005，24（16）：2882－2888.

[4]Lekhnitskii，S.G.Theory of Elasticity of an Anisotropic Body[M].Trranslated by Fem P.San Francisco，Holden－Day lne，1963.

[5]Amadei B.In situ stress measurements in anisotropic rock[J].International Journal of Rock Mechanics and Mining Science&Geomechanics Abstracts，1984，21（6）：327－338.

[6]Aadnoy B S，Chenevert M E.Stability of highly inclined boreholes[J].SPE Drilling Engineering，1987，2（4）：364－374.

[7]Ong S H.Borehole stability[M].Norman：University of Oklahoma，1994：201－210.

[8]衡帥，楊春和，張保平，等.頁巖各向異性特征的試驗研究[J].巖土力學，2015，36（3）：609－616.

[9]王漢，陳平，張智.頁巖力學參數各向異性對井壁應力分布的影響[J].長江大學學報（自然版），2015，12（20）：41－46.

[10]王倩，王鵬，等.頁巖力學參數各向異性研究[J].天然氣工業，2012，32（12）：62－65.

[11]BHARGAVA R，SHARMA K.A study of finite size effects on cracked 2-D piezoelectric media using extended finite elementmethod [J].ComputationalMaterialsScience，2011，50（6）：1834-1845.

[12]Mclean，M.R.，&Addis，M.A.Wellbore stability：The effect of strength criteria on mud weight recommendations[J].SPE Annual Technical Conference and Exhibition，1990：23－26.

[13]Yi，X.，Russell，J.E.，&Ong，S.H.Improving Borehole Stability Analysis by Quantifying the Effects of Intermediate Principal Stress using Polyaxial Rock Strength Test Data[M].The 40th U.S.Symposium on Rock Mechanics（USRMS），Anchorage，AK，2005.