Что учитывает коэффициент сжимаемости газа

Коэффициент сжимаемости

Коэффициент сжимаемости не следует путать со сжимаемостью (также известной как коэффициент сжимаемости или изотермической сжимаемости ) материала, которая является мерой относительного изменения объема жидкости или твердого тела в ответ на изменение давления.

Содержание

Определение и физическое значение [ править ]

Коэффициент сжимаемости часто определяется в термодинамике и технике как:

В статистической механике описание таково:

Fugacity [ править ]

Коэффициент сжимаемости связан с летучестью соотношением

f = P exp ⁡ ( ∫ Z − 1 P d P ) <\displaystyle f=P\exp \left(\int <\frac

>dP\right)>

Обобщенные графики коэффициента сжимаемости для чистых газов [ править ]

Что касается сжимаемости газов, принцип соответствующих состояний указывает, что любой чистый газ при одинаковой пониженной температуре и пониженном давлении должен иметь одинаковый коэффициент сжимаемости. T r <\displaystyle T_> P r <\displaystyle P_>

Пониженные температура и давление определяются как

Вот и известны как критическая температура и критическое давление газа. Они являются характеристиками каждого конкретного газа с температурой, выше которой невозможно сжижить данный газ, и минимальным давлением, необходимым для сжижения данного газа при его критической температуре. Вместе они определяют критическую точку жидкости, выше которой не существует отдельных жидкой и газовой фаз данной жидкости. T c <\displaystyle T_> P c <\displaystyle P_> T c <\displaystyle T_> P c <\displaystyle P_>

Данные «давление-объем-температура» (PVT) для реальных газов варьируются от одного чистого газа к другому. Однако, когда коэффициенты сжимаемости различных однокомпонентных газов отображаются в зависимости от давления вместе с изотермами температуры, многие из графиков демонстрируют аналогичные формы изотерм.

Чтобы получить обобщенный график, который можно использовать для многих различных газов, значения пониженного давления и температуры, и используются для нормализации данных коэффициента сжимаемости. На рисунке 2 показан пример обобщенного графика коэффициента сжимаемости, полученного из сотен экспериментальных точек данных PVT для 10 чистых газов, а именно метана, этана, этилена, пропана, н-бутана, изопентана, н-гексана, азота, диоксида углерода и пар. P r <\displaystyle P_> T r <\displaystyle T_>

Как читать обобщенную диаграмму сжимаемости [ править ]

Чтобы прочитать диаграмму сжимаемости, необходимо знать пониженное давление и температуру. Если пониженное давление или температура неизвестны, необходимо найти уменьшенный удельный объем. В отличие от пониженного давления и температуры, уменьшенный удельный объем не определяется с помощью критического объема. Приведенный удельный объем определяется как

ν R = ν actual R T cr / P cr <\displaystyle \nu _=<\frac <\nu _<\text>>>/P_<\text>>>>

Наблюдения, сделанные с помощью обобщенной диаграммы сжимаемости [ править ]

Глядя на обобщенную диаграмму сжимаемости, можно сделать три наблюдения. Вот эти наблюдения:

Теоретические модели [ править ]

Вириальное уравнение особенно полезно для описания причин неидеальности на молекулярном уровне (очень немногие газы являются одноатомными), поскольку оно выводится непосредственно из статистической механики:

Z = 1 + B V m + C V m 2 + D V m 3 + … <\displaystyle Z=1+<\frac >>>+<\frac >^<2>>>+<\frac >^<3>>>+\dots >

Где коэффициенты в числителе известны как вириальные коэффициенты и являются функциями температуры.

Коэффициент сжимаемости связан с потенциалом межмолекулярной силы φ следующим образом:

Z = 1 + 2 π N A V m ∫ 0 ∞ ( 1 − exp ⁡ ( φ k T ) ) r 2 d r <\displaystyle Z=1+2\pi <\frac >>>>\int _<0>^<\infty >\left(1-\exp \left(<\frac <\varphi >>\right)\right)r^<2>dr>

В статье о реальном газе представлены более теоретические методы вычисления коэффициентов сжимаемости.

