劉暉，向諄，周恩澤，周洋卓君，黃宇，向坤軒，周游

（1.長(zhǎng)沙理工大學(xué)電氣與信息工程學(xué)院智能電網(wǎng)運(yùn)行與控制湖南省重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)沙 410114；2.廣東電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，廣州 510080）

0 引言

近年來，隨著“退耕還林”政策的實(shí)施，以及人民生產(chǎn)生活用火習(xí)俗和各種極端天氣等因素影響，輸電走廊附近易爆發(fā)大范圍山火災(zāi)害，嚴(yán)重影響電網(wǎng)安全穩(wěn)定運(yùn)行［1-3］。山火分布規(guī)律研究表明，輸電走廊山火具有時(shí)空分布不均衡特性［4-6］。在時(shí)間上呈現(xiàn)較強(qiáng)的季節(jié)和時(shí)節(jié)性特點(diǎn)，每年9 月至翌年4 月是相對(duì)高發(fā)期，并隨著春節(jié)、清明等民俗節(jié)氣因素而呈現(xiàn)逐日增長(zhǎng)的爆發(fā)性趨勢(shì)；在空間上具有明顯地域性特征，靠近人類生活和交通運(yùn)輸區(qū)域的輸電走廊區(qū)段通常發(fā)生山火的概率也會(huì)高于偏遠(yuǎn)地區(qū)。開展輸電走廊山火時(shí)空分布研究，確定特定時(shí)段的山火高發(fā)生風(fēng)險(xiǎn)區(qū)域，有助于運(yùn)維單位開展針對(duì)性的防山火工作，提升電網(wǎng)安全穩(wěn)定運(yùn)行性。

早期氣象和林業(yè)部門根據(jù)天氣因素和歷史火災(zāi)數(shù)據(jù)建立森林火險(xiǎn)氣象指數(shù)，預(yù)測(cè)大范圍區(qū)域火險(xiǎn)等級(jí)［7-10］，但是該類方法不適用于目標(biāo)需求更加精細(xì)化的輸電走廊山火預(yù)測(cè)。文獻(xiàn)［11］在森林火險(xiǎn)氣象指數(shù)的基礎(chǔ)上，引入了地表燃燒因子和歷史火情因子，構(gòu)建圖模型評(píng)估輸電線路山火風(fēng)險(xiǎn)等級(jí)。文獻(xiàn)［12］綜合考慮氣象、衛(wèi)星監(jiān)測(cè)熱點(diǎn)、工農(nóng)業(yè)用火和隱患點(diǎn)要素，通過賦予不同的權(quán)重計(jì)算山火預(yù)警值實(shí)現(xiàn)輸電走廊山火風(fēng)險(xiǎn)評(píng)估和預(yù)警。但是上述方法中的各要素權(quán)重由專家主觀確定，受人為影響較大。為降低專家主觀經(jīng)驗(yàn)不足造成的山火風(fēng)險(xiǎn)誤判，文獻(xiàn)［13］結(jié)合電網(wǎng)歷史監(jiān)測(cè)數(shù)據(jù)，利用BP 神經(jīng)網(wǎng)絡(luò)評(píng)估輸電線路山火風(fēng)險(xiǎn)。文獻(xiàn)［14-15］選取人為、氣象、地形和植被因子，構(gòu)建Logistic 回歸模型實(shí)現(xiàn)了輸電走廊山火風(fēng)險(xiǎn)評(píng)估。文獻(xiàn)［16-17］采用樸素貝葉斯網(wǎng)絡(luò)進(jìn)行山火風(fēng)險(xiǎn)評(píng)估，繪制了山火風(fēng)險(xiǎn)分布圖。但是這些模型通常要求因子和要素之間相互獨(dú)立，且所選變量為年平均值，未考慮氣象等因子隨時(shí)間和季節(jié)的動(dòng)態(tài)變化對(duì)山火發(fā)生的影響。綜上，目前輸電線路山火風(fēng)險(xiǎn)評(píng)估和預(yù)測(cè)存在兩方面問題，一是構(gòu)建模型過程中變量對(duì)山火發(fā)生事件的影響差異主要以主觀或客觀權(quán)重來體現(xiàn)，沒有考慮山火影響因子之間的相互耦合關(guān)系；二是模型變量通常選擇的是當(dāng)下時(shí)刻的值或者年平均分布情況，然而山火發(fā)生的概率除了和當(dāng)下的氣象條件有關(guān)以外，還和過去幾天的氣象和近期的人類活動(dòng)變化有關(guān)，現(xiàn)有的因子變量無法體現(xiàn)山火發(fā)生的時(shí)空聚集性。

