黃光球， 苗 丹

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

近年來(lái)，以高濃度氮氧化物、二氧化硫、PM2.5為核心污染物的有毒霧霾大范圍地、頻繁地在我國(guó)工業(yè)化進(jìn)程較快的城市爆發(fā)，尤其是2013年，多個(gè)城市空氣污染指數(shù)達(dá)到500 μg/m3上限[1]，有毒霧霾已經(jīng)成為當(dāng)前不容忽視的問(wèn)題.對(duì)此，國(guó)內(nèi)外學(xué)者主要從有毒霧霾的危害性評(píng)價(jià)、成因、特征、治理方法等方面進(jìn)行了大量研究，其中對(duì)于有毒霧霾危害性評(píng)價(jià)的研究主要運(yùn)用文字描述、實(shí)驗(yàn)監(jiān)測(cè)、創(chuàng)建模型三類方法.然而，這三類方法在評(píng)價(jià)過(guò)程中都存在一些缺失：文字描述無(wú)法觀測(cè)到有毒霧霾危害產(chǎn)生的過(guò)程，也無(wú)法得到危害產(chǎn)生的大小程度，只能單純地描述有毒霧霾對(duì)人體健康、交通、農(nóng)業(yè)產(chǎn)生的危害[2-5]；實(shí)驗(yàn)監(jiān)測(cè)往往僅針對(duì)霧霾顆粒中重金屬的破壞作用[6]，或者有毒霧霾對(duì)植物光合作用的影響[7]，亦或是有毒霧霾與農(nóng)作物病害發(fā)生率之間的關(guān)系[8]，忽略了有毒霧霾危害的全面性和生態(tài)環(huán)境的整體性；創(chuàng)建的物元可拓模型[9]、故障樹(shù)模型[10]等對(duì)于有毒霧霾危害性的評(píng)價(jià)都是從霧霾產(chǎn)生的角度出發(fā)，這樣的評(píng)價(jià)方式趨向于原因?qū)虻脑u(píng)價(jià)，難以對(duì)已經(jīng)發(fā)生的有毒霧霾危害性進(jìn)行估量評(píng)價(jià).

Petri網(wǎng)既有直觀的圖形表示，又有堅(jiān)實(shí)的數(shù)學(xué)基礎(chǔ)[11-12]，具有很好的評(píng)價(jià)功能，但是卻很少用于危害性評(píng)價(jià)領(lǐng)域.本文通過(guò)發(fā)揮Petri網(wǎng)的特性，用形式化的方法描述有毒霧霾中的毒性物質(zhì)在地表生態(tài)環(huán)境危害產(chǎn)生的過(guò)程，且融入危害性評(píng)價(jià)指標(biāo)的計(jì)算，得出危害性評(píng)價(jià)指標(biāo)的量化值，對(duì)有毒霧霾的危害性進(jìn)行定量評(píng)價(jià)，并結(jié)合危害等級(jí)標(biāo)準(zhǔn)，在定量評(píng)價(jià)的同時(shí)進(jìn)行合理的定性評(píng)價(jià).

1 函數(shù)Petri網(wǎng)的概念

1.1 函數(shù)Petri網(wǎng)的定義

1962年，Carl Adam Petri首次提出了適用于描述離散、異步、并發(fā)系統(tǒng)模型的Petri網(wǎng)[11]，函數(shù)Petri網(wǎng)是在Petri網(wǎng)的基礎(chǔ)上進(jìn)化來(lái)的.一個(gè)九元組FPN=(S，T，W，M，N，D，λ，V，ω)為函數(shù)Petri網(wǎng)，其中S={s1，s2，…，sm}稱為FPN的庫(kù)所集，T={t1，t2，…，tm}稱為FPN的變遷集.W為PN上的權(quán)函數(shù)，在函數(shù)Petri網(wǎng)中，W=1.M表示FPN的標(biāo)識(shí)，它是一個(gè)列向量，其中第i個(gè)元素表示第i個(gè)庫(kù)所中的托肯數(shù)目，M0為初始標(biāo)識(shí).N是事件累積發(fā)生次數(shù)的有限級(jí)，N={1，2，…，k}，k≤365.D是事件次數(shù)，00，n>0，m>D，a、n、m均為常數(shù);② 單調(diào)遞減的變遷速率，λ=a*exp[(m-D)/n]，a>0，n>0，m

