"ГОСТ Р 8.662-2009 (ИСО 20765-1:2005). Национальный стандарт Российской Федерации. Государственная система обеспечения единства измерений. Газ природный. Термодинамические свойства газовой фазы. Методы расчетного определения для целей транспортирования и распределения газа на основе фундаментального уравнения состояния AGA8"
(утв. и введен в действие Приказом Ростехрегулирования от 15.12.2009 N 1043-ст)

Цели и принципы стандартизации в Российской Федерации установлены Федеральным законом от 27 декабря 2002 г. N 184-ФЗ "О техническом регулировании", а правила применения национальных стандартов Российской Федерации - ГОСТ Р 1.0-2004 "Стандартизация в Российской Федерации. Основные положения".

1. Подготовлен Федеральным государственным унитарным предприятием "Российский научно-технический центр информации по стандартизации, метрологии и оценке соответствия" (ФГУП "СТАНДАРТИНФОРМ") на основе собственного аутентичного перевода на русский язык стандарта, указанного в пункте 4.

2. Внесен Техническим комитетом по стандартизации ТК 180 "Государственная служба стандартных справочных данных".

3. Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 15 декабря 2009 г. N 1043-ст.

4. Настоящий стандарт является модифицированным по отношению к международному стандарту ИСО 20765-1:2005 "Газ природный. Расчетное определение термодинамических свойств. Часть 1. Свойства газовой фазы для целей транспортирования и распределения газа" (ISO 20765-1:2005 "Natural gas - Calculation of thermodynamic properties - Part 1: Gas phase properties for transmission and applications") путем изменения его структуры для приведения в соответствие с правилами, установленными в ГОСТ Р 1.5 (подразделы 4.2 и 4.3), а также путем изменения отдельных фраз (слов, значений, показателей, ссылок), которые выделены в тексте полужирным курсивом.

Внесение указанных технических отклонений направлено на учет особенностей объекта стандартизации, характерных для Российской Федерации, а также связано с целесообразностью использования ссылочных национальных и межгосударственных стандартов вместо ссылочных международных стандартов.

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5-2004 (подраздел 3.5).

Сравнение структуры настоящего стандарта со структурой указанного международного стандарта приведено в дополнительном Приложении ДА.

5. Введен впервые.

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

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

Стандарт распространяется на подготовленные для транспортирования по магистральным газопроводам газы в диапазонах давления p и температуры t, при которых на практике осуществляют транспортирование и распределение газов.

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

ГОСТ 8.417-2002 Государственная система обеспечения единства измерений. Единицы величин

ГОСТ 31369-2008 Газ природный. Вычисление теплоты сгорания, плотности, относительной плотности, числа Воббе на основе компонентного состава (ИСО 6976-1:1995, MOD)

ГОСТ 31371.1-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 1. Руководство по проведению анализа (ИСО 6974-1:2000, MOD)

Примечание. При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодно издаваемому информационному указателю "Национальные стандарты", который опубликован по состоянию на 1 января текущего года, и по соответствующим ежемесячно издаваемым информационным указателям, опубликованным в текущем году. Если ссылочный стандарт заменен (изменен), то при пользовании настоящим стандартом следует руководствоваться заменяющим (измененным) стандартом. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, применяется в части, не затрагивающей эту ссылку.

3.1. Термины и определения

В настоящем стандарте применены термины по ГОСТ 8.417, ГОСТ 31371.1, с учетом [1] - [6], а также следующие термины с соответствующими определениями:

3.1.1. Неидеальная составляющая свойства (residual property): составляющая термодинамического свойства, которая определяется неидеальным (реальным) поведением газа или гомогенной газовой смеси, то есть представляет собой разницу между термодинамическим свойством реального газа или газовой смеси и идеально-газовой составляющей этого свойства газа или газовой смеси при тех же самых температуре и плотности.

3.1.2. Калорическое свойство (caloric property): характеристика газа или гомогенной газовой смеси, которую можно рассчитать по фундаментальному уравнению состояния (с привлечением данных об идеально-газовых значениях теплоемкости, энтальпии и энтропии).

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

3.1.3. Термическое свойство (volumetric property): характеристика газа или гомогенной газовой смеси, которые можно рассчитать по термическому уравнению состояния.

Примечание. Настоящий стандарт распространяется на такие термические свойства, как фактор сжимаемости и плотность; перевод термина дан в соответствии с [6].

