[go: up one dir, main page]

CN106525180B - Method of Calculating Compressibility Factors of Natural Gas - Google Patents

Method of Calculating Compressibility Factors of Natural Gas Download PDF

Info

Publication number
CN106525180B
CN106525180B CN201611085801.8A CN201611085801A CN106525180B CN 106525180 B CN106525180 B CN 106525180B CN 201611085801 A CN201611085801 A CN 201611085801A CN 106525180 B CN106525180 B CN 106525180B
Authority
CN
China
Prior art keywords
gas
natural gas
makings
temperature
pressure
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
CN201611085801.8A
Other languages
Chinese (zh)
Other versions
CN106525180A (en
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.)
Dewen Apparatus Instrument Shanghai Co ltd
Original Assignee
Ningbo University
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 Ningbo University filed Critical Ningbo University
Priority to CN201611085801.8A priority Critical patent/CN106525180B/en
Publication of CN106525180A publication Critical patent/CN106525180A/en
Application granted granted Critical
Publication of CN106525180B publication Critical patent/CN106525180B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F15/00—Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
    • G01F15/02—Compensating or correcting for variations in pressure, density or temperature
    • G01F15/04—Compensating or correcting for variations in pressure, density or temperature of gases to be measured
    • G01F15/043—Compensating or correcting for variations in pressure, density or temperature of gases to be measured using electrical means
    • G01F15/046—Compensating or correcting for variations in pressure, density or temperature of gases to be measured using electrical means involving digital counting

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Measuring Volume Flow (AREA)

Abstract

The invention discloses a kind of Method of Calculating Compressibility Factors of Natural Gas, characterized by the following steps: 1) determines the natural gas quality in pipeline；2) sampling obtains the makings temperature t and makings pressure p of the natural gas quality in pipeline, and wherein the range of t is 0 DEG C~60 DEG C, and the range of p is 50KPa~4000KPa；3) according to formula z=s1t2+s2pt+s3t+s4P+c obtains gas deviation factor z, wherein s4It is the parameter determined by the natural gas quality, c is constant, s1=A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is constant.

Description

