CN106525180B - Method of Calculating Compressibility Factors of Natural Gas - Google Patents
Method of Calculating Compressibility Factors of Natural Gas Download PDFInfo
- 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
- makings
- natural gas
- pressure
- temperature
- 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
- 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
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
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. a kind of Method of Calculating Compressibility Factors of Natural Gas, characterized by the following steps:
1) natural gas quality in pipeline is determined；
2) in pipeline sampling obtain the makings temperature t and makings pressure p of the natural gas quality, wherein the range of t be 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 the natural gas gas
The parameter that matter determines, 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, 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；
Coefficient s4It obtains as follows:
1) natural gas quality is calculated in selected standard makings temperature first with existing gas deviation factor formula
Spend t0, 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:
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 (1)
| 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 |
Citations (4)
| 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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7354291B2 (en) * | 2006-03-01 | 2008-04-08 | Panduit Corp. | Plug locking assembly |
-
2016
- 2016-11-30 CN CN201611085801.8A patent/CN106525180B/en active Active
Patent Citations (4)
| 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)
| 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 |
|---|---|---|
| CN104713606B (en) | The flow-measuring method and device of multicomponent gas | |
| EP3032230B1 (en) | Flow meter and a method of calibration | |
| CN112067304B (en) | Method for measuring inlet flow of compressor in engine whole machine test | |
| CN103969614B (en) | A kind of calibration steps of digital multimeter | |
| CN109358095B (en) | Quantitative measurement method and system for gas concentration of each component of mixed gas | |
| CN102768049A (en) | Intelligent differential pressure type flow rate sensing device and design method of intelligent differential pressure type flow rate sensing device | |
| CN109100051A (en) | The temperature correction and device of the dynamic response of temperature sensor | |
| CN106525180B (en) | Method of Calculating Compressibility Factors of Natural Gas | |
| CN112964834A (en) | Calibration method of dynamic calibrator for fixed pollution source | |
| CN209673463U (en) | Fluorine gas, hex and nitrogen hybrid standard sample compounding system | |
| CN106643934B (en) | Natural gas velocity of sound calculation method | |
| CN114184648B (en) | Moisture content calibration method for resistance-capacitance humidity sensor | |
| CN106918377A (en) | Calibrating installation, sensitivity determining module and correlation method for virtual flowmeter | |
| RU2721900C2 (en) | Method and instrument for determining gas properties by correlation | |
| CN112098457B (en) | Polynomial regression calibration method of thermal conductivity measuring instrument | |
| RU2682540C9 (en) | The method of setting the flow measurement channel with a restriction device | |
| CN109781193A (en) | Turbine flowmeter calibrated in situ method applied to testing ground | |
| CN103162771B (en) | The scaling method of apparatus for metering fluids | |
| CN117951897A (en) | Conversion method of kinematic viscosity of phosphate fire-resistant oil at different temperatures | |
| Teplukh et al. | Design of linear capillary measuring transducers for low gas flow rates | |
| CN116380339B (en) | Film gauge vacuum gauge calibration method and related equipment | |
| CN109990855A (en) | A kind of ultrasonic gas flowmeter mark condition velocity of sound calculation method | |
| TWI416619B (en) | Methods for performing actual flow verification | |
| CN116879513B (en) | Verification method, device, equipment and storage medium of gas analysis system | |
| Marshall et al. | Advances in Flow Measurement Using a Frictional Pressure Drop |
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 |
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 |
|
| TR01 | Transfer of patent right |