Межгосударственный стандарт ГОСТ 30319.3-2015 "Газ природный. Методы расчета физических свойств. Вычисление физических свойств на основе данных о компонентном составе" (введен в действие приказом Федерального агентства по техническому регулированию и метрологии от 30 ноября 2015 г. N 2075-ст)

Natural gas. Methods of calculation of physical properties. Calculation of physical properties on base information on component composition

Дата введения - 1 января 2017 г.

Взамен 30319.3-96

Предисловие

Цели, основные принципы и порядок проведения работ по межгосударственной стандартизации установлены в ГОСТ 1.0-92 "Межгосударственная система стандартизации. Основные положения" и ГОСТ 1.2-2009 "Межгосударственная система стандартизации. Стандарты межгосударственные, правила и рекомендации по межгосударственной стандартизации. Правила разработки, принятия, применения, обновления и отмены"

Сведения о стандарте

1 Разработан Обществом с ограниченной ответственностью "Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ", Техническим комитетом по стандартизации ТК 52 "Природный и сжиженные газы"

2 Внесен Федеральным агентством по техническому регулированию и метрологии (Росстандарт)

3 Принят Межгосударственным советом по стандартизации, метрологии и сертификации (протокол от 27 августа 2015 г. N 79-П)

За принятие проголосовали:

Краткое наименование страны по МК (ИСО 3166) 004-97	Код страны по МК (ИСО 3166) 004-97	Сокращенное наименование национального органа (по управлению строительством) по стандартизации
Армения Беларусь Казахстан Киргизия Молдова Россия Украина	AM BY KZ KG MD RU UA	Минэкономики Республики Армения Госстандарт Республики Беларусь Госстандарт Республики Казахстан Кыргызстандарт Молдова-Стандарт Росстандарт Минэкономразвития Украины

4 Приказом Федерального агентства по техническому регулированию и метрологии от 30 ноября 2015 г. N 2075-ст межгосударственный стандарт ГОСТ 30319.3-2015 введен в действие в качестве национального стандарта Российской Федерации с 1 января 2017 г.

5 Взамен 30319.3-96

1 Область применения

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

1.2 Настоящий стандарт применяют для расчета указанных в 1.1 физических свойств природного газа при давлениях до 30 МПа включительно и температурах от 250 до 350 К.

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

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

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

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

ГОСТ 31371.2-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 2. Характеристики измерительной системы и статистические оценки данных

ГОСТ 31371.3-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 3. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов до с использованием двух насадочных колонок

ГОСТ 31371.4-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 4. Определение азота, диоксида углерода и углеводородов и в лаборатории и с помощью встроенной измерительной системы с использованием двух колонок

ГОСТ 31371.5-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 5. Определение азота, диоксида углерода и углеводородов и в лаборатории и при непрерывном контроле с использованием трех колонок

ГОСТ 31371.6-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 6. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов с использованием трех капиллярных колонок

ГОСТ 31371.7-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 7. Методика выполнения измерений молярной доли компонентов

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

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

3 Термины, определения и обозначения

3.1 В настоящем стандарте применены термины и определения по ГОСТ 30319.1.

3.2 Основные условные обозначения величин, принятые в стандарте, приведены в таблице 2 ГОСТ 30319.1.

4 Методы расчета физических свойств природного газа

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

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

,

(1)

где - приведенное давление;

- приведенная температура;

- безразмерный комплекс (см. 4.1.1.2).

4.1.1.1 Приведенные давление () и температуру () рассчитывают по формулам:

,

(2)

,

(3)

где - параметр приведения для давления, МПа;

- параметр приведения для температуры, равный 1 К.

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

,

(4)

где - смесевой параметр размера, ;

R - универсальная газовая постоянная (см. таблицу 2 ГОСТ 30319.1).

Смесевой параметр размера () рассчитывают по формуле

,

(5)

где Nc - число компонентов природного газа;

{} и {} - параметры размера компонентов природного газа, значения которых приведены в таблице А.1 (приложение А);

{} - параметры бинарного взаимодействия, значения которых приведены в таблице А.2 (приложение А).

4.1.1.2 Безразмерный комплекс () рассчитывают по формуле

,

(6)

где {}, {}, {}, {}, и {} - коэффициенты и показатели степеней, значения которых приведены в таблице А.3 (приложение А);

{}, {} - функции молярных долей компонентов природного газа.

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

(7)

.

(8)

Вспомогательные функции ( и ) рассчитывают по формулам:

,

(9)

,

(10)

,

(11)

,

(12)

,

(13)

,

(14)

,

(15)

,

(16)

,

(17)

где {}, {}, {}, {}, {} - параметры, значения которых приведены в таблице А.3 (приложение А);

{}, {}, {}, {}, {}, {} - параметры компонентов природного газа, значения которых приведены в таблице А.1 (приложение А);