1.2 函數(shù)Petri網(wǎng)的圖形表達(dá)

在Petri網(wǎng)中，庫(kù)所用圓圈表示，變遷用長(zhǎng)方形表示，托肯由庫(kù)所中的小黑點(diǎn)表示.為了方便表達(dá)起過(guò)渡作用的、沒(méi)有實(shí)際意義的、為滿足某些條件而人為添加的庫(kù)所和變遷稱為虛庫(kù)所和虛變遷，用虛線圓圈表示虛庫(kù)所，用虛線長(zhǎng)方形表示虛變遷.為了區(qū)別函數(shù)形式不同的變遷速率，將單調(diào)遞增的變遷速率用λ↑表示，單調(diào)遞減的變遷率用λ↓表示，恒定不變的變遷速率用λ表示.影響變遷觸發(fā)的前集庫(kù)所對(duì)后繼庫(kù)所的貢獻(xiàn)值ω若為0或者-1時(shí)，在庫(kù)所指向變遷的箭線末端表示，沒(méi)有表示的均默認(rèn)為1.

1.3 函數(shù)Petri網(wǎng)變遷觸發(fā)規(guī)則

九元組FPN=(S，T，W，M，N，D，λ，V，ω)是一個(gè)函數(shù)Petri網(wǎng)，不同庫(kù)所的標(biāo)識(shí)和狀態(tài)值通過(guò)變遷的觸發(fā)來(lái)改變，而任何變遷的觸發(fā)都必須滿足一定的使能條件，即?s∈?t，M(s)≥W(s，t)，即M(s)≥1時(shí)，變遷t∈T才能使能.若一個(gè)庫(kù)所的標(biāo)識(shí)同時(shí)滿足兩個(gè)變遷使能條件，用空心長(zhǎng)方形來(lái)表示觸發(fā)優(yōu)先級(jí)別高的變遷，即優(yōu)先觸發(fā)空心長(zhǎng)方形表示的變遷.變遷觸發(fā)規(guī)則：

若s∈?t-t?，M′(s)取M(s)-1，若s∈t?-?t，M′取M(s)，其他情況下M′(s)取M(s).

如圖1(a)所示，?ti為si，M(si)=0<1，則ti不滿足使能條件；?tj為sj，M(sj)=1，tj滿足使能條件，因此變遷tj觸發(fā)，觸發(fā)后結(jié)果如圖1(b).

如圖2(a)所示，?ti為si和sj，M(si)=1且M(sj)=1，則ti滿足使能條件，tj也滿足使能條件.sj庫(kù)所的標(biāo)識(shí)同時(shí)滿足兩個(gè)變遷的使能條件，則優(yōu)先觸發(fā)空心長(zhǎng)方形表示的變遷ti.ti觸發(fā)后的結(jié)果如圖2(b)，仍有變遷滿足使能條件，此時(shí)變遷tj觸發(fā)，觸發(fā)結(jié)果如圖2(c).

1.4 函數(shù)Petri網(wǎng)庫(kù)所狀態(tài)值的函數(shù)表達(dá)式

用函數(shù)Petri網(wǎng)描述一種因果關(guān)系，庫(kù)所的狀態(tài)值與輸入有關(guān)，與輸出無(wú)關(guān).輸入包括輸入變遷以及與輸入變遷相連的前集庫(kù)所的狀態(tài)值，為了更好地進(jìn)行計(jì)算，將庫(kù)所狀態(tài)值的改變規(guī)則分為四大類.

(1) 單個(gè)實(shí)變遷輸入的庫(kù)所的狀態(tài)值：如圖3所示，庫(kù)所sj只有通過(guò)變遷ti觸發(fā)產(chǎn)生的狀態(tài)值輸入，?ti為si，當(dāng)ti觸發(fā)，sj的狀態(tài)值函數(shù)：VD(sj)=VD-1(sj)+VD(si)*λD，i*ω(si，ti).

