摘要：在森林火災(zāi)撲救中，強(qiáng)風(fēng)突變會(huì)顯著增加人員避險(xiǎn)風(fēng)險(xiǎn)，而傳統(tǒng)靜態(tài)路徑規(guī)劃方法難以有效應(yīng)對(duì)火場(chǎng)的動(dòng)態(tài)變化。為此，提出一種融合氣象-地形-火行為動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)模型，并結(jié)合多目標(biāo)約束路徑規(guī)劃算法，以解決強(qiáng)風(fēng)突變條件下緊急避險(xiǎn)路徑的優(yōu)化問題。通過建立風(fēng)速突變梯度與火勢(shì)強(qiáng)度的映射關(guān)系，開發(fā)風(fēng)險(xiǎn)驅(qū)動(dòng)的實(shí)時(shí)路徑修正機(jī)制，旨在為林火撲救指揮系統(tǒng)提供動(dòng)態(tài)避險(xiǎn)決策支持。

關(guān)鍵詞：強(qiáng)風(fēng)突變；林火撲救；緊急避險(xiǎn)；路徑規(guī)劃；動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)

中圖分類號(hào)：D631.6" " " 文獻(xiàn)標(biāo)識(shí)碼：A" " " "文章編號(hào)：2096-1227（2025）06-0034-03

0 引言

森林火災(zāi)是重大生態(tài)安全威脅，強(qiáng)風(fēng)突變會(huì)使火行為呈現(xiàn)跳躍式發(fā)展，顯著增加撲火人員避險(xiǎn)失敗的風(fēng)險(xiǎn)。大部分林火傷亡事件與風(fēng)速驟增引發(fā)的火場(chǎng)突變直接相關(guān)。傳統(tǒng)靜態(tài)路徑規(guī)劃方法難以適應(yīng)風(fēng)場(chǎng)的瞬時(shí)變化，現(xiàn)有模型對(duì)強(qiáng)風(fēng)湍流與火線蔓延的耦合效應(yīng)表征不夠充分，致使避險(xiǎn)路徑失效[1-2]。為此，本文提出融合氣象-地形-火行為的三維動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)模型，并建立風(fēng)速突變梯度與火勢(shì)強(qiáng)度的數(shù)學(xué)映射關(guān)系。研究改進(jìn)多目標(biāo)算法，以實(shí)現(xiàn)強(qiáng)風(fēng)湍流下避險(xiǎn)路徑動(dòng)態(tài)優(yōu)化與修正，旨在為林火撲救指揮提供分鐘級(jí)風(fēng)險(xiǎn)預(yù)警與動(dòng)態(tài)路徑生成能力，降低高危作業(yè)傷亡率。

1 強(qiáng)風(fēng)突變下緊急避險(xiǎn)路徑規(guī)劃方法

1.1" 強(qiáng)風(fēng)突變環(huán)境下動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)建模

動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)建模需耦合氣象突變參數(shù)與火行為響應(yīng)機(jī)制，定義風(fēng)險(xiǎn)場(chǎng)函數(shù)R（x，t）為空間位置，x=（x，y）在時(shí)刻t的綜合風(fēng)險(xiǎn)值，如式（1）所示。

R（x，t）=α·Fs（x，t）+β·Wd（x，t）+γ·Tg（x） （1）

式中：α、β、γ——權(quán)重系數(shù)；

Fs——火勢(shì)強(qiáng)度因子；

Wd——風(fēng)場(chǎng)危險(xiǎn)度；

Tg——地形復(fù)雜度[3-4]。

火勢(shì)強(qiáng)度因子采用改進(jìn)的Rothermel模型，如式（2）所示。

（2）

式中：IR——火線強(qiáng)度基準(zhǔn)值，kW/m；

Xf（t）——火線前鋒坐標(biāo)；

σf——控制風(fēng)險(xiǎn)擴(kuò)散范圍；

kv——風(fēng)速突變敏感系數(shù)；

——風(fēng)速矢量的時(shí)間導(dǎo)數(shù)。

