潘科，關(guān)守安

(1.大連交通大學(xué) 土木與安全工程學(xué)院，遼寧 大連 116028； 2．遼寧省安全科學(xué)研究院，遼寧 沈陽(yáng) 110004)

0　引言

地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)在保障地鐵車站及地鐵列車的安全運(yùn)營(yíng)方面具有重要的作用.近年來(lái)，關(guān)于地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)體系及其方法的研究發(fā)展迅速.文獻(xiàn)[1]分析了《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》在權(quán)重確定及評(píng)價(jià)方法方面的不足，嘗試在地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)中引入變權(quán)理論和相對(duì)差異函數(shù)，以彌補(bǔ)個(gè)別指標(biāo)因權(quán)重過(guò)小導(dǎo)致的評(píng)價(jià)結(jié)果失真問(wèn)題.文獻(xiàn)[2]針對(duì)地鐵建設(shè)和運(yùn)營(yíng)特點(diǎn)建立了地鐵運(yùn)營(yíng)安全評(píng)價(jià)的多級(jí)可拓評(píng)價(jià)模型.文獻(xiàn)[3]為評(píng)估地鐵安全運(yùn)營(yíng)水平，建立了基于DT 法的地鐵設(shè)備設(shè)施動(dòng)態(tài)分級(jí)模型.文獻(xiàn)[4]基于Petri網(wǎng)絡(luò)建立了地鐵火災(zāi)應(yīng)急響應(yīng)模型.地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)體系的完備性、權(quán)重計(jì)算及評(píng)價(jià)方法的科學(xué)性是影響風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的主要因素，在現(xiàn)有的指標(biāo)權(quán)重確定方法中，專家經(jīng)驗(yàn)法和AHP法不同程度上且有一定的主觀性，本文基于地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)的兩大基本特征，即各指標(biāo)對(duì)于地鐵運(yùn)營(yíng)事故發(fā)生的可能性(RP)和后果的嚴(yán)重程度(RI)兩方面考慮其影響，并采用三角模糊理論來(lái)解決風(fēng)險(xiǎn)影響因素的不確定性問(wèn)題[5- 7]，從而確定各風(fēng)險(xiǎn)指標(biāo)的權(quán)重；考慮到可變模糊集理論解決模糊問(wèn)題的合理性[8]，論文提出了地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)的可變模糊評(píng)價(jià)方法，該方法能夠科學(xué)合理地確定影響地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)各評(píng)價(jià)指標(biāo)對(duì)各風(fēng)險(xiǎn)等級(jí)的相對(duì)隸屬度，并且可以通過(guò)變化模型及其參數(shù)，合理地確定出地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)等級(jí)，提高其結(jié)果的可信度.

1　地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)體系的構(gòu)建

基于系統(tǒng)安全的原理，影響系統(tǒng)安全的關(guān)鍵因素是“人-車(設(shè)備)-環(huán)境-管理”4個(gè)方面，對(duì)于地鐵運(yùn)營(yíng)，將“人”的因素歸納為運(yùn)營(yíng)組織及管理層，針對(duì)《城市軌道交通技術(shù)規(guī)范》(GB 50490-2009)在“運(yùn)營(yíng)”、“限界”、“土建工程”、“車輛”、“機(jī)電設(shè)備”等方面的要求，并結(jié)合《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》(GB/T 50438-2007)及文獻(xiàn)[4]的成果，建立并完善了地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系，如表1所示.

表1　地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系及取值

注：(1)表1中1級(jí)指標(biāo)、2級(jí)指標(biāo)的有關(guān)內(nèi)容主要參考《城市軌道交通技術(shù)規(guī)范》(GB 50490-2009)地鐵運(yùn)營(yíng)的有關(guān)要求，重點(diǎn)參考《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》(GB/T 50438-2007)的要求.(2)表1中的指標(biāo)分值為《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》附錄A-P中規(guī)定的該項(xiàng)總分值.考慮到計(jì)算的需要，將附錄規(guī)定的該項(xiàng)總分值與實(shí)際取值轉(zhuǎn)化為百分制，即第6列及第7列所示，如100(10)表示附錄中該項(xiàng)總分值為10分，轉(zhuǎn)變后為100分，85(8.5)表示該項(xiàng)實(shí)際取值為8.5分，將其轉(zhuǎn)變成百分制后為85分.(3)對(duì)于《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》附錄A-P中規(guī)定的有關(guān)定性定量指標(biāo)分值，除依據(jù)該條款外，還應(yīng)參照現(xiàn)行的地鐵運(yùn)營(yíng)規(guī)范，如對(duì)于附錄N土建評(píng)價(jià)表中“地下、高架結(jié)構(gòu)與車站建筑”，其N1項(xiàng)要考慮《城市軌道交通技術(shù)規(guī)范》7.3.18條關(guān)于地鐵換乘車站防火分區(qū)的有關(guān)要求，N11項(xiàng)要考慮進(jìn)站口的安檢系統(tǒng).

