CN110359906A - 一种基于短期生产数据的地层压力计算方法 - Google Patents

Abstract

本发明涉及油气田开发技术领域，具体涉及一种基于短期生产数据的地层压力计算方法。一种基于短期生产数据的地层压力计算方法，通过选取拟稳态阶段生产时间间隔小于1个月的任意两个生产点的地层压力和累计产气量，结合原始地层压力、原始条件下的偏差系数、压力≥13MPa条件下偏差系数与压力之间关系式，建立了任意生产点的地层压力计算方法，代入任意生产点的累计产气量，即可计算对应的地层压力。本发明相较于实测法和物质物质平衡法，降低了关井和测试带来的经济损失；相较于流动物质平衡法避免了近似计算，结果更加符合实际；本发明仅需两个生产点的相关参数，计算过程更加简单。

Description

一种基于短期生产数据的地层压力计算方法

技术领域

本发明涉及油气田开发技术领域，具体涉及一种基于短期生产数据的地层压力计算方法。

背景技术

地层压力作为气藏开发的动力，决定着气井的产气量和经济效益，时刻掌握地层压力的变化，对保障气井和气藏的增产、稳产和后期开发调整至关重要。现有地层压力获取技术包括实测法和计算法两大类：实测法最准确但需要较长时间的关井，长时间关井会给有生产任务的气井会带来经济上的损失；计算法简单快速，物质平衡法和流动物质平衡法为目前最常用的方法，但都具有参数获取困难和计算精度较差等劣势。

物质平衡法采用物质平衡方程计算地层压力，但原始地质储量难以准确获取。实际应用过程中，通过气井多次的长时间关井，测得气藏在不同条件下的地层压力值，基于累计产气量与视地层压力成直线关系，进一步外推得到原始地质储量。物质平衡法要求气井多次关井，对于低渗气藏而言并不现实，不适合大规模应用。

流动物质平衡法认为气井稳定生产一段时间后井底流压曲线和地层压力曲线平行，获取得到曲线平行起点的地层压力，即可计算任意点的地层压力。实际应用过程中，曲线平行起点的地层压力没法获取，因此往往采用原始地层压力代替曲线平行起点的地层压力，进行近似计算，计算得到的地层压力与实际地层压力之间存在较大误差。

发明内容

针对背景技术存在的问题，本发明提出一种简单快速且计算精准的基于短期生产数据的地层压力计算方法。

本发明的技术方案在于：

对于定容气驱且p≥13MPa的气藏，采用定容气驱气藏物质平衡方程进行描述：

式中：p_ri和p_i分别为任意生产点的地层压力和原始地层压力，单位为MPa；Z_ri和Z_i分别为p_ri对应的偏差系数和p_i对应的偏差系数；G_pi和G分别为任意生产点的累积采出气量和原始地质储量，单位为10⁴m³。

当气井进入拟稳态生产阶段，对于p≥13MPa的气藏，地层压力和井底流压的下降速率相同，选取处于拟稳态阶段且时间间隔小于1个月的任意两个生产点，存在：

Δp_r＝Δp_wf＝p_r2-p_r1＝p_wf2-p_wf1 (2)

式中：p_r1和p_r2分别为任意两个生产点对应的地层压力，单位为MPa；p_wf1和p_wf2分别为任意两个生产点对应的井底流压，单位为MPa；Δp_r和Δp_wf为任意两个生产点对应的地层压力差值和井底流压差值，单位为MPa。

拟稳态阶段，任意两个生产点的物质平衡方程存在如下关系：

式中：Z_r1为p_r1的偏差系数，Z_r2为p_r2的偏差系数。

当两个生产点之间的时间间隔小于1个月时，由于时间间隔较短，Z_r1与Z_r2之间相差非常小，可近似认为Z_r1＝Z_r2，取用替代Z_r1与Z_r2，得到：

式中：为平均偏差系数。

对于p≥13MPa的气藏，Z和p可采用线性关系式进行描述：

Z＝ap+b (5)

式中：a为斜率，b为截距。

进而得到：

Z_ri-Z_wfi＝a(p_ri-p_wfi) (6)

式中：Z_wfi为p_wfi对应的偏差系数；p_wfi为p_ri对应的井底流压，单位为MPa。

