一种基于流型转换界限的水平气井持液率计算方法技术

本发明专利技术提供一种基于流型转换界限的水平气井持液率计算方法，包括：根据持液率随气流速变化趋势构建三参数方程式，将计算出的一定井筒液流速时常压/常温条件下垂直管中环状流

【技术实现步骤摘要】
一种基于流型转换界限的水平气井持液率计算方法


[0001]本专利技术属于气井排水采气
，具体涉及一种基于流型转换界限的水平气井持液率计算方法。

技术介绍

[0002]气藏开发到中后期产水明显上升，及时将液体携带出井筒是维持气井正常生产、延长生命周期的重要举措。由持液率主导的重力降是井筒压力梯度的重要组成部分，同时也是气井连续携液设计中的关键参数。因此，准确预测持液率是明确气井井筒压降、实现气井连续携液的关键。
[0003]目前，关于持液率的计算方法众多，有机理模型法、经验公式法、实验测试法等。其中，经验或半经验模型通常以实验测试为基础，仅适用于实验参数范围内的情况，计算时需要先进行模型评价优选。机理模型则通过流型识别后选择相应的计算方法，然而不同学者建立的流型判别标准各异，导致计算结果精度差异大，工程适用性差。
[0004]此外，上述计算方法多基于垂直井建立，而水平井分为垂直段、倾斜段及水平段，井斜角由90
°
到0
°
连续变化，流型分布复杂多变，受井斜角、管径变化等因素影响更明显，目前针对水平井一般采用描述水平管、倾斜管和垂直管的3个流型图来分段处理，各流型图实验条件差异大，适应范围也不同，尤其是产水量极小的气井，气液比极高，易超出工程常用气液两相管流流型图的坐标值范围，导致预测水平气井气液两相流型误差大。由此可见，目前缺乏简便实用且针对水平井全井筒的统一持液率计算模型。
[0005]本专利技术提供一种根据气相表观流速计算水平井持液率的三参数模型，无需大量的模型比对优选及流型判断计算，适用于工程计算，为简便、准确预测水平井全井筒持液率及井筒压降提供思路。

技术实现思路

[0006]本专利技术针对现有计算方法的不足，基于流型转换界限，通过构建持液率随气流速变化的三参数式方程提供一种水平气井持液率计算方法。
[0007]为实现上述目的，本专利技术提供如下技术方案：
[0008]步骤一：获取气井井筒某一位置的压力、温度、倾斜角度、气流速、液流速、管径、气体相对密度等参数；
[0009]步骤二：构建持液率随气流速变化的三参数式方程；
[0010]H
L
＝Av
SGB
+C
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0011]式中，H
L
为垂直段持液率，无量纲；v
SG
为表观气流速，m/s；A、B、C均为模型系数，无量纲。
[0012]步骤三：计算一定气井井筒液流速时常压/常温条件下(0.1MPa，20℃)垂直管中环状流
‑
搅动流转变点、搅动流
‑
段塞流转变点和段塞流
‑
泡状流转变点的持液率及对应气流速值；
[0013](1)环状流
‑
搅动流转变点对应气流速：
[0014][0015]式中，v
A
‑
C
——环状流
‑
搅动流转变点对应气流速，m/s；f
i
——内摩擦因子，无量纲；ρ
G,标况
、ρ
L,标况
——标况(0.1MPa，20℃)下气体、液体密度，kg/m3；g——重力加速度，9.8；r0——油管半径，m；δ——液膜厚度，m。
[0016]采用Fore模型计算内摩擦因子：
[0017][0018]式中，Re
G
——气体雷诺数，无量纲。
[0019]液膜厚度计算公式：
[0020][0021]式中Q
L
——液流量，m3/s；μ
L
——液体粘度，mPa
·
s。
[0022]环状流
‑
搅动流转变点对应持液率：
[0023][0024]式中，H
A
‑
C
——环状流
‑
搅动流转变点对应持液率，无量纲；D——油管直径，m。
[0025]那么可以得到环状流
‑
搅动流转变点持液率与气流速的关系式：
[0026]H
A
‑
C
＝Av
A
‑
CB
+C
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0027](2)搅动流
‑
段塞流转变点、段塞流
‑
泡状流转变点对应持液率分别为
[0028]H
C
‑
S
＝0.22
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0029]H
S
‑
B
＝0.75
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0030]式中，H
C
‑
S
、H
S
‑
B
分别为搅动流
‑
段塞流转变点和段塞流
‑
泡状流转变点对应持液率，无量纲。
[0031]利用Hasan
‑
Kabir模型计算搅动流
‑
段塞流转变点气流速：
[0032]v
C
‑
S
＝[0.0051(ρ
L,标况
v
SL2
)
1.7
/ρ
G,标况
]0.5
(当ρ
L,标况
v
SL2
＜74.4时)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0033]或
[0034]v
C
‑
S
＝[(25.41gρ
L,标况
v
SL2
‑
38.9)/ρ
G,标况
]0.5
(当ρ
L,标况
v
SL2
＞74.4时)
ꢀꢀꢀꢀ
(10)
[0035]式中，v
C
‑
S
——搅动流
‑
段塞流转变点对应气流速，m/s；v
SL
——表观液流速，m/s；σ——气水界面张力，N/m；
[0036]采用漂移模型计算段塞流
‑
泡状流转变点气流速：
[0037][0038]式中，v
S
‑
B
——段塞流
‑
泡状流转变点对应气流速，m/s；v
SL
——表观液流速，m/s；
σ——气水界面张力，N/m；
[0039]那么可以得到搅动流
‑
段塞流转变点和段塞流
‑
泡状流转变点持液率与气流速的关系式：
[0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于流型转换界限的水平气井持液率计算方法，其特征在于，所述方法主要包括以下步骤：步骤一：获取气井井筒某一位置的压力、温度、倾斜角度、气流速、液流速、管径、气体相对密度等参数；步骤二：构建持液率随气流速变化的三参数式方程；步骤三：计算一定气井井筒液流速时常压/常温条件下(0.1MPa，20℃)垂直管中环状流
‑
搅动流转变点、搅动流
‑
段塞流转变点和段塞流
‑
泡状流转变点的持液率及对应气流速值；步骤四：计算三参数式方程的系数A、B、C；步骤五：采用角度修正关系式对上述垂直段持液率计算结果进行角度修正得到该倾斜角的持液率计算新模型；步骤六：基于流动相似准则利用无因次准数将高压下气流速转化为常压下对应气流速；步骤七：将转化后的常压下对应气流速v
SG,标况
代入步骤五所得的新模型即可得到一定压力条件下水平气井持液率。2.根据权利要求1所述的一种基于流型转换界限的水平气井持液率计算方法，其特征在于：步骤二构建持液率随气流速变化的三参数式方程具体为：H
L
＝Av
SGB
+C
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式中，H
L
为垂直段持液率，无量纲；v
SG
为表观气流速，m/s；A、B、C均为模型系数，无量纲。3.根据权利要求1所述的一种基于流型转换界限的水平气井持液率计算方法，其特征在于：步骤四中对模型系数B进行迭代求解，具体如下：将环状流
‑
搅动流转变点、搅动流
‑
段塞流转变点和段塞流
‑
泡状流转变点持液率与气流速的关系式进行组合消去系数A和C可得：对式(2)进行变形...

【专利技术属性】
技术研发人员：罗程程，高丽蓉，刘永辉，杨建英，王强，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：