為此，本文提出了一種基于動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)（dynamic Bayesian network，DBN）的山火風(fēng)險(xiǎn)時(shí)空分布評(píng)估方法。首先，綜合選取和收集了人為、氣象、植被和地理四大類共計(jì)17 個(gè)山火因子，利用隨機(jī)森林算法篩選最優(yōu)因子集。然后利用靜態(tài)貝葉斯網(wǎng)絡(luò)消除因子之間耦合關(guān)系的基礎(chǔ)上，建立了動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)模型，對(duì)比月時(shí)間尺度下靜態(tài)貝葉斯網(wǎng)絡(luò)和動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)模型對(duì)山火風(fēng)險(xiǎn)評(píng)估的影響，探索引入過去的因素對(duì)于山火風(fēng)險(xiǎn)評(píng)估的影響，并研究不同時(shí)間尺度對(duì)的動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)性能的影響。最后，利用最優(yōu)山火評(píng)估模型繪制輸電走廊山火風(fēng)險(xiǎn)分布圖驗(yàn)證了模型的工程適用性。

1 輸電走廊山火風(fēng)險(xiǎn)時(shí)空評(píng)估框架

為了更好地應(yīng)對(duì)突發(fā)性的山火災(zāi)害，本文提出了一種輸電走廊山火風(fēng)險(xiǎn)時(shí)空評(píng)估方法，主要分為數(shù)據(jù)收集與預(yù)處理、因子優(yōu)選和山火風(fēng)險(xiǎn)分布評(píng)估3個(gè)部分，總體框架如圖1所示。

圖1 輸電走廊山火風(fēng)險(xiǎn)時(shí)空評(píng)估框架Fig.1 Framework for assessing the spatial-temporal distribution of wildfire risk in transmission line corridors

1）數(shù)據(jù)收集與預(yù)處理。首先，從氣象部門和電力部門收集得到氣象、人為、植被和地理數(shù)據(jù)，采用箱線圖檢測(cè)數(shù)據(jù)集中的異常值并對(duì)空缺值進(jìn)行填充，得到山火因子數(shù)據(jù)集。然后，為進(jìn)行不同時(shí)間的山火風(fēng)險(xiǎn)評(píng)估，將歷史火點(diǎn)按照一定的時(shí)間間隔進(jìn)行劃分，得到不同時(shí)間的火點(diǎn)分布。最后，將研究區(qū)域按照經(jīng)緯度進(jìn)行網(wǎng)格劃分，其中火點(diǎn)所在網(wǎng)格作為火點(diǎn)樣本點(diǎn)，其他網(wǎng)格則作為非火點(diǎn)樣本點(diǎn)，獲得不同空間上的火點(diǎn)分布。利用地理信息軟件將山火特征數(shù)據(jù)提取至火點(diǎn)和非火點(diǎn)樣本點(diǎn)，構(gòu)成原始的山火樣本數(shù)據(jù)集。

2）因子優(yōu)選。原始的山火樣本數(shù)據(jù)集因子眾多，過多的冗余因子可能會(huì)帶來不必要的噪聲影響模型性能。因此采用隨機(jī)森林算法進(jìn)行因子優(yōu)選，通過Bootstrap重抽樣技術(shù)將數(shù)據(jù)集劃分為訓(xùn)練數(shù)據(jù)集和袋外數(shù)據(jù)集，在訓(xùn)練集上同步訓(xùn)練決策樹。然后逐一對(duì)因子施加擾動(dòng)并計(jì)算袋外誤差，得到各因子的重要度排序，得到最優(yōu)因子集。

3）山火風(fēng)險(xiǎn)分布評(píng)估。將篩選得到的山火因子數(shù)據(jù)集離散化，作為動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)模型的輸入，采用極大似然估計(jì)法學(xué)習(xí)初始網(wǎng)絡(luò)和轉(zhuǎn)移網(wǎng)絡(luò)的條件概率。然后以受試者工作特性曲線（receiver operating characteristic，ROC）和曲線下方面積（area under curve，AUC）作為評(píng)價(jià)指標(biāo)，評(píng)估模型在測(cè)試集上的預(yù)測(cè)效果。

2 山火影響因子分析與因子優(yōu)選

2.1 山火影響因子分析

可燃物、助燃物和著火源是火災(zāi)發(fā)生必備的三要素［18］，其中可燃物和助燃物為自然因子，著火源包括人為火源和自然火源。早期的森林火災(zāi)預(yù)測(cè)預(yù)報(bào)研究［19］主要從自然因子出發(fā)，考慮氣象和植被等因子，建立綜合評(píng)價(jià)指標(biāo)對(duì)大范圍的森林區(qū)域進(jìn)行籠統(tǒng)的火險(xiǎn)等級(jí)預(yù)測(cè)，該方法準(zhǔn)確性主要依賴于氣象數(shù)據(jù)的精確度。而輸電線路大多穿梭于人類活動(dòng)較為頻繁的山區(qū)和田地，線路周邊居民生產(chǎn)生活以及祭祖用火頻繁、且地形和植被等因子相互關(guān)聯(lián)，山火爆發(fā)隨機(jī)性強(qiáng)，僅考慮自然因子的林火預(yù)報(bào)難以實(shí)現(xiàn)對(duì)小范圍輸電線路山火的準(zhǔn)確預(yù)測(cè)。2017年國(guó)家電網(wǎng)公司頒布《架空輸電線路山火分布圖繪制導(dǎo)則》［20］，采用了歷史火點(diǎn)密度表征人為火源分布情況，同時(shí)考慮植被燃燒危害等級(jí)評(píng)估輸電線路山火風(fēng)險(xiǎn)。但該方法未考慮氣象和地形等助燃物要素，難以保證評(píng)估時(shí)效性。因此在評(píng)估輸電走廊山火時(shí)，應(yīng)該從火災(zāi)三要素出發(fā)，綜合考慮影響山火發(fā)生的氣象、人為、植被和地形等因子。