气井在拟稳态阶段生产压差保持稳定，故Z_r和Z_wf的差值为常数值C：

C＝a(p_r1-p_wf1)＝a(p_r2-p_wf2)＝a(p_ri-p_wfi) (7)

Z_ri＝Z_wfi+C (8)

将(8)式代入到(4)式，得到：

式中：为平均偏差系数，Z_wf1和Z_wf2分别为p_wf1和p_wf2对应的偏差系数。

将(9)式代入到(3)式，得到：

将(10)式代入到(1)式，得到：

将(5)式代入到(11)式，得到：

通过迭代，可计算任意生产点的地层压力p_ri。

一种基于短期生产数据的地层压力计算方法，包括如下步骤：

1)获取气井的原始地层压力p_i和原始条件下的偏差系数Z_i；获取气井处于拟稳态阶段且时间间隔小于1个月的任意两个生产点的井底流压(p_wf1，p_wf2)和累计产气量(G_p1，G_p2)；获取p≥13MPa条件下的偏差系数Z与压力p的关系式Z＝ap+b，得到斜率a和截距b；

2)计算第一个生产点对应的相关参数：

p_wf1和p_wf2对应的偏差系数Z_wf1和Z_wf2的平均值

地层压力p_r1的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r1：

设定相对误差ε≤1％，判断p_r1是否满足满足则采用下述公式计算p_r1对应的C₁：

C₁＝a(p_r1-p_wf1)；

3)计算第二个生产点对应的相关参数：

地层压力p_r2的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r2：

设定相对误差ε≤1％，判断p_r2是否满足满足则采用下述公式计算p_r2对应的C₂；

C₂＝a(p_r2-p_wf2)；

4)采用下述公式计算任意生产点的地层压力p_ri：

其中，C为C₁与C₂的平均值，G_pi为任意生产点的累计产气量。

优选地，若p_r1不满足即时，则进行迭代：以p_r1作为计算及p_r1，若再次判断重复迭代直至获取最终的p_r1，并计算p_r1对应的C₁；

若p_r2不满足即时，则进行迭代：以p_r2作为计算及p_r2，若再次判断重复迭代直至获取最终的p_r2，并计算p_r2对应的C₂；

计算任意生产点的地层压力p_ri。

本发明的技术效果在于：

本发明相较于实测法和物质物质平衡法，降低了关井和测试带来的经济损失；本发明相较于流动物质平衡法避免了近似计算，结果更加符合实际；本发明仅需两个生产点的相关参数，计算过程更加简单。

附图说明

图1为本发明具体实验例的偏差系数Z与压力p的关系图。

具体实施方式

一种基于短期生产数据的地层压力计算方法，包括如下步骤：

1)获取气井的原始地层压力p_i和原始条件下的偏差系数Z_i；获取气井处于拟稳态阶段且时间间隔小于1个月的任意两个生产点的井底流压(p_wf1，p_wf2)和累计产气量(G_p1，G_p2)；获取p≥13MPa条件下的偏差系数Z与压力p的关系式Z＝ap+b，得到斜率a和截距b；

2)计算第一个生产点对应的相关参数：

p_wf1和p_wf2对应的偏差系数Z_wf1和Z_wf2的平均值

地层压力p_r1的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r1：

设定相对误差ε≤1％，判断p_r1是否满足满足则采用下述公式计算p_r1对应的C₁：

C₁＝a(p_r1-p_wf1)；

3)计算第二个生产点对应的相关参数：

地层压力p_r2的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r2：

设定相对误差ε≤1％，判断p_r2是否满足满足则采用下述公式计算p_r2对应的C₂；

C₂＝a(p_r2-p_wf2)；

4)采用下述公式计算任意生产点的地层压力p_ri：

其中，C为C₁与C₂的平均值，G_pi为任意生产点的累计产气量。

优选地，若p_r1不满足即时，则进行迭代：以p_r1作为计算及p_r1，若再次判断重复迭代直至获取最终的p_r1，并计算p_r1对应的C₁；

若p_r2不满足即时，则进行迭代：以p_r2作为计算及p_r2，若再次判断重复迭代直至获取最终的p_r2，并计算p_r2对应的C₂；

