陸秋琴，王 勃，黃光球

(西安建筑科技大學(xué) 管理學(xué)院，陜西 西安 710055)

0 引言

燃?xì)庖呀?jīng)處于城市能源組成部分的核心地位，在IEA成員國(guó)中，有16個(gè)國(guó)家的城市燃?xì)?即天然氣、液化氣)消耗量超過一次能源消耗總量的20%，其中7個(gè)國(guó)家超過30%[1]。我國(guó)現(xiàn)在正處于城市化、工業(yè)化的高速發(fā)展階段，城市的燃?xì)馄占奥室呀?jīng)高達(dá)98.42%，每年的燃?xì)夤?yīng)總量約2 803億m3。隨著近年來(lái)城市燃?xì)夤┬枞笨谌找鏀U(kuò)大、極端氣候、突發(fā)事故等方面的影響，導(dǎo)致城市的“氣荒”情況越來(lái)越嚴(yán)重，可能造成工業(yè)生產(chǎn)和電力供應(yīng)停滯，社會(huì)經(jīng)濟(jì)活動(dòng)和日常生活難以正常運(yùn)行；同時(shí)也可能引起社會(huì)的恐慌，短時(shí)間內(nèi)激發(fā)消費(fèi)需求，引起消費(fèi)劇增、價(jià)格暴漲，對(duì)社會(huì)穩(wěn)定造成嚴(yán)重的影響[2]。因此，針對(duì)城市愈演愈烈的“氣荒”狀況，通過對(duì)城市燃?xì)舛倘睉?yīng)急體系的評(píng)估，找出影響城市燃?xì)舛倘睉?yīng)急能力的主要因素，提出科學(xué)、有效的應(yīng)對(duì)策略，對(duì)提高城市突發(fā)性燃?xì)舛倘钡膽?yīng)急能力具有一定的現(xiàn)實(shí)意義，同時(shí)也有助于減少由城市突發(fā)性燃?xì)舛倘笔鹿室鸬母黝惤?jīng)濟(jì)損失。

針對(duì)以上不足，本文提出一種基于模糊Petri網(wǎng)的城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估模型，該模型基于應(yīng)急能力影響因素間的關(guān)聯(lián)過程，評(píng)價(jià)城市突發(fā)性燃?xì)舛倘睉?yīng)急能力。該評(píng)估模型具有以下特點(diǎn)：

1)該模型能夠結(jié)合評(píng)估體系中一些模糊或不確定的影響因素，根據(jù)層次化的模擬仿真對(duì)影響因素眾多且因素間信息反饋不明確的復(fù)雜城市突發(fā)性燃?xì)舛倘睉?yīng)急能力進(jìn)行評(píng)估。

2)該評(píng)估模型系統(tǒng)全面地從應(yīng)急的多個(gè)方面來(lái)構(gòu)建應(yīng)急能力影響因素體系，同時(shí)可以從不同層次間動(dòng)態(tài)的應(yīng)急反饋過程，形象、直觀地顯示城市突發(fā)性應(yīng)急過程中的邏輯因果關(guān)系與機(jī)理，圖形化、科學(xué)化地對(duì)應(yīng)急過程進(jìn)行仿真分析。

3)該模型的仿真過程通過采用極大代數(shù)運(yùn)算來(lái)提高整體模型運(yùn)算效率，同時(shí)采用AHP來(lái)減少應(yīng)急影響因素的評(píng)估主觀性，可以更加合理、有效地對(duì)城市突發(fā)性燃?xì)舛倘睉?yīng)急能力進(jìn)行評(píng)估。

1 應(yīng)急能力評(píng)估指標(biāo)體系

1.1 指標(biāo)體系構(gòu)建準(zhǔn)則及步驟