2　三角模糊理論及其權(quán)重確定方法的構(gòu)建

風(fēng)險(xiǎn)為事故發(fā)生的可能性及后果嚴(yán)重程度的綜合度量，而不同的指標(biāo)對(duì)于事故發(fā)生的可能性及后果嚴(yán)重程度是不一樣的，且有一定的不確定性，因此科學(xué)的權(quán)重確定方法是考慮風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)對(duì)事故發(fā)生后果及嚴(yán)重程度的影響，并采用數(shù)學(xué)方法來(lái)解決其賦值過(guò)程中的不確定問(wèn)題，而三角模糊理論可以有效解決這一問(wèn)題.

(1)

式中，在地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)中，bl為指標(biāo)賦值區(qū)間中的最小值;br為指標(biāo)賦值區(qū)間中的最大值;bq為指標(biāo)賦值區(qū)間中的最可能值.

(2)

式中，al，ar，aq的意義同bl，br，bq.

根據(jù)地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)影響因素及三角模糊理論，其權(quán)重確定步驟如下.

(1)定義風(fēng)險(xiǎn)因素(RF)

根據(jù)地鐵運(yùn)營(yíng)的有關(guān)特點(diǎn)及風(fēng)險(xiǎn)的兩大基本特征，即事故發(fā)生的可能性(RP)以及后果的嚴(yán)重程度(RI).依據(jù)文獻(xiàn)[5]及文獻(xiàn)[6]，表2、表3給出RP和RI模糊數(shù)的定義，包括其相應(yīng)的語(yǔ)言變量、地鐵運(yùn)營(yíng)中關(guān)于可能性及后果嚴(yán)重程度的參考值和對(duì)應(yīng)的三角模糊數(shù).

表2　事故發(fā)生的可能性(RP)描述

表3　后果的嚴(yán)重程度(RI)描述

注：事故運(yùn)營(yíng)等級(jí)分類參照《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》(GB/T 50438-2007)第19.2.1條.

(2)確定風(fēng)險(xiǎn)因素(RF)

按照風(fēng)險(xiǎn)的定義，得出風(fēng)險(xiǎn)層次結(jié)構(gòu)的底層指標(biāo)的風(fēng)險(xiǎn)因素RF的計(jì)算公式如下：

RF=RP?RI

(3)

(3)解模糊化過(guò)程

使用重心解模糊法[7]，其具體轉(zhuǎn)化公式如下：

(4)

(4)歸一化處理

1級(jí)指標(biāo)和2級(jí)指標(biāo)的權(quán)重ωi和ωij可由式(8)求得：

(5)

3　地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)的可變模糊評(píng)價(jià)模型

設(shè)地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)因素集為{y1,y2,…,yn}，y1,y2,…,yn為1級(jí)指標(biāo)，每個(gè)1級(jí)指標(biāo)有m個(gè)級(jí)指標(biāo)，則有待評(píng)價(jià)對(duì)象的特征值矩陣為X=(xij).xij的值主要由被評(píng)價(jià)對(duì)象的實(shí)際情況與《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》附錄A-P中規(guī)定的有關(guān)定性定量指標(biāo)分值取得，如對(duì)于安全管理評(píng)價(jià)中的“事故應(yīng)急救援體系(20分)”包括“應(yīng)急救援機(jī)構(gòu)(3)”、“預(yù)案制定情況(4分)”、“預(yù)案管理(1分)”、“應(yīng)急設(shè)備及人員配備情況及維護(hù)體系(6分)”、“事故應(yīng)急培訓(xùn)及應(yīng)急演練(4分)”、“當(dāng)年緊急事故處置(2分)”五項(xiàng)內(nèi)容，而每項(xiàng)內(nèi)容又根據(jù)其分項(xiàng)組成進(jìn)行打分，其分值結(jié)果相對(duì)符合實(shí)際情況.