Method of Calculating Compressibility Factors of Natural Gas
Technical field
The present invention relates to natural gas flow measuring, especially a kind of Method of Calculating Compressibility Factors of Natural Gas.
Background technique
Day with the Large scale construction of the rapid development of gas industry, especially gas line network, for Trade Measures Right throughput instrument is increasing, and operating pressure is continuously improved, and range of flow also constantly increases.Natural gas flow measuring is more The continuous measurement of parameter, multicomponent gas, magnitude measurement have not reversibility, and accuracy of measurement is by many factors shadow It rings.
Only under low pressure, high temperature real gas just can approximation be counted as perfect gas, due to real gas and ideal The difference of gas, so that be difficult to evaluate to the accuracy and reliability in gas flow measurement, especially low temperature, pressure piping gas The measurement of body flow.In this case, the tested makings in pipeline cannot be described with The Ideal-Gas Equation And processing.
This deviation of real gas and perfect gas can illustrate that this ratio is claimed using the ratio of PV and RT For compressibility factor, compressibility factor is defined as: authorized pressure and at a temperature of, the volume of any mass gas is with the gas in phase The ratio of the gas volume calculated under the conditions of by perfect gas law.Compressibility factor is indicated with zed, if Z > 1, The volume of real gas is greater than the volume of the perfect gas of equivalent at that the same temperature and pressure；If Z < 1, identical Temperature and pressure under real gas volume be less than equivalent perfect gas volume.
Compressibility factor is calculated there are many method, the AGA8 report that American Gas association delivers is calculated by gas component The compressibility factor of natural gas and other related hydrocarbons gases, AGA8-92DC method are exactly one of very important method, should Calculation method is relatively complicated, and design repeatedly asks partial derivative and integral calculation, simultaneously because the cost control of gas flowmeter, stream Monolithic microcomputer kernel computing capability built in meter is general, inefficiency when requiring in face of high-intensitive, high-frequency calculating.
Summary of the invention
The technical problem to be solved by the present invention is in view of the above-mentioned problems of the prior art, providing a kind of simple, fast Prompt and accurate Method of Calculating Compressibility Factors of Natural Gas.
The technical scheme of the invention to solve the technical problem is: a kind of Method of Calculating Compressibility Factors of Natural Gas, It is characterized by comprising following steps:
1) natural gas quality in pipeline is determined；
2) sampling obtains the makings temperature t and makings pressure p of the natural gas quality in pipeline, and wherein the range of t is 0 DEG C~60 DEG C, the range of p is 50KPa~4000KPa；
3) according to formula z=s1t2+s2pt+s3t+s4P+c obtains gas deviation factor z, wherein s4It is by natural gas gas The parameter that matter determines, c is constant, s1=A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is Constant.
In order to enable the calculating error of gas deviation factor is minimum, s1=9.707 × 10-2s4-2.944×10-7, s2 =-8.115 × 10-3s4+3.136×10-8, s3=-5.007s4+3.214×10-5, thus (9.707 × 10 z=-2s4-2.944 ×10-7)t2+(-8.115×10-3s4+3.136×10-8)pt+(-5.007s4+3.214×10-5)t+s4p+0.9986。
Coefficient s4It obtains as follows:
1) the natural gas quality temperature is calculated in selected mark first with existing gas deviation factor formula Quasi- makings temperature, standard makings pressure p0Under standard gas deviation factor value z0；
2) the standard gas deviation factor value z that will be calculated0And corresponding standard makings pressure p0With standard makings Temperature t0Substitute into coefficient s4Calculation formula and obtain:
Compared with the prior art, the advantages of the present invention are as follows: it is suitable for domestic gas flow rate calculation, it can be accurate, efficient Gas deviation factor instantly is calculated, applies in the gas flowmeter of mainstream instantly, working efficiency can be improved.
Detailed description of the invention
Fig. 1 is the flow chart of calculation method of the invention.
Specific embodiment
The present invention will be described in further detail below with reference to the embodiments of the drawings.
A kind of Method of Calculating Compressibility Factors of Natural Gas, it is contemplated that the Main change model of makings temperature is chosen in actual application Enclosing is 0 DEG C to 60 DEG C, in the environment of the Main change range of makings pressure is 50KPa to 4000KPa, provides a kind of easy day Right gas compressibility factor calculation formula overcomes the problems, such as that the existing solution of the AGA8-92DC equation proposed in No. AGA8 report is cumbersome, And on the basis of existing AGA8-92DC equation, the relative error of formula calculated result of the present invention is controlled within 0.5%.
Specifically, gas deviation factor z is the quadratic polynomial about makings temperature and makings pressure, wherein containing only One coefficient s to be calibrated4, the form of formula are as follows: z=s1t2+s2pt+s3t+s4p+c.Wherein p is the makings pressure of natural gas to be measured Power, t are the makings temperature of natural gas to be measured, and coefficient c is the constant term unrelated with natural gas quality to be measured, coefficient s4Be with it is to be measured The related amount of natural gas quality.It further illustrates, coefficient s1、s2、s3Available and coefficient s4Related linear relation expression, s1 =A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is constant.
In order to enable the calculating error of gas deviation factor is minimum, and in the present embodiment, c=0.9986, s1=9.707 ×10-2s4-2.944×10-7, s2=-8.115 × 10-3s4+3.136×10-8, s3=-5.007s4+3.214×10-5。
Coefficient s4Determination can be in the following manner: under the conditions of known to the makings constituent, assisted using American Gas The AGA8-92DC equation calculations that No. 8 report AGA8 of meeting are proposed are in ethos pressure p0With ethos temperature t0Under natural gas Compressibility factor size z0, the z that will be calculated0And corresponding p0、t0S can be obtained by substituting into following formula4The value of coefficient.
Therefore, the concrete form of calculation method of the present invention are as follows:
Z=(9.707 × 10-2s4-2.944×10-7)t2+(-8.115×10-3s4+3.136×10-8)pt+(-5.007s4 +3.214×10-5)t+s4p+0.9986
Embodiment one
By taking the gas sample 1 that No. 8 reports of American Gas association are previously mentioned as an example, the main component molar percentage of the gas sample is shown in attached Table 1.
The various makings main component molar percentages (%) of table 1
Meanwhile obtaining the visible Fig. 1 of process for the compressibility factor calculation formula for being directed to the gas sample 1:
1) go out first with the AGA8-92DC equation calculation that No. 8 reports of American Gas association propose in makings temperature t0For 10 DEG C, makings pressure p0For the standard gas deviation factor value z under 1000KPa0.(calculate z0When, choose other makings temperature Degree with makings pressure as standard reference value also), obtain gas deviation factor value z0=0.97744；In this step, Standard gas deviation factor value z can also be obtained with other Method of Calculating Compressibility Factors of Natural Gas0；
2) the gas deviation factor value z that will be calculated0And corresponding makings pressure p0With makings temperature t0Substitute into system Number s4Calculation formula:
To which coefficient s be calculated4=-2.478 × 10-5；
3) coefficient s has been determined4Value after, by coefficient s4It substitutes into and corrects compressibility factor calculation formula in formula of the present invention:
Z=(9.707 × 10-2s4-2.944×10-7)t2+(-8.115×10-3s4+3.136×10-8)pt+(-5.007s4 +3.214×10-5)t+s4p+0.9986
The compressibility factor calculation formula for being directed to gas sample 1 can be obtained as a result, it may be assumed that
Z=-2.699 × 10-6t2+2.324×10-7pt+1.562×10-4t-2.478×10-5p+0.9986
At this point, relevant error can be carried out in order to verify to obtain the feasibility of the compressibility factor calculation formula about gas sample 1 Analysis.In the case where choosing a series of makings temperature and makings pressure condition, proposed with No. 8 report AGA8 of American Gas association On the basis of the calculated result of AGA8-92DC equation, the relative error size of formula of the present invention is calculated.Concrete condition can be shown in Table 2.
The calculating of 2 gas sample 1 of table is compared
Embodiment two
Feasibility in order to further illustrate the present invention has chosen No. 8 reports of 5 kinds of American Gas associations again here and is previously mentioned Gas sample, the specific Component molar percentage of this 5 kinds of gas samples can be shown in Table 1.
The compressibility factor value of different gas samples is sampled first, it is contemplated that the practical application of compressibility factor, therefore take makings Temperature DEG C variation from 0 DEG C to 60,5 DEG C are step-length；Makings pressure takes to be changed from 50KPa to 2000KPa, takes wherein that 50KPa is extremely 300KPa using 50KPa as step-length, 300KPa to 2000KPa using 200KPa as step-length, using AGA8 report in the AGA8- that proposes 92DC equation carries out the compressibility factor value at different makings temperature and makings pressure to 5 kinds of gas samples and samples.Due to data sampling point Compare more, table 3 to table 7 has intercepted compressibility factor value of 5 kinds of gas samples on a part of sampled point respectively.
Meanwhile the step of passing through Fig. 1, (we select the determining compressibility factor calculation formula for corresponding to different gas samples here Ethos pressure is 1000KPa, and ethos temperature is 10 DEG C).Above-mentioned 5 kinds of gas samples are calculated each according to obtained formula Compressibility factor value on sampled point, also due to data sampling point is more, table 3 to table 7 has intercepted 5 kinds of gas samples in a part respectively Compressibility factor value on sampled point.
On the basis of the AGA8-92DC equation that existing AGA8 report proposes, the compressibility factor value by the equation is obtained Calculated result z1, and calculate the calculated result z of gas deviation factor calculation formula proposed by the present invention2.And according to opposite mistake Poor calculation formulaCorresponding relative error size is calculated.
Table 3 to table 7 has intercepted relative error size of 5 kinds of gas samples on a part of sampled point respectively, it is seen that relative error Control is 0.5% or less, it was demonstrated that invention formula has stronger feasibility.
The calculating of 3 gas sample 2 of table is compared
The calculating of 4 gas sample 3 of table is compared
The calculating of 5 gas sample 4 of table is compared
The calculating of 6 gas sample 5 of table is compared
The calculating of 7 gas sample 6 of table is compared
Above illustrate one group of coefficient s1、s2、s3, and coefficient s1、s2、s3There can also be other alternative solutions, meet phase Error is controlled in 0.5% requirement below: such as
s1=9.711 × 10-2s4-2.931×10-7
s2=-8.125 × 10-3s4+3.105×10-8
s3=-5.018s4+3.178×10-5
Z=(9.711 × 10-2s4-2.931×10-7)t2+(-8.125×10-3s4+3.105×10-8)pt+(-5.018s4 +3.178×10-5)t+s4p+0.9986
For another example:
s1=9.7 × 10-2s4-2.955×10-7
s2=-8.114 × 10-3s4+3.139×10-8
s3=-4.998s4+3.229×10-5
Z=(9.7 × 10-2s4-2.955×10-7)t2+(-8.114×10-3s4+3.139×10-8)pt+(-4.998s4+ 3.229×10-5)t+s4p+0.9986
For another example:
s1=9.696 × 10-2s4-3.056×10-7
s2=-8.11 × 10-3s4+3.184×10-8
s3=-5.000s4+3.283×10-5
Z=(9.696 × 10-2s4-3.056×10-7)t2+(-8.11×10-3s4+3.184×10-8)pt+(-5.000s4+ 3.283×10-5)t+s4p+0.9986

