CN110017873A - 一种基于界面波的气液两相流流量测量方法 - Google Patents

一种基于界面波的气液两相流流量测量方法 Download PDF

Info

Publication number
CN110017873A
CN110017873A CN201910144260.9A CN201910144260A CN110017873A CN 110017873 A CN110017873 A CN 110017873A CN 201910144260 A CN201910144260 A CN 201910144260A CN 110017873 A CN110017873 A CN 110017873A
Authority
CN
China
Prior art keywords
liquid
flow
apparent velocity
velocity
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.)
Granted
Application number
CN201910144260.9A
Other languages
English (en)
Other versions
CN110017873B (zh
Inventor
陈威
朱晓玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua Research Institute Of Pearl River Delta
Shenzhen Star Technology Co Ltd
Original Assignee
Tsinghua Research Institute Of Pearl River Delta
Shenzhen Star Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua Research Institute Of Pearl River Delta, Shenzhen Star Technology Co Ltd filed Critical Tsinghua Research Institute Of Pearl River Delta
Priority to CN201910144260.9A priority Critical patent/CN110017873B/zh
Publication of CN110017873A publication Critical patent/CN110017873A/zh
Application granted granted Critical
Publication of CN110017873B publication Critical patent/CN110017873B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/704Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F7/00Volume-flow measuring devices with two or more measuring ranges; Compound meters

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明涉及一种基于界面波的气液两相流流量测量方法,其主要包括以下步骤:1)获取气‑液界面的波动信息,得到两个特征参数,所述两个特征参数为平均液膜厚度δ和波速uw;2)分别建立平均液膜厚度δ和波速uw与气相表观流速usg、液相表观流速usl之间的关系式;3)联立关系式求解,得到气相表观流速usg、液相表观流速usl;4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql。本发明可以同时测量气相流量和液相流量,利用本发明的方法,只需获取气‑液界面的波动信息中的两个特征参数,便可得到气液两相流流量,方法简单,在得到同样结果的前提下,所需要设备更少,可以减少投入成本。

Description