3.1.4. Термодинамическое свойство (thermodynamic property): термическое или калорическое свойство.

3.1.5. Уравнение состояния (equation of state): математическое выражение взаимосвязи между параметрами состояния газа или гомогенной газовой смеси.

Примечание. В этой части следует учитывать различия между двумя видами уравнения состояния, а именно: 1) термическим уравнением состояния, которое устанавливает взаимосвязь между такими параметрами состояния, как давление, температура и объем, занимаемый заданным количеством вещества; 2) фундаментальным уравнением состояния, которое устанавливает взаимосвязь между плотностью, температурой и свободной энергией Гельмгольца.

3.2. Обозначения

3.2.1. Условные обозначения

Условные обозначения величин приведены в таблице 1.

3.2.2 Подстрочные индексы

Индексы в условных обозначениях величин означают следующее:

i - для i-го компонента (от i = 1 до N);

j - для второго (j-го) компонента в бинарном взаимодействии (от j = 2 до N);

n - номер коэффициента в уравнении состояния (от n = 1 до 58);

0 - для идеально-газового состояния;

r - для неидеальной составляющей;

- частная производная по относительной молярной плотности;

- для стандартного состояния (, );

- частная производная по обратной относительной температуре.

4.1. Сущность метода

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

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

В настоящем методе используют результаты детального анализа молярного состава смеси, то есть должны быть указаны все компоненты, присутствующие в смеси в количествах, превышающих 0,00005 молярной доли (50 миллионных частей молярной доли) <*>. Для типичного природного газа такой анализ должен включать алканы, примерно до или , совместно с азотом, диоксидом углерода и гелием. Обычно изомеры алканов свыше объединяют с нормальными изомерами с учетом молярной массы.

--------------------------------

<*> - несистемная единица.

Для некоторых природных газов следует учитывать дополнительные компоненты, такие как углеводороды и , водяной пар и сероводород. Для коммунально-бытовых газов необходимо учитывать наличие водорода и окиси углерода.

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

4.2. Фундаментальное уравнение состояния для свободной энергии Гельмгольца

4.2.1. Исходные положения

Примечание. Исходные положения для фундаментального уравнения состояния для свободной энергии Гельмгольца приведены в Приложении A.

4.2.2. Свободная энергия Гельмгольца

Свободную энергию Гельмгольца f гомогенной газовой смеси при заданных давлении и температуре выражают как сумму идеально-газовой составляющей и неидеальной составляющей , характеризующей поведение реального газа, согласно уравнению

которое при безразмерном представлении свободной энергии в виде преобразуют в уравнение

где X - вектор молярных долей смеси;

- обратная (безразмерная) относительная температура, связанная с температурой T по соотношению

где L = 1 K.

Поскольку, согласно статистической термодинамике, свободная энергия Гельмгольца представляет собой непосредственный результат числа и типов молекулярных взаимодействий в смеси и является тем самым непосредственной функцией молярной плотности и молярных долей молекул в смеси, в уравнения (1) и (2) в качестве входных переменных записывают соответственно - молярную плотность и - относительную плотность, а не давление p.

Относительную плотность связывают с молярной плотностью по уравнению

где K - смесевой параметр размера.

Идеально-газовую составляющую относительной свободной энергии Гельмгольца получают из уравнений для идеально-газовой изобарной теплоемкости (4.2.2) и избыточной составляющей из уравнения состояния AGA8 (4.2.3).

4.2.3. Идеально-газовая составляющая свободной энергии Гельмгольца

Идеально-газовую составляющую свободной энергии Гельмгольца через энтальпию и энтропию выражают в соответствии с уравнением

В свою очередь, энтальпию и энтропию выражают через идеально-газовую изобарную теплоемкость по уравнениям (6) и (7), где в качестве пределов интегрирования устанавливают температуры и T:

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

Неидеальную составляющую свободной энергии Гельмгольца определяют, используя уравнения (6) и (7), как функцию , и X, по уравнению

(более подробную информацию см. в Приложении B).

4.2.4. Неидеальная составляющая свободной энергии Гельмгольца

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

где B - второй вириальный коэффициент;

- функция молярных долей компонентов;

, , , - коэффициенты уравнения состояния и функции молярных долей компонентов.

Фактор сжимаемости Z связывают с неидеальной составляющей относительной свободной энергии Гельмгольца согласно уравнению

где - частная производная от по относительной плотности при постоянных и X.