Claims (1)

1.一种天然气压缩因子计算方法，其特征在于：包括如下步骤：1. a natural gas compression factor calculation method, is characterized in that: comprise the steps: 1)确定管道内的天然气气质；1) Determine the natural gas quality in the pipeline; 2)在管道内采样得到所述天然气气质的气质温度t和气质压力p，其中t的范围为0℃～60℃，p的范围为50KPa～4000KPa；2) Sampling in the pipeline to obtain the gas temperature t and gas pressure p of the natural gas gas, where t ranges from 0°C to 60°C and p ranges from 50KPa to 4000KPa; 3)根据公式z＝s1t2+s2pt+s3t+s4p+c得到天然气压缩因子z，其中s4是由所述天然气气质决定的参量，s1＝9.707×10-2s4-2.944×10-7，s2＝-8.115×10-3s4+3.136×10-8，s3＝-5.007s4+3.214×10-5，c＝0.9986，由此3) According to the formula z=s 1 t 2 +s 2 pt+s 3 t+s 4 p+c, the natural gas compression factor z is obtained, where s 4 is a parameter determined by the natural gas quality, s 1 =9.707×10 − 2 s 4 -2.944×10 -7 , s 2 =-8.115×10 -3 s 4 +3.136×10 -8 , s 3 =-5.007s 4 +3.214×10 -5 , c=0.9986, thus z＝(9.707×10-2s4-2.944×10-7)t2+(-8.115×10-3s4+3.136×10-8)pt+(-5.007s4+3.214×10-5)t+s4p+0.9986；z=(9.707×10 -2 s 4 -2.944×10 -7 )t 2 +(-8.115×10 -3 s 4 +3.136×10 -8 )pt+(-5.007s 4 +3.214×10 -5 )t +s 4 p+0.9986; 系数s4通过如下步骤得到：The coefficient s4 is obtained by the following steps: 1)首先利用现有的天然气压缩因子公式计算出所述天然气气质在选定的标准气质温度t0，标准气质压力p0下的标准天然气压缩因子值z0；1) First, use the existing natural gas compression factor formula to calculate the standard natural gas compression factor value z 0 of the natural gas gas at the selected standard gas temperature t 0 and standard gas pressure p 0 ; 2)将计算得到的标准天然气压缩因子值z0以及对应的标准气质压力p0与标准气质温度t0代入系数s4的计算公式而得到：2) Substitute the calculated standard natural gas compression factor value z 0 and the corresponding standard gas pressure p 0 and standard gas temperature t 0 into the calculation formula of the coefficient s 4 to obtain:
CN201611085801.8A 2016-11-30 2016-11-30 Method of Calculating Compressibility Factors of Natural Gas Active CN106525180B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611085801.8A CN106525180B (en) 2016-11-30 2016-11-30 Method of Calculating Compressibility Factors of Natural Gas

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611085801.8A CN106525180B (en) 2016-11-30 2016-11-30 Method of Calculating Compressibility Factors of Natural Gas