Физическая причина зависимости температуры и давления [ править ]

На рисунке справа показан обзор, охватывающий широкий диапазон температур. При низкой температуре (100 K) кривая имеет характерную форму галочки, восходящая часть кривой почти прямо пропорциональна давлению. При промежуточной температуре (160 К) наблюдается плавная кривая с широким минимумом; хотя участок высокого давления снова почти линейный, он больше не прямо пропорционален давлению. Наконец, при высокой температуре (400 К) Z больше единицы при всех давлениях. Для всех кривых Z приближается к единице для идеального газа при низком давлении и превышает это значение при очень высоком давлении.

Чтобы лучше понять эти кривые, на втором рисунке более подробно показано поведение при низких температурах и давлениях. Все кривые начинаются с Z, равного единице при нулевом давлении, а Z сначала уменьшается по мере увеличения давления. В этих условиях N ₂ представляет собой газ, поэтому расстояние между молекулами велико, но становится меньше по мере увеличения давления. Это увеличивает притягивающее взаимодействие между молекулами, сближая молекулы и делая объем меньше, чем у идеального газа при той же температуре и давлении. Более высокая температура снижает эффект притягивающих взаимодействий, и газ ведет себя почти идеальным образом.

На последних рисунках показано поведение при температурах, значительно превышающих критические. На отталкивающие взаимодействия по существу не влияет температура, но на притягивающее взаимодействие влияет все меньше и меньше. Таким образом, при достаточно высокой температуре отталкивающие взаимодействия преобладают при всех давлениях.

Это можно увидеть на графике, показывающем поведение при высоких температурах. По мере увеличения температуры начальный наклон становится менее отрицательным, давление, при котором Z является минимальным, становится меньше, а давление, при котором начинают преобладать отталкивающие взаимодействия, т.е. когда Z изменяется от менее единицы к значению больше единицы, становится меньше. При температуре Бойля (327 К для N ₂ ) эффекты притяжения и отталкивания нейтрализуют друг друга при низком давлении. Тогда Z остается на уровне идеального газа, равном единице, до давлений в несколько десятков бар. Выше температуры Бойля коэффициент сжимаемости всегда больше единицы и медленно, но неуклонно увеличивается с увеличением давления.

Экспериментальные значения [ править ]

Сжимаемость воздуха [ править ]

Источник

КОЭФФИЦИЕНТ СЖИМАЕМОСТИ ГАЗА

Смотреть что такое «КОЭФФИЦИЕНТ СЖИМАЕМОСТИ ГАЗА» в других словарях:

