一种基于广义流度的油气藏复杂井型井控动半径计算方法技术

本发明专利技术涉及一种基于广义流度的油气藏复杂井型井控动半径计算方法，未投产井进行预测，首先基于广义流度思想，建立求解产能方程或通过现场实测确定产能方程，计算绘制流入动态曲线，建立井筒管流方程计算绘制流出动态曲线，计算预测单井合理产能；建立单井压力不稳定流动模型，并求解得到不稳定流动井底压力解；结合注采过程中对井口压力的约束条件

【技术实现步骤摘要】
一种基于广义流度的油气藏复杂井型井控动半径计算方法


[0001]本专利技术涉及油气藏工程
，具体涉及一种基于广义流度并考虑地层压力约束及井口压力约束的复杂井型井控动半径计算方法
。

技术介绍

[0002]在油气藏开发设计和调整工作中，井控半径的研究一直是受到重视的课题，其关系着井网部署及合理井网密度确定以及井控制范围及波及范围确定，指导井网设计调整
、
生产措施实施及井网控制储量评价等
。
[0003]目前在井控半径研究方面，有多种不同的定义以及计算方法，其一为探测半径法，主要是计算压力脉冲波传播的距离，对此，李传亮等人
(2002)
探讨了探测半径的计算公式，王新海等人
(2008)
给出了非达西渗流
、
非牛顿流体等几种特殊的探测半径计算方法，崔迪生等人
(2005)
推导了径向复合油藏的探测半径计算方法，但该种方法主要为压力传播距离，实际与井控半径相差较大，压力虽然传播到了但实际无法控制到远处的储量，探测半径的计算值通常偏大；其二为生产数据分析方法，主要是通过生产流量及压力数据，结合现代产量递减等分析理论，计算动用范围及动储量等参数，廖伟等人
(2019)
基于产量不稳定分析方法评价了气井动用储量及井控半径，但该方法需要丰富的日生产数据及井底压力数据，对未开采井及无法获取井底压力数据的井来说无法应用分析；其三为临界压力梯度法，主要针对低渗透油气藏存在的启动压力梯度，地层中压力梯度小于启动压力梯度时，外围储层不渗流，从而获得井控半径，傅春梅等人
(2009)
考虑低渗气藏的启动压力梯度研究了气井有效供气半径，此类方法存在一定局限性，启动压力梯度值的确定难度较大；其四为稳定流法，通过稳定流动达西方程反推井控半径，但仅限于常规直井均质模型且仅适用于稳定流动条件，不适用于不稳定流动和复杂井型及储层模型；同时还有其它不同的方法，苟宏刚
(2005)
提出了考虑流量比的影响半径计算方法，但其流量标准难以确定；唐立根
(2018)
结合不稳定渗流方程及物质平衡方程建立了一套储气库井控半径计算方法，方法针对于储气库且其中压力流量等参数需人为确定，难度较大
。
[0004]在实际应用过程中，井控动半径与生产流量
、
地层压力
、
井底压力
、
生产时间等众多参数相关，尤其是压力参数对井控动半径影响较大，现有井控半径方法均集中于某一种参数影响，考虑因素不全面，因此需进一步研究一种考虑地层压力
、
井底压力等参数约束以及考虑到合理产能等其它因素的井控动半径的计算方法
。

技术实现思路