Publications (2)

Publication Number Publication Date
CN106525180A CN106525180A (en) 2017-03-22
CN106525180B true CN106525180B (en) 2019-05-21

Family

ID=58354215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611085801.8A Active CN106525180B (en) 2016-11-30 2016-11-30 Method of Calculating Compressibility Factors of Natural Gas

Country Status (1)

Country Link
CN (1) CN106525180B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738579A (en) * 2018-07-02 2020-01-31 中国石油化工股份有限公司天然气分公司计量研究中心 Method for calculating precision of compression factor of natural gas flow computer and application thereof
CN114186165B (en) * 2020-09-14 2024-11-26 中国石油天然气股份有限公司 A method for obtaining natural gas compression factor
CN114372369A (en) * 2022-01-11 2022-04-19 西南石油大学 Method for calculating ultrahigh pressure natural gas compression factor

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
UA76071C2 (en) * 2005-07-06 2006-06-15 Affiliated Company Ukrtransgaz Method for determining flow rate of compressed air in operation of a compressor unit
CN103335682A (en) * 2013-07-15 2013-10-02 成都千嘉科技有限公司 Measuring method for gas flow of natural gas
CN104713606A (en) * 2015-03-12 2015-06-17 新奥科技发展有限公司 Method and device for measuring flow of multi-component gas
CN105842273A (en) * 2016-03-16 2016-08-10 华北电力大学（保定） Acquiring method and system for compression factor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7354291B2 (en) * 2006-03-01 2008-04-08 Panduit Corp. Plug locking assembly

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
UA76071C2 (en) * 2005-07-06 2006-06-15 Affiliated Company Ukrtransgaz Method for determining flow rate of compressed air in operation of a compressor unit
CN103335682A (en) * 2013-07-15 2013-10-02 成都千嘉科技有限公司 Measuring method for gas flow of natural gas
CN104713606A (en) * 2015-03-12 2015-06-17 新奥科技发展有限公司 Method and device for measuring flow of multi-component gas
CN105842273A (en) * 2016-03-16 2016-08-10 华北电力大学（保定） Acquiring method and system for compression factor

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
天然气压缩因子影响因素的研究;李师瑶等;《当代化工》;20161031;第45卷(第10期);第2378-2381页
天然气压缩因子的计算与体积计量;秦朝葵等;《天然气工业》;20031130(第06期);第130-134页
天然气压缩因子计算方法的评价;韩洪升;《油气储运》;19941231(第01期);第13-21页