коэффициент сжимаемости газа — — [http://slovarionline.ru/anglo russkiy slovar neftegazovoy promyishlennosti/] Тематики нефтегазовая промышленность EN gas compressibilitygas deviation factorgas compressibility factor … Справочник технического переводчика

Коэффициент сжимаемости — отношение относительной вертикальной деформации (изменения коэффициента пористости) к давлению, вызвавшему эту деформацию. Источник: ГОСТ 30416 96: Грунты. Лабораторные испытания. Общие положения оригинал документа … Словарь-справочник терминов нормативно-технической документации

коэффициент изотермической сжимаемости газа — — [http://slovarionline.ru/anglo russkiy slovar neftegazovoy promyishlennosti/] Тематики нефтегазовая промышленность EN coefficient of isothermal compresibility of gas … Справочник технического переводчика

коэффициент сжимаемости — Коэффициент Z = pV/RT, характеризующий отклонение свойств данного вещества от свойств идеального газа … Политехнический терминологический толковый словарь

Коэффициент вязкости — Механика сплошных сред Сплошная среда Классическая механика Закон сохранения массы · Закон сохранения импульса … Википедия

ГОСТ Р 8.740-2011: Государственная система обеспечения единства измерений. Расход и количество газа. Методика измерений с помощью турбинных, ротационных и вихревых расходомеров и счетчиков — Терминология ГОСТ Р 8.740 2011: Государственная система обеспечения единства измерений. Расход и количество газа. Методика измерений с помощью турбинных, ротационных и вихревых расходомеров и счетчиков оригинал документа: 3.2.1 вспомогательные… … Словарь-справочник терминов нормативно-технической документации

МИ 3082-2007: Государственная система обеспечения единства измерений. Выбор методов и средств измерений расхода и количества потребляемого природного газа в зависимости от условий эксплуатации на узлах учета. Рекомендации по выбору рабочих эталонов для их поверки — Терминология МИ 3082 2007: Государственная система обеспечения единства измерений. Выбор методов и средств измерений расхода и количества потребляемого природного газа в зависимости от условий эксплуатации на узлах учета. Рекомендации по выбору… … Словарь-справочник терминов нормативно-технической документации

теплофизические характеристики газа — 3.2 теплофизические характеристики газа: Величины, характеризующие теплофизические свойствагаза. Примечание В настоящем документе в качестве теплофизических характеристик газа приняты плотность, коэффициент сжимаемости, вязкость, показатель… … Словарь-справочник терминов нормативно-технической документации

Сжимаемость электронного газа — термодинамический коэффициент характеризующий изменение давления электронного газа при изменении его объёма. Для электронного газа, по аналогии с обычным идеальным газом можно вводят понятие сжимаемости K, обратная величина которой определяется… … Википедия

форма — 3.2 форма (form): Документ, в который вносятся данные, необходимые для системы менеджмента качества. Примечание После заполнения форма становится записью. Источник: ГОСТ Р ИСО/ТО 10013 2007: Менеджмент организации. Руководство по документированию … Словарь-справочник терминов нормативно-технической документации

Источник

1 Назначение и область применения

Настоящий стандарт устанавливает четыре метода определения коэффициента сжимаемости природного газа: при неизвестном полном компонентном составе природного газа (два метода) и известном компонентном составе.

Стандарт устанавливает предпочтительные области применения каждого метода по измеряемым параметрам (давление, температура, плотность природного газа при стандартных условиях и компонентный состав природного газа), однако не запрещает использование любого из методов и в других областях.

Допускается применять любые другие методы расчета коэффициента сжимаемости, однако погрешность расчета коэффициента сжимаемости по этим методам не должна превышать погрешностей, приведенных в настоящем стандарте (см. 3.2.1).

Используемые в настоящем стандарте определения и обозначения приведены в соответствующих разделах ГОСТ 30319.0.

2 Нормативные ссылки

В настоящем стандарте использованы ссылки на следующие стандарты:

ГОСТ 30319.0-96 Газ природный. Методы расчета физических свойств. Общие положения

ГОСТ 30319.1-96 Газ природный. Методы расчета физических свойств. Определение физических свойств природного газа, его компонентов и продуктов его переработки

3 Определение коэффициента сжимаемости

3.1 Общие положения

Коэффициент сжимаемости вычисляют по формуле

Рабочие условия характеризуются такими давлениями и температурами, которые определяются измерениями в процессе добычи, переработки и транспортирования природного газа. Давление p_c и температура T_c при стандартных условиях приведены в ГОСТ 30319.0.

3.2 Методы расчета коэффициента сжимаемости

3.2.1 Пределы применимости методов расчета и область их применения и погрешности расчета коэффициента сжимаемости

В таблице 1 приведены общие результаты апробации методов расчета. Апробация проведена на обширном массиве высокоточных экспериментальных данных о факторе сжимаемости природного газа 1.

Погрешность данных не превышает 0,1 %.

Область применения и погрешность метода расчета

Отклонения от экспериментальных данных

Для расчета коэффициента сжимаемости природного газа при определении его расхода и количества рекомендуется применять:

Метод NX19 мод. и уравнение состояния GERG-91 мод. могут быть использованы при неизвестном полном компонентном составе природного газа, расчет по этим методам не требует применения ЭВМ.

Расчет по уравнениям состояния AGA8-92DC и ВНИЦ СМВ может быть осуществлен только при наличии ЭВМ и известном полном компонентном составе природного газа, при этом должны быть выдержаны следующие диапазоны концентраций компонентов (в мол. %):

пропан £ 3,5 бутаны £ 1,5

азот £ 15 диоксид углерода £ 15

сероводород £ 30 (УС ВНИЦ СМВ) и £ 0,02 (УС AGA8-92DC)

Выбор конкретного метода расчета коэффициента сжимаемости допускается определять в контракте между потребителем природного газа и его поставщиком с учетом требований настоящего стандарта.

В таблице 1 приняты следующие обозначения:

, (2)

, (3)

, (4)

, (5)

Погрешность расчета коэффициента сжимаемости d приведена в таблице 1 без учета погрешности исходных данных.

3.2.2 Модифицированный метод NX19 мод.

В соответствии с требованиями стандарта Германии [17] расчет фактора сжимаемости по модифицированному методу NX19 мод. основан на использовании уравнения следующего вида

где , (7)

, (8)

, (9)

, (10)

, (11)

Корректирующий множитель F в зависимости от интервалов параметров р_а и D Т_а вычисляют по формулам:

, (12)

, (13)

Параметры р_a и Т_a определяются по следующим соотношениям:

, (15)

, (16)

В формулах (17), (18) вместо молярных долей диоксида углерода и азота допускается применять их объемные доли (r_y и r_a).

Европейская группа газовых исследований на базе экспериментальных данных, собранных в [12], и уравнения состояния вириального типа [18], разработала и опубликовала в [13, 14] УС

Коэффициенты уравнения состояния определяют из следующих выражений:

, (20)

, (21)

, (24)

, (25)

, (26)

, (28)

, (29)

, (30)

, (31)

, (32)

. (33)

В формулах (23), (27) Н рассчитывают по выражению

, (34)

В выражении (35) молярную долю эквивалентного углеводорода (x_э) рассчитывают с использованием формулы (22), а фактор сжимаемости при стандартных условиях (z_с) рассчитывают по формуле (24) ГОСТ 30319.1, а именно

После определения коэффициентов уравнения состояния (19) В_m и С_m рассчитывают фактор сжимаемости при заданных давлении (р, МПа) и температуре (Т, К) по формуле

, (38)

, (39)

, (40)

, (41)

Коэффициент сжимаемости природного газа рассчитывают по формуле (1), а именно

, (44)

Фактор сжимаемости при стандартных условиях z_с рассчитывают по формуле (36).

В проекте стандарта ISO/TC 193 SC1 № 62 [15] Американской Газовой Ассоциацией для расчета фактора сжимаемости предложено использовать уравнение состояния

Если состав газа задан в объемных долях, то молярные доли рассчитываются по формуле (12) ГОСТ 30319.1.

Приведенную плотность определяют по формуле

, (46)

Коэффициенты УС рассчитывают из следующих соотношений:

где N — количество компонентов в природном газе.

Бинарные параметры <E_ij, G_ij> и параметры <U, G, K_m, Q, F> рассчитывают с использованием следующих уравнений:

, (49)

, (50)

, (51)

, (52)

, (54)

, (55)

Параметры бинарного взаимодействия, которые не приведены в этой таблице, а также при i = j, равны единице.

Для расчета фактора сжимаемости по уравнению состояния (45) необходимо определить плотность r _м при заданных давлении (р, МПа) и температуре (Т, К).

Плотность r _м из УС (45) определяют по методу Ньютона в следующем итерационном процессе:

1) начальную плотность определяют по формуле

, (56)

где приведенное давление вычисляют из выражения

, (57)

2) плотность на k-м итерационном шаге определяют из выражений

