11
где ρci - плотность i-го компонента при стандартных условиях (см. таблицу Б.2);
gi - массовая доля i-го компонента;
N - количество основных компонентов.
Для расчета фактора сжимаемости по уравнению состояния (62) необходимо определить плотность ρм при заданных давлении (р, МПа) и температуре (Т, К).
Плотность ρм из УС (62) определяют по методу Ньютона в следующем итерационном процессе:
1) начальную плотность определяют по формуле
, (75)
где приведенное давление вычисляют из выражений
, (76)
, (77)
а псевдокритические плотность (ρпк), температуру (Тпк) и фактор Питцера (Ω) рассчитывают по формулам (64), (66) и (69);
2) плотность на k-м итерационном шаге определяется из выражений
, (78)
, (79)
где z(k-1) рассчитывают из УС (62) при плотности на итерационном шаге (k-1), т.е. при ρм(k-1), a безразмерный комплекс A1 определяют из выражения
, (80)
4) критерий завершения итерационного процесса.
, (81)
если критерий (81) не выполняется, то необходимо продолжить итерационный процесс, начиная с пункта 2) алгоритма.
После определения фактора сжимаемости при рабочих и стандартных условиях по формуле (1) рассчитывают коэффициент сжимаемости. Допускается рассчитывать фактор сжимаемости при стандартных условиях по формуле (20) ГОСТ 30319.1.
4 Влияние погрешности исходных данных на погрешность расчета коэффициента сжимаемости
При измерении расхода и количества природного газа, транспортируемого в газопроводах, давление (р), температуру (T), плотность при стандартных условиях (ρc) и состав (хi) измеряют с определенной погрешностью. Перечисленные параметры являются исходными данными для расчета коэффициента сжимаемости.
В соответствии с рекомендациями ИСО 5168 [16] погрешность расчета коэффициента сжимаемости, которая появляется в связи с погрешностью измерения исходных данных, определяют по формуле
, (82)