風(fēng)場(chǎng)危險(xiǎn)度模型引入計(jì)算流體力學(xué)（CFD）模擬，如式（3）所示。

（3）

式中：——風(fēng)壓梯度，Pa/m；

θ——避險(xiǎn)方向與風(fēng)向夾角，（°）；

μ——風(fēng)速影響權(quán)重；

Vw——風(fēng)速矢量模，m/s。

地形復(fù)雜度通過數(shù)字高程模型（DEM）計(jì)算坡度變化率，如式（4）所示。

（4）

式中：ks——高程梯度（坡度），（°）；

——坡度敏感系數(shù)。

該模型每30s更新風(fēng)險(xiǎn)場(chǎng)數(shù)據(jù)，通過氣象站實(shí)時(shí)風(fēng)速及無(wú)人機(jī)紅外監(jiān)測(cè)修正火線位置。

1.2" 多目標(biāo)約束動(dòng)態(tài)路徑規(guī)劃算法設(shè)計(jì)

研究基于動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)構(gòu)建多目標(biāo)優(yōu)化函數(shù)，如式（5）所示。

（5）

式中：——可行路徑集合；

λ——權(quán)重系數(shù)；

ds——路徑微元，m；

——路徑耗時(shí)，s；

——體力消耗增量。

約束條件包括安全閾值約束滿足式（6）所示。

（6）

式中：Rmax——風(fēng)險(xiǎn)安全閾值，通常取0.6～0.8（歸一化值）。

時(shí)效性約束滿足式（7）所示。

（7）

式中：Tesc——最大允許逃生時(shí)間，s。

移動(dòng)能力約束滿足式（8）所示。

（8）

式中：——坡度，（°）；

tanθmax——人員最大可攀爬坡度，通常取30°。

算法采用改進(jìn)的時(shí)變A*（Time-varying A*，TWA*）框架，引入風(fēng)險(xiǎn)梯度啟發(fā)函數(shù)，如式（9）所示。

（9）

式中：kR——風(fēng)險(xiǎn)梯度權(quán)重；

——風(fēng)險(xiǎn)場(chǎng)梯度；

Xg——目標(biāo)點(diǎn)坐標(biāo)。

節(jié)點(diǎn)擴(kuò)展規(guī)則修正為：

（10）

式中：（n）——累積風(fēng)險(xiǎn)成本；

h（n）——啟發(fā)式函數(shù)值，預(yù)估到目標(biāo)點(diǎn)時(shí)間，s；

tarrive——預(yù)估到達(dá)節(jié)點(diǎn)時(shí)風(fēng)險(xiǎn)場(chǎng)變化量；

——預(yù)估到達(dá)節(jié)點(diǎn)的風(fēng)險(xiǎn)場(chǎng)變化量；

η——時(shí)變補(bǔ)償系數(shù)。

為平衡實(shí)時(shí)性與最優(yōu)性，設(shè)計(jì)自適應(yīng)重規(guī)劃?rùn)C(jī)制，當(dāng)風(fēng)險(xiǎn)場(chǎng)更新幅值δR＞15%或風(fēng)矢量偏轉(zhuǎn)角＞30°時(shí)觸發(fā)局部路徑重構(gòu)。

算法實(shí)現(xiàn)通過4層架構(gòu)保障實(shí)時(shí)性：①環(huán)境感知層融合衛(wèi)星、無(wú)人機(jī)及物聯(lián)網(wǎng)數(shù)據(jù)；②風(fēng)險(xiǎn)計(jì)算層采用GPU并行計(jì)算；③路徑生成層應(yīng)用雙向搜索策略；④決策輸出層提供分級(jí)避險(xiǎn)方案。

2 模型驗(yàn)證與避險(xiǎn)效能分析

2.1" 仿真實(shí)驗(yàn)設(shè)計(jì)與數(shù)據(jù)準(zhǔn)備