城市突發(fā)性燃?xì)舛倘睉?yīng)急能力是指從城市整體應(yīng)急系統(tǒng)的角度出發(fā)，通過政府部門的指揮決策、法規(guī)制度、應(yīng)急預(yù)案以及企業(yè)、人員、設(shè)備積極響應(yīng)等各方面因素參與下，對(duì)城市突發(fā)性燃?xì)舛倘睉?yīng)急全過程的緊急應(yīng)對(duì)與事后恢復(fù)能力。在城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估過程中，評(píng)估指標(biāo)體系的構(gòu)建是最關(guān)鍵的一個(gè)環(huán)節(jié)，它將直接影響應(yīng)急能力評(píng)估結(jié)果的科學(xué)性、全面性、可比性及有效性，但是影響城市突發(fā)性燃?xì)舛倘睉?yīng)急能力的各類因素中存在大量的主觀性與不確定性、耦合性的問題，因此必須用合適的評(píng)估準(zhǔn)則進(jìn)行深入地剖析來(lái)確定應(yīng)急能力評(píng)估體系。本文從應(yīng)急能力評(píng)估的目標(biāo)出發(fā)，綜合考慮了法律制度、人員、設(shè)備、管理、指揮決策等因素，在篩選具體指標(biāo)時(shí)以系統(tǒng)性、科學(xué)性、可行性為指導(dǎo)原則，遵循了(P)計(jì)劃-(D)實(shí)施-(C)檢驗(yàn)-(A)改進(jìn)的構(gòu)建程序[14]，建立能夠綜合反映城市突發(fā)性燃?xì)舛倘睉?yīng)急能力狀況的應(yīng)急指標(biāo)體系。構(gòu)建步驟如圖1所示。

圖1 城市突發(fā)性燃?xì)舛倘睉?yīng)急評(píng)估指標(biāo)體系構(gòu)建步驟Fig.1 Construction steps of index system for emergency assessment of urban sudden gas shortage

1.2 評(píng)估指標(biāo)體系的構(gòu)建

城市突發(fā)性燃?xì)舛倘睉?yīng)急是一個(gè)多階段且復(fù)雜的反饋綜合應(yīng)急網(wǎng)絡(luò)體系，由監(jiān)測(cè)預(yù)警體系、資源保障體系、應(yīng)急響應(yīng)體系、應(yīng)急組織體系、后期恢復(fù)評(píng)估體系5個(gè)子體系共同構(gòu)成，各個(gè)子體系在城市突發(fā)性燃?xì)舛倘笔录l(fā)生過程中，共同對(duì)城市突發(fā)性燃?xì)舛倘睉?yīng)急能力發(fā)揮著至關(guān)重要的作用。通過對(duì)應(yīng)急系統(tǒng)影響因素內(nèi)容以及文獻(xiàn)[12]分析，建立包括5個(gè)二級(jí)指標(biāo)和26個(gè)三級(jí)指標(biāo)的城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估體系，如表1所示。

表1 城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估指標(biāo)體系Table 1 Index system for emergency capability assessment of urban sudden gas shortage

1.3 應(yīng)急能力等級(jí)評(píng)判標(biāo)準(zhǔn)

表2 應(yīng)急能力等級(jí)評(píng)價(jià)標(biāo)準(zhǔn)Table 2 Assessment criteria of emergency capability grade

2 FPN評(píng)估模型及推理算法

2.1 模糊Petri網(wǎng)模型的定義

模糊Petri網(wǎng)是根據(jù)傳統(tǒng)Petri網(wǎng)不能處理評(píng)價(jià)指標(biāo)的不確定性和模糊性的特征擴(kuò)展而來(lái)，通過模糊集理論形成模糊Petri網(wǎng)評(píng)價(jià)體系來(lái)描述具有模糊行為的并發(fā)系統(tǒng)，因此被廣泛應(yīng)用于各個(gè)復(fù)雜系統(tǒng)的模擬仿真等實(shí)際問題。城市突發(fā)性燃?xì)舛倘睉?yīng)急是一個(gè)復(fù)雜動(dòng)態(tài)過程，模糊Petri網(wǎng)可以很好地展現(xiàn)出應(yīng)急行為變遷與反饋過程，還可以通過推理分析出影響因素之間的因果關(guān)系。

為了更明確地表達(dá)城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)價(jià)的模糊Petri網(wǎng)(USGSFPN)，將其定義為1個(gè)8元組，如式(1)所示。

USGSFPN=(P,T,D,IN,OUT,V,λ,M)

(1)

2.2 產(chǎn)生式規(guī)則的模糊FPN表示