. (58)

, (59)

, (60)

4) критерий завершения итерационного процесса

если критерий (61) не выполняется, то необходимо продолжить итерационный процесс, начиная с пункта 2) алгоритма.

После определения фактора сжимаемости при рабочих и стандартных условиях по формуле (1) рассчитывают коэффициент сжимаемости.

3.2.5 Уравнение состояния ВНИЦ СМВ

Во Всероссийском научно-исследовательском центре стандартизации, информации и сертификации сырья, материалов и веществ (ВНИЦ СМВ) для расчета фактора сжимаемости природного газа разработано уравнение состояния

Коэффициенты УС определяют по формуле

, (63)

Псевдокритические параметры природного газа и его фактор Питцера вычисляют по формулам:

где , (65)

(; )

где , (67)

; (68)

(; )

где , (70)

Если заданный компонентный состав природного газа включает, кроме основных, другие компоненты (но не более 1 % в сумме), то молярные доли этих компонентов прибавляют к соответствующим долям основных компонентов следующим образом:

— ацетилен и этилен к этану;

— пропилен к пропану;

— углеводороды от н-пентана и выше к н-бутану;

— прочие компоненты к азоту.

Если состав газа задан в объемных долях, то молярные доли рассчитывают по формуле (12) ГОСТ 30319.1.

