王倩倩，孟繁宇，2，曾益萍，張少標，吳國華，楊麗麗

(1.南方科技大學(xué) 統(tǒng)計與數(shù)據(jù)科學(xué)系,廣東 深圳 518055； 2.南方科技大學(xué) 前沿與交叉科學(xué)研究院,廣東 深圳 518055;3.深圳市城市公共安全技術(shù)研究院,廣東 深圳 518055)

0 引言

城市火災(zāi)風(fēng)險評估一直是政府和民眾關(guān)心的問題，如何對火災(zāi)事故風(fēng)險進行分析并量化，是目前我國相關(guān)學(xué)者正在研究的重要課題。國內(nèi)外學(xué)者對城市區(qū)域火災(zāi)風(fēng)險評估進行一系列的研究：張一先等[1]采用指數(shù)法對蘇州古城區(qū)的火災(zāi)危險性進行分級；魯鈺雯等[2]通過ArcGIS軟件對廈門市火災(zāi)風(fēng)險進行評估；馬礪等[3]采用主成分分析法對火災(zāi)成因進行分析和量化；李犇[4]在辨識研究區(qū)域火災(zāi)危險源的基礎(chǔ)上，篩選出相應(yīng)的火災(zāi)風(fēng)險評估指標，進行火災(zāi)風(fēng)險評估；賈進章等[5]采用突變理論與模糊數(shù)學(xué)相結(jié)合的評價方法求得火災(zāi)危險性評價結(jié)果；劉紀華[6]提出基于FTA，ETA的綜合分析；王元鵬等[7]引用“社群”相關(guān)概念和算法，建立火災(zāi)社區(qū)之間拓撲關(guān)系的分析模型。 文獻[1-7]雖然可以快速地對火災(zāi)風(fēng)險進行評估，但存在過度依賴專業(yè)的判斷和過去的經(jīng)驗的問題，基于數(shù)據(jù)運用聚類算法建立風(fēng)險評估模型的研究較為少見。因此，本文在前人研究的基礎(chǔ)上，基于數(shù)據(jù)思想運用DBSCAN聚類算法，并結(jié)合層次分析法構(gòu)建火災(zāi)風(fēng)險評估的驗證模型，同時加入火災(zāi)原因因子通過線性變換進行量化，綜合得出風(fēng)險評估結(jié)果。

1 研究方法及數(shù)據(jù)來源

1.1 研究方法

1.1.1 DBSCAN聚類算法

DBSCAN(Density-Based Spatial Clustering of Applications with Noise)是1種經(jīng)典的基于密度空間的聚類方法[8]，其將簇定義為密度相連的點的最大集合，能夠把足夠密度的區(qū)域劃分為簇，并可以在含有噪點的空間中發(fā)現(xiàn)任意形狀的簇。

本文使用該算法識別深圳市鹽田區(qū)火災(zāi)事故的空間聚集特征，通過多次改變參數(shù)值進行多次迭代最終得到最優(yōu)聚類效果。在運用DBSCAN算法進行聚類時，需要計算2點之間的距離，使用haversine方法計算2個經(jīng)緯度之間的距離，最后運用DBSCAN算法進行聚類。其中haversine計算公式如式(1)所示：

(1)

式中：haversine為半正矢函數(shù)；d為所求距離，km；R為地球半徑，km；φ1,φ2為2點的緯度，(°)；Δλ為2點經(jīng)度的差值，(°)。

1.1.2 層次分析法

層次分析法[9,14]將決策過程中的信息采用定性與定量相結(jié)合的方式進行優(yōu)化處理。

AHP具體步驟如下：

1)建立系統(tǒng)的層次結(jié)構(gòu)模型

按照因素間的相互關(guān)系影響以及隸屬關(guān)系，將因素聚集為不同層次，本文將建立3層結(jié)構(gòu)模型。DBSCAN算法對火災(zāi)事故的地理位置聚類產(chǎn)生的結(jié)果集內(nèi)的中心點作為第1層數(shù)據(jù)；火災(zāi)事故場景構(gòu)成層次結(jié)構(gòu)模型的第2層數(shù)據(jù)；火災(zāi)事故原因作為層次結(jié)構(gòu)模型的第3層數(shù)據(jù)。

2)構(gòu)造判斷矩陣

層次分析結(jié)構(gòu)模型確立后，問題將轉(zhuǎn)化為層次中各因素和上層因素相對重要性的排序問題。本文將采用1～9及其倒數(shù)的標度方法[9]對每1層的權(quán)重進行賦值，構(gòu)造判斷矩陣。若A層元素中ak與下1層中B1,B2,…,Bn有聯(lián)系，則判斷矩陣見表1。其中aii=1,aji=1/aij(i,j=1,2,…,n)。

表1 判斷矩陣Table 1 Judgment matrix

3)計算判斷矩陣的最大特征根和特征向量

②計算最大特征值λmax，如式(2)所示：

(2)

式中：λmax為最大特征值；(BW)i為BW的第i個分量；BW為特征向量。

1.1.3 區(qū)域火災(zāi)綜合風(fēng)險值計算

1)線性變換

