基于風險場的區域突發性環境風險評價方法研究

邢永健,王 旭*,可 欣,吳天輝,滕文超(.沈陽航空航天大學安全工程學院,遼寧 沈陽 036；.沈陽航空航天大學能源與環境學院,遼寧 沈陽 036)

?

基于風險場的區域突發性環境風險評價方法研究

邢永健1,王 旭1*,可 欣2,吳天輝1,滕文超2(1.沈陽航空航天大學安全工程學院,遼寧 沈陽 110136；2.沈陽航空航天大學能源與環境學院,遼寧 沈陽 110136)

摘要：在大尺度空間,風險源、風險受體和風險傳播途徑的多樣性和相互作用的復雜性是區域環境風險評價的難點.同時,現有評價方法很少關注風險因子的釋放規律,釋放后在空間中的分布格局,以及風險受體受到損害的途徑和程度.本文應用風險場理論,分析描述風險場形成和對風險受體作用的機制.結合南京化工園區實例進行環境風險源識別,采用集對分析等方法構建了各風險源產生的環境風險場,再分析處于風險場中的環境風險受體,最后得到區域環境風險水平分布,同時將其劃分為5個等級.評價結果顯示,位于風險源周圍或河道下游的人口稠密、生態環境敏感地區環境風險值R≥6,風險水平處于極高、高等級,與園區現狀具有較好的一致性.建議依據環境風險水平的分級、分類結果建立起以預防為主的環境風險管理體系.

關鍵詞：風險場；區域環境風險；風險評價；風險受體；集對分析法

* 責任作者, 教授, wangx1960@126.com

區域風險評價(Regional Risk Assessment)是一種估計和比較環境影響在大尺度地理區域特征的方法[1],其特征表現在多源釋放產生的多重壓力、風險傳播途徑以及風險受體的多樣性和相互作用的復雜性.

近年來,國內外研究對區域環境風險評價進行了多方面的探索.通過構建風險模型對環境風險后果進行分析,進而提出基于后果的環境風險評價模型,如考慮危險物質特性、環境遷移規律和環境危害特性的評價模型[2],以風險受體脆弱性與暴露的乘積表示風險后果的評價模型[3],考慮重大損失的全部重要類型來評估整體后果的評價模型[4].這些模型的實質是分析風險受體以及受損規律來評價風險源,但未能足夠闡述風險源對風險受體的作用途徑.

同時,對環境風險源和受體的綜合評價模型以及針對環境風險受體的評價模型有了進一步的發展.賈倩等[5]構建了基于危險物質、生產工藝、設備設施、企業管理、企業布局的突發環境風險綜合評價指標體系和相應的風險評價模型;Pizzol等[6]針對空間信息變化問題,提出了一種基于空間分析和相對風險的評價方法,對污染場地的風險進行大小排序,為風險管理提供依據;Chen等[7]開發了特定受體風險分布地圖以更好地對污染土地管理進行情景分析; Giubilato 等[8]提出了一種在區域尺度上基于證據權重的環境化學風險源的分級方法,支持管理者識別優先環境污染物和優先管理的區域; Zabeo等[9]通過應用多目標決策分析和空間分析技術評價區域環境風險受體的脆弱性;薛鵬麗等[10]從環境風險受體敏感性和適應力兩方面構建了脆弱性概念模型,建立了上海市環境污染事故風險受體綜合脆弱性評價指標體系.總體而言,現有評價模型和方法很少關注風險因子如何釋放,釋放后在空間中的分布格局,以及風險受體如何暴露于風險因子并受到損害的程度.

本文探索應用區域環境風險系統理論中的環境風險場理論[11],并對其進行實際應用的可操作性優化,以期彌補相關研究的不足.

1　理論分析

1.1 環境風險場

風險因子在環境空間中形成某種分布格局是風險危害發生的前提,我們稱這種分布格局為“環境風險場”(Environmental Risk Field, ERF).由于環境風險是一種不確定性事件,因此,環境風險場的形成也具有一定的概率.在環境風險場形成過程中,風險因子的狀態、時空特性都將發生變化.

在實際的風險評價中,關注的是ERF對風險受體能產生危害的能力,即ERF的強度,這里將ERF的強度稱為環境風險場強度,簡稱場強,用E表示.ERF中的任一處(x, y)質點的風險場強度將決定于該點處風險因子的暴露水平和這種暴露出現的概率.公式表示如下:

其中:

式中:Ex,y為(x, y)處質點的風險場強;Ex,yi為質點(x, y)處,風險因子i的風險場強;Px,yi為風險因子i出現在質點(x, y)處的概率;Cx,yj為質點(x, y)處,風險因子i的暴露水平;P′x,yi為出現在質點(x, y)處的風險因子i的釋放概率;P″x,yi為控制機制失效的概率;P″′x,yi為風險因子i轉運到質點(x, y)處的概率;i為質點(x, y)處可能出現的風險因子,i= 1,2,...,n.

1.2 環境風險受體

風險受體是指風險因子可能危害的人、有價值物體、自然環境及社會系統[12].風險受體易損性可以被定義為受體對風險因子的敏感性,它是風險受體固有的特性,與風險因子種類和受體本身有關[9].如果要考慮質點(x, y)處所有可能出現的風險受體,那么某種風險受體與其出現的概率必須綜合考量.因此,質點(x, y)處風險受體易損性可表示為

式中:Vx,y為質點(x, y)處的風險受體易損性;Px,yj為風險受體j出現在質點(x, y)處的概率;Vx,yj為質點(x, y)處,風險受體j的易損性;j為質點(x, y)處可能出現的風險受體,j=1,2,...,k.

根據環境風險系統理論,風險場與風險源之間有著密切的因果關系,風險場中的風險因子源于風險源;而風險場與風險受體二者有相對獨立性,并不相互依存.鑒于此,在具體的環境風險評價中,為了實現可操作性,我們可先假設風險受體處于勻強風險場中,這樣只需根據風險受體本身的特性就可得出受體易損性.

1.3 風險計算模型

空間中質點(x,y)處的風險與該處可能出現的風險場強和風險受體易損性共同決定,計算模

型為

式中:Rx,y為質點(x, y)處的風險;Ex,yi為質點(x, y) 處,風險因子i的風險場強;Px,yj為風險受體j出現在質點(x, y)處的概率;Vx,yj為質點(x, y)處,風險受體j的易損性;i為質點(x, y)處可能出現的風險因子,i=1,2,...,n;j為質點(x, y)處可能出現的風險受體,j=1,2,...,k.

由于風險場(風險因子)與風險受體的多樣性以及相互作用的復雜性,為簡化計算,質點(x, y)處的風險可表示為

式中:Rx,y為質點(x, y)處的風險;Ex,y為質點(x, y)處的風險場強;Vx,y為質點(x, y)處的風險受體易損性.

綜上所述,把質點(x, y)處的Ex,y和Vx,y作為Rx,y的變量,但是Ex,y和Vx,y如何計算更具科學性,目前還沒有定論.劉桂友等[13]參考USEPA規定(職業人群可接受風險值為10-5～ 10-4p/a,非職業人群可接受風險值為10-7～10-6p/a),對風險值在10-8～10-3p/a進行分級.同時,為方便在風險信息矩陣中顯示風險值,把指數形式的風險值轉換為熟悉的小數形式.類似的,我們將二者相乘后轉化為小數形式(式(6))進行表達.

2　案例分析

南京化工園區位于長江北岸,目前規劃面積45km2(包括長蘆片區26km2,玉帶片區19km2),實際開發面積29.2km2.規劃區內以岳子河為界,岳子河以西為長蘆片,岳子河以東為玉帶片.南京化工園區是以石油化工和合成材料、精細化工等石油深加工為主的綜合性化工園,園區內及外圍分布若干居民區、生態保護區,且河網密布,一旦發生突發環境事件,極易造成重大環境損失.

2.1 環境風險場的構建

2.1.1 區域網格化 區域網格化是解決區域風險評價的復雜性問題的方法之一[13-14].把評價區域的二維空間用500m等步長劃分為相同大小的正方形區域,然后用一個62行65列的矩陣來表示評價區域,矩陣中的元素Ex,y代表二維空間中對應的正方形區域的環境風險場強度,通常用正方形中心點的場強來代表該區域.利用網格法的思想,可把研究區域分割形成一系列正方形區域,整個區域環境風險場強度可以用矩陣形式(式(7))表示,此時矩陣內各單元的環境風險場強度為0.

