ГОСТ 28656-2019. Межгосударственный стандарт. Газы углеводородные сжиженные. Расчетный метод определения плотности и давления насыщенных паров

(введен в действие Приказом Росстандарта от 03.04.2019 № 119-ст)

ИСТОРИЯ ИЗМЕНЕНИЙ НОРМАТИВНО-ТЕХНИЧЕСКИХ АКТОВ

flowchart LR
  n0["ГОСТ 28656-90"]:::rst
  click n0 "https://bukva.ai/docs/popravka-popravki-gost-32115-2013" "ГОСТ 28656-90. Газы углеводородные сжиженные. Расчетный метод определения плотности и давления насыщенных паров (утв. Постановлением Госстандарта СССР от 04.09.1990 N 2499)"
  n1["ГОСТ 28656-2019"]:::doc
  click n1 "https://bukva.ai/docs/standart-gost-28656-2019" "ГОСТ 28656-2019. Межгосударственный стандарт. Газы углеводородные сжиженные. Расчетный метод определения плотности и давления насыщенных паров (введен в действие Приказом Росстанда"
  n0 --> n1
classDef rst fill:#FFEBEB,stroke:#8B0000,color:#000;
classDef doc fill:#EBFFEB,stroke:#006400,color:#000;
classDef notEnt fill:#FFFFEB,stroke:#FF8C00,color:#000;

Введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 3 апреля 2019 г. № 119-ст

МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ

ГАЗЫ УГЛЕВОДОРОДНЫЕ СЖИЖЕННЫЕ

РАСЧЕТНЫЙ МЕТОД ОПРЕДЕЛЕНИЯ ПЛОТНОСТИ И ДАВЛЕНИЯ НАСЫЩЕННЫХ ПАРОВ

Liquefied hydrocarbon gases. Calculation method for determination of saturated vapour density and pressure

ГОСТ 28656-2019

МКС 75.160.30

Дата введения 1 января 2020 года

Предисловие

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

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

1 РАЗРАБОТАН Акционерным обществом «Волжский научно-исследовательский институт углеводородного сырья» (АО «ВНИИУС»)

2 ВНЕСЕН Межгосударственным техническим комитетом по стандартизации МТК 52 «Природный и сжиженные газы»

3 ПРИНЯТ Межгосударственным советом по стандартизации, метрологии и сертификации (протокол от 30 января 2019 г. № 115-П)

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

4 Приказом Федерального агентства по техническому регулированию и метрологии от 3 апреля 2019 г. № 119-ст межгосударственный стандарт ГОСТ 28656-2019 введен в действие в качестве национального стандарта Российской Федерации с 1 января 2020 г.

5 ВЗАМЕН ГОСТ 28656-90

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

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

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

1.2 Настоящий метод применяют для определения плотности СУГ в диапазоне температур от минус 50 °C до плюс 50 °C и избыточного давления насыщенных паров СУГ в интервале от 0,06 до 2,0 МПа при температурах минус 35 °C, минус 30 °C, минус 20 °C, плюс 45 °C.

1.3 Настоящий стандарт предназначен для вычисления плотности и давления насыщенных паров СУГ, в которых диапазон массовой доли компонентов составляет от 0,005% до 99,80%.

Примечания

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

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

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

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

ГОСТ 10679-2019 Газы углеводородные сжиженные. Метод определения углеводородного состава

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

ГОСТ 33012-2014 (ISO 7941:1988) Пропан и бутан товарные. Определение углеводородного состава методом газовой хроматографии

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

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

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

3.1.1 сжиженные углеводородные газы; СУГ: Смесь углеводородов (пропана, пропилена, бутанов, бутиленов и бутадиенов с присутствием метана, этана, этилена и (или) пентанов и пентенов), преобразованная в жидкое состояние.

3.1.2 плотность сжиженного углеводородного газа: Масса СУГ, заключенная в единице его объема при определенных значениях давления и температуры.

3.1.3 давление насыщенных паров: Давление, при котором жидкость находится в равновесном состоянии со своей газовой фазой; давление насыщенных паров складывается из избыточного давления и атмосферного давления.

