"ГОСТ Р 70927-2023. Национальный стандарт Российской Федерации. Газ природный. Методы расчета физических свойств. Вычисление коэффициента сжимаемости в области низких температур"
(утв. и введен в действие Приказом Росстандарта от 21.09.2023 N 899-ст)

1 РАЗРАБОТАН Обществом с ограниченной ответственностью Центр Метрологии "СТП" (ООО ЦМ "СТП"), Обществом с ограниченной ответственностью "Газпром межрегионгаз" (ООО "Газпром межрегионгаз") - Управляющей организации АО "Газпром газораспределение"

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 052 "Природный и сжиженные газы"

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 21 сентября 2023 г. N 899-ст

4 ВВЕДЕН ВПЕРВЫЕ

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

1.2 Настоящий стандарт применяют для расчета коэффициента сжимаемости природного газа, находящегося в однофазном газовом состоянии при избыточном давлении до 0,6 МПа (включительно) и температуре от 220 до 250 К (исключительно), с содержанием компонентов согласно таблице 1.

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

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

ГОСТ 5542 Газ природный промышленного и коммунально-бытового назначения. Технические условия

ГОСТ 17310 Газы. Пикнометрический метод определения плотности

ГОСТ 20060 Газ природный. Определение температуры точки росы по воде

ГОСТ 20061 Газ природный. Определение температуры точки росы по углеводородам

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

ГОСТ 30319.2 Газ природный. Методы расчета физических свойств. Вычисление физических свойств на основе данных о плотности при стандартных условиях и содержании азота и диоксида углерода

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

ГОСТ 31371.3 (ИСО 6974-3:2000) Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 3. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов до C₈ с использованием двух насадочных колонок

ГОСТ 31371.4 (ИСО 6974-4:2000) Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 4. Определение азота, диоксида углерода и углеводородов C₁ - C₅ и C₆₊ в лаборатории и с помощью встроенной измерительной системы с использованием двух колонок

ГОСТ 31371.5 (ИСО 6974-5:2014) Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 5. Определение азота, диоксида углерода и углеводородов C₁ - C₅ и C₆₊ изотермическим методом

ГОСТ 31371.6 (ИСО 6974-6:2002) Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 6. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов C₁ - C₈ с использованием трех капиллярных колонок

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

ГОСТ 34721 Газ природный. Определение плотности пикнометрическим методом

ГОСТ Р 8.654 Государственная система обеспечения единства измерений. Требования к программному обеспечению средств измерений. Основные положения

ГОСТ Р 53762 Газы горючие природные. Определение температуры точки росы по углеводородам

ГОСТ Р 53763 Газы горючие природные. Определение температуры точки росы по воде

ГОСТ Р 56333 Газы горючие природные. Стандартные условия измерения и вычисления физико-химических свойств

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

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

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

H - теплота сгорания, МДж/м³;

M - молярная масса, кг/кмоль;

T - термодинамическая температура, К;

p - абсолютное давление, МПа;

- молярная плотность, кмоль/м³;

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

z - коэффициент сжимаемости;

- относительная погрешность, %;

- плотность, кг/м³.

Нижние индексы:

а - азот;

с - физическая величина при стандартных условиях по ГОСТ Р 56333 (p_с = 0,101325 МПа, T_с = 293,15 К);

д - диоксид углерода;

э - эквивалентный углеводород.

Метод расчета коэффициента сжимаемости z природного газа основан на уравнении состояния следующего вида:

где B_m - второй вириальный коэффициент уравнения состояния, м³/кмоль;

C_m - третий вириальный коэффициент уравнения состояния, м⁶/кмоль²;

- молярная плотность, кмоль/м³.

Вириальные коэффициенты уравнения состояния B_m, м³/кмоль, и C_m, м⁶/кмоль², вычисляют по формулам:

где x_э, x_а, x_д - молярные доли эквивалентного углеводорода, азота и диоксида углерода соответственно.

Молярную долю эквивалентного углеводорода x_э и входящие в формулы (2) и (3) коэффициенты вычисляют по формулам:

Значение теплоты сгорания эквивалентного углеводорода H_э, МДж/м³, вычисляют по формуле

Значение молярной массы эквивалентного углеводорода M_э, кг/кмоль, вычисляют по формуле

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

4.2.1 Используя измеренные значения температуры, плотности природного газа при стандартных условиях, а также значения молярных долей азота и диоксида углерода вычисляют значения вириальных коэффициентов уравнения состояния B_m и C_m по формулам (2) - (18).

Примечания

1 Плотность при стандартных условиях определяют по ГОСТ 17310, ГОСТ 34721 или вычисляют по ГОСТ 31369.

2 Молярные доли азота и диоксида углерода измеряют по ГОСТ 31371.3 - ГОСТ 31371.7.