根據(jù)USGSFPN的定義，每1個(gè)變遷都可以被看成是1個(gè)運(yùn)行規(guī)則，那么變遷的輸入庫(kù)所(d1,d2,d3,,dn)和輸出庫(kù)所(dg)分別表示運(yùn)行規(guī)則的前提條件和結(jié)論命題，可以通過邏輯聯(lián)接符號(hào)“and”或“or”進(jìn)行因果表示。通過對(duì)城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估體系分析，認(rèn)為應(yīng)急系統(tǒng)中的某1個(gè)影響因素可能就會(huì)導(dǎo)致整體應(yīng)急能力下降。即IFd1(w1) ord2(w2)or,,ordn(wn)，THENdg(wg1)，變遷發(fā)生規(guī)則以及觸發(fā)規(guī)則如圖2所示。變遷的可信度數(shù)學(xué)表達(dá)式如式(2)所示[14]：

wg1=max(w1μ1，w2μ2,,wnμn)

(2)

圖2 USGSPFN模型推理規(guī)則Fig.2 Inference rules of USGSPFN model

2.3 模糊矩陣推理算法及步驟

1)模糊推理算法

在USGSFPN模糊推理之前，為了清晰、準(zhǔn)確、簡(jiǎn)潔地表達(dá)推理矩陣運(yùn)算，推理過程中引入以下算子：設(shè)A,B,C均為n×m維矩陣，D為n維向量[16]，則有：

①加法算子⊕:A⊕B=C,則有cij=max(aij,bij),(i=1,2,,n;j=1,2,,m)

②乘法算子?:D?B=C,則有cij=di×bij,(i=1,2，,n;j=1,2,,m)

2)USGSFPN模糊推理步驟

根據(jù)相關(guān)算子規(guī)則進(jìn)行模糊推理運(yùn)算。首先將城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估體系轉(zhuǎn)化為USGSFPN網(wǎng)模型，根據(jù)模糊Petri網(wǎng)的USGSFPN定義分別確定應(yīng)急能力影響因素集IN與應(yīng)急行為集OUT，其中|P|=n，|T|=m，相應(yīng)矩陣形式分別為IN={αij}，OUT={γji}；然后根據(jù)歷史數(shù)據(jù)與應(yīng)急專家知識(shí)庫(kù)數(shù)據(jù)評(píng)估標(biāo)準(zhǔn)得到初狀態(tài)M(0)；根據(jù)層次分析法可得到應(yīng)急系統(tǒng)中各個(gè)應(yīng)急影響因素的權(quán)重值ωi，同時(shí)得到模糊變遷的置信度向量U={Ui}={ωi}。具體模糊推理計(jì)算步驟如下所示，其中，輸入包括：初始化輸入矩陣IN；輸出矩陣OUT；置信度向量U；應(yīng)急能力評(píng)判矩陣Q和初始狀態(tài)矩陣M(0)。輸出包括：各個(gè)庫(kù)所應(yīng)急能力評(píng)估值。

①步驟一：令初始狀態(tài)k=0，變遷發(fā)生時(shí)，初始應(yīng)急能力影響因素輸入庫(kù)所的可信度INTM(K)；變遷觸發(fā)后，應(yīng)急行為輸出庫(kù)所的可信度U?OUT。

②步驟二：計(jì)算變遷觸發(fā)生后的下一個(gè)狀態(tài)，應(yīng)急能力影響因素行為變更的可信度為M(k+1)=M(k)⊕(U?OUT)[INTM(k)]。

③步驟三：若M(k+1)≠M(fèi)(k)，則令k=k+1，轉(zhuǎn)到步驟二；若M(k+1)=M(k)，則模糊推理結(jié)束。

④步驟四：計(jì)算各個(gè)庫(kù)所應(yīng)急能力等級(jí)F=M(k)QT，F(xiàn)為n×1階矩陣，表示n個(gè)庫(kù)所的應(yīng)急能力評(píng)估等級(jí)值,該矩陣的最后1個(gè)元素對(duì)應(yīng)的等級(jí)標(biāo)準(zhǔn)即為城市突發(fā)性燃?xì)舛倘睉?yīng)急能力指標(biāo)綜合等級(jí)。

3 實(shí)例仿真

3.1 構(gòu)建FPN模型

為了驗(yàn)證上述應(yīng)急能力模型和模糊推理算法的可行性和科學(xué)性，本文以西安市突發(fā)性燃?xì)舛倘睉?yīng)急為例進(jìn)行模擬仿真分析。通過基本的模糊規(guī)則邏輯關(guān)系將城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估指標(biāo)體系(表1)轉(zhuǎn)化為FPN模型，如圖3所示。

圖3 城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)估體系FPN模型Fig.3 FPN model of emergency capability assessment system for urban sudden gas shortage