{}, {}, {} - параметры бинарного взаимодействия компонентов природного газа, значения которых приведены в таблице А.2 (приложение А).

4.1.2 Решение уравнения (1) осуществляют в итерационном процессе методом Ньютона; значение начального приближения () рассчитывают, используя заданные значения температуры, давления и молярных долей природного газа (см. 5.2.2). После вычисления приведенной плотности () в итерационном процессе плотность смеси рассчитывают по формуле

,

(18)

где - молярная масса смеси, кг/кмоль;

- смесевой параметр размера (см. формулу (5)).

Молярную массу смеси рассчитывают по формуле

,

(19)

где - молярная масса i-го компонента природного газа, значения которой для каждого компонента приведены в таблицах А.1 и А.9 (приложение А);

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

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

,

(20)

где - безразмерный комплекс (см. 4.1.1.2).

Примечание - Безразмерный комплекс в формуле (20) рассчитывают при заданных значениях (Т, {}) и найденном в результате решения уравнения (1) значении приведенной плотности ().

4.2 Методы расчета показателя адиабаты и скорости звука

4.2.1 Показатель адиабаты и скорость звука рассчитывают по следующим формулам:

,

(21)

,

(22)

где , и - безразмерные комплексы (см. 4.2.2);

- безразмерная изобарная теплоемкость природного газа в идеально-газовом состоянии (см. 4.2.3).

4.2.2 Безразмерные комплексы , , и рассчитывают по следующим формулам:

(23)

,

(24)

.

(25)

Коэффициенты, показатели степеней, параметры, функции, входящие в формулы расчета (23) - (25), те же самые, которые входят в формулы расчета безразмерного комплекса (см. 4.1.1.2).

4.2.3 Безразмерную изобарную теплоемкость природного газа в идеально-газовом состоянии () рассчитывают по формуле

,

(26)

где {} - безразмерные изобарные теплоемкости компонентов природного газа в идеально-газовом состоянии;

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

Значения величин {} рассчитывают по формуле

,

(27)

где .

Коэффициенты {}, {}, {}, {}, {}, {}, {}, {}, {}, формулы (27) приведены в таблице А.4 (приложение А).

4.3 Метод расчета коэффициента динамической вязкости

4.3.1 Вязкость природного газа рассчитывают по формуле

,

(28)

где - вязкость природного газа в разреженном состоянии;

- молярная масса природного газа (см. формулу (19));

- псевдокритическое давление природного газа;

- псевдокритическая температура природного газа (см. формулу (37));

- избыточная составляющая вязкости природного газа.

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

,

(29)

где R - универсальная газовая постоянная (см. таблицу 2 ГОСТ 30319.1);

- псевдокритическая молярная плотность природного газа (см. формулу (36));

- ацентрический фактор Питцера i-го компонента природного газа, значения {} для компонентов приведены в таблице А.5 (приложение А);

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

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

,

(30)

где и - соответственно вязкость i-го и j-го компонентов природного газа в разреженном состоянии;

и - соответственно молярная масса i-го и j-го компонентов природного газа, значения которых для каждого компонента приведены в таблицах А.1 и А.9 (приложение А).

Вязкость компонентов природного газа в разреженном состоянии () вычисляют по формуле

,

(31)

где {} - коэффициенты, значения которых для каждого компонента приведены в таблице А.6 приложения А;

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

4.3.3 Избыточную составляющую вязкости рассчитывают по формуле

,

(32)

где {}, {}, {} - коэффициенты и показатели степеней, значения которых приведены в таблице А.7 (приложение А);

,..., - параметры преобразований для приведенных значений плотности и температуры природного газа (см. формулу (33));

и - приведенные плотность и температура природного газа (см. формулы (34), (35)).

Параметры преобразований для приведенных значений плотности и температуры природного газа рассчитывают по формуле

, i = 1,2,...,6,

(33)

где {} и {} - коэффициенты, значения которых приведены в таблице А.8 (приложение А);

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

Приведенные плотность () и температуру () природного газа рассчитывают по формулам:

,

(34)

,

(35)

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

Псевдокритическую молярную плотность () и температуру () вычисляют по следующим формулам:

,

(36)

,

(37)

где {, }, {, } и {, } - критические плотности, молярные массы и критические температуры для компонентов (i, j) природного газа соответственно;

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

Значения критических параметров {}, {} и молярной массы {} для компонентов природного газа приведены в таблицах А.5 и А.1, А.9 (приложение А) соответственно.

5 Алгоритм расчета физических свойств природного газа

5.1 Исходные данные

5.1.1 Исходными данными для расчета физических свойств природного газа являются:

- молярные доли компонентов природного газа {};

- абсолютное давление природного газа;

- температура природного газа.