Путем исключения Z, используя уравнения (9) и (10) и интегрирование по относительной плотности, приходят к уравнению для неидеальной составляющей относительной свободной энергии Гельмгольца

(более подробную информацию см. в Приложении C).

4.2.5. Относительная свободная энергия Гельмгольца

Исходя из возможности рассчитать все термодинамические свойства по фундаментальному уравнению состояния (2) для относительной свободной энергии Гельмгольца , путем использования уравнения (8) для идеально-газовой составляющей и уравнения (11) для неидеальной составляющей , относительную свободную энергию Гельмгольца представляют в виде уравнения

4.3. Термодинамические свойства, получаемые из свободной энергии Гельмгольца

4.3.1. Исходные положения

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

Каждую производную определяют в виде суммы идеально-газовой составляющей (Приложение B) и неидеальной составляющей (Приложение C). Вводят подстановки, приведенные в уравнениях, для упрощения получения требуемых аналитических выражений:

(более детальные выражения для , , , и см. в Приложении C; необходимые общие выражения для различных термодинамических свойств приведены в 4.3.2.1 - 4.3.2.9 [уравнения (17) - (26)]: в уравнениях (19) - (24) нижние символы относятся к молярным величинам (то есть отнесенным к 1 киломолю), и соответствующие верхние символы относятся к массовым величинам (то есть отнесенным к 1 килограмму); переход от молярных переменных к массовым достигается делением на молярную массу M).

Примечание. В этих уравнениях R - молярная газовая постоянная; соответственно R/M - массовая газовая постоянная.

Молярную массу M смеси получают на основе компонентного состава X и молярных масс компонентов по уравнению

Значения молярных масс компонентов приведены в [7] и [8], в которых эти значения являются идентичными значениям, приведенным в ГОСТ 31369.

Примечание. Значения, приведенные для молярных масс, в большинстве случаев не идентичны современным значениям, применяемым в международной метрологической практике. Они представляют собой значения, которые были использованы в основном при получении уравнения состояния AGA8 и в этой связи оставлены без изменений; расхождения во всех случаях не превышают 0,001 кг/кмоль.

В уравнениях (20), (21) и (23) - (26) основные выражения для свойств h, s, , , k и w представляют в нескольких вариантах для того, чтобы полученные значения свойств можно было использовать для упрощения последующих вычислений. Такой подход можно применять в случаях, когда необходимо определять несколько или все термодинамические свойства. В каждом подпункте сначала представлены основные термодинамические соотношения, а далее вспомогательные выражения.

4.3.2. Уравнения для термодинамических свойств

4.3.2.1. Фактор сжимаемости и плотность

Выражение для фактора сжимаемости Z представляют согласно уравнению

где - производная от свободной энергии Гельмгольца по относительной молярной плотности [см. также уравнение (10)]. Молярную плотность и удельную (массовую) плотность D связывают с давлением согласно уравнению

Значения фактора сжимаемости Z, рассчитываемые по настоящему стандарту, как правило, не должны отличаться от значений, рассчитываемых по стандарту [1]. Для этого выход из итерационного процесса при расчете плотности на шаге 4 блок-схемы расчета термодинамических свойств (рисунок F.1, Приложение F) должен проводиться по критерию разности абсолютных значений рассчитанного и заданного (входного) давления; эта разность должна быть меньше .

4.3.2.2. Внутренняя энергия

Выражение для внутренней энергии u представляют согласно уравнению

4.3.2.3. Энтальпия

Выражение для энтальпии h представляют согласно уравнению

4.3.2.4. Энтропия

Выражение для энтропии s представляют согласно уравнению

4.3.2.5. Изохорная теплоемкость

Выражение для изохорной теплоемкости представляют согласно уравнению

4.3.2.6. Изобарная теплоемкость

Изобарную теплоемкость определяют по уравнению

4.3.2.7. Коэффициент Джоуля-Томсона

Выражение для коэффициента Джоуля-Томсона представляют уравнением

4.3.2.8. Показатель адиабаты

Выражение для показателя адиабаты k представляют согласно уравнению

4.3.2.9. Скорость звука

Выражение для скорости звука w представляют согласно уравнению

5.1. Входные переменные

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

Переход от давления к относительной плотности проводят по 5.2.