3.1.4 абсолютное давление: Истинное давление, отсчитываемое от абсолютного нуля (давление абсолютного вакуума).

3.1.5 избыточное давление: Давление, равное разности между абсолютным и атмосферным давлением.

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

3.1.7 идеальный газ: Газ, подчиняющийся законам идеального газа.

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

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

4 Метод определения плотности сжиженных углеводородных газов

4.1 Определение плотности сжиженных углеводородных газов

4.1.1 Значение плотности СУГ \rho_{t} , кг/м3, вычисляют на основе закона аддитивности по данным измеренного компонентного состава, определенного хроматографическим методом и значениям плотности индивидуальных углеводородов, входящих в состав СУГ, при заданной температуре по формуле

\rho_{t} = 100/ \sum_{i = 1}^{n} \frac{w_{i}}{\rho_{it}} , (1)

где n — число компонентов сжиженного газа;

wi — массовая доля i-го компонента, %;

\rho_{it} — плотность i-го компонента при данной температуре t, кг/м3.

4.1.2 Если компонентный состав измерен в молярных долях, то плотность вычисляют по формуле

\rho_{t} = \sum_{i = 1}^{n} x_{i} \cdot \rho_{it} , (2)

где xi — молярная доля i-го компонента, доли единицы.

Компонентный состав определяют по ГОСТ 10679 или ГОСТ 33012.

4.1.3 Плотность индивидуальных углеводородов в жидком состоянии в зависимости от температуры приведена в таблице А.1 (приложение А).

Значения молярных масс индивидуальных компонентов приведены в приложении Б.

4.1.4 Если в таблице А.1 (приложение А) отсутствует значение плотности компонента при конкретной температуре измерений, то ее значение вычисляют интерполированием табличных значений плотностей, соответствующих температурам, ближайшим к данной.

4.1.5 Примеры расчета плотности СУГ приведены в таблицах В.1, В.2 (приложение В).

4.2 Оформление результатов вычисления плотности сжиженных углеводородных газов

4.2.1 За результат вычисления плотности СУГ при данной температуре принимают значение единичного определения.

4.2.2 Результат вычисления плотности СУГ записывают в виде

\rho_{t} \pm U_{(\rho_{t})} , (3)

где U_{(\rho_{t})} — расширенная неопределенность результата вычисления плотности для данной температуры t, кг/м3, при коэффициенте охвата k = 2, [1] - [3].

U_{(\rho_{t})} вычисляют по таблице 1. Вычисленные значения плотности СУГ и расширенной неопределенности (абсолютной погрешности) округляют до первого десятичного знака.

Таблица 1

Расширенная неопределенность результатов вычисления плотности сжиженных углеводородных газов

4.3 Требования к показателям точности метода

Метод обеспечивает получение результатов вычисления плотности СУГ по измеренному компонентному составу со значением расширенной неопределенности U_{(\rho_{t})} , не превышающей значений, приведенных в таблице 1, при доверительной вероятности P = 0,95.

5 Метод определения давления насыщенных паров

5.1 Давление насыщенных паров СУГ вычисляют по углеводородному составу, определенному методом газовой хроматографии в молярных долях, и значениям летучести углеводородов, входящих в состав СУГ, соответствующим заданной температуре измерений.

Углеводородный состав, определенный в массовых долях, пересчитывают в молярные доли xi по формуле

x_{i} = \frac{w_{i}}{M_{i}} / \sum_{i = 1}^{n} \frac{w_{i}}{M_{i}} , (4)

где Mi — молярная масса i-го компонента по таблице Б.1 приложения Б, кг/кмоль.

5.2 Абсолютное давление насыщенных паров СУГ P, МПа, вычисляют методом последовательного приближения, задавая произвольные значения двух ближайших значений давления насыщенных паров при данной температуре (приложение Г), по формуле

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} , (5)

где P_{z}\prime — меньшее выбранное значение абсолютного давления СУГ, МПа, по таблицам Г.1 - Г.₈ (приложение Г);

P_{z}\prime\prime — большее выбранное значение абсолютного давления СУГ, МПа, по таблицам Г.1 - Г.8 (приложение Г).