5.1.2 Молярные доли компонентов природного газа определяют хроматографическим анализом по ГОСТ 31371.1 - ГОСТ 31371.7. Измерения молярных долей компонентов могут выполняться как потоковыми, так и лабораторными хроматографами. Если измерены объемные доли компонентов природного газа, то для перевода их в молярные доли используют следующую формулу

, i = 1,..., Nc,

(38)

где - объемная доля i-го компонента природного газа;

- коэффициент сжимаемости i-го компонента природного газа при стандартных условиях, значения которого приведены в таблице А.1 приложения А;

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

5.1.3 Избыточное давление природного газа измеряют с применением соответствующих средств измерений. Для расчета абсолютного давления и перевода его в МПа применяют следующую формулу

,

(39)

где и - переводные коэффициенты, значения которых приведены в таблице 1;

- избыточное давление природного газа;

- атмосферное давление.

Таблица 1 - Переводные коэффициенты и

Единица измерения	Коэффициенты и
1
1
1 МПа	1
1 бар
1 мм рт. ст.

Пример - Перевод давления Р, МПа, при заданных = 10 ; = 750 мм рт. ст. По таблице 1 находим значения коэффициентов: ; , затем рассчитываем абсолютное давление: МПа.

5.1.4 Температуру природного газа измеряют с применением соответствующих средств измерений, как правило, в градусах Цельсия. Для перевода измеренной температуры t, °С в температуру Т, К применяют следующую формулу

T = t + 273,15.

(40)

5.2 Алгоритм расчета

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

- смесевой параметр размера () по формуле (5);

- давление нормировки () по формуле (4);

- молярную массу () по формуле (19);

- функции молярных долей компонентов и по формулам (7) - (17).

5.2.2 Расчет приведенной плотности () осуществляется в результате решения уравнения (1).

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

,

(41)

где R - универсальная газовая постоянная (см. таблицу 1 ГОСТ 30319.1).

Окончательное значение приведенной плотности () определяется по методу Ньютона в следующем итерационном процессе:

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

,

,

(42)

где безразмерные комплексы , рассчитывают по формулам (6) и (23) при плотности на итерационном шаге (к-1), т.е. при ;

б) условие завершения итерационного процесса

,

(43)

где приведенное давление рассчитывают по формуле

,

(44)

где безразмерный комплекс рассчитывают по формуле (6) при плотности на итерационном шаге (k), т.е. при .

Если условие (43) не выполняется, то продолжают итерационный процесс, возвращаясь к пункту а) итерационного процесса. Если условие (43) выполняется, то уравнение (1) считается решенным. После этого рассчитывают плотность по формуле (18) и коэффициент сжимаемости (z) по формуле (20) при , т.е. при найденном решении уравнения (1).

5.2.3 Расчет показателя адиабаты и скорости звука выполняют по формулам (21) и (22) при заданных () и () и найденном значении .

5.2.4 Расчет вязкости осуществляется по формулам (28) - (37) при заданных значениях (Т) и () и найденному значению молярной плотности

,

(45)

где - плотность, рассчитанная по формуле (18) при значении .

Блок-схема и примеры расчета физических свойств природного газа по представленным в стандарте методам приведены, соответственно, на рисунке 1 и в приложении Б.

6 Диапазон применения и погрешности расчета физических свойств природного газа

6.1 Диапазон применения и погрешности методов расчета физических свойств

6.1.1 Методы расчета, приведенные в настоящем стандарте, предназначены для расчета физических свойств природного газа в следующих диапазонах параметров:

- по температуре - от 250 до 350 К включительно;

- по давлению - от 0,1 до 30,0 МПа включительно.

При этом молярные доли компонентов природного газа не должны выходить за диапазоны, которые приведены в таблице 2.

6.1.2 Погрешности методов расчета физических свойств природного газа с диапазонами молярных долей компонентов, которые представлены в таблице 2, и во всем диапазоне температур и давлений, приведенном в 6.1.1, находятся в следующих пределах:

;

;

;

.

"Рисунок 1 - Блок-схема расчета физических свойств природного газа"

Таблица 2 - Компоненты природного газа и диапазоны молярных долей компонентов