由于山火災(zāi)害還具有一定的時(shí)間變化特性，上一時(shí)刻的氣象信息和發(fā)生火災(zāi)情況會(huì)對(duì)當(dāng)前時(shí)刻山火發(fā)生造成一定影響。因此，本文綜合考慮氣象、人為、植被和地形4 類因子，通過靜態(tài)因子和動(dòng)態(tài)因子相結(jié)合來刻畫山火時(shí)間相關(guān)性。其中選取了歷史火點(diǎn)密度、人口密度、GDP、線路到公路距離、線路到居民點(diǎn)距離、植被類型、NDVI、可燃物載量、土地利用類型、海拔、坡度和坡向12 個(gè)基本短時(shí)間內(nèi)不隨著時(shí)間變化的因子作為靜態(tài)因子；并選擇目標(biāo)時(shí)段內(nèi)的連續(xù)無雨日、平均溫度、平均降水量、平均濕度和平均風(fēng)速5 個(gè)隨時(shí)間變化的動(dòng)態(tài)因子來反映時(shí)間變化對(duì)山火發(fā)生風(fēng)險(xiǎn)影響。構(gòu)建DBN 模型，通過證據(jù)推理計(jì)算山火發(fā)生的概率，評(píng)估山火風(fēng)險(xiǎn)。

2.2 基于隨機(jī)森林的因子優(yōu)選

通常為全面挖掘山火發(fā)生規(guī)律，在構(gòu)建模型時(shí)傾向于采用較多的因子。而引入過多的因子會(huì)帶來“維數(shù)災(zāi)”問題，導(dǎo)致模型復(fù)雜度增大，訓(xùn)練時(shí)間成本增加，甚至引入噪聲降低模型精度；與此同時(shí)，隨著變量的增加，貝葉斯網(wǎng)絡(luò)的參數(shù)學(xué)習(xí)所需樣本量也會(huì)急劇增加，給模型參數(shù)評(píng)估帶來困難。為此，本文采用隨機(jī)森林特征選擇算法［21-22］從原始的高維山火數(shù)據(jù)集中篩選出最優(yōu)因子集，提高模型的計(jì)算效率和擬合精度。

隨機(jī)森林是一種以決策樹為基學(xué)習(xí)器的集成監(jiān)督學(xué)習(xí)模型。隨機(jī)森林通過Bootstrap自助采樣法和節(jié)點(diǎn)隨機(jī)分裂技術(shù)建立決策樹，以少數(shù)服從多數(shù)的投票方式集成決策樹輸出結(jié)果。通過計(jì)算每個(gè)因子在每棵樹上的貢獻(xiàn)程度，可以評(píng)估因子對(duì)山火發(fā)生事件的重要程度，對(duì)應(yīng)的算法步驟如下。

步驟1：從原始樣本集中隨機(jī)有放回地抽取K個(gè)訓(xùn)練數(shù)據(jù)集{D1，D2，…，DK}和K個(gè)對(duì)應(yīng)的袋外數(shù)據(jù)集{B1，B2，…，BK}；

步驟2：設(shè)置k=1，在訓(xùn)練集Dk上建立決策樹，使用決策樹對(duì)Bk進(jìn)行預(yù)測(cè)，統(tǒng)計(jì)袋外數(shù)據(jù)預(yù)測(cè)誤差，記為；

步驟3：對(duì)Bk中的因子Xi，i=1，2，…，n施加擾動(dòng)，統(tǒng)計(jì)袋外數(shù)據(jù)預(yù)測(cè)誤差；

步驟4：對(duì)于k=2，3，…，K，重復(fù)步驟2、3；

步驟5：利用式（1）度量因子Xi的重要度；

3 基于動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)的山火風(fēng)險(xiǎn)評(píng)估模型

3.1 貝葉斯網(wǎng)絡(luò)和動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)

貝葉斯網(wǎng)絡(luò)（Bayesian network，BN）是一種用于概率推理的圖形化數(shù)學(xué)模型，它利用有向無環(huán)圖描述變量間的依賴關(guān)系，在給定觀測(cè)樣本數(shù)據(jù)或?qū)＜蚁闰?yàn)知識(shí)下，BN 能夠利用貝葉斯概率理論對(duì)目標(biāo)變量進(jìn)行概率推理，在含有多變量的不確定性問題中得到廣泛應(yīng)用。DBN 在BN 的基礎(chǔ)上引入了時(shí)序概念，通過刻畫隨機(jī)變量的時(shí)序轉(zhuǎn)移和因果關(guān)系，對(duì)動(dòng)態(tài)時(shí)變隨機(jī)過程進(jìn)行推理，在負(fù)荷預(yù)測(cè)、數(shù)據(jù)鏈路預(yù)測(cè)、海底管道失效預(yù)測(cè)［23-25］等動(dòng)態(tài)推理預(yù)測(cè)領(lǐng)域具有成熟應(yīng)用。