一种基于界面波的气液两相流流量测量方法
技术领域
本发明涉及一种基于界面波的气液两相流流量测量方法。
背景技术
油气能源是一个国家工业和经济发展的重要能源,关系着国家的经济命脉,其重要性不言而喻。我国是世界上石油生产和需求大国,据统计,2017年我国石油产量为1.92亿吨,天然气产量为1330.07亿立方米,位居世界前列。油气能源产量巨大,使得对油气能源的准确测量和高效管理显得愈发重要。实现智慧油田,对油田开采进行数字化管理,成为了时代发展的趋势,同时也给多相流测量技术带来了巨大的机遇和挑战。
目前,油田现场对于流量的计量仍然以传统的分离式方法为主,体积庞大、价格昂贵的弊端严重阻碍了该方法的使用和发展,因此研发成本低廉、安装方便、在线测量的不分离式多相流测量技术变得越来越紧迫。经过几十年的发展,现有的测量方法/设备是多种多样的,但是往往一种测量方法/设备只能测得一种流量,难以满足同时测量多相流量的需求,而将多种测量方法/设备加以组合来实现多相流量的测量又会导致设备体积庞大、成本高昂、安装困难等弊端,背离了不分离式测量方法的研发初衷。因此,研发能够同时测得多相流量的测量方法/设备,是不分离式多相流测量技术的发展趋势。
发明内容
本发明提出一种基于界面波的气液两相流流量测量方法,其解决了现有多相流测量方法难以同时测量气相流量和液相流量的技术问题。
本发明解决上述问题的技术方案是:一种基于界面波的气液两相流流量测量方法,其特殊之处在于,包括以下步骤:
1)获取气-液界面的波动信息,得到两个特征参数,所述两个特征参数为平均液膜厚度δ和波速uw
2)分别建立平均液膜厚度δ和波速uw与气相表观流速usg、液相表观流速usl之间的关系式;
3)联立关系式求解,得到气相表观流速usg、液相表观流速usl
4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql
进一步地,上述步骤1)中:
获取气-液界面的波动信息,是指通过界面波测量设备测得气-液界面的瞬时波动信息。
进一步地,上述步骤2)中:
平均液膜厚度δ与气相表观流速usg、液相表观流速usl的关系式为:
Relf=(1-E)Rel (4)
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,l表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
其中,C表示权重系数。
进一步地,上述步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl具体为:
3.1)联立关系式(1)~(8):
3.2)利用迭代算法求解:
3.21)根据经验设定usl的搜索范围[usl0,usl1],令usl=usl0
3.22)通过关系式(1)~(6)迭代,求得usg收敛值usg1
3.23)通过关系式(7)和(8)迭代,求得usg收敛值usg2
3.24)根据测量精度选取步长,在搜索范围内遍历usl,重复3.22)和3.23)两个步骤,计算每一种usl对应的usg1和usg2
3.25)找出usg1和usg2的所有组合中使得|usg1-usg2|最小的那一组,然后将二者的平均值作为usg的预测值,该组对应的usl作为usl的预测值。
进一步地,上述步骤4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql,具体为:
Qg=Ausg (9)
Ql=Ausl (10)
其中A表示管道的横截面积。
进一步地,上述步骤1)中:测量瞬时气液界面的仪器具体为超声波脉冲发射接收仪。
本发明的优点:
本发明提出一种基于界面波的气液两相流流量测量方法,可以同时测量气相流量和液相流量,利用本发明的方法,只需获取气-液界面的波动信息中的两个特征参数,便可得到气液两相流流量,方法简单,在得到同样结果的前提下,所需要设备更少,可以减少投入成本。
附图说明
图1是本发明实施例中的基于界面波的气液两相流流量测量方法原理图。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
参见图1,一种基于界面波的气液两相流流量测量方法,包括以下步骤:
1)获取气-液界面的波动信息,得到两个特征参数,所述两个特征参数为平均液膜厚度δ和波速uw
2)分别建立平均液膜厚度δ和波速uw与气相表观流速usg、液相表观流速usl之间的关系式;
3)联立关系式求解,得到气相表观流速usg、液相表观流速usl
4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql
进一步地,上述步骤1)中:
获取气-液界面的波动信息,是指通过界面波测量设备测得气-液界面的瞬时波动信息。这里,界面波测量设备可以是任意一种用于测量瞬时气液界面的仪器,最常用的有电导探针、高速摄像头、超声设备等,优选的,为广东汕头超声电子股份有限公司生产的超声波脉冲发射接收仪,产品型号:CTS-8077PR。
进一步地,上述步骤2)中:
平均液膜厚度δ与气相表观流速usg、液相表观流速usl的关系式为:
Relf=(1-E)Rel (4)
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,1表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
其中,C表示权重系数。
进一步地,上述步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl,具体为:
3.1)联立关系式(1)~(8):
3.2)利用迭代算法求解:
3.21)根据经验设定usl的搜索范围[usl0,usl1],令usl=usl0
3.22)通过关系式(1)~(6)迭代,求得usg收敛值usg1
3.23)通过关系式(7)和(8)迭代,求得usg收敛值usg2
3.24)根据测量精度选取步长,在搜索范围内遍历usl,重复3.22)和3.23)两个步骤,计算每一种usl对应的usg1和usg2
3.25)找出usg1和usg2的所有组合中使得|usg1-usg2|最小的那一组,然后将二者的平均值作为usg的预测值,该组对应的usl作为usl的预测值。
进一步地,上述步骤4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Q1,具体为:
Qg=Ausg (9)
Q1=Ausl (10)
其中A表示管道的横截面积。
以上所述仅为本发明的实施例,并非以此限制本发明的保护范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的系统领域,均同理包括在本发明的保护范围内。

Claims (6)