2.1.2 環境風險源識別 風險場由風險源產生,合理辨識環境風險源是構建環境風險場的前提.因研究區域范圍較大,可視一家企業為一個環境風險源單元.突發性環境風險源危險性分析不是基于危害后果進行的,而是基于風險源自身的危險性進行的,突發環境風險源危險的大小與化學物質的理化性質、危險性和物質數量的多少密切相關.環境風險源中有毒有害物質在發生泄漏后會對環境和人群帶來影響,但有些易燃易爆物質發生燃燒爆炸時可能引發多米諾效應和次生、衍生事故.因此,參照《企業突發環境事件風險評估指南(試行)》[15]確定環境風險源指數Q.

式中:qi為每種危險化學品實際存在量,t;Qi為與各危險化學品相對應的臨界量,t.

不同的風險物質有不同的理化性質,對環境風險受體的作用機制和損害程度也不同,該方法應用歸一化的思想,把不同物質的風險量度同一化.同時,環境風險物質臨界量的確定遵循危害等值的原則,這為后續的風險疊加提供了可行性.

由式(8)計算得出園區各環境風險源指數.根據系統安全優先次序,也為優化評價方法,篩選風險源指數最大的前25家企業作為評價對象.表1所列是園區部分企業的環境風險源信息以及在矩陣中的位置.

表1　南京化工園區部分企業環境風險源信息Table 1　Environmental risk sources information of partial enterprises in Nanjing Chemical Industry Park

表2　南京化工園區近20年年均風頻統計(%)Table 2 The average of annual wind frequency in Nanjing Chemical Industry Park in past 20years (%)

2.1.3 環境風險場強度計算 環境風險場按風險因子傳播途徑可分為大氣環境風險場、水環境風險場和土壤環境風險場.土壤環境風險場因其時間跨度大,在研究突發性環境風險時不予考慮.

對大氣環境風險因子分布可按梯形模糊關系進行簡化計算[13],這為構建大氣環境風險場強度計算模型提供了依據.同時,環境風險場的形成存在確定與不確定的因素,如風險因子的載體大氣是確定的,但大氣在事故狀態下如何流動是不確定的,而且確定與不確定的演變是一個連續的、動態的過程,由于時空的變化,環境風險場在不斷改變,也是一個動態的過程.因此,本文應用集對分析法[16]來構建大氣環境風險場強度計算模型.評價區域地處長江三角洲,在此,假設該區域地勢平坦開闊,且忽略人工建筑對氣體擴散的影響.那么,區域某單元單一風險因子的暴露水平可表示為[17]

式中:Cx,y為計算點的單一風險因子的暴露水平;Q為風險源點的環境風險源指數;l為計算點與風險源點的距離;i、j分別為差異系數、對立系數,取i1=0.5、i2=-0.5、j≡-1;s1、s2、s3、s4分別取1000,3000,5000,10000m.

水環境風險因子擴散基本沿河道遷移,所影響區域僅限河道周邊[13],形成的水環境風險場呈現明顯的流域性特征.在有限的區域內有毒有害液體濃度仍遠大于容許濃度值[13],故對水環境風險因子暴露水平做等值擴散處理.

式(2)中的P′、P″可通過歷史數據統計得到,P′′′則需具體分析.《環境風險評價實用技術、方法和案例》中有關石化行業的事故概率統計為10-5/a級別[18].在研究大氣環境風險場時,P′P″可取10-5/a,P′′′等于研究區域年平均風頻,見表2;在研究水環境風險場時,因其較大氣環境風險場難以形成,P′P″P′′′可取10-6/a.不同企業的P不盡相同,可根據《企業突發環境事件風險評估指南(試行)》中企業生產工藝過程與環境風險控制水平進行修正[15],P的修正系數見表3.因缺少各企業資料,本文未作修正.

表3　企業生產工藝與環境風險控制水平以及P的修正系數Table 3　Enterprise production process and environmental risk control levels and correction factors of P

分別計算25個風險源所形成大氣、水環境風險場強度,并由式(1)得到區域各單元大氣、水環境風險場強度.若某單元同時存在大氣、水環境風險場,取其大者為最終環境風險場強度.因篇幅有限,僅取(17, 26)單元周圍1km各單元風險場強度矩陣為例示之,如下.

2.2 環境風險受體易損性計算

同2.1.1,得到整個區域環境風險受體易損性指數矩陣(式(10)),此時矩陣內各單元的環境風險受體易損性指數為0.

