摘要
本文讨论了与森林灭火建模相关的问题。提出了森林火灾与水障相互作用的数学公式。论证了火和植被结构的演变对蒸发过程的影响。得出的结论是,在存在茂密的低植被层的情况下使用水屏障对阻止火灾蔓延无效。
简介
对世界各地大型森林火灾的动态分析表明,现有的处理方法存在缺陷。成功扑灭森林火灾需要充分预测火灾的动态,并考虑到为扑灭火灾而采取的行动。
森林火灾有许多不同层次的数学模型。经验和半经验模型基于统计和实验数据的概括。有几篇论文致力于基于半经验模型对消防进行建模,如[7]。
森林火灾的物理模型提供了基于力学反应介质定律的详细情况,因此这些模型需要大量的计算。与经验模型不同,物理模型允许模拟火灾与灭火用水之间的详细交互作用,并根据正在进行的过程的物理规律找到最有效的灭火策略。
灭火效率取决于向火区提供足够数量的水以湿润和冷却森林燃料的能力。灭火场景主要取决于灭火设备在最大可能距离内将水输送至燃烧目标的能力。水流中液滴的初始尺寸和初始速度对于在燃烧区获得水至关重要。在中讨论了液滴蒸发和飞行轨迹的建模。本文致力于研究水幕的屏蔽特性。
根据《建筑物和隧道内火灾动力学》中讨论的使用自动喷水灭火系统在封闭空间和隧道内灭火与开放空间内的火灾动力学有很大不同。因此,在森林火灾期间,会发生明显的垂直对流。手稿显示,在水雾灭火的情况下,90-95%以上的水雾不会产生预期效果。
尽管有大量的论文,但考虑到物理和化学过程以及空气动力学的完整模型,扑灭森林火灾的问题仍未得到解决。在现有的扑灭森林火灾的方法中,应强调以下方法:从上方排放水,通过地面方式向火灾区域供水,并形成水屏障。现有的森林灭火规则通常基于计算所需水量的工程方法和观测结果。近年来扑灭大型火灾的经验表明,火灾的动态往往超出预测,并导致灾难性后果。因此,考虑关键物理、化学和水动力过程的森林火灾数值模拟问题仍然是一个重要和相关的问题。
本研究的目的是调查森林燃料密度和地形对水屏障有效性的影响。
方法
森林火灾的动态取决于可燃物质的密度分布和地形。本文研究了森林结构的两种变体:一级和两级。在这两种情况下,森林高度均为4m。一级森林被建模为整个高度上密度为2kg/m3的均质可燃介质。
两级森林的底层凋落物、干草和小灌木,高度为40cm,燃料密度为8kg/m3。两级森林上层高度为3.6m,密度为2kg/m3。问题的物理模型如图1所示。图1(a)显示了火灾在山区蔓延的情况下问题的几何结构。两级森林的结构如图1(b)所示。第3节的结果还显示了水障和地形。假设水障的宽度为40 cm,高度为4 m,这与森林的高度相对应。
水是扑灭森林火灾最常用的方法之一。剧烈沸腾发生在进入燃烧区时,这导致热能的吸收和水蒸气的释放。假设在高于水沸点的温度下,整个供应的能量都用于沸腾过程,流体本身的温度几乎保持恒定。因此,沸腾过程不可能像其他物理和化学过程那样将质量速度定义为水温的函数。
ρ
4
∂
φ
4
∂
t
=
−
R
4
+
f
(
x
,
z
,
t
)
(1)
\tag{1} \rho_4 \frac{\partial \varphi_4}{\partial t}=-R_4+f(x, z, t)
ρ4∂t∂φ4=−R4+f(x,z,t)(1)
我们假设自由水储备可以通过煮沸和补充来消耗,例如,在给定区域使用水炮的结果
根据森林火灾的物理模型(不考虑自由水沸腾的过程),温度由
∂
(
(
∑
i
=
1
4
ρ
i
φ
ℓ
c
p
t
+
ρ
5
c
p
5
)
T
)
∂
t
+
∂
(
ρ
5
c
p
5
U
T
)
∂
x
+
∂
(
ρ
5
c
p
5
W
T
)
∂
z
=
=
∂
∂
x
(
λ
i
∂
T
∂
x
)
+
∂
∂
z
(
λ
i
∂
T
∂
z
)
−
q
2
R
2
+
q
3
R
3
+
q
5
R
5
+
k
s
(
c
U
R
−
4
σ
(
T
4
−
T
e
4
)
)
(2)
\tag{2} \begin{aligned} & \frac{\partial\left(\left(\sum_{i=1}^4 \rho_i \varphi_{\ell} c_{p t}+\rho_5 c_{p 5}\right) T\right)}{\partial t}+\frac{\partial\left(\rho_5 c_{p 5} U T\right)}{\partial x}+\frac{\partial\left(\rho_5 c_{p 5} W T\right)}{\partial z}= \\ & =\frac{\partial}{\partial x}\left(\lambda_i \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial z}\left(\lambda_i \frac{\partial T}{\partial z}\right)-q_2 R_2+q_3 R_3+q_5 R_5+k_s\left(c U_R-4 \sigma\left(T^4-T_e^4\right)\right) \end{aligned}
∂t∂((∑i=14ρiφℓcpt+ρ5cp5)T)+∂x∂(ρ5cp5UT)+∂z∂(ρ5cp5WT)==∂x∂(λi∂x∂T)+∂z∂(λi∂z∂T)−q2R2+q3R3+q5R5+ks(cUR−4σ(T4−Te4))(2)
经过简单的变换和替换后,方程(2)将采用
T
′
(
∑
i
=
1
4
ρ
i
φ
i
c
p
i
+
ρ
5
c
p
s
)
+
T
∂
∂
t
(
∑
i
=
1
4
ρ
i
φ
i
c
p
i
+
ρ
5
c
p
5
)
+
∂
(
ρ
5
c
p
5
U
T
)
∂
x
+
∂
(
ρ
5
c
p
5
W
T
)
∂
z
=
=
∂
∂
x
(
λ
t
∂
T
∂
x
)
+
∂
∂
z
(
λ
t
∂
T
∂
z
)
−
q
2
R
2
+
q
3
R
3
+
q
5
R
5
+
k
s
(
c
U
R
−
4
σ
(
T
4
−
T
e
4
)
)
.
(3)
\tag{3} \begin{aligned} & T^{\prime}\left(\sum_{i=1}^4 \rho_i \varphi_i c_{p i}+\rho_5 c_{p s}\right)+T \frac{\partial}{\partial t}\left(\sum_{i=1}^4 \rho_i \varphi_i c_{p i}+\rho_5 c_{p 5}\right)+\frac{\partial\left(\rho_5 c_{p 5} U T\right)}{\partial x}+\frac{\partial\left(\rho_5 c_{p 5} W T\right)}{\partial z}= \\ & =\frac{\partial}{\partial x}\left(\lambda_t \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial z}\left(\lambda_t \frac{\partial T}{\partial z}\right)-q_2 R_2+q_3 R_3+q_5 R_5+k_s\left(c U_R-4 \sigma\left(T^4-T_e^4\right)\right) . \end{aligned}
T′(i=1∑4ρiφicpi+ρ5cps)+T∂t∂(i=1∑4ρiφicpi+ρ5cp5)+∂x∂(ρ5cp5UT)+∂z∂(ρ5cp5WT)==∂x∂(λt∂x∂T)+∂z∂(λt∂z∂T)−q2R2+q3R3+q5R5+ks(cUR−4σ(T4−Te4)).(3)
在方程(3)中, T ′ T^{\prime} T′的值决定了在没有蒸发过程的情况下的温度变化率。只要水的温度不低于沸点并且具有足够的能量来维持该过程,水就会沸腾。因此,理论上,以下情况是可能的:
a)存在非零水供应,沸腾的质量流率受到能量通量密度的限制,温度保持在Tb水平。在这种情况下,沸腾速率由下式确定。
R
4
=
1
q
4
T
′
(
∑
i
=
1
4
ρ
i
φ
i
c
p
i
+
ρ
5
c
p
5
)
.
(4)
\tag{4} R_4=\frac{1}{q_4} T^{\prime}\left(\sum_{i=1}^4 \rho_i \varphi_i c_{p i}+\rho_5 c_{p 5}\right) .
R4=q41T′(i=1∑4ρiφicpi+ρ5cp5).(4)
b)温度高于沸点,所有进入的水都被蒸发掉。在这种情况下,沸腾速率可以表示为
R 4 = f ( x , z , t ) (5) \tag{5} R_4=f(x, z, t) R4=f(x,z,t)(5)
c)如果温度低于沸点T<Tb或没有游离水,则不会发生沸腾过程。
R
4
=
0
(6)
\tag{6} R_4=0
R4=0(6)
该模型不考虑温度或浓度突然增加的情况,因为在这种情况下,模型中的沸腾速率将减小为δ函数。与[[15]]不同,本文使用了一种改进的自由水沸腾率计算方法
∂
(
(
∑
i
=
1
4
ρ
i
φ
i
c
p
i
+
ρ
5
c
p
5
)
T
)
∂
t
+
∂
(
ρ
5
c
p
5
U
T
)
∂
x
+
∂
(
ρ
5
c
p
5
W
T
)
∂
z
=
=
∂
∂
x
(
λ
t
∂
T
∂
x
)
+
∂
∂
z
(
λ
t
∂
T
∂
z
)
+
q
1
R
1
−
q
2
R
2
+
q
3
R
3
−
q
4
R
4
+
q
5
R
5
+
k
s
(
c
U
R
−
4
σ
(
T
4
−
T
e
4
)
)
(7)
\tag{7} \begin{aligned} & \frac{\partial\left(\left(\sum_{i=1}^4 \rho_i \varphi_i c_{p i}+\rho_5 c_{p 5}\right) T\right)}{\partial t}+\frac{\partial\left(\rho_5 c_{p 5} U T\right)}{\partial x}+\frac{\partial\left(\rho_5 c_{p 5} W T\right)}{\partial z}= \\ & =\frac{\partial}{\partial x}\left(\lambda_t \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial z}\left(\lambda_t \frac{\partial T}{\partial z}\right)+q_1 R_1-q_2 R_2+q_3 R_3-q_4 R_4+q_5 R_5+k_s\left(c U_R-4 \sigma\left(T^4-T_e^4\right)\right) \end{aligned}
∂t∂((∑i=14ρiφicpi+ρ5cp5)T)+∂x∂(ρ5cp5UT)+∂z∂(ρ5cp5WT)==∂x∂(λt∂x∂T)+∂z∂(λt∂z∂T)+q1R1−q2R2+q3R3−q4R4+q5R5+ks(cUR−4σ(T4−Te4))(7)
Q
=
(
1
−
α
c
)
R
1
+
R
2
+
R
4
+
M
c
M
1
R
3
(8)
\tag{8} Q=\left(1-\alpha_c\right) R_1+R_2+R_4+\frac{M_c}{M_1} R_3
Q=(1−αc)R1+R2+R4+M1McR3(8)
为了模拟添加自由水阶段,应引入方程(7)-(8)
为了确定屏障中水的初始体积分数,使用以下公式:
φ 4 b = φ 1 ρ 1 ρ 4 w b (9) \tag{9} \varphi_{4 b}=\varphi_1 \frac{\rho_1}{\rho_4} w_b φ4b=φ1ρ4ρ1wb(9)
根据上述模型,模拟了森林火灾期间水的补充和消耗过程。在这项工作中,假设存在一个水屏障,表示含有游离水的区域,该区域在加热到沸腾温度后,将所有提供给它的能量用于蒸发过程。在计算开始时,屏障的密度由内部水的质量分数确定。具体密度值在相应的计算中给出。
考虑分裂两个过程的交错模式的大粒子方法[16],以及用于解决该问题的空间坐标[17]。
大颗粒方法的优点是其稳定性和实施简单。这种方法的缺点是一级精度和显著的人工粘度。使用交错模式,可以模拟与起伏均匀性相关的角点。模式示意图如图2所示。在中心确定了与相应单元格相关的数量。
结果和讨论
根据上述模型,模拟了火与水屏障的相互作用。图-3.1-3.2中的温度显示为灰度颜色;箭头表示速度场;水屏障显示为实心黑色区域,其厚度对应于屏障中剩余水的比例沿整个森林高度的分布。
图-3.1显示了燃烧前缘与障碍物碰撞时一级森林中火灾的动态。平原地形的护栏密度为1.6 kg/m3,丘陵地形的护栏为2.5 kg/m3。可以清楚地看到,在非均质地形的情况下,火到达屏障的速度稍快。从计算结果中可以看出,在火灾穿过平原的传播过程中,火灾的前部在通过屏障后得到恢复,但底层不参与这一过程。当火势沿着山坡蔓延时,在到达背风坡之前,火势无法再次增加,因为火势主要在炉膛内水平燃烧,大部分受热物质从林冠中排出。如果火灾在7.6秒时蔓延到平原,可以看到屏障几乎没有受到影响。
在8秒时,可以看到燃烧区的温度显著降低,屏障部分沸腾。一旦屏障停止存在,火锋在8.4秒时再次升温。然而,屏障阻止了热量传播,所以森林燃料中的水分并没有蒸发。因此,火灾传播将在8秒时,可以看到燃烧区的温度显著降低,屏障部分沸腾。一旦屏障停止存在,火锋在8.4秒时再次升温。然而,屏障阻止了热量传播,所以森林燃料中的水分并没有蒸发。因此,火灾传播需要额外的能量来蒸发水分。缺乏能量导致火灾停止。
在图3.1的右栏中,我们可以看到火与山丘高原上的水屏障的碰撞动力学。与平原相比,高原火灾的典型坡度更大,因此在8.4秒内,障碍物不会完全消失。尽管火的能量足以蒸发掉所有的水,但火也会熄灭。
图-3.2显示了森林火灾在两级森林中分别通过丘陵和平原蔓延的演变。在这种情况下,防潮层的密度为20 kg/m3。在火灾初期,蔓延模式大致相同,不取决于地形。在这两种情况下,在10秒的时间内,被加热的气流流过屏障,几乎不会蒸发。在向风斜坡上蔓延的过程中,上锋的速度会增加。在平原上蔓延的情况下,在14秒的瞬间,火焰到达屏障,这导致其在上前方蒸发时的功耗。
数值模拟结果表明,对流效应和外部速度场的共同作用使得屏障的下部蒸发较弱,从而导致较低密度层的燃烧停止。
当考虑到火灾在山坡上蔓延时,也会有类似的效果,尽管由于屏障的位置,这会稍晚一些。在这两种情况下,由于材料密度较高和风速较低,火灾的下边缘都在上边缘之后。因此,几乎整个加热气流都被引导到水屏障。在22秒的瞬间,当水屏障蒸发时,显示出热气流。此时,主动燃烧区位于距离障碍物约2-4米的位置,这提供了足够的时间来克服障碍物。
结论
研究结果表明,在两层森林燃料的情况下,来自底层燃烧炉膛的对流气体流被引向水障,导致其沸腾,因此,水障灭火变得无效。在一级森林的情况下没有观察到这种影响。底层对水屏障效率的影响大于山丘。