Пример расчета давления насыщенных паров методом последовательного приближения приведен в приложении Д.

Значения \Delta P_{z}\prime и \Delta P_{z}\prime\prime вычисляют по формулам:

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime ; (6)

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime , (7)

где P_{0}\prime и P_{0}\prime\prime — значения абсолютного давления насыщенных паров, МПа, вычисленные по формулам:

P_{0}\prime = \Sigma x_{i}f_{i}\prime ; (8)

P_{0}\prime\prime = \Sigma x_{i}f_{i}\prime\prime , (9)

где f_{i}\prime и f_{i}\prime\prime — значения летучести (фугитивности) i-го компонента СУГ при абсолютных давлениях P_{z}\prime и P_{z}\prime\prime , МПа, приведенные в таблицах Г.1 - Г.8 (приложение Г).

В результате вычисления должно соблюдаться условие P_{0}\prime>P_{z}\prime . Если P_{0}\prime \leq P_{z}\prime , то расчет прекращают, задают следующую пару значений давления насыщенных паров и повторяют процедуру приближения.

5.3 Избыточное давление насыщенных паров СУГ Pизб, МПа, вычисляют по формуле

Pизб = P — Pатм, (10)

где P — абсолютное давление насыщенных паров СУГ, МПа;

Pатм — атмосферное давление, МПа; Pатм = 101,3 кПа (0,1 МПа).

5.4 Примеры расчета давления насыщенных паров СУГ приведены в таблицах Е.1 - Е.4 (приложение Е).

5.5 Оформление результатов вычисления давления насыщенных паров СУГ

5.5.1 За результат вычисления значения давления насыщенных паров СУГ при данной температуре принимают значение единичного определения.

5.5.2 Результат вычисления давления насыщенных паров СУГ Pизб, МПа, записывают в виде

Pизб +/- U(Pизб), (11)

где U(Pизб) — расширенная неопределенность результата вычисления давления насыщенных паров для данной температуры t, МПа, при коэффициенте охвата k = 2;

U(Pизб) вычисляют по таблице 2.

Таблица 2

Расширенная неопределенность результатов вычисления давления насыщенных паров СУГ

Вычисленные значения давления насыщенных паров СУГ и расширенной неопределенности (абсолютной погрешности) округляют до второго десятичного знака.

5.6 Требования к показателям точности метода

Метод обеспечивает получение результатов вычисления избыточного давления насыщенных паров СУГ по измеренному компонентному составу со значением расширенной неопределенности U(Pизб), не превышающей значений, приведенных в таблице 2 при доверительной вероятности P = 0,95.

Приложение А (обязательное)

ЗНАЧЕНИЯ ПЛОТНОСТИ УГЛЕВОДОРОДОВ В ЖИДКОМ СОСТОЯНИИ

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

Таблица А.1

Значения плотности углеводородов в жидком состоянии в зависимости от температуры

Продолжение таблицы А.1

Продолжение таблицы А.1

Окончание таблицы А.1

Приложение Б (обязательное)

ЗНАЧЕНИЯ МОЛЯРНЫХ МАСС ИНДИВИДУАЛЬНЫХ КОМПОНЕНТОВ

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

Таблица Б.1

Значения молярных масс индивидуальных компонентов

Приложение В (рекомендуемое)

ПРИМЕР РАСЧЕТА ПЛОТНОСТИ

В таблицах В.1, В.2 приведены примеры расчета плотности при различных температурах.

Таблица В.1

Пример расчета плотности СУГ при температуре 20 °C через массовые доли

Таблица В.2

Пример расчета плотности СУГ при температуре 20 °C через молярные доли

Приложение Г (обязательное)

ЗНАЧЕНИЯ ЛЕТУЧЕСТИ (ФУГИТИВНОСТИ) КОМПОНЕНТОВ СЖИЖЕННЫХ УГЛЕВОДОРОДНЫХ ГАЗОВ

В таблицах Г.1 - Г.8 приведены значения летучести (фугитивности) компонентов СУГ при различных температурах.

Таблица Г.1

