CN103674139A - 基于质量流量计参数统计特征的两相流流量测量方法 - Google Patents
基于质量流量计参数统计特征的两相流流量测量方法 Download PDFInfo
- Publication number
- CN103674139A CN103674139A CN201310670292.5A CN201310670292A CN103674139A CN 103674139 A CN103674139 A CN 103674139A CN 201310670292 A CN201310670292 A CN 201310670292A CN 103674139 A CN103674139 A CN 103674139A
- Authority
- CN
- China
- Prior art keywords
- flow
- density
- mass flowmeter
- gas
- liquid
- 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
Links
Images
Abstract
本发明属于流体测量技术领域,公开了一种基于质量流量计参数统计特征的两相流流量测量方法,该方法应用科氏质量流量计测得的瞬时质量流量和瞬时密度及其统计特征,结合存储于计算机中的关联关系及各关系式系数,实现气液两相流含率及流量的测量。它具体包括以下步骤:(1)建立流量和密度的统计特征与气液两相流含气率的关联关系;(2)建立密度降与质量降的关联关系;(3)建立科氏质量流量计密度示值与气液两相流混合密度的关联关系;(4)根据建立的各关联关系估计气液两相流含气率、气相质量流量和液相质量流量。本发明的有益效果在于:该测量方法无需采用高效气液分离器进行气液分离,成本低、流量检测精度高、实时性好、易于实现。适用于气液两相流的含气率及分相流量的测量。
Description
技术领域
本发明属于流体测量技术领域,具体涉及到一种基于质量流量计参数统计特征的两相流流量测量方法。
背景技术
气液两相流广泛存在于石油、化工、冶金、能源、动力、轻工等部门,其应用的普遍性及重要性促使该领域研究工作的迅速发展。流量是气液两相流系统中具有重要科学和工程意义的参数,对气液两相流应用系统的计量、控制、可靠性和效率等均具有重大影响,长期以来,一直是两相流领域的重要研究方向。
气液两相流量测量是一个亟待解决的问题,国内外学者对其进行了大量的研究,其中主要方法有分离计量法、单相流量计法、节流法、相关测量法等。
分离计量法采用庞大的分离设备,把两相流体流量测量转化成单相流体流量测量,具有工作可靠、测量精度高、测量范围宽且不受气液两相流流型变化影响等优点。这种方法的最大缺点是分离设备体积庞大,价格昂贵,并需要建立专门的计量站和测试管线,这在很大程度上增加了成本。
单相流量计法是将单相流流量测量仪表应用到两相流流量测量中的方法,由于这些单相流量计在理论研究和实际应用上都比较成熟,使得该方法在工业应用中更容易被接受。根据单相流量计组合的不同,该方法可以分为两个单相流量计组合法、单相流量计与密度计组合法和波动信号特征值法等。
应用节流法测量气液两相流量时,主要通过节流式仪表测量原理与两相流模型来实现。按照不同的假设条件,国内外研究者建立了均相流模型、分相流模型、Murdock关系式、Chisholm关系式、林宗虎关系式、漂移通量模型等数理模型。部分模型中的参数需要通过实验数据确定,当实验装置或应用条件不同时,模型中的参数也略有不同。
相关测量法以相关技术为基础测量两相流流量。该方法测量流速范围宽,不需要预先建立数理模型,无需通过实验确定参数,因此,相关流量计法为解决两相流量测量提供了一种强有力的技术手段。该技术的优点是可以采用不同的传感器来获得流体流动信号,再结合信号及数据处理技术,构成流体流量测量系统,实现非接触式测量。但相关流量测量技术目前仍存在一些问题需要进一步探讨,例如相关速度的物理意义、互相关函数峰值的确定方法、相关流量计的标定等。
科氏质量流量计能够直接测量流体的质量流量而无需测量流体的密度,同时,科氏质量流量计测量单相流体介质时,测量精度很高,可靠性好,维护费用低,因此在一些计量领域得到了很好的应用。然而,在很多应用领域中都存在着大量的气液两相流。目前研究者普遍认为在较高的含气率情况下,管子震动不稳定,很难得到稳定的信号。鉴于科氏质量流量计测量气液两相介质时测量误差较大,测量结果不能很好的满足石油石化及其它工业生产领域的要求,科氏质量流量计很少被考虑应用于气液两相介质流量的测量。
当待测介质含有气体时,介质密度将减小,同时,待测介质中含有的气体影响科氏质量流量计测量管的震动,从而导致流量计解算的介质流量偏小。可通过待测介质密度及流量测量值的变化估计待测介质的实际流量。
发明内容
本发明的目的在于根据科氏质量流量计测量气液两相流的流量和密度波动的统计特征,实现气液两相流分相流量的测量。一种基于质量流量计参数统计特征的两相流流量测量方法,具体内容包括以下步骤:
步骤A:建立流量和密度的统计特征与气液两相流含气率的关联关系;
步骤B:建立密度降与质量降的关联关系;
步骤C:应用流量和密度的统计特征与气液两相流含气率的关联关系估计气液两相流含气率α;
步骤D:应用密度降与质量降的关联关系估计管道内气相质量流量Mg和液相质量流量Ml。
上述基于质量流量计参数统计特征的两相流流量测量方法,在步骤A建立流量和密度的统计特征与气液两相流含气率的关联关系中,具体包括以下步骤:
A1、根据计算科氏质量流量计密度示值的平均值,根据计算科氏质量流量计密度示值的标准差,根据计算科氏质量流量计流量示值的平均值,根据计算科氏质量流量计流量示值的标准差,根据 计算科氏质量流量计流量示值的峰度,式中,N为采样点数,在实验中设置,ρi为科氏质量流量计测得的密度瞬时值,mi为科氏质量流量计测得的流量瞬时值;
A2、建立关联关系 式中,α为气液两相流含气率,在实验中根据 计算得到,ρl为液相密度,ρg为气相密度,ml为液相流量参考值,在实验中设置,mg为气相流量参考值,在实验中设置,λ1、λ2、λ3、λ4、γ1、γ2、γ3 和γ4为关系式系数,R1、R2、R3和R4为建立上述关联关系的回归系数,关系式系数和回归系数根据实验数据,由最小二乘法确定。
上述基于质量流量计参数统计特征的两相流流量测量方法,在步骤B建立密度降与质量降的关联关系中,具体包括以下步骤:
B1、建立密度降Δρ与科氏质量流量计质量降Δm之间的关联关系Δm=c×Δρ+d,式中, ρt为气液两相流混合密度,由ρt=ρg×α+(1-α)×ρl计算得到,α在步骤A2中得到,mapp为科氏质量流量计流量示值,由科氏质量流量计测量得到,c和d为关系式系数,根据实验数据,由最小二乘法确定;
B2、建立科氏质量流量计密度示值ρapp与气液两相流混合密度ρt之间的关系ρapp=a×ρt+b,式中,ρapp为科氏质量流量计密度示值,由科氏质量流量计测量得到,a和b为关系式系数,根据实验数据,由最小二乘法确定。
上述基于质量流量计参数统计特征的两相流流量测量方法,在步骤C应用流量和密度的统计特征与气液两相流含气率的关联关系估计气液两相流含气率α中,具体包括以下步骤:
上述基于质量流量计参数统计特征的两相流流量测量方法,在步骤D应用密度降与质量降的关联关系估计管道内气相质量流量Mg和液相质量流量Ml中,具体包括以下步骤:
D1、根据计算气液两相流混合密度ρt,式中,a和b在步骤B2中得到,ρapp 由科氏质量流量计在线测量得到;
D3、根据Δm=c×Δρ+d计算质量降Δm,式中,c和d在步骤B1中得到;
上述基于质量流量计参数统计特征的两相流流量测量方法,其特征在于根据实验数据,由最小二乘法离线确定a、b、c、d、λ1、λ2、λ3、λ4、γ1、γ2、γ3、γ4、R1、R2、R3和R4,并存储于计算机中。
本发明的有益效果在于:无需采用高效气液分离器进行气液分离,应用科氏质量流量计测得的瞬时质量流量和瞬时密度及其统计特征,结合存储于计算机中的关联关系及各关系式系数,实现气液两相流流量的测量。利用离线实验数据建立流量和密度的统计特征与气液两相流含气率的关联关系、建立密度降—质量降关联关系,利用已建立的各关联关系估计气液两相流含气率、气相质量流量和液相质量流量。该测量方法成本低、流量检测精度高、实时性好、易于实现。适用于两相流含气率及流量的测量。
附图说明
图1为科氏质量流量计测量含气率8%气液两相流时的密度瞬时波动图;
图2为科氏质量流量计测量含气率8%气液两相流时的流量瞬时波动图;
图3为根据气液两相流实验数据建立的测量值统计特征关联关系;
图4为密度降与科氏质量流量计质量降之间的关联关系图;
图5为科氏质量流量计密度示值与气液两相流混合密度之间的关联关系图;
图6为含气率预测结果图;
图7为气相质量流量预测结果图;
图8为液相质量流量预测结果图。
具体实施方式
本发明首先利用离线实验数据建立流量和密度的统计特征与气液两相流含气率的关联关系、并建立科氏质量流量计密度降—质量降关联关系,然后,应用科氏质量流量计测得的瞬时质量流量和瞬时密度及其统计特征,结合已建立的存储于计算机中的各关联关系及各关系式系数,估计气液两相流含气率、气相质量流量和液相质量流量。下面结合附图与具体实施方式对本发明作进一步详细描述。
本实施例实现气相体积流量从0.25m3/h到3.25m3/h,液相体积流量从3m3/h到10m3/h 的气液两相流的含气率、气相质量流量和液相质量流量的测量,其中,气相为空气,液相为实验柴油。
步骤A:建立流量和密度的统计特征与气液两相流含气率的关联关系。首先,根据 计算科氏质量流量计密度示值的平均值,根据计算科氏质量流量计密度示值的标准差,根据计算科氏质量流量计流量示值的平均值,根据计算科氏质量流量计流量示值的标准差,根据 计算科氏质量流量计流量示值的峰度,式中,N为采样点数,在实验中设置,ρi为科氏质量流量计测得的密度瞬时值,mi为科氏质量流量计测得的流量瞬时值。
在本实施例中,气液两相流的含气率变化范围为2.4%到52%,采样点数N设置为1000,图1和图2分别为科氏质量流量计测量含气率8%气液两相流时的密度瞬时波动和流量瞬时波动。
然后,建立关联关系 式中,α为气液两相流含气率,在实验中根据 计算得到,ρl为液相密度,ρg为气相密度,ml为液相流量参考值,在实验中设置,mg为气相流量参考值,在实验中设置,λ1、λ2、λ3、λ4、γ1、γ2、γ3和γ4为关系式系数,R1、R2、R3和R4为建立上述关联关系的回归系数,关系式系数和回归系数根据实验数据,由最小二乘法确定。
若关联关系为y=β0+βx,那么,按照最小二乘法原理,关系式系数计算式如下:
在本实施例中,根据含气率从2.4%到52%,气相体积流量从0.25m3/h到3.25m3/h,液相体积流量从3m3/h到10m3/h的气液两相流实验数据建立的上述关联关系如图3所示,其中,图3(a)为关联关系图3(b)为关联关系sdρ=λ2×α+γ2,图3(c)为关联关系sdm=λ3×α+γ3,图3(d)为关联关系kurm=λ4×α+γ4。根据实验数据,由最小二乘法确定的关系式系数为 回归系数为
步骤B:建立密度降与质量降的关联关系。首先,建立密度降Δρ与科氏质量流量计质量降Δm之间的关联关系Δm=c×Δρ+d,式中, ρt为气液两相流混合密度,由ρt=ρg×α+(1-α)×ρl计算得到,α在步骤A中得到,mapp为科氏质量流量计流量示值,由科氏质量流量计测量得到,c和d为关系式系数,根据实验数据,由最小二乘法确定。
在本实施例中,建立的密度降Δρ与科氏质量流量计质量降Δm之间的关联关系如图4所示。关系式系数为c=1.5124,d=2.4229。
然后,建立科氏质量流量计密度示值ρapp与气液两相流混合密度ρt之间的关系ρapp=a×ρt+b,式中,ρapp为科氏质量流量计密度示值,由科氏质量流量计测量得到,a和b为关系式系数,根据实验数据,由最小二乘法确定。
在本实施例中,建立的科氏质量流量计密度示值与气液两相流混合密度之间的关联关系如图5所示。关系式系数为a=0.6689,b=0.2815。
步骤C:应用流量和密度的统计特征与气液两相流含气率的关联关系估计气液两相流含气率α。首先,根据计算估计值α1,根据计算估计值α2,根据计算估计值α3,根据计算估计值α4。然后,计算式中,R1、R2、R3、R4分别为建立关联关系sdρ=λ2×α+γ2、sdm=λ3×α+γ3、kurm=λ4×α+γ4的回归系数。最后,根据 计算气液两相流含气率αmea。
含气率测量相对误差err计算式如下:
图6为含气率预测结果。
步骤D:应用密度降与质量降的关联关系估计管道内气相质量流量Mg和液相质量流量Ml。首先,根据计算气液两相流混合密度ρt,式中,a和b在步骤B中得到,ρapp由科氏质量流量计在线测量得到。然后,根据计算密度降Δρ,根据Δm=c×Δρ+d计算质量降Δm,式中,c和d在步骤B中得到。最后,根据计算管道内液相质量流量Ml,根据计算管道内气相质量流量Mg,式中,mapp由科氏质量流量计在线测量得到,αmea在步骤C中得到。
图7和图8分别为气相质量流量和液相质量流量的预测结果。
Claims (6)
1.一种基于质量流量计参数统计特征的两相流流量测量方法,其特征包括如下步骤:
步骤A:建立流量和密度的统计特征与气液两相流含气率的关联关系;
步骤B:建立密度降与质量降的关联关系;
步骤C:应用流量和密度的统计特征与气液两相流含气率的关联关系估计气液两相流含气率α;
步骤D:应用密度降与质量降的关联关系估计管道内气相质量流量Mg和液相质量流量Ml。
2.根据权利要求1所述的一种基于质量流量计参数统计特征的两相流流量测量方法,其特征在于上述步骤A建立流量和密度的统计特征与气液两相流含气率的关联关系中,具体包括以下步骤:
A1、根据计算科氏质量流量计密度示值的平均值,根据计算科氏质量流量计密度示值的标准差,根据计算科氏质量流量计流量示值的平均值,根据计算科氏质量流量计流量示值的标准差,根据 计算科氏质量流量计流量示值的峰度,式中,N为采样点数,在实验中设置,ρi为科氏质量流量计测得的密度瞬时值,mi为科氏质量流量计测得的流量瞬时值;
3.根据权利要求1所述的一种基于质量流量计参数统计特征的两相流流量测量方法,其特征在于上述步骤B建立密度降与质量降的关联关系中,具体包括以下步骤:
B1、建立密度降Δρ与科氏质量流量计质量降Δm之间的关联关系Δm=c×Δρ+d,式中, ρt为气液两相流混合密度,由ρt=ρg×α+(1-α)×ρl计算得到,α在步骤A2中得到,mapp为科氏质量流量计流量示值,由科氏质量流量计测量得到,c和d为关系式系数,根据实验数据,由最小二乘法确定;
B2、建立科氏质量流量计密度示值ρapp与气液两相流混合密度ρt之间的关系ρapp=a×ρt+b,式中,ρapp为科氏质量流量计密度示值,由科氏质量流量计测量得到,a和b为关系式系数,根据实验数据,由最小二乘法确定。
6.根据权利要求1所述的一种基于质量流量计参数统计特征的两相流流量测量方法,其特征在于根据实验数据,由最小二乘法离线确定a、b、c、d、λ1、λ2、λ3、λ4、γ1、γ2、γ3、γ4、R1、R2、R3和R4,并存储于计算机中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310670292.5A CN103674139B (zh) | 2013-12-11 | 2013-12-11 | 基于质量流量计参数统计特征的两相流流量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310670292.5A CN103674139B (zh) | 2013-12-11 | 2013-12-11 | 基于质量流量计参数统计特征的两相流流量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103674139A true CN103674139A (zh) | 2014-03-26 |
CN103674139B CN103674139B (zh) | 2016-04-06 |
Family
ID=50312396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310670292.5A Expired - Fee Related CN103674139B (zh) | 2013-12-11 | 2013-12-11 | 基于质量流量计参数统计特征的两相流流量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103674139B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108412481A (zh) * | 2018-03-05 | 2018-08-17 | 中国石油天然气股份有限公司 | 超声波多普勒多相流分相流量测井资料解释方法和装置 |
CN115096383A (zh) * | 2022-07-15 | 2022-09-23 | 海默新宸水下技术(上海)有限公司 | 基于等效密度测算多相流中气相流量的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10281846A (ja) * | 1997-04-09 | 1998-10-23 | Oval Corp | コリオリ流量計を利用したパターン認識法による多相流量計 |
CN1351709A (zh) * | 1999-03-26 | 2002-05-29 | 微动公司 | 具有统计学优选技术的流量计校准系统 |
CN101819056A (zh) * | 2005-09-19 | 2010-09-01 | 微动公司 | 用于流量计的校验诊断的仪表电子器件和方法 |
CN102346058A (zh) * | 2010-07-29 | 2012-02-08 | 中国石油大学(华东) | 科氏质量流量计测量夹气液体流量的模型法 |
CN102589628A (zh) * | 2003-02-10 | 2012-07-18 | 因万西斯系统股份有限公司 | 多相哥氏流量计 |
-
2013
- 2013-12-11 CN CN201310670292.5A patent/CN103674139B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10281846A (ja) * | 1997-04-09 | 1998-10-23 | Oval Corp | コリオリ流量計を利用したパターン認識法による多相流量計 |
CN1351709A (zh) * | 1999-03-26 | 2002-05-29 | 微动公司 | 具有统计学优选技术的流量计校准系统 |
CN102589628A (zh) * | 2003-02-10 | 2012-07-18 | 因万西斯系统股份有限公司 | 多相哥氏流量计 |
CN101819056A (zh) * | 2005-09-19 | 2010-09-01 | 微动公司 | 用于流量计的校验诊断的仪表电子器件和方法 |
CN102346058A (zh) * | 2010-07-29 | 2012-02-08 | 中国石油大学(华东) | 科氏质量流量计测量夹气液体流量的模型法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108412481A (zh) * | 2018-03-05 | 2018-08-17 | 中国石油天然气股份有限公司 | 超声波多普勒多相流分相流量测井资料解释方法和装置 |
CN115096383A (zh) * | 2022-07-15 | 2022-09-23 | 海默新宸水下技术(上海)有限公司 | 基于等效密度测算多相流中气相流量的方法 |
CN115096383B (zh) * | 2022-07-15 | 2022-11-22 | 海默新宸水下技术(上海)有限公司 | 基于等效密度测算多相流中气相流量的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103674139B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2760930C (en) | Multi-phase fluid measurement apparatus and method | |
CN102759383B (zh) | 基于单节流元件的气液两相流气相流量在线测量方法及装置 | |
CN100437046C (zh) | 基于截面测量的气液两相流测量方法及装置 | |
CN101178347B (zh) | 狭缝文丘里节流装置及气液两相流测量系统 | |
GB2432425A (en) | Isokinetic sampling method and apparatus for multiphase flow from subterranean wells | |
CN201241692Y (zh) | 湿天然气多相流量计 | |
CN101900589B (zh) | 基于质量流量计的夹气液体流量测量方法 | |
CN107843297B (zh) | 基于v锥的低含气率气液两相流液相流量在线测量装置及方法 | |
Hua et al. | Wet gas meter based on the vortex precession frequency and differential pressure combination of swirlmeter | |
CN102346058B (zh) | 科氏质量流量计测量夹气液体流量的模型法 | |
CN201837418U (zh) | 高精度宽量程一体化节流装置 | |
CN101666770A (zh) | 低含气率原油测量装置及方法 | |
CN105486358A (zh) | 基于文丘里管双差压的气液两相流参数测量方法 | |
CN102207398B (zh) | 用于燃油终端结算的超声波流量测量的装置和方法 | |
CN103063320B (zh) | 一种高温高压复杂天然气节流温度和气液比的测定方法 | |
US7447599B2 (en) | Method and system for generating an uncertainty value | |
CN103674139B (zh) | 基于质量流量计参数统计特征的两相流流量测量方法 | |
CN201032457Y (zh) | 基于标准文丘里管与涡街流量计的低含气率气液两相流测量装置 | |
CN201327400Y (zh) | 湿蒸汽两相流流量计 | |
CN202748069U (zh) | 湿气流量测量装置 | |
Wang et al. | Mass flow measurement of two-phase carbon dioxide using coriolis flowmeters | |
RU114338U1 (ru) | Устройство для измерения продукции нефтедобывающих скважин | |
CN114547892A (zh) | 基于液膜流动参数建模的涡街湿气分相流量测量方法 | |
CN103353319A (zh) | 基于直通型气体超声波流量计的湿气流量测量方法 | |
CN103323066A (zh) | 一种低含液率气液两相流测量方法及测量系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20160406 Termination date: 20161211 |