DBN 可被定義為(B1，B→)，B1為初始時(shí)刻的BN，定義了變量之間的因果關(guān)系和概率分布，B→為轉(zhuǎn)移網(wǎng)絡(luò)，它將前后相鄰的時(shí)間片的BN 聯(lián)結(jié)起來，描述了2 個(gè)相鄰時(shí)刻變量間的轉(zhuǎn)移概率P(Xt|Xt-1)，即：

圖2 DBN網(wǎng)絡(luò)拓?fù)銯ig.2 Topology of DBN network

通常為簡(jiǎn)化動(dòng)態(tài)模型推理的復(fù)雜度，假設(shè)DBN 滿足一階馬爾科夫和轉(zhuǎn)移概率不變兩個(gè)前提假設(shè)［26］。其中一階馬爾科夫假設(shè)要求節(jié)點(diǎn)之間的有向邊必須位于同一時(shí)間片或相鄰時(shí)間片之間，不能跨越時(shí)間片，即因子變化只與上一個(gè)時(shí)間片的狀態(tài)有關(guān)，與更早之前的時(shí)間片的狀態(tài)無關(guān)；轉(zhuǎn)移概率不變假設(shè)是指B→的轉(zhuǎn)移概率P(Xt|Xt-1)不隨時(shí)間變化?；谏鲜黾僭O(shè)，可將DBN 拓展至第T個(gè)時(shí)間片，得到DBN 任一節(jié)點(diǎn)的聯(lián)合概率，如式（3）所示。

3.2 基于DBN的山火風(fēng)險(xiǎn)評(píng)估

構(gòu)建DBN 山火風(fēng)險(xiǎn)評(píng)估模型時(shí)，包含網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)和參數(shù)學(xué)習(xí)兩個(gè)步驟。網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)指根據(jù)專家經(jīng)驗(yàn)或者樣本數(shù)據(jù)搜索初始網(wǎng)絡(luò)結(jié)構(gòu)B1，并對(duì)B1進(jìn)行時(shí)間延拓得到DBN 結(jié)構(gòu)，本文采用專家經(jīng)驗(yàn)法，通過分析山火影響因子之間的因果關(guān)系獲得合適的山火風(fēng)險(xiǎn)評(píng)估DBN 結(jié)構(gòu)；然后基于時(shí)間序列樣本訓(xùn)練得到DBN 各節(jié)點(diǎn)條件概率；最后通過DBN推理算法計(jì)算樣本點(diǎn)的山火發(fā)生概率。

3.2.1 模型節(jié)點(diǎn)參數(shù)的計(jì)算

DBN 節(jié)點(diǎn)參數(shù)即為各節(jié)點(diǎn)的條件概率表（conditional probability table，CPT），反映了各影響因子之間的概率依賴關(guān)系。CPT 通?？捎蓪＜医?jīng)驗(yàn)直接給出或者基于樣本數(shù)據(jù)訓(xùn)練得到。由于山火災(zāi)害爆發(fā)隨機(jī)性強(qiáng)，專家經(jīng)驗(yàn)給定的參數(shù)難以適用于各因素相互耦合的山火預(yù)測(cè)。因此本文基于歷史時(shí)間序列樣本學(xué)習(xí)節(jié)點(diǎn)的CPT，充分挖掘節(jié)點(diǎn)之間的依賴關(guān)系。

在山火歷史數(shù)據(jù)集觀測(cè)完全的條件下，采用統(tǒng)計(jì)學(xué)中的極大似然法［27］估計(jì)網(wǎng)絡(luò)節(jié)點(diǎn)的CPT。設(shè)樣本集合為D={d1，d2，…，dm}建立對(duì)數(shù)似然函數(shù)為：

式中：m為樣本數(shù)；n為時(shí)間片內(nèi)節(jié)點(diǎn)數(shù)；T為時(shí)間片數(shù)。定義θi，j，k，t=，則式（4）可化為：

式中j、k為節(jié)點(diǎn)取值，本文中取離散值1、2、3、4。Ii，j，k，l，t為：

利用拉格朗日乘數(shù)法求取L的極大值，得到參數(shù)θi，j，k，t的估計(jì)值為：

3.2.2 山火發(fā)生概率推理

確定模型網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)后，通過DBN 推理對(duì)山火發(fā)生概率進(jìn)行預(yù)測(cè)。DBN 推理方法可分為精確推理和近似推理2 種方法，本文采用精確推理中的直接計(jì)算推理算法［28］，該方法無需進(jìn)行復(fù)雜的圖形變換，根據(jù)貝葉斯公式和條件獨(dú)立性假設(shè)即可實(shí)現(xiàn)后驗(yàn)概率推理，推理計(jì)算效率較高。