Значения летучести (фугитивности) компонентов СУГ при температуре плюс 45 °C

Таблица Г.2

Значения летучести (фугитивности) компонентов СУГ при температуре минус 20 °C

Таблица Г.3

Значения летучести (фугитивности) компонентов СУГ при температуре минус 30 °C

Таблица Г.4

Значения летучести (фугитивности) компонентов СУГ при температуре минус 35 °C

Таблица Г.5

Значения летучести (фугитивности) непредельных углеводородов СУГ при температуре плюс 45 °C

Таблица Г.6

Значения летучести (фугитивности) непредельных углеводородов СУГ при температуре минус 20 °C

Таблица Г.7

Значения летучести (фугитивности) непредельных углеводородов СУГ при температуре минус 30 °C

Таблица Г.8

Значения летучести (фугитивности) непредельных углеводородов СУГ при температуре минус 35 °C

Приложение Д (рекомендуемое)

ПРИМЕР РАСЧЕТА ДАВЛЕНИЯ НАСЫЩЕННЫХ ПАРОВ ПРИ ТЕМПЕРАТУРЕ ПЛЮС 45 °C МЕТОДОМ ПОСЛЕДОВАТЕЛЬНОГО ПРИБЛИЖЕНИЯ

Задают произвольные значения абсолютного давления насыщенных паров P_{z}\prime и P_{z}\prime\prime .

Принимают P_{z}\prime = 1,5\text{МПа} и P_{z}\prime\prime = 2,0\text{МПа} .

При выбранных значениях давления насыщенных паров согласно данным таблиц Г.1, Г.5 (приложение Г) выбирают значения летучести f_{i}\prime и рассчитывают P_{0}\prime по формуле P_{0}\prime = \Sigma x_{i}f_{i}\prime .

Расчет P_{0}\prime приведен в таблице Д.1.

Таблица Д.1

Расчет P_{0}\prime при P_{z}\prime = 1,5\text{МПа} и P_{z}\prime\prime = 2,0\text{МПа}

Так как при P_{z}\prime = 1,5\text{МПа} получают P_{0}\prime = 0,54\text{МПа} , следовательно, условие P_{0}\prime>P_{z}\prime не выполняется, расчет прерывают.

Задают следующую пару значений абсолютного давления насыщенных паров P_{z}\prime и P_{z}\prime\prime и повторяют процедуру.

Принимают пару значений P_{z}\prime = 1,0\text{МПа} и P_{z}\prime\prime = 1,5\text{МПа} и возобновляют расчет (см. таблицу Д.2).

Таблица Д.2

Расчет P_{0}\prime при P_{z}\prime = 1,0\text{МПа} и P_{z}\prime\prime = 1,5\text{МПа}

По результатам расчета из таблицы Д.2 получают P_{0}\prime = 0,51\text{МПа} менее P_{z}\prime = 1,0\text{МПа} , следовательно, условие P_{0}\prime>P_{z}\prime не выполняется, расчет прекращают, задают следующую пару значений абсолютного давления насыщенных паров P_{z}\prime и P_{z}\prime\prime и повторяют процедуру.

Принимают пару значений P_{z}\prime = 0,5\text{МПа} и P_{z}\prime\prime = 1,0\text{МПа} и возобновляют расчет (см. таблицу Д.3).

Таблица Д.3

Расчет P_{0}\prime при P_{z}\prime = 0,5\text{МПа} и P_{z}\prime\prime = 1,0\text{МПа}

По результатам расчета получают P_{0}\prime = 0,47\text{МПа} менее P_{z}\prime = 0,5\text{МПа} , следовательно, условие P_{0}\prime>P_{z}\prime не выполняется, расчет прекращают, задают следующую пару значений абсолютного давления насыщенных паров P_{z}\prime и P_{z}\prime\prime и повторяют процедуру.

Принимают пару значений P_{z}\prime = 0,1\text{МПа} и P_{z}\prime\prime = 0,5\text{МПа} .

При выбранных значениях P_{z}\prime и P_{z}\prime\prime рассчитывают P_{0}\prime и P_{0}\prime\prime (см. таблицу Д.4).

