重复压裂地应力场数值模拟计算

1、重复压裂地应力场数值模拟计算张红妮1 何 坚21吉林油田分公司采油工艺研究院 吉林松原 （138000）2吉林油田公司信息管理部 吉林松原 （138000） 摘要： 经过水力压裂后的油气井，由于受压裂工艺、材料的限制，规模欠小等原因会导致水力裂缝导流能力大幅降低而逐渐失去作用，为了获得稳产和经济的开采效益必须进行重复压裂。在己经压裂改造过储层中，由于人工裂缝的存在及油气井长期的生产活动，导致地应力发生变化。通过建立数学模型模拟求解，得出初次压裂和采油注水产生的诱导应力的大小，且在原最小水平主应力方向产生的诱导应力远远大于在最大主应力方向，应力发生反转，重复压裂产生新缝。关键词：重复压裂 诱导应

2、力 模拟计算 应力转向中图分类号：TE3488引言初次人工裂缝产生后，油气井长期的生产活动将在井眼和初次人工裂缝周围的椭圆形区域内导致局部孔隙压力重新分布，改变油藏中的应力分布状况。因此，当重复压裂井中的诱导应力差足以改变地层中的初始应力差时，则在井筒和初次人工裂缝周围的椭圆形区域内发生应力重定向3。根据弹性力学理论和岩石破裂准则，裂缝总是沿着垂直于最小水平主应力方向起裂，那么重复压裂产生新裂缝将可能垂直于初次裂缝缝长方向起裂和延伸，一直延伸到椭圆形的应力重定向边界处(应力各向同性点)，超过应力各向同性点后，应力场方向恢复到初始应力状态，重复压裂新裂缝将逐渐重新转向到平行于初次裂缝缝长方向继续

3、延伸。1、地应力模型建立与计算为了研究裂缝转向，需要得到井眼周围有效应力场分布。井眼周围应力场可由三部分组成：（1）远场引发应力的原地应力场； （2）人工裂缝产生的诱导应力场；（3）流体流动引发的孔隙压力诱导应力场。1.1 原始地应力的计算1.1.1 建立模型垂向地应力是由上覆地层重力引起的，它是随着地层密度和深度的变化的，在深度H出的垂向应力为： （1-1） （1-2） （1-3）1.1.2 模拟计算在地应力计算公式中，构造应力系数和是未知的，必须建立一种方法来确定它们。具体的做法是用用压裂资料计算某一深度的地应力，再根据求出的地应力值反算构造应力系数和。 （1-4） 假设： 一口井所对应的

4、压裂数据可以列一个方程，但其中有两个未知数，如果能够获得两个以上井的破裂压力数据，由方程（14）可以构成一个方程组，其形式为： （1-5）求解方程（2-5）需要至少两组测量数据，当测量数据大于未知数个数时，将构成一个超定方程组，可用最小二乘法解，使解的误差最小化。最终可以求得最大和最小地应力。将确定好的和 代入下式，即可确定和值。 （1-6）将求到的该地区的构造应力常数和目的井的参数代入（2-2）和（2-3）式，即可求得该井的原始最大最小主应力。1.2 人工裂缝引起的诱导应力1.2.1 建立数学力学模型 假设无限大储层中含有一条对称双翼的垂直裂缝，水力裂缝人工诱导地应力场属于平面应变问题，根据

5、弹性力学理论，建立平面应变问题的平衡微分方程、平面问题的几何方程和平面问题的本构方程为数学模型。模型的边界条件为：在0，处： ，0；在0，处： ，在处： ，。1.2.2模拟计算对上述模型根弹性力学按应力求界平面应变问题并进行傅立叶积分变化，求得二维垂直裂缝所诱导的应力场为： （2-10） （2-11） （2-12） 式中：为裂缝的半长。1.3 流体流动引起的诱导应力1.3.1建立应力应变模型假设岩石为具有小变形的弹性介质，储集层和流体等温。应力应变模型基于弹性理论的三个基本关系式组成：应力平衡关系、应变位移关系和应变应力压力关系。用i, j, k分别代表坐标系统中的x, y, z轴。应力平衡方

6、程： （2-13）式中：1代表方向，2代表方向，3代表方向。应变位移方程： （2-14）应变应力压力方程： （2-15） i¹j 式中：为应变增量，为有效应力增量，为流体压力增量，为位移矢量增量，为Boit孔隙弹性常量，为杨氏模量，是剪切模量，为泊松比。1.3.2模拟计算对于一维情况下，在r方向的平衡方程（位移方程）： （2-16） 边界条件： 当 当 式中： 由井眼压力和流动压力引起的有效径向应力和有效切向应力。有效径向应力： （2-17） 有效切向应力： （2-18）1.4初次人工裂缝储层总应力场分布叠加前面提到的三组应力场（原地应力、人工裂缝诱导应力和流动有效应力），得到含有初

7、次人工裂缝的总应力场。垂直裂缝井重复压裂前时间和空间上的应力场为： 如果井壁处的初始最大水平应力方向变为当前的最小水平应力方向，初始最小水平应力方向变为当前的最大水平应力方向，则重复压裂将产生新裂缝。2 实例计算2.1 区块的基本状况新民油田开采的主力油层为扶余油层及杨大城子油层，油层埋藏深度在10501400米。新民油田属岩性断块油藏，断层为近北东向和南北向同生断层，主力砂体呈北东南西和北西南东向延伸。根据上述原因，同时考虑到注采井网的完整性，对于新民油田采用了初期采油速度较高的反九点法面积注水的注采方式。新民油田井排方向的选择为垂直断层，沿着砂体，避开裂缝，使各项水井见水相对均匀。选取新民