對(duì)于T個(gè)時(shí)間片的網(wǎng)絡(luò)，假定影響山火發(fā)生的人為、氣象、植被和地理因子分別為變量R、Q、Z和D，山火為變量F，則在觀測(cè)值｛R1：T，Q1：T，Z1：T，D1：T，F(xiàn)1：T-1｝下F在時(shí)間片T上的后驗(yàn)概率為：

由條件獨(dú)立性假設(shè)，聯(lián)合分布概率為：

結(jié)合式（8）—（9）可以實(shí)現(xiàn)第T個(gè)時(shí)間片山火發(fā)生概率預(yù)測(cè)。

3.3 模型評(píng)估指標(biāo)

山火預(yù)測(cè)是一種二分類問題，因此采用混淆矩陣［29］來度量模型泛化性能，如表1 所示。真正類數(shù)（true positive，TP）記作NTP，表示實(shí)際有火網(wǎng)格樣本中被正確預(yù)測(cè)為有火的數(shù)量；假負(fù)類數(shù)（false negative，F(xiàn)N）記作NFN，表示實(shí)際有火網(wǎng)格樣本中被錯(cuò)誤預(yù)測(cè)為無火的數(shù)量；假正類數(shù)（false positive，F(xiàn)P）記作NFP，表示實(shí)際無火網(wǎng)格樣本中被錯(cuò)誤預(yù)測(cè)為有火的數(shù)量；真負(fù)類數(shù)（true negative，TN）記作NTN，表示實(shí)際無火網(wǎng)格樣本中被正確預(yù)測(cè)為無火的個(gè)數(shù)。

表1 混淆矩陣Tab.1 Definition of confusion matrix

NTP和NTN為模型預(yù)測(cè)正確的樣本數(shù)，其占比越高，預(yù)測(cè)效果越好。為綜合評(píng)估預(yù)測(cè)性能，引入準(zhǔn)確率（accuracy）和召回率（recall），如式（10）-（11）所示：

式中：A為全體樣本中被正確預(yù)測(cè)的比例；R為實(shí)際有火樣本中被正確預(yù)測(cè)的比例，R越高意味著有火預(yù)測(cè)越準(zhǔn)確。

采用ROC 曲線和AUC 值來衡量模型的分類性能［30］。ROC 曲線繪制方法為：將樣本預(yù)測(cè)值降序排列，分別以每個(gè)樣本的預(yù)測(cè)值作為分類閾值，將預(yù)測(cè)值大于或等于分類閾值的樣本判斷為正類，反之為負(fù)類；然后根據(jù)式（12）和式（13）計(jì)算真正類率（true positive rate，TPR）記作RTP和假正類率（false positive rate，F(xiàn)PR）記作RFP并分別以其為橫坐標(biāo)和縱坐標(biāo)繪制曲線。ROC曲線越靠近圖形左上方，表明模型分類效果越好。

AUC 值為ROC 曲線在區(qū)間［0，1］的面積，表示模型根據(jù)預(yù)測(cè)值將正樣本排在負(fù)樣本前面的概率，反映了模型對(duì)樣本的排序分類能力。當(dāng)AUC大于0.5 時(shí)，表示模型具備分類能力，AUC 大于0.8時(shí)，模型具有優(yōu)良的分類性能。

4 算例分析

4.1 研究區(qū)域及數(shù)據(jù)預(yù)處理

本文選擇華南地區(qū)廣東省作為研究區(qū)域。2015—2020 年間廣東省發(fā)生森林火災(zāi)3 584 起，空間分布如圖3 所示。全省總面積17.97 萬平方千米，地形以山地、丘陵為主。廣東省大部分地區(qū)屬亞熱帶季風(fēng)氣候，全年雨熱同期，適宜植被生長(zhǎng)。近年來森林覆蓋不斷擴(kuò)張，已達(dá)58.7%。2021 年，廣東常住人口總量達(dá)12 684 萬人。密集的人口分布，適宜的氣候條件和豐富的植被使其易受山火災(zāi)害侵襲。

圖3 2015—2020年廣東省森林火災(zāi)空間分布Fig.3 Spatial distribution of wildfires in Guangdong Province,2015—2020

為建立研究區(qū)域山火預(yù)測(cè)模型，收集了廣東省2010—2020年的歷史火點(diǎn)數(shù)據(jù)，包括火災(zāi)發(fā)生時(shí)間和經(jīng)緯度信息。由于超過90%的山火災(zāi)害由人為野外用火造成［31］，因此將2010—2014 年的歷史火點(diǎn)用于計(jì)算歷史火點(diǎn)密度，綜合反映研究區(qū)域人為用火分布狀況。為全面考慮山火影響因素，從氣象、人為、植被和地理4個(gè)方面收集了17個(gè)山火因子數(shù)據(jù)，如表2 所示。其中氣象因子以及可燃物載量數(shù)據(jù)由國(guó)家衛(wèi)星氣象中心提供，人為、植被和地理因子數(shù)據(jù)從中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（https：//www.resdc.cn/）下載得到，分辨率均為1 km×1 km。