Если вместо давления p в качестве входной переменной приводят удельную (массовую) плотность D, значения получают непосредственно без проведения процедур по 5.2, а именно как , где M - молярная масса, определяемая по уравнению (16). Компонентный состав необходимо определять для следующего 21 компонента: азота, диоксида углерода, метана, этана, пропана, н-бутана, изобутана (2-метилпропана), н-пентана, изопентана (2-метилбутана), н-гексана, н-гептана, н-октана, н-нонана, н-декана, водорода, кислорода, углерода, воды, сероводорода, гелия и аргона (о допустимых пределах молярных долей компонентов, устанавливаемых рассматриваемым методом, - см. 6.2). Какой-либо следовой компонент, присутствующий в газе, но не идентифицируемый в числе перечисленного 21 компонента, можно приблизительно учитывать в совокупности с подходящим из перечисленных выше (Приложение E). Сумма молярных долей компонентов должна быть равна единице.

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

Примечание. Если состав газа известен в объемных долях, его необходимо пересчитать в молярные доли, используя метод по ГОСТ 31369.

5.2. Переход от давления к относительной плотности

Комбинируя уравнения (4), (9) и (18), получают уравнение

Если входные переменные выражены через давление, обратную относительную температуру и компонентный состав, используя уравнение (27), можно определить приведенную молярную плотность . Значения параметров , , K(X) и коэффициентов , , и уравнения (27) можно получить из уравнений и таблиц, представленных в Приложении D [уравнений (D.1), (D.6) и (D.11), и таблицы D.1 соответственно] для этих характеристик. Численные значения параметров для всех индивидуальных компонентов и их бинарного взаимодействия, также требуемые для решения уравнений (D.1), (D.6) и (D.11), представлены в таблицах D.2 и D.3 соответственно.

Решение можно получить подходящим численным методом. На практике применяют стандартный алгоритм решения уравнения состояния относительно плотности как наиболее удобный и достаточно корректный. В таких алгоритмах обычно используют начальное приближение для значения плотности (часто идеально-газовое приближение) и продолжают расчеты путем итерационных вычислений p и для того, чтобы найти значение , которое воспроизводит известное значение p в пределах заданного уровня точности. В рассматриваемом методе в качестве допустимого критерия устанавливают условие, что давление, рассчитываемое из вычисленных значений относительной молярной плотности , должно воспроизводить исходное значение p в пределах .

5.3. Проведение основных вычислений

Для проведения вычислений в целом необходимо выполнение нижеследующих условий.

Требуемый набор входных переменных имеется в наличии. С этим выверенным набором входных переменных: относительной плотностью , обратной относительной температурой и компонентным составом X используют фундаментальное уравнение состояния для вычисления относительной свободной энергии Гельмгольца и других термодинамических свойств. Согласно уравнению (12) определяют относительную свободную энергию Гельмгольца как . Согласно уравнению (11) определяют неидеальную составляющую свободной энергии Гельмгольца как функцию относительной плотности , обратной относительной температуры и молярного компонентного состава X. Идеально-газовую составляющую , определяемую уравнением (8), получают согласно уравнению (B.3) Приложения B таким образом, чтобы выразить как представлено в уравнении

Значения всех коэффициентов , и от до для идеального газа даны в Приложении B для всех из 21 возможного компонента газа.

Производные от по (относительной) плотности и (обратной приведенной) температуре, которые необходимы для вычислений различных термодинамических свойств, можно получить из уравнений (C.2) - (C.6) (Приложение C). Наконец, различные термодинамические свойства можно рассчитать по уравнениям (17) - (26). Значения коэффициентов , , и и параметров , являющихся функциями молярных долей компонентов, приведены в Приложении D.

Более подробное изложение процедур вычислений представлено в Приложении F.

6.1. Давление и температура

Метод настоящего стандарта применим только к газам, подготовленным для транспортирования по магистральным газопроводам (6.2), в диапазонах давлений и температур, при которых обычно осуществляют транспортирование и распределение газа. Соответствующие диапазоны давлений и температур представлены в таблице 2. Настоящий метод применим только к смесям в газообразном состоянии, и на него распространяется специальное ограничение, что он не должен применяться для каких-либо условий, когда расчетное значение фактора сжимаемости может быть менее 0,5.

6.2. Газ, подготовленный для транспортирования по магистральным газопроводам

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

Возможные следовые компоненты природного газа и конкретные рекомендации по их учету - см. Приложение F. Общее суммарное содержание следовых компонентов не должно превышать 0,0005 молярной доли.

7.1. Неопределенности расчета свойств для газа, подготовленного для транспортирования по магистральным газопроводам