Для расчета фактора сжимаемости по уравнению состояния (62) необходимо определить плотность r _м при заданных давлении (р, МПа) и температуре (Т, К).

Плотность r _м из УС (62) определяют по методу Ньютона в следующем итерационном процессе:

1) начальную плотность определяют по формуле

, (75)

где приведенное давление вычисляют из выражений

, (76)

, (77)

а псевдокритические плотность ( r _пк), температуру (Т_пк) и фактор Питцера ( W ) рассчитывают по формулам (64), (66) и (69);

2) плотность на k-м итерационном шаге определяется из выражений

, (78)

, (79)

, (80)

4) критерий завершения итерационного процесса.

если критерий (81) не выполняется, то необходимо продолжить итерационный процесс, начиная с пункта 2) алгоритма.

После определения фактора сжимаемости при рабочих и стандартных условиях по формуле (1) рассчитывают коэффициент сжимаемости.

4 Влияние погрешности исходных данных на погрешность расчета коэффициента сжимаемости

При измерении расхода и количества природного газа, транспортируемого в газопроводах, давление (р), температуру (T), плотность при стандартных условиях ( r _c) и состав (х_i) измеряют с определенной погрешностью. Перечисленные параметры являются исходными данными для расчета коэффициента сжимаемости.

В соответствии с рекомендациями ИСО 5168 [16] погрешность расчета коэффициента сжимаемости, которая появляется в связи с погрешностью измерения исходных данных, определяют по формуле

, (84)

q_k — условное обозначение k-го параметра исходных данных (р. Т, r _c, х_i,);