為減小所得數(shù)據(jù)的差異敏感度，對數(shù)據(jù)進行線性變換，此處采用底數(shù)為2的對數(shù)對其轉(zhuǎn)化，計算第3層風(fēng)險值。線性變換及風(fēng)險值計算公式分別如式(3)～(4)所示：

(3)

(4)

式中：ri為所需要轉(zhuǎn)化的數(shù)據(jù)ai通過對數(shù)轉(zhuǎn)化后的數(shù)值；a和b為線性變換上下限，a=50，b=100；RL3為第3層風(fēng)險值。

2)綜合風(fēng)險值計算

由1.1.2節(jié)計算得知分類后第1，2層的權(quán)重以及第3層的風(fēng)險值，對于聚類算法產(chǎn)生的每個區(qū)域的火災(zāi)風(fēng)險值，可通過每層的數(shù)據(jù)及其風(fēng)險權(quán)重計算得到。應(yīng)用線性加權(quán)方法確定火災(zāi)風(fēng)險的模型計算如式(5)～(7)所示：

(5)

(6)

(7)

1.2 數(shù)據(jù)來源

本文研究的案例為深圳市鹽田區(qū)2011—2019年的火災(zāi)事故數(shù)據(jù)，字段主要有火災(zāi)發(fā)生時間、地點、火災(zāi)發(fā)生原因。通過數(shù)據(jù)分析，探究在不同火災(zāi)事故場景下火災(zāi)發(fā)生起數(shù)隨時間、空間以及火災(zāi)類型之間的關(guān)系，確定所需要的模型變量；刪除冗雜數(shù)據(jù)，剩余火災(zāi)427起。

將火災(zāi)事故地點按照《關(guān)于城市消防管理工作的規(guī)定》[10-11]對消防安全重點單位實施分級管理和劃分，并結(jié)合對數(shù)據(jù)的認識及需要將場所的管轄范圍分為10類。賓館、飯館、旅店、招待所、咖啡館、酒吧等娛樂場所為第1類P1；室內(nèi)、室外農(nóng)貿(mào)市場、商場，三合一、多合一場所為第2類P2；候車(機、船)室、公共交通工具為第3類P3；體育場(館)、游泳場(館)、公園、醫(yī)院、博物館、美術(shù)館、動物園、圖書館、書店、報刊亭為第4類P4；住宅、職工宿舍、學(xué)生宿舍為第5類P5；辦公場所，公司為第6類P6；工廠、拆遷工地、建筑工地、建筑施工場所、倉庫為第7類P7；變壓器、露天框架為第8類P8；花卉、林木(含果園)、垃圾及廢棄物、苗圃、農(nóng)田、養(yǎng)殖場為第9類P9；其他、誤報為第10類P10。

2 深圳市鹽田區(qū)火災(zāi)風(fēng)險計算與分析

2.1 研究步驟

為篩選出合適變量進入模型，本文將深圳市鹽田區(qū)火災(zāi)事故數(shù)據(jù)進行清洗與分析，分析中發(fā)現(xiàn)火災(zāi)事故發(fā)生場所與眾多因素關(guān)聯(lián)性很強，因此分別探究在不同場景下火災(zāi)發(fā)生起數(shù)隨時間(年份、月份、小時、周末)、空間和火災(zāi)原因的分布特性，最終確定地理位置、火災(zāi)事故場所及事故原因3個指標來計算城市區(qū)域火災(zāi)風(fēng)險值。研究步驟如圖1所示。

2.2 區(qū)域火災(zāi)風(fēng)險的聚類模型

DBSCAN聚類算法

本文經(jīng)過數(shù)次迭代調(diào)換參數(shù)，最終確定Eps值為0.5，MmPts值為10，聚類結(jié)果為4類，分別標記為A，B，C，D類，聚類結(jié)果圖如圖2所示。

將每類數(shù)據(jù)匯總，計算標準偏差，觀察每1類內(nèi)部數(shù)據(jù)的波動性大小，具體見表2。

由圖2及表2可知，對算法聚類產(chǎn)生的每1類的數(shù)據(jù)計算標準偏差，值接近于0，說明每1簇內(nèi)部的點比較集中，波動較小，結(jié)果表明參數(shù)值選取的相對合理。

圖2 DBSCAN聚類結(jié)果Fig.2 DBSCAN clustering results

表2 DBSCAN算法聚類結(jié)果匯總Table 2 Summary of DBSCAN algorithm clustering results

2.3 層次分析模型構(gòu)建及驗證

2.3.1 分類方法的分層劃分

根據(jù)聚類結(jié)果構(gòu)建層次分析法的第1層數(shù)據(jù)，由于每類的數(shù)量反映火災(zāi)的發(fā)生起數(shù)，即數(shù)值越大則在此區(qū)域發(fā)生的火災(zāi)越多，反之則越少。因此根據(jù)每類的發(fā)生起數(shù)作為指標值的參考值構(gòu)造判斷矩陣，見表3。

表3 第1層判斷矩陣數(shù)據(jù)表示及計算Table 3 Data representation and calculation of first-level judgment matrix