Also Published As

Publication number Publication date
CN106525180A (en) 2017-03-22

Similar Documents

Publication Publication Date Title
CN112067304B (en) Method for measuring inlet flow of compressor in engine whole machine test
CN110514390B (en) A three-hole pressure probe measurement uncertainty evaluation method for two-dimensional flow field
EP3032230B1 (en) Flow meter and a method of calibration
CN105043511B (en) Vortex-shedding meter coefficient modification method based on Reynolds number
CN106525180B (en) Method of Calculating Compressibility Factors of Natural Gas
CN101603985A (en) Method for measuring sine signal with high accuracy
CN109358095B (en) Quantitative measurement method and system for gas concentration of each component of mixed gas
CN108152325A (en) A kind of method based on Guarded hot plate calibration heat-flow meter method conductometer
CN203376111U (en) Pressure gauge online calibration equipment for petroleum drilling operation
CN107340038A (en) Mass flowmenter automatic Calibration and measuring and testing device system and detection method
CN102768049A (en) Intelligent differential pressure type flow rate sensing device and design method of intelligent differential pressure type flow rate sensing device
CN107152958A (en) Mass flowmeter calibration detecting system and method based on multinomial algorithm
CN116086732A (en) High-precision leakage meter comparison system and method
CN104749544A (en) Electromagnetic compatibility current probe transmission resistance temperature error correction method
CN112098457B (en) Polynomial regression calibration method of thermal conductivity measuring instrument
RU2721900C2 (en) Method and instrument for determining gas properties by correlation
CN112964834A (en) Calibration method of dynamic calibrator for fixed pollution source
CN106643934B (en) Natural gas velocity of sound calculation method
CN102384816A (en) Suction resistance meter differential pressure sensor calibration system and method
CN116182969A (en) Quick correction method for error of ultrasonic gas meter
CN109781193B (en) In-situ calibration method for turbine flowmeter applied to test site
CN109060456B (en) A calibration device and method for a fixed source dilution channel sampler
RU2682540C9 (en) The method of setting the flow measurement channel with a restriction device
CN114459572B (en) A flow rate collection method for flow meter
CN117951897A (en) Conversion method of kinematic viscosity of phosphate fire-resistant oil at different temperatures

Legal Events

Date Code Title Description
C06 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: 20221206

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