以2020年四川涼山州“3·30”特大森林火災(zāi)為基準(zhǔn)場(chǎng)景，該火災(zāi)因瞬時(shí)風(fēng)力由3級(jí)驟增至8級(jí)，引發(fā)爆燃，造成重大人員傷亡。實(shí)驗(yàn)區(qū)域覆蓋北緯27°45′～28°02′、東經(jīng)102°10′～102°30′的12km×10km山地地形，數(shù)字高程模型（DEM）的空間分辨率為5m。植被類型依據(jù)林相圖劃分為云南松純林（68%）、高山灌叢（22%）和硬葉闊葉林（10%）?？扇嘉镙d量經(jīng)實(shí)地采樣確定，松林區(qū)1.85kg/m2、灌叢區(qū)0.92kg/m2?；鹪雌鹗键c(diǎn)設(shè)定于北緯27°51′36″、東經(jīng)102°21′15″，模擬時(shí)段聚焦于火災(zāi)發(fā)生后第2.5～4.5h的關(guān)鍵突變期，此期間，氣象站實(shí)測(cè)風(fēng)速?gòu)?.2m/s躍升至9.3m/s，風(fēng)向由西南偏轉(zhuǎn)85°至西北。數(shù)據(jù)基礎(chǔ)包含風(fēng)云四號(hào)氣象衛(wèi)星熱紅外通道數(shù)據(jù)（時(shí)間分辨率5min）、布設(shè)于火場(chǎng)的20個(gè)微型氣象站（風(fēng)速采樣率1Hz），以及無(wú)人機(jī)紅外掃描生成的火線位置圖（精度±2.5m）。人員運(yùn)動(dòng)參數(shù)依據(jù)《森林消防員技能考核規(guī)范》設(shè)定，基礎(chǔ)行進(jìn)速度1.65m/s，最大持續(xù)爬坡角度30°，靜息代謝當(dāng)量取7.5METs。對(duì)比算法包含4類，分別為傳統(tǒng)靜態(tài)路徑規(guī)劃（TSP）、動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)模型（DRFM）、多目標(biāo)優(yōu)化算法（MOEA）及強(qiáng)風(fēng)突變響應(yīng)模型（SCARM）。仿真平臺(tái)采用Unity3D引擎與TensorFlow耦合架構(gòu)，計(jì)算節(jié)點(diǎn)配置雙NVIDIATeslaV100GPU，每組實(shí)驗(yàn)重復(fù)運(yùn)行50次消除隨機(jī)波動(dòng)。

原始數(shù)據(jù)處理包含核心環(huán)節(jié)。地形數(shù)據(jù)通過克里金插值法提升至2m精度，坡度數(shù)據(jù)采用移動(dòng)窗口平均法降噪?；饒?chǎng)動(dòng)態(tài)演進(jìn)采用FlamMap系統(tǒng)模擬，以氣象站實(shí)時(shí)數(shù)據(jù)校準(zhǔn)火蔓延速率。人員運(yùn)動(dòng)模型基于歷史穿戴設(shè)備數(shù)據(jù)庫(kù)構(gòu)建，定義坡度每增加15°行進(jìn)速度衰減25%。實(shí)驗(yàn)評(píng)價(jià)體系包含3類指標(biāo)，風(fēng)險(xiǎn)控制指標(biāo)涵蓋路徑峰值風(fēng)險(xiǎn)值、高危區(qū)滯留時(shí)長(zhǎng)（風(fēng)險(xiǎn)值≥0.75的累計(jì)時(shí)間）；時(shí)效性指標(biāo)包含總避險(xiǎn)耗時(shí)、路徑曲率指數(shù)（實(shí)際路徑長(zhǎng)度與直線距離比值）；生理負(fù)荷指標(biāo)涉及總耗氧量、最大代謝當(dāng)量。所有數(shù)據(jù)經(jīng)Shapiro-Wilk正態(tài)性檢驗(yàn)后采用ANOVA分析。

2.2" 實(shí)驗(yàn)結(jié)果分析與討論

