9
где {Eij*, Gij*, Uij*, Kij*} - параметры бинарного взаимодействия, которые даны в таблице А.3.
Для расчета фактора сжимаемости по уравнению состояния (45) необходимо определить плотность ρм при заданных давлении (р, МПа) и температуре (Т, К).
Плотность ρм из УС (45) определяют по методу Ньютона в следующем итерационном процессе:
1) начальную плотность определяют по формуле
, (56)
где приведенное давление вычисляют из выражения
, (57)
2) плотность на k-м итерационном шаге определяют из выражений
. (58)
, (59)
где z(k-1) - рассчитывают из УС (45) при плотности на итерационном шаге (k-1), т.е. при ρм(k-1), а безразмерный комплекс А1 определяют из выражения
, (60)
при этом ρп = Кт3ρм(k-1);
4) критерий завершения итерационного процесса
, (61)
если критерий (61) не выполняется, то необходимо продолжить итерационный процесс, начиная с пункта 2) алгоритма.
После определения фактора сжимаемости при рабочих и стандартных условиях по формуле (1) рассчитывают коэффициент сжимаемости. Допускается рассчитывать фактор сжимаемости при стандартных условиях по формуле (20) ГОСТ 30319.1.
3.2.5 Уравнение состояния ВНИЦ СМВ
Во Всероссийском научно-исследовательском центре стандартизации, информации и сертификации сырья, материалов и веществ (ВНИЦ СМВ) для расчета фактора сжимаемости природного газа разработано уравнение состояния
, (62)
где ckl - коэффициенты УС;
ρп = ρм/ρпк - приведенная плотность;
Тп = Т/Тпк - приведенная температура;
ρм - молярная плотность, кмоль/м3;
ρпк и Тпк - псевдокритические параметры природного газа.
Коэффициенты УС определяют по формуле
, (63)