软泥岩井眼弹塑性变形的拉格朗日元法计算.pdfVIP

三些 趔丛 ，石得勤。，张 利‘ 。石油大学，北京 昌平 102200~ 。中国海洋石油生产研究中心，北京 101149； 。大庆油田钻井所，黑龙江 大庆 163413； ‘国家专利局，北京 100083． 摘 要：首次将大变形原理引人井眼及其围岩弹塑性变形计算，应用拉格朗日元涪对 塑性泥页岩地层的井眼缩径问题进行分析求解。拉格朗日元法能够计算材料非线性 和几何非线性问题，采用常用的摩尔一库仑强度准则对所选择的计算模型进行拉格朗 日元舟析计算，求解出井眼的收缩变形、井眼围岩位移场分布以及井眼围岩塑性区的 萎主盏急型；苎量塑!巫壁望丝銮形；；鳞咄 存s—E-。-- 分类号关键词：整韭兰；苎量塑!垂!壁望丝立丑；‰ 扦 ． 雾衫∥矛渺计 ：0344．7 文献标识码：A 、 f ’ ’ 』 1 拉格朗日元法计算原理 拉格朗日元法(FastLagrangianAnalysisofContinua)计算程序…，即连续介质的快速拉 格朗日元分析，简称FLAC，用于解决工程力学计算的有限差分程序，适合于计算材料大变形 问题，而岩石力学问题也多为非线性大变形问题。倒如，石油钻井过程中遇到的塑性软泥岩地 层的井眼缩径问题就是大变形问题。下面以一维问题为倒推导拉格朗日元法差分方程的建立 过程。假定一个杆为弹性，侧面不受限制，两端受拉，杆的密度为P，弹性模量为E，a。为应力， “，为位移，如图1所示。由本构方程：％一￡警，和杆的运动平衡方程：p磬=孥，可以建 立一维拉格朗日元法的位移“，的差分方程：“：(f+△f)一“：(f)+“：(f+拿)△f。 对上式用显式法迭代求解。对于某一个结点而言，在每一时刻它受到来自其周围区域的合 力的影响，如果合力不等于零，也就是具有了失稳力，根据牛顿定律，结点就在产生加速度，进 收稿日期：1998一02—15 作者简升：邓金根(1963一)，男，石油大学(北京)石油工程系教授、博士，中国岩石力学与工程学会理事．主要从事石油工 程岩石力学方面的研究。 33 而可以在一个时步中求得速度和位移的增量。对于每一个区域而言，可以根据其周围结点的运 动速度求得它的应变率，然后可以根据材料的本构关系求得应力的增量，再撼加下一级载荷。 如此循环计算直至问题收敛。 2 强度准则 钻井实践和实验室研究表明，多 数岩石的破坏一般表现为剪切破 坏。设剪切面上的剪应力为r，法向 应力为一，根据库仑等人的研究。]， 认为岩石失稳破坏时，剪切面上的一 剪切应力必须克服岩石的固有强度 c值，也称作粘聚力，再加上作用于 剪切面上的摩擦力(P·a)。即： r= C + p·口 其中：P为岩石的内摩擦系数，P一 由于岩石是孔隙介质，根据孔隙弹性理论，应采用有效应力进行井眼围岩弹塑性变形计 算。由地层孔隙压力实测数据可知：地层孔隙压力梯度G =0．0104 MPa／m。对于塑性泥页 岩，我们取Biot有效应力系数为0．5进行计算，可以求得有效最大水平地应力 一和最小水平 地应力 ： 一日一 H aPP— GHH — aGppH 一 28．2MPa h= Oh— aPp— GhH — aGppH = 19．8MPa 衰1 南海油田乐东9-1·1井岩石的 帖聚力c和内摩接角 Table 1 Cohesive strength and Internal friction angle of the rock in Ledong 8-1—1 well in the Nanhai sea oilfield 图3 计算的