ISO 20765-2:2015
(Main)Natural gas — Calculation of thermodynamic properties — Part 2: Single-phase properties (gas, liquid, and dense fluid) for extended ranges of application
Natural gas — Calculation of thermodynamic properties — Part 2: Single-phase properties (gas, liquid, and dense fluid) for extended ranges of application
ISO 20765-2:2015 specifies a method to calculate volumetric and caloric properties of natural gases, manufactured fuel gases, and similar mixtures, at conditions where the mixture may be in either the homogeneous (single-phase) gas state, the homogeneous liquid state, or the homogeneous supercritical (dense-fluid) state.
Gaz naturel — Calcul des propriétés thermodynamiques — Partie 2: Propriétés des phases uniques (gaz, liquide, fluide dense) pour une gamme étendue d'applications
L'ISO 20765-2:2015 spécifie une méthode de calcul des propriétés volumétriques et calorifiques des gaz naturels, des gaz naturels manufacturés, et des mélanges similaires, dans des conditions telles que le mélange peut exister à l'état gazeux homogène (phase unique), à l'état liquide homogène, ou à l'état supercritique homogène (fluide dense).
General Information
Buy Standard
Standards Content (Sample)
INTERNATIONAL ISO
STANDARD 20765-2
First edition
2015-01-15
Natural gas — Calculation of
thermodynamic properties —
Part 2:
Single-phase properties (gas, liquid,
and dense fluid) for extended ranges
of application
Gaz naturel — Calcul des propriétés thermodynamiques —
Partie 2: Propriétés des phases uniques (gaz, liquide, fluide dense)
pour une gamme étendue d’applications
Reference number
ISO 20765-2:2015(E)
©
ISO 2015
---------------------- Page: 1 ----------------------
ISO 20765-2:2015(E)
COPYRIGHT PROTECTED DOCUMENT
© ISO 2015
All rights reserved. Unless otherwise specified, no part of this publication may be reproduced or utilized otherwise in any form
or by any means, electronic or mechanical, including photocopying, or posting on the internet or an intranet, without prior
written permission. Permission can be requested from either ISO at the address below or ISO’s member body in the country of
the requester.
ISO copyright office
Case postale 56 • CH-1211 Geneva 20
Tel. + 41 22 749 01 11
Fax + 41 22 749 09 47
E-mail copyright@iso.org
Web www.iso.org
Published in Switzerland
ii © ISO 2015 – All rights reserved
---------------------- Page: 2 ----------------------
ISO 20765-2:2015(E)
Contents Page
Foreword .v
1 Scope . 1
2 Normative references . 2
3 Terms and definitions . 2
4 Thermodynamic basis of the method . 4
4.1 Principle . 4
4.2 The fundamental equation based on the Helmholtz free energy . 4
4.2.1 Background. 4
4.2.2 The Helmholtz free energy . 5
4.2.3 The reduced Helmholtz free energy . 5
4.2.4 The reduced Helmholtz free energy of the ideal gas . 6
4.2.5 The pure substance contribution to the residual part of the reduced
Helmholtz free energy . 6
4.2.6 The departure function contribution to the residual part of the reduced
Helmholtz free energy . 7
4.2.7 Reducing functions . 8
4.3 Thermodynamic properties derived from the Helmholtz free energy . 8
4.3.1 Background. 8
4.3.2 Relations for the calculation of thermodynamic properties in the
homogeneous region . 9
5 Method of calculation .11
5.1 Input variables .11
5.2 Conversion from pressure to reduced density .11
5.3 Implementation .12
6 Ranges of application .13
6.1 Pure gases .13
6.2 Binary mixtures .14
6.3 Natural gases .17
7 Uncertainty of the equation of state .18
7.1 Background .18
7.2 Uncertainty for pure gases . .18
7.2.1 Natural gas main components.18
7.2.2 Secondary alkanes .19
7.2.3 Other secondary components .21
7.3 Uncertainty for binary mixtures .21
7.4 Uncertainty for natural gases .23
7.4.1 Uncertainty in the normal and intermediate ranges of applicability of
natural gas .24
7.4.2 Uncertainty in the full range of applicability, and calculation of properties
beyond this range .25
7.5 Uncertainties in other properties .25
7.6 Impact of uncertainties of input variables .25
8 Reporting of results .25
Annex A (normative) Symbols and units .27
Annex B (normative) The reduced Helmholtz free energy of the ideal gas .29
Annex C (normative) Values of critical parameters and molar masses of the pure components .35
Annex D (normative) The residual part of the reduced Helmholtz free energy .36
Annex E (normative) The reducing functions for density and temperature .48
Annex F (informative) Assignment of trace components .55
© ISO 2015 – All rights reserved iii
---------------------- Page: 3 ----------------------
ISO 20765-2:2015(E)
Annex G (informative) Examples .57
Bibliography .60
iv © ISO 2015 – All rights reserved
---------------------- Page: 4 ----------------------
ISO 20765-2:2015(E)
Foreword
ISO (the International Organization for Standardization) is a worldwide federation of national standards
bodies (ISO member bodies). The work of preparing International Standards is normally carried out
through ISO technical committees. Each member body interested in a subject for which a technical
committee has been established has the right to be represented on that committee. International
organizations, governmental and non-governmental, in liaison with ISO, also take part in the work.
ISO collaborates closely with the International Electrotechnical Commission (IEC) on all matters of
electrotechnical standardization.
The procedures used to develop this document and those intended for its further maintenance are
described in the ISO/IEC Directives, Part 1. In particular the different approval criteria needed for the
different types of ISO documents should be noted. This document was drafted in accordance with the
editorial rules of the ISO/IEC Directives, Part 2 (see www.iso.org/directives).
Attention is drawn to the possibility that some of the elements of this document may be the subject of
patent rights. ISO shall not be held responsible for identifying any or all such patent rights. Details of
any patent rights identified during the development of the document will be in the Introduction and/or
on the ISO list of patent declarations received (see www.iso.org/patents).
Any trade name used in this document is information given for the convenience of users and does not
constitute an endorsement.
For an explanation on the meaning of ISO specific terms and expressions related to conformity
assessment, as well as information about ISO’s adherence to the WTO principles in the Technical Barriers
to Trade (TBT) see the following URL: Foreword - Supplementary information
The committee responsible for this document is ISO/TC 193, Natural Gas, Subcommittee SC 1, Analysis
of Natural Gas.
ISO 20765 consists of the following parts, under the general title Natural gas — Calculation of
thermodynamic properties:
— Part 1: Gas phase properties for transmission and distribution applications
— Part 2: Single-phase properties (gas, liquid, and dense fluid) for extended ranges of application
— Part 3: Two-phase properties (vapour-liquid equilibria)
© ISO 2015 – All rights reserved v
---------------------- Page: 5 ----------------------
INTERNATIONAL STANDARD ISO 20765-2:2015(E)
Natural gas — Calculation of thermodynamic properties —
Part 2:
Single-phase properties (gas, liquid, and dense fluid) for
extended ranges of application
1 Scope
This part of ISO 20765 specifies a method to calculate volumetric and caloric properties of natural gases,
manufactured fuel gases, and similar mixtures, at conditions where the mixture may be in either the
homogeneous (single-phase) gas state, the homogeneous liquid state, or the homogeneous supercritical
(dense-fluid) state.
NOTE 1 Although the primary application of this document is to natural gases, manufactured fuel gases,
and similar mixtures, the method presented is also applicable with high accuracy (i.e., to within experimental
uncertainty) to each of the (pure) natural gas components and to numerous binary and multi-component mixtures
related to or not related to natural gas.
For mixtures in the gas phase and for both volumetric properties (compression factor and density)
and caloric properties (for example, enthalpy, heat capacity, Joule-Thomson coefficient, and speed of
sound), the method is at least equal in accuracy to the method described in Part 1 of this International
Standard, over the full ranges of pressure p, temperature T, and composition to which Part 1 applies. In
some regions, the performance is significantly better; for example, in the temperature range 250 K to
275 K (–10 °F to 35 °F). The method described here maintains an uncertainty of ≤ 0,1 % for volumetric
properties, and generally within 0,1 % for speed of sound. It accurately describes volumetric and
caloric properties of homogeneous gas, liquid, and supercritical fluids as well as those in vapour-liquid
equilibrium. Therefore its structure is more complex than that in Part 1.
NOTE 2 All uncertainties in this document are expanded uncertainties given for a 95 % confidence level
(coverage factor k = 2).
The method described here is also applicable with no increase in uncertainty to wider ranges of
temperature, pressure, and composition for which the method of Part 1 is not applicable. For example, it
is applicable to natural gases with lower content of methane (down to 0,30 mole fraction), higher content
of nitrogen (up to 0,55 mole fraction), carbon dioxide (up to 0,30 mole fraction), ethane (up to 0,25 mole
fraction), and propane (up to 0,14 mole fraction), and to hydrogen-rich natural gases. A practical usage is
the calculation of properties of highly concentrated CO mixtures found in carbon dioxide sequestration
2
applications.
The mixture model presented here is valid by design over the entire fluid region. In the liquid and
dense-fluid regions the paucity of high quality test data does not in general allow definitive statements
of uncertainty for all sorts of multi-component natural gas mixtures. For saturated liquid densities of
LNG-type fluids in the temperature range from 100 K to 140 K (–280 °F to –208 °F), the uncertainty is
≤(0,1 – 0,3) %, which is in agreement with the estimated experimental uncertainty of available test data.
The model represents experimental data for compressed liquid densities of various binary mixtures
to within ±(0,1 – 0,2) % at pressures up to 40 MPa (5800 psia), which is also in agreement with the
estimated experimental uncertainty. Due to the high accuracy of the equations developed for the binary
subsystems, the mixture model can predict the thermodynamic properties for the liquid and dense-fluid
regions with the best accuracy presently possible for multi-component natural gas fluids.
© ISO 2015 – All rights reserved 1
---------------------- Page: 6 ----------------------
ISO 20765-2:2015(E)
2 Normative references
The following documents, in whole or in part, are normatively referenced in this document and are
indispensable for its application. For dated references, only the edition cited applies. For undated
references, the latest edition of the referenced document (including any amendments) applies.
ISO 7504, Gas Analysis — Vocabulary
ISO 14532, Natural gas — Vocabulary
ISO 20765-1, Natural gas — Calculation of thermodynamic properties — Part 1: Gas phase properties for
transmission and distribution applications
ISO 80000-5:2007, Quantities and units — Part 5: Thermodynamics
3 Terms and definitions
For the purposes of this document, the terms and definitions in ISO 80000-5:2007 and/or ISO 20765-1,
ISO 7504, ISO 14532, and the following apply.
NOTE 1 See Annex A for the list of symbols and units used in this part of ISO 20765.
NOTE 2 Figure 1 is a schematic representation of the phase behaviour of a typical natural gas as a function of
pressure and temperature. The positions of the bubble and dew lines depend upon the composition. This phase
diagram may be useful in understanding the definitions below.
14
SUPERCRITICAL
cricondenbar
13
DENSE FLUID
12
STATE
11
10
9
critical point dew
8
line
7
LIQUID PHASE
6
TWO-PHASE
5
cricondentherm
bubble
VAPOUR-
LIQUID
4
line
3
2
GAS
1
PHASE
0
100 150 200 250 300 350 400
Figure 1 — Phase diagram for a typical natural gas
3.1
bubble pressure
pressure at which an infinitesimal amount of vapour is in equilibrium with a bulk liquid for a
specified temperature
2 © ISO 2015 – All rights reserved
Pressure/MPa
---------------------- Page: 7 ----------------------
ISO 20765-2:2015(E)
3.2
bubble temperature
temperature at which an infinitesimal amount of vapour is in equilibrium with a bulk liquid for a
specified pressure
Note 1 to entry: The locus of bubble points is known as the bubble line.
Note 2 to entry: More than one bubble temperature may exist at a specific pressure. Moreover, more than one
bubble pressure may exist at a specified temperature, as explained in the example given in 3.6.
3.3
cricondenbar
maximum pressure at which two-phase separation can occur
3.4
cricondentherm
maximum temperature at which two-phase separation can occur
3.5
critical point
unique saturation point along the two-phase vapour-liquid equilibrium boundary where both the vapour
and liquid phases have the same composition and density
Note 1 to entry: The critical point is the point at which the dew line and the bubble line meet.
Note 2 to entry: The pressure at the critical point is known as the critical pressure and the temperature as the
critical temperature.
Note 3 to entry: A mixture of given composition may have one, more than one, or no critical points. In addition,
the phase behaviour may be quite different from that shown in Fig. 1 for mixtures (including natural gases)
containing, e.g., hydrogen or helium.
3.6
dew pressure
pressure at which an infinitesimal amount of liquid is in equilibrium with a bulk vapour for a
specified temperature
Note 1 to entry: More than one dew pressure may exist at the specified temperature. For example, isothermal
compression at 300 K with a gas similar to that shown in Figure 1: At low pressure the mixture is a gas. At just
above 2 MPa (the dew pressure), a liquid phase initially forms. As pressure increases more liquid forms in the
two-phase region, but a further increase in pressure reduces the amount of liquid (retrograde condensation) until
at about 8 MPa where the liquid phase disappears at the upper dew pressure, and the mixture is in the dense gas
phase. In the two-phase region, the overall composition is as specified, however the coexisting vapour and liquid
will have different compositions.
3.7
dew temperature
temperature at which an infinitesimal amount of liquid is in equilibrium with a bulk vapour for a
specified pressure
Note 1 to entry: More than one dew temperature may exist at a specified pressure, similar to the example given in 3.6.
Note 2 to entry: The locus of dew points is known as the dew line.
3.8
supercritical state
dense phase region above the critical point (often considered to be a state above the critical temperature
and pressure) within which no two-phase separation can occur
© ISO 2015 – All rights reserved 3
---------------------- Page: 8 ----------------------
ISO 20765-2:2015(E)
4 Thermodynamic basis of the method
4.1 Principle
The method is based on the concept that natural gas or any other type of mixture can be completely
characterized in the calculation of its thermodynamic properties by component analysis. Such an
analysis, together with the state variables of temperature and density, provides the necessary input
data for the calculation of properties. In practice, the state variables available as input data are generally
temperature and pressure, and it is thus necessary to first iteratively determine the density using the
equations provided here.
These equations express the Helmholtz free energy of the mixture as a function of density, temperature,
and composition, from which all other thermodynamic properties in the homogeneous (single-phase)
gas, liquid, and supercritical (dense-fluid) regions may be obtained in terms of the Helmholtz free energy
and its derivatives with respect to temperature and density. For example, pressure is proportional to
the first derivative of the Helmholtz energy with respect to density (at constant temperature).
NOTE These equations are also applicable in the calculation of two-phase properties (vapour-liquid
equilibria). Additional composition-dependent derivatives are required and are presented in Part 3 of this
International Standard.
The method uses a detailed molar composition analysis in which all components present in amounts
exceeding 0,000 05 mole fraction (50 ppm) are specified. For a typical natural gas, this might include
alkane hydrocarbons up to about C or C together with nitrogen, carbon dioxide, and helium. Typically,
7 8
isomers for alkanes C and higher may be lumped together by molar mass and treated collectively as the
6
normal isomer.
For some fluids, additional components such as C , C , water, and hydrogen sulfide may be present and
9 10
need to be taken into consideration. For manufactured gases, hydrogen, carbon monoxide, and oxygen
may also be present in the mixture.
More precisely, the method uses a 21-component analysis in which all of the major and most of the minor
components of natural gas are included (see Clause 6). Any trace component present but not identified as one
of the 21 specified components may be assigned appropriately to one of these 21 components (see Annex F).
4.2 The fundamental equation based on the Helmholtz free energy
4.2.1 Background
[1]
The GERG-2008 equation was published by the Lehrstuhl für Thermodynamik at the Ruhr-Universität
Bochum in Germany as a new wide-range equation of state for the volumetric and caloric properties of
[2] [1]
natural gases and other mixtures. It was originally published in 2007 and later updated in 2008.
[3]
The new equation improves upon the performance of the AGA-8 equation for gas phase properties and
in addition is applicable to the properties of the liquid phase, to the dense-fluid phase, to the vapour-
liquid phase boundary, and to properties for two-phase states. The ranges of temperature, pressure,
and composition to which the GERG-2008 equation of state applies are much wider than the AGA-8
equation and cover an extended range of application. The Groupe Européen de Recherches Gazières
(GERG) supported the development of this equation of state over several years.
The GERG-2008 equation is explicit in the Helmholtz free energy, a formulation that enables all
thermodynamic properties to be expressed analytically as functions of the free energy and of its
derivatives with respect to the state conditions of temperature and density. There is generally no need
for numerical differentiation or integration within any computer program that implements the method.
4 © ISO 2015 – All rights reserved
---------------------- Page: 9 ----------------------
ISO 20765-2:2015(E)
4.2.2 The Helmholtz free energy
The Helmholtz free energy a of a fluid mixture at a given mixture density ρ, temperature T, and molar
o r
composition x can be expressed as the sum of a describing the ideal gas behaviour and a describing
the residual or real-gas contribution, as follows:
or
aT(,ρρ,)xa=+(,Tx,) aT(,ρ ,)x (1)
4.2.3 The reduced Helmholtz free energy
The Helmholtz free energy is often used in its dimensionless form α=a/(RT) as
or
αδ(,τα,)xT=+(,ρα,)xx(,δτ,) (2)
In this equation, the reduced (dimensionless) mixture density δ is given by
ρ
δ = (3)
ρ ()x
r
and the inverse reduced (dimensionless) mixture temperature τ is given by
Tx()
r
τ = (4)
T
where
ρ and Τ are reducing functions for the mixture density and mixture temperature (see 4.2.7) depending
r r
on the molar composition of the mixture only.
r
The residual part α of the reduced Helmholtz free energy is given by
r rr
αδ(,τα,)xx=+(,δτ,) Δαδ(,τ,)x (5)
o
r
In this equation, the first term on the right-hand side α describes the contribution of the residual parts
o
of the reduced Helmholtz free energy of the pure substance equations of state, which are multiplied by
the mole fraction of the corresponding substance, and calculated at the reduced mixture variables δ and
r
τ (see equation (8)). The second term Δα is the departure function, which is the double summation over
all binary specific and generalized departure functions developed for the respective binary mixtures
(see equation (10)).
© ISO 2015 – All rights reserved 5
---------------------- Page: 10 ----------------------
ISO 20765-2:2015(E)
4.2.4 The reduced Helmholtz free energy of the ideal gas
o
The reduced Helmholtz free energy α represents the properties of the ideal-gas mixture at a given
mixture density ρ, temperature Τ, and molar composition x according to
N
o o
αρ(,Tx,)=+xT[(αρ,) lnx ] (6)
∑ iio i
i=1
o
In this equation, the term ∑x lnx is the contribution from the entropy of mixing, and αρ(,T) is the
i i
oi
dimensionless form of the Helmholtz free energy in the ideal-gas state of component i, as given by
∗
T T T
ρ R
o o o c,i o c,,i o o c,i
αρ(,T)l= nl++nn +n n + n ln sinh ϑ
oi oii,,12o oi,3 ∑ oik,,oik
ρ R T T T
c,i
k=46,
(7)
T
o o c,i
− n lnncosh ϑ
∑ oik, oik,
T
k=57,
where
ρ and Τ are the critical parameters of the pure components (see Annex C).
c,i c,i
o o
The values of the coefficients n and the parameters ϑ for all 21 components are given in Annex B.
oik, oik,
NOTE 1 The method prescribed is taken without change from the method prescribed in Part 1 of this
International Standard. The user should however be aware of significant differences that result inevitably from
the change in definition of the inverse reduced temperature τ between Part 1 and Part 2.
-1 -1 [4]
NOTE 2 R = 8,314 472 J·mol ·K was the internationally accepted standard for the molar gas constant at the
time of development of the equation of state. Equation (7) results from the integration of the equations for the
ideal-gas heat capacities taken from [5], where a different molar gas constant was used than the one adopted in
-1 -1
the mixture model presented here. The ratio R*/R with R*=8,314 51 J·mol ·K takes into account this difference
and therefore leads to the exact solution of the original equations for the ideal-gas heat capacity.
4.2.5 The pure substance contribution to the residual part of the reduced Helmholtz free energy
The contribution of the residual parts of the reduced Helmholtz free energy of the pure substance
r
equations of state α to the residual part of the reduced Helmholtz free energy of the mixture is
o
N
r r
αδ(,τα,)xx= (,δτ) (8)
o ∑ iio
i=1
6 © ISO 2015 – All rights reserved
---------------------- Page: 11 ----------------------
ISO 20765-2:2015(E)
where
r
αδ(,τ) is the residual part of the reduced Helmholtz free energy of component i (i.e., the residual part
oi
of the respective pure substance equation of state listed in Table 2) and is given by
K KK+
PolE,,iixp
Pol,i
cc
oik,
dt dt
r −δ
ooik,,ik ooik,,ik
αδ(,τδ)=+nnτδ τ e (9)
oii∑ oo,k ∑ ik,
k=1 kK=+1
Pol,i
r
The equations for α use the same basic structure as further detailed in Annex D.2. The values of the
oi
coefficients n and the exponents d , t and c for all 21 components are given in Annex D.2.2.
oi,k oi,k oi,k oi,k
4.2.6 The departure function contribution to the residual part of the reduced Helmholtz free
energy
The purpose of the departure function is to further improve the accuracy of the mixture model in the
description of thermodynamic properties in addition to fitting the parameters of the reducing functions
(see 4.2.7) when sufficiently accurate experimental data are available to characterize the properties of the
r
mixture. The departure function Δα of the multi-component mixture is the double summation over all
binary specific and generalized departure functions developed for the binary subsystems and is given by
N−1 N
rr
ΔΔαδ(,τα,)xx= (,δτ,) (10)
ij
∑ ∑
i=1 ji=+1
with
rr
Δαδ(,τα,)xx= xF (,δτ) (11)
ij ij ij ij
r r
In this equation, the function α (δ,τ) is the part of the departure function Δα (δ,τ,x ) that depends only
ij ij
on the reduced mixture variables δ and τ, as given by
K
Pol,ij
dt
r
ij,,kijk
αδ(,τδ)= n τ
ij ∑ ij,k
k=1
(12)
KK+
2
PolE,,ij xp ij
dt −ηηδ()−−εβ ()δγ−
ij,,kijk ij,,kijk ij,,kijk
+ neδτ
∑ ij,k
kK=+1
Pol,ij
where
r
αδ(,τ) was developed either for a specific binary mixture (a binary specifi
...
NORME ISO
INTERNATIONALE 20765-2
Première édition
2015-01-15
Gaz naturel — Calcul des propriétés
thermodynamiques —
Partie 2:
Propriétés des phases uniques (gaz,
liquide, fluide dense) pour une gamme
étendue d’applications
Natural gas — Calculation of thermodynamic properties —
Part 2: Single-phase properties (gas, liquid, and dense fluid) for
extended ranges of application
Numéro de référence
ISO 20765-2:2015(F)
©
ISO 2015
---------------------- Page: 1 ----------------------
ISO 20765-2:2015(F)
DOCUMENT PROTÉGÉ PAR COPYRIGHT
© ISO 2015
Droits de reproduction réservés. Sauf indication contraire, aucune partie de cette publication ne peut être reproduite ni utilisée
sous quelque forme que ce soit et par aucun procédé, électronique ou mécanique, y compris la photocopie, l’affichage sur
l’internet ou sur un Intranet, sans autorisation écrite préalable. Les demandes d’autorisation peuvent être adressées à l’ISO à
l’adresse ci-après ou au comité membre de l’ISO dans le pays du demandeur.
ISO copyright office
Case postale 56 • CH-1211 Geneva 20
Tel. + 41 22 749 01 11
Fax + 41 22 749 09 47
E-mail copyright@iso.org
Web www.iso.org
Publié en Suisse
ii © ISO 2015 – Tous droits réservés
---------------------- Page: 2 ----------------------
ISO 20765-2:2015(F)
Sommaire Page
Avant-propos .v
1 Domaine d’application . 1
2 Références normatives . 2
3 Termes et définitions . 2
4 Base thermodynamique de la méthode . 4
4.1 Principe . 4
4.2 Équation fondamentale basée sur l’énergie libre de Helmholtz . 4
4.2.1 Rappel . 4
4.2.2 Énergie libre de Helmholtz . 5
4.2.3 Énergie libre réduite de Helmholtz. 5
4.2.4 Énergie libre réduite de Helmholtz du gaz parfait . 6
4.2.5 Contribution de la substance pure à la partie résiduelle de l’énergie libre
réduite de Helmholtz . . 6
4.2.6 Contribution de la fonction de départ à la partie résiduelle de l’énergie
libre réduite de Helmholtz . 7
4.2.7 Fonctions réductrices . 8
4.3 Propriétés thermodynamiques établies à partir de l’énergie libre de Helmholtz . 8
4.3.1 Rappel . 8
4.3.2 Relations pour le calcul des propriétés thermodynamiques dans la
phase homogène . 9
5 Méthode de calcul .10
5.1 Variables d’entrée .10
5.2 Conversion de la pression en masse volumique réduite .11
5.3 Application .12
6 Plages d’application .12
6.1 Gaz purs .12
6.2 Mélanges binaires .14
6.3 Gaz naturels .17
7 Incertitude de l’équation d’état .18
7.1 Rappel .18
7.2 Incertitude inhérente aux gaz purs .18
7.2.1 Principaux constituants du gaz naturel .18
7.2.2 Alcanes secondaires .18
7.2.3 Autres constituants secondaires .20
7.3 Incertitude inhérente aux mélanges binaires .21
7.4 Incertitude inhérente aux gaz naturels .22
7.4.1 Incertitude associée aux plages normales et intermédiaires d’applicabilité
du gaz naturel .23
7.4.2 Incertitude dans le la plage d’application complète, et calcul des
propriétés au-delà de cette plage .25
7.5 Incertitudes associées à d’autres propriétés .25
7.6 Impact des incertitudes des variables d’entrée .25
8 Consignation des résultats dans un rapport .25
Annexe A (normative) Symboles et unités .27
Annexe B (normative) Énergie libre réduite de Helmholtz du gaz parfait .29
Annexe C (normative) Valeurs des paramètres critiques et des masses molaires des
constituants purs .35
Annexe D (normative) Partie résiduelle de l’énergie libre réduite de Helmholtz .36
Annexe E (normative) Fonctions réductrices pour la masse volumique et la température .48
© ISO 2015 – Tous droits réservés iii
---------------------- Page: 3 ----------------------
ISO 20765-2:2015(F)
Annexe F (informative) Attribution des constituants à l’état de traces .56
Annexe G (informative) Exemples .58
Bibliographie .61
iv © ISO 2015 – Tous droits réservés
---------------------- Page: 4 ----------------------
ISO 20765-2:2015(F)
Avant-propos
L’ISO (Organisation internationale de normalisation) est une fédération mondiale d’organismes
nationaux de normalisation (comités membres de l’ISO). L’élaboration des Normes internationales est
en général confiée aux comités techniques de l’ISO. Chaque comité membre intéressé par une étude
a le droit de faire partie du comité technique créé à cet effet. Les organisations internationales,
gouvernementales et non gouvernementales, en liaison avec l’ISO participent également aux travaux.
L’ISO collabore étroitement avec la Commission électrotechnique internationale (IEC) en ce qui concerne
la normalisation électrotechnique.
Les procédures utilisées pour élaborer le présent document et celles destinées à sa mise à jour sont
décrites dans les Directives ISO/CEI, Partie 1. Il convient, en particulier de prendre note des différents
critères d’approbation requis pour les différents types de documents ISO. Le présent document a été
rédigé conformément aux règles de rédaction données dans les Directives ISO/IEC, Partie 2 (voir www.
iso.org/directives).
L’attention est appelée sur le fait que certains des éléments du présent document peuvent faire l’objet de
droits de propriété intellectuelle ou de droits analogues. L’ISO ne saurait être tenue pour responsable
de ne pas avoir identifié de tels droits de propriété et averti de leur existence. Les détails concernant les
références aux droits de propriété intellectuelle ou autres droits analogues identifiés lors de l’élaboration
du document sont indiqués dans l’Introduction et/ou dans la liste des déclarations de brevets reçues par
l’ISO (voir www.iso.org/brevets).
Les appellations commerciales éventuellement mentionnées dans le présent document sont données pour
information, par souci de commodité, à l’intention des utilisateurs et ne sauraient constituer un engagement.
Pour une explication de la signification des termes et expressions spécifiques de l’ISO liés à l’évaluation de
la conformité, ou pour toute information au sujet de l’adhésion de l’ISO aux principes de l’OMC concernant
les obstacles techniques au commerce (OTC), voir le lien suivant: Avant-propos — Informations
supplémentaires.
Le comité chargé de l’élaboration du présent document est l’ISO 193, Gaz naturel, sous-comité SC 1,
Analyse du gaz naturel.
L’ISO 20765 comprend les parties suivantes, présentées sous le titre général Gaz naturel — Calcul des
propriétés thermodynamiques:
— Partie 1: Propriétés de la phase gazeuse pour des applications de transport et de distribution
— Partie 2: Propriétés des fluides monophasiques (gaz, liquide, fluide dense) pour un domaine d’application
étendu (ce document)
— Partie 3: Propriétés des milieux diphasiques (équilibres liquide-vapeur)
© ISO 2015 – Tous droits réservés v
---------------------- Page: 5 ----------------------
NORME INTERNATIONALE ISO 20765-2:2015(F)
Gaz naturel — Calcul des propriétés thermodynamiques —
Partie 2:
Propriétés des phases uniques (gaz, liquide, fluide dense)
pour une gamme étendue d’applications
1 Domaine d’application
La présente partie de l’ISO 20765 spécifie une méthode de calcul des propriétés volumétriques et
calorifiques des gaz naturels, des gaz naturels manufacturés, et des mélanges similaires, dans des
conditions telles que le mélange peut exister à l’état gazeux homogène (phase unique), à l’état liquide
homogène, ou à l’état supercritique homogène (fluide dense).
NOTE 1 Bien que ce document vise principalement les gaz naturels, les gaz naturels manufacturés, et les
mélanges similaires, la méthode présentée est également applicable avec une grande précision (c’est-à-dire, à
l’incertitude expérimentale près) à chacun des composants du gaz naturel (pur) et à de nombreux mélanges
binaires ou à plusieurs constituants liés au gaz naturel.
Pour les mélanges en phase gazeuse et pour les propriétés volumétriques (facteur de compressibilité et
masse volumique) et les propriétés calorifiques (par exemple, enthalpie, capacité thermique, coefficient
de Joule-Thomson, et vitesse du son), la méthode offre une précision au moins égale à la méthode décrite
dans la Partie 1 de la présente Norme internationale, pour toutes les plages de pression p, de température
T, et de composition pour lesquelles la Partie 1 est applicable. Dans certaines régions, la performance est
significativement meilleure; par exemple, dans la plage de température allant de 250 K à 275 K (–10 °F
à 35 °F). La méthode décrite ici conserve une incertitude ≤ 0,1 % pour les propriétés volumétriques,
et généralement de l’ordre de 0,1 % pour la vitesse du son. Elle décrit précisément les propriétés
volumétriques et calorifiques des fluides gazeux, liquides, et supercritiques homogènes ainsi que celles
de l’équilibre liquide-vapeur. Sa structure est par conséquent plus complexe que celle de la Partie 1.
NOTE 2 Toutes les incertitudes citées dans le présent document sont des incertitudes élargies données pour un
niveau de confiance de 95 % (facteur d’élargissement k = 2).
La méthode décrite ici est également applicable, sans augmentation de l’incertitude, à des plages plus
larges de température, de pression, et de composition pour lesquelles la méthode de la Partie 1 n’est pas
applicable. Par exemple, elle est applicable aux gaz naturels ayant une teneur plus faible en méthane
(fraction molaire inférieure à 0,30), une teneur supérieure en azote (fraction molaire inférieure ou égale
à 0,55), en dioxyde de carbone (fraction molaire inférieure ou égale à 0,30), en éthane (fraction molaire
inférieure ou égale à 0,25), et en propane (fraction molaire inférieure ou égale à 0,14), et aux gaz naturels
riches en hydrogène. Un usage pratique est le calcul des propriétés des mélanges à forte concentration
de CO rencontrés dans les applications de captage et stockage du dioxyde de carbone.
2
Le modèle de mélange présenté ici est applicable à toute la phase fluide. Dans les phases liquides et fluides
denses, le peu de données d’essai de bonne qualité ne permet pas en général d’estimer précisément
l’incertitude de toutes les sortes de mélanges de gaz naturels à plusieurs constituants. Pour la masse
volumique du liquide saturé de type GNL et dans la plage de température allant de 100 K à 140 K
(–280 °F à –208 °F), l’incertitude est ≤ (0,1 – 0,3) %, ce qui est en accord avec l’incertitude expérimentale
estimée des données d’essai disponibles. Le modèle représente les données expérimentales concernant
les masses volumiques des liquides comprimés de divers mélanges binaires à ± (0,1 à 0,2) % près à des
pressions inférieures ou égales à 40 MPa (5800 psia), ce qui est également en accord avec l’incertitude
expérimentale estimée. En raison de la grande précision des équations élaborées pour les sous-systèmes
binaires, le modèle de mélange peut prédire les propriétés thermodynamiques pour les phases liquides
et fluides denses avec la meilleure précision actuellement possible pour des fluides de type gaz naturel
à plusieurs constituants.
© ISO 2015 – Tous droits réservés 1
---------------------- Page: 6 ----------------------
ISO 20765-2:2015(F)
2 Références normatives
Les documents suivants, en totalité ou en partie, sont référencés de manière normative dans le présent
document et sont indispensables pour son application. Pour les références datées, seule l’édition citée
s’applique. Pour les références non datées, la dernière édition du document de référence s’applique (y
compris les éventuels amendements).
ISO 7504, Analyse des gaz — Vocabulaire
ISO 14532, Gaz naturel — Vocabulaire
ISO 20765-1, Gaz naturel — Calcul des propriétés thermodynamiques — Partie 1: Propriétés de la phase
gazeuse pour des applications de transport et de distribution
ISO 80000-5:2007, Grandeurs et unités — Partie 5: Thermodynamique
3 Termes et définitions
Pour les besoins de ce document, les termes et définitions de l’ISO 80000-5:2007 et/ou de l’ISO 20765-1,
l’ISO 7504, l’ISO 14532, et les suivants s’appliquent.
NOTE 1 Voir l’Annexe A pour la liste des symboles et unités utilisés dans la présente partie de l’ISO 20765.
NOTE 2 La Figure 1 est une représentation schématique du comportement des phases d’un gaz naturel typique
en fonction de la pression et de la température. Les positions des courbes de rosée et de bulle dépendent de la
composition. Ce diagramme de phase peut être utile à la compréhension des définitions suivantes.
Figure 1 — Diagramme de phase d’un gaz naturel type
3.1
pression de bulle
pression à laquelle une quantité infinitésimale de vapeur est en équilibre avec un liquide (phase continue)
à une température spécifiée
2 © ISO 2015 – Tous droits réservés
---------------------- Page: 7 ----------------------
ISO 20765-2:2015(F)
3.2
température de bulle
température à laquelle une quantité infinitésimale de vapeur est en équilibre avec un liquide (phase
continue) à une pression spécifiée
Note 1 à l’article: L’ensemble des points de bulle est dénommé « courbe de bulle ».
Note 2 à l’article: Plusieurs températures de bulle peuvent exister à une pression spécifique. De plus, plusieurs
pressions de bulle peuvent exister à une température spécifiée, comme expliqué dans l’exemple donné en 3.6.
3.3
cricondenbar
pression maximale à laquelle la séparation des deux phases peut se produire
3.4
cricondentherm
température maximale à laquelle la séparation des deux phases peut se produire
3.5
point critique
point de saturation unique sur la courbe d’équilibre liquide-vapeur où la phase vapeur et la phase liquide
ont des compositions et des masses volumiques identiques
Note 1 à l’article: Le point critique est le point auquel la courbe de rosée et la courbe de bulle se rencontrent.
Note 2 à l’article: La pression au point critique est appelée « pression critique » et la température « température
critique ».
Note 3 à l’article: Un mélange de composition donnée peut avoir un ou plusieurs points critiques ou n’en avoir
aucun. De plus, le comportement des phases peut être différent de celui représenté à la Figure 1 pour les mélanges
(y compris les gaz naturels) contenant par exemple de l’hydrogène ou de l’hélium.
3.6
pression de rosée
pression à laquelle une quantité infinitésimale de liquide est en équilibre avec une vapeur (phase
continue) à une température spécifiée
Note 1 à l’article: Plusieurs pressions de rosée peuvent exister à une température spécifique. Par exemple, la
compression isotherme à 300 K d’un gaz similaire à celui représenté à la Figure 1. A basse pression, le mélange
est gazeux. À une pression légèrement supérieure à 2 MPa (la pression de rosée), une phase liquide se forme.
A mesure que la pression augmente, davantage de liquide se forme dans la zone diphasique, mais une nouvelle
augmentation de pression réduit la quantité de liquide (condensation rétrograde) jusqu’à environ 8 MPa où la
phase liquide disparaît à une pression supérieure à la pression de rosée, et le mélange est en phase gazeuse
dense. Dans la zone diphasique, la composition globale est comme spécifiée, toutefois les phases vapeur et liquide
coexistantes auront des compositions différentes, la vapeur étant au point de rosée, et le liquide au point de bulle.
3.7
température de rosée
température à laquelle une quantité infinitésimale de liquide est en équilibre avec une vapeur (phase
continue) à une pression spécifiée
Note 1 à l’article: Plusieurs température de rosée peuvent exister à une pression spécifiée, d’une façon similaire à
l’exemple donné en 3.6.
Note 2 à l’article: L’ensemble des points de rosée est dénommé « courbe de rosée ».
3.8
état supercritique
zone de phase dense au-dessus du point critique (souvent considérée comme étant un état au-dessus de
la température critique et de la pression critique) dans laquelle une séparation en deux phases ne peut
pas se produire
© ISO 2015 – Tous droits réservés 3
---------------------- Page: 8 ----------------------
ISO 20765-2:2015(F)
4 Base thermodynamique de la méthode
4.1 Principe
La méthode est fondée sur le principe selon lequel le gaz naturel ou tout autre type de mélange peut
être totalement caractérisé par le calcul de ses propriétés thermodynamiques grâce à l’analyse des
constituants. Ce type d’analyse fournit, conjointement avec les variables d’état de température et de
masse volumique, les données d’entrée requises pour le calcul des propriétés. Dans la pratique, les
variables d’état disponibles sous la forme de données d’entrée sont généralement la température et la
pression; il est donc nécessaire de d’abord déterminer itérativement la masse volumique à l’aide des
équations fournies ici.
Ces équations expriment l’énergie libre de Helmholtz du mélange comme une fonction de la masse
volumique, de la température et de la composition, à partir de laquelle toutes les autres propriétés
thermodynamiques dans les phases gazeuse homogène (fluide monophasique), liquide homogène et
supercritique homogène (fluide dense) peuvent être obtenues en termes d’énergie libre de Helmholtz
et de ses dérivées par rapport à la température et à la masse volumique. Par exemple, la pression est
proportionnelle à la première dérivée de l’énergie de Helmholtz par rapport à la masse volumique (à
température constante).
NOTE Ces équations sont également applicables au calcul des propriétés de la zone diphasique (équilibres
liquide-vapeur). Des dérivées supplémentaires dépendantes de la composition sont requises et sont présentées
dans l’ISO 20765-3.
La méthode utilise une analyse détaillée de la composition molaire pour laquelle tous les constituants
présents en des quantités supérieures à une fraction molaire de 0,00005 (50 ppm) sont spécifiés. Pour un
gaz naturel type, cela peut inclure des hydrocarbures alcanes avec jusqu’à C ou C , ainsi que l’azote, le
7 8
dioxyde de carbone et l’hélium. En règle générale, les isomères pour les alcanes en C et au-delà peuvent
6
être regroupés par masse molaire et traités collectivement comme l’isomère normal.
Pour certains fluides, des constituants supplémentaires – les hydrocarbures en C et C , l’eau et le
9 10
sulfure d’hydrogène, par exemple – peuvent être présents et doivent être pris en compte. Pour les gaz
manufacturés, de l’hydrogène, du monoxyde de carbone, et de l’oxygène peuvent également être présents
dans le mélange.
La méthode utilise plus précisément une analyse de 21 constituants dans laquelle l’ensemble des
constituants majeurs et mineurs du gaz naturel sont inclus (voir Article 6). Tout constituant à l’état de
traces présent mais non identifié comme l’un des 21 constituants spécifiés peut être attribué de manière
appropriée à l’un de ces 21 constituants (voir Annexe F).
4.2 Équation fondamentale basée sur l’énergie libre de Helmholtz
4.2.1 Rappel
[1]
L’équation GERG-2008 a été publiée par le « Lehrstuhl für Thermodynamik at the Ruhr-Universität
Bochum » en Allemagne comme une nouvelle équation d’état à large plage pour les propriétés
volumétriques et calorifiques des gaz naturels et d’autres mélanges. Elle a été publiée initialement
[2] [1]
en 2007 et mise à jour ultérieurement en 2008. La nouvelle équation améliore la performance de
[3]
l’équation AGA-8 pour les propriétés de la phase gazeuse et par ailleurs est applicable aux propriétés
de la phase liquide, de la phase fluide dense, à la frontière des phases liquide-vapeur, et aux propriétés
des états diphasiques. Les plages de température, de pression, et de composition auxquelles l’équation
d’état GERG-2008 s’applique sont beaucoup plus larges que l’équation AGA-8 et couvrent une plage
d’application étendue. Le Groupe Européen de Recherches Gazières (GERG) a soutenu l’élaboration de
cette équation d’état pendant plusieurs années.
L’équation GERG-2008 établit explicitement l’énergie libre de Helmholtz, une formulation qui permet
d’exprimer analytiquement toutes les propriétés thermodynamiques en fonction de l’énergie libre de
Helmholtz et de ses dérivées par rapport aux conditions d’état de température et de masse volumique. Il
4 © ISO 2015 – Tous droits réservés
---------------------- Page: 9 ----------------------
ISO 20765-2:2015(F)
n’est généralement pas nécessaire de procéder à une différentiation ou à une intégration numérique au
moyen d’un quelconque programme informatique qui utilise cette méthode.
4.2.2 Énergie libre de Helmholtz
L’énergie libre de Helmholtz a d’un mélange fluide à une masse volumique ρ donnée, une température T
o
donnée, et une composition molaire x donnée peut être exprimée comme la somme de a décrivant le
r
comportement d’un gaz parfait et de a décrivant la contribution d’un gaz résiduel ou réel, comme suit:
or
aT(,ρρ,)xa=+(,Tx,) aT(,ρ ,)x (1)
4.2.3 Énergie libre réduite de Helmholtz
L’énergie libre de Helmholtz est souvent utilisée sous une forme adimensionnelle α = a/(RT) par
or
αδ(,τα,)xT=+(,ρα,)xx(,δτ,) (2)
Dans cette équation, la masse volumique réduite (adimensionnelle) du mélange δ est donnée par
ρ
δ= (3)
ρ ()x
r
et la température inverse réduite (adimensionnelle) τ est donnée par
Tx()
r
τ = (4)
T
où
ρ et Τ sont des fonctions réductrices de la masse volumique du mélange et la température du mélange
r r
(voir 4.2.7) qui dépendent uniquement de la composition molaire du mélange.
r
La partie résiduelle α de l’énergie libre réduite de Helmholtz est donnée par
r rr
αδ(,τα,)xx=+(,δτ,) Δαδ(,τ,)x (5)
o
r
Dans cette équation, le premier terme à partir de la droite α décrit la contribution des parties résiduelles
o
de l’énergie libre réduite de Helmholtz des équations d’état de la substance pure, qui sont multipliées par
la fraction molaire de la substance correspondante et calculées aux variables réduites δ et τ du mélange
r
(voir Équation (8)). Le second terme Δα est la fonction de départ, qui est la double somme de toutes
les fonctions binaires spécifiques et généralisées élaborées pour les mélanges binaires respectifs (voir
Équation (10)).
© ISO 2015 – Tous droits réservés 5
---------------------- Page: 10 ----------------------
ISO 20765-2:2015(F)
4.2.4 Énergie libre réduite de Helmholtz du gaz parfait
o
L’énergie libre réduite de Helmholtz α représente les propriétés du mélange de gaz parfait à une masse
volumique ρ donnée, une température Τ donnée, et une composition molaire x donnée du mélange
conformément à:
N
o o
αρ(,Tx,)=+xT[(αρ,) lnx ] (6)
∑ iio i
i=1
o
Dans cette équation, le terme ∑x lnx est la contribution de l’entropie de mélange, et αρ(,T) est la
i i
oi
forme adimensionnelle de l’énergie libre de Helmholtz à l’état de gaz parfait du constituant i:
∗
T T T
ρ R
o o o c,i o c,,i o o c,i
αρ(,T)l= nl++nn +n n + n ln sinh ϑ
oi oii,,12o oi,3 ∑ oik,,oik
ρ R T T T
c,i
k=46,
(7)
T
o o c,i
− n lnncosh ϑ
∑ oik, oik,
T
k=57,
où
ρ et Τ sont les paramètres critiques des constituants purs (voir Annexe C).
c,i c,i
o o
Les valeurs des coefficients n et des paramètres ϑ pour les 21 constituants sont données
oik, oik,
dans l’Annexe B.
La méthode spécifiée provient sans modification de la méthode spécifiée dans l’ISO 20765-1. Il convient
que l’utilisateur tienne compte des différences significatives qui résultent inévitablement de la
modification de la définition de la température réduite inverse, τ, entre l’ISO 20765-1 et l’ISO 20765-2.
−1 −1
NOTE 2 R = 8,314 472 J·mol ·K est la valeur standard internationalement reconnue en vigueur pour la
[4]
constante molaire universelle des gaz parfaits. L’Équation (7) résulte de l’intégration des équations des
[5]
capacités thermiques des gaz parfaits provenant de, où la constante molaire des gaz utilisée est différente
−1 −1
de celle adoptée dans le modèle de mélange présenté ici. L
...
Questions, Comments and Discussion
Ask us and Technical Secretary will try to provide an answer. You can facilitate discussion about the standard in here.