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

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

Info

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
Application number
CN201910144260.9A
Other languages
English (en)
Other versions
CN110017873A (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.)
Research Institute Of Tsinghua Pearl River Delta
Shenzhen Leengstar Technology Co ltd
Original Assignee
Research Institute Of Tsinghua Pearl River Delta
Shenzhen Leengstar 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 Research Institute Of Tsinghua Pearl River Delta, Shenzhen Leengstar Technology Co ltd filed Critical Research Institute Of Tsinghua 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

Images

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的关系式为:
Figure BDA0001979491990000021
Figure BDA0001979491990000022
Figure BDA0001979491990000023
Relf=(1-E)Rel (4)
Figure BDA0001979491990000024
Figure BDA0001979491990000025
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,l表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
Figure BDA0001979491990000026
Figure BDA0001979491990000027
其中,C表示权重系数。
进一步地,上述步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl具体为:
3.1)联立关系式(1)~(8):
Figure BDA0001979491990000031
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的关系式为:
Figure BDA0001979491990000051
Figure BDA0001979491990000052
Figure BDA0001979491990000053
Relf=(1-E)Rel (4)
Figure BDA0001979491990000054
Figure BDA0001979491990000055
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,1表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
Figure BDA0001979491990000056
Figure BDA0001979491990000057
其中,C表示权重系数。
进一步地,上述步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl,具体为:
3.1)联立关系式(1)~(8):
Figure BDA0001979491990000061
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之间的关系式:
Figure FDA0002562621220000011
Figure FDA0002562621220000012
Figure FDA0002562621220000013
Relf=(1-E)Rel (4)
Figure FDA0002562621220000014
Figure FDA0002562621220000015
其中D表示管道内径,Re表示雷诺数,ρ表示密度,μ表示粘度,E表示液滴夹带率,We表示韦伯数,σ表示表面张力,下标lf表示液膜,g表示气相,l表示液相;
波速uw与气相表观流速usg、液相表观流速usl的关系式为:
Figure FDA0002562621220000016
Figure FDA0002562621220000017
其中,C表示权重系数;
3)联立关系式求解,得到气相表观流速usg、液相表观流速usl
4)通过气相表观流速usg、液相表观流速usl计算出气相体积流量Qg、液相体积流量Ql
2.根据权利要求1所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤1)中:
获取气-液界面的波动信息,是指通过界面波测量设备测得气-液界面的瞬时波动信息。
3.根据权利要求1所述的一种基于界面波的气液两相流流量测量方法,其特殊之处在于:步骤3)中联立关系式求解,得到气相表观流速usg、液相表观流速usl,具体为:
3.1)联立关系式(1)~(8):
Figure FDA0002562621220000021
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)中测量瞬时气液界面的仪器具体为超声波脉冲发射接收仪。
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 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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101883967B (zh) * 2007-12-05 2012-11-28 普拉德研究及开发股份有限公司 超声波夹钳式多相流量计
CN106932444B (zh) * 2017-03-14 2019-10-18 中国石油大学(北京) 一种基于电导率的混油界面检测实验装置及实验方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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