7.1.1. Диаграммы неопределенностей для метана

В качестве примера для неопределенностей, которые можно ожидать для смесей, содержащих молярную долю метана, близкую к единице, на рисунках 1 - 3 представлены соответственно диаграммы неопределенностей для фактора сжимаемости, скорости звука и энтальпии чистого метана. Неопределенность, выраженная с 95%-ной доверительной вероятностью, данная для каждой области, представляет наибольшее значение в пределах суммы неопределенности в заслуживающих доверия справочных данных и разницы между справочными данными и значениями свойства, рассчитанного методом настоящего стандарта. Использованные справочные данные рассчитаны по уравнениям, приведенным в [10].

7.1.2. Диаграммы неопределенностей для природного газа

В качестве примера неопределенностей, которые можно ожидать для природного газа, на рисунках 4 - 6 представлены соответственно диаграммы неопределенностей для фактора сжимаемости, скорости звука и энтальпии. Неопределенность, выраженная с 95%-ной доверительной вероятностью, дана для каждой области наибольших отклонений (для широкого спектра природных газов) между измеренным значением свойств и значением, рассчитанным с использованием метода настоящего стандарта. Для фактора сжимаемости экспериментальные значения взяты из [10]; для скорости звука - из [11]; а для энтальпии - из [12].

Примечание 1. Для всех газов диаграмма неопределенности для плотности идентична по форме диаграмме для фактора сжимаемости.

Примечание 2. Для всех газов неопределенность для показателя адиабаты приблизительно в два раза больше неопределенности для скорости звука.

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

Для газов при низких давлениях (ниже, как правило, 1 МПа), ведущих себя подобно идеальному газу (фактор сжимаемости более чем 0,95), можно ожидать, что все калорические свойства должны быть получены с малой неопределенностью. Это следует из того факта, что в рассматриваемом частном случае большая часть каждого свойства определяется из идеально-газовой части свободной энергии Гельмгольца, которая будучи рассчитана прямо из высокоточных данных по изобарной теплоемкости идеального газа (Приложение B) сама является высокоточной. В таком случае неопределенности для плотности, фактора сжимаемости, скорости звука, изохорной и изобарной теплоемкости, показателя адиабаты и коэффициента Джоуля-Томпсона находятся в пределах 0,1%.

7.2. Влияние неопределенностей входных переменных

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

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

Для проверки не следует использовать лишние цифры (см. пример в Приложении G).

Уравнение состояния [8] AGA8 было опубликовано в 1992 г. Комитетом по измерениям при транспортировке газа Американской газовой ассоциации как метод для высокоточных вычислений фактора сжимаемости. В этой связи оно является объектом рассмотрения. Для вычислений на основе уравнения состояния AGA8 всех термодинамических свойств необходимо выполнить два основных требования:

a) уравнение должно быть математически преобразовано к виду, где избыточная свободная энергия Гельмгольца представляется в явном виде, то есть в виде фундаментального уравнения состояния для свободной энергии Гельмгольца. На основе соответствующих соотношений [8] можно рассчитывать все избыточные термодинамические свойства из избыточной свободной энергии Гельмгольца и ее производных по температуре и плотности как параметров состояния;

b) для расчета калорических свойств необходимо наличие идеально-газовой составляющей свободной энергии Гельмгольца как функции температуры. В этом случае целесообразно использовать соответствующие соотношения [13], [14] для идеально-газовой изобарной теплоемкости. При этом также требуются производные свободной энергии Гельмгольца по параметрам состояния.

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

Предлагаемый метод вычислений следует применять в прикладных программах, в частности в программах, которые можно использовать при транспортировании и распределении газа.

B.1. Вычисление идеально-газовой составляющей свободной энергии Гельмгольца

a) Идеально-газовую изобарную теплоемкость индивидуального компонента по [13], [15] можно записать в следующем виде

Уравнение (B.1) можно обобщить на случай N-компонентной смеси, используя следующее уравнение

b) Это уравнение для идеально-газовой изобарной теплоемкости можно включить в уравнение (8) для относительной свободной энергии Гельмгольца . В результате получают уравнение

В уравнении (B.3) константы и связывают с константами интегрирования и уравнения (8) в соответствии с уравнениями:

Важно отметить, что в уравнении (B.3) является функцией относительной молярной плотности реального газа, а не относительной молярной плотности идеального газа и, следовательно, не может быть полностью рассчитана до тех пор, пока значение неизвестно (см. Приложение D). Уравнение (B.3) используют для расчета относительной свободной энергии Гельмгольца идеального газа.

