8
, (24)
в формулах (23) и (24) N - количество тестовых точек в массивах
, (25)
где Qрасч и Qтест - условное обозначение, соответственно, расчетного по аттестуемым методам и рассчитанного в перечислении 2) алгоритма тестового значений физического свойства природного газа (ρ, к, и, μ);
4) определяется погрешность расчета свойства Q по аттестуемым методам согласно ИСО 5168 [5]
, (26)
где δQ - погрешность расчета физических свойств по УС (1) и по уравнению для вязкости (15), значение которой для каждого свойства приведено в таблице 1.
Если для аттестуемых методов в качестве исходных данных используют плотность смеси природного газа при стандартных условиях (ρc), ее значение для тестовых смесей необходимо рассчитывать по УС (1). Допускается также рассчитывать плотность ρc по формуле (16) ГОСТ 30319.1 (см. п. 3.3.2).
ПРИЛОЖЕНИЕ А
(рекомендуемое)
Листинг программы расчета физических свойств природного газа
Расчет физических свойств природного газа по уравнению состояния (1) и по уравнению для вязкости (15) реализован на ПЭВМ, совместимых с IBM PC/AT/XT, на языке программирования ФОРТРАН-77.
C **************************************************************
C * *
С * Программа расчета физических свойств (плотности, показателя *
С * адиабаты, скорости звука и вязкости) природного газа по *
С * уравнению состояния ВНИЦ СМВ. *
С * *
C **************************************************************
IMPLICIT REAL*8(A-H,O-Z)
CHARACTER*26 AR
DIMENSION PI(100),TI(100),ROP(100,100),PAP(100,100),
*WP(100,100),ETAP(100,100)
COMMON/P/P/T/T/RON/RON/YI/YC(25)/NPR/NPR/Z/Z/TS/RO,PA,W
*/ETA/ETA/AR/AR(25)
200 WRITE(*,300)
300 FORMAT(18(/))
WRITE(*,400)
400 FORMAT(
*’ Расчет физических свойств природного газа’/
*’ по уравнению состояния’/////)
WRITE(*,1)
1 FORMAT(’ Введите исходные данные для расчета.’/)
WRITE(*,35)
35 FORMAT(’ Введите 0, если состав задан в молярных долях’/
*’ или 1, если состав задан в объемных долях ’\)