表2 輸電走廊山火因子Tab.2 Wildfire factors of transmission line corridors

根據(jù)山火災(zāi)害在時(shí)間上呈現(xiàn)月度分布不均衡這一特性，以月為時(shí)間單位構(gòu)建山火樣本集。首先將2015—2020年的火點(diǎn)以月為單位進(jìn)行時(shí)間劃分，生成逐月火點(diǎn)時(shí)間序列數(shù)據(jù)集。然后為實(shí)現(xiàn)網(wǎng)格化收集數(shù)據(jù)，利用地理信息軟件將研究區(qū)域按照經(jīng)緯度劃分成0.05 °×0.05 °的網(wǎng)格，并按照目標(biāo)時(shí)間間隔計(jì)算落在網(wǎng)格內(nèi)的火點(diǎn)數(shù)，將火點(diǎn)數(shù)大于等于1 的網(wǎng)格作為有火網(wǎng)格，否則為無火網(wǎng)格。由于氣象數(shù)據(jù)為廣東省氣象站點(diǎn)數(shù)據(jù)，本文采用克里金插值法成1 km×1 km 分辨率柵格數(shù)據(jù)。最后，逐一提取網(wǎng)格山火因子數(shù)據(jù)，生成時(shí)間序列的有火和無火網(wǎng)格樣本集。

4.2 山火因子篩選

按照2.2 節(jié)所述方法對(duì)17 個(gè)山火因子進(jìn)行篩選，得到各因子貢獻(xiàn)度排序以及累計(jì)貢獻(xiàn)度如圖4所示。前9 個(gè)山火因子累計(jì)貢獻(xiàn)度高達(dá)95.07%［32］，后8 個(gè)因子貢獻(xiàn)度小，可視為冗余因子，因此選取前9個(gè)因子作為模型輸入數(shù)據(jù)。

圖4 各因子貢獻(xiàn)度及累計(jì)貢獻(xiàn)度Fig.4 Various factors VIM and cumulative VIM

4.3 靜態(tài)模型與動(dòng)態(tài)模型對(duì)比

4.3.1 網(wǎng)絡(luò)構(gòu)建

在短期內(nèi)，氣象因子時(shí)序關(guān)聯(lián)性強(qiáng)，隨時(shí)間變化范圍大，其在上一時(shí)刻和當(dāng)前時(shí)刻的狀態(tài)會(huì)持續(xù)對(duì)植被可燃性產(chǎn)生影響，具有時(shí)間累計(jì)效應(yīng)。因此在模型中引入時(shí)間因素，將氣象、植被等動(dòng)態(tài)因子上下兩個(gè)時(shí)刻的狀態(tài)結(jié)合起來。DBN 模型一方面基于概率理論和圖論耦合氣象、人為、植被和地理等多種因子，另一方面利用有向邊將動(dòng)態(tài)變量在時(shí)間軸上進(jìn)行連接，通過轉(zhuǎn)移概率挖掘時(shí)序關(guān)聯(lián)性，從而融入時(shí)間因素對(duì)山火發(fā)生概率進(jìn)行推理。本文構(gòu)建2時(shí)間片的DBN模型，基于各因子在上一時(shí)刻和當(dāng)前時(shí)刻的狀態(tài)對(duì)山火發(fā)生概率進(jìn)行推理，以反映當(dāng)前時(shí)刻山火發(fā)生的難易程度。

為探索引入時(shí)間因素對(duì)山火發(fā)生風(fēng)險(xiǎn)評(píng)估準(zhǔn)確性的影響，構(gòu)建靜態(tài)BN 模型，與2 個(gè)相鄰時(shí)間片的DBN 模型形成對(duì)比。首先以月為時(shí)間尺度，構(gòu)建靜態(tài)BN 網(wǎng)絡(luò)結(jié)構(gòu)，然后將BN 拓展至2時(shí)間片的DBN模型。

基于篩選得到的9 個(gè)山火因子構(gòu)建靜態(tài)BN。首先，對(duì)各因子進(jìn)行皮爾遜相關(guān)性分析，繪制相關(guān)性熱力圖，如圖5所示。

圖5 山火影響因子的相關(guān)性熱力圖Fig.5 Wildfire factors correlation heat map

相關(guān)系數(shù)絕對(duì)值大于0.1 的兩個(gè)因子之間存在一定的相關(guān)關(guān)系［33］。例如海拔與人口密度、溫度和火點(diǎn)密度的相關(guān)系數(shù)分別為-0.24、-0.18 和0.33，這是因?yàn)楹０卧礁撸鄳?yīng)的溫度越低，因此人類通常生活在海拔較低的地區(qū)，人口密度越大，則由人類活動(dòng)產(chǎn)生的歷史火點(diǎn)密度越大，該地區(qū)越容易發(fā)生山火。利用有向箭頭將這些因果關(guān)系進(jìn)行連接，得到靜態(tài)BN結(jié)構(gòu)如圖6所示。