Таблица Д.4

Расчет P_{0}\prime и P_{0}\prime\prime при P_{z}\prime = 0,1\text{МПа} и P_{z}\prime\prime = 0,5\text{МПа}

В результате расчета при P_{z}\prime = 0,1\text{МПа} и P_{z}\prime\prime = 0,5\text{МПа} получают P_{0}\prime = 0,4254\text{МПа} более P_{z}\prime = 0,1\text{МПа} . Так условие P_{0}\prime>P_{z}\prime выполняется, то продолжают расчет.

При выбранных значениях давления насыщенных паров по таблицам Г.1, Г.5 (приложение Г) выбирают значения летучести f_{i}\prime\prime , рассчитывают P_{0}\prime\prime по формуле P_{0}\prime\prime = \Sigma x_{i}f_{i}\prime\prime

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime = 0,425 - 0,1 = 0,325 ;

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime = 0,470 - 0,5 = - 0,030 ;

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} = 0,1 + (0,5 - 0,1) \frac{0,325}{0,325 - (- 0,030)} = 0,47\text{МПа} .

В результате методом последовательного приближения получают P = 0,47 МПа, следовательно

Pизб = 0,47 — 0,1 = 0,37 МПа;

Pизб = (0,37 +/- 0,07) МПа.

Приложение Е (рекомендуемое)

ПРИМЕРЫ РАСЧЕТА ДАВЛЕНИЯ НАСЫЩЕННЫХ ПАРОВ

В таблицах Е.1 - Е.4 приведены примеры расчета давления насыщенных паров СУГ при различных температурах.

Таблица Е.1

Пример расчета давления насыщенных паров при температуре 45 °C

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime = 1,2646 - 1,0 = 0,2646 ;

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime = 1,3329 - 1,5 = - 0,1671 ;

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} = 1,0 + (1,5 - 1,0) \frac{0,2646}{0,2646 - (- 0,1671)} = 1,31\text{МПа} ;

Pизб = 1,31 — 0,1 = 1,21 МПа;

Pизб = (1,21 +/- 0,14) МПа.

Таблица Е.2

Пример расчета давления насыщенных паров при температуре минус 20 °C

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime = 0,2585 - 0,1 = 0,1585 ;

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime = 0,2677 - 0,5 = - 0,2323 ;

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} = 1,0 + (0,5 - 0,1) \frac{0,1585}{0,1585 - (- 0,2323)} = 0,262\text{МПа} ;

Pизб = 0,262 — 0,1 = 0,16 МПа;

Pизб = (0,16 +/- 0,04) МПа.

Таблица Е.3

Пример расчета давления насыщенных паров при температуре минус 30 °C

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime = 0,1976 - 0,1 = 0,0976 ;

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime = 0,2063 - 0,5 = - 0,2937 ;

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} = 0,1 + (0,5 - 0,1) \frac{0,0976}{0,0976 - (- 0,2937)} = 0,1998\text{МПа} ;

P_{\text{изб}} = 0,1998 - 0,1 = 0,0998 \approx 0,10\text{МПа} ;

Pизб = (0,10 +/- 0,02) МПа.

Таблица Е.4

Пример расчета давления насыщенных паров при температуре минус 35 °C

\Delta P_{z}\prime = P_{0}\prime - P_{z}\prime = 0,1765 - 0,1 = 0,0765 ;

\Delta P_{z}\prime\prime = P_{0}\prime\prime - P_{z}\prime\prime = 0,184 - 0,5 = - 0,3160 ;

P = P_{z}\prime + (P_{z}\prime\prime - P_{z}\prime) \frac{\Delta P_{z}\prime}{\Delta P_{z}\prime - \Delta P_{z}\prime\prime} = 0,1 + (0,5 - 0,1) \frac{0,0765}{0,0765 - (- 0,316)} = 0,178\text{МПа} .

P_{\text{изб}} = 0,178 - 0,1 = 0,078 \approx 0,08\text{МПа} ;

Pизб = (0,08 +/- 0,02) МПа.

БИБЛИОГРАФИЯ