(2) 虛變遷輸入的虛庫(kù)所的狀態(tài)值：如圖4所示，虛庫(kù)所sj只通過(guò)虛變遷ti觸發(fā)產(chǎn)生的狀態(tài)值輸入，?ti為si，當(dāng)ti觸發(fā)，sj的狀態(tài)值函數(shù)VD(sj)=VD(si)*λD，i*ω(si,ti).

(3) 兩個(gè)輸入變遷，變遷屬于順序關(guān)系：如圖5所示，sj的狀態(tài)值受到變遷ti和tj的共同影響，而且tj觸發(fā)后，ti才滿足觸發(fā)的條件，當(dāng)兩個(gè)變遷相繼觸發(fā)后，sj的狀態(tài)值函數(shù)表達(dá)式為：

VD(sj)=VD(sk)*ω(sk,tk)*λD，k+

[VD(sk)*ω(sk,tk)*λD，k*ω(sj,ti)+VD(si)*ω(si,ti)]*λD，i.

(4) 多個(gè)輸入變遷，變遷屬于并列關(guān)系：如圖6所示，sj的狀態(tài)值受到變遷t1，t2,…,tk的共同影響，當(dāng)變遷觸發(fā)，sj的狀態(tài)值函數(shù)表達(dá)式為：

2 函數(shù)Petri網(wǎng)危害性評(píng)價(jià)模型

2.1 危害性評(píng)價(jià)模型

(1) 有毒霧霾對(duì)地表自然生態(tài)環(huán)境危害性評(píng)價(jià)模型的基本假設(shè)：① 有毒霧霾的危害具有區(qū)域性，危害在某一區(qū)域內(nèi)是全面覆蓋的.② 忽略有毒霧霾污染物濃度對(duì)評(píng)價(jià)指標(biāo)的影響.③ 忽略毒性物質(zhì)在生態(tài)環(huán)境子系統(tǒng)間的相互作用.④ 變遷速率的函數(shù)關(guān)系式是可以通過(guò)監(jiān)測(cè)或者數(shù)據(jù)的規(guī)律性得到的.

表1 模型中庫(kù)所和變遷名稱、庫(kù)所初始狀態(tài)值表

(2) 有毒霧霾對(duì)地表自然生態(tài)環(huán)境危害：有毒霧霾的每一次出現(xiàn)都會(huì)攜帶各種毒性物質(zhì)，比如氮氧化物、重金屬、二氧化硫等，這些物質(zhì)通過(guò)沉降、降雨等方式進(jìn)入地表的各個(gè)自然生態(tài)環(huán)境子系統(tǒng)，直接影響生物的生長(zhǎng)發(fā)育，重者會(huì)造成生物的徹底滅亡，使各個(gè)子系統(tǒng)受到不同程度的危害.氮氧化物進(jìn)入水體，造成水體富營(yíng)養(yǎng)化，增加水生生物的死亡率，降低水的利用價(jià)值.重金屬和二氧化硫?qū)χ参锏挠绊戄^大：進(jìn)入草地，就會(huì)降低草的抵抗力，導(dǎo)致草地面積退化；污染森林土壤，影響樹(shù)木正常發(fā)育，造成樹(shù)木枯萎死亡；土地退化和樹(shù)木的減少導(dǎo)致水土大面積流失；污染農(nóng)田土壤，抑制農(nóng)作物吸收養(yǎng)分，造成農(nóng)作物減產(chǎn)絕收.地表自然生態(tài)系統(tǒng)是一個(gè)復(fù)雜且各子系統(tǒng)相互聯(lián)系的大系統(tǒng)，在有毒霧霾中的毒性物質(zhì)進(jìn)入自然生態(tài)環(huán)境的同時(shí)，生態(tài)環(huán)境會(huì)表現(xiàn)出一定的修復(fù)能力，能進(jìn)行自我保護(hù)，以此來(lái)削弱有毒霧霾的危害程度.

(3) 有毒霧霾對(duì)地表自然生態(tài)環(huán)境的危害性評(píng)價(jià)體系：通過(guò)對(duì)有毒霧霾危害產(chǎn)生過(guò)程的分析，將富營(yíng)養(yǎng)化湖泊面積、草地退化面積、森林退化面積、農(nóng)作物減產(chǎn)量、水土流失面積作為危害性評(píng)價(jià)的子指標(biāo)，地表自然生態(tài)環(huán)境危害損失為綜合指標(biāo).

