石油地质与工程２０１９年５月ＰＥＴＲＯＬＥＵＭＧＥＯＬＯＧＹＡＮＤＥＮＧＩＮＥＥＲＩＮＧ第３３卷第３期文章编号：１６７３—８２１７（２０１９）０３—００８２—０５非均匀地应力下井壁Ｖｏｎｉｉｓｅｓ应力及稳定性分析李冬梅１，窦益华２王小增３，李明飞２，蒲红斌４（１．中国石化西北油田分公司石油工程技术研究院，新疆乌鲁木齐８３００１ｌ；２．西安石油大学机械工程学院，陕西西安７１００６５；３．嘉应学院电子信息工程学院，广东梅州５１４０１５；４．中国石化河南油田分公司采油一厂，河南南阳４７３１３２）摘要：为分析非均匀地应力作用下井壁ＶｏｎＭｉｓｅｓ（米塞斯）应力及其对井壁稳定性的影响，建立了非均匀地应力作用下井壁附近岩石的力学分析模型。根据弹性力学中应力与应力函数之间的关系，应用叠加原理，导出了非均匀地应力下井壁径向应力、环向应力和剪切应力计算公式，分析了非均匀地应力对井壁稳定性的影响，确定了以钻井液密度为自变量的井壁ＶｏｎＭｉｓｅｓ应力的二次拟合公式。采用拟合公式计算的井壁ＶｏｎＭｉｓｅｓ应力误差小于２．５％，利用该公式可以确定井壁岩石不会破裂的钻井液密度。关键词：非均匀地应力；Ｙｏｎ中图分类号：ＴＥ２６Ｍｉｓｅｓ应力；井壁稳定性；叠加原理；钻井液密度文献标识码：ＡＶｏｎＭｉｓｅｓｓｔｒｅｓｓａｎｄｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｗｅｌｌｂｏｒｅｗａｌｌｕｎｄｅｒｎｏｎ－ｕｎｉｆｏｒｍｉｎ－ｓｉｔｕｓｔｒｅｓｓＬＩＤｏｎｇｍｅｉｌ，ＤＯＵＹｉｈｕａ２，ＷＡＮＧＸｉａｏｚｅｎ９１，ＬＩＭｉｎｇｆｃｉ２，ＰＵＨｏｎｇｂｉｎ４ｆ１．ＰｅｔｒｏｌｅｕｍＥｎｇｉｎｅｅｒｉｎｇＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＮｏｒｔｈｗｅｓｔＯｉｌｆｉｅｌｄＣｏｍｐａｎｙ，ＳＩＮＯＰＥＣ，Ｕｒｕｍｕｑｉ，Ｘｉｎｊｉａｎｇ８３００１１；２．ＳｃｈｏｏｌｏｆＭｅｃｈａｎｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ’ａｎＰｅｔｒｏｌｅｕｍＵｎｉｖｅｒｓｉｔｙ，Ｘｉ’ａｎ，Ｓｈａｎｎｘｉ７１００６５，Ｃｈｉｎａ；３．ＣｏｌｌｅｇｅｏｆＥｌｅｃｔｒｏｎｉｃＩｎｆｏｒｍａｔｉｏｎＯｉｌＰｒｏｄｕｃｔｉｏｎＰｌａｎｔＥｎｇｉｎｅｅｒｉｎｇ，ＪｉａｙｉｎｇＵｎｉｖｅｒｓｉｔｙ，Ｍｅｉｚｈｏｕ，Ｇｕａｎｇｄｏｎｇ５１４０１５，Ｃｈｉｎａ；４．ＴａｈｅＰｒｏｊｅｃｔＤｅｐａｒｔｍｅｎｔ，Ｎｏ．１ｏｆＨｅｎａｎＯｉｌｆｉｅｌｄＣｏｍｐａｎｙ，Ｓ１ＮＯＰＥＣ，Ｎａｎｙａｎｇ，Ｈｅｎａｎ４７３１３２，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：ＩｎｏｒｄｅｒｔｏａｎａｌｙｚｅｔｈｅｉｎｆｌｕｅｎｃｅｏｆＶｏｎＭｉｓｅｓｓｉｔｕｓｔｒｅｓｓａｎｄｉｔｓｅｆｌｅｅｔｏｎｓｔｒｅｓｓｏｎｔｈｅｂｏｒｅｈｏｌｅｗａｌｌｕｎｄｅｒｎｏｎ－ｕｎｉｆｏｒｉＢｉｎ－ｎｅａｒｔｈｅｓｔａｂｉｌｉｔｙｏｆｔｈｅｂｏｒｅｈｏｌｅｗａｌｌ，ａｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｍｏｄｅｌｏｆｔｈｅｒｏｃｋｓｔｒｅｓｓｔｈｅｗｅｌｌｂｏｒｅｗａｌｌｕｎｄｅｒｎｏｎ—ｕｎｉｆｏｒｌｎｉｎ—ｓｉｔｕａｎｄｓｔｒｅｓｓｗａｓｅｓｔａｂｌｉｓｈｅｄ．Ａｃｃｏｒｄｉｎｇｔｏｔｈｅｒｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｓｔｒｅｓｓｏｎｆｕｎｅｔｉｏｎｉｎｅｌａｓｔｉｃｍｅｃｈａｎｉｃｓ．ｂａｓｅｄｓｔｒｅｓｓｔｈｅｓｕｐｅｒｐｏｓｉｔｉｏｎｐｒｉｎｃｉｐｌｅｓ，ｔｈｅｆｏｒｍｕｌａｓｏｆｒａｄｉａｌｓｔｒｅｓｓｏｎｓｔｒｅｓｓ，ｃｉｒｃｕｍｆｅｒｅｎｔｉａｌａｎｄｓｈｅａｒｓｔｒｅｓｓａｒｅｄｅｒｉｖｅｄ，ｔｈｅｉｎｆｌｕｅｎｃｅｏｆｎｏｎ—ｕｎｉｆｏｒｍｉｎ—ｓｉｔｕｓｔｒｅｓｓｏｎｗｅｌｌｂｏｒｅｗａＵｓｔａｂｉｌｉｔｖｉｓａｎａｌｙｚｅｄ．ｔｈｅｑｕａｄｒａｔｉｃｆｉｔｔｉｎｇｆｏｒｍｕｌａｏｆｔｈｅＶｏｎＭｉｓｅｓｄｅｎｓｉｔｖｏｆｄｒｉｌｌｉｎｇｆｌｕｉｄａｓｔｈｅｗｅｌｌｈｏｒｅｗａｌｌｗｉｔｈｔｈｅｓｔｒｅｓｓｔｈｅｉｎｄｅｐｅｎｄｅｎｔｖａｒｉａｂｌｅｉｓｄｅｔｅｒｍｉｎｅｄ．ＴｈｅｃａｎｅｒｒｏｒｏｆＶｏｎＭｉｓｅｓｃａｌｃｕｌａｔｅｄｂｙｔｈｅｆｉｔｔｉｎｇｆｏｒｍｕｌａｉｓｌｅｓｓｔｈａｎ２．５％．ａｎｄｔｈｉｓｆｏｒｍｕｌａｔｈｅｗｅｌｌｂｏｒｅｗａｌｌｒｏｃｋｗｉｔｈｏｕｔｆｒａｃｔｕｒｅｓ．Ｋｅｙｂｅｕｓｅｄｔｏｄｅｔｅｒｍｉｎｅｔｈｅｄｅｎｓｉｔｙｏｆｄｒｉｌｌｉｎｇｆｌｕｉｄｉｎｗｏｒｄｓ：ｎｏｎ－ｕｎｉｆｏｒｍｉｎ——ｓｉｔｕｓｔｒｅｓｓ；ＶｏｎＭｉｓｅｓｓｔｒｅｓｓ；ｗｅｌｌｂｏｒｅｓｔａｂｉｌｉｔｙ；ｓｕｐｅｒｐｏｓｉｔｉｏｎｐｒｉｎｃｉｐｌｅ；ｄｒｉｌｌｉｎｇｆｌｕｉｄｄｅｎｓｉｔｙ钻井过程中的井壁稳定性是事关钻进速度和钻井安全的重要问题，井眼附近岩石应力大小及其分布规律是影响井壁稳定性的主要因素之一，负压钻井和钻遇枯竭地层时，井壁更易失稳且导致钻井事故。因井壁失稳，国内油田在钻井过程中发生过不象ｐ１；姚直书等考虑井壁的空间轴对称，运用叠加原收稿日期：２０１９０１—２４作者简介：李冬梅，高级工程师，１９７３年生，１９９６年毕业于中国地质大学（武汉）油藏工程专业，现从事完井测试与试油评价技术及管理工作。基金项目：国家自然科学基金项目“页岩气水平井压裂与生产套管变形机理及其控制机制研究”（ＮＳＦＣ５１６７４１９９ｏ同程度的井漏、井壁垮塌和卡钻等复杂事故”。２１，为此，人们对井壁稳定性问题开展了针对性研究。周国庆等结合模拟试验，分析了华东地区井壁破坏现李冬梅等．非均匀地应力下井壁ＶｏｎＭｉｓｅｓ应力及稳定性分析理导出了井壁应力计算公式【４】；李军陋１、王建军等陋ｌ研究了地层各向异性和多孔弹性易坍塌地层井壁围岩的应力分布规律；针对泥页岩地层钻井易发生井壁失稳、影响钻井作业效率等问题，冉小丰ｒ”、光新军等【８１基于连续损伤力学理论，将塑性损伤演化及渗流相互耦合的概念引入摩尔一库伦破坏准则，建立了基于等效塑性应变的损伤演化模型，给出了泥岩的渗透性演化方程；何朋立等一１研究表明，温度作用也是导致井壁内应力增大，诱发井壁破裂的重要因素；欧焕农”叫利用内外边界条件，给出了井壁环向应力分布的解析解和数值解；郭莉霞等”ｕ研究表明，最大主应力与最小主应力之差越过莫尔圆包络线将使井壁发生剪切破坏，诱发平行于井眼轴线的纵向裂缝；王渭明等”２】导出了正交各向异性多层复合井壁的应力、位移计算公式。可以看出，现有研究多针对均匀地应力下井壁的稳定性，而非均匀地应力作用下井壁应力、井壁稳定性问题有待深入研究。因此，本文建立了非均匀地应力下井壁岩石力学分析模型，根据弹性力学中应力与应力函数之间的关系，应用叠加原理，推导了非均匀地应力下井壁径向应力、环向应力和剪切应力公式，分析了非均匀地应力对井壁稳定性的影响，确定了以钻井液密度为自变量的井壁岩石等效应力拟合公式，利用该公式可以计算井壁岩石不会破裂的钻井液密度。１非均匀地应力下井壁力学模型地层钻开后，没有下人套管的井眼内壁受到水平方向上的地层压力和井眼内液体压力。井壁及周围岩层受力如图１ａ所示，ｐｈ是水平方向最小地应力，阳是水平方向最大地应力，ｐｆ是井眼内壁静液柱压力，单位均为ＭＰａ。为确定非均匀地应力作用下井眼内壁应力分布，可以应用叠加原理将井壁载荷分解为只有内压作用、两向压缩和拉压共同作用三种情况，分别如图ｌｂ～ｄ所示。井眼钻开后，非均匀地应力作用于距离井壁中心无穷远处，井壁周围岩石为各向同性连续介质。２非均匀地应力下井壁等效应力分析２．１内压作用下井壁应力分析内压作用下井壁受力模型如图１ｂ所示。根据Ｌａｍｅ厚壁筒理论，力学模型中地层宽度远大于井眼图１非均匀地应力作用下井壁力学分析模型直径，相当于Ｌａｍｅ公式中外径趋近于无穷大。由Ｌａｍｅ公式可知，内压下半径ｒ处井壁径向和环向应力分别为：ｄ２ｄ２听１２一万仇，∞１２历Ｐｉ／．、ＬＪ，式中：ｄ为井眼直径，ｍ２．２均匀压缩地应力下井壁应力分析均匀压缩地应力下井壁受力模型如图２所示。水平方向地应力相同，均为（ＰＨ＋Ｐｈ）／２。采用极坐标法求解井壁附近应力，将模型周围直边界变换为圆边界。选取直径Ｄ远大于ｄ的圆周，根据圣维南原理，该圆周上Ａ点的应力与无孑ＬＨ，－ｊ的应力相同。直角坐标系下该点的应力分量为：ａｘ＝一＿ＰＨ＋Ｐｈ，ｏ－ｙ＝一＿ＰＨ＋Ｐｈ，ｈｙｍｍ０（２）由应力直角坐标向极坐标转换的公式：听＝叹ＣＯＳ２９＋ｑｓｉｎ２日＋２ｒｚｙｓｉｎ臼ＣＯＳ９％２ｑＣＯＳ２日＋呶ｓｉｎ２臼一２Ｔｘｖｓｉｎ日ＣＯＳｐＴｒ０＝（ｏ－ｙ—ａｘ）ｓｉｎ臼ＣＯＳ８＋Ｔｘｖ（ＣＯＳ２臼一ｓｉｎ２８）（３）可得极坐标下直径＿Ｄ处的应力：ａｒｍｍ一—ＰＨ＿＋Ｐｈ，∞２一—ＰＨ＿＋Ｐｈ，Ｔｒ日＝０（４）式中：日为极坐标中井眼任一点和原点连线与ｘ轴的夹角，ｒａｄ。井壁附近径向和环向应力分布相当于内压０和外压（ｐＨ＋Ｐｈ）／２作用下，内径为ｄ和外径为Ｄ的圆筒的径向和环向应力，由Ｌａｍｅ公式和ｄ／Ｄ趋于０可得：％：一半（１一筹）石油地质与工程２０１９年第３期ａ０２－－一ＴＰＨ＋Ｐｈ＼、（１＋匀Ｔｒ日２＝０（５）图２均匀压缩地应力下井壁受力模型２．３拉伸、压缩地应力下井壁应力分析拉伸、压缩地应力共同作用下井壁受力模型如图３所示。一个方向上的拉伸地应力为（ｐＨ—ｐｈ）／２，其垂直方向上的压缩地应力为（ｐＨ—Ｐｈ）／２。采用极坐标求解井壁附近应力，将周围直边界变换为圆边界。选取直径Ｄ远大于ｄ的圆周，根据圣维南原理，该圆周上Ａ点的应力与无孑Ｌ时的应力相同，即：吼＝下Ｐｈ－－ＰＨ，ｑ＝一下Ｐｈ－－ＰＨ，ｋｙ＝０（６）由应力的直角坐标向极坐标转换公式可得直径Ｄ处的应力分布：听＝ＴＰｈ－－ＰＨＣＯＳ２０，ａｏ＝ＴＰｈ－－ＰＨＴｒ丹＝一Ｐｈ－－．ＰＨｓｉｎ２０（７）ｏ～－ｐＯｎ图３拉伸、压缩地应力下井壁受力模型直径为Ｄ处的应力边界条件为：ａｒ＝０，ＴｒＯ＝０（８）由半逆解法，可设田为ｒ的某个函数乘Ｕ２ｃｏｓ２ｐ，设Ｔｒ日为ｒ的另一个函数乘１２２ｓｉｎ２０，则应力函数为：咖＝厂（ｒ）ＣＯＳ２０（９）带人极坐标下的相容方程（甭０２＋；昙＋嘉杀）２咖＝０（１０）求解可得：咖（ｒ，８）＝ｃ。ｓＰＨ２８（一Ｐｈ－－４ｒ２＋ｄ２—Ｐｈ－－—ＰＨ—ｄ４—Ｐｈ－－—ＰＨｌｌ）（１Ｉ）根据应力函数与应力之间的关系：＝；石ａ，ｂ十一１塑，２；石十一一，Ｔ－－ｒａＯｒＴ，等ｃｒｏ－ｒｒ２ａｒｚ＝２雨，ｒＯ＝一杀Ｇ翔２一石ｂ剖（１２）可得：一一Ｐｈ－－ＰＨｃｏｓ２８（１一祟）（１一再３ｄ２）Ｊ，铂。＝～Ｐｈ－－：ＰＨｓｉｎ２０（１一景）（１＋４３ｄ万２）（１３）一一一Ｐｈ－－ＰＨ％３２一—了一厶ｃｏｓ２０ＣＯＳ２（１＋羔Ｊｋ）Ｉ、ｌ＋瓦ｉｕＪ２．４非均匀地应力作用下井壁应力分析假设井壁在上述３种载荷单独作用下变形较小，在弹性范围内井壁应力可以采用叠加原理，得到非均匀地应力作用下井壁径向应力、环向应力和剪应力的分布公式：听＝一ＴＰｈ＋ＰＨ＼（１一杀）＋ＴＰｈ－－ＰＨｃｏｓ２臼×１一景）（１一丝４ｒ２、，／一巩景ａｏ＝一ＴＰｈ＋ＰＨ＼（１＋翱一ＴＰｈ－－ＰＨｃｏｓ２ｐｆ３ｄ４＼ｄ２×Ｉ１＋再Ｊ＋仇甭２一——ｉ—ｚｓｉｎＦｌ＿ｌ一蕊Ｊ、＾■），×（１＋甭３ｄ２）（１４）当，－－－９００时，在距离井眼内壁无穷远处的径向和环向剪切应力分别为：，２熙听２一—丁一＋——■厂一．Ｐｈ＋ＰＨ（ｐｈ—ＰＨ）ＣＯＳ２０盯ｒｒ＿÷∞ＺＺ盯ｐ２ｒｌｉ．．ｍ。ａｏｒ－÷∞Ｚ‘Ｔ’ｒｐ＝一—（Ｐｈ—－－ＰｉＨ）一ｓｉｎ２一——ｒ一—————ｉ———一２０（１５）当ｒ＝ａ／２时，井眼内壁的径向和环向的剪切席力分别为：李冬梅等．非均匀地应力下井壁ＹｏｎＭｉｓｅｓ应力及稳定性分析盯，，ｒ＝Ｐｆ，ｆ”阳＝０盯”日＝（Ｐｆ—ｐｈ—ＰＨ）一２（ｐｈ—ＰＨ）×ｃｏｓ２０（１６）山！础３非均匀地应力下井壁等效应力实例分析设井眼直径为２１５．９ｍｍ，水平最大地应力１１０ＭＰａ，最小地应力１００ＭＰａ，井眼内压力３０ＭＰａ。井壁附近岩石处于三向压应力状态，可采用第四强度理论进行校核，以下计算的等效应力均为ＶｏｎＭｉｓｅｓ应力。三向应力包括径向应力、环向应力和垂直方向地应力。根据平面应变理论，垂直方向地应力ａｚ＝肛（听＋ａ０）。图４为应用本文前分析公式，得到的非均匀地应力下井壁附近岩石的三维等效应力分布。Ｒ目较神图６不同角度的ＶｏｎＭｉ￥ｅＳ应力井壁内压和钻井液密度成正比，钻井液密度越大，井壁内压也越大。因此，钻井液密度会影响井壁岩石的等效应力。根据非均匀地应力作用下井壁径向应力、环向应力和剪应力的分布公式（１４），可以得到以钻井液密度为自变量的井壁最大Ｖｏｎ应力的二次拟合公式为：Ｐｅｍａｘ＝１７．４２４２ｙ２—１１３．１２１２ｙ＋２２０．７９３９Ｍｉｓｅｓ（１７）式中：ｙ为钻井液密度，ｇ／ｃｍ３。假设钻井液密度为１．９９／ｃｍ３，采用公式（１７）计孟Ｗ赣淞捌鞍算得到井壁最大等效应力为６８．８ＭＰａ，采用应力理∞∞加∞蚰∞加啪论公式计算的等效应力为６７．１ＭＰａ，误差为２．５％。因此，根据钻井液密度，采用拟合公式计算井壁最大等效应力是可行的。根据得到的井壁岩石抗压强图４三维ＶｏｎＭｉｓｅｓ应力度Ｐｒｏｃｋ和钻进安全系数孔，可得到许可井壁最大ＶｏｎＭｉｓｅｓ应力Ｐ。ｍａｘ的计算公式：Ｐｅｍａｘ图５是不同半径处ＶｏｎＭｉｓｅｓ应力，图６是不同角度处ＶｏｎＭｉｓｅｓ应力。可以看出，离井壁越远，ＶｏｎＭｉｓｅｓ应力越小；９０。方向即垂直于最大水平地应力方向上的ＶｏｎＭｉｓｅｓ应力最大，０度方向上的ＶｏｎＭｉｓｅｓ应力最小；距离井壁越近，ＶｏｎＭｉｓｅｓ应力在不同方向上的差距越大；随着距离井眼中心距离的增加，各个方向的ＶｏｎＭｉｓｅｓ应力趋于相同。最大Ｖｏｎ２＿Ｐｒｏ－ｃｋ（１８）将许可井壁最大ＶｏｎＭｉｓｅｓ应力带人公式（１７），可以计算出不会导致井壁岩石破裂的钻井液密度Ｙｍａｘ０４结论（１）井壁最大ＶｏｎＭｉｓｅｓ应力与水平方向最大地应力和最小地应力及井筒内压有关，随距井眼中心距离的增加，不同方向ＶｏｎＭｉｓｅｓ应力指数减少。（２）垂直于最大地应力方向上的井壁岩石ＶｏｎＭｉｓｅｓ应力最大，０度方向最小；距离井壁越近，不Ｍｉｓｅｓ应力发生在井壁９０。方向上，数值为１５８ＭＰａ；如果该值小于井壁岩石的抗压强度，则井壁是安全的，不会发生破裂。乱皇日同方向上ＶｏｎＭｉｓｅｓ应力的差距越大；随距井眼中心距离的增加，各个方向的ＶｏｎＭｉｓｅｓ应力趋于相同。（３）采用以钻井液密度为自变量的井壁应力拟合公式计算的ＶｏｎＭｉｓｅｓ应力误差不超过２．５％，利用该公式可以确定钻井过程中井壁岩石不会破裂的钻井液密度。图５不同半径处的ＶｏｎＭｉＲ倒赣神ＳＯＳ应力（下转第８９页）罗帅．华北地区ＳＹｌ井钻井技术难点及对策·８９·［２］张金川，金之钧，袁明生．页岩气成藏机理和分布【Ｊ］．天４结论与建议【３］３然气工业，２００４，２４（７）：１５—１８．刘洪林，王红岩，刘人和，等．中国页岩气资源极其勘探潜力分析［Ｊ］．地质学报，２０１０，８４（９）：１３７４—１３７８．（１）钻头选型对缩短建井周期、降低施工成本意义重大。ＳＹｌ井上部地层全面钻进过程中机械钻速虽然较高，但是在钙质泥岩及硅质胶结砂岩段等可［４］肖金裕，杨兰平，李茂森，等．有机盐聚合醇钻井液在页岩气井中的应用［Ｊ］．钻井液与完井液，２０１１，２８（６）ｌ２１—２３．钳｜生差的地层，钻头选型应进一步提高其攻击性。（２）井壁失稳是页岩气探井施工的主要难题。应将防止黏卡、缩径、垮塌等复杂情况作为钻井液设计的核心，提高钻井液的抑制性，尤其是在长裸眼段长时间钻进中应着重提高钻井液的防塌能力。如条件允许，可改为三开井身结构，或在上部非取心井段尝试使用气体钻进技术。（３）取心阶段机械钻速尚有进一步提升空间。取心过程中对比了多种型号取心钻头，机械钻速多维持在０．３０一ｏ．６０ｍ／ｈ，与地矿勘查单位相比取心速度不高，可借鉴其成熟的取心工艺，能进一步提高取心段的钻进速度。参考文献［１］邹才能，董大忠，王社教，等．中国页岩气形成机理、地质特征及资源潜力［Ｊ］．石油勘探与开发，２０１０，３７（６）：６４１—６５３．［５］魏晓亮，张金川，党伟，等．ＳＹｌ井海陆过渡相页岩发育特征及其含气性［Ｊ］．科学技术与工程，２０１６，１６（２６）：４２—５０．［６】６王志刚．涪陵页岩气勘探开发重大突破与启示［Ｊ］．石油与天然气地质，２０１５，３６（１）：１—６．［７】郭少斌，付娟娟，高丹，等．中国海陆相交互页岩气研究现状与展望［Ｊ］．石油实验地质，２０１５，３７（５）：５３５—５４０．［８］王社教，李登华，李建忠，等．鄂尔多斯盆地页岩气勘探潜力分析［Ｊ］．天然气工业，２０１１，３１（１２）：４０—４６．［９】李中明，张栋，张古彬，等．豫西地区海陆过渡相含气页岩层系优选方法及有利区预测［Ｊ］．地学前缘，２０１６，２３（２）：３９—４７．【１０］孙少亮．中长筒保形取心技术在页岩气井中的应用［Ｊ］．钻采工艺，２０１３，３６（５）：１１ｌ一１１３．编辑：赵川喜（上接第８５页）４５—４６．［７］冉小丰，王越之，贾善坡，等．基于渗流一应力一损伤参考文献耦合模型的泥页岩井壁稳定性研究［Ｊ］．中国科技论文，２０１５，１０（３）：３７０３７４．［１】胡勇科，李丹，邱元瑞，等．南堡３号构造井壁失稳问题分析与对策叨．石油地质与工程，２０１４，２８（６）：１０１—１０３．【８］光新军，陈曾伟，刘修善，等．三向应力状态下层理性泥页岩地层井壁稳定性研究［Ｊ】．科学技术与工程，２０１５，１５（５）：５４—５７．［２］蒋晓红，刘其明，衡俊鹏．川西大邑地区气体钻井井壁稳定性评价研究叨．石油地质与工程，２０１１，２５（２）ｉ９５—９７．【９】何朋立，郭力，王在泉．深厚表土层井壁温度竖向附加应力分析【Ｊ］．金属矿山，２０１５，（４）：１１９—１２２．［１０］欧焕农．应力状态对井壁失稳安全性的影响分析［Ｊ］．科学技术与工程，２０１３，１４（６）：１５８８１５９１．［３］周国庆，程锡禄．特殊地层中的井壁应力计算问题［Ｊ］．中国矿业大学学报，１９９５，２４（４）：２４—２５．【４】姚直书，李瑞君．考虑竖向附加力时井壁应力计算方法［Ｊ］．东北煤炭技术，１９９７，４：３—４．［５】李军，柳贡慧，陈勉．正交各向异性地层井壁围岩应力新模型［Ｊ】．岩石力学与工程学报，２０１１，３０（１２）：２４８１—２４８２．【１１］郭莉霞，蒋晓红，黄振华．完整性与非完整性井壁应力模拟分析研究［Ｊ］．钻采工艺，２０１１，３４（５）：１０一１２．［１２］王渭明，路林海．正交各向异性复合井壁应力变形分析与应用［Ｊ】．力学与实践，２００９，３１（１）：５２—５６．［６］王建军，林凯，王新虎，等．易坍塌地层椭圆形井眼内套管应力的有限元分析闭．石油大学学报，２００４，２８（２）：编辑：蒲洪果