圖6 輸電走廊山火風(fēng)險(xiǎn)靜態(tài)貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)Fig.6 Static Bayesian network structure of transmission line corridor wildfire risk

對(duì)靜態(tài)BN 在時(shí)間軸上進(jìn)行延拓，生成含有2個(gè)時(shí)間片的DBN，如圖7 所示。DBN 網(wǎng)絡(luò)中溫度、降水量、風(fēng)速、濕度和連續(xù)無雨日為動(dòng)態(tài)變量，因此利用有向邊將其相鄰時(shí)間片節(jié)點(diǎn)連接，以反映變量的時(shí)間變化，其余靜態(tài)變量則無需在時(shí)間片間進(jìn)行連接。

圖7 輸電走廊山火風(fēng)險(xiǎn)動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)Fig.7 The dynamic Bayesian network structure of transmission line corridor wildfire risk

4.3.2 建模結(jié)果分析

由于DBN 處理離散化的數(shù)據(jù)效率更高，且模型結(jié)果更加穩(wěn)定［34］，因此對(duì)連續(xù)型數(shù)據(jù)離散化處理。為體現(xiàn)有火與無火網(wǎng)格樣本分級(jí)離散化后的差異，逐一統(tǒng)計(jì)各因子的頻數(shù)分布，根據(jù)數(shù)據(jù)分布差異性進(jìn)行離散化處理，如表3所示。

為保障模型結(jié)果的穩(wěn)定性和真實(shí)性，采用5 折交叉驗(yàn)證法分別對(duì)BN 和DBN 進(jìn)行訓(xùn)練和測(cè)試，即通過分層采樣的方式將數(shù)據(jù)集劃分為5 個(gè)同等大小的互斥子集，每次以其中4 個(gè)子集作為訓(xùn)練集，余下的1 個(gè)子集作為測(cè)試集。利用極大似然估計(jì)法學(xué)習(xí)BN 和DBN 的節(jié)點(diǎn)參數(shù)，然后對(duì)測(cè)試集樣本進(jìn)行山火概率評(píng)估計(jì)算，將預(yù)測(cè)概率大于等于0.5 的樣本判識(shí)為有火，而小于0.5 的樣本判識(shí)為無火。將5 次測(cè)試集預(yù)測(cè)結(jié)果取平均值，使用ROC 曲線、AUC 值、預(yù)測(cè)準(zhǔn)確率和召回率對(duì)兩個(gè)模型進(jìn)行評(píng)價(jià)對(duì)比，結(jié)果如表4和圖8所示。

表4 模型評(píng)價(jià)指標(biāo)對(duì)比Tab.4 Comparison of Model Evaluation Indexes

圖8 DBN和BN模型的ROC曲線Fig.8 ROC curves of DBN and BN model

從表中可看出，DBN模型的準(zhǔn)確率為78.47%，高于BN 的70.99%，說明DBN 對(duì)于樣本集的擬合度更高。DBN 召回率為65.89%，優(yōu)于BN 的49.18%，表明DBN 對(duì)于有火網(wǎng)格樣本的預(yù)測(cè)更準(zhǔn)確，能有效降低有火的誤判或漏判。從圖8 可見，DBN 的ROC 曲線處于BN 的左上方，其AUC 值更高，說明DBN具有更好的分類效果。

綜上所述，在月時(shí)間尺度下，基于相鄰2 月的樣本構(gòu)建的DBN 山火風(fēng)險(xiǎn)評(píng)估模型明顯優(yōu)于靜態(tài)BN 模型，這意味著引入上一個(gè)階段的因素后有助于貝葉斯網(wǎng)絡(luò)進(jìn)行山火風(fēng)險(xiǎn)評(píng)估，并降低對(duì)有火的誤判或漏判率。

4.4 不同時(shí)間尺度下的山火預(yù)測(cè)模型對(duì)比分析

由4.3 節(jié)可知，動(dòng)態(tài)貝葉斯模型山火風(fēng)險(xiǎn)評(píng)估效果優(yōu)于靜態(tài)模型。但是在月時(shí)間尺度下，DBN的準(zhǔn)確率和召回率僅達(dá)78.47%和65.89%，并非足夠理想。為此，本文進(jìn)一步細(xì)化時(shí)間尺度，通過重新處理動(dòng)態(tài)特征數(shù)據(jù)，研究DBN 在15 d、10 d、7 d、3 d 等不同時(shí)間尺度下的山火風(fēng)險(xiǎn)評(píng)估效果。