Компонент	Диапазоны молярных долей
Метан
Этан
Пропан
Бутаны в сумме
Пентаны в сумме
Гексан
Азот
Диоксид углерода
Гелий
Водород
Примечания 1 Молярные доли остальных компонентов не превышают суммарно 0,0015. 2 Если в природном газе молярная доля гелия не превышает 0,0005, то при расчете физических свойств можно принять , а измеренную по ГОСТ 31371.1 - ГОСТ 31371.7 молярную долю гелия суммировать с молярной долей азота. 3 Если в природном газе молярная доля водорода не превышает 0,0005, то при расчете физических свойств можно принять , а измеренную по ГОСТ 31371.1 - ГОСТ 31371.7 молярную долю водорода суммировать с молярной долей азота. 4 Если измерены по ГОСТ 31371.1 - ГОСТ 31371.7 молярные доли кислорода и аргона, то при расчете физических свойств их значения необходимо суммировать с молярной долей азота. 5 Если измерены по ГОСТ 31371.1 - ГОСТ 31371.7 молярные доли н-гептана и н-октана, то при расчете физических свойств их значения необходимо суммировать с молярной долей н-гексана. 6 Для исключения возникновения дополнительной погрешности расчета физических свойств необходимо молярную массу смеси рассчитывать по формуле (19) с учетом всех компонентов, молярная доля которых не равна нулю (молярные массы кислорода, аргона, н-гептана и н-октана приведены в таблице А.9 приложения А).

Погрешности методов расчета физических свойств природного газа, соответствующие конкретным диапазонам температуры и давления, приведены в таблицах 3 - 5.

Таблица 3 - Погрешности методов расчета плотности и коэффициента сжимаемости (с доверительной вероятностью 95%)

Т, К	р, МПа	, , %
От 250,0 до 267,0 включ.	От 0,1 до включ.	0,1
	Св. до включ.	0,2
	Св. до 30,0 включ.	0,4
Св. 267,0 до 280,0 включ.	От 0,1 до включ.	0,1
	Св. до 30,0 включ.	0,2
Св. 280,0 до 295,0 включ.	От 0,1 до 30,0 включ.	0,1
Св. 295,0 до 310,0 включ.	От 0,1 до включ.	0,1
	Св. до 30,0 включ.	0,2
Св. 310,0 до 350,0 включ.	От 0,1 до включ.	0,1
	Св. до 30,0 включ.	0,2
Примечания 1 = 0,32353Т - 78,882. 2 = 0,94118Т - 221,29. З = 1,7308Т - 454,62. 4 = -1,2000Т + 384,00. 5 = 0,30000Т - 81,000.

Таблица 4 Погрешности методов расчета скорости звука и показателя адиабаты (с доверительной вероятностью 95 %)

Т, К	р, МПа	, %	, %
От 250,0 до 350,0 включ.	От 0,1 до включ.	0,2	0,5
	Св. до включ.	0,8	1,8
	Св. до 30,0 включ.	2,0	4,4
Примечания 1 = 0,06T - 9,0. 2 = 0,20T - 40,0.

Таблица 5 - Погрешности методов расчета коэффициента динамической вязкости (с доверительной вероятностью 95%)

Т, К	р, МПа	, %
От 250 до 350 включ.	От 0,1 до 1,0 включ.	0,6
	Св. 1,0 до 10,0 включ.	1,9
	Св. 10,0 до 20,0 включ.	2,6
	Св. 20,0 до 30,0 включ.	4,0

6.2 Учет погрешности измерения давления, температуры и молярных долей компонентов природного газа

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

,

(46)

,

(47)

,

(48)

,

(49)

,

(50)

где , , , , и - погрешности методов расчета соответственно коэффициента сжимаемости, плотности, скорости звука, показателя адиабаты и вязкости, значения которых приведены в таблицах 3 - 5;

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

6.2.2 Погрешности , , , , и вычисляют по следующим формулам:

,

(51)

,

(52)

,

(53)

,

(54)

,

(55)

где Nc - число компонентов природного газа;

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

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

, , , , - соответственно, коэффициент сжимаемости, плотность, скорость звука, показатель адиабаты и вязкость, алгоритм расчета которых приведен в 6.2.3;

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

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

В случае бинарной смеси формула (51) приобретает следующий вид:

,

(56)

где z - коэффициент сжимаемости, значение которого рассчитано при измеренных значениях , , и ;

- коэффициент сжимаемости, значения которого рассчитывают:

- для k = 1 при , , и ;

- для k = 2 при , , и ;

- для k = 3 при , , и ;

- для k = 4 при , , и ;

- коэффициент сжимаемости, значения которого рассчитывают:

- для k = 1 при , , и ;

- для k = 2 при , , и ;

- для k = 3 при , , и ;

- для k = 4 при , , и .

При этом значения давления, температуры и молярных долей компонентов с нижними индексами, включающими плюс и минус, рассчитывают по формулам:

,

(57)

,

(58)

,

(59)

,

(60)

,

(61)

,

(62)

,

(63)

,

(64)

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

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

Международный стандарт ISO 20765-1:2005(Е)*	Natural gas - Calculation of thermodynamic propertis - Part 1: Gas phase properties for transmission and distribution applications
ГОСТ Р 8.770-2011	Государственная система обеспечения единства измерений. Газ природный. Коэффициент динамической вязкости сжатого газа с известным компонентным составом. Метод расчетного определения

_____________________________

* С указанным стандартом можно ознакомиться в ФГУП "Стандартинформ".