基于测井资料的岩石力学参数有限元反演分析.pdfVIP

长江大学学报 (自然科学版) 2011年6月 第8卷第6期 JournalofYangtzeUniversity(NatSciEdit) Jun．2011，Vo1．8No．6 基于测井资料的岩石力学参数有限元反演分析 姬 辉 ， ． 楼一珊，翟晓鹏，史文专，秦杜( 兰 。) 史连超 (安东石油技术集团有限公司，北京1001o2) [摘要]提 出了一种利用有限元反演套变地层岩石力学参数 的分析方法。将地层、水泥环、套管作为研究 对象，建立套变地层 中套管变形与岩石力学参数的关系模型，利用多臂井径测井资料进行参数反演 ，结 合摩尔库伦准则确保反演结果的唯一性 。将该方法应用到中原油 田卫 95—86井建立相应的有限元模型进行 反演计算，反演结果与实际相符。 [关键词]岩石力学参数 ；测井资料 ；套管变形；摩尔库伦准则；反演计算 [中图分类号]TE21 [文献标识码]A [文章编号]l673 油井套管在生产过程中由于承受地层应力的作用 ，极易产生套管变形。套管变形与套变地层 的岩石 力学参数有着密切的联系。然而至今还没有能直接获取套变地层岩石力学参数的方法 ，只能通过模拟计 算得到。常见 的计算方法有最小二乘法 、遗传算法 等 ，但计算方法复杂 ，求解困难。为此 ，笔者提出 了一种利用有 限元反演计算岩石力学参数的方法 。 1 套管变形与地应力和岩石力学参数的关系 套管变形往往是 由于套管受不均匀挤压造成 的_2]，表 现为套管在最大水平地应力方 向为压缩变形，在最小水平地 应力方 向为拉伸变形 ，最终由圆形变为椭圆状。由于岩层一 般较厚 ，根据厚壁圆筒理论 ，将此问题处理成平面应变 问题 加以分析 ，取图 1所示地层一水泥环一套管系统为研究对象 ， 为最大水平地应力 ； 为最小水平地应力。假定套管及水 泥环都为理想 圆形，厚度均匀 ，r为套管内半径 Ir2为套管外 半径 ；F3为水泥环外半径。 根据弹性力学厚壁圆筒平面应变问题 的求解方法 ，极坐 标系中的应力分量表达式为： 图1 套管一水泥环一地层力学模型 一 +等一[一a2一a2()4D(a，一)。]cos2 一 碧一[2A+6B()+ ()。]cos2 ㈩ r一一 + 一[2A一6B()6C(a一)2D／a，]sin2 式 中， 为内边界周向应力分量 ； 为内边界径 向应力分量 ； 为内边界切向应力分量 ；r为井眼周围岩体 任一点到0点的距离 ；为井眼周围岩体任一点与最大水平地应力所成的角 ；为应力函数，一声 ；A、B、 c、D为与弹性参数和几何参数相关的常数 ；n为边界半径。 [收稿 日期]2011—04～27 [作者简介]姬辉 ，男，硕士生，现主要从事岩石力学和套管损伤方面的研究工作。 第 8卷 第 6期 姬辉等：基于测井资料的岩石力学参数有限元反演分析 对于平面应力 问题本构关 系为 ： fe一car--yaO {e一吉(a8一 ) (2) f79．： = 式中，e，e，7 分别为周 向应变 、径向应变和切向应变。 根据平面应变问题的平衡方程 、物理方程和几何方程 ，以及地层外边界 r一 。。时应力边界条件，可 得套管的应力分量为： f口，一 (2A1+Az厂 )一 (A5+ 4A3r- 一 3A6r-)cos20 lde一 (2A1一A2r)+(A5+12A4r。一3A6r-)cos2O (3) r【= (A5+6A4r。一2A3r-+3A6厂)sin20 套管位移垂直分量和水平分量分别为 ： I一 E(k一1)A1r—A2r-]+E(