МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ УДК 622.279:681.5 ИНТЕРПОЛЯЦИОННЫЕ МОДЕЛИ ЗАВИСИМОСТИ ФАКТОРА СЖИМАЕМОСТИ ПРИРОДНОГО ГАЗА ОТ ДАВЛЕНИЯ И ТЕМПЕРАТУРЫ В.А. <...> В расчетах движения реального газа часто приходится применять уравнение его состояния, которое, как известно, записывается с использованием фактора сжимаемости z(р, Т) (в отечественной литературе чаще называемого коэффициентом сверхсжимаемости) в виде [1] ρ 10 M pR T z, 3 (1) где – плотность газа, кг/м3; р – абсолютное давление, МПа; R = 8,31451 – универсальная газовая постоянная, кДж/(кмоль·К); Т – термодинамическая температура, К; М – молярная масса газа, кг/кмоль, которая в зависимости от состава газовой смеси вычисляется по формуле M xM . i 1 ii n (2) В (2) xi – объемные доли компонентов в составе газа, а Mi – их молярные массы, которые можно найти, например, из табл. <...> Кроме расчета плотности реального газа фактор сжимаемости z(р, Т) приходится применять и при обработке данных газогидродинамических исследований скважин, и в оценках остаточных запасов газа по p/z-диаграммам и т. д. <...> Поскольку фактор сжимаемости зависит не только от состава газа, но и от условий, в которых газ находится, то невозможно получить какую-то усредненную величину или таблицу значений z, не допустив при этом большой погрешности. <...> В частности, в ОАО "Газпром" был разработан межгосударственный стандарт ГОСТ 0319-2-96 [3], устанавливающий обязательные рамки погрешностей при аналитическом расчете фактора сжимаемости и предлагающий четыре метода расчета, удовлетворяющие установленным требованиям. <...> Кроме того, известен [1, 4] аналитический метод расчета на основе уравнения Пенга-Робинсона. <...> Однако применение в вычислительных схемах громоздких расчетов по алгоритмам [1, 3] сильно замедляет работу программ для ЭВМ. <...> Желательно для конкретного состава газа иметь аппроксимационную формулу, позволяющую с достаточной точностью находить z в довольно широких диапазонах <...>