第2層數(shù)據(jù)由火災(zāi)事故場景構(gòu)成。根據(jù)各類場所的性質(zhì)、功能類型及相關(guān)文獻[12-13]，得到每類場所的指標值排序：P6

基于第1層數(shù)據(jù)構(gòu)造第2層數(shù)據(jù)的判斷矩陣，每類包含的場所不同，因此構(gòu)造判斷矩陣不同，利用1.1.2節(jié)的計算方法，得到各項指標的權(quán)重、最大特征值及隨機一致性比例。由于篇幅有限，本文將以A區(qū)域為例進行展示每層的判斷矩陣計算值，見表4。

第3層數(shù)據(jù)是由火災(zāi)事故原因構(gòu)成，根據(jù)不同類型火災(zāi)發(fā)生起數(shù)[3]作為指標的參考值，由1.1.3中方法計算得到第3層風(fēng)險值，風(fēng)險矩陣見表5。

表4 第2層A類區(qū)域判斷矩陣數(shù)據(jù)表示及計算Table 4 Data representation and calculation of second-level judgment matrix for area A

表5 第3層A類區(qū)域風(fēng)險矩陣判斷矩陣數(shù)據(jù)表示及計算Table 5 Data representation and calculation of third-level judgment matrix for area A

2.3.2 區(qū)域火災(zāi)綜合風(fēng)險值計算

通過1.1.3方法對模型計算，得到每個區(qū)域的火災(zāi)風(fēng)險值，以區(qū)域A為例，見表6。

根據(jù)區(qū)域火災(zāi)實際防控設(shè)定量化范圍[50，100]及風(fēng)險分為4級[14-15]，見表7。

表6 區(qū)域A風(fēng)險值計算Table 6 Calculation of risk value for area A

表7 火災(zāi)風(fēng)險等級量化Table 7 Quantification of fire risk levels

2.4 結(jié)果分析

將各區(qū)域風(fēng)險計算結(jié)果匯總，見表8。4類區(qū)域地理分布如圖3所示。

表8 區(qū)域風(fēng)險值匯總Table 8 Summary of regional risk values

圖3 模型結(jié)果區(qū)域分布Fig.3 Regional distribution of model results

聚類區(qū)域結(jié)果討論：

1)A類區(qū)域中住宅類最多，占比39.3%，其次為垃圾廢棄物、花卉類，占比25.2%。起火原因類型中電氣火災(zāi)占比最高達33.1%，間接驗證了基于聚類并結(jié)合場所、原因?qū)馂?zāi)風(fēng)險計算得到的區(qū)域風(fēng)險值的合理性。因此，在針對此區(qū)域的消防安全中，應(yīng)加強社區(qū)安全消防的宣傳，定期組織疏散滅火演習(xí)，及時發(fā)現(xiàn)和整改火災(zāi)隱患，降低區(qū)域火災(zāi)風(fēng)險。

2)B類區(qū)域位于高速交叉處，地形復(fù)雜，為事故多發(fā)地，且周圍工業(yè)園區(qū)較多，一旦發(fā)生火災(zāi)，造成的影響將會很大。在B區(qū)域中交通類型火災(zāi)占比較大，為24.04%，垃圾廢棄物類占比為12.57%。由第3部分的計算得到B類區(qū)域風(fēng)險最高，達96.25。建議B區(qū)域?qū)S房、倉庫，要逐一簽訂防火安全責(zé)任書，嚴查嚴管，防患未然。

3)C類區(qū)域毗鄰梧桐山，有部分居民區(qū)，地點較為偏僻，在2011—2019年間共發(fā)生22起火災(zāi)，其中住宅類火災(zāi)為15起，火災(zāi)原因多為電氣火災(zāi)，為火災(zāi)事故低風(fēng)險區(qū)，建議及時更換老舊設(shè)施，排查隱患，提高用火安全意識。

4)D類區(qū)域位于大梅沙景區(qū)，周圍酒店居多，有部分小區(qū)和學(xué)校。其中住宅類和交通類火災(zāi)較多，合計達50%，火災(zāi)原因多為電氣火災(zāi)，占比38.9%。建議D類區(qū)域嚴格按照消防安全規(guī)定合理設(shè)置消防設(shè)施，提高酒店住宿的防火等級。

3 結(jié)論

1)基于DBSCAN聚類算法計算城市區(qū)域火災(zāi)風(fēng)險，聚類結(jié)果產(chǎn)生4類，由于各區(qū)域功能不同及地理位置不同，區(qū)域火災(zāi)風(fēng)險值有明顯差異。

2)通過結(jié)合層次分析法進行驗證，模型計算出的4個區(qū)域結(jié)果與實際火災(zāi)情況具有一致性，結(jié)果表明B區(qū)域的消防安全狀況相對于鹽田區(qū)其他區(qū)域風(fēng)險較高，C區(qū)域火災(zāi)風(fēng)險較低，消防安全系數(shù)較高。

3)本文區(qū)域火災(zāi)風(fēng)險值計算的嘗試，可為開展城市區(qū)域火災(zāi)風(fēng)險評估提供參考和借鑒。