CN106643934B - 天然气音速计算方法 - Google Patents
天然气音速计算方法 Download PDFInfo
- Publication number
- CN106643934B CN106643934B CN201611005642.6A CN201611005642A CN106643934B CN 106643934 B CN106643934 B CN 106643934B CN 201611005642 A CN201611005642 A CN 201611005642A CN 106643934 B CN106643934 B CN 106643934B
- Authority
- CN
- China
- Prior art keywords
- natural gas
- makings
- sound
- gas 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
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/66—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 measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种天然气音速计算方法,天然气音速为v,其特征在于:计算天然气音速简化公式为:v=s1t2+s2pt+s3p+s4t+c,其中p是待测天然气的气质压力,t是待测天然气的气质温度,s1、s2是常数,s3=Ac+B,s4=A1c+B1,其中A、B、A1、B1均为常数,c是由待测天然气气质决定的参量。
Description
技术领域
本发明涉及天然气流量计量,尤其是一种天然气音速计算方法。
背景技术
随着天然气工业的飞速发展,特别是输气管网的大规模建设,用于贸易计量的天然气流量仪表日益增多,其工作压力不断提高,流量范围也不断增大。天然气流量计量是多参数、多组分气体的连续测量,其量值测量具有不可回复性,其测量准确度受众多因素影响。
现有的计算天然气音速方法主要是根据美国燃气协会发表的AGA10号报告中提出的计算公式来计算的,该计算方法较为繁琐,设计多次求偏导数以及积分计算,同时由于气体流量计的成本控制,流量计内置的单片机核心计算能力一般,面对高强度、高频率的计算要求时效率低下。
发明内容
本发明所要解决的技术问题是针对上述现有技术存在的问题,提供一种计算简单而高效的天然气音速计算方法。
本发明解决上述技术问题所采用的技术方案为:一种天然气音速计算方法,其特征在于:1)确定管道内的天然气气质;2)根据公式v=s1t2+s2pt+s3p+s4t+c,得到管道内采样点的天然气音速v,其中p是待测天然气的气质压力,t是待测天然气的气质温度,s1、s2是常数,s3=Ac+B,s4=A1c+B1,其中A、B、A1、B1均为常数,c是由待测天然气气质决定的参量,并且t的范围为-20℃~60℃,p的范围为50KPa~2000KPa。
为了使得天然气音速的计算误差最小,s1=-1.128×10-3,s2=5.8398×10-5,s3=1.26×10-5c-1.018×10-2,s4=1.814×10-3c-4.139×10-2,由此v=-1.128×10-3t2+5.8398×10-5pt+(1.26×10-5c-1.018×10-2)p+(1.814×10-3c-4.139×10-2)t+c。
系数c通过如下步骤得到:
1)首先利用现有的天然气音速计算公式计算出待测天然气气质在选定的标准气质温度t0,标准气质压力p0下的标准天然气音速值v0;
2)将计算得到的标准天然气音速值v0、对应的标准气质压力p0与标准气质温度t0代入系数c的计算公式而得到系数c:
与现有技术相比,本发明的优点在于:适用于家用燃气流量计算,能精确、高效的计算出当下的天然气音速,运用到当下主流的气体流量计当中,能提高工作效率。
附图说明
图1为天然气音速计算流程图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
一种天然气音速计算方法,考虑到城镇燃气管道的真实工作环境,气质温度的主要变化范围为-20℃至60℃,气质压力的主要变化范围为50KPa至2000KPa的环境下,因此,本发明的计算方法所涉及到的气质温度与气质压力均在上述所提及的变化范围之内。本发明的计算方法可以克服AGA10号报告中提出的天然气音速计算公式存在的求解繁琐的问题,并以现行的AGA10报告提出的天然气音速计算方式为基准,将本发明计算结果的相对误差控制在0.5%以内。
具体地,天然气音速v是关于气质温度与气质压力的二次多项式:
v=s1t2+s2pt+s3p+s4t+c
其中系数s1、s2是与气质无关的常数项,系数c是与气质有关的参量,p是气质压力(可测得),t是气质温度(可测得),而系数s3、系数s4均可用与系数c有关的线性关系式表示,上述公式中只有一个待定标量c。s3=Ac+B,s4=A1c+B1,其中A、B、A1、B1均为常数。
为了使得天然气音速的计算误差最小,系数s1为常数项,且值为s1=-1.128×10-3;系数s2为常数项,且值为s2=5.8398×10-5;系数s3可以用系数c的线性表达式表示,该式为:s3=1.26×10-5c-1.018×10-2;系数s4可以用与系数c有关的线性关系式来表示,系数s4的计算公式为:s4=1.814×10-3c-4.139×10-2,只与待标定量c有关。
而系数c的确定可通过以下方式:在气质组成成分已知的条件下,利用美国燃气协会10号报告AGA10提出的音速计算公式计算在任一气质压力p0与任一气质温度t0下的天然气音速大小v0,将计算得到的v0以及对应的p0、t0代入下式中便可得到c的值:
因此,本发明公式的具体形式为:
v=-1.128×10-3t2+5.8398×10-5pt+(1.26×10-5c-1.018×10-2)p+(1.814×10- 3c-4.139×10-2)t+c
实施例1:
以美国燃气协会10号报告所提到的气质GulfCoast例,该气质的主要成分摩尔百分比见表1。
表1各种气质主要成分摩尔百分比
同时,得到针对于该气质的音速计算公式的流程可见图1:
1)首先利用美国燃气协会AGA10号报告提出的天然气音速计算公式计算出在标准气质温度t0为20℃,标准气质压力p0为1000KPa下的标准天然气音速值v0(计算v0时,选取其他的气质温度与气质压力作为标准参考值也均可),得到标准天然气音速值v0=430.104514m/s;在此步骤中,也可以用其他天然气音速计算方法得到标准天然气音速值v0;
2)将计算得到的标准天然气音速值v0以及对应的标准气质压力p0与标准气质温度t0代入系数c的计算公式:
从而计算得到系数c=419.1554;
3)确定了系数c的值后,将系数c代入本发明公式中修正音速公式:
v=-1.128×10-3t2+5.8398×10-5pt+(1.26×10-5c-1.018×10-2)p+(1.814×10- 3c-4.139×10-2)t+c
便可以得到针对于气质Gulf Coast的音速简化计算公式,即:
v=-1.128×10-3t2+5.8398×10-5pt-4.898×10-3p+0.7190t+419.1554
此时,为了验证得到关于气质Gulf Coast的音速简化计算公式的可行性,可进行相关的误差分析。在选取一系列气质温度与气质压力情况下,以美国燃气协会10号报告AGA10提出的计算公式为基准,计算本发明的相对误差大小。具体情况可见表2。
表2气质GulfCoast计算比较
实施例2:为了进一步说明本发明的可行性,这里又选取了3种美国燃气协会10号报告所提到的气质,即Amarillo、High N2、High CO2,这三种气质的具体成分摩尔百分比可见表1。
首先对3种天然气的音速值进行采样,考虑到城镇燃气管道的真实工作环境,故取气质温度从-20℃至60℃变化,5℃为步长;气质气压取从50KPa到2000KPa变化,取其中50KPa至300KPa以50KPa为步长,300KPa至2000KPa以200KPa为步长,利用AGA10报告中提出的音速计算公式对3种气质进行不同气质温度与气质压力下的音速值采样。由于数据采样点比较多,表3~表5分别截取了3种气质在一部分采样点上的音速值。
同时,通过图1所示的步骤,利用AGA10报告提出的音速计算公式计算出20摄氏度、1000Kpa下不同天然气的音速值v0,从而确定对应于不同气质的音速计算公式。根据得到的公式分别计算上述3种气质在各个采样点上的音速值,同样由于数据采样点较多,表3~表5分别截取了3种气质在一部分采样点上的音速值。
以AGA10报告提出的天然气音速计算方法为基准,得到通过AGA10报告计算出的音速值v1,并计算出本发明提出的天然气音速简化计算公式的计算结果v2。并根据相对误差计算公式计算得到对应的相对误差大小。
表3气质Amarillo计算比较
表4气质High N2计算比较
表5气质High CO2计算比较
上述的表3~表5分别截取了3种气质在一部分采样点上的相对误差大小,由此可见相对误差均控制在了0.5%以下,证明本发明的计算公式具有较强的可行性。
上文中例举了一组系数s1、s2、s3、s4,而系数s1、s2、s3、s4还可以有其他替代方案,满足相对误差均控制在0.5%以下的要求:如
s1=-1.130×10-3
s2=5.607×10-5
s3=-6.800×10-6c-2.080×10-3
s4=1.752×10-3c-1.553×10-2
v=-1.130×10-3t2+5.607×10-5pt+(-6.800×10-6c-2.080×10-3)p+(1.752×10- 3c-1.553×10-2)t+c
又如:
s1=-1.142×10-3
s2=5.864×10-5
s3=2.540×10-5c-1.549×10-2
s4=1.864×10-3c-6.21×10-2
v=-1.142×10-3t2+5.864×10-5pt+(2.540×10-5c-1.549×10-2)p+(1.864×10- 3c-6.21×10-2)t+
又如:
s1=-1.121×10-3
s2=5.828×10-5
s3=1.450×10-5c-1.097×10-2
s4=1.825×10-3c-4.584×10-2
v=-1.121×10-3t2+5.828×10-5pt+(1.450×10-5c-1.097×10-2)p+(1.825×10- 3c-4.584×10-2)t+c。
Claims (1)
1.一种天然气音速计算方法,其特征在于:
1)确定管道内的天然气气质;
2)根据公式v=s1t2+s2pt+s3p+s4t+c,得到管道内采样点的天然气音速v,其中p是待测天然气的气质压力,t是待测天然气的气质温度,t的范围为-20℃~60℃,p的范围为50KPa~2000KPa,s1=-1.128×10-3,s2=5.8398×10-5,s3=1.26×10-5c-1.018×10-2,s4=1.814×10-3c-4.139×10-2,由此v=-1.128×10-3t2+5.8398×10-5pt+(1.26×10-5c-1.018×10-2)p+(1.814×10-3c-4.139×10-2)t+c,c是由待测天然气气质决定的参量,通过如下步骤得到:
1)首先利用现有的天然气音速计算公式计算出待测天然气气质在选定的标准气质温度t0,标准气质压力p0下的标准天然气音速值v0;
2)将计算得到的标准天然气音速值v0、对应的标准气质压力p0与标准气质温度t0代入系数c的计算公式而得到系数c:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611005642.6A CN106643934B (zh) | 2016-11-15 | 2016-11-15 | 天然气音速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611005642.6A CN106643934B (zh) | 2016-11-15 | 2016-11-15 | 天然气音速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106643934A CN106643934A (zh) | 2017-05-10 |
CN106643934B true CN106643934B (zh) | 2019-03-29 |
Family
ID=58806026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611005642.6A Active CN106643934B (zh) | 2016-11-15 | 2016-11-15 | 天然气音速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106643934B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109990855A (zh) * | 2019-04-09 | 2019-07-09 | 宁波大学 | 一种超声波气体流量计标况音速计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1077796A (zh) * | 1991-12-23 | 1993-10-27 | 蒂特达斯屈托拉仪器股份公司 | 监测气体流量,特别是天然气流量的方法和装置 |
EP2887027A1 (en) * | 2013-12-18 | 2015-06-24 | Siemens Energy, Inc. | Active measurement of gas flow velocity or simultaneous measurement of velocity and temperature, including in gas turbine combustors |
CN105910696A (zh) * | 2016-06-24 | 2016-08-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种可压缩气体声速测量方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS564019A (en) * | 1979-06-25 | 1981-01-16 | Toshiba Corp | Acoustic velocity measuring device |
-
2016
- 2016-11-15 CN CN201611005642.6A patent/CN106643934B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1077796A (zh) * | 1991-12-23 | 1993-10-27 | 蒂特达斯屈托拉仪器股份公司 | 监测气体流量,特别是天然气流量的方法和装置 |
EP2887027A1 (en) * | 2013-12-18 | 2015-06-24 | Siemens Energy, Inc. | Active measurement of gas flow velocity or simultaneous measurement of velocity and temperature, including in gas turbine combustors |
CN105910696A (zh) * | 2016-06-24 | 2016-08-31 | 中国航空工业集团公司西安飞机设计研究所 | 一种可压缩气体声速测量方法及装置 |
Non-Patent Citations (1)
Title |
---|
《渐缩形管道中亚音速气流的物态分析》;张立红等;《山东建筑工程学院学报》;20050330(第01期);第85-87页 |
Also Published As
Publication number | Publication date |
---|---|
CN106643934A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110514391B (zh) | 一种五孔压力探针测量三维流场不确定度评定方法 | |
EP2711678B1 (en) | Correction of pressure signals measured during supersonic wind tunnel testing | |
CN110514390B (zh) | 一种三孔压力探针测量二维流场不确定度评定方法 | |
CN204085645U (zh) | 带有自校准结构的气体流量标准装置 | |
CN201402160Y (zh) | 具有零点标定功能的皮托管流量计 | |
CN105571653B (zh) | 磨煤机风量测量方法和系统 | |
EP1102967A1 (en) | Method and apparatus to determine a property of a multiphase fluid | |
EP2124019A3 (en) | Real time error determination for inertial instruments | |
CN102226904A (zh) | 电站大型锅炉一次风风量的软测量方法 | |
CN105045973B (zh) | 一种变弧长自适应采样方法 | |
CN109840671A (zh) | 运行改进效果计算设备、运行改进效果计算方法和记录介质 | |
CN107631773A (zh) | 用于运行流量测量仪器的方法和流量测量仪器 | |
CN106643934B (zh) | 天然气音速计算方法 | |
CN102879061A (zh) | 一种基于拟合方程的水表误差校正方法 | |
RU2009117712A (ru) | Способ мониторинга коррозии трубопровода и устройство для его осуществления | |
CN107289985A (zh) | 在线标定和补偿的动态测量系统及方法 | |
CN103605323B (zh) | 化工生产的离散控制方法及装置 | |
CN106918377B (zh) | 用于虚拟流量计的校准装置、灵敏度确定模块及相应方法 | |
CN106093905B (zh) | 基于apFFT相位差的精密测量雷达频偏分类测频方法 | |
CN106525180B (zh) | 天然气压缩因子计算方法 | |
CN105181049B (zh) | 管道内流体流量的测量方法、测量系统及辅助测量系统 | |
CN203534679U (zh) | 孔板流量计的校正系统 | |
CN109635421B (zh) | 一种基于威布尔模型的一般压力表检测周期动态优化方法 | |
CN109990855A (zh) | 一种超声波气体流量计标况音速计算方法 | |
CN109632029B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221212 Address after: 200120 floor 1-2, north end of Building 1, No. 51, Mindong Road, Pudong New Area, Shanghai Patentee after: DEWEN APPARATUS INSTRUMENT (SHANGHAI) Co.,Ltd. Address before: 315211, Fenghua Road, Jiangbei District, Zhejiang, Ningbo 818 Patentee before: Ningbo University |