油井围岩应力与孔压场的数值分析.pdfVIP

l工程力学》增刊 1999年 ·669· 油井围岩应力和孑L压场的数值分析 金国东 赵忠良 中匿科学院力学新7童l∞0柏) (大虎油田蕈八采油厂规捌设计研姗 摘要：本文基于比奥的同结耦合理论，在考虑岩层的非线性特性及流体特性的情况下，对深井周围岩 层的应力和孔压的分布进行了有限元分析，提出了适合实际应用的数值方法．通过数值解与理论解的比 较以及模拟现场实际问题，说明了数值模型是正确有效的．结果表明：模型能够精确地预舅藏逮低于临 界流速条件下应力和孔隙压力的分布，同时也能预测流速超过理论临界值时地层的全面的拉伸破坏．为 进一步分析和了解在渗流条件下固体飘粒通过未固结可变形的多孔介质漉人井眼的机理以及岩层破坏的 机理提供了相应的数值基础． 关键词：有限元法同结应力孔压岩层 在石油天然气工业中，井壁稳定性一直是人们所关注的问题．在井黜q钻进过程中， 井壁的不稳定会引起井壁的坍塌、破裂，甚至井眼的报废，造成钻井生产费用的增加；在 采油阶段，地层的不稳定会引起砂的流动或出砂．而出砂又会加剧套管的磨损。在严重的 情况下，会导致油井生产的停止或油井的报废．另外，在整个采油过程中，井眼压力需保 持一个最小值，以便达到产量最大，除非采油是在对涉及相关的物理过程有一个充分了解 的基础上进行的，否则，也会引起井壁的失稳及油井出砂．尽管这些稳定问题在石油工业 中经常遇到，然而有关这些问题的机理和过程，人们至今尚未弄清楚．这是因为人们对井 眼周围岩层的应力和孔压分布的实际情况缺乏足够的了解，同时缺乏对井壁稳定问题作出 综合评价所需的足够详细的定量观察数据． Rjsnes，Bra廿j＆H跏d…等人提出了在稳定渗流的条件下，弹性材料中圆形洞附近 的应力和孔隙压力的分布的理论公式．这些公式解释了塑性剪切破坏的形成，在拉伸破坏 开始之前，结果是正确的．同时黜mes等人也提出了在整个地层中引起拉伸应力所需的Il缶 界流动速率的表达式．他们认为：在速率到达临界值之前，尽管地层中可能有塑性破坏存 在，地层仍能保持稳定．一旦流动速率超过临界值，整个地层就变得不稳定，固体颗粒开 始流动．根据Rjsnes等人的理论公式，开始出砂要求流速超过某一临界值qc．如果井眼 没有下套管，地层能够支持的最大流速为： n。珈扣乜州s+警，Z ㈣ Ⅳ 式中D为产层厚度，中为摩擦角，c为粘聚力，k为固有渗透率，“为流体粘度． 盒国东男19729出生顿士助研 ·670· 《工程力学》增刊 1999年 方程(1)描述的稳定准则只适用于稳态流动条件．在大多数现场问题中，由于流体流动 的短暂性，井壁附近的流体压力梯度开始可能大大超过长期稳态值．这会导致早期的失 稳，而这在方程(1)中则没有考虑．随周围材料的坍塌，坍塌区域的渗透率会增大，假设没 有其它变化，压力梯度也会减小．当压力梯度减小到临界值的下限时，出砂就会停止．稳 态条件到达时，地层特性可能已发生重要的变化． 本文基于比奥的固结耦合理论“1，在考虑岩层的非线性特性及流体特性的情况下，对 深井周围岩层的应力和孔压的分布进行了有限元分析，为进一步分析和了解井眼周围岩层 的失稳机理，及油井出砂机理等提供了可靠的数值方法． 基本方程 在深井周围岩层中取一该分体，只考虑重力．z坐标向下为正，应力以压为正．假设 岩层均匀，各向同性，流体可压缩． a)平衡方程：疗u．j—Fi=O (2) 式中盯u是总应力矢量，Fi是体积力矢量． b)有效应力原理：盯；=盯u—u民 (3) 式中盯u是有效应力矢量·u是孔隙流体的压力，占u是Kmn∞k盯’符号． c)应力应变关系：盯；=D睦l毛I (4) 式中skl是应变矢量，D哇I由切线刚度模量E。、B。等构成．模量变化如下【3】…： 切线变形模量：n。=。^。[詈]吨[，一号爱笔主≥端]2