8、某区块的*井，模拟计算该井地应力大小。*井是2005年投产新井，本井11、12小层砂体厚度10.0米、静态发育较好，为主力油层。初产液9.5吨、产油6.0吨、含水36.84%，目前日产液1.3吨、日产油1.0吨、含水23.08% 。产液量下降了7.3倍，且注水见效程度不明显。2.2 目的井地应力计算2.2.1原始地应力计算 确定上覆地层应力根据前面给出的上覆地应力的确定方法，有某区块的密度测井曲线确定上覆地层应力。测量井段为在11001400m，每隔1m取一个深度上对应的密度值。各深度下的密度数据离散点，以及对些数据点进行回归，其回归关系公式： （3-1） 式中：为随井深变化的岩石密度，g/c

9、m3； 为地层深度，m。对式（3-1）进行积分，可以得到不同井深的上覆岩层应力的计算公式。 （3-2）式中：为上覆岩层应力，MPa；为地层深度，m。确定静态泊松比根据民* *井的岩心进行的三轴试验，将声波测井得到的动态泊松比转换成静态泊松比。由于测井数据全井的连续，所得到的静态泊松比也是连续的，其回归关系公式如下： （3-3） 式中：为岩石静态泊松比；新民油田某区块地应力计算将新民油田某区块进行16网格剖分，根据该区块压裂井资料，计算该区快最大水平构造系数 1.187，最小水平构造应力系数0.773，地应力计算公式为 （3-4） （3-5） 民* *井目的层段是1108.4-1119.2米，原

10、始地应力场计算结果如下表： 表1 新油田民* *井原始地应力表2.2.2人工裂缝诱导应力本井2004年11月压裂投产，初次加砂20m3 ,预计半缝长为79米左右，穿透比0.75（裂缝高度21米左右）民* *井初次裂缝人工诱导应力场计算结果如下表：表2 新油田民* *井人工裂缝诱导应力增量表最大水平主应sHmax （MPa）21.2最小水平主应sHmin （MPa）15.2井眼半径 rw （m）0.10原始地层压力 PO （MPa）11.4地层渗透率 K （mdc）5.4杨氏模量 E （MPa）28000泊松比 µ0.24孔隙度 0.25压缩系数 Ct0.0008流体黏度 v （mp.

11、s） 7.5流体黏度 v （mp.s） 7.5最大水平主应sHmax （MPa）21.2最小水平主应sHmin （MPa）15.2井眼半径 rw （m）0.10注水压力 P （MPa）25.0原始地层压力 PO （MPa）11.4地层渗透率 K （mdc）5.4杨氏模量 E （MPa）28000泊松比 µ0.24孔隙度 0.25压缩系数 Ct0.0008流体黏度 v （mp.s） 7.5表3 民* *井计算参数表 表4 4口注水井计算参数表 2.2.3注水压力诱导应力计算选取新民* *井作为研究区域，周围有注水井4口，产油井为定压生产，采用的生产井和注水井计算参数如表4、5根据两表中

12、参数算出孔隙压力在裂缝方向上（x）轴和过井轴垂直于裂缝方向(y轴)上变化的孔隙压力分布，以及相应的最大主应力方向的应力变化和最小主应力方向的应力变化。民* *井注水压力诱导应力计算表5 新油田民* *井区注水采油诱导应力增量民*井总应力场力学计算结果如下表：表6 新油田民* *井总诱导应力表民* *井在开发一段时间之后，由于初次压裂造缝和注水采油均产生了诱导应力，且原最小主应力方向的诱导应力大于原最大主应力方向的诱导应力，使原最小主应力方向的应力大于原最大主应力方向的应力，两相应力差值在井筒附近为3Mpa左右，离井筒越远差值越小。裂缝总是沿着最大主应力方向延伸，所以理论上民* *井重复压裂将产

13、生于原裂缝有一定夹角的新裂缝3、几点认识（1）初次人工裂缝诱导应力在最小水平主应力方向（垂直裂缝方向）诱导水平应力大，在最大水平主应力方向（裂缝延伸方向）诱导水平应力小；张开裂缝诱导应力在裂缝面上最大，离缝越远诱导应力越小。（2）长期注水采油导致储层应力场发生变化，并且这种变化随着生产时间增加而增加，在空间上离井筒越近变化量越大，而后随距离的增加变化减小，在一定的远处地应力不发生变化。（3）在垂直初次裂缝方向和裂缝延伸方向，最大主应力的房的应力变化要最小主应力方向的应力变化大，这种变化的差异随距离井筒增加，越来越小，在离裂缝足够远处应力场几乎不受影响。 参考蔡美峰，等主编岩石力学与工程第四章岩

14、石本构关系与强度理论（185-191页）The mumerical simulate and calculate of terrestrial stress field when refracturingZhanghongni He jianPeturleum Production Technology Research Institute ofJilin Oilfield Company Abstrct : Because of the limit of fracturing technics and materials , the flow conductivity of the frac

15、ture will decline day by day after the stimulate of hydralic fracturing .For the sake of stable yileds and ecnomical explorition benefits,we must take to hadraulic refracturing. Crustal stress of The reservoir which had hydraulic fractured will change , as the artificial fracture and the long time o

16、il production .The induced stess of the first time fracture, manufacture and flooding will be calculated with the establishing of mathmatics modle. Also ,the induced stress in the aspect of primary mini-horizon pracinpal stress is more bigger than in the aspect of primary major horizon pracinpal stress . The new fracture will be build up when refracture as the stress turn invert . Key words : refracture