[0005]本专利技术旨在针对上述问题，提出一种基于广义流度的油气藏复杂井型井控动半径计算方法
。
[0006]本专利技术的技术方案在于：一种基于广义流度的油气藏复杂井型井控动半径计算方法，其特征在于：方法如下：
S1
：针对已有开发数据井，通过初始地层压力
、
目前地层压力测试及生产和注入过程中的流量数据，结合物质平衡方程计算已有开发数据井的井控动半径；
S2
：针对未投产井进行预测，首先基于广义流度思想，建立求解产能方程或通过现场实测确定产能方程，计算绘制流入动态曲线，同时建立井筒管流方程计算绘制流出动态曲线，根据节点分析法结合流入动态曲线及流出动态曲线计算预测单井合理产能；
S3
：基于广义流度，建立单井压力不稳定流动模型，并求解得到不稳定流动井底压力解；
S4
：结合注采过程中对井口压力的约束条件
、
预测单井合理产能及单井压力不稳定流动模型计算井底压力，从而迭代计算确定不同时刻合理注采流量；
S5
：结合地层压力约束条件及物质平衡方程迭代计算最优井控动半径
。
[0007]进一步的，步骤
S1
的具体过程为：根据物质平衡方程反推井控动半径方程，根据井控动半径方程代入参数直接计算得到井控动半径，井控动半径方程为：气藏井控动半径方程为：水井井控动半径方程为：式中：
r
e
‑
井控动半径，
m
；
N
f
‑
累计油流量
(
采出为正，注入为负
)
，
m3；
B
o
‑
油体积系数，
m3/m3；
φ
‑
孔隙度，小数；
h
‑
有效厚度，
m
；
C
t
‑
综合压缩系数，
1/Pa
；
B
oi
‑
初始条件下的原油体积系数，
m3/m3；
p
i
‑
原始地层压力，
Pa
；
p
‑
压力，
Pa
；
G
f
‑
累计气流量
(
采出为正，注入为负
)
，
m3；
Z
i
‑
原始地层压力对应的天然气偏差因子，无量纲；
Z
‑
天然气偏差因子，无量纲；
W
f
‑
累计水流量
(
采出为正，注入为负
)
，
m3；
B
w
‑
水体积系数，
m3/m3；
B
wi
‑
初始条件下的水体积系数，
m3/m3。
[0008]方程中的累计流量下标符号统一采用
f
，既能表示采出也能表示注入，也可以按业界常用的采用下标
p
表示采出，下标
i
表示注入，此时方程中的地层压力差
p
i
‑
p
需取绝对值
。
[0009]进一步的，步骤
S2
的具体过程为：
S21
：基于广义流动思想表征流体的流动规律，对任意一种流体运动方程，将其写作
v
＝
﹣
λ
▽
p
的广义流度模型形式；其中，
λ
为广义流度，
m2/(Pa
·
s)
；
v
表示流体流速，
m/s
；
▽
p
表示压力梯度，
Pa/m
；将流体的流动规律采用如上广义流度模型形式进行表征；
S22
：通过产能方程计算确定地层压力与井底压力及流量的关系，进而绘制井底压力与流量的流入动态曲线；
S23本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种基于广义流度的油气藏复杂井型井控动半径计算方法，其特征在于：方法如下：
S1
：针对已有开发数据井，通过初始地层压力
、
目前地层压力测试及生产和注入过程中的流量数据，结合物质平衡方程计算已有开发数据井的井控动半径；
S2
：针对未投产井进行预测，首先基于广义流度思想，建立求解产能方程或通过现场实测确定产能方程，计算绘制流入动态曲线，同时建立井筒管流方程计算绘制流出动态曲线，根据节点分析法结合流入动态曲线及流出动态曲线计算预测单井合理产能；
S3
：基于广义流度，建立单井压力不稳定流动模型，并求解得到不稳定流动井底压力解；
S4
：结合注采过程中对井口压力的约束条件
、
预测单井合理产能及单井压力不稳定流动模型计算井底压力，从而迭代计算确定不同时刻合理注采流量；
S5
：结合地层压力约束条件及物质平衡方程迭代计算最优井控动半径
。2.
根据权利要求1所述基于广义流度的油气藏复杂井型井控动半径计算方法，其特征在于：所述
S1
的具体过程为：根据物质平衡方程反推井控动半径方程，根据井控动半径方程代入参数直接计算得到井控动半径，油井井控动半径方程为：气井井控动半径方程为：水井井控动半径方程为：式中：
r
e
‑
井控动半径，
m
；
N
f
‑
累计油流量，
m3；
B
o
‑
油体积系数，
m3/m3；
φ
‑
孔隙度，小数；
h
‑
有效厚度，
m
；
C
t
‑
综合压缩系数，
1/Pa
；
B
oi
‑
初始条件下的原油体积系数，
m3/m3；
p
i
‑
原始地层压力，
Pa
；
p
‑
压力，
Pa
；
G
f
‑
累计气流量，
m3；
Z
i
‑
原始地层压力对应的天然气偏差因子，无量纲；
Z
‑
天然气偏差因子，无量纲；
W
f
‑
累计水流量，
m3；
B
w
‑
水体积系数，
m3/m3；
B
wi
‑
初始条件下的水体积系数，
m3/m3。3.
根据权利要求2所述基于广义流度的油气藏复杂井型井控动半径计算方法，其特征在于：所述
S2
的具体过程为：
S21
：基于广义流动思想表征流体的流动规律，对任意一种流体运动方程，将其写作
v
＝
﹣
λ
▽
p
的广义流度模型形式；其中，
λ
为广义流度，
m2/(Pa
·
s)
；
v
表示流体流...

【专利技术属性】
技术研发人员：林加恩，李尧，
申请(专利权)人：西安华线石油科技有限公司，
类型：发明
国别省市：