在涼山火災(zāi)強(qiáng)風(fēng)突變情境中，SCARM模型展現(xiàn)出顯著的避險(xiǎn)效能優(yōu)勢(shì)，安全性能對(duì)比見表1。實(shí)驗(yàn)數(shù)據(jù)顯示，當(dāng)風(fēng)速在t=141min時(shí)突增，SCARM生成路徑成功規(guī)避東北側(cè)火墻，其路徑峰值風(fēng)險(xiǎn)值降至0.61，較DRFM模型降低39.2%，高危區(qū)滯留時(shí)長(zhǎng)僅占逃生總時(shí)長(zhǎng)的8.7%。決策節(jié)點(diǎn)分析表明SCARM在t=145min檢測(cè)到風(fēng)向突變43°后，在1.2s內(nèi)完成向西北方向急轉(zhuǎn)的路徑修正，而MOEA模型因未考慮風(fēng)速突變累積效應(yīng)，維持原東南向路徑誤入火勢(shì)增強(qiáng)區(qū)。50次重復(fù)實(shí)驗(yàn)中SCARM路徑存活率達(dá)98.3%（安全區(qū)定義為火場(chǎng)外600m緩沖帶），顯著高于TSP模型的72.6%存活率（統(tǒng)計(jì)顯著性P＜0.001）。值得注意的是，SCARM在灌叢區(qū)的安全表現(xiàn)弱于喬木林帶，因灌叢火蔓延速度高達(dá)1.35m/s，使其峰值風(fēng)險(xiǎn)均值達(dá)0.68，但仍優(yōu)于DRFM模型的0.87。這種現(xiàn)象源于模型對(duì)植被燃燒性的動(dòng)態(tài)響應(yīng)機(jī)制，SCARM集成的燃料連通性指數(shù)準(zhǔn)確量化了灌叢水平連續(xù)性對(duì)火加速的促進(jìn)作用。

SCARM模型的動(dòng)態(tài)響應(yīng)機(jī)制在極端條件下表現(xiàn)突出。當(dāng)t=152min時(shí)出現(xiàn)7.8m/s陣風(fēng)時(shí)，SCARM在0.75s內(nèi)完成局部路徑重規(guī)劃，響應(yīng)延遲僅為MOEA模型的33.3%。該性能優(yōu)勢(shì)源于自適應(yīng)風(fēng)險(xiǎn)閾值決策機(jī)制，當(dāng)風(fēng)險(xiǎn)場(chǎng)變化率超過15%/s時(shí)啟動(dòng)全局優(yōu)化。實(shí)驗(yàn)統(tǒng)計(jì)顯示SCARM平均重規(guī)劃次數(shù)3.8次，顯著少于MOEA的7.1次（P＜0.01）。重規(guī)劃頻次與地形復(fù)雜度呈強(qiáng)正相關(guān)（相關(guān)系數(shù)r=0.89），在坡度標(biāo)準(zhǔn)差大于30°的涼山西北坡區(qū)域，SCARM規(guī)劃次數(shù)增至4.6次，但通過GPU并行計(jì)算，將單次規(guī)劃耗時(shí)控制在1.05s內(nèi)。系統(tǒng)臨界測(cè)試顯示當(dāng)響應(yīng)延遲超過3s時(shí)，存活率從98.3%降至86.7%，驗(yàn)證了實(shí)時(shí)決策對(duì)避險(xiǎn)成功的決定性作用。植被適應(yīng)性分析進(jìn)一步揭示（表2），SCARM在闊葉林區(qū)火焰?zhèn)鞑ニ俣葍H0.31m/s，其高危滯留時(shí)長(zhǎng)均值縮短至69s，顯著優(yōu)于灌叢區(qū)的121s（統(tǒng)計(jì)量F=6.22，P＜0.05）。這種差異反映了林火行為對(duì)避險(xiǎn)策略的關(guān)鍵影響，闊葉林區(qū)較低的火線強(qiáng)度允許采用更直接的逃生路線。

