Таблица N 1. Значения коэффициента "a"

Таблица N 1

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

3. Угловой коэффициент облучения является геометрической характеристикой, зависящей от взаимного расположения и формы поверхностей источника излучения и объекта. Угловой коэффициент облучения единичной площадки F₂ от факела с видимой поверхностью F₁ имеет следующий вид

где F₁ - излучающая поверхность факела, видимая со стороны облучаемой площадки;

- угол между нормалью к поверхности факела F₁ и направлением на облучаемую площадку F₂;

- угол между нормалью к облучаемой площадке F₂ и направлением на излучающую поверхность факела F₁;

r - расстояние между поверхностью факела F₂ и облучаемой площадкой F₂.

Приближенное значение интеграла (5), т.е. значения в конкретной точке территории с координатами (x, y), вычисляется численными методами при предварительном разбиении поверхности факела на элементарные площадки (рисунок 1) по формуле

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

Рисунок 1. Схема (пример) аппроксимации усеченного конуса

пламени факела усеченной пирамидой

При этом рекомендуется использовать аппроксимацию поверхности пламени следующими приближениями:

приближение лежащим полуцилиндром;

приближение прямым или наклонным цилиндром;

приближение усеченным конусом.

Рекомендуемые к использованию при расчете аналитические выражения коэффициента для указанных форм факела при расчетах q(x, y) приведены в пунктах 4, 5 настоящего приложения.

4. Расчет геометрических размеров пламени и теплового излучения при горении газа по сценариям из группы "Пожар в котловане"

Для сценариев указанной группы пламя пожара моделируется в виде цилиндрического твердого теплового излучателя, вертикального или наклонного (рисунок 2)

Рисунок 2. Схематическое представление пламени пожара

на газопроводе при сценариях из группы "Пожар в котловане"

в виде вертикального цилиндра

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

где Q_ф - общее тепловыделение пожара, кВт; - низшая теплота сгорания метана, кДж/кг; G, кг/с - суммарный массовый расход газа при его аварийном истечении из двух концов разрушенного газопровода на заданный момент времени (отсчет времени - от момента разрушения газопровода) или его осредненное значение за заданный промежуток времени .

Коэффициент излучения струевого пламени принимается . Значение удельной интенсивности излучения пламени E_f при сценариях группы "Пожар в котловане" не может превышать 120 кВт/м². Если расчеты по формуле (4) дают более высокие значения E_f, то для последующих оценок теплового излучения, воспринимаемого объектом-приемником, используют максимально возможное значение E_f = 120 кВт/м².

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

Рисунок 3. Модель излучателя в виде наклонного цилиндра

Формулы для угловых коэффициентов облучения от наклонного цилиндра в направлении ветра:

Формулы для излучения от наклонного цилиндра в направлении, перпендикулярном ветру:

где - угол отклонения оси пламени (цилиндра) от вертикали под действием ветра, рассчитываемый по формуле

где - средняя в пределах высоты пламени скорость ветра, м/с; m - массовая скорость поступления топлива в зону реакции с единицы поверхности очага горения, кг/м²·с; - плотность паров топлива при температуре поверхности раздела фаз (можно принимать равной плотности газа при нормальных условиях), кг/м³; g - ускорение силы тяжести, м/с²; D_эф - эффективный (видимый) диаметр очага горения, м.

Для других углов расположения площадок (не по направлению ветра, как в формулах 9а - б, и не перпендикулярно ветру, как в формулах 9в - г) значения и рассчитываются численным интегрированием по формуле (5).

Значения и при можно получить как частный случай общего решения.

Расчет величины для различных значений (с подветренной и наветренной стороны) и характерного отношения L_ф / D_эф = 3 показан на рисунке 4.

Рисунок 4. Угловой коэффициент облучения единичной

площадки на уровне поверхности земли от наклонного цилиндра

(R_эк = D_эф / 2)

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

Для вертикальных поверхностей стены (например, узла сопряжения с крышей)

Верхние индексы (+ и -) означают соответственно наклон пламени в сторону единичной площадки и от нее.

5. Расчет геометрических размеров пламени и теплового излучения при горении газа по сценариям из группы "Струевое пламя"

При расчете принимается:

- длина видимой части пламени настильной (прилегающей к земле) струи увеличивается на 20 - 25% по сравнению с длиной свободного пламени;

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

начальный ("слабосветящийся") участок факела длиной также оказывает незначительное влияние на характеристики теплового излучения;

для сценариев группы "Струевое пламя" факел моделируется как твердый тепловой излучатель в виде приподнятого или лежачего усеченного наклонного конуса (рисунок 6).

