基于返排数据计算页岩气井 压裂有效体积的方法及应用.pdf

引用本文：王益民，李继庆，万云强，等基于返排数据计算页岩气井压裂有效体积的方法及应用石油实验地质，（）：，（）：基于返排数据计算页岩气井压裂有效体积的方法及应用王益民，李继庆，万云强，刘 莉，张 谦，朱朝光，汤亚顽中国石化 江汉油田分公司 勘探开发研究院，武汉；中国石化 江汉油田分公司，湖北 潜江 摘要：页岩气多段压裂水平井的缝网参数通常利用产气数据解释获得。虽然对产气数据进行解释能得到裂缝参数，但很难获得压裂有效体积。为快捷地获取页岩气井压裂有效体积，该研究建立了页岩气多段压裂水平井压裂液返排数学模型，并结合水相物质平衡方程和渗流方程，推导获得了页岩气井多段压裂水平井压裂有效体积表达式。通过分析模型解发现，在边界控制流阶段产量规整化压力与物质平衡时间双对数曲线为单位斜率直线，利用该阶段产量规整化压力及物质平衡时间数据，可以计算页岩气井压裂有效体积，从而形成了一套基于压裂液返排数据计算页岩气井压裂有效体积的方法。实例应用表明：（）构建的页岩气多段压裂水平井压裂液返排模型能快速地计算页岩气井压裂有效体积，且计算结果可靠；（）计算页岩气井压裂有效体积不能忽略试气期间的压裂液返排数据，否则计算结果偏小；（）形成的页岩气井压裂有效体积计算方法还能识别邻井压裂干扰，并定量化表征邻井压裂干扰对页岩气井压裂有效体积的影响。该研究成果为油田现场估算页岩气井压裂有效体积提供了一种新的方法，同时也为油田现场识别邻井压裂干扰提供了新的思路和方法。关键词：边界控制流；双对数诊断图版；返排数据；压裂有效体积；页岩气中图分类号：文献标识码：，；，：，：（），收稿日期（）：；修订日期（）：；出版日期（）：。作者简介：王益民（），男，博士，副研究员，从事非常规油气田开发研究。：。基金项目：中国石油化工股份有限公司科技部项目“页岩气开发调整区动态分析方法与生产规律研究”（）资助。第 卷第 期 年 月 石 油 实 验 地 质 ，；（），；（），：；页岩气藏具有“低孔隙度、低渗透率”的特点，需要通过水力压裂才能进行经济开采。在大多数情况下，页岩气水平井经过多段压裂后，沿着水平井方向会形成复杂缝网。表征裂缝参数对评价压裂效果具有重要意义。目前大多数学者利用产气数据解释得到裂缝参数。虽然利用产气数据能有效地解释获得裂缝参数，但是很难表征完整的裂缝网络，这是因为改造缝网最初充满的是压裂液而不是气体。鉴于改造缝网在放喷前被压裂液填充，研究认为利用压裂液返排数据可以获取压裂缝网信息。利用压裂液返排数据可以解释获得人工缝网导流能力、裂缝渗透率及裂缝半长，还可以解释获得压裂有效体积。等构建了水相渗流模型，该模型能利用返排液数据计算压裂有效体积，因模型未考虑气体压缩性对返排液量的影响，且该模型计算压裂有效体积需要确定裂缝渗透率，造成使用该模型计算得到的压裂有效体积误差较大。和 假设地层流体为单相水流动，结合流动物质平衡法，给出了利用返排液数据计算压裂有效体积的方法。但该计算方法也未考虑气体压缩性对返排液量的影响，造成计算的压裂有效体积偏离实际情况。本文基于页岩气水平井压裂返排特征，考虑页岩中气体压缩性对返排液量的影响，建立并求解了页岩气多段压裂水平井压裂液返排模型，并绘制了产水量规整化压力与物质平衡时间双对数曲线，形成了一套基于双对数诊断图版计算页岩气井压裂有效体积的方法。此外，尽管页岩压裂改造效果受地质条件的影响，但其改造后形成的压裂有效体积均可利用本文提出的方法进行计算，即本文形成的页岩气井压裂有效体积计算方法对不同地质条件下页岩储层具有普遍适用性。本文的研究成果为油田现场估算页岩气井压裂有效体积提供了一种新的方法，同时也为油田现场识别邻井压裂干扰提供了新的思路和方法。压裂液返排机理页岩气井压裂完成后，缝网内充满压裂液。根据能否形成气体流动通道，将裂缝分为两类：有效裂缝，气体能进入该裂缝，并在其中流动；无效裂缝，因裂缝闭合，气体无法进入其中。在有效裂缝内，气体进入其中，并驱替出压裂液。故利用返排数据可以计算得到有效裂缝体积。在返排初期，有效裂缝内只有单相水流动，压裂液处于不稳定流动阶段。随着时间的推移，部分有效裂缝被气体突破，在被突破的裂缝内压裂液流动进入边界控制流阶段，而在未被突破的裂缝内压裂液仍处于不稳定流动阶段；当所有的有效裂缝被气体突破后，此时压力波传播到有效裂缝边界，在 改造区内压裂液流动进入边界控制流阶段。多段压裂水平井压裂液返排模型 物理模型对页岩气藏中一口水平井进行多段水力压裂，其物理模型如图 所示，模型假设如下：（）水平井被 条不可变形裂缝穿透，裂缝半长为、缝宽为，缝高为（图）；（）压裂改造形成的裂缝面为矩形，且裂缝内为单相水线性流动（图）；（）页岩气井压裂有效体积内空间视为恒温定容体；02xfhfwxyz图 页岩气多段压裂水平井物理模型 石 油 实 验 地 质 ：第 卷 2xfhf图 裂缝内流体流动示意 （）忽略重力的影响，地层岩石及流体均微可压缩，流体作等温达西渗流；（）不考虑基质向裂缝的流体补给。数学模型 水相物质平衡方程在没有能量补给条件下，地层主要依靠岩石骨架和流体的弹性能驱替压裂液，故累产水的地下体积等于地层压力下降引起的地下气体、水的膨胀量以及气藏孔隙体积减小量之和。由于 与、方向相反，以体积增大方向为参考，则有：（）式中：为产水速率，；为压力 下水相的体积系数；为页岩气井压裂有效体积，；为气体体积，；为水的体积，。根据水、气及岩石等温压缩系数的定义，可分别得到如下公式：（）（）（）式中：、分别为水、气体、岩石的等温压缩系数，。联立公式（）（）得：()（）对上式进行整理得：（）其中，对公式（）进行积分，则得：（）式中：为原始地层压力，；为平均地层压力，；为累计产水量，；为总压缩系数，。水相渗流方程模型假定裂缝为矩形，压裂液在裂缝中的流线为一组互相平行的直线，则页岩气多段压裂水平井水相渗流微分方程为：（）式中：为水的黏度，；为裂缝渗透率，；为裂缝孔隙度；为裂缝条数；为裂缝半长，；为裂缝宽度，；为裂缝高度，。井底边界条件：（）裂缝末端边界条件：（）联立公式（）（），获得裂缝压力分布表达式：，()（）式中：（，）为裂缝内 时刻 处的压力，；为井底流压，。为构建平均地层压力与井底流压之间的关系，首先计算平均地层压力。模型假设裂缝形状为矩形，则平均地层压力可以表示为：（）（）式中：为单翼缝有效体积，。联立公式（）（），则得：（）对上式进行整理，则平均地层压力与井底流压 第 期 王益民，等 基于返排数据计算页岩气井压裂有效体积的方法及应用的关系如下：（）模型的解联立公式（）和公式（），可以获得原始地层压力与井底流压关系，如下所示：（）对式（）进行变形，得：（）定义 个变量，产量规整化压力：；物质平衡时间：，则式（）可表示为：（）式中：为产量规整化压力，（）；为物质平衡时间，。从公式（）中可以看出，在边界控制流阶段，产量规整化压力与物质平衡时间双对数曲线为单位斜率直线；直角坐标系中，在边界控制流阶段产量规整化压力与物质平衡时间是线性关系，利用直线的截距可以获得裂缝渗透率、裂缝半长及裂缝条数之间的关系，同时利用直线的斜率可以计算页岩气井压裂有效体积，则压裂有效体积表达式为：（）式中：为直角坐标系中产量规整化压力与物质平衡时间关系曲线的斜率。实例应用 页岩气井压裂有效体积计算方法利用本文构建的页岩气多段压裂水平井压裂液返排模型，拟合实测压裂液返排液数据和压力，可以解释得到页岩气多段压裂水平井压裂有效体积，具体计算方法如下：（）收集产水数据及压力数据；（）绘制产量规整化压力（）及物质平衡时间（）双对数诊断图版；（）在双对数坐标中做单位斜率直线，并与实测 双对数曲线拟合，识别边界控制流；（）提取边界控制流阶段产量规整化压力及物质平衡时间数据，在直角坐标系中绘制 曲线，通过线性回归确定最佳匹配的斜率；（）计算综合压缩系数；（）计算页岩气多段压裂水平井压裂有效体积。实例应用以涪陵页岩气田 口典型井为例，利用返排液数据计算页岩气多段压裂水平井的压裂有效体积。实例 分考虑试气期间压裂液返排数据和不考虑试气期间压裂液返排数据 种情况，讨论试气期间返排数据对页岩气井压裂有效体积的影响，并将计算结果与微地震结果对比，验证模型的可靠性。实例 中 井周围有新的压裂井，通过实例 展示利用页岩气井压裂有效体积计算方法识别邻井压裂干扰，并计算邻井压裂干扰对页岩气井压裂有效体积的影响。实例 井为涪陵页岩气田一口多段压裂水平井，该井于 年 月 日完成压裂，试气期间返排液量为 。井于 年 月 日正式投产，水平井长度 ，气层中部垂深 ，初始地层压力 ，总压缩系数为 。分考虑试气期间压裂液返排数据和不考虑试气期间压裂液返排数据 种情况，利用 井产水数据及压力数据绘制产量规整化压力及物质平衡时间双对数曲线（图）。从图 中可以看出，种情况下的产量规整化压力与物质平衡时间双对数曲线均出现了斜率为 的直线段，说明 种情况下压裂液都进入了边界控制流阶段。此外，种情况0.010.11101001101001 00010 000100 000?MPa/(m/d)3/(?d/?图 涪陵页岩气田 井产量规整化压力与物质平衡时间双对数诊断图版 石 油 实 验 地 质 ：第 卷 下的产量规整化压力与物质平衡时间双对数曲线出现了明显的偏移，说明试气期间压裂液返排数据对计算页岩气井压裂有效体积存在较大的影响。提取边界控制流阶段的产量规整化压力与物质平衡时间数据，在直角坐标系中绘制产量规整化压力与物质平衡时间关系曲线（图）。图 中离散数据点为现场实测数据，虚线（趋势线）为拟合结果。从图 可以看出，拟合结果很好，因此可以将趋势线看成产量规整化压力与物质平衡时间关系曲线，利用直线的斜率可以计算页岩气井压裂有效体积。利用公式（）计算不考虑试气期间压裂液返排数据情况下 井压裂有效体积为 ，而考虑试气期间压裂液返排数据情况下 井压裂有效体积计算结果为 。对比 井微地震监测结果（）发现，在不考虑试气期间压裂液返排数据情况下，模型计算的页岩气井压裂有效体积结果偏小，且误差较大（）；而在考虑试气期间压裂液返排数据情况下，模型计算的页岩气井压裂有效体积结果可靠，误差仅为。以上结果表明：计算页岩气多段压裂水平井的压裂有效体积需要考虑试气期间的压裂液返排数据，若忽略该阶段的返排数据，模型计算结果偏小；在考虑试气期间压裂液返排数据情况下，利用本文提出的数学模型可较可靠地计算页岩气井压裂有效体积。实例 井为涪陵页岩气田一口多段压裂水平井，该井于 年 月 日完成压裂，试气期间返排液量为 。井于 年 月 日正式投产，水平井长度 ，气层中部垂深 ，初始地层压力，总压缩系数为 。02468101202 0004 0006 0008 00010 00012 000yx=0.001+0.026 4R=0.995 8?d/?MPa/(m/d)3/(?yx=0.001 7+0.0130 6R=0.991图 涪陵页岩气田 井边界控制流阶段产量规整化压力与物质平衡时间关系曲线 年 月 日，井周边压裂了一口新井 井，井试气期间返排液量为 ，井于 年 月 日投产。对 井的渗流过程分为 个阶段展开研究：阶段 为从 井投产到 井压裂前；阶段 以 井开始压裂为起始点。利用 井 个阶段的产水数据及压力数据，在双对数坐标中分别绘制阶段 和阶段 的产量规整化压力及物质平衡时间双对数诊断图版（图）。从图 中可以看出，双对数曲线均出现了斜率为 的直线段，说明 个阶段的压裂液都进入了边界控制流阶段。此外，阶段 的单位斜率直线向右下方偏离阶段 的单位斜率直线，说明 井压裂对 井产生了压裂干扰，井形成的裂缝与 井的裂缝形成了连通，造成 井的压裂有效体积增大。根据现场动态监测资料显示，井压裂的确对 井造成了压裂干扰。井压裂对 井压裂有效体积的影响可以利用图 中 条单位斜率线上的数据进行量化。提取边界控制流阶段的产量规整化压力与物质平衡时间数据，在直角坐标系中绘制产量规整化压力与物质平衡时间关系曲线（图）。图 中离散数据点为现场实测数据，虚线（趋势线）为拟合结果。从图 可以看出，拟合结果很好，因此可以将趋势线看成产量规整化压力与物质平衡时间关系曲线，利用直线的斜率可以计算页岩气井压裂有效体积。在 井压裂前、后，利用公式（）计算 井的压裂有效体积分别为 和 。受 井压裂影响，井的压裂有效体积在原始体积上增加了。通过实例 不难发现，利用页岩气井压裂有效体积计算方法，可以识别邻井压裂干扰，并定量计算出邻井压裂干扰对页岩气井压裂有效体积的影响。0.010.11101001101001 00010 000100 000?1?2?d/?MPa/(m/d)3/(图 涪陵页岩气田 井产量规整化压力与物质平衡时间双对数诊断图版 第 期 王益民，等 基于返排数据计算页岩气井压裂有效体积的方法及应用01234567891002 0004 0006 0008 00010 000?1?2yx=0.000 9+0.018 6R=0.999 9yx=0.001 1+0.019 5R=0.997 9?d/?MPa/(m/d)3/(图 涪陵页岩气田 井边界控制流阶段产量规整化压力与物质平衡时间关系曲线 结论（）水相进入边界控制流阶段后，产量规整化压力与物质平衡时间呈线性关系，在直角坐标系中，利用拟合线的斜率可以计算页岩气井压裂有效体积，拟合线的截距反映了裂缝渗透率、裂缝半长及裂缝条数之间的关系。（）计算页岩气井压裂有效体积需要考虑试气期间的压裂液返排数据，若忽略该阶段的返排数据，模型计算的页岩气井压裂有效体积会比实际值小。（）通过实例应用，将模型计算的页岩气井压裂有效体积与微地震监测结果对比，误差仅为，验证了本文构建的数学模型计算页岩气井压裂有效体积的可靠性。（）利用页岩气井压裂有效体积计算方法，可以识别邻井压裂干扰，并计算邻井压裂干扰对页岩气井压裂有效体积的影响。利益冲突声明 所有作者声明不存在利益冲突。作者贡献 王益民参与模型建立及推导；王益民、李继庆、万云强、刘莉、张谦、朱朝光、汤亚顽参与论文写作和修改。所有作者均阅读并同意最终稿件的提交。，参考文献：王永胜页岩气水平井产能评价方法研究成都：西南石油大学，：，纪国法，张公社，许冬进，等页岩气体积压裂支撑裂缝长期导流能力研究现状与展望 科学技术与工程，（）：，（）：，：，（）：，：，：，：，：（）：，：，（）：，：，：，：，：，：，（编辑徐文明）石 油 实 验 地 质 ：第 卷