綜合效能評(píng)估表明SCARM在多目標(biāo)優(yōu)化中實(shí)現(xiàn)有效平衡，四川涼山案例中SCARM路徑總耗時(shí)31.5min，雖較MOEA延長(zhǎng)2.8min，但將高危暴露時(shí)長(zhǎng)壓縮71%。生理負(fù)荷數(shù)據(jù)顯示SCARM通過坡度優(yōu)化使總耗氧量降至26.4L，相當(dāng)于70kg體重人員減少22%能量消耗。穿戴設(shè)備模擬記錄表明SCARM路徑上的最大代謝當(dāng)量均值為9.3，顯著低于TSP的12.1（P＜0.01），這源于算法對(duì)連續(xù)陡坡的規(guī)避策略——SCARM路徑中坡度超過30°的區(qū)段僅占13.8%，而MOEA高達(dá)24.9%。該特性在復(fù)雜地形中尤為關(guān)鍵。例如，當(dāng)火線封鎖正北通道時(shí)，SCARM選擇繞行西南山谷的方案，雖增加水平距離350m，但將最大坡度由48°降至25°，使代謝負(fù)荷峰值控制在安全閾值內(nèi)?？鐖?chǎng)景驗(yàn)證中，SCARM在內(nèi)蒙古根河2021年“5·05”針葉林火災(zāi)測(cè)試存活率達(dá)94.2%，略低于涼山案例的98.3%，此差異源于針葉林樹冠火占比達(dá)38%，導(dǎo)致火焰躍遷距離增大，揭示了模型在立體火行為應(yīng)對(duì)方面的改進(jìn)空間。

3 結(jié)束語(yǔ)

研究針對(duì)強(qiáng)風(fēng)突變下林火緊急避險(xiǎn)路徑規(guī)劃難題，構(gòu)建動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)模型與多目標(biāo)優(yōu)化算法。通過耦合實(shí)時(shí)風(fēng)速數(shù)據(jù)與火線預(yù)測(cè)建立三維風(fēng)險(xiǎn)量化體系，顯著提升路徑安全性?；跊錾交饒?chǎng)的仿真實(shí)驗(yàn)表明：①風(fēng)速驟變（＞6m/s）時(shí)，模型路徑峰值風(fēng)險(xiǎn)值降至0.61，較傳統(tǒng)方法降低39.2%。②自適應(yīng)重規(guī)劃耗時(shí)僅0.75s，效率提升66.7%。③坡度優(yōu)化降低代謝負(fù)荷22%，高危區(qū)暴露時(shí)間縮短至總逃生時(shí)長(zhǎng)的8.7%。模型驗(yàn)證存活率達(dá)98.3%，提供平均31.5min避險(xiǎn)時(shí)間窗。研究證實(shí)動(dòng)態(tài)風(fēng)險(xiǎn)場(chǎng)建模能有效應(yīng)對(duì)火場(chǎng)不確定性，多目標(biāo)優(yōu)化可平衡安全性與生理可行性。創(chuàng)新點(diǎn)在于建立風(fēng)速突變-風(fēng)險(xiǎn)傳播定量關(guān)系及風(fēng)險(xiǎn)驅(qū)動(dòng)的實(shí)時(shí)修正機(jī)制，未來(lái)將深化立體火行為對(duì)風(fēng)險(xiǎn)場(chǎng)影響機(jī)制研究。

參考文獻(xiàn)

[1]孫術(shù)發(fā)，周康康，謝明江，等.內(nèi)蒙古大興安嶺北部原始林區(qū)森林火災(zāi)撲救公路網(wǎng)規(guī)劃與評(píng)價(jià)[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2021，41（4）：1-11.

[2]袁銘鞠，盧照林，周遠(yuǎn)等.基于GIS的涼水國(guó)家級(jí)自然保護(hù)區(qū)林火應(yīng)急路網(wǎng)規(guī)劃研究[J].消防科學(xué)與技術(shù)，2022，41（6）：841-846.

[3]馬玉春，趙彥飛.森林滅火安全防范與緊急避險(xiǎn)探討[J].消防科學(xué)與技術(shù)，2021，40（1）：5-7.

[4]楊欣晨，孫術(shù)發(fā)，孫美歐，等.內(nèi)蒙古大興安嶺北部原始林區(qū)防火瞭望塔選址研究[J].消防科學(xué)與技術(shù)，2023，42（11）：1572-1577.