3.1　相對(duì)差異函數(shù)模型及其指標(biāo)特征值的模糊化

3.1.1相對(duì)隸屬函數(shù)及相對(duì)差異函數(shù)

設(shè)U為論域，u為U的任意元素，元素u的一對(duì)對(duì)立模糊概念，或u對(duì)立的兩種基本模糊屬性：A和Ac.分別賦予A和Ac處于共維差異中介過(guò)渡的兩個(gè)端點(diǎn)以1、0與0、1的區(qū)間數(shù).分別在1到0和0到1的數(shù)軸上構(gòu)成一對(duì)[1，0]與[0，1]閉區(qū)間數(shù)的連續(xù)統(tǒng).對(duì)于U中的元素u，都在該連續(xù)統(tǒng)的任一點(diǎn)上指定了一對(duì)數(shù)fA和fAc，將fA和fAc分別稱為u對(duì)A和Ac的相對(duì)隸屬度.定義如下式的映射[11]

μA(μ):u|→fA∈[0，1]

μAc(μ):u|→fAc∈[0，1]

(6)

μA(u)和μAc(u)分別為u對(duì)A、Ac的相對(duì)隸屬函數(shù).

設(shè)fd=fA-fAc，fd為u對(duì)A的相對(duì)差異度.映射

D(u):u|→fd∈[-1,1]

(7)

D(u)為u對(duì)A的相對(duì)差異函數(shù).

3.1.2利用相對(duì)差異函數(shù)確定相對(duì)隸屬度的方法

對(duì)于地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)因素集為{y1,y2,…,yn}，假設(shè)每個(gè)因素的風(fēng)險(xiǎn)論域可劃分為k個(gè)等級(jí).根據(jù)標(biāo)準(zhǔn)規(guī)范或者評(píng)價(jià)目的，可將第i個(gè)1級(jí)指標(biāo)中第j個(gè)因素的第h個(gè)風(fēng)險(xiǎn)等級(jí)的標(biāo)準(zhǔn)值區(qū)間設(shè)定為[aijh,bijh]，該因素第h個(gè)風(fēng)險(xiǎn)等級(jí)的上下界區(qū)間為[cijh,dijh].而mijh為區(qū)間中D(u)=1的點(diǎn)值，即該因素在第h個(gè)風(fēng)險(xiǎn)等級(jí)的最可能取值.xij為論域上任意點(diǎn)的量值.

根據(jù)可變模糊集理論，地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)的可變模糊集合的吸引域Iab、范圍域Icd和完全隸屬于矩陣Iab的點(diǎn)值矩陣M分別如式(8)～(10)所示：

Iab=([aijh，bijh])

(8)

Icd=([cijh，dijh])

(9)

M=(mijh)

(10)

式中,M的第i行、第h列元素為mijh.根據(jù)h確定mijh：當(dāng)h=1時(shí)mi1=ai1；當(dāng)h=k時(shí)mik=bik；當(dāng)1

若xij落在mijh值的左側(cè)，xij相對(duì)于風(fēng)險(xiǎn)等級(jí)h的隸屬函數(shù)為：

(11)

反之則：

(12)

根據(jù)式(11)、(12)，確定評(píng)價(jià)對(duì)象特征值矩陣X對(duì)各個(gè)風(fēng)險(xiǎn)等級(jí)的相對(duì)隸屬度矩陣，如式(13)所示.

iU=(μAh(xij))

(13)

3.2　基于相對(duì)隸屬函數(shù)的模糊可變?cè)u(píng)價(jià)模型

根據(jù)文獻(xiàn)[8]，地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)的模糊可變?cè)u(píng)價(jià)模型如式(14)所示：

(14)

將式(14)計(jì)算結(jié)果歸一化后得到第i個(gè)1級(jí)指標(biāo)的隸屬度矩陣iUh：

iUh=(iuh)

(15)

式中，iuh為第i個(gè)1級(jí)指標(biāo)歸一化后的相對(duì)于風(fēng)險(xiǎn)等級(jí)h(h=1,2,…,k)的隸屬度；iUh的第i行、第h列元素為iuh.