在環境風險評價中,很少全面考慮和定量分析受體易損性.區域尺度的環境風險受體易損性差異明顯,暴露于環境風險場時后果各異,風險管理也顯著不同.因此,合理量化環境風險受體易損性意義重大.由于暴露受體的敏感性越強,易損性就越大;而適應力越強,則易損性越低[10],那么環境風險受體易損性計算模型可表示為

式中:V、S和A分別為風險受體易損性指數、敏感性指數和適應力指數.

表4　S和A量化和計算模型Table 4　S, A quantification and calculation model

根據表4,通過式(3) (受體做靜態處理時, Px,yj=0或1)分別計算區域各單元環境風險受體易損性指數,下面是(17, 26)單元周圍1km各單元環境風險受體易損性指數矩陣.

2.3 環境風險值計算

根據式(6),計算區域各單元的環境風險值,并通過去余取整的方法得到區域環境風險值矩陣,(17, 26)單元周圍1km各單元環境風險值矩陣如下.

應用ArcGIS進行可視化處理,繪制環境風險地圖(圖1).

根據園區環境風險值計算和圖1表征結果,區域環境風險水平也可以劃分為5個等級水平:極高(R≥7)、高(R=6)、中(R=5)、低(R=4)、極低(R≤3).

圖1　南京化工園區及周邊環境風險地圖Fig.1　Environmental risk map for Nanjing Chemical Industry Park and the surrounding

3　環境風險管理

依據行政區進行環境風險管理是一種有效的手段[6-7,9,12].根據式(12)計算南京化工園區及周邊行政區的環境風險值(四舍五入),并進行環境風險分級,同時通過分析環境風險源和受體分布確定行政區的顯著風險類型,見表5.

式中:n0、n1、n2、…為行政區內各單元相同風險值的數量.

大廠街道、葛塘街道為高風險地區.該地區分布大量居民區,人口密集,有馬汊河-長江生態公益林二級管控區,且位于化工園區盛行風向的下風向.應做好環境應急人員與隊伍建設,制定詳細的人員疏散方案和救援工作,完善應急監測及預警決策指揮等.

評價區域環境風險源主要位于長蘆街道,該地區多個環境風險源瀕臨長江,一旦有毒有害物質流入長江,極易對下游飲用水水源地、生態濕地造成破壞性影響,因此,要做好各環境風險源的風險評價,優化產業布局,并制定有針對性的應急預案,防止火災、爆炸、泄漏等事故及可能引起的次生、衍生環境污染及人員傷亡事故(有毒有害氣體擴散,消防水、物料泄漏物及反應生成物,從雨水排口、清凈下水排口、污水排口、廠門或圍墻排出,污染生態環境等).

表5　南京化工園區及周邊行政區的環境風險值Table 5　Environmental risk values for districts of Nanjing Chemical Industry Park and the surrounding

其他街道緊鄰化工園區,人口密集、生態環境敏感,應根據區域環境風險評價結果優化城市規劃與土地利用規劃,制定突發環境事件應急預案,預防為主,完善環境風險管理體系.

4　結語

在工業園區規劃建設如火如荼發展的今天,開展區域突發性環境風險評價方法研究是具有現實意義的.本文根據風險場理論分析描述風險場的形成及其對風險受體的損害機制.結合實例將研究區域網格化,進行環境風險源識別,運用集對分析等方法構建環境風險場強度、環境風險受體易損性指數和環境風險值計算模型,繼而確定了區域環境風險水平分布.研究表明,風險場理論運用于區域環境風險評價是可行的,其量化后的計算結果與園區現狀具有較好的一致性.根據結果判定園區及周邊宜采取分類、分級環境風險管理策略.

參考文獻：

[1] Hunsaker C T, Graham R L, Suter G W II, et al. Assessing ecological risk on a regional scale [J]. Environmental Management, 1990,14(3):325-332.

[2] 趙 肖,郭振仁.基于環境后果評價的環境風險源分級模型研究[J]. 安全與環境學報, 2010,10(2):105-108.

[3] 周 薇,錢 瑜,李林子,等.基于環境風險后果評價的危險化學品碼頭選址 [J]. 環境保護科學, 2012,38(5):59-64.

[4] Arunraj N S, Maiti J. A methodology for overall consequence modeling in chemical industry [J]. Journal of Hazardous Materials, 2009,169(1):556-574.

[5] 賈 倩,黃 蕾,袁增偉,等.石化企業突發環境風險評價與分級方法研究 [J]. 環境科學學報, 2010,30(7):1510-1517.

