CN115096383B - 基于等效密度测算多相流中气相流量的方法 - Google Patents
基于等效密度测算多相流中气相流量的方法 Download PDFInfo
- Publication number
- CN115096383B CN115096383B CN202210830795.3A CN202210830795A CN115096383B CN 115096383 B CN115096383 B CN 115096383B CN 202210830795 A CN202210830795 A CN 202210830795A CN 115096383 B CN115096383 B CN 115096383B
- Authority
- CN
- China
- Prior art keywords
- historical
- flow
- gas
- calculating
- real
- 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/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
-
- 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
- G01F1/40—Details of construction of the flow constriction devices
- G01F1/44—Venturi tubes
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于等效密度测算多相流中气相流量的方法,首先利用历史数据计算历史等效气密度;其次对历史等效气密度进行归一化处理,得到归一化历史等效气密度;再次,对归一化历史等效气密度和历史含气率进行拟合,以线性函数表达二者关系;最后,根据实时数据反算实时等效气密度,由此算出实时气相体积流量。采用本发明的显著效果是,避免了采用经验公式计算时,由于多变量的影响而造成的滑脱比、干度和气相流量计量误差。本方法以一种全新的方式测算多相流中的气相流量,计量结果误差较小,具有实际应用前景。
Description
技术领域
本发明涉及石油生产计量,具体涉及多相流的流量测定方法。
背景技术
在原油开采过程中,为了确定各油井的油、气产量,了解地层油气含量及地层结构的变化,需要对油井产出流体中各组分的体积流量或质量流量进行连续的计量并提供实时计量数据,以优化生产参数,提高采收率。近年来,随着全球能源问题的日益严重和石油开采技术的不断发展,对原油计量提出了更高的要求。相对于传统三相分离器,文丘里伽马多相流量计(能实时、连续测量流量、相分率等数据)对被测油气水混合物不用进行相分离,现场安装工艺简洁,结构紧凑,占空间小;基本上可以做到无人值守,尤其在海洋石油和油井测试中具有很大的经济效益。
众所周知,在低压力操作条件下,用文丘里流量计测量气液两相流时,由于气液之间存在显著的相滑脱,通常会造成较大的气流量系统偏差,这需要引入适当的模型来提高气量测量精度。常用的方式是从伯努利方程出发,基于均相流模型计算多相流的总质量流量,但在计算气相流量时必须引入“干度”数据。干度无法进行实时、直接测量,其需要知道气液滑脱比才能计算。而滑脱比也无法直接测量,必须用经验公式来描述,为提升经验公式的合理性和计算准确性,国内外众多学者都致力于相关研究;但不论经验公式如何修正、调整,其涉及到的变量太多,实际运用效果仍不理想。因此通过总流量和干度确定气相流量的方法在实际应用时还是存在一些问题。
发明内容
有鉴于以上问题,本发明目的在于提出一种新颖的、简单可行的相流量计量方法,其采用文丘里流量计和伽马射线收发机构来解决多相流中气量计量问题。初步研究分析表明此方法中的变量较少,与文丘里的结构和介质物理性质的关系比较弱。特别是对于低操作压力下的气相流量测量,这一方法有较好的应用前景。
所采用的主要技术方案如下:
一种基于等效密度测算多相流中气相流量的方法,按以下步骤进行:
公式①中:
其中ε为可膨胀系数;
C 为文丘里流出系数;
A 为文丘里喉部截面积;
d 为文丘里喉部内直径;
D 为文丘里入口内直径;
公式③中:
a、b为常数,通过多组数据拟合得到;
附图说明
具体实施方式
以下结合实施例和附图对本发明作进一步说明。
实施例1:
一种基于等效密度测算多相流中气相流量的方法,按以下步骤进行:
所述文丘里伽马流量计包括文丘里流量测定模块、伽马相分率测定模块、信号收
发模块、计算单元等,其具体结构在现有技术中多有记载,在此不做赘述;具体来讲,由文丘
里流量测定模块测定所述历史差压,由伽马相分率测定模块测定所述历史含气率,
由所述计算单元输出所述历史气相体积流量;而历史工况气密度则是标定得到,同理
历史工况液相密度也是标记得到;
公式①中:
K是与文丘里相关的参数,均可通过现有技术手段测定得到;
其中:
ε为可膨胀系数;
C 为文丘里流出系数;
A 为文丘里喉部截面积;
d 为文丘里喉部内直径;
D 为文丘里入口内直径;
公式③中:
a是拟合直线的斜率;
b为拟合直线的纵截距;
步骤一中所述的历史数据组可以是历史计量时留存的数据;也可以是为了进行以上步骤二、三、四专门进行的测试数据,例如在实验室或检测单位的配备的多相流试验环道上进行测试。主要目的是为线性拟合提供原始数据,从而得到a、b的值。
公式组④中个参数的定义与以上描述保持一致;在此不再赘述。
步骤五中用到的文丘里伽马流量计与步骤一中用到的文丘里伽马流量计可以是相同参数型号,也可以是不同参数型号。
实施例2:
公式⑤中:
ρ为多相流的实时工况混合密度;
实施例3:
一种多相流气气液滑脱比计量方法,包括按照实施例2所述的方法计算干度x;还包括按以下公式⑥计算气液滑脱比S;
试验例:
试验介质有:介质①“原油-天然气-产出水”以及介质②“变压器油-空气-水”。
文丘里均竖直安装,DN50(0.5)的伽马传感器安装在文丘里的上游入口,其它所有DN80、DN50文丘里均在其喉部安装单能伽玛传感器来测量截面含气率。
流态稳定后,每个试验点测试时间均为15-30分钟。
实验参数范围:
操作压力:1.5-6.0Bar;
DN80液体流量范围:7.0-74.0m3/h;
DN50液体流量范围:4.5-27.0m3/h;
DN80气体流量范围:7.0-380.0m3/h;
DN50气体流量范围:4.6-82.0m3/h;
干度范围:0-0.09;
截面含气率范围:0-96%;
体积含气率范围:0-98% 。
油、气、水介质属性见表1。
表1、介质①、②属性表
介质① | 介质② | |
油密度@20℃(kg/m<sup>3</sup>) | 856.0 (原油) | 875.2 (变压器油) |
水密度@20℃(kg/m<sup>3</sup>) | 1001.5 (地层水) | 998.5 (自来水) |
气密度@20℃(kg/m<sup>3</sup>) | 0.902 (天然气) | 1.204 (空气) |
从图1可以看出,所有的数据点(归一化历史等效气密度的自然对数)都基本分布
在与历史含气率成直线关系的一条带上;因而,将归一化历史等效气密度的自然对数与对应的历史含气率进行线性拟合是合理的;而作为反向计算,可以预见的是,
采用公式④计算实时气相体积流量也是可行的。
以DN80文丘里管为例:
将采用DN80文丘里管的测算得到的归一化历史等效气密度的自然对数和
对应的历史含气率绘制在平面直角坐标系上,如图2所示。从图2中可以看出,当具体到
DN80文丘里管时,数据点的离散性更小,归一化历史等效气密度的自然对数和对应
的历史含气率之间的线性关系更明确;可预料的是:采用公式④计算的实时气相体积流
量应该会更准确。拟合DN80的所有数据,按照公式③进行线性回归得到的归一化历史等效
气密度的自然对数和对应的历史含气率的具体关系是:
由此确定参数a=-8.514;b=10.184。
约20%数据量的相对误差在±5%以内(含5%);
约50%数据量的相对误差在±10%以内(含10%);
约75%数据量的相对误差在±15%以内(含15%);
约95%数据量的相对误差在±20%以内(含20%)。
由以上结果可知:就DN80文丘里管的测算结果而言,本发明方法没有表现出对文丘里的喉径比和介质物理性质的强的依赖性。虽然部分数据的相对误差表现还不够完美,但这样的误差已经具有实用意义。
以DN50文丘里管为例:
将采用DN50文丘里管的测算得到的归一化历史等效气密度的自然对数和
对应的历史含气率绘制在平面直角坐标系上,如图3所示。从图3中可以看出,DN50文丘里
管运行测算后,数据点的离散性较DN80大一些。拟合DN50文丘里管测算的所有数据,按照公
式③进行线性回归得到的归一化历史等效气密度的自然对数和对应的历史含气
率的具体关系是:
由此确定参数a=-9.1227;b=10.218。
约25%数据量的相对误差在±5%以内(含5%);
约65%数据量的相对误差在±10%以内(含10%);
约75%数据量的相对误差在±15%以内(含15%);
约95%数据量的相对误差在±20%以内(含20%)。
由以上结果可知:就DN50文丘里管的测算结果而言,本发明方法没有表现出对文丘里的喉径比和介质物理性质的强的依赖性。虽然部分数据的相对误差表现还不够完美,但这样的误差已经具有实用意义。
有益效果:采用本发明的方法,避免了采用经验公式计算时,由于多变量的影响而造成的滑脱比、干度和气相流量计量误差。本方法以一种全新的方式测算多相流中的气相流量,计量结果误差较小,具有实际应用前景。
最后需要说明的是,上述描述仅仅为本发明的优选实施例,本领域的普通技术人员在本发明的启示下,在不违背本发明宗旨及权利要求的前提下,可以做出多种类似的表示,这样的变换均落入本发明的保护范围之内。
Claims (4)
1.一种基于等效密度测算多相流中气相流量的方法,其特征在于按以下步骤进行:
公式①中:
其中ε为可膨胀系数;
C 为文丘里流出系数;
A 为文丘里喉部截面积;
d 为文丘里喉部内直径;
D 为文丘里入口内直径;
公式③中:
a、b为常数,通过多组数据拟合得到;
公式组④中:
4.根据权利要求1、2或3所述的基于等效密度测算多相流中气相流量的方法,其特征在于:所述文丘里伽马流量计包括文丘里流量测定模块和伽马相分率测定模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210830795.3A CN115096383B (zh) | 2022-07-15 | 2022-07-15 | 基于等效密度测算多相流中气相流量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210830795.3A CN115096383B (zh) | 2022-07-15 | 2022-07-15 | 基于等效密度测算多相流中气相流量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115096383A CN115096383A (zh) | 2022-09-23 |
CN115096383B true CN115096383B (zh) | 2022-11-22 |
Family
ID=83297291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210830795.3A Active CN115096383B (zh) | 2022-07-15 | 2022-07-15 | 基于等效密度测算多相流中气相流量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115096383B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115420342B (zh) * | 2022-11-03 | 2023-03-24 | 海默新宸水下技术(上海)有限公司 | 一种基于含气率拟合的湿天然气计量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0372598A2 (en) * | 1988-12-03 | 1990-06-13 | Schlumberger Limited | Impedance cross correlation logging tool |
CN101517379A (zh) * | 2006-07-21 | 2009-08-26 | 因万西斯系统股份有限公司 | 多相科里奥利流量计 |
CN101839738A (zh) * | 2010-04-30 | 2010-09-22 | 西北工业大学 | 一种湿蒸汽流量仪及测量方法 |
CN102749111A (zh) * | 2012-07-24 | 2012-10-24 | 兰州海默科技股份有限公司 | 湿气流量测量方法及其装置 |
CN103674139A (zh) * | 2013-12-11 | 2014-03-26 | 中国石油大学(华东) | 基于质量流量计参数统计特征的两相流流量测量方法 |
CN110411521A (zh) * | 2019-06-28 | 2019-11-05 | 中国石油大学(北京) | 一种基于双喷嘴的油井多相流分相含率在线计量方法 |
CN110726444A (zh) * | 2019-08-22 | 2020-01-24 | 无锡洋湃科技有限公司 | 一种基于科氏力质量流量计的湿气流量计量方法及其装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2399641B (en) * | 2003-03-18 | 2005-08-31 | Schlumberger Holdings | Method and apparatus for determining the gas flow rate of a gas-liquid mixture |
GB2432425B (en) * | 2005-11-22 | 2008-01-09 | Schlumberger Holdings | Isokinetic sampling method and system for multiphase flow from subterranean wells |
CN102435245B (zh) * | 2012-01-06 | 2014-01-15 | 兰州海默科技股份有限公司 | 一种蒸汽流量计量装置及计量方法 |
CN103292849B (zh) * | 2013-03-25 | 2015-12-23 | 兰州海默科技股份有限公司 | 一种在线测量水平管中湿气的气相流量和液相流量的装置及方法 |
GB2572836B (en) * | 2018-09-13 | 2020-09-02 | M-Flow Tech Ltd | Void fraction calibration method |
-
2022
- 2022-07-15 CN CN202210830795.3A patent/CN115096383B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0372598A2 (en) * | 1988-12-03 | 1990-06-13 | Schlumberger Limited | Impedance cross correlation logging tool |
CN101517379A (zh) * | 2006-07-21 | 2009-08-26 | 因万西斯系统股份有限公司 | 多相科里奥利流量计 |
CN101839738A (zh) * | 2010-04-30 | 2010-09-22 | 西北工业大学 | 一种湿蒸汽流量仪及测量方法 |
CN102749111A (zh) * | 2012-07-24 | 2012-10-24 | 兰州海默科技股份有限公司 | 湿气流量测量方法及其装置 |
CN103674139A (zh) * | 2013-12-11 | 2014-03-26 | 中国石油大学(华东) | 基于质量流量计参数统计特征的两相流流量测量方法 |
CN110411521A (zh) * | 2019-06-28 | 2019-11-05 | 中国石油大学(北京) | 一种基于双喷嘴的油井多相流分相含率在线计量方法 |
CN110726444A (zh) * | 2019-08-22 | 2020-01-24 | 无锡洋湃科技有限公司 | 一种基于科氏力质量流量计的湿气流量计量方法及其装置 |
Non-Patent Citations (5)
Title |
---|
A sifting method of the void fraction and the bubble density with bubble size in gas-liquid two phase flow;Sekoguchi K.等;《Technology Reports of the Kyushu University》;19780101;第51卷(第1期);第85-91页 * |
利用γ射线测量原油含水率和含气率的数学算法分析;白秋果 等;《核电子学与探测技术》;20020520;第22卷(第03期);第34-36页 * |
基于射线技术的低压竖直管气液两相流模型研究;潘艳芝 等;《工程热物理学报》;20201130;第41卷(第11期);第2756-2761页 * |
基于狭缝文丘里的气液两相流检测技术研究;方立德;《中国博士学位论文全文数据库基础科学辑》;20090715(第07期);第A004-14页 * |
应用文丘里管和空泡份额传感器测量油气两相流的实验研究;贾志海 等;《高校化学工程学报》;20090615;第23卷(第03期);第51-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115096383A (zh) | 2022-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101413817B (zh) | 双差压节流湿气测量方法 | |
CN100434870C (zh) | 气液多相流分离整流装置及其测量装置 | |
CN107843297B (zh) | 基于v锥的低含气率气液两相流液相流量在线测量装置及方法 | |
US20150134275A1 (en) | Steam flow rate metering device and metering method therefor | |
CN109506729B (zh) | 一种气液两相流参数在线检测方法及装置 | |
CN110514257B (zh) | 一种基于文丘里的低含液率湿气两相流量测量装置及方法 | |
CN115096383B (zh) | 基于等效密度测算多相流中气相流量的方法 | |
CN110411521B (zh) | 一种基于双喷嘴的油井多相流分相含率在线计量方法 | |
CN110726444B (zh) | 一种基于科氏力质量流量计的湿气流量计量方法及其装置 | |
CN107290008A (zh) | 一种孔板涡街气液计量装置及其计算方法 | |
CN101666770A (zh) | 低含气率原油测量装置及方法 | |
CN205477586U (zh) | 管式油井气液两相流量计 | |
CN112710703B (zh) | 一种带有导电特性补偿的电导网格传感器三相流成像方法 | |
CN107939367B (zh) | 一种压裂水马力确定方法 | |
CN112129363A (zh) | 基于电容传感器的气液两相流测量装置 | |
CN111323087B (zh) | 基于强制环状流技术的湿气流量实时测量装置及方法 | |
CN108254034B (zh) | 一种油气水分相流量在线计量装置及其计量方法 | |
CN103323066B (zh) | 一种低含液率气液两相流测量方法及测量系统 | |
CN115855187A (zh) | 一种基于折算滑速比拟合的湿天然气计量方法 | |
CN206959918U (zh) | 一种孔板涡街气液计量装置 | |
CN115790758A (zh) | 基于温度补偿的伽马传感器计数矫正方法 | |
CN209446110U (zh) | 一种气液两相流参数在线检测装置 | |
CN111504407A (zh) | 一种产液气藏井中气液两相流量的在线测量方法 | |
CN111442809B (zh) | 一种基于竖直文丘里管的气液段塞流流量在线测量方法 | |
CN103674139A (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 |