采用與月時(shí)間尺度相同的方法，建立不同時(shí)間尺度下的動(dòng)態(tài)數(shù)據(jù)集，分別構(gòu)建DBN，利用5 折交叉驗(yàn)證法進(jìn)行訓(xùn)練和測(cè)試，結(jié)果如表5 和圖9 所示。隨著時(shí)間尺度的減小，模型準(zhǔn)確率和召回率不斷上升。其中3 d 模型結(jié)果最優(yōu)，準(zhǔn)確率達(dá)86.39%，召回率達(dá)80.33%。ROC 曲線隨著時(shí)間減小逐漸靠近左上方，對(duì)應(yīng)的AUC 值也逐漸增大，表明模型分類能力不斷提升。

表5 不同時(shí)間尺度下DBN模型評(píng)價(jià)指標(biāo)對(duì)比Tab.5 Comparison of evaluation indexes of DBN model under different time scales

圖9 不同時(shí)間尺度下DBN模型ROC曲線Fig.9 ROC curves of DBN model under different time scales

圖10 2022年清明節(jié)期間廣東省山火風(fēng)險(xiǎn)分布Fig.10 Distribution of wildire risk in Guangdong Province during Qingming Festival in 2022

上述結(jié)果表明，山火預(yù)測(cè)準(zhǔn)確率隨著時(shí)間尺度減小而提高，其原因主要為：根據(jù)山火因子重要度排序，氣象因子為重要因子，一方面，隨著時(shí)間尺度變小，溫度、降水量等氣象因子數(shù)據(jù)越來越精確，更加能夠反映當(dāng)時(shí)山火的發(fā)生容易程度；另一方面，DBN 將上下2 個(gè)時(shí)間片的動(dòng)態(tài)因子聯(lián)系起來，融入了時(shí)間信息，使得模型的擬合度不斷提升，進(jìn)而提高模型預(yù)測(cè)準(zhǔn)確度。

4.5 模型工程化應(yīng)用

為實(shí)現(xiàn)山火風(fēng)險(xiǎn)空間分布評(píng)估的工程化應(yīng)用，收集了2022年廣東省清明節(jié)期間（4月1日—4月6日）山火因子數(shù)據(jù)，采用訓(xùn)練好的3天DBN模型計(jì)算4月4日至4 月6 日山火發(fā)生概率。并分別以25%、50%、75%為間斷點(diǎn)，將山火發(fā)生概率劃分為低風(fēng)險(xiǎn)，中等風(fēng)險(xiǎn)，較高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)四級(jí)風(fēng)險(xiǎn)等級(jí)，利用地理信息軟件將山火風(fēng)險(xiǎn)可視化，繪制成1 km×1 km分辨率的廣東省山火風(fēng)險(xiǎn)分布圖，如圖11所示。

從分布圖可以發(fā)現(xiàn)，廣東省高風(fēng)險(xiǎn)區(qū)域主要集中在北部和東部地區(qū)，應(yīng)當(dāng)重點(diǎn)關(guān)注。根據(jù)當(dāng)?shù)仉娏咎峁┑膹V東省2022 年4 月4 日至4 月6 日期間共計(jì)84 起輸電走廊山火經(jīng)緯度信息，提取火點(diǎn)位置對(duì)應(yīng)的風(fēng)險(xiǎn)等級(jí)。經(jīng)驗(yàn)證，其中68 個(gè)火點(diǎn)落在較高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)區(qū)域，預(yù)測(cè)準(zhǔn)確率達(dá)80.95%，驗(yàn)證了分布圖的合理性和適用性。

5 結(jié)論

本文考慮時(shí)間因素對(duì)山火發(fā)生的影響，提出了一種基于動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)的輸電走廊山火風(fēng)險(xiǎn)時(shí)空評(píng)估方法。

1）根據(jù)輸電走廊山火月度分布不均衡特性，分別建立基于靜態(tài)貝葉斯網(wǎng)絡(luò)和動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)的山火風(fēng)險(xiǎn)評(píng)估模型，發(fā)現(xiàn)引入上一個(gè)時(shí)間的因子的動(dòng)態(tài)模型提高了評(píng)估準(zhǔn)確性。

2）在月時(shí)間尺度基礎(chǔ)上，研究了時(shí)間尺度對(duì)DBN 模型效果的影響，分別建立了15 d、10 d、7 d、3 d 動(dòng)態(tài)貝葉斯網(wǎng)絡(luò)山火風(fēng)險(xiǎn)評(píng)估模型。隨著時(shí)間尺度變小，模型評(píng)估效果逐步提升，其中以3 d為時(shí)間尺度的DBN 模型測(cè)試集準(zhǔn)確率達(dá)86.39%，召回率達(dá)80.33%，這是因?yàn)殡S時(shí)間細(xì)化，氣象因子數(shù)據(jù)越來越精確，提高了模型預(yù)測(cè)準(zhǔn)確性。

3）利用3 d 為時(shí)間尺度的DBN 模型評(píng)估了2022年清明節(jié)期間廣東省山火風(fēng)險(xiǎn)分布，繪制了1 km×1 km 分辨率山火風(fēng)險(xiǎn)分布圖，80.95%真實(shí)火點(diǎn)被預(yù)測(cè)成功，驗(yàn)證了本文所提方法的有效性。