1.一种基于界面波的气液两相流流量测量方法,其特殊之处在于,包括以下步骤:
1)获取气-液界面的波动信息,得到两个特征参数,所述两个特征参数为平均液膜厚度δ和波速uw
2)分别建立平均液膜厚度δ和波速uw与气相表观流速usg、液相表观流速usl之间的关系式;
3)联立关系式求解,得到气相表观流速usg、液相表观流速usl
4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql
2.根据权利要求1所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤1)中:
获取气-液界面的波动信息,是指通过界面波测量设备测得气-液界面的瞬时波动信息。
3.根据权利要求1或2所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤2)中:
平均液膜厚度δ与气相表观流速usg、液相表观流速usl的关系式为:
Relf=(1-E)Rel (4)
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,l表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
其中,C表示权重系数。
4.根据权利要求3所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl,具体为:
3.1)联立关系式(1)~(8):
3.2)利用迭代算法求解:
3.21)根据经验设定usl的搜索范围[usl0,usl1],令usl=usl0
3.22)通过关系式(1)~(6)迭代,求得usg收敛值usg1
3.23)通过关系式(7)和(8)迭代,求得usg收敛值usg2
3.24)根据测量精度选取步长,在搜索范围内遍历usl,重复3.22)和3.23)两个步骤,计算每一种usl对应的usg1和usg2
3.25)找出usg1和usg2的所有组合中使得|usg1-usg2|最小的那一组,然后将二者的平均值作为usg的预测值,该组对应的usl作为usl的预测值。
5.根据权利要求4所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql,具体为:
Qg=Ausg (9)
Ql=Ausl (10)
其中A表示管道的横截面积。
6.根据权利要求2所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤1)中测量瞬时气液界面的仪器具体为超声波脉冲发射接收仪。
CN201910144260.9A 2019-02-27 2019-02-27 一种基于界面波的气液两相流流量测量方法 Active CN110017873B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910144260.9A CN110017873B (zh) 2019-02-27 2019-02-27 一种基于界面波的气液两相流流量测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910144260.9A CN110017873B (zh) 2019-02-27 2019-02-27 一种基于界面波的气液两相流流量测量方法

Publications (2)

Publication Number Publication Date
CN110017873A true CN110017873A (zh) 2019-07-16
CN110017873B CN110017873B (zh) 2020-09-04

Family

ID=67189121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910144260.9A Active CN110017873B (zh) 2019-02-27 2019-02-27 一种基于界面波的气液两相流流量测量方法

Country Status (1)

Country Link
CN (1) CN110017873B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110794114A (zh) * 2019-10-18 2020-02-14 河海大学 一种基于电磁感应原理测量卵石推移质输沙率的方法
CN111351541A (zh) * 2020-02-12 2020-06-30 清华大学深圳国际研究生院 基于界面波和差压流量计的气液两相流流量测量方法
CN112945318A (zh) * 2021-01-29 2021-06-11 天津大学 基于液膜厚度测量和涡街流量计的湿气分相流量测量方法
CN113049047A (zh) * 2021-01-29 2021-06-29 天津大学 基于牛顿迭代的涡街湿气分相流量测量方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100299088A1 (en) * 2007-12-05 2010-11-25 Schlomberger Technology Corporation Ultrasonic clamp-on multiphase flowmeter
CN103038611A (zh) * 2010-06-30 2013-04-10 普拉德研究及开发股份有限公司 多相流量计和液体膜测量方法
CN103424406A (zh) * 2013-09-03 2013-12-04 上海理工大学 管道内气液两相流的图像法测量装置及方法
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
CN104502240A (zh) * 2014-12-17 2015-04-08 中国科学院工程热物理研究所 液体参数测量系统
CN106932444A (zh) * 2017-03-14 2017-07-07 中国石油大学(北京) 一种基于电导率的混油界面检测实验装置及实验方法
CN107270980A (zh) * 2017-07-21 2017-10-20 西安交通大学 一种垂直管内气液两相环状流液膜厚度及流量的测量方法
CN108090321A (zh) * 2017-12-20 2018-05-29 西南石油大学 一种水平气井倾斜管液膜携带临界气流量计算方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100299088A1 (en) * 2007-12-05 2010-11-25 Schlomberger Technology Corporation Ultrasonic clamp-on multiphase flowmeter
CN103038611A (zh) * 2010-06-30 2013-04-10 普拉德研究及开发股份有限公司 多相流量计和液体膜测量方法
US20140013857A1 (en) * 2010-06-30 2014-01-16 Alexandre Lupeau Multiphase Flowmeter and Liquid Film Measurement Method
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
CN103424406A (zh) * 2013-09-03 2013-12-04 上海理工大学 管道内气液两相流的图像法测量装置及方法
CN104502240A (zh) * 2014-12-17 2015-04-08 中国科学院工程热物理研究所 液体参数测量系统
CN106932444A (zh) * 2017-03-14 2017-07-07 中国石油大学(北京) 一种基于电导率的混油界面检测实验装置及实验方法
CN107270980A (zh) * 2017-07-21 2017-10-20 西安交通大学 一种垂直管内气液两相环状流液膜厚度及流量的测量方法
CN108090321A (zh) * 2017-12-20 2018-05-29 西南石油大学 一种水平气井倾斜管液膜携带临界气流量计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
AMUMUL AHSAN: "《Two phase flow,phase change and numerical modeling》", 30 September 2011, INTECH *
KUMAR R 等: ""Film Thickness and Wave Velocity Measurements in a Vertical Duct"", 《JOURNAL OF FLUIDS ENGINEERING》 *
梁法春 等: ""预测水平管气-液两相环状流液膜分布的简化模型"", 《长江大学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110794114A (zh) * 2019-10-18 2020-02-14 河海大学 一种基于电磁感应原理测量卵石推移质输沙率的方法
CN110794114B (zh) * 2019-10-18 2021-09-28 河海大学 一种基于电磁感应原理测量卵石推移质输沙率的方法
CN111351541A (zh) * 2020-02-12 2020-06-30 清华大学深圳国际研究生院 基于界面波和差压流量计的气液两相流流量测量方法
CN111351541B (zh) * 2020-02-12 2021-06-25 清华大学深圳国际研究生院 基于界面波和差压流量计的气液两相流流量测量方法
CN112945318A (zh) * 2021-01-29 2021-06-11 天津大学 基于液膜厚度测量和涡街流量计的湿气分相流量测量方法
CN113049047A (zh) * 2021-01-29 2021-06-29 天津大学 基于牛顿迭代的涡街湿气分相流量测量方法
CN112945318B (zh) * 2021-01-29 2022-12-06 天津大学 基于液膜厚度测量和涡街流量计的湿气分相流量测量方法