(4) 有毒霧霾對(duì)地表自然生態(tài)環(huán)境危害性評(píng)價(jià)模型：結(jié)合函數(shù)Petri網(wǎng)的定義，將有毒霧霾對(duì)地表自然生態(tài)環(huán)境產(chǎn)生的危害過(guò)程，環(huán)境的自我修復(fù)影響過(guò)程，以及危害性評(píng)價(jià)體系，通過(guò)函數(shù)Petri網(wǎng)的網(wǎng)系統(tǒng)表達(dá)出來(lái)，建立有毒霧霾對(duì)地表自然生態(tài)環(huán)境的危害性評(píng)價(jià)研究模型，如圖7所示.圖7中庫(kù)所和變遷的具體名稱，如表1所示.

2.2 危害度定義

為了更好地表達(dá)有毒霧霾對(duì)地表自然生態(tài)環(huán)境的危害性，將危害性用危害度表示.考慮到霧霾具有區(qū)域性，霧霾對(duì)生態(tài)環(huán)境產(chǎn)生的危害是過(guò)程性的，積累性的，因此，本文采用公式(1)來(lái)計(jì)算危害度：

H=EL/Gi,

(1)

式中：H為危害度；EL是年有毒霧霾造成的地表自然生態(tài)環(huán)境的經(jīng)濟(jì)損失，即為危害性評(píng)價(jià)模型中D=N時(shí)的VD(S23)，Gi為研究地區(qū)研究年份國(guó)內(nèi)生產(chǎn)總值.

根據(jù)危害度計(jì)算結(jié)果，將有毒霧霾對(duì)地表自然生態(tài)環(huán)境的危害度劃分為極嚴(yán)重(H>0.1%)、嚴(yán)重(0.05%

2.3 危害性評(píng)價(jià)的基本步驟

(1) 明確隨危害事件次數(shù)變化的變遷速率滿足的函數(shù)關(guān)系式.(2) 設(shè)最大迭代次數(shù)d即為事件累積發(fā)生次數(shù)k，當(dāng)前迭代次數(shù)d=0，事件次數(shù)D=0.(3) 令迭代次數(shù)d=d+1，事件次數(shù)D=d.通過(guò)步驟(1)得到相應(yīng)的變遷速率，結(jié)合函數(shù)Petri網(wǎng)的變遷發(fā)生規(guī)則運(yùn)行Petri網(wǎng)，變遷發(fā)生后，通過(guò)庫(kù)所狀態(tài)值的函數(shù)表達(dá)式計(jì)算出各庫(kù)所的狀態(tài)值.(4) 若d

3 實(shí)例說(shuō)明與分析

某縣區(qū)，占地面積1 279.42 km2，由于周邊多礦產(chǎn)開(kāi)發(fā)，工業(yè)廠房較多，汽車(chē)數(shù)量增長(zhǎng)較快等因素，致使當(dāng)?shù)仡l繁受到有毒霧霾的侵襲，近50年來(lái)霧霾年發(fā)生次數(shù)在20～170次之間.接連不斷的霧霾不僅給當(dāng)?shù)鼐用竦娜粘Ｉ钤斐闪藰O大困擾，更使得當(dāng)?shù)氐淖匀簧鷳B(tài)環(huán)境受到了極大的破壞.2016年，該縣霧霾發(fā)生97次，生產(chǎn)總值為513.8億元.各變遷速率的函數(shù)關(guān)系式如表2所示.

表2 各變遷速率的函數(shù)關(guān)系式表