c) Наиболее достоверные из известных в научной литературе данные об идеально-газовой изобарной теплоемкости [14] были использованы в качестве исходных данных для получения численных значений констант с до для каждого из 21 индивидуального компонента, рассматриваемого в настоящем стандарте.

Значения констант и , а также констант с до , используемые в уравнении (B.3), приведены в таблице B.1.

B.2. Производные от идеально-газовой составляющей свободной энергии Гельмгольца

Для вычисления ряда термодинамических свойств необходимо иметь первые и вторые частные производные от идеально-газовой составляющей свободной энергии Гельмгольца [уравнение (B.3)] по обратной относительной температуре . Необходимые для этого математические выражения представляют в виде уравнений:

C.1. Вычисление свободной энергии Гельмгольца

Относительную свободную энергию Гельмгольца для газовой смеси заданного состава рассчитывают по уравнению

где рассчитывают в соответствии с Приложением B. Ту часть уравнения (C.1), которая учитывает неидеальную составляющую свободной энергии Гельмгольца (вся правая часть уравнения, за исключением ) определяют из уравнения (11). Уравнения для расчета функций , и K(X) в уравнении (C.1) приведены в D.1. Значения различных констант в уравнении (C.1) даны в D.2.

C.2. Производные от свободной энергии Гельмгольца

Для вычисления совокупности термодинамических свойств требуются первые и/или вторые частные производные по обратной относительной температуре и/или по относительной плотности от относительной свободной энергии Гельмгольца. Необходимые для этого математические выражения записывают в виде уравнений:

Уравнения и константы, необходимые для расчета по уравнениям (C.2) - (C.6), приведены в Приложении D.

D.1. Вспомогательные уравнения

Для расчета по уравнению (C.1) необходимо определить значения второго вириального коэффициента функций и параметра размера K(X). Кроме того, так как уравнение (C.1) включает в качестве аргумента относительную молярную плотность , необходимо получить ее значение при задании давления в качестве входной переменной. Эти шаги рассматриваются в перечислениях a) - d):

a) второй вириальный коэффициент рассчитывают по уравнениям (D.1) - (D.5) с использованием необходимых констант из таблиц D.1 - D.3:

где

где

b) функции (для n от 13 до 58) рассчитывают по уравнениям (D.6) - (D.10) с использованием необходимых констант из таблиц D.1 - D.3:

где

В таблице D.2 равны нулю для всех компонентов, за исключением водорода, для которого , а равны нулю для всех компонентов, за исключением воды, для которой . Многие из параметров бинарного взаимодействия в таблице D.3 равны единице;

c) смесевой параметр размера K(X) рассчитывают по уравнению (D.11) с использованием необходимых констант из таблиц D.2 и D.3

Многие значения равны единице;

d) относительную молярную плотность определяют в результате решения уравнения (D.12) (это уравнение получено из уравнения (27) настоящего стандарта, преобразованного к явному виду для давления)

Уравнение (D.12) решают с использованием стандартных алгоритмов расчета плотности по уравнению состояния.

D.2. Значения констант для уравнения состояния

Настоящий раздел содержит значения всех констант, необходимых для использования уравнения состояния AGA8 в уравнениях (C.1) и (D.1) - (D.12). В таблице D.1 даны значения констант самого уравнения состояния. В таблице D.2 представлены значения характеристических параметров, определяющих свойства индивидуальных компонентов смеси. В таблице D.3 даны значения параметров бинарного взаимодействия между молекулами различных компонентов смеси.

Параметры взаимодействия для любой пары компонентов, не приведенные в таблице D.3, равны 1,0.

Для того чтобы рассчитать, используя метод, устанавливаемый в настоящем стандарте, термодинамические свойства природного газа или подобной смеси, которая содержит следы одного или более компонентов, не включенных в таблицу 3, необходимо включить каждый такой следовой компонент в один из 21 основного и второстепенного компонента, для которого было разработано уравнение AGA8. Рекомендации по такому включению даны в таблице E.1.

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

Примечание. Набор следовых компонентов по таблице E.1 и таблице 3 настоящего стандарта соответствует ГОСТ 31369.