3 Допускается принимать значения плотности при стандартных условиях, молярные доли азота и диоксида углерода из паспорта качества природного газа по ГОСТ 5542.

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

Примечание - Если подкоренное выражение получилось отрицательным, то его принимают равным нулю.

4.2.3 Примеры расчета коэффициента сжимаемости природного газа по настоящему методу приведены в приложении А.

Алгоритм пересчета объема природного газа, приведенного к стандартным условиям, с использованием настоящего метода приведен в приложении Б.

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

- температура от 220 до 250 К (исключительно) [от минус 53,15 °C до минус 23,15 °C (исключительно)];

- избыточное давление от 0 до 0,6 МПа;

- плотность при стандартных условиях от 0,66 до 1,05 кг/м³;

- молярная доля азота и диоксида углерода не более 0,20 соответственно;

- высшая теплота сгорания от 20 до 48 МДж/м³.

Примечание - Значение высшей теплоты сгорания H, МДж/м³, вычисляют по ГОСТ 31369 или по формуле

- однофазное газовое состояние.

5.1.2 Контроль однофазности природного газа выполняют путем определения температур точек росы по воде в соответствии с ГОСТ 20060 или ГОСТ Р 53763 и по углеводородам в соответствии с ГОСТ 20061 или ГОСТ Р 53762 на газораспределительных станциях и сопоставления их значений с величиной температуры природного газа в точке отбора пробы.

5.2.1 Программное обеспечение, реализующее настоящий метод расчета коэффициента сжимаемости, а также расчет фактического объемного расхода и объема (количества) природного газа, приведенных к стандартным условиям, при температурах от 220 до 250 К (исключительно) (далее - ПО), должно соответствовать требованиям ГОСТ Р 8.654.

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

5.2.3 В эксплуатационной документации на ПО должны быть указаны следующие идентификационные данные:

- идентификационное наименование;

- номер версии (идентификационный номер);

- цифровой идентификатор (контрольные суммы метрологически значимых частей);

- алгоритмы вычислений цифрового идентификатора.

5.2.4 При разделении ПО его метрологически значимая часть не должна поддаваться влиянию метрологически незначимой части согласно ГОСТ Р 8.654.

5.2.5 Относительная погрешность вычислений, обусловленная программной реализацией, не должна превышать +/- 0,02%.

5.3.1 Доверительные границы относительной погрешности (при доверительной вероятности 0,95) коэффициента сжимаемости , обусловленной методом расчета, изложенным в настоящем стандарте, не превышают +/- 0,11%.

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

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

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

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

z_qk+, z_qk- - значения коэффициента сжимаемости, алгоритм расчета которого приведен в 5.3.4.

5.3.4 Алгоритм расчета значений коэффициента сжимаемости z_qk+, z_qk-.

Значения коэффициента сжимаемости z_qk+ вычисляют:

- для k = 1 при p_и+, T_и, , x_а.и и x_д.и;

- для k = 2 при p_и, T_и+, , x_а.и и x_д.и;

- для k = 3 при p_и, T_и, , x_а.и и x_д.и;

- для k = 4 при p_и, T_и, , x_а.и+ и x_д.и;

- для k = 5 при p_и, T_и, , x_а.и и x_д.и+.

Значения коэффициента сжимаемости z_qk- вычисляют:

- для k = 1 при p_и-, T_и, , x_а.и и x_д.и;

- для k = 2 при p_и, T_и-, , x_а.и и x_д.и;

- для k = 3 при p_и, T_и, , x_а.и и x_д.и;

- для k = 4 при p_и, T_и, , x_а.и- и x_д.и;

- для k = 5 при p_и, T_и, , x_а.и и x_д.и-.

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

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

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

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

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

Объем природного газа, приведенный к стандартным условиям, V_с, м³, при применении расходомеров и счетчиков, измеряющих объемный расход (объем) газа в рабочих условиях, вычисляют по формуле на основании усредненных за час значений параметров расхода газа:

Значение коэффициента сжимаемости в рабочих условиях вычисляют по ГОСТ 30319.2 при температуре природного газа более или равной 250 К (минус 23,15 °C) и по настоящей методике при температуре менее 250 К (минус 23,15 °C).

При отсутствии в вычислителе (корректоре) реализованного алгоритма по настоящей методике и температуре природного газа менее 250 К (минус 23,15 °C) объем природного газа, приведенный к стандартным условиям, вычисляют по формуле на основании усредненных за час значений параметров расхода газа:

где V_{с_выч} - значение объема природного газа, приведенного к стандартным условиям, м³, рассчитанное вычислителем (корректором) по формуле

где z_выч - значение коэффициента сжимаемости, вычисленное по ГОСТ 30319.2, с использованием подстановочного значения температуры;

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