3.2 FPN應(yīng)急能力評(píng)估參數(shù)求解

根據(jù)USGSFPN的定義，可得到城市突發(fā)性燃?xì)舛倘睉?yīng)急FPN模型的輸入矩陣IN和輸出矩陣OUT，通過查閱《西安市統(tǒng)計(jì)年鑒2018》、相關(guān)歷史數(shù)據(jù)以及城市應(yīng)急專家知識(shí)庫(kù)等,以及分析歷史數(shù)據(jù)和專家評(píng)判法給出的庫(kù)所評(píng)估向量取值量化規(guī)則可得到初始評(píng)估向量矩陣M(0)，如式(3)所示。

(3)

IN={αij},α1,1=1,α2,2=1,α3,3=1,α4,4=1,α5,5=1,α6,6=1,α7,7=1,α8,8=1,α9,9=1,α10,10=1,α11,11=1,α12,12=1,α13,13=1,α14,14=1,α15,15=1,α16,16=1,α17,17=1,α18,18=1,α19,19=1,α20,20=1,α21,21=1,α22,22=1,α23,23=1,α24,24=1,α25,25=1,α26,26=1,α27,27=1,α28,28=1,α29,29=1,α30,30=1,α31,31=1,其余的輸入值αij=0,(i=1,2,,n;j=1,2,,m)。

OUT={γij},γ1,27=1,γ2,27=1,γ3,27=1,γ4,27=1,γ5,27=1,γ6,28=1,γ7,28=1,γ8,28=1,γ9,28=1,γ10,28=1,γ11,29=1,γ12,29=1,γ13,29=1,γ14,29=1,γ15,29=1,γ16,29=1,γ17,29=1,γ18,30=1,γ19,30=1,γ20,30=1,γ21,30=1,γ22,30=1,γ23,31=1,γ24,31=1,γ25,31=1,γ26,31=1,γ27,32=1,γ28,32=1,γ29,32=1,γ30,32=1,γ31,32=1,其余的輸出值γij=0,(i=1,2,,n;j=1,2,,m)。

3.3 影響因素的權(quán)重值計(jì)算

計(jì)算各影響因素的權(quán)重值(即變遷的置信度向量U)，根據(jù)城市突發(fā)性燃?xì)舛倘睉?yīng)急評(píng)估體系，以二級(jí)指標(biāo)為例，通過層次分析法從應(yīng)急預(yù)警(B1)、資源保障(B2)、應(yīng)急響應(yīng)(B3)、應(yīng)急組織(B4)、后期恢復(fù)評(píng)估(B5)等5個(gè)方面確定指標(biāo)權(quán)重，根據(jù)指標(biāo)重要程度采用1～9標(biāo)度法得到判斷矩陣A-Bi。經(jīng)歸一化處理后的權(quán)重值如表3所示。判斷矩陣A-Bi的λmax=5.353 239，CI=0.088 31,RI=1.12，進(jìn)而求得CR=0.078 8<0.1,滿足一致性檢驗(yàn)。

表3 A-Bi層判斷矩陣及指標(biāo)權(quán)重Table 3 A-Bi layer judgment matrix and index weights

同理可得其他各變遷的置信度向量為：U={0.30,0.43,0.14,0.08,0.05,0.46,0.18,0.23,0.09,0.04,0.49,0.14,0.05,0.16,0.16,0.28,0.36,0.05,0.10,0.11,0.07,0.03,0.08,0.48,0.20,0.19,0.04,0.53,0.27,0.15,0.05}。

3.4 模糊推理運(yùn)算

根據(jù)USGSFPN模糊推理步驟，利用MATLAB軟件進(jìn)行模糊運(yùn)算，由推理算法程序進(jìn)行迭代得到最終的M(k)為：

M(2)=M(3)=[

根據(jù)最后一個(gè)元素為5.424,表明了西安市突發(fā)性燃?xì)舛倘睉?yīng)急能力等級(jí)為中等；二級(jí)指標(biāo)中，監(jiān)測(cè)預(yù)警能力為6.078，表明監(jiān)測(cè)預(yù)警方面能力等級(jí)較高;資源保障能力為6.204，表明資源保障方面應(yīng)急能力等級(jí)較高;應(yīng)急響應(yīng)能力為4.974，表明應(yīng)急響應(yīng)能力等級(jí)中等;應(yīng)急組織能力為5.464，表明應(yīng)急組織能力等級(jí)中等;后期處理能力5.572，表明后期恢復(fù)評(píng)估能力等級(jí)中等。三級(jí)指標(biāo)應(yīng)急能力等級(jí)限于篇幅不詳細(xì)列出。