Настоящее Приложение представляет общую последовательность действий, необходимых для вычисления термодинамических свойств в соответствии с настоящим стандартом. Метод достаточно сложный, поэтому для его применения необходимо использовать компьютерные программы. В настоящем Приложении обсуждаются вычислительные аспекты метода. На рисунке F.1 приведена блок-схема процесса вычислений термодинамических свойств с использованием компьютерных программ.

Шаг 1

Для расчета необходимо ввести молярную долю каждого компонента природного газа, абсолютную температуру T и абсолютное давление p.

В рассматриваемом методе температуру задают в градусах Кельвина (K), а давление в мегапаскалях (МПа). Если используют какие-либо другие единицы измерения, их нужно перевести в градусы Кельвина или мегапаскали. Для этого необходимо использовать переводные коэффициенты на основе положений ГОСТ 8.417 с учетом стандарта [3].

Иногда вместо давления p в качестве входной величины используют массовую плотность D. В таком случае все вычисления выполняют как обычно, но с учетом содержания последней ступени шага 4.

Шаг 2

a) Для заданных значений , используя уравнения (D.2) - (D.5), рассчитывают значения .

b) Таким же образом для заданных значений , используя уравнения (D.6) - (D.10), рассчитывают значения .

c) По уравнению (D.11) рассчитывают значения смесевого параметра размера K(X).

d) Рассчитывают следующие значения:

- относительной температуры по уравнению (3) при ;

- относительной плотности идеального газа при стандартных условиях по уравнению (4) при .

Плотность идеального газа при стандартных условиях определена в пункте 4.2.3.

Шаг 3

a) Для заданного значения T по уравнению (3) рассчитывают , а затем, используя значения [шаг 2, перечисление a)], по уравнению (D.1) рассчитывают значения второго вириального коэффициента .

b) Таким же образом, для заданных значений T, используя значения [шаг 2, перечисление b)], рассчитывают значения для n от 13 до 58.

Шаг 4

a) Подставляют рассчитанные значения , и K(X) [шаг 3, перечисления a), b) и шаг 2, перечисление c)] соответственно в уравнение (D.12). Это уравнение для p, которое было задано в шаге 1, теперь содержит одно неизвестное - относительную плотность .

b) Решают это уравнение относительно . Требуемое решение может быть получено с использованием подходящего численного метода, но на практике наиболее подходящим алгоритмом может быть стандартный алгоритм нахождения плотности из уравнения состояния. Такие алгоритмы обычно используют для первоначальной оценки плотности (часто в приближении идеального газа) и итеративным путем находят такое значение , которое воспроизводит значение p с заранее заданной точностью. Подходящим критерием в настоящем случае является тот, при котором рассчитанное давление при заданной молярной плотности отличается от заданного значения давления p менее чем на .

Если в качестве входной величины вместо давления p используют массовую плотность D, то вычисляют непосредственно без итераций как , где M(X) - молярная масса, рассчитанная из уравнения (16) настоящего стандарта.

Шаг 5

Подставляют и в уравнение (B.3) вместе с различными константами и функциями обратной относительной температуры для расчета идеально-газовой части относительной свободной энергии Гельмгольца. Используют также уравнения (B.6) и (B.7) для расчета первой и второй частных производных по обратной относительной температуре.

Шаг 6

Используют уравнение (B.1) для расчета общей относительной свободной энергии Гельмгольца . Рассчитывают по уравнениям (C.2) и (C.4) первые частные производные и по обратной относительной температуре и относительной плотности соответственно. Уравнение (C.3) используют для расчета второй частной производной по обратной относительной температуре, уравнение (C.5) - для расчета и уравнение (C.6) - для расчета .

Шаг 7

a) Используют уравнения (17) - (26) настоящего стандарта для расчета полного набора термодинамических свойств: фактора сжимаемости Z; плотности ; внутренней энергии u; энтальпии h; энтропии s; изохорной теплоемкости ; изобарной теплоемкости ; коэффициента Джоуля-Томсона ; показателя адиабаты k и скорости звука w.

b) На основе уравнений (18) - (23) можно представлять термодинамические свойства, рассчитанные на основе значений молярной (при использовании символов нижнего регистра) или массовой (при использовании символов верхнего регистра) доли. Для преобразования молярных величин в массовые используют молярную массу смеси M(X), которую рассчитывают по уравнению (16).

Шаг 8

Возврат на шаг 1 (выбор новых входных переменных) или окончание вычислений.

Следующие примеры расчетов приведены для целей проверки программных решений (таблицы G.1 - G.7).