CN105387895B - 基于流动密度的长喉颈文丘里管湿气流量测量方法 - Google Patents
基于流动密度的长喉颈文丘里管湿气流量测量方法 Download PDFInfo
- Publication number
- CN105387895B CN105387895B CN201510674413.2A CN201510674413A CN105387895B CN 105387895 B CN105387895 B CN 105387895B CN 201510674413 A CN201510674413 A CN 201510674413A CN 105387895 B CN105387895 B CN 105387895B
- Authority
- CN
- China
- Prior art keywords
- formula
- calculated
- flow
- calculate
- gas
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
- G01F1/34—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure
- G01F1/36—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure the pressure or differential pressure being created by the use of flow constriction
Landscapes
- Measuring Volume Flow (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
Abstract
本发明涉及一种基于流动密度的长喉颈文丘里管湿气流量测量方法:计算气相流量初值;计算气相流速;计算气相Froude数:计算流出系数;计算L‑M参数;计算干度初值;计算截面含气率;求虚高;计算出气相单独流过文丘里管时的差压;计算得到两相混合质量流量;计算两相流的流动密度;计算出体积含气率;计算出新的质量含气率和新的截面含气率,判断是否收敛,如果不收敛,将计算出的截面含气率作为初值重新计算;计算虚高;判断是否收敛,如果不收敛,将计算出的虚高作为下一次迭代的初值重新执行1~13步的计算,直至收敛为止;得到湿气的气相和液相流量。本发明可以消除两相介质的速度差对测量结果的影响。
Description
技术领域
本发明属于两相流量测量技术领域,涉及一种湿气流量测量方法。
背景技术
湿气通常是指气相为连续相,液相为离散相的气液两相流。美国机械工程师学会ASME将其界定为Lockhart-Martinelli参数(简称L-M参数)X小于0.3的气液两相流。湿气普遍存在于自然界与工业现场中,如常规天然气田井口产出气(湿天然气)、煤层气、湿饱和蒸汽以及页岩气等均属湿气。随着工业生产的发展,对湿气不分离计量的要求越来越高。然而已有研究表明,湿气流量测量问题非常复杂,检测难度大,利用传统的单相气体流量计直接测量湿气流量,会因液相的影响而产生一定的误差,甚至无法正常工作。因此湿气测量在相当长的一段时间内都将是流量测量领域的一个重要研究方向和研究热点。
文献1提出的“基于截面含气率的文丘里管湿气虚高双迭代模型”[1],利用节流原件测量原理得到两相介质的流动密度。截面含气率指的是某一瞬间气相所占截面与总流道截面之比[2],此定义为瞬时量。而在实际的工程应用时,主要关心的是两相流的体积流量,即流经流量计气相和液相一段时间内的累积体积流量。由于在实际流动中气液两相之间是存在速度差的,截面含气率并不等于体积含气率。
专利2012104654424给出了一种基于长喉颈文丘里的双差压湿气流量测量装置。
参考文献:
[1]Chao Yuan,Ying Xu,Tao Zhang,Xili Ba,Huaxiang Wang,Dual-iterativemodel for gas condensate measurementbased onvoid fraction,Journal ofNaturalGas Science and Engineering,2015,24:330-336.
[2]阎昌琪.哈尔滨工程大学出版社,2009.
发明内容
本发明的目的是基于上述长喉颈文丘里湿气流量测量装置,提供一种可以消除两相介质的速度差对测量结果的影响的湿气流量测量方法。本发明的技术方案如下:
一种基于流动密度的长喉颈文丘里管湿气流量测量方法,该方法利用长喉颈文丘里差压式流量装置进行湿气流量测量,流量计算方法如下:
①设虚高初值OR0为1,流出系数初值C0为1,带入公式(1)得到气相流量初值,并令气相流量Wg=Wg0:
其中:A为管道的截面积,β=d/D为流量计的节流比,D为管道直径,d为喉部直径,ρg为被测气体密度,Δpg为文丘里的前差压,ε为可膨胀系数,
其中:为压力比,即文丘里前端压力P与前差压ΔP1的差值与前端压力P的比值,ΔP1为前差压即文丘里上游和喉部取压位置之间的压力差,κ为等熵指数,取κ=1.4;
②根据公式(2)计算出气相流速:
③根据公式(3)计算出气相Froude数:
其中,ρl为被测气液两相流中液相密度;
④根据公式(4)计算出流出系数,并令C0=C:
其中,L为长喉颈文丘里管的喉部长度;
⑤根据式(5)计算出L-M参数:
其中:ΔP2为后差压即文丘里喉部和下游取压位置之间的压力差;
⑥根据公式(6)计算干度初值,并令质量含气率x=x0:
⑦根据公式(7)计算截面含气率:
其中:αH为均相流条件下的截面含气率;
⑧根据公式(8)求出虚高;
⑨由公式(9)计算出气相单独流过文丘里管时的差压;
⑩根据式(10)计算得到两相混合质量流量:
根据公式(11)计算出两相流的流动密度;
根据公式(12)计算出体积含气率;
根据公式(13)计算出新的质量含气率x,之后根据公式(7)计算出新的截面含气率α,并将新计算出的截面含气率与第7步计算的截面含气率比较,两者的相对误差小于0.1%时,则认为结果已收敛,否则,将计算出的截面含气率作为初值重新执行8~13步的计算,直至收敛为止,将该值带入下一步计算;
根据公式(8)计算出虚高;
将计算出的虚高与上一次虚高迭代值比较,两者的相对误差小于0.1%时,则认为结果已收敛,不再计算;若大于0.1%,将计算出的虚高作为下一次迭代的初值重新执行1~13步的计算,直至收敛为止;
得到湿气的气相和液相流量。
本发明直接利用体积含气率,从而消除两相介质的速度差对测量结果的影响。
附图说明
图1为长喉颈文丘里流量计的结构原理图。
图2迭代结果的相对误差。
具体实施方式
如图1所示,长喉颈文丘里是一种差压式流量计,对于差压式流量计,根据伯努利方程和连续性方程的推导,即可通过测量出的差压值计算单相流的质量流量:
其中:C为差压式流量计的流出系数,D为管道直径,β=d/D为流量计的节流比,d为喉部直径,为被测气体密度,为文丘里的前差压,为可膨胀系数,
其中:τ为压力比,为等熵指数,是在基本可逆绝热(等熵)转换的条件下,压力的相对变化量与密度的相对变化量的比值。会随着气体的性质及其温度和压力的变化而变化,在本发明中取=1.4。
液相的存在使得在应用差压式流量计测量湿气的过程中,测得的差压高于等量气体单独流过流量计时的差压,称为“虚高”。此时湿气中的气相流量的计算公式为
其中:为湿气流过差压式流量计时产生的差压,测得的差压高于等量气体单独流过流量计时的差压。
为了求出真实的气相流量,引入一个量纲为1的修正系数虚高OR,其表达式为
因此欲准确测量湿气两相流的气相流量,则必须准确知道虚高OR的计算模型。本专利提出了一种采用双迭代方法的虚高OR的计算模型,计算过程如下:
1、设虚高初值OR0为1,流出系数初值C0为1,带入公式(1)得到气相流量初值,并令气相流量Wg=Wg0:
其中:A为管道的截面积,β=d/D为流量计的节流比,D为管道直径,d为喉部直径,ρg为被测气体密度,Δpg为文丘里的前差压,ε为可膨胀系数,
其中:为压力比即文丘里前端压力P与前差压ΔP1的差值与前端压力P的比值,ΔP1为前差压即文丘里上游和喉部取压位置之间的压力差,κ为等熵指数,取κ=1.4;
2、根据公式(2)计算出气相流速:
3、根据公式(3)计算出气相Froude数:
其中,ρl为被测气液两相流中液相密度;
4、根据公式(4)计算出流出系数,并令C0=C:
其中,L为长喉颈文丘里管的喉部长度;
5、根据式(5)计算出L-M参数:
其中:ΔP2为后差压即文丘里喉部和下游取压位置之间的压力差。
6、根据公式(6)计算干度初值,并令质量含气率x=x0:
7、根据公式(7)计算截面含气率:
其中:αH为均相流条件下的截面含气率。
8、根据公式(8)求出虚高;
9、由公式(9)计算出气相单独流过文丘里管时的差压;
10、根据式(10)计算得到两相混合质量流量:
11、根据公式(11)计算出两相流的流动密度;
12、根据公式(12)计算出体积含气率;
13、根据公式(13)计算出新的质量含气率x,之后根据公式(7)计算出新的截面含气率α,并将新计算出的截面含气率与第7步计算的截面含气率比较,两者的相对误差小于0.1%时,则认为结果已收敛,否则,将计算出的截面含气率作为初值重新执行8~13步的计算,直至收敛为止,将该值带入下一步计算。
14、根据公式(8)计算出虚高;
15、将计算出的虚高与上一次虚高迭代值比较,两者的相对误差小于0.1%时,则认为结果已收敛,不再计算。若大于0.1%,将计算出的虚高作为下一次迭代的初值重新执行1~13步的计算,直至收敛为止。
通过以上的迭代过程可以得到虚高值,依此推算出湿气中的气相和液相流量,从而实现湿气流量的测量。在测量过程中,所测量的体积流量受文丘里管取压口之间的差压影响,而差压值受取压口之间的所有流体的影响,并非某一截面信息所能表示,因此在迭代过程中,使用流动密度更加符合实际的工程应用,也与测量原理形成统一。综上所述,此专利模型更加合理。
经实验数据验证,图2给出第一、三、五次迭代后的虚高预测相对误差,可明显看出迭代结果的速度和精度都有了明显的提高。
Claims (1)
1.一种基于流动密度的长喉颈文丘里管湿气流量测量方法,该方法利用长喉颈文丘里差压式流量装置进行湿气流量测量,流量计算方法如下:
①设虚高初值OR0为1,流出系数初值C0为1,带入公式(1)得到气相流量初值,并令气相流量Wg=Wg0:
其中:A为管道的截面积,β=d/D为流量计的节流比,D为管道直径,d为喉部直径,ρg为被测气体密度,ΔPg为文丘里的前差压,ε为可膨胀系数,
其中:为压力比,即文丘里前端压力P与前差压ΔP1的差值与前端压力P的比值,ΔP1为前差压即文丘里上游和喉部取压位置之间的压力差,κ为等熵指数,取κ=1.4;
②根据公式(2)计算出气相流速:
③根据公式(3)计算出气相Froude数:
其中,ρl为被测气液两相流中液相密度;
④根据公式(4)计算出流出系数,并令C0=C:
其中,L为长喉颈文丘里管的喉部长度;
⑤根据式(5)计算出L-M参数:
其中:ΔP2为后差压即文丘里喉部和下游取压位置之间的压力差;
⑥根据公式(6)计算干度初值,并令质量含气率x=x0:
⑦根据公式(7)计算截面含气率:
其中:αH为均相流条件下的截面含气率;
⑧根据公式(8)求出虚高;
⑨由公式(9)计算出气相单独流过文丘里管时的差压;
⑩根据式(10)计算得到两相混合质量流量:
根据公式(11)计算出两相流的流动密度;
根据公式(12)计算出体积含气率;
根据公式(13)计算出新的质量含气率x,之后根据公式(7)计算出新的截面含气率α,并将新计算出的截面含气率与第7步计算的截面含气率比较,两者的相对误差小于0.1%时,则认为结果已收敛,否则,将计算出的截面含气率作为初值重新执行8~13步的计算,直至收敛为止,将该值带入下一步计算;
根据公式(8)计算出虚高;
将计算出的虚高与上一次虚高迭代值比较,两者的相对误差小于0.1%时,则认为结果已收敛,不再计算;若大于0.1%,将计算出的虚高作为下一次迭代的初值重新执行1~13步的计算,直至收敛为止;
得到湿气的气相和液相流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510674413.2A CN105387895B (zh) | 2015-10-16 | 2015-10-16 | 基于流动密度的长喉颈文丘里管湿气流量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510674413.2A CN105387895B (zh) | 2015-10-16 | 2015-10-16 | 基于流动密度的长喉颈文丘里管湿气流量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105387895A CN105387895A (zh) | 2016-03-09 |
CN105387895B true CN105387895B (zh) | 2018-06-19 |
Family
ID=55420442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510674413.2A Expired - Fee Related CN105387895B (zh) | 2015-10-16 | 2015-10-16 | 基于流动密度的长喉颈文丘里管湿气流量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105387895B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105675070A (zh) * | 2016-03-28 | 2016-06-15 | 中国海洋石油总公司 | 异型文丘里流量计和用其测量多相流中气液相流量的方法 |
CN106123983A (zh) * | 2016-06-28 | 2016-11-16 | 浙江汽轮成套技术开发有限公司 | 空气压缩机筒状进气管的进气流量测量方法 |
CN106768117B (zh) * | 2017-02-20 | 2019-09-20 | 天津大学 | 一种湿蒸汽干度流量双参数测量方法 |
CN108168612A (zh) * | 2017-12-27 | 2018-06-15 | 中国石油大学(华东) | 基于差压信号波动的气液两相流体积含气率测量方法 |
CN108871511B (zh) * | 2018-06-29 | 2020-06-09 | 南京罕华流体技术有限公司 | 一种联动型工业流量精确计量方法 |
CN109635405B (zh) * | 2018-12-05 | 2020-12-18 | 中国石油大学(北京) | 一种基于空间守恒的多相流动瞬态计算方法及系统 |
CN109945936B (zh) * | 2019-03-25 | 2020-06-09 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种基于节流装置测量的水及蒸汽流量计算方法 |
CN110514257B (zh) * | 2019-08-29 | 2020-08-18 | 西安交通大学 | 一种基于文丘里的低含液率湿气两相流量测量装置及方法 |
CN111155981B (zh) * | 2019-12-20 | 2023-05-09 | 中煤科工集团西安研究院有限公司 | 一种煤层气井双煤层合层排采分层产气量监测方法 |
CN111506122B (zh) * | 2020-04-21 | 2023-09-19 | 中国舰船研究设计中心 | 一种洗消通道压力流量控制方法 |
CN112747796A (zh) * | 2020-12-23 | 2021-05-04 | 一汽解放汽车有限公司 | 一种柴油机进气流量的测量装置及计算方法 |
CN112964316B (zh) * | 2021-02-01 | 2024-03-22 | 深圳市联恒星科技有限公司 | 基于长喉颈文丘里管湿气测量及不确定度预测方法、系统 |
CN113375741A (zh) * | 2021-04-29 | 2021-09-10 | 安徽中控仪表有限公司 | 基于三差压数据拟合模型的湿气两相流量计量装置及方法 |
CN115420342B (zh) * | 2022-11-03 | 2023-03-24 | 海默新宸水下技术(上海)有限公司 | 一种基于含气率拟合的湿天然气计量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101413817A (zh) * | 2008-12-05 | 2009-04-22 | 天津大学 | 双差压节流湿气测量方法 |
CN103048019A (zh) * | 2012-11-15 | 2013-04-17 | 天津大学 | 基于长喉颈文丘里的双差压湿气流量测量装置 |
CN202994216U (zh) * | 2012-11-15 | 2013-06-12 | 天津大学 | 基于长喉颈文丘里的双差压湿气流量测量装置 |
-
2015
- 2015-10-16 CN CN201510674413.2A patent/CN105387895B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101413817A (zh) * | 2008-12-05 | 2009-04-22 | 天津大学 | 双差压节流湿气测量方法 |
CN103048019A (zh) * | 2012-11-15 | 2013-04-17 | 天津大学 | 基于长喉颈文丘里的双差压湿气流量测量装置 |
CN202994216U (zh) * | 2012-11-15 | 2013-06-12 | 天津大学 | 基于长喉颈文丘里的双差压湿气流量测量装置 |
Non-Patent Citations (2)
Title |
---|
Dual-iterative model for gas condensate measurement based on void fraction;Chao Yuan 等;《Journal of Natural Gas Science and Engineering》;20150408;第24卷;全文 * |
长喉颈文丘里喉部取压位置对湿气测量模型影响的试验研究;徐英 等;《机械工程学报》;20140630;第50卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105387895A (zh) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105387895B (zh) | 基于流动密度的长喉颈文丘里管湿气流量测量方法 | |
CN103292849B (zh) | 一种在线测量水平管中湿气的气相流量和液相流量的装置及方法 | |
Lien et al. | The entrance length for fully developed turbulent channel flow | |
Zhou et al. | Experimental measurement of wind-sand flux and sand transport for naturally mixed sands | |
CN104453861B (zh) | 一种高压气井井筒温度分布的确定方法以及系统 | |
WO2013102312A1 (zh) | 一种蒸汽流量计量装置及计量方法 | |
CN105675070A (zh) | 异型文丘里流量计和用其测量多相流中气液相流量的方法 | |
CN101802563A (zh) | 两相流量计 | |
Hua et al. | Wet gas meter based on the vortex precession frequency and differential pressure combination of swirlmeter | |
Sun | Mass flow measurement of gas–liquid bubble flow with the combined use of a Venturi tube and a vortex flowmeter | |
US20180245463A1 (en) | Method for determining a water cut of an oil-water mixture produced from an oil well | |
CN103134834A (zh) | 一种湿蒸汽干度测量装置及方法 | |
CN103438942B (zh) | 双文丘里管 | |
CN106593408A (zh) | 油井流速的获取方法与装置 | |
Zheng et al. | Wet gas measurements of long-throat Venturi Tube based on forced annular flow | |
Xu et al. | An overreading model for nonstandard Venturi meters based on H correction factor | |
CN205483094U (zh) | 倒u型管与靶式流量计组合式天然气湿气流量测量系统 | |
Adefila et al. | Flow measurement of wet CO2 using an averaging pitot tube and coriolis mass flowmeters | |
Hardie | Developing measurement facilities for carbon capture and storage | |
CN102590057B (zh) | 吸阻标准棒吸阻值的测量方法和用于验证吸阻标准棒校准值的方法 | |
Wang et al. | High pressure steam–water two-phase flow measurements by flow division and separation method | |
Kelly et al. | Turbulent interchange in triangular array bare rod bundles | |
CN206470284U (zh) | 管道内气体流速的测量装置 | |
Wibowo et al. | Production test analysis of XYZ-Well at Deing geothermal field using horizontal discharge lip pressure method with Russel James equation and Hiriart equation | |
Li et al. | Mass flowrate measurement of wet steam using combined V-cone and vortex flowmeters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180619 Termination date: 20181016 |