3.5 分析及建議

在USGSFPN評(píng)估模型中，由于上層指標(biāo)支配著下層指標(biāo)，同時(shí)下層指標(biāo)共同制約著上層指標(biāo)，因此使用下層指標(biāo)評(píng)估等級(jí)均值來(lái)作為下層指標(biāo)對(duì)上層指標(biāo)的制約相關(guān)性閥值，若下層指標(biāo)評(píng)估等級(jí)值低于閥值，則對(duì)上層指標(biāo)的制約相關(guān)性越強(qiáng)。根據(jù)上述方法和仿真計(jì)算結(jié)果顯示，西安市突發(fā)性燃?xì)舛倘睉?yīng)急能力水平為中等，制約綜合應(yīng)急能力的二級(jí)指標(biāo)主要有應(yīng)急響應(yīng)能力、應(yīng)急組織能力、后期恢復(fù)評(píng)估能力(閥值=5.658 4)。同理，通過數(shù)值分析可得對(duì)上述3個(gè)二級(jí)指標(biāo)制約相關(guān)性較強(qiáng)的三級(jí)指標(biāo)，并根據(jù)其相應(yīng)的不足對(duì)提升西安市突發(fā)性燃?xì)舛倘睉?yīng)急能力提出建議。

1)在應(yīng)急響應(yīng)能力方面，應(yīng)建立應(yīng)急專項(xiàng)領(lǐng)導(dǎo)小組加強(qiáng)城市突發(fā)性事件應(yīng)急指揮決策反應(yīng)能力，總結(jié)、積累應(yīng)急事件經(jīng)驗(yàn)，不斷完善突發(fā)性事件應(yīng)急預(yù)案的完整性，建立應(yīng)急響應(yīng)的預(yù)演機(jī)制、加強(qiáng)宣傳力度，使得政府和民眾應(yīng)急響應(yīng)能力進(jìn)一步提升。

2)在應(yīng)急組織能力方面，精簡(jiǎn)應(yīng)急單位、部門的數(shù)量并且做到明確應(yīng)急組織的責(zé)任劃分，建立城市區(qū)域間的聯(lián)動(dòng)協(xié)調(diào)機(jī)制，不斷提升應(yīng)急組織效率的同時(shí)，也提高應(yīng)急組織能力的快速化和實(shí)效性。

3)在后期恢復(fù)評(píng)估能力方面，需要加強(qiáng)對(duì)應(yīng)急專業(yè)隊(duì)伍的建設(shè)和后期恢復(fù)資金的投入，來(lái)提高事后評(píng)估的準(zhǔn)確性、減少恢復(fù)的時(shí)間，加強(qiáng)突發(fā)性事件后期恢復(fù)評(píng)估能力的建設(shè)也是為應(yīng)急能力的不斷提升打下堅(jiān)實(shí)的基礎(chǔ)。

4 結(jié)論

1)本文提出的模糊Petri網(wǎng)方法綜合考慮了系統(tǒng)中的各種影響因素，通過模糊推理和算法得到應(yīng)急能力等級(jí)評(píng)價(jià)指標(biāo)，在某種程度上避免了影響因素分析的不確定性和主觀性問題，使應(yīng)急能力評(píng)估結(jié)果更加準(zhǔn)確、有效。

2)通過對(duì)模糊Petri網(wǎng)模型及其推理算法在城市突發(fā)性燃?xì)舛倘睉?yīng)急能力評(píng)價(jià)中的應(yīng)用，表明該方法不僅能夠評(píng)估城市應(yīng)急系統(tǒng)綜合的應(yīng)急能力等級(jí)，而且還能評(píng)估各個(gè)子系統(tǒng)的影響因素以及中間層的應(yīng)急能力等級(jí)，具有可行性、科學(xué)性和準(zhǔn)確性的特點(diǎn)。

3)以西安市突發(fā)性燃?xì)舛倘睉?yīng)急為例進(jìn)行模擬仿真分析，最終根據(jù)模擬仿真結(jié)果，在應(yīng)急響應(yīng)能力、應(yīng)急組織能力、后期恢復(fù)評(píng)估能力等3個(gè)方面提出措施建議。