應(yīng)用級(jí)別特征值公式，1級(jí)指標(biāo)i的級(jí)別特征值為hi，令H=(h1,h2,…,hn)，則：

H=(1,2,…,k)·jUh

(16)

式中，k表示對(duì)應(yīng)的風(fēng)險(xiǎn)等級(jí);“·”表示矩陣相乘.

根據(jù)H對(duì)待評(píng)價(jià)對(duì)象進(jìn)行綜合評(píng)價(jià)得到其級(jí)別特征值：

h0=(ω1,ω2,…ωn)·HT

(17)

式中，ωi為第i個(gè)1級(jí)指標(biāo)的權(quán)重.

4　實(shí)例分析

1)根據(jù)表2、表3確定各指標(biāo)的RP和RI評(píng)價(jià)值，利用式(3)、(4)確定RF值及RFT值，確定利用式(5)確定1級(jí)指標(biāo)及2級(jí)指標(biāo)的權(quán)重，分別如表1的第3，5列所示，如1級(jí)指標(biāo)的權(quán)重確定過(guò)程如表4所示.

表4　1級(jí)指標(biāo)RP和RI評(píng)價(jià)值及其權(quán)重確定過(guò)程

2)根據(jù)該地鐵運(yùn)營(yíng)的實(shí)際情況，采用集值統(tǒng)計(jì)法獲得各指標(biāo)分值，如表2的第6列所示.結(jié)合評(píng)價(jià)對(duì)象，參考文獻(xiàn)[8]及《地鐵運(yùn)營(yíng)安全評(píng)價(jià)標(biāo)準(zhǔn)》中關(guān)于地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)及相關(guān)風(fēng)險(xiǎn)級(jí)別的分級(jí)標(biāo)準(zhǔn)，確定風(fēng)險(xiǎn)等級(jí)論域?yàn)閂={不可接受，基本可接受需改進(jìn)，可接受，可忽略}，以數(shù)值1，2，3和4分別表示這4個(gè)風(fēng)險(xiǎn)等級(jí).在該體系中，取所有指標(biāo)的風(fēng)險(xiǎn)分值區(qū)間均一致，則風(fēng)險(xiǎn)評(píng)價(jià)的可變模糊集合的吸引域?yàn)镮ab=[0,60)[60,80)[80,90)[90,100])，范圍域?yàn)镮cd=([0,80)[0,90)[60,100)[80,100]).根據(jù)指標(biāo)的吸引域區(qū)間及相對(duì)隸屬函數(shù)對(duì)于點(diǎn)值矩陣的要求，確定各指標(biāo)4個(gè)等級(jí)的點(diǎn)值矩陣為(mi1,mi2,mi3,mi4)=(0,65,85,100).

根據(jù)以上的各參數(shù)取值，對(duì)每個(gè)指標(biāo)值應(yīng)用式(11)～(13)計(jì)算，并將結(jié)果歸一化，可得每個(gè)指標(biāo)的相對(duì)隸屬度矩陣，如表1第8～11列.根據(jù)2級(jí)指標(biāo)的權(quán)重及相對(duì)隸屬度，應(yīng)用式(14)～(16)可求得各1級(jí)指標(biāo)的級(jí)別特征值，如表5第2～5行所示；1級(jí)指標(biāo)的權(quán)重及級(jí)別特征值利用式(17)可得到目標(biāo)層的級(jí)別特征值，如表5第6行所示.由表5可知，4種模型下的評(píng)價(jià)結(jié)果的離差只有0.062，評(píng)價(jià)結(jié)果均為3級(jí)(可接受)，偏向于2級(jí)(基本可接受需改進(jìn)).

表5　1級(jí)指標(biāo)4組模型參數(shù)的評(píng)價(jià)結(jié)果表

將該地鐵運(yùn)營(yíng)線路經(jīng)基礎(chǔ)分值評(píng)價(jià)法與可變模糊評(píng)價(jià)所得的兩種評(píng)價(jià)結(jié)果進(jìn)行對(duì)比，如表6.

表6　2種評(píng)價(jià)結(jié)果的對(duì)比