[6] Pizzol L, Critto A, Agostini P, et al. Regional risk assessment for contaminated sites Part 2: Ranking of potentially contaminated sites [J]. Environment International, 2011,37(8):1307-1320.

[7] Chen I C, Ng S, Wang G S, et al. Application of receptor-specific risk distribution in the arsenic contaminated land management [J]. Journal of Hazardous Materials, 2013,262(15):1080-1090.

[8] Giubilato E, Zabeo A, Critto A, et al. A risk-based methodology for ranking environmental chemical stressors at the regional scale [J]. Environment International, 2014,65:41-53.

[9] Zabeo A, Pizzol L, Agostini P, et al. Regional risk assessment for contaminated sites Part 1: Vulnerability assessment by multicriteria decision analysis [J]. Environment International, 2011,37(8):1295-1306.

[10] 薛鵬麗,曾維華.上海市環境污染事故風險受體脆弱性評價 [J].環境科學學報, 2011,31(11):2556-2561.

[11] 畢 軍,楊 潔,李其亮.區域環境風險分析和管理 [M]. 北京:中國環境科學出版社, 2006.

[12] Albanese S, Vivo B D, Lima A, et al. Prioritising environmental risk at the regional scale by a GIS aided technique: The Zambian Copperbelt Province case study [J]. Journal of Geochemical Exploration, 2014,144:433-442.

[13] 劉桂友,徐琳瑜.一種區域環境風險評價方法——信息擴散法[J]. 環境科學學報, 2007,27(9):1549-1556.

[14] 謝元博,李 巍,郝芳華.基于區域環境風險評價的產業布局規劃優化研究 [J]. 中國環境科學, 2013,33(3):560-568.

[15] 環境保護部.企業突發環境事件風險評估指南(試行) [S/OL]. http://www.zhb.gov.cn/gkml/hbb/bgt/201404/W02014041554313 9322805.pdf, 2014-04-03.

[16] 趙克勤.集對分析及其初步應用 [M]. 杭州:浙江科學技術出版社, 2000.

[17] 邢永健,王 旭,杜 航.集對分析在區域大氣環境風險評價中的應用研究 [J]. 中國環境科學, 2016,33(2):634-640.

[18] 胡二邦.環境風險評價實用技術、方法和案例 [M]. 北京:中國環境科學出版社, 2009.

Method of regional acute environmental risk assessment based on risk field.

XING Yong-jian1, WANG Xu1*, KE Xin2, WU Tian-hui1, TENG Wen-chao2(1.School of Safety Engineering, Shenyang Aerospace University, Shenyang 110136, China；2.College of Energy and Environment, Shenyang Aerospace University, Shenyang 110136, China). China Environmental Science, 2016,36(4)：1268～1274

Abstract：In large geographic area, multiple sources release multiple stressors and affect multiple habits and endpoints. The diversity and complexity of the interaction between sources, stressors and habits is one of the challenges of Regional Environmental Risk Assessment (RERA) for large geographic areas. Meanwhile, the existed assessment methods do not pay enough attention on the release rules of risk factors, the distribution patterns of risk factors in space, and the path and the degree of the damage on receptors caused by risk factors. The theory of Risk Field (RF) was applied to analyze the generation mechanism of risk field and the influence mechanism on risk receptors. The environmental risk of Nanjing Chemical Industry Park was investigated as the example. At first, the environmental risk sources were identified. Then, methods, like Set Pair Analysis, were used to construct environmental risk field of various risk sources. After analyzing the environmental risk receptors in the risk field, the distribution of regional environmental risk level could be obtained. The densely populated areas or ecologically sensitive areas located around the risk sources or the river downstream had higher environmental risk level (environmental risk value R≥6), which was in good agreement with the actual situation. This paper suggested establishing the prevention-oriented environmental risk management system based on the classification results of the environmental risks.

Key words：risk field；regional environmental risk；risk assessment；risk receptor；Set Pair Analysis

作者簡介：邢永健(1989-),男,江蘇徐州人,碩士研究生,主要從事系統安全理論與應用、環境風險評價與管理研究.

基金項目：國家水體污染控制與治理科技重大專項(2012ZX07202-004-02)

收稿日期：2015-07-07

中圖分類號：X820.4

文獻標識碼：A

文章編號：1000-6923(2016)04-1268-07