Also Published As

Publication number Publication date
CN110017873B (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
CN110017873A (zh) 一种基于界面波的气液两相流流量测量方法
CN101354273B (zh) 复合式气体流量测量方法及其装置
CN1134651C (zh) 分流型流量计
RU2499229C2 (ru) Способ и устройство для определения состава и расхода влажного газа
CN201241692Y (zh) 湿天然气多相流量计
CN105675070A (zh) 异型文丘里流量计和用其测量多相流中气液相流量的方法
CN103038609A (zh) 用于测量多相流体混合物的至少一个特征值的装置
Sun Mass flow measurement of gas–liquid bubble flow with the combined use of a Venturi tube and a vortex flowmeter
CN106643945B (zh) 一种同质气液混合介质质量流量测试装置及方法
CN107328447A (zh) 井口气液两相流量计量方法及装置
CN101900589B (zh) 基于质量流量计的夹气液体流量测量方法
CN102288730A (zh) 雾滴含量测量的方法及其雾滴含量测量装置
CN201032457Y (zh) 基于标准文丘里管与涡街流量计的低含气率气液两相流测量装置
CN102628702B (zh) 油水两相流部分分离在线计量的装置的应用方法
CN201242456Y (zh) 一种基于涡街、均速管、旁路管的气体流量组合测量装置
CN208721144U (zh) 基于多孔节流与声发射技术的两相流参数检测装置
CN207439632U (zh) 可视化先导钴靶件钴棒流致振动实验装置
CN103323066B (zh) 一种低含液率气液两相流测量方法及测量系统
CN202471152U (zh) 一种蒸汽流量计量装置
Falcone et al. ANUMET: A novel wet gas flowmeter
CN107830979A (zh) 可视化先导钴靶件钴棒流致振动实验装置
CN105444828B (zh) 超声波流量计的整流器的设计方法
CN103353319A (zh) 基于直通型气体超声波流量计的湿气流量测量方法
CN207798137U (zh) 一种地面试油过程中提高测量精度的装置
CN203037301U (zh) 天然气贸易流量计

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