ГОСТ Р 8.769—2011
Приложение F
(справочное)
Подпрограмма SGERG.FOR на языке FORTRAN
Комментарий к программе.
Вириальное уравнение состояния GERG—88.
Подпрограмма SGERG
SGERG рассчитывает фактор сжимаемости природного газа с использованием неполного компонентного
состава газа.
Легальные копии этой программы могут быть получены учленов рабочей группы GERG по определению фак
тора сжимаемости природного газа в той версии, которая представлена в Технической монографии ТМ 5
(GERG,1991).
Расчеты основаны на следующих исходных данных:
от 0.0 до 0.3;
от 20 до 48;
от 0.55 до
0.9; от 0.0 до
0.1.
Допустимые интервалы:
- 1 — ХЗ; молярная доля С02
- 2 — HS: теплота сгорания МДж
- 3 — RM; относительная плотность*
- 4 — Х5; молярная доля Н2
Далее вводят входные переменные:
р — давление.бар;
Т — температура. °С
от 0 до 120;
от -23 до *65.
В программе рассчитывают Х2 — молярную долю азота и выводят ее на печать.
Результаты вычислений:
Z — фактор сжимаемости;
D — молярная плотность, моль • м~3.
СИНТАКСИС. CALL SGERG (Х2. ХЗ. HS. RM. Х5. Р. ТС. Z. D).
Значения коэффициентов, используемых в этой программе, соответствуют значениям, приведенным в под
программе GAS682
ot
20 июля 1988 г.,приложенной всообщении 8807 Лаборатории Ван-дер-Ваальса.Амстердам
J. Р. J. Michels. J. A. Schouten.
August 16.1991
Значения газовой постоянной, молярных масс, значения теплоты сгорания и плотности воздуха соответству
ют приведенным в ГОСТ 31369.
С
ссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссссс
с
SUBROUTINE SGERG(X2,X3.HS,RM.X5.P,TC.Z.D)
IMPLICIT DOUBLE PRECISION (A-H.O-Z)
IF( P.LT. 0.0 .OR. P GT.120.0)STOP ’PRESSURE OUT OF RANGE’
IF(TC.LT.-23.0 .OR. TC.GT .65.0) STOP 1TEMPERATURE OUT OF RANGE’
CALL SGERG 1(P.TC,X2,X3.X5.HS.RM.Z.D)
RETURN
END
C ...................................................................................................* ........................* ..............................................................................
SUBROUTINE SGERG1(P.TC,Q2.Q3.Q5.QM,RM.Z,D)
IMPLICIT DOUBLE PRECISION (A-H.O-Z)
COMMON /RBLOK/ AMOL.HS
COMMON /XBLOK/ Х1 .X2.X3.X11 .X12.X13.X22.X23.X33
>.X5.X7.X15.X17.X25.X55.X77
COMMON/MBLOK/GM1RO.GM1R1.GM2.GM3.GMS.GM7.FA.FB.RL.TO.H5.H7.R
HS = QM
X3 = Q3
X5 = Q5
* Измерения количества при Г * 0.0 “С; р = 1.01325 бар; температура сгорания: 25.0 °С.
25