從兩種評(píng)價(jià)方法的評(píng)價(jià)結(jié)果來(lái)看，其總體風(fēng)險(xiǎn)水平都偏向“可接受”，但對(duì)于可變模糊評(píng)價(jià)方法，還能得出期1級(jí)指標(biāo)所處于的風(fēng)險(xiǎn)水平，如對(duì)于“安全管理評(píng)價(jià)”1級(jí)指標(biāo)，其風(fēng)險(xiǎn)處于“需改進(jìn)”狀態(tài)，因此可以通過(guò)完善其所屬風(fēng)險(xiǎn)較高的2級(jí)指標(biāo)來(lái)改善其風(fēng)險(xiǎn)水平，如對(duì)于“事故應(yīng)急預(yù)案及應(yīng)急體系”，該地鐵運(yùn)營(yíng)公司雖然對(duì)應(yīng)急救援預(yù)案進(jìn)行了編制，但預(yù)案的演練方式較為單一，并且演練缺乏相應(yīng)的計(jì)劃、總結(jié)及評(píng)估.因此建議該地鐵運(yùn)營(yíng)公司應(yīng)定期組織多種形式的應(yīng)急演練，并對(duì)其內(nèi)容進(jìn)行總結(jié)、評(píng)估，根據(jù)地鐵運(yùn)營(yíng)發(fā)展現(xiàn)狀對(duì)預(yù)案進(jìn)行維護(hù)更新.

5　結(jié)論

(1)根據(jù)地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)體系中各指標(biāo)對(duì)于風(fēng)險(xiǎn)發(fā)生的可能性及后果嚴(yán)重程度來(lái)確定指標(biāo)權(quán)重的思想是合理性的，采用的三角模糊理論在處理地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)評(píng)價(jià)的不確定性問(wèn)題是可行的；

(2)應(yīng)用可變模糊集理論，用相對(duì)差異函數(shù)來(lái)確定地鐵運(yùn)營(yíng)風(fēng)險(xiǎn)隸屬度向量是可行的.可變模糊評(píng)價(jià)模型不僅權(quán)重向量可變，且模型參數(shù)α、p可變，該方法突破傳統(tǒng)方法只用一種模型對(duì)評(píng)價(jià)對(duì)象進(jìn)行評(píng)價(jià)的作法，可以提高評(píng)價(jià)的可靠性，而從計(jì)算結(jié)果可以看出通過(guò)變換模型參數(shù)得到的級(jí)別特征值較穩(wěn)定.

參考文獻(xiàn):

[1]潘科，石劍云. 變權(quán)和相對(duì)差異函數(shù)在地鐵運(yùn)營(yíng)安全評(píng)價(jià)中的應(yīng)用[J].鐵道學(xué)報(bào)，2009，31(3)：20- 25.

[2]潘科，王洪德，石劍云. 多級(jí)可拓評(píng)價(jià)方法在地鐵運(yùn)營(yíng)安全評(píng)價(jià)中的應(yīng)用[J].鐵道學(xué)報(bào)，2011，33(5):14- 19.

[3]胥旋，鐘茂華，史聰靈，等. 基于動(dòng)態(tài)分級(jí)法的地鐵設(shè)備設(shè)施風(fēng)險(xiǎn)評(píng)價(jià)研究[J]. 2013,23(7)：61- 66.

[4]LI QIMING,DENG YONGLIANG,LIU CONG,et al.Modeling and analysis of subway fire emergency response:An empirical study [J].Safety Science,2016,84：171- 180.

[5]范英，李辰，晉民杰，等.三角模糊數(shù)和層次分析法在風(fēng)險(xiǎn)評(píng)價(jià)中的應(yīng)用研究[J].2014,24(7)：70- 74.

[6]WAN SHUPING,WANG FENG,LIN LIPIAN,et al.Some new generalized aggregation operators for triangular intuitionistic fuzzy numbers and application to multi-attribute group decision making[J].Computers & Industrial Engineering,2016,93：286- 301.

[7]NIETO-MOROTE A,RUZ- VILA F.A fuzzy approach to construction project risk assessment[J].International Journal of Project Managerment,2011,29(2):220- 231.

[8]PAN KE，JIANG QIU-YUAN，YAN SHAN-YU.Variable Fuzzy Assessment Method and Its Application in the Evaluation of Chemical Industry Park Emergency Response Capability[C]//Progress in safety science and technology,Beijing：Science Press，2010:171- 176.