CN113343607B - 一种基于管壁液膜运动规律的携液临界气流速计算方法 - Google Patents
一种基于管壁液膜运动规律的携液临界气流速计算方法 Download PDFInfo
- Publication number
- CN113343607B CN113343607B CN202110748685.8A CN202110748685A CN113343607B CN 113343607 B CN113343607 B CN 113343607B CN 202110748685 A CN202110748685 A CN 202110748685A CN 113343607 B CN113343607 B CN 113343607B
- Authority
- CN
- China
- Prior art keywords
- liquid
- liquid film
- calculating
- carrying
- flow rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于管壁液膜运动规律的携液临界气流速计算方法,包括:定义气井积液临界点;建立液膜沿管壁的一维速度剖面方程;计算向下流动液膜厚度;计算液膜总厚度;计算层流层液膜厚度;计算气芯拖曳力;计算携液临界气流速。本发明优点在于:充分考虑了液膜的流动规律,定义了气井积液新的判断标准,基于液膜力平衡分析推导得到新的携液临界气流速模型。本发明方法适用范围更广,也具有更准确的预测结果。
Description
技术领域
本发明涉及油气田排水采气技术领域,特别涉及一种基于管壁液膜运动规律的携液临界气流速计算方法。
背景技术
天然气作为一种安全可靠、利用率高、应用广泛的清洁能源,在国家的能源战略地位与日俱增。而国内气井开采中后期占比高,导致气井积液问题频发,严重影响气井正常生产。准确预测气井积液并及时实施排水采气工艺措施是气井稳产的关键。气井积液最直接的判断方法之一是进行压力梯度测试,判断气井井筒中是否存在压力梯度拐点,但并不是所以气井都有条件进行压力测试。因此,建立符合现场实际情况的气井积液预测理论模型,提前判断气井是否积液并及时介入排水采气工艺措施是实现气井稳产的一种重要手段。
工程中常用预测气井积液的模型可分为液滴反转模型和液膜反转模型。其中,液滴反转模型是基于井筒中单个液滴力平衡分析所建,但大量实验和研究表明,井筒中液滴数量少且只有极少量液滴发生反转,说明液膜的向下流动是造成气井积液的主要因素。液膜反转模型则是假定液膜反转时是气井积液的起始点,但液膜反转时对应气量非常大,大量工程实践也表明液膜反转模型计算值远远大于真实气井积液时对应气量。因此,目前缺乏有效能表征气井积液的准确模型。
发明内容
本发明针对现有技术的缺陷,提供了一种基于管壁液膜运动规律的携液临界气流速计算方法。
为了实现以上发明目的,本发明采取的技术方案如下:
一种基于管壁液膜运动规律的携液临界气流速计算方法,包括以下步骤:
步骤一:定义气井积液临界点;定义气井积液临界点为向下液膜流量为总产液量的一半,且向上净流量全部以波动形式向上携带。将液膜层划分为层流层及波动层。层流层的液膜总流量为零,即层流层液膜向下流动流量=层流层液膜向上流动流量=1/2QL,液膜以波动形式向上流量为QL,向上流动液膜总流量则为3/2QL。
步骤二:建立液膜沿管壁的一维速度剖面方程;
步骤三:计算向下流动液膜厚度;
步骤四:计算液膜总厚度;
步骤五:计算层流层液膜厚度;
步骤六:计算气芯拖曳力;
步骤七:计算携液临界气流速。
进一步地,步骤二中建立液膜沿管壁的一维速度剖面方程为:液膜沿管壁流动时遵循Navier-Stokes方程:
式中,μL为液体粘度,Pa·s;ρL为液体密度,kg/m3;v为液膜沿轴向速度分布,m/s;g为重力加速度,m/s2;p为压力,Pa。
忽略压力项并对方程(1)进行积分,分别得到关于速度变化率及速度的表达:
当气井处于积液临界点时,有边界条件如下:
v|x=0=0 (5)
式中,δ0为向下流动液膜厚度,m。
分别将式(4)和(5)代入式(2)、(3),可得速度剖面方程:
进一步地,步骤三中计算向下流动液膜厚度具体为:当气井处于积液临界点时,对向下流动液膜进行速度积分:
由式(7)可得:
进一步地,步骤四中计算液膜总厚度具体为:当气井处于积液临界点时,对向上流动液膜进行速度积分:
由式(9)可得:
Aδ3+Bδ2+C=0 (10)
C=QL (13)
式中A,B,C为系数。具体求解步骤为:
(1)根据给定位置处液体密度ρL、液体粘度μL和计算的向下流动液膜厚度δ0计算系数A、B、C;
(3)利用下式求解得到式(6)中一元三次方程的实根,即液膜总厚度:
进一步地,步骤五中计算层流层液膜厚度具体为:当气井处于积液临界点时,对层流层液膜进行速度积分:
式中,δ'为层流层液膜厚度,m。
由式(16)可得:
进一步地,步骤六中计算气芯拖曳力具体为:当气井处于积液临界点时,波动层液膜以波动形式向上携带,不由气芯提供拖曳力。层流层液膜向上流动由气芯提供拖曳力,力平衡关系为:
τi=ρLgδ0 (18)
式中,τi——内切应力,Pa。
进一步地,步骤七中计算携液临界气流速具体为:采用Wallis模型计算携液临界气流速:
式中,vSG为携液临界气流速,m/s;fi为内摩擦因子,无量纲;ρG为气体密度,kg/m3;ro为油管半径,m。
与现有技术相比,本发明的优点在于:
充分考虑了液膜的流动规律,定义了气井积液新的判断标准,进而基于液膜力平衡分析推导得到新的携液临界气流速模型。因此,该方法适用范围更广,也具有更准确的预测结果。
附图说明
图1是本发明实施例一维液膜速度剖面示意图;
图2是本发明实施例液膜波动部分示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
气井的生产总会有液体产出,在建立气井积液模型过程中,液膜流动规律应将作为主要的影响因素。假定QL为气井稳定生产的产液量且全部以液膜形式流动,如图1所示。
步骤一:定义气井积液临界点为向下液膜流量为总产液量的一半,且向上净流量全部以波动形式向上携带。将液膜层划分为层流层及波动层,液膜厚度如图2所示。层流层的液膜总流量为零,即层流层液膜向下流动流量=层流层液膜向上流动流量=1/2QL,液膜以波动形式向上流量为QL,向上流动液膜总流量则为3/2QL。
步骤二:建立液膜沿管壁的一维速度剖面方程。液膜沿管壁流动时遵循Navier-Stokes方程:
式中,μL为液体粘度,Pa·s;ρL为液体密度,kg/m3;v为液膜沿轴向速度分布,m/s;g为重力加速度,m/s2;p为压力,Pa。
忽略压力项并对方程(1)进行积分,分别得到关于速度变化率及速度的表达式:
当气井处于积液临界点时,有边界条件如下:
v|x=0=0 (5)
式中,δ0为向下流动液膜厚度,m。
分别将式(4)和(5)代入式(2)、(3),可得速度剖面方程:
步骤三:计算向下流动液膜厚度。当气井处于积液临界点时,对向下流动液膜进行速度积分:
由式(7)可得:
步骤四:计算液膜总厚度。当气井处于积液临界点时,对向上流动液膜进行速度积分:
由式(9)可得:
Aδ3+Bδ2+C=0 (10)
C=QL (13)
式中A,B,C为系数。具体求解步骤为:
(1)根据给定位置处液体密度ρL、液体粘度μL和计算的向下流动液膜厚度δ0计算系数A、B、C;
(3)利用下式求解得到式(6)中一元三次方程的实根,即液膜总厚度:
步骤五:计算层流层液膜厚度。当气井处于积液临界点时,对层流层液膜进行速度积分:
式中,δ'为层流层液膜厚度,m。
由式(16)可得:
步骤六:计算气芯拖曳力。当气井处于积液临界点时,波动层液膜以波动形式向上携带,不由气芯提供拖曳力。层流层液膜向上流动由气芯提供拖曳力,力平衡关系为:
τi=ρLgδ0 (18)
式中,τi为内切应力,Pa。
步骤七:计算携液临界气流速。采用Wallis模型计算携液临界气流速:
式中,vSG为携液临界气流速,m/s;fi为内摩擦因子,无量纲;ρG为气体密度,kg/m3;ro为油管半径,m。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (7)
1.一种基于管壁液膜运动规律的携液临界气流速计算方法,其特征在于,包括以下步骤:
步骤一:定义气井积液临界点;定义气井积液临界点为向下液膜流量为总产液量的一半,且向上净流量全部以波动形式向上携带;将液膜层划分为层流层及波动层;层流层的液膜总流量为零,即层流层液膜向下流动流量=层流层液膜向上流动流量=1/2QL,液膜以波动形式向上流量为QL,向上流动液膜总流量则为3/2QL;
步骤二:建立液膜沿管壁的一维速度剖面方程;
步骤三:计算向下流动液膜厚度;
步骤四:计算液膜总厚度;
步骤五:计算层流层液膜厚度;
步骤六:计算气芯拖曳力;
步骤七:计算携液临界气流速。
6.根据权利要求1所述的携液临界气流速计算方法,其特征在于:步骤六中计算气芯拖曳力具体为:当气井处于积液临界点时,波动层液膜以波动形式向上携带,不由气芯提供拖曳力;层流层液膜向上流动由气芯提供拖曳力,力平衡关系为:
τi=ρLgδ0 (18)
式中,τi——内切应力,Pa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110748685.8A CN113343607B (zh) | 2021-07-02 | 2021-07-02 | 一种基于管壁液膜运动规律的携液临界气流速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110748685.8A CN113343607B (zh) | 2021-07-02 | 2021-07-02 | 一种基于管壁液膜运动规律的携液临界气流速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113343607A CN113343607A (zh) | 2021-09-03 |
CN113343607B true CN113343607B (zh) | 2022-04-01 |
Family
ID=77482261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110748685.8A Active CN113343607B (zh) | 2021-07-02 | 2021-07-02 | 一种基于管壁液膜运动规律的携液临界气流速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113343607B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110878688A (zh) * | 2018-08-28 | 2020-03-13 | 中国石油化工股份有限公司 | 一种高气液比水平井临界携液流量的计算方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106401570B (zh) * | 2015-07-30 | 2019-05-07 | 中国石油化工股份有限公司 | 页岩气井产水的确定方法、积液的确定方法及排液方法 |
CN110866321B (zh) * | 2018-08-14 | 2022-05-31 | 中国石油化工股份有限公司 | 一种适用于油气水同产水平井积液判别方法 |
NO345508B1 (en) * | 2019-05-24 | 2021-03-15 | Ledaflow Tech Da | Method and tool for planning and dimensioning subsea pipelines for produced fluids |
-
2021
- 2021-07-02 CN CN202110748685.8A patent/CN113343607B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110878688A (zh) * | 2018-08-28 | 2020-03-13 | 中国石油化工股份有限公司 | 一种高气液比水平井临界携液流量的计算方法 |
Non-Patent Citations (1)
Title |
---|
基于紊流条件下的气井临界携液流量计算模型;明瑞卿等;《地质科技情报》;20180515(第03期);第254-258页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113343607A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110334431B (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
Peng et al. | Analysis of particle deposition in a new‐type rectifying plate system during shale gas extraction | |
CN106600442B (zh) | 一种深水气井测试临界携液条件优化方法 | |
CN108266176B (zh) | 一种基于井筒模型的天然气井口流量计算方法 | |
CN104847314B (zh) | 高温高压油气直井单相流射孔完井参数优化方法 | |
CN107977513B (zh) | 一种基于路径搜索的天然气动态流动温度预测方法 | |
Luan et al. | A new model for the accurate prediction of liquid loading in low-pressure gas wells | |
CN110344818A (zh) | 柱塞气举井柱塞上行阶段液体段塞界面跟踪方法 | |
Chappell et al. | Aerodynamic performance of suction-side gill region film cooling | |
CN114819347B (zh) | 多相内流激励下的输流管道振动响应预测方法 | |
CN114580100B (zh) | 压裂水平井全井筒压力计算方法、设备和计算机可读储存介质 | |
CN113779901B (zh) | 一种多井统一计量下的单井动态储量计算方法 | |
CN113343607B (zh) | 一种基于管壁液膜运动规律的携液临界气流速计算方法 | |
Zhao et al. | A case study of hydraulic collecting a single spherical particle | |
CN111553096B (zh) | 一种可变喷油率的柴油喷雾贯穿距预测方法 | |
CN108090321B (zh) | 一种水平气井倾斜管液膜携带临界气流量计算方法 | |
CN116838327A (zh) | 一种气水油三相流井筒压降预测方法 | |
CN105184061A (zh) | 产气井温度压力分布的数值模拟方法 | |
CN112800588B (zh) | 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 | |
CN207990210U (zh) | 天然气地下储气库单井注采同管双向计量装置 | |
Zahedi et al. | Liquid film thickness prediction in elbows for annular flows | |
Tang et al. | Effect of air injector on the airlift performance in air-water-solid three-phase flow | |
CN113969779B (zh) | 一种注气井井筒压力分布的确定方法 | |
CN116856908B (zh) | 一种页岩气井携砂临界流速的实验确定方法 | |
CN114580611B (zh) | 一种基于pinn的制冷剂多相流充注流量获取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |