CN106525180A - Natural gas compression factor calculation method - Google Patents
Natural gas compression factor calculation method Download PDFInfo
- Publication number
- CN106525180A CN106525180A CN201611085801.8A CN201611085801A CN106525180A CN 106525180 A CN106525180 A CN 106525180A CN 201611085801 A CN201611085801 A CN 201611085801A CN 106525180 A CN106525180 A CN 106525180A
- Authority
- CN
- China
- Prior art keywords
- gas
- natural gas
- makings
- 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.)
- Granted
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 natural gas compression factor calculation method. The method is characterized by including the following steps that firstly, natural gas quality in a pipeline is determined; secondly, quality temperature t and quality pressure p of the natural gas quality are obtained through sampling in the pipeline, wherein t ranges from 0 DEG C to 60 DEG C and p ranges from 50 KPa to 4000 KPa; thirdly, a natural gas compression factor z is obtained according to the formula that z equals to s1t2 plus s2pt plus s3t plus s4p plus c, wherein s4 is a parameter determined according to the natural gas quality, c is a constant, s1 equals to A1s4 plus B1, s2 equals to A2s4 plus B2, s3 equals to A3s4 plus B3, and A1, B1, A2, B2, A3 and B3 are all constants.
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 technology
With the Large scale construction of developing rapidly for gas industry, particularly gas line network, for the day of Trade Measures
So throughput instrument is increasing, and its operating pressure is improved constantly, and range of flow also constantly increases.Natural gas flow measuring is many
The continuous measurement of parameter, multicomponent gases, with not reversibility, its accuracy of measurement receives many factors shadow for its value measurement
Ring.
Only under low pressure, high temperature, real gas just can approximately be counted as ideal gases, 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, particularly 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 process.
Real gas and this deviation of ideal gases, can illustrate using the ratio of PV and RT, this ratio is claimed
For compressibility factor, the definition of compressibility factor is：In authorized pressure with a temperature of, the volume of any mass gas is with the gas in phase
The ratio of the gas volume calculated by perfect gas law with the conditions of.Compressibility factor is represented with zed, if Z ＞ 1,
The volume of real gas is more than the volume of the ideal gases of equivalent at that the same temperature and pressure；If Z ＜ 1, identical
Temperature and pressure under real gas volume less than equivalent ideal gases volume.
Calculating compressibility factor has various methods, the AGA8 reports that American Gas association delivers, and is calculated by gas component
The compressibility factor of natural gas and other related hydrocarbons gases, AGA8-92DC methods are exactly the very important method of one of which, should
Computational methods are relatively complicated, and partial derivative and integral and calculating are repeatedly asked in design, simultaneously because the cost control of gas flowmeter, stream
The built-in monolithic microcomputer kernel computing capability of gauge is general, inefficiency when requiring in the face of high intensity, high-frequency calculating.
The content of the invention
The technical problem to be solved is the problem existed for above-mentioned prior art, there is provided a kind of simple, fast
Prompt and accurate Method of Calculating Compressibility Factors of Natural Gas.
The present invention solve the technical scheme that adopted of above-mentioned technical problem for：A kind of Method of Calculating Compressibility Factors of Natural Gas,
It is characterized in that：Comprise the steps：
1) determine the natural gas quality in pipeline；
2) sample in pipeline and obtain makings temperature t and makings pressure p of the natural gas quality, wherein the scope of t is 0
DEG C～60 DEG C, the scope 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 is determined, c is constant, s1=A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is
Constant.
In order that the calculation 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 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 s4Obtain as follows：
1) the natural gas quality temperature is calculated in the mark selected first with existing gas deviation factor formula
Quasi- makings temperature, standard makings pressure p0Under standard gas deviation factor value z0；
2) by calculated standard gas deviation factor value z0And corresponding standard makings pressure p0With standard makings
Temperature t0Substitute into coefficient s4Computing formula and obtain：
Compared with prior art, it is an advantage of the current invention that：Suitable for domestic gas flow rate calculation, can be accurate, efficient
Gas deviation factor instantly is calculated, is applied in the middle of the gas flowmeter of main flow instantly, work efficiency can be improved.
Description of the drawings
Fig. 1 is the flow chart of the computational methods of the present invention.
Specific embodiment
The present invention is described in further detail below in conjunction with accompanying drawing embodiment.
A kind of Method of Calculating Compressibility Factors of Natural Gas, it is contemplated that actual application, chooses the Main change model of makings temperature
Enclose for 0 DEG C to 60 DEG C, in the environment of the Main change scope of makings pressure is 50KPa to 4000KPa, there is provided a kind of easy day
Right gas compressibility factor computing formula, the loaded down with trivial details problem of the solution that the AGA8-92DC equations proposed in overcoming No. AGA8 report are present,
And on the basis of existing AGA8-92DC equations, the relative error of formula result of calculation of the present invention is controlled within 0.5%.
Specifically, gas deviation factor z is the quadratic polynomial with regard to makings temperature with makings pressure, is wherein comprised only
One coefficient s to be calibrated4, the form of formula is：Z=s1t2+s2pt+s3t+s4p+c.Wherein p is the makings pressure of natural gas to be measured
Power, t is 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 relevant amount of natural gas quality.Further illustrate, coefficient s1、s2、s3Can use and coefficient s4Relevant linear relation represents, s1
=A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is constant.
In order that the calculation 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 makings constituent, assisted using American Gas
The AGA8-92DC Equation for Calculating that No. 8 report AGA8 of meeting are proposed are in ethos pressure p0With ethos temperature t0Under natural gas
Compressibility factor size z0, by calculated z0And corresponding p0、t0Just s is obtained in substituting into following formula4The value of coefficient.
Therefore, the concrete form of computational methods of the present invention is：
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.
Table 1 various makings main component molar percentage (%)
Meanwhile, the visible Fig. 1 of flow process for the compressibility factor computing formula for obtaining being directed to the gas sample 1：
1) the AGA8-92DC Equation for Calculating proposed first with No. 8 reports of American Gas association goes out in makings temperature t0For
10 DEG C, makings pressure p0For standard gas deviation factor value z under 1000KPa0.Z is calculated (0When, 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) by calculated gas deviation factor value z0And corresponding makings pressure p0With makings temperature t0Substitute into system
Number s4Computing formula：
So as to be calculated coefficient s4=-2.478 × 10-5；
3) determine coefficient s4Value after, by coefficient s4Compressibility factor computing formula is corrected in substituting into 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
Thus, just can obtain being directed to the compressibility factor computing formula of gas sample 1, i.e.,：
Z=-2.699 × 10-6t2+2.324×10-7pt+1.562×10-4t-2.478×10-5p+0.9986
Now, the error of correlation can in order to verify the feasibility for obtaining the compressibility factor computing formula with regard to gas sample 1, be carried out
Analysis.In the case where a series of makings temperature and makings pressure condition are chosen, proposed with No. 8 report AGA8 of American Gas association
On the basis of the result of calculation of AGA8-92DC equations, the relative error size of formula of the present invention is calculated.Concrete condition can be shown in Table 2.
2 gas sample of table 1 is calculated and is compared
Embodiment two
In order to further illustrate the feasibility of the present invention, have chosen No. 8 reports of 5 kinds of American Gas associations again here and be previously mentioned
Gas sample, the concrete Component molar percentage ratio of this 5 kinds of gas samples can be shown in Table 1.
First the compressibility factor value of different gas samples is sampled, it is contemplated that the practical application of compressibility factor, therefore take makings
Temperature changes from 0 DEG C to 60 DEG C, and 5 DEG C is step-length；Makings pressure takes and changes from 50KPa to 2000KPa, takes wherein 50KPa extremely
300KPa with 50KPa as step-length, 300KPa to 2000KPa with 200KPa as step-length, using the AGA8- that proposes in AGA8 reports
92DC equations are carried out different makings temperature and are sampled with the compressibility factor value under makings pressure to 5 kinds of gas samples.Due to data sampling point
Comparison is more, and 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 by Fig. 1, it is determined that the compressibility factor computing formula corresponding to different gas samples (we select here
Ethos pressure is 1000KPa, and ethos temperature is 10 DEG C).Above-mentioned 5 kinds of gas samples are calculated at each according to the formula for obtaining
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 respectively in a part
Compressibility factor value on sampled point.
On the basis of the AGA8-92DC equations proposed by existing AGA8 reports, the compressibility factor value by the equation is obtained
Result of calculation z1, and calculate result of calculation z of gas deviation factor computing formula proposed by the present invention2.And missed according to relative
Difference computing formulaIt is calculated corresponding relative error size.
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 less than 0.5%, it was demonstrated that invention formula has stronger feasibility.
3 gas sample of table 2 is calculated and is compared
4 gas sample of table 3 is calculated and is compared
5 gas sample of table 4 is calculated and is compared
6 gas sample of table 5 is calculated and is compared
7 gas sample of table 6 is calculated and is compared
A system number s is illustrated above1、s2、s3, and coefficient s1、s2、s3There can also be other replacement schemes, meet phase
Requirement below 0.5% is controlled to error：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
And for 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
And for 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 (3)
1. a kind of Method of Calculating Compressibility Factors of Natural Gas, it is characterised in that：Comprise the steps：
1) determine the natural gas quality in pipeline；
2) sampling obtains makings temperature t and makings pressure p of the natural gas quality in the pipeline, wherein the scope of t be 0 DEG C～
60 DEG C, the scope 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 is determined, c is constant, s1=A1s4+B1, s2=A2s4+B2, s3=A3s4+B3, wherein A1、B1、A2、B2、A3、B3It is
Constant.
2. Method of Calculating Compressibility Factors of Natural Gas as claimed in claim 1, it is characterised in that：
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
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。
3. Method of Calculating Compressibility Factors of Natural Gas as claimed in claim 2, it is characterised in that：Coefficient s4As follows
Arrive：
1) natural gas quality is calculated in the standard makings temperature selected first with existing gas deviation factor formula
Degree t0, standard makings pressure p0Under standard gas deviation factor value z0；
2) by calculated standard gas deviation factor value z0And corresponding standard makings pressure p0With standard makings temperature
t0Substitute into coefficient s4Computing 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 true CN106525180A (en) | 2017-03-22 |
| CN106525180B 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) |
Cited By (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 (5)
| 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 |
| EP1994608A2 (en) * | 2006-03-01 | 2008-11-26 | Panduit Corporation | Plug locking assembly |
| 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 |
-
2016
- 2016-11-30 CN CN201611085801.8A patent/CN106525180B/en active Active
Patent Citations (5)
| 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 |
| EP1994608A2 (en) * | 2006-03-01 | 2008-11-26 | Panduit Corporation | Plug locking assembly |
| 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 |
|---|
| 李师瑶等: "天然气压缩因子影响因素的研究", 《当代化工》 * |
| 秦朝葵等: "天然气压缩因子的计算与体积计量", 《天然气工业》 * |
| 韩洪升: "天然气压缩因子计算方法的评价", 《油气储运》 * |
Cited By (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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106525180B (en) | 2019-05-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104713606B (en) | The flow-measuring method and device of multicomponent gas | |
| US5648605A (en) | Flowmeter calibration method | |
| JPWO2011078242A1 (en) | Mass flow controller system | |
| CN101846572B (en) | Method for decreasing basic error of pressure sensor | |
| CN105043511B (en) | Vortex-shedding meter coefficient modification method based on Reynolds number | |
| CN202522295U (en) | System for calibrating turbine flow meter under condition of large flow of low-temperature medium | |
| US20160169730A1 (en) | Flow meter and a method of calibration | |
| CN103969614B (en) | A kind of calibration steps of digital multimeter | |
| CN107655552A (en) | A kind of Ultrasonic water meter flow velocity modification method based on the non-measured section of propagation time difference and drift error | |
| CN103616313A (en) | On-line gas density measuring method | |
| CN210036908U (en) | High-pressure natural gas flow primary standard device with adjustable sonic nozzle | |
| CN109358095A (en) | A kind of method for quantitative measuring and system of mixed gas each component gas concentration | |
| CN109307573B (en) | Air leakage rate testing method for air preheater | |
| CN106525180A (en) | Natural gas compression factor calculation method | |
| CN205388506U (en) | Circulating high -pressure gas flow rate standard facility | |
| CN209673463U (en) | Fluorine gas, hex and nitrogen hybrid standard sample compounding system | |
| CN104964729B (en) | A kind of calibrating installation of fluid metering instrument | |
| CN116182969A (en) | Quick correction method for error of ultrasonic gas meter | |
| EP3287751A1 (en) | Gas meter system and heating-value derivation method | |
| RU2446005C1 (en) | Method for preparation of multi-component gas mixes | |
| CN104390831A (en) | Method and device for preparing standard methacrylic acid gas sample | |
| CN106643934A (en) | Calculation method for gas sound velocity | |
| CN113514621A (en) | Method for testing dynamic performance of gas sensor | |
| CN109781193A (en) | Turbine flowmeter calibrated in situ method applied to testing ground | |
| CN104280199A (en) | Mass spectrograph linear section testing method |
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 |