Приложение 6

Приложение № 6
к пункту 12 Методики

Порядок проведения расчета и математические модели для определения времени блокирования путей эвакуации опасными факторами пожара

I. Порядок проведения расчета

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

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

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

Скорость выгорания для этих случаев определяется формулами:

            _

           | Ψуд × π × v2 × t2  - для кругового распространения пожара

           | Ψуд × 2 × v × t × b  - для линейного распространения пожара   

Ψ = <                          t                                                                                             , (П6.1)

           | Ψуд × F × √ -----  - для неустановившегося горения ГЖ

           |_                    tст

где Ψуд - удельная скорость выгорания (для жидкостей установившаяся), кг/(с × м2);
v - скорость распространения пламени, м/с;
b - ширина полосы горючей нагрузки, м;
tст - время стабилизации горения горючей жидкости, с;
F - площадь очага пожара, м2.

 

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

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

Предельно допустимые значения по каждому из опасных факторов пожара составляют:
по повышенной температуре - 70оС;
по тепловому потоку - 1400 Вт/м2;
по потере видимости - 20 м;
по пониженному содержанию кислорода - 0,226 кг/м3;
по каждому из токсичных газообразных продуктов горения
(СО2 - 0,11 кг/м3; СО - 1,16 × 10-3 кг/м3; HCL - 23 × 10-6 кг/м3).

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

Определяется время блокирования tбл:

tбл = min {tп.в.кр, tTкр, tт.г.кр, tO2кр, tт.п.кр}.                                                    (П6.2)

 

II. Классификация и область применения методов
математического моделирования пожара

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

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

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

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

для помещений, в которых один из геометрических размеров гораздо больше (меньше) остальных (тоннели, закрытые автостоянки большой площади и.т.д.);

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

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

 

III. Интегральная математическая модель расчета газообмена в здании, при пожаре

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

Уравнения движения, связывающие значения перепадов давлений на проемах с расходами газов через проемы, имеют вид:

                                           ________ 

Gji = sign (ΔPji) × µ × F ×√2×ρ×|ΔPji|  ,                                                                   (П6.3)

где Gji - расход газов через проем между двумя (j-м и i-м) смежными помещениями, кг/с;
µ - коэффициент расхода проема (µ = 0,8 для закрытых проемов и µ = 0,64 для открытых);
F - площадь сечения проема, м2;
ρ - плотность газов, проходящих через проем, кг/м3;
ΔPji - средний перепад полных давлений между j-м и i-м помещением, Па.

 

Направление (знак) расхода определяется знаком разности давлений ΔPji. В зависимости от этого плотность  ρ принимает различные значения.
Знак расхода газов (входящий в помещение расход считается положительным, выходящий - отрицательным) и значение ρ зависят от знака перепада давлений:

                           _ 

                          | -1, ρ = ρj, при ΔP < 0

sign (ΔP), ρ = <                                          .                                                               (П6.4)

                          | +1, ρ = ρi, при ΔP > 0

                           _

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

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

            Δp     

G = √ -----  ,                                                                                                                (П6.5)

            S

                           1 

где S = -------------------------  - характеристика гидравлического сопротивления на границе зон;

              2 × g × ρ × k × F2

F - площадь поперечного сечения шахты;
k - коэффициент (допускается принимать равным 0,05 с2/м);
g = 9,81 м/с2 - ускорение свободного падения;
Δp - перепад давлений между узлами.

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

Уравнение баланса массы выражается формулой:

d (ρj × Vj) / dt = Ψ + ΣkGk - ΣiGi ,                                                                               (П6.6)

где Vi - объем помещения, м3;
t - время, с;
ΣkGk - сумма расходов, входящих в помещение, кг/с;

ΣiGi - сумма расходов, выходящих из помещения, кг/с;

Ψ - скорость выгорания пожарной нагрузки, кг/c.

 

Уравнение сохранения энергии выражается формулой:

d(Cv × ρj × Vj × Tj) / dt = Cp × Σk (Tk  × Gk) - Cp × Tj × ΣiGi + Qг - Qw,                    (П6.7) 

где Сv, Сp - удельная изохорная и изобарная теплоемкости, кДж/(кг×K);
Ti, Tj - температуры газов в i- м и j-м помещениях, K;
QГ - количество тепла, выделяемого в помещении при горении, кВт;
Qw - тепловой поток, поглощаемый конструкциями и излучаемый через проемы, кВт.

Для помещения очага пожара величина QГ определяется по формуле:

Qг = η × Qн + I ,

где η - коэффициент полноты горения;
Qн - низшая теплота сгорания, кДж/кг;
I - энтальпия газифицированной горючей нагрузки.
Для остальных помещений QГ = 0.

Коэффициент полноты горения η определяется по формуле:

                                      Xox,a × Ga

η = ηa × K + (1 - K) × ----------------  ,                                                                           (П6.8)

                                       Ψ × Lox

где ηa - коэффициент полноты горения в режиме пожара, регулируемом горючей нагрузкой, определяемый формулой:

ηa = 0,63 + 0,2 × Xox,a + 1500 × X6ox,a .                                                                  (П6.9)

 

Коэффициент K рассчитывается по формуле:

        _           _            _       _                 _

      |   Xox,m  |B         |         |        Xox,m  ||

K = |  --------- |   exp  | B ×  | 1 - ---------  ||  ,                                                                (П6.10)

      |   Xox,a   |           |         |         Xox,a  ||

       _            _             _        _                 _

              _                          _

             |          Xox,a           |2

где B = | ---------------------- |   ;

             |      Xox,a - 0,01     |

             _                           _

Xox,a - начальная концентрация кислорода в помещении очага пожара;
Xox,m - текущая концентрация кислорода в помещении очага пожара.

 

Уравнение Менделеева-Клайперона выражается формулой:

                        R 

Pj = ρj × Tj × -----  ,                                                                                                      (П6.11)

                        M

где Pj - давление газа в j-м помещении, Па;
Tj - температура газа в j-м помещении, K;
R = 8,31 - универсальная газовая постоянная, Дж/(моль×К);
M - молярная масса газа, моль.

 

Параметры газа в помещении определяются из уравнения баланса масс отдельных компонентов продуктов горения и кислорода и уравнения баланса оптической плотности дыма.
Уравнение баланса масс отдельных компонентов продуктов горения и кислорода:

d(XLj × ρj × Vj) / dt = Ψ × LL + Σk(XL,k × Gk) - XL,j × ΣiGi  ,                                      (П6.12)

где XL,i, XL,j - концентрация L-го компонента продуктов горения в i- м и j-м помещениях, кг/кг;
LL - количество L-го компонента продуктов горения (кислорода), выделяющегося (поглощающегося) при сгорании одного килограмма пожарной нагрузки, кг/кг.

 

Уравнение баланса оптической плотности дыма:

Vj × dμj / dt = Ψ × Dm + Σkk × Gk) / ρk - μj × ΣiGi /  ρj  ,                                       (П6.13)

где μi, μj - оптическая плотность дыма в i- м и j-м помещениях, Нп×м-1;
Dm - дымообразующая способность пожарной нагрузки, Нп×м2/кг.

Оптическая плотность дыма при обычных условиях связана с расстоянием предельной видимости в дыму формулой:

lпр = 2,38/μ .                                                                                                                  (П6.14)

 

Для помещений без источника тепла система уравнений (П6.6), (П6.7) и (П6.8) упрощается и представляется в виде:

    _                                                    _______________

   |Gij (τ) = sign (ΔPij (τ)) × μ × F × √2 × ρk × | ΔPij(τ) |

<                 Gij(τ)                                                                    ,                                       (П6.15)

   |Σni=1 -----------  = 0

   |                ρk

    _

                 1

где ρk = ---- × [ρi + sign (ΔPij) × ρi + ρj - sign (ΔPij) × ρj]  .

                 2