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

?

基于風(fēng)險場的區(qū)域突發(fā)性環(huán)境風(fēng)險評價方法研究

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

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

關(guān)鍵詞：風(fēng)險場；區(qū)域環(huán)境風(fēng)險；風(fēng)險評價；風(fēng)險受體；集對分析法

* 責(zé)任作者, 教授, wangx1960@126.com

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

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

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

本文探索應(yīng)用區(qū)域環(huán)境風(fēng)險系統(tǒng)理論中的環(huán)境風(fēng)險場理論[11],并對其進行實際應(yīng)用的可操作性優(yōu)化,以期彌補相關(guān)研究的不足.

1　理論分析

1.1 環(huán)境風(fēng)險場

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

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

其中:

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

1.2 環(huán)境風(fēng)險受體

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

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

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

1.3 風(fēng)險計算模型

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

型為

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

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

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

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

2　案例分析

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

2.1 環(huán)境風(fēng)險場的構(gòu)建

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

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

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

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

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

表1　南京化工園區(qū)部分企業(yè)環(huán)境風(fēng)險源信息Table 1　Environmental risk sources information of partial enterprises in Nanjing Chemical Industry Park

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

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

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

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

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

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

表3　企業(yè)生產(chǎn)工藝與環(huán)境風(fēng)險控制水平以及P的修正系數(shù)Table 3　Enterprise production process and environmental risk control levels and correction factors of P

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

2.2 環(huán)境風(fēng)險受體易損性計算

同2.1.1,得到整個區(qū)域環(huán)境風(fēng)險受體易損性指數(shù)矩陣(式(10)),此時矩陣內(nèi)各單元的環(huán)境風(fēng)險受體易損性指數(shù)為0.

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

式中:V、S和A分別為風(fēng)險受體易損性指數(shù)、敏感性指數(shù)和適應(yīng)力指數(shù).

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

根據(jù)表4,通過式(3) (受體做靜態(tài)處理時, Px,yj=0或1)分別計算區(qū)域各單元環(huán)境風(fēng)險受體易損性指數(shù),下面是(17, 26)單元周圍1km各單元環(huán)境風(fēng)險受體易損性指數(shù)矩陣.

2.3 環(huán)境風(fēng)險值計算

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

應(yīng)用ArcGIS進行可視化處理,繪制環(huán)境風(fēng)險地圖(圖1).

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

圖1　南京化工園區(qū)及周邊環(huán)境風(fēng)險地圖Fig.1　Environmental risk map for Nanjing Chemical Industry Park and the surrounding

3　環(huán)境風(fēng)險管理

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

式中:n0、n1、n2、…為行政區(qū)內(nèi)各單元相同風(fēng)險值的數(shù)量.

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

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

表5　南京化工園區(qū)及周邊行政區(qū)的環(huán)境風(fēng)險值Table 5　Environmental risk values for districts of Nanjing Chemical Industry Park and the surrounding

其他街道緊鄰化工園區(qū),人口密集、生態(tài)環(huán)境敏感,應(yīng)根據(jù)區(qū)域環(huán)境風(fēng)險評價結(jié)果優(yōu)化城市規(guī)劃與土地利用規(guī)劃,制定突發(fā)環(huán)境事件應(yīng)急預(yù)案,預(yù)防為主,完善環(huán)境風(fēng)險管理體系.

4　結(jié)語

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

參考文獻：

[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] 趙 肖,郭振仁.基于環(huán)境后果評價的環(huán)境風(fēng)險源分級模型研究[J]. 安全與環(huán)境學(xué)報, 2010,10(2):105-108.

[3] 周 薇,錢 瑜,李林子,等.基于環(huán)境風(fēng)險后果評價的危險化學(xué)品碼頭選址 [J]. 環(huán)境保護科學(xué), 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] 賈 倩,黃 蕾,袁增偉,等.石化企業(yè)突發(fā)環(huán)境風(fēng)險評價與分級方法研究 [J]. 環(huán)境科學(xué)學(xué)報, 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] 薛鵬麗,曾維華.上海市環(huán)境污染事故風(fēng)險受體脆弱性評價 [J].環(huán)境科學(xué)學(xué)報, 2011,31(11):2556-2561.

[11] 畢 軍,楊 潔,李其亮.區(qū)域環(huán)境風(fēng)險分析和管理 [M]. 北京:中國環(huán)境科學(xué)出版社, 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] 劉桂友,徐琳瑜.一種區(qū)域環(huán)境風(fēng)險評價方法——信息擴散法[J]. 環(huán)境科學(xué)學(xué)報, 2007,27(9):1549-1556.

[14] 謝元博,李 巍,郝芳華.基于區(qū)域環(huán)境風(fēng)險評價的產(chǎn)業(yè)布局規(guī)劃優(yōu)化研究 [J]. 中國環(huán)境科學(xué), 2013,33(3):560-568.

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

[16] 趙克勤.集對分析及其初步應(yīng)用 [M]. 杭州:浙江科學(xué)技術(shù)出版社, 2000.

[17] 邢永健,王 旭,杜 航.集對分析在區(qū)域大氣環(huán)境風(fēng)險評價中的應(yīng)用研究 [J]. 中國環(huán)境科學(xué), 2016,33(2):634-640.

[18] 胡二邦.環(huán)境風(fēng)險評價實用技術(shù)、方法和案例 [M]. 北京:中國環(huán)境科學(xué)出版社, 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-),男,江蘇徐州人,碩士研究生,主要從事系統(tǒng)安全理論與應(yīng)用、環(huán)境風(fēng)險評價與管理研究.

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

收稿日期：2015-07-07

中圖分類號：X820.4

文獻標(biāo)識碼：A

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