分析案例內(nèi)容，N=97，G2016=513.8億元.為了使結(jié)果不僅能評(píng)價(jià)當(dāng)年的危害性，還可以直觀的凸顯危害性與有毒霧霾發(fā)生次數(shù)的動(dòng)態(tài)關(guān)系.考慮到當(dāng)?shù)亟?0年來(lái)霧霾年發(fā)生的次數(shù)，將有毒霧霾對(duì)地表自然生態(tài)環(huán)境危害性評(píng)價(jià)模型和函數(shù)Petri網(wǎng)庫(kù)所狀態(tài)值的函數(shù)關(guān)系式結(jié)合起來(lái)，通過(guò)MATLAB軟件計(jì)算得到有毒霧霾累積發(fā)生次數(shù)在185次內(nèi)的危害性評(píng)價(jià)指標(biāo)所有數(shù)據(jù)，其中霧霾發(fā)生97次時(shí)，富營(yíng)養(yǎng)化湖泊面積為0.937 39 km2，草地退化面積為20.337 45 km2，森林退化面積為14.278 28 km2，農(nóng)作物減產(chǎn)量為97.362 62 t，水土流失面積為187.356 9 km2，地表自然生態(tài)環(huán)境危害損失為2 268.131 82萬(wàn)元，并將有毒霧霾發(fā)生次數(shù)在185次以內(nèi)的所有危害度用圖8表示.

從圖8可以得到兩點(diǎn)結(jié)論.第一，按照危害性等級(jí)標(biāo)準(zhǔn)，年發(fā)生次數(shù)在0～28次時(shí)，有毒霧霾的危害性等級(jí)為較輕；年發(fā)生次數(shù)在29～106次時(shí)，危害性等級(jí)為較嚴(yán)重；年發(fā)生次數(shù)在107～168時(shí)，危害性等級(jí)為嚴(yán)重；年發(fā)生次數(shù)大于等于169次時(shí)，危害性等級(jí)為極嚴(yán)重.該縣霧霾發(fā)生次數(shù)為97次，危害性為較嚴(yán)重.第二，圖8準(zhǔn)確地反應(yīng)出有毒霧霾對(duì)生態(tài)環(huán)境產(chǎn)生的危害具有邊際效應(yīng)遞增特征：有毒霧霾不僅給生態(tài)環(huán)境造成損害，還減弱了生物的生存能力、抵抗能力以及自我修復(fù)能力，使霧霾產(chǎn)生的危害越來(lái)越大.

通過(guò)上述得到的結(jié)果，給減輕有毒霧霾危害提供了新的思路：① 在明確有毒霧霾造成的各生態(tài)子系統(tǒng)的破壞量后，可以有針對(duì)性地進(jìn)行生態(tài)環(huán)境保護(hù)和生態(tài)再造.② 對(duì)于報(bào)道相關(guān)氣象災(zāi)害的媒體或者組織，應(yīng)當(dāng)實(shí)時(shí)報(bào)道霧霾發(fā)生的年累計(jì)次數(shù)，讓公眾更清晰地了解到問(wèn)題的嚴(yán)重性，增強(qiáng)環(huán)保意識(shí).③ 不僅可以從源頭上控制霧霾的毒性，也可以通過(guò)生態(tài)建設(shè)增強(qiáng)生物自我修復(fù)能力的方式減小有毒霧霾對(duì)地表自然生態(tài)環(huán)境的危害，比如選擇種性優(yōu)良、產(chǎn)量高、品質(zhì)好、抗逆性強(qiáng)的品種進(jìn)行林業(yè)和草業(yè)種植.

4 結(jié) 論

本文結(jié)合Petri網(wǎng)的特點(diǎn)，提出了一種針對(duì)霧霾中毒性物質(zhì)對(duì)地表自然生態(tài)環(huán)境的危害性定量評(píng)價(jià)的方法.重點(diǎn)分析了有毒霧霾產(chǎn)生危害的過(guò)程，預(yù)見(jiàn)了研究地區(qū)研究年份地表自然生態(tài)環(huán)境在有毒霧霾的侵蝕下造成經(jīng)濟(jì)損失的累積趨勢(shì)，定量化評(píng)估有毒霧霾造成的地表自然生態(tài)子系統(tǒng)環(huán)境損失量，為政府生態(tài)修復(fù)提供依據(jù)；定量化評(píng)估有毒霧霾造成的地表自然生態(tài)環(huán)境經(jīng)濟(jì)損失，使公眾直觀地認(rèn)識(shí)到環(huán)境問(wèn)題的嚴(yán)重性；危害度的定量定性化評(píng)價(jià)有助于科學(xué)判別環(huán)境政策的有效性，為政府完善相關(guān)政策提供決策依據(jù).