Рисунок 5. Встречное струевое горение. Излучающая

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

Рисунок 6. Схематическое представление струевого

пламени при сценариях из группы "Струевое пламя" (показана

одна струя)

Общая длина струи определяется по формуле (13) (с учетом заданной интенсивности истечения G₁ (или G₂))

где - диаметр аварийного газопровода;

Отрыв факела пламени от среза разорванной трубы (отверстия) h, м, определяется по формуле

Длина видимой части пламени R_L, м, определяется по формуле

При моделировании струи в виде усеченного конуса ширина малого основания W₁, м, определяется по формуле

где и - плотности газа соответственно в ударной плоскости и в критическом сечении, рассчитанные по приложению 9, кг/м³.

Ширина большого основания усеченного конуса W₂, м, определяется по формуле

При моделировании струи в виде цилиндра (на ранней стадии истечения) ширина основания W₀, м, определяется по формуле

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

Остальные геометрические размеры настильной струи рассчитываются по формулам (14) - (18).

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

При необходимости уточнения величины рекомендуется использовать следующую формулу , где u_стр - скорость достигаемая в струе после ее расширения до давления окружающей среды 1 атм. Величина u_стр рассчитывается в соответствии с приложением N 9 настоящего Руководства.

Значение удельной интенсивности излучения пламени E_f при таком горении не может превышать 200 кВт/м². Если расчеты по формуле (2) дают более высокие значения E_f, то для последующих оценок теплового излучения, воспринимаемого объектом-приемником, следует использовать максимально возможное значение E_f = 200 кВт/м².

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

Угловой коэффициент для одиночного настильного факела может быть рассчитан аналитически, если принять, что излучающая поверхность факела аппроксимируется боковой поверхностью полуцилиндра, лежащего на поверхности земли с характерной длиной и диаметром , где L_max = 1,25L_ф. При консервативных оценках длина полуцилиндра принимается равной 1,25L_ф. В этом случае форма пламени имеет вид лежащего на земле полуцилиндра, который имеет три излучающие поверхности: две концевых поверхности, через которые проходит ось цилиндра, S_C1 и S_C2, и боковую поверхность S₃ (рисунок 7). Вблизи места аварии размещаются два полуцилиндра, каждый из них аппроксимирует тепловое излучение от соответствующего пламени факела.

Рисунок 7. Модель излучателя в виде двух полуцилиндров

Расчет угловых коэффициентов излучения проводится для точек поверхности земли, расположенных на лучах Л₁ ... Л₆ от боковой поверхности пламени, на лучах Л_C1 и Л_C2 от торцевых поверхностей цилиндра. Лучи Л₁, Л₃, Л₄ и Л₆ проводятся по нормали к оси цилиндра через боковую поверхность цилиндра, лучи Л₂ и Л₅ - по нормали к оси цилиндра на половине его длины. Лучи Л_C1 и Л_C2 проводятся по оси цилиндра в направлении от торцевых поверхностей, начало этих лучей лежит на торцевых поверхностях цилиндра.

Для лучей Л₁ ... Л₆ угловой коэффициент излучения от лежащего полуцилиндра для вертикально ориентированных площадок-приемников излучения и для горизонтально расположенных площадок , расположенных на поверхности грунта, рассчитывается по формулам

где R_экв - эффективный радиус факела пламени, определяется по формуле - расстояние по горизонтали от оси цилиндра до приемника излучения, м; L_ф1 - длина полуцилиндра пламени с одной стороны от луча, проведенного по нормали к оси цилиндра, м; L_ф2 - длина полуцилиндра пламени с другой стороны от луча, проведенного по нормали к оси цилиндра, м.

Для лучей Л_C1 и Л_C2 расчет угловых коэффициентов излучения при x > R_'экв проводится по формуле (полуокружность для точек, лежащих на оси)

В силу симметрии течения расчеты проводят только для точек расположенных на лучах Л_1, Л₂ и Л_C1. Для луча Л₁ L_ф1 = 0 и L_ф2 = 0.8 · L_B. Для луча Л₂ L_ф1 = L_ф2 = 0.4 · L_B. Результаты, полученные для луча Л₁ переносят на лучи Л₃, Л₄ и Л₆; результаты для луча Л₂ - на луч Л₅; результаты для луча Л_C1 - на луч Л_C2.

На поверхности земли под нижней поверхностью полуцилиндра предполагается прямой огневой контакт с пламенем при тепловых нагрузках на поверхность объектов не менее 200 кВт/м².

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