ГОСТ 30319.2-2015. Межгосударственный стандарт. Газ природный. Методы расчета физических свойств. Вычисление физических свойств на основе данных о плотности при стандартных условиях и содержании азота и диоксида углерода
(введен в действие Приказом Росстандарта от 10.11.2015 № 1744-ст)
Примечание
Документ введен в действие с 1 января 2017 года. Взамен ГОСТ 30319.1-96, ГОСТ 30319.2-96.
ИСТОРИЯ ИЗМЕНЕНИЙ НОРМАТИВНО-ТЕХНИЧЕСКИХ АКТОВ
flowchart LR
n0["ГОСТ 30319.1-96"]:::rst
click n0 "https://bukva.ai/docs/gost-30319-1-96" "ГОСТ 30319.1-96. Межгосударственный стандарт. Газ природный. Методы расчета физических свойств. Определение физических свойств природного газа, его компонентов и продуктов его... ("
n1["ГОСТ 30319.2-96"]:::rst
click n1 "https://bukva.ai/docs/gost-30319-2-96" "ГОСТ 30319.2-96. Межгосударственный стандарт. Газ природный. Методы расчета физических свойств. Определение коэффициента сжимаемости (введен в действие Постановлением Госстандарта "
n2["ГОСТ 30319.2-2015"]:::doc
click n2 "https://bukva.ai/docs/gost-30319-2-2015" "ГОСТ 30319.2-2015. Межгосударственный стандарт. Газ природный. Методы расчета физических свойств. Вычисление физических свойств на основе данных о плотности при стандартных условия"
n1 --> n2
n0 --> n2
classDef rst fill:#FFEBEB,stroke:#8B0000,color:#000;
classDef doc fill:#EBFFEB,stroke:#006400,color:#000;
classDef notEnt fill:#FFFFEB,stroke:#FF8C00,color:#000;Введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 10 ноября 2015 г. № 1744-ст
МЕТОДЫ РАСЧЕТА ФИЗИЧЕСКИХ СВОЙСТВ
ВЫЧИСЛЕНИЕ ФИЗИЧЕСКИХ СВОЙСТВ НА ОСНОВЕ ДАННЫХ О ПЛОТНОСТИ ПРИ СТАНДАРТНЫХ УСЛОВИЯХ И СОДЕРЖАНИИ АЗОТА И ДИОКСИДА УГЛЕРОДА
Natural gas. Methods of calculation of physical properties. Calculation of physical properties on base information on density of standards conditions and nitrogen and carbon dioxide contents
Дата введения 1 января 2017 года
Предисловие
Цели, основные принципы и порядок проведения работ по межгосударственной стандартизации установлены ГОСТ 1.0-92 «Межгосударственная система стандартизации. Основные положения» и ГОСТ 1.2-2009 «Межгосударственная система стандартизации. Стандарты межгосударственные, правила и рекомендации по межгосударственной стандартизации. Правила разработки, принятия, применения, обновления и отмены»
Сведения о стандарте
1 РАЗРАБОТАН Обществом с ограниченной ответственностью «Научно-исследовательский институт природных газов и газовых технологий — Газпром ВНИИГАЗ», Техническим комитетом по стандартизации ТК 52 «Природный и сжиженные газы»
2 ВНЕСЕН Федеральным агентством по техническому регулированию и метрологии Российской Федерации (Росстандарт)
3 ПРИНЯТ Межгосударственным советом по стандартизации, метрологии и сертификации по переписке (протокол от 27 августа 2015 г. № 79-П)
4 Приказом Федерального агентства по техническому регулированию и метрологии от 10 ноября 2015 г. № 1744-ст межгосударственный стандарт ГОСТ 30319.2-2015 введен в действие в качестве национального стандарта Российской Федерации с 1 января 2017 г.
5 ВЗАМЕН ГОСТ 30319.1-96, ГОСТ 30319.2-96
Информация об изменениях к настоящему стандарту публикуется в ежегодном информационном указателе «Национальные стандарты», а текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячном информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет
1. Область применения
1.1 Настоящий стандарт предназначен для расчета коэффициента сжимаемости, плотности, показателя адиабаты, коэффициента динамической вязкости природного газа и скорости распространения звука в среде природного газа по измеренным значениям давления, температуры, плотности при стандартных условиях, содержания азота и диоксида углерода.
1.2 Настоящий стандарт применяют для расчета указанных в 1.1 физических свойств природного газа при давлениях до 7,5 МПа включительно и температурах от 250 до 350 К.
1.3 Методы расчета физических свойств, приведенные в настоящем стандарте, могут быть использованы при разработке программного обеспечения вычислителей расхода природного газа.
2. Нормативные ссылки
В настоящем стандарте использованы ссылки на следующие межгосударственные стандарты:
ГОСТ 31369-2008 Газ природный. Вычисление теплоты сгорания, плотности, относительной плотности и числа Воббе на основе компонентного состава
ГОСТ 31371.1-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 1. Руководство по проведению анализа
ГОСТ 31371.2-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 2. Характеристики измерительной системы и статистические оценки данных
ГОСТ 31371.3-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 3. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов до C8 с использованием двух насадочных колонок
ГОСТ 31371.4-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 4. Определение азота, диоксида углерода и углеводородов C1 — C5 и C6+ в лаборатории и с помощью встроенной измерительной системы с использованием двух колонок
ГОСТ 31371.5-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 5. Определение азота, диоксида углерода и углеводородов C1 — C5 и C6+ в лаборатории и при непрерывном контроле с использованием трех колонок
ГОСТ 31371.6-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 6. Определение водорода, гелия, кислорода, азота, диоксида углерода и углеводородов C1 — C8 с использованием трех капиллярных колонок
ГОСТ 31371.7-2008 Газ природный. Определение состава методом газовой хроматографии с оценкой неопределенности. Часть 7. Методика выполнения измерений молярной доли компонентов
ГОСТ 30319.1-2015 Газ природный. Методы расчета физических свойств. Общие положения
ГОСТ 30319.3-2015 Газ природный. Методы расчета физических свойств. Вычисление физических свойств на основе данных о компонентном составе
Примечание - При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодному информационному указателю «Национальные стандарты», который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя «Национальные стандарты» за текущий год. Если ссылочный стандарт заменен (изменен), то при пользовании настоящим стандартом следует руководствоваться заменяющим (измененным) стандартом. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, применяется в части, не затрагивающей эту ссылку.
3. Термины, определения и обозначения
3.1 В настоящем стандарте применены термины и определения по ГОСТ 30319.1.
3.2 Основные условные обозначения величин, принятые в стандарте, приведены в таблице 2 ГОСТ 30319.1.
4. Определение коэффициента сжимаемости и плотности
4.1 Метод расчета
Метод расчета коэффициента сжимаемости и плотности основан на применении уравнения состояния, которое приведено в международном стандарте [1], и имеет следующий вид:
z = 1 + B_{m}\tilde{\rho} + C_{m}\tilde{\rho}^{2}, (1)
где Bm и Cm — коэффициенты уравнения состояния.
Коэффициенты уравнения состояния рассчитывают по следующим формулам:
где xэ — молярная доля эквивалентного углеводорода;
xу — молярная доля диоксида углерода.
Молярную долю эквивалентного углеводорода и входящие в формулы (2) и (3) коэффициенты рассчитывают по следующим формулам:
B2 = -0,1446 + 7,4091·10-4T — 9,1195·10-7T2, (6)
B23 = -0,339693 + 1,61176·10-3T — 2,04429·10-6T2, (7)
B3 = -0,86834 + 4,0376·10-3T — 5,1657·10-6T2, (8)
C2 = 7,8498·10-3 — 3,9895·10-5T + 6,1187·10-8T2, (10)
C3 = 2,0513·10-3 + 3,4888·10-5T — 8,3703·10-8T2, (11)
C223 = 5,52066·10-3 — 1,68609·10-5T + 1,57169·10-8T2, (12)
C233 = 3,58783·10-3 + 8,06674·10-6T — 3,25798·10-8T2, (13)
B* = 0,72 + 1,875·10-5(320 — T)2, (14)
C* = 0,92 + 0,0013·(T — 270). (15)
В формулах (5), (9) значение теплоты сгорания эквивалентного углеводорода (Hэ) рассчитывают по формуле
где Mэ — молярная масса эквивалентного углеводорода, значение которой рассчитывают по формуле
M_{\text{э}} = (24,05525 \cdot z_{c}\rho_{c} - 28,0135x_{a} - 44,01x_{y})/x_{\text{э}}, (17)
где xэ — молярная доля эквивалентного углеводорода, см. формулу (4);
xу — молярная доля диоксида углерода.
В формуле (17) коэффициент сжимаемости при стандартных условиях рассчитывают по формуле
z_{c} = 1 - (0,0741\rho_{c} - 0,006 - 0,063x_{a} - 0,0575x_{y})^{2}. (18)
4.2 Алгоритм расчета
4.2.1 Используя измеренные значения температуры, плотности при стандартных условиях, а также измеренные по ГОСТ 31371.1 — ГОСТ 31371.7 значения молярных долей азота и диоксида углерода, рассчитывают значения коэффициентов уравнения состояния (1) Bm и Cm по формулам (2) - (18).
Примечание — Допускается рассчитывать плотность при стандартных условиях по ГОСТ 31369.
4.2.2 Рассчитывают коэффициент сжимаемости при измеренных значениях температуры и давления по формуле
A_{2} = [A_{0} - (A_{0}^{2} - A_{1}^{3})^{0,5}]^{1/3}, (21)
4.2.3 Плотность природного газа рассчитывают по формуле
\rho = 10^{3}M_{m}p/(RTz), (26)
где Mm — молярная масса природного газа, которую рассчитывают по следующей формуле
M_{m} = 24,05525 \cdot z_{c}\rho_{c}. (27)
В формуле (27) коэффициент сжимаемости при стандартных условиях рассчитывают по формуле (18).
4.3 Диапазон применения и погрешности расчета
4.3.1 Диапазон применения и погрешности метода расчета коэффициента сжимаемости и плотности по уравнению состояния (1) приведены в таблице 1.
Диапазон применения и погрешности метода расчета коэффициента сжимаемости и плотности
| Диапазон применения и погрешности метода расчета | Отклонения от экспериментальных данных | ||||
| Диапазон применения | \rho_{c} , кг/м3 | p, МПа | \delta_{\text{м}} , % | \delta_{\text{сист}} , % | \delta_{\text{ст}} , % |
| 20 <= H, МДж/м3 <= 48 0,66 \leq \rho_{\text{с}} , кг/м3 <= 1,05 0 <= xа <= 0,20 0 <= xу <= 0,20 250 <= T, К <= 350 0,1 <= p, МПа <= 7,5 | От 0,66 до 0,70 включ. | От 0,1 до 3,0 включ. | 0,11 | 0,002 | 0,02 |
| Св. 3,0 до 7,5 включ. | 0,16 | 0,01 | 0,06 | ||
| Св. 0,70 до 0,75 включ. | От 0,1 до 3,0 включ. | 0,13 | -0,005 | 0,04 | |
| Св. 3,0 до 7,5 включ. | 0,17 | -0,005 | 0,07 | ||
| Св. 0,75 до 1,05 включ. | От 0,1 до 3,0 включ. | 0,12 | 0,006 | 0,03 | |
| Св. 3,0 до 7,5 включ. | 0,16 | -0,005 | 0,06 | ||
| Примечания 1 Оценка значений погрешностей метода расчета коэффициента сжимаемости и плотности выполнена в соответствии с А.1 (приложение А). 2 H — высшая теплота сгорания, рассчитанная по ГОСТ 31369 или по следующей формуле: H = 92,819(0,51447\rho_{c} + 0,05603 - 0,65689x_{a} - x_{y}) . 3 \delta_{\text{м}} — погрешность метода расчета коэффициента сжимаемости и плотности. 4 \delta_{\text{сист}} и \delta_{\text{ст}} — соответственно систематическое и стандартное отклонения расчетных значений коэффициента сжимаемости от экспериментальных данных. | |||||
4.3.2 Погрешности расчета коэффициента сжимаемости (\delta_{z}) и плотности (\delta_{\rho}) по уравнению состояния (1) с учетом погрешности измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода (исходных данных для расчета) вычисляют по формулам
\delta_{z} = (\delta_{\text{м}}^{2} + \delta_{z\text{ид}}^{2})^{0,5}, (28)
\delta_{\rho} = (\delta_{\text{м}}^{2} + \delta_{\rho\text{ид}}^{2})^{0,5}, (29)
где \delta_{\text{м}} — погрешность метода расчета коэффициента сжимаемости и плотности, значения которой приведены в таблице 1;
\delta_{z\text{ид}} — погрешность расчета коэффициента сжимаемости, которая появляется дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода;
\delta_{\rho\text{ид}} — погрешность расчета плотности, которая появляется дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода.
Формулы расчета \delta_{z\text{ид}} и \delta_{\rho\text{ид}} приведены в А.3 (приложение А).
5. Определение показателя адиабаты
5.1 Показатель адиабаты применяют при расчете коэффициента расширения природного газа.
5.2 Показатель адиабаты зависит от параметров состояния газа (давления и температуры), а в случае природного газа и от его компонентного состава.
5.3 Показатель адиабаты природного газа вычисляют по формуле Кобза [2], которая имеет следующий вид:
5.4 Диапазон применения и погрешности метода расчета показателя адиабаты по формуле (30) приведены в таблице 2.
5.5 Погрешность расчета показателя адиабаты (\delta_{k}) с учетом погрешности измерения давления, температуры, плотности при стандартных условиях и молярных долей азота (исходных данных для расчета) вычисляют по формуле
\delta_{k} = (\delta_{\text{м}}^{2} + \delta_{k\text{ид}}^{2})^{0,5}, (31)
где \delta_{\text{м}} — погрешность метода расчета показателя адиабаты, значения которой приведены в таблице 2;
\delta_{k\text{ид}} — погрешность расчета показателя адиабаты, которая появляется дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях и молярных долей азота.
Формулы расчета \delta_{k\text{ид}} приведены в А.3 (приложение А).
Диапазон применения и погрешности метода расчета показателя адиабаты
| Диапазон применения и погрешности метода расчета | Отклонения от данных ГОСТ 30319.3 | ||||
| Диапазон применения | T, К | p, МПа | \delta_{\text{м}} , % | \delta_{\text{сист}} , % | \delta_{\text{ст}} , % |
| 20 <= H, МДж/м3 <= 48 0,66 \leq \rho_{\text{с}} , кг/м3 <= 1,05 0 <= xа <= 0,20 0 <= xу <= 0,20 250 <= T, К <= 350 0,1 <= p, МПа <= 7,5 | От 250 до 260 включ. | От 0,1 до 6,0 включ. | 2,5 | -0,17 | 1,2 |
| Св. 6,0 до 7,5 включ. | 4,9 | -1,3 | 2,2 | ||
| Св. 260 до 350 включ. | От 0,1 до 7,5 включ. | 2,5 | -0,49 | 1,2 | |
| Примечания 1 Оценка значений погрешностей метода расчета показателя адиабаты выполнена в соответствии с А.2 (приложение А). 2 H — высшая теплота сгорания, рассчитанная по ГОСТ 31369 или по следующей формуле: H = 92,819(0,51447\rho_{c} + 0,05603 - 0,65689x_{a} - x_{y}) . 3 \delta_{\text{м}} — погрешность метода расчета показателя адиабаты. 4 \delta_{\text{сист}} и \delta_{\text{ст}} — соответственно систематическое и стандартное отклонения расчетных значений показателя адиабаты от данных, рассчитанных (полученных) по методу ГОСТ 30319.3. | |||||
6. Определение скорости звука
6.1 Скорость звука применяют при определении поправочного множителя показаний вибрационных плотномеров и при контроле точности результатов измерений расхода природного газа ультразвуковыми преобразователями расхода.
6.2 Скорость звука зависит от параметров состояния газа (давления и температуры), а в случае природного газа и от его компонентного состава.
6.3 Скорость звука вычисляют по формуле
u = 18,591(\frac{Tkz}{z_{c}\rho_{c}})^{0,5}. (32)
Входящие в правую часть формулы (32) свойства рассчитывают:
- показатель адиабаты — по формуле (30);
- коэффициент сжимаемости — по методу раздела 4;
- коэффициент сжимаемости при стандартных условиях — по формуле (18).
Формула (32) получена из уравнений термодинамики для скорости звука и показателя адиабаты [3].
6.4 Диапазон применения и погрешности метода расчета скорости звука по формуле (32) приведены в таблице 3.
Диапазон применения и погрешности метода расчета скорости звука
| Диапазон применения и погрешности метода расчета | Отклонения от данных ГОСТ 30319.3 | ||||
| Диапазон применения | T, К | p, МПа | \delta_{\text{м}} , % | \delta_{\text{сист}} , % | \delta_{\text{ст}} , % |
| 20 <= H, МДж/м3 <= 48 0,66 \leq \rho_{\text{с}} , кг/м3 <= 1,05 0 <= xа <= 0,20 0 <= xу <= 0,20 250 <= T, К <= 350 0,1 <= p, МПа <= 7,5 | От 250 до 260 включ. | От 0,1 до 6,0 включ. | 1,3 | -0,06 | 0,63 |
| Св. 6,0 до 7,5 включ. | 2,7 | -0,62 | 1,2 | ||
| Св. 260 до 350 включ. | От 0,1 до 7,5 включ. | 1,3 | -0,25 | 0,61 | |
| Примечания 1 Оценка значений погрешностей метода расчета скорости звука выполнена в соответствии с А.2 (приложение А). 2 H — высшая теплота сгорания, рассчитанная по ГОСТ 31369 или по следующей формуле: H = 92,819(0,51447\rho_{c} + 0,05603 - 0,65689x_{a} - x_{y}) . 3 \delta_{\text{м}} — погрешность метода расчета скорости звука. 4 \delta_{\text{сист}} и \delta_{\text{ст}} — соответственно систематическое и стандартное отклонения расчетных значений скорости звука от данных, рассчитанных (полученных) по методу ГОСТ 30319.3. | |||||
6.5 Погрешность расчета скорости звука (\delta_{u}) с учетом погрешности измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода (исходных данных для расчета) вычисляют по формуле
\delta_{u} = (\delta_{\text{м}}^{2} + \delta_{u\text{ид}}^{2})^{0,5}, (33)
где \delta_{\text{м}} — погрешность метода расчета скорости звука, значения которой приведены в таблице 3;
\delta_{u\text{ид}} — погрешность расчета скорости звука, которая появляется дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода.
Формулы расчета \delta_{u\text{ид}} приведены в А.3 (приложение А).
7. Определение коэффициента динамической вязкости
7.1 Вязкость применяют для вычисления числа Рейнольдса, которое является одной из важнейших характеристик течения вязкой среды и определяется отношением инерционных сил к силам вязкости. Число Рейнольдса применяют для определения коэффициента истечения.
7.2 Вязкость природного газа рассчитывают по формуле
\mu = \frac{3,24(T^{0,5} + 1,37 - 9,09\rho_{c}^{0,125})}{\rho_{c}^{0,5} + 2,08 - 1,5(x_{a} + x_{y})} [1 + \frac{p_{\text{п}}^{2}}{30(\text{Т}_{\text{п}} - 1)}], (34)
xу — молярная доля диоксида углерода;
pп и Tп — приведенные давление и температура, соответственно.
Приведенные давление и температуру вычисляют по формулам
где pпк и Tпк — псевдокритические давление в МПа и температура в К, которые рассчитывают по формулам
\text{р}_{\text{пк}} = 2,9585(1,608 - 0,05994\rho_{\text{с}} + \text{х}_{\text{у}} - 0,392\text{х}_{\text{а}}), (37)
\text{Т}_{\text{пк}} = 88,25(0,9915 + 1,759\rho_{\text{с}} - \text{х}_{\text{у}} - 1,681\text{х}_{\text{а}}). (38)
7.3 Диапазон применения и погрешности метода расчета вязкости по формуле (34) приведены в таблице 4.
Диапазон применения и погрешности метода расчета вязкости
| Диапазон применения и погрешности метода расчета | Отклонения от данных ГОСТ 30319.3 | |||
| Диапазон применения | p, МПа | \delta_{\text{м}} , % | \delta_{\text{сист}} , % | \delta_{\text{ст}} , % |
| 20 <= H, МДж/м3 <= 48 0,66 \leq \rho_{\text{с}} , кг/м3 <= 1,05 0 <= xа <= 0,20 0 <= xу <= 0,20 250 <= T, К <= 350 0,1 <= p, МПа <= 7,5 | От 0,1 до 1,0 включ. | 5,6 | -3,1 | 2,3 |
| Св. 1,0 до 7,5 включ. | 6,7 | -4,2 | 2,2 | |
| Примечания 1 Оценка значений погрешностей метода расчета вязкости выполнена в соответствии с А.2 (приложение А). 2 H — высшая теплота сгорания, рассчитанная по ГОСТ 31369 или по следующей формуле: H = 92,819(0,51447\rho_{c} + 0,05603 - 0,65689x_{a} - x_{y}) . 3 \delta_{\text{м}} — погрешность метода расчета вязкости. 4 \delta_{\text{сист}} и \delta_{\text{ст}} — соответственно систематическое и стандартное отклонения расчетных значений вязкости от данных, рассчитанных (полученных) по методу ГОСТ 30319.3. | ||||
7.4 Погрешность расчета вязкости (\delta_{\mu}) с учетом погрешности измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода (исходных данных для расчета) вычисляют по формуле
\delta_{\mu} = (\delta_{\text{м}}^{2} + \delta_{\mu\text{ид}}^{2})^{0,5}, (39)
где \delta_{\text{м}} — погрешность метода расчета вязкости, значения которой приведены в таблице 4;
\delta_{\mu\text{ид}} — погрешность расчета вязкости, которая появляется дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода.
Формулы расчета \delta_{\mu\text{ид}} приведены в А.3 (приложение А).
Примеры расчета физических свойств природного газа по представленным в стандарте методам приведены в приложении Б.
ОЦЕНКА ПОГРЕШНОСТЕЙ МЕТОДОВ РАСЧЕТА ФИЗИЧЕСКИХ СВОЙСТВ ПРИРОДНОГО ГАЗА
А.1 Погрешность метода расчета коэффициента сжимаемости и плотности
Погрешность метода расчета коэффициента сжимаемости вычисляют по формуле
\delta_{\text{м}} = [\delta_{\text{сист}}^{2} + (2\delta_{\text{ст}})^{2} + \delta_{\text{эксп}}^{2}]^{0,5}, (А.1)
где \delta_{\text{сист}} — систематическое отклонение расчетных значений коэффициента сжимаемости от экспериментальных данных;
\delta_{\text{ст}} — стандартное отклонение расчетных значений коэффициента сжимаемости от экспериментальных данных;
\delta_{\text{эксп}} — погрешность экспериментальных данных (0,1%).
Систематическое и стандартное отклонения рассчитывают по следующим формулам:
\delta_{\text{сист}} = \frac{1}{N} \sum_{i = 1}^{N} \delta_{i}, (А.2)
\delta_{\text{ст}} = [\frac{1}{N - 1} \sum_{i = 1}^{N} (\delta_{i} - \delta_{\text{сист}})^{2}]^{0,5}, (А.3)
где N — число экспериментальных данных;
\delta_{i} — отклонение расчетного значения коэффициента сжимаемости от его экспериментального значения в i-й точке экспериментальных данных, которое рассчитывают по формуле
\delta_{i} = 100[\frac{z_{\text{расч},i} - z_{\text{эксп},i}}{z_{\text{эксп},i}}] (i = 1, 2, ..., N — 1, N), (А.4)
где zрасч,i и zэксп,i — соответственно, расчетное и экспериментальное значения коэффициента сжимаемости в i-й точке экспериментальных данных.
Примечание - Численно значения погрешности метода расчета коэффициента сжимаемости и плотности равны.
А.2 Погрешности методов расчета показателя адиабаты, скорости звука и коэффициента динамической вязкости
Значения погрешностей рассчитаны по формуле
\delta_{\text{м}} = [\delta_{\text{сист}}^{2} + (2\delta_{\text{ст}})^{2} + \delta_{\text{ГОСТ}}^{2}]^{0,5}, (А.5)
где \delta_{\text{сист}} — систематическое отклонение расчетных значений показателя адиабаты, скорости звука и вязкости от данных, рассчитанных (полученных) по методу ГОСТ 30319.3;
\delta_{\text{ст}} — стандартное отклонение расчетных значений показателя адиабаты, скорости звука и вязкости от данных, полученных по методу ГОСТ 30319.3;
\delta_{\text{ГОСТ}} — погрешность данных, полученных по методу ГОСТ 30319.3, о показателе адиабаты, скорости звука и вязкости.
Систематическое и стандартное отклонения рассчитаны по формулам (А.2) и (А.3), соответственно, за исключением вычисления отклонения расчетного значения показателя адиабаты, скорости звука и вязкости от их значений в i-ой точке полученных по методу ГОСТ 30319.3 данных (\delta_{i}) , которые рассчитывают по следующим формулам:
\delta_{i} = 100[\frac{k_{\text{расч},i} - k_{\text{ГОСТ},i}}{k_{\text{ГОСТ},i}}] (i = 1, 2, ..., N — 1, N), (А.6)
\delta_{i} = 100[\frac{u_{\text{расч},i} - u_{\text{ГОСТ},i}}{u_{\text{ГОСТ},i}}] (i = 1, 2, ..., N — 1, N), (А.7)
\delta_{i} = 100[\frac{\mu_{\text{расч},i} - \mu_{\text{ГОСТ},i}}{\mu_{\text{ГОСТ},i}}] (i = 1, 2, ..., N — 1, N), (А.8)
где kрасч,i и kГОСТ,i — соответственно, расчетное значение показателя адиабаты и его значение, полученное по методу ГОСТ 30319.3;
uрасч,i и uГОСТ,i — соответственно, расчетное значение скорости звука и ее значение, полученное по методу ГОСТ 30319.3;
\mu_{\text{расч},i} и \mu_{\text{ГОСТ},i} — соответственно, расчетное значение вязкости и ее значение, полученное по методу ГОСТ 30319.3;
N — число данных, полученных по методу ГОСТ 30319.3.
А.3 Расчет дополнительных погрешностей методов расчета коэффициента сжимаемости, плотности, показателя адиабаты, скорости звука и коэффициента динамической вязкости
А.3.1 Погрешности расчета коэффициента сжимаемости (\delta_{z\text{ид}}) , плотности (\delta_{\rho\text{ид}}) , показателя адиабаты (\delta_{k\text{ид}}) , скорости звука (\delta_{u\text{ид}}) и вязкости (\delta_{\mu\text{ид}}) , которые появляются дополнительно в связи с погрешностью измерения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода, вычисляют по следующим формулам:
\delta_{z\text{ид}} = \frac{100}{z} [\sum_{k = 1}^{5} (z_{q_{k} +} - z_{q_{k} -})^{2}]^{0,5}, (А.9)
\delta_{\rho\text{ид}} = \frac{100}{\rho} [\sum_{k = 1}^{5} (\rho_{q_{k} +} - \rho_{q_{k} -})^{2}]^{0,5}, (А.10)
\delta_{k\text{ид}} = \frac{100}{k} [\sum_{k = 1}^{5} (k_{q_{k} +} - k_{q_{k} -})^{2}]^{0,5}, (А.11)
\delta_{u\text{ид}} = \frac{100}{u} [\sum_{k = 1}^{5} (u_{q_{k} +} - u_{q_{k} -})^{2}]^{0,5}, (А.12)
\delta_{\mu\text{ид}} = \frac{100}{\mu} [\sum_{k = 1}^{5} (\mu_{q_{k} +} - \mu_{q_{k} -})^{2}]^{0,5}, (А.13)
где qk — условное обозначение k-го параметра применяемых для расчета исходных данных, т.е. измеренные значения давления (pи), температуры (Tи), плотности при стандартных условиях (\rho_{\text{си}}) , молярных долей азота (xаи) и диоксида углерода (xуи) природного газа;
z, \rho , k, u и \mu — соответственно, коэффициент сжимаемости, плотность, показатель адиабаты, скорость звука и вязкость, значения которых рассчитывают при измеренных значениях давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода природного газа;
z_{q_{k} +} , \rho_{q_{k} +} , k_{q_{k} +} , u_{q_{k} +} и \mu_{q_{k} +} — соответственно, коэффициент сжимаемости, плотность, показатель адиабаты, скорость звука и вязкость, алгоритм расчета которых приведен в А.3.2;
z_{q_{k} -} , \rho_{q_{k} -} , k_{q_{k} -} , u_{q_{k} -} и \mu_{q_{k} -} — соответственно, коэффициент сжимаемости, плотность, показатель адиабаты, скорость звука и вязкость, алгоритм расчета которых приведен в А.3.2.
А.3.2 Ниже приведен алгоритм расчета значений коэффициента сжимаемости ( z_{q_{k} +} и z_{q_{k} -} ). Расчет аналогичных значений плотности ( \rho_{q_{k} +} и \rho_{q_{k} -} ), показателя адиабаты ( k_{q_{k} +} и k_{q_{k} -} ), скорости звука ( u_{q_{k} +} и u_{q_{k} -} ) и вязкости ( \mu_{q_{k} +} и \mu_{q_{k} -} ) осуществляют так же, как и для коэффициента сжимаемости z_{q_{k} +} и z_{q_{k} -} .
Значения коэффициента сжимаемости z_{q_{k} +} рассчитывают:
- для k = 1 при pи+, Tи, \rho_{\text{си}} , xаи и xуи,
- для k = 2 при pи, Tи+, \rho_{\text{си}} , xаи и xуи,
- для k = 3 при pи, Tи, \rho_{\text{си} +} , xаи и xуи,
- для k = 4 при pи, Tи, \rho_{\text{си}} , xаи+ и xуи,
- для k = 5 при pи, Tи, \rho_{\text{си}} , xаи и xуи+.
Значения коэффициента сжимаемости z_{q_{k} -} рассчитывают:
- для k = 1 при pи-, Tи, \rho_{\text{си}} , xаи и xуи,
- для k = 2 при pи, Tи-, \rho_{\text{си}} , xаи и xуи,
- для k = 3 при pи, Tи, \rho_{\text{си} -} , xаи и xуи,
- для k = 4 при pи, Tи, \rho_{\text{си}} , xаи- и xуи,
- для k = 5 при pи, Tи, \rho_{\text{си}} , xаи и xуи-.
При этом значения давления, температуры, плотности при стандартных условиях, молярных долей азота и диоксида углерода с нижними индексами, включающими плюс и минус, рассчитывают по формулам:
\text{р}_{\text{и} +} = \text{р}_{\text{и}}(1 + 0,005\delta_{\text{р}}), (А.14)
\text{р}_{\text{и} -} = \text{р}_{\text{и}}(1 - 0,005\delta_{\text{р}}), (А.15)
\text{Т}_{\text{и} +} = \text{Т}_{\text{и}}(1 + 0,005\delta_{\text{Т}}), (А.16)
\text{Т}_{\text{и} -} = \text{Т}_{\text{и}}(1 - 0,005\delta_{\text{Т}}), (А.17)
\rho_{\text{си} +} = \rho_{\text{си}}(1 + 0,005\delta_{\rho\text{с}}), (А.18)
\rho_{\text{си} -} = \rho_{\text{си}}(1 - 0,005\delta_{\rho\text{с}}), (А.19)
\text{х}_{\text{аи} +} = \text{х}_{\text{аи}}(1 + 0,005\delta_{\text{ха}}), (А.20)
\text{х}_{\text{аи} -} = \text{х}_{\text{аи}}(1 - 0,005\delta_{\text{ха}}), (А.21)
\text{х}_{\text{уи} +} = \text{х}_{\text{уи}}(1 + 0,005\delta_{\text{ху}}), (А.22)
\text{х}_{\text{уи} -} = \text{х}_{\text{уи}}(1 - 0,005\delta_{\text{ху}}), (А.23)
где \delta_{p} , \delta_{T} , \delta_{\rho c} , \delta_{xa} и \delta_{xy} — соответственно, погрешности измерения pи, Tи, \rho_{\text{си}} , xаи и xуи, численные значения которых определяют в соответствии с применяемыми методиками или средствами их измерений.
ПРИМЕРЫ РАСЧЕТА ФИЗИЧЕСКИХ СВОЙСТВ ПРИРОДНОГО ГАЗА
Б.1 Примеры расчета рекомендуется использовать в качестве тестовых данных при программной реализации методов расчета физических свойств природного газа, которые даны в настоящем стандарте.
Б.2 Примеры расчета приведены в форме таблиц. При этом в таблице Б.1 даны исходные данные для расчета физических свойств смесей, имитирующих природный газ, а в таблицах Б.2 и Б.3 приведены расчетные значения физических свойств этих смесей при соответствующих температурах и давлениях.
Исходные данные для расчета физических свойств смесей, имитирующих природный газ
| Исходные данные | Смеси | |
| № 1 | № 2 | |
| Плотность при стандартных условиях, кг/м3 | 0,7000 | 0,8263 |
| Молярная доля азота | 0,003 | 0,057 |
| Молярная доля диоксида углерода | 0,006 | 0,076 |
Расчетные значения физических свойств для смеси № 1
Расчетные значения физических свойств для смеси № 2
БИБЛИОГРАФИЯ
| [1] | Международный стандарт ISO 12213-3:2006 (E) <*> | Natural gas — Calculation of compression factor — Part 3: Calculation using physical properties |
| [2] | Кобза З., Добровольски Б., Гонтарек Я. (Польская высшая инженерная школа) Анализ влияния неточности определения показателя адиабаты природных газов на погрешность расчета расхода | |
| [3] | Шпильрайн Э.Э., Кессельман П.М. Основы теории теплофизических свойств веществ. М. «Энергия», 1977. 248 с | |
<*> С указанным стандартом можно ознакомиться в ФГУП «Стандартинформ».
(2)
(3)
(5)
(9)
(30)