CN110017873B - 一种基于界面波的气液两相流流量测量方法 - Google Patents
一种基于界面波的气液两相流流量测量方法 Download PDFInfo
- Publication number
- CN110017873B CN110017873B CN201910144260.9A CN201910144260A CN110017873B CN 110017873 B CN110017873 B CN 110017873B CN 201910144260 A CN201910144260 A CN 201910144260A CN 110017873 B CN110017873 B CN 110017873B
- Authority
- CN
- China
- Prior art keywords
- liquid
- gas
- phase
- flow velocity
- velocity
- 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
- 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/704—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F7/00—Volume-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 (5)
1.一种基于界面波的气液两相流流量测量方法,其特殊之处在于,包括以下步骤:
1)获取气-液界面的波动信息,得到两个特征参数,所述两个特征参数为平均液膜厚度δ和波速uw;
2)分别建立平均液膜厚度δ和波速uw与气相表观流速usg、液相表观流速usl之间的关系式:
Relf=(1-E)Rel (4)
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,l表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
其中,C表示权重系数;
3)联立关系式求解,得到气相表观流速usg、液相表观流速usl;
4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql。
2.根据权利要求1所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤1)中:
获取气-液界面的波动信息,是指通过界面波测量设备测得气-液界面的瞬时波动信息。
3.根据权利要求1所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤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.根据权利要求3所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql,具体为:
Qg=Ausg (9)
Ql=Ausl (10)
其中A表示管道的横截面积。
5.根据权利要求2所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤1)中测量瞬时气液界面的仪器具体为超声波脉冲发射接收仪。
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 CN110017873A (zh) | 2019-07-16 |
CN110017873B true 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) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794114B (zh) * | 2019-10-18 | 2021-09-28 | 河海大学 | 一种基于电磁感应原理测量卵石推移质输沙率的方法 |
CN111351541B (zh) * | 2020-02-12 | 2021-06-25 | 清华大学深圳国际研究生院 | 基于界面波和差压流量计的气液两相流流量测量方法 |
CN112945318B (zh) * | 2021-01-29 | 2022-12-06 | 天津大学 | 基于液膜厚度测量和涡街流量计的湿气分相流量测量方法 |
CN113049047B (zh) * | 2021-01-29 | 2022-12-02 | 天津大学 | 基于牛顿迭代的涡街湿气分相流量测量方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424406A (zh) * | 2013-09-03 | 2013-12-04 | 上海理工大学 | 管道内气液两相流的图像法测量装置及方法 |
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 |
CN104502240A (zh) * | 2014-12-17 | 2015-04-08 | 中国科学院工程热物理研究所 | 液体参数测量*** |
CN107270980A (zh) * | 2017-07-21 | 2017-10-20 | 西安交通大学 | 一种垂直管内气液两相环状流液膜厚度及流量的测量方法 |
CN108090321A (zh) * | 2017-12-20 | 2018-05-29 | 西南石油大学 | 一种水平气井倾斜管液膜携带临界气流量计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101883967B (zh) * | 2007-12-05 | 2012-11-28 | 普拉德研究及开发股份有限公司 | 超声波夹钳式多相流量计 |
CN106932444B (zh) * | 2017-03-14 | 2019-10-18 | 中国石油大学(北京) | 一种基于电导率的混油界面检测实验装置及实验方法 |
-
2019
- 2019-02-27 CN CN201910144260.9A patent/CN110017873B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 中国科学院工程热物理研究所 | 液体参数测量*** |
CN107270980A (zh) * | 2017-07-21 | 2017-10-20 | 西安交通大学 | 一种垂直管内气液两相环状流液膜厚度及流量的测量方法 |
CN108090321A (zh) * | 2017-12-20 | 2018-05-29 | 西南石油大学 | 一种水平气井倾斜管液膜携带临界气流量计算方法 |
Non-Patent Citations (2)
Title |
---|
"Film Thickness and Wave Velocity Measurements in a Vertical Duct";Kumar R 等;《Journal of Fluids Engineering》;20021231;第124卷;第634页摘要、第640页第1栏、第641页第2栏、第642页 * |
"预测水平管气-液两相环状流液膜分布的简化模型";梁法春 等;《长江大学学报》;20081231;第5卷(第3期);第37-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110017873A (zh) | 2019-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110017873B (zh) | 一种基于界面波的气液两相流流量测量方法 | |
US10472957B2 (en) | Entropy based multiphase flow detection | |
CN103513270B (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
RU2499229C2 (ru) | Способ и устройство для определения состава и расхода влажного газа | |
CN104729972B (zh) | 确定碳酸盐岩气水相对渗透率的方法及装置 | |
CN101556234A (zh) | 气水两相流多参数测量方法和装置 | |
CN105486358B (zh) | 基于文丘里管双差压的气液两相流参数测量方法 | |
CN206038586U (zh) | 气藏水侵物理模拟实验装置 | |
CN105004837B (zh) | 天然气水合物沉积物多测量单元分析方法 | |
CN104389586A (zh) | 油井动液面深度测量装置及方法 | |
CN107655533A (zh) | 一种基于回波能量积分的气体超声波流量计信号处理方法和*** | |
CN101900589B (zh) | 基于质量流量计的夹气液体流量测量方法 | |
CN203965422U (zh) | 一种井***水在线分析装置 | |
CN104237107A (zh) | 地层中低渗透率储层的视渗透率解释方法及*** | |
CN102538913B (zh) | 单声道超声水表流量测量特性校正方法 | |
CN111351541B (zh) | 基于界面波和差压流量计的气液两相流流量测量方法 | |
US20140090835A1 (en) | Systems and methods for the determination of gas permeability | |
Spina | Organized structures in a supersonic turbulent boundary layer | |
CN106767693A (zh) | 一种新型气压计及其高度测量方法 | |
US20150160057A1 (en) | Systems and methods for determining mass flow measurements of fluid flows | |
Qi et al. | Engineering simulation tests on multiphase flow in middle-and high-yield slanted well bores | |
CN108871152A (zh) | 一种管道内径测量设备及方法 | |
KR101522249B1 (ko) | 초음파를 이용한 질량유량측정 프로그램 및 이를 이용한 측정장치 | |
CN107816345A (zh) | 一种油井套管气气量计量的装置和方法 | |
CN203534679U (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 |