计算任意生产点的地层压力p_ri。

具体实施例

实施过程如下：

以鄂尔多斯盆地东南部S4井为例，先进行稳产压降测试，然后进行为期三个月的关井压力恢复测试，选取稳产压降测试阶段的两个生产点进行计算，其中第二个生产点为稳产压降测试阶段的终点，第一个生产点与第二个生产点的时间间隔为21天。

实施过程如下：

1)获取S4井的：原始地层压力p_i＝26.315MPa；原始条件下的偏差系数Z_i＝0.906；第一个生产点的井底流压p_wf1＝15.757MPa，第二个生产点的井底流压p_wf2＝15.505MPa；第一个生产点的累计产气量G_p1＝1154.43×10⁴m³，第二个生产点的累计产气量G_p2＝1218.90×10⁴m³；p≥13MPa条件下的偏差系数Z与压力p的关系式Z＝0.0014p+0.8828，得到a＝0.0014，b＝0.8828；

2)计算第一个生产点对应的相关参数：

p_wf1和p_wf2对应的偏差系数Z_wf1和Z_wf2的平均值：

地层压力p_r1的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r1：

设定相对误差ε＝0.1％，不满足的条件，进行第一次迭代；以p_r1＝21.8846MPa作为第一次迭代近似值计算得到p_r1＝21.9303MPa；再次判断，进行第二次迭代；将p_r1＝21.9303MPa作为第二次迭代近似值计算得到p_r1＝21.9322MPa；再次判断，则认为21.9322MPa为最终的p_r1值；

计算p_r1对应的C₁＝0.0086；

3)计算第二个生产点对应的相关参数：

地层压力p_r2的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r2：

设定相对误差ε＝0.1％，不满足的条件，进行第一次迭代；以p_r2＝21.6242MPa作为第一次迭代近似值计算得到p_r2＝21.6696MPa；再次判断，进行第二次迭代；将p_r2＝21.6696MPa作为第二次迭代近似值计算得到p_r2＝21.6714MPa；再次判断，则认为21.6714MPa为最终的p_r2值；

计算p_r2对应的C₂＝0.0086；

4)计算C＝(C₁+C₂)/2＝0.0086；

5)得到S4井任意生产点的地层压力p_ri计算公式：

结果验证：

根据本发明提供的方法计算得到第二个生产点的地层压力p_r2为21.6714MPa，通过实测法得到第二个生产点的地层压力p_r2为22.5331MPa，相对误差绝对值为3.82％，满足工程计算要求。

Claims (2)

1.一种基于短期生产数据的地层压力计算方法，其特征在于：包括如下步骤：

1)获取气井的原始地层压力p_i和原始条件下的偏差系数Z_i；获取气井处于拟稳态阶段且时间间隔小于1个月的任意两个生产点的井底流压(p_wf1，p_wf2)和累计产气量(G_p1，G_p2)；获取p≥13MPa条件下的偏差系数Z与压力p的关系式Z＝ap+b，得到斜率a和截距b；

2)计算第一个生产点对应的相关参数：

p_wf1和p_wf2对应的偏差系数Z_wf1和Z_wf2的平均值

地层压力p_r1的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r1：

设定相对误差ε≤1％，判断p_r1是否满足满足则采用下述公式计算p_r1对应的C₁：

C₁＝a(p_r1-p_wf1)；

3)计算第二个生产点对应的相关参数：

地层压力p_r2的初始近似值

对应的常数值

对应的偏差系数

采用下述公式计算p_r2：

设定相对误差ε≤1％，判断p_r2是否满足满足则采用下述公式计算p_r2对应的C₂：

C₂＝a(p_r2-p_wf2)；

4)采用下述公式计算任意生产点的地层压力p_ri：

其中，C为C₁与C₂的平均值，G_pi为任意生产点的累计产气量。

2.根据权利要求1所述一种基于短期生产数据的地层压力计算方法，其特征在于：

所述步骤2)中，若判断则进行迭代：以p_r1作为计算及p_r1，若再次判断重复迭代直至

所述步骤3)中，若判断则进行迭代：以p_r2作为计算及p_r2，若再次判断重复迭代直至