Коэффициент сжимаемости ` К (среднее значение) рассчитывают по выбранному рекомендуемому методу расчета при средних параметрах q_k.

Общую погрешность расчета коэффициента сжимаемости определяют по формуле

, (85)

Для методов NX19 мод. и УС GERG-91 мод. допускается рассчитывать погрешность d _и.д по формуле

— при расчете К по методу NX19 мод.

, (87)

, (88)

, (89)

, (90)

, (91)

— при расчете К по методу GERG-91

, (92)

, (93)

, (94)

, (95)

. (96)

5 Программная и техническая реализация расчета коэффициента сжимаемости

Расчет коэффициента сжимаемости природного газа по указанным в стандарте методам реализован на ПЭВМ, совместимых с IBM PC/AT/XT, на языке программирования ФОРТРАН-77. Листинг программы приведен в приложении В.

В приложениях Г и Д приведены примеры расчета соответственно коэффициента сжимаемости и погрешности вычисления коэффициента сжимаемости, которая вызвана погрешностью определения исходных данных.

ПРИЛОЖЕНИЕ А

Таблицы констант и параметров уравнения состояния AGA8-92DC

Параметры бинарного взаимодействия

ПРИЛОЖЕНИЕ Б

Таблицы коэффициентов и параметров уравнения состояния ВНИЦ СМВ

1 Плотность ( r _ki ), температура ( T _ki ) в критической точке и фактор Питцера ( W _i ) отличаются от литературных данных и применимы только для уравнения состояния ВНИЦ СМВ.

ПРИЛОЖЕНИЕ В

(рекомендуемое)

Листинг программы расчета коэффициента сжимаемости природного газа

С * Программа расчета коэффициента сжимаемости природного газа *

С * (основной модуль) *

DATA AR/’ метана (СН4)’,’ этана (С2Н6)’,’ пропана (С3Н8)’,

*’ диоксида углерода (СO2)’,’ сероводорода (H2S)’,

*’ ацетилена (С2Н2)’,’ этилена (С2Н4)’,’ пропилена (С3Н6)’,

*’ нео-пентана (нео-С5Н12)’,’ н-гексана (н-С6Н14)’,

*’ бензола (С6Н6)’,’ н-гептана (н-С7Н16)’,’ толуола (С7Н8)’,

*’ н-декана (н-С10Н22)’,’ гелия (Не)’,’ водорода (Н2)’,

*’ моноксида углерода (СО)’,’ кислорода (О2)’/

IF(NVAR.EQ.5) GO TO 134

1 FORMAT(’ Введите исходные данные для расчета.’/)

*’ Плотность при 293.15 К и 101.325 кПа, в кг/куб.м ’

53 FORMAT(’ Введите 0, если состав азота и диоксида углерода’,

*’ задан в молярных долях’/

*’ или 1, если состав этих компонентов задан’,

3 FORMAT (’ Значение молярной доли, в мол. %’)

33 FORMAT(’ Значение объемной доли, в об. %’)

WRITE(*,’(A\)’) ’ диоксида углерода (С02) ’

35 FORMAT(’ Введите 0, если состав задан в молярных долях’/

*’ или 1, если состав задан в объемных долях ’\)

*’ Введите количество точек по давлению: ’

*’ Введите количество точек по температуре: ’

*’ Введите значения давлений в МПа: ’

*’ Введите значения температур в К: ’

*’ Ввод исходных данных завершен. ’

GO TO (10,20,30,40) NVAR

20 CALL GERG2(ICALC,YA,YY)

30 CALL AGA8DC(ICALC)

IF(NVAR.EQ.l) CALL NX19(YA,YY)

IF(NVAR.EQ.2) CALL GERG2(ICALC,YA,YY)

IF(NVAR.EQ.3) CALL AGA8DC(ICALC)

IF(NVAR.EQ.4) CALL VNIC(ICALC)

IF (NRANGE) 134,134,200

*10X,’ Расчет коэффициента сжимаемости природного газа’//

*10Х,’ 1. Модифицированный метод NX 19 ’/

*10Х,’ 2. Уравнение состояния GERG-91 ’/

*10Х,’ 3. Уравнение состояния AGA8-92DC ’/

*10Х,’ 4. Уравнение состояния ВНИЦ СМВ ’/

*’Введите порядковый номер метода расчета или 5 для выхода в ДОС’,

*’ Выбранная Вами методика при заданных параметрах «не работает»’/

CHARACTER*26 AR(25), FNAME

*’(модифицированный метод NX19)’,

*’(уравнение состояния GERG-91)’,

*’(уравнение состояния AGA8-92DC)’,

*’(уравнение состояния ВНИЦ СМВ)’/

IF(NYST.EQ.2) WRITE(*,’(A\)’) ’ Введите имя файла ’

*’ Включите принтер, вставьте бумагу и нажмите ’

*13X,’Коэффициент сжимаемости природного газа.’/

1 FORMAT(’ Плотность при 293.15 К и 101.325 кПа ’,F6.4,’ кг/куб.м’)

3 FORMAT(’ Содержание в мол. %’)

33 FORMAT(’ Содержание в об.%’)

IF(YC(I).NE.0D0.AND.I.NE.25) GO TO 9

PAUSE ’ Для продолжения вывода нажмите ’

PAUSE ’ Для продолжения вывода нажмите ’

*’ Для продолжения вывода вставьте бумагу и нажмите ’

21 FORMAT(’ p, МПа ’,6А9)

23 FORMAT(10X,6(:,’ | ’,F6.2))

IF(NLN.EQ.NL) GO TO 133

IF(J.EQ.NP.AND.I+NL-1.EQ.NT) GO TO 29

PAUSE ’ Для продолжения вывода нажмите ’

IF(J.EQ.NP.AND.I+NL-1.EQ.NT) GO TO 29

*’ Для продолжения вывода вставьте бумагу и нажмите ’

PAUSE ’ Вывод завершен, для продолжения работы нажмите ’

IF(NBOLB.EQ.l) GO TO 22

С * Подпрограмма расчета коэффициента сжимаемости природного *

С * газа по модифицированному методу NX19. *

IF(NCONT.EQ.l) GO TO 134

С * Подпрограмма расчета коэффициента сжимаемости природного *

С * газа по модифицированному уравнению состояния GERG-91. *

IF(ICALC.EQ.2) GO TO 3

IF(Z.EQ.0D0) GO TO 133

С Расчет теплоты сгорания эквивалентного углеводорода (Н)

IF(Z.EQ.0D0) GO TO 133

IF(Z.EQ.0D0) GO TO 133

IF(Z.EQ.0D0) GO TO 133

С Подпрограмма, реализующая схему Кардано для определения

С фактора сжимаемости из уравнения состояния

С * Подпрограмма расчета коэффициента сжимаемости природного *

С * газа по уравнению состояния AGA8-92DC. *

IF(ICALC.NE.l) GO TO 3

IF(Z.EQ.0D0) GO TO 133

IF(NPR.EQ.0D0) GO TO 5

IF(YI(I).EQ.0D0) GO TO 11

IF(Z.EQ.0D0) GO TO 134

С Подпрограмма, реализующая итерационный процесс определения

С плотности из уравнения состояния (метод Ньютона)

IF(ITER.GT.10) GO TO 4

IF(DABS(DR/X).GT.1.D-6) GO TO 1

С * Подпрограмма расчета коэффициента сжимаемости природного *

С * газа по уравнению состояния ВНИЦ СМВ. *

IF(ICALC.NE.1) GO TO 1

IF(Z.EQ.0D0) GO TO 133

IF(NPR.EQ.1) GO TO 333

IF(NPR.EQ.0) GO TO 555

IF(YI(I).EQ.0D0) GO TO 155

IF(Z.EQ.0D0) GO TO 134

С Подпрограмма, реализующая итерационный процесс определения

С плотности из уравнения состояния (метод Ньютона)

IF(ITER.GT.10) GO TO 4

IF(DABS(DR/X).GT.1.D-6) GO TO 1

4 CALL COMPL(X,T,NPRIZ)

IF(NPRIZ.NE.0) GO TO 7

ПРИЛОЖЕНИЕ Г

Примеры расчета коэффициента сжимаемости природного газа

Г.1 Модифицированный метод NX19

Плотность при 0,101325 МПа и 293,15 К: 0,6799 кг/м 3

Г.2 Уравнение состояния GERG-91

Плотность при 0,101325 МПа и 293,15 К: 0,6799 кг/м 3

Г.3 Уравнение состояния AGA8-92DC

Состав природного газа в молярных процентах:

Плотность при 0,101325 МПа и 293,15 К: 0,6799 кг/м 3

Г.4 Уравнение состояния ВНИЦ СМВ

Состав природного газа в молярных процентах:

Плотность при 0,101325 МПа и 293,15 К: 0,7675 кг/м3

ПРИЛОЖЕНИЕ Д

Влияние погрешности исходных данных на погрешность расчета коэффициента сжимаемости природного газа (примеры расчета)

Д.1 Модифицированный метод NX19

Исходные данные (заданные параметры)

Плотность, кг/м 3 (0,101325 МПа, 293,15 К)

диоксида углерода (СО₂)

Д.2 Уравнение состояния GERG-91

Исходные данные (заданные параметры)

Плотность, кг/м 3 (0,101325 МПа, 293,15 К)

диоксида углерода (СО₂)

Д.3 Уравнение состояния AGA8-92DC

Исходные данные (заданные параметры)

диоксида углерода (СО₂)

Д.4 Уравнение состояния ВНИЦ СМВ

Исходные данные (заданные параметры)

диоксида углерода (СО₂)

ПРИЛОЖЕНИЕ Е

Библиография

Ключевые слова: природный газ, методы расчета коэффициента сжимаемости, давление, температура, плотность при стандартных условиях, компонентный состав, молярные и объемные доли, коэффициент сжимаемости, фактор сжимаемости, плотность, погрешность, уравнение состояния, итерационный процесс, листинг программы

Источник