李國偉,周廷剛，曾 波，沈敬偉，李宗峰

(西南大學(xué) a 三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,b 地理科學(xué)學(xué)院，重慶 400715)

近年來，三峽庫區(qū)水環(huán)境形勢不容樂觀，長江流域突發(fā)性水污染事件也有逐年增加的趨勢。葛洲壩庫區(qū)黃柏河黃磷污染事件、萬州航空油泄漏事件、沱江污染事件、重慶墊江縣英特化工有限公司爆炸事故等，這些突發(fā)性水污染事故的發(fā)生給生態(tài)環(huán)境、人群健康和社會經(jīng)濟(jì)發(fā)展造成極大的損失和影響[1]。事故型水環(huán)境污染最突出的特點(diǎn)就是事發(fā)地點(diǎn)和事發(fā)時間的不確定性。因此，除了需要考慮污染事故對水體本身及生態(tài)系統(tǒng)的風(fēng)險外，還需要關(guān)注污染事故對其他敏感目標(biāo)如水源、取水口的風(fēng)險。

國際上對突發(fā)性水環(huán)境污染事故的研究，最初主要集中在對海域突發(fā)性污染事故的風(fēng)險分析和評估上。Emerson等[2]分析了海域移動風(fēng)險源(如船舶)及海岸線固定風(fēng)險源(如化工廠、港口)的安全狀況，探討了突發(fā)性水污染事件對海岸生態(tài)環(huán)境的影響。美國頒布的《反恐法》明確規(guī)定，水源是美國重要的基礎(chǔ)設(shè)施，該法第5部分對涵蓋美國90%以上人口的8 000多個供水人口在3 300人以上的水源地供水系統(tǒng)的飲用水源保障和安全做出了專門規(guī)定，要求對這些水源供水系統(tǒng)進(jìn)行易損性評價并制定應(yīng)急對策計(jì)劃[3]。法國開發(fā)出一個稱為“SEANS”的軟件包，可為突發(fā)性水污染事故提供應(yīng)急決策；歐盟在萊茵河、多瑙河、易北河先后建立了多瑙河水質(zhì)突發(fā)事故預(yù)警系統(tǒng)(AEWPS-Danube Accident Emergency Warning and Prevention System)，旨在開發(fā)跨界河流污染事故早期信息發(fā)布系統(tǒng)，建立流域國際合作程序，以便相關(guān)國家及時應(yīng)對[4]。

我國最近10年來對水環(huán)境污染的研究更多集中在固定點(diǎn)源、農(nóng)業(yè)面源的污染控制，流域生態(tài)管理，土地利用調(diào)整等方面[5-7]，事故型水環(huán)境污染風(fēng)險問題沒有得到應(yīng)有的重視，相關(guān)研究內(nèi)容相對較少，主要是一些較為單一、分散的專題研究，如污染物質(zhì)在水體中擴(kuò)散行為的模型模擬[8]、突發(fā)性水環(huán)境污染事故后的水質(zhì)健康風(fēng)險評估[9]、河流突發(fā)污染事件的環(huán)境風(fēng)險評價模型和事故泄露行為的模擬分析[10]、水環(huán)境健康風(fēng)險評價模型[11]等。也有學(xué)者從區(qū)域角度建立了水環(huán)境預(yù)警、分析系統(tǒng)，如張羽[12]從水源地安全的角度，針對水源地突發(fā)性水污染事件進(jìn)行了風(fēng)險評價理論及方法的實(shí)證研究；張鈞[13]對江河水源地突發(fā)事故的預(yù)警體系與模型進(jìn)行了研究；呂連宏等[14]對沿江化工園突發(fā)水污染事故的預(yù)警指標(biāo)體系進(jìn)行了研究；袁潔等[15]建立了城市飲用水源地地理信息系統(tǒng)；張艷軍[16]建立了“基于SOA的三峽庫區(qū)水環(huán)境安全預(yù)警平臺集成”。以上研究為保障飲用水源、水污染防治和動態(tài)監(jiān)控等提供了強(qiáng)大的軟硬件環(huán)境和技術(shù)支撐。但總體來看，對事故型水環(huán)境污染風(fēng)險評價的研究還不足，關(guān)于水環(huán)境風(fēng)險源的辨識、分級和標(biāo)準(zhǔn)、定性定量評估等都還缺乏研究，在理論和技術(shù)方法上也有待進(jìn)一步分析。目前對三峽庫區(qū)事故型水環(huán)境風(fēng)險的研究還非常薄弱。

本研究從水環(huán)境污染風(fēng)險源及敏感目標(biāo)信息辨識、水環(huán)境風(fēng)險源及敏感目標(biāo)分級評估、分區(qū)評估等方面進(jìn)行綜合分析，以.NET為開發(fā)平臺，采用瀏覽器/服務(wù)器(Browser/Server，B/S)架構(gòu)，建立三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)，以期為政府相關(guān)部門對三峽庫區(qū)水污染事故的風(fēng)險評估及決策提供支持，也為整個三峽庫區(qū)水環(huán)境的可持續(xù)發(fā)展提供有效的信息參考。

1 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的總體設(shè)計(jì)

1.1 系統(tǒng)的總體架構(gòu)設(shè)計(jì)

以Microsoft Visual Studio 2010 和ArcGIS Server 9.3.1為開發(fā)平臺，以SQL Server2008數(shù)據(jù)庫進(jìn)行數(shù)據(jù)管理與分析，基于C#程序設(shè)計(jì)語言，搭載上富客戶端Silverlight技術(shù)，設(shè)計(jì)三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)。系統(tǒng)運(yùn)用B/S模式的3層體系架構(gòu)，從系統(tǒng)開發(fā)的業(yè)務(wù)流程上看，系統(tǒng)平臺可分為數(shù)據(jù)層、功能層和應(yīng)用層3個層次(圖1)。

1.2 系統(tǒng)的數(shù)據(jù)庫設(shè)計(jì)

三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的數(shù)據(jù)庫主要由空間數(shù)據(jù)庫、業(yè)務(wù)數(shù)據(jù)庫和用戶數(shù)據(jù)庫3部分構(gòu)成。其中空間數(shù)據(jù)庫主要是指三峽庫區(qū)的基礎(chǔ)地理數(shù)據(jù)，主要包括行政區(qū)界、道路、主要居民區(qū)、河流水資源利用和生態(tài)環(huán)境等數(shù)據(jù)；業(yè)務(wù)數(shù)據(jù)庫主要是風(fēng)險源和敏感目標(biāo)等數(shù)據(jù)，比如名稱、經(jīng)緯度、類型、所在區(qū)縣、暴露人口、?；访Q、現(xiàn)存量(t)、臨界量(t)、允許限值(mg/L)、事故發(fā)生概率、化學(xué)品校正系數(shù)、化學(xué)品濃度、跨界影響系數(shù)等；用戶數(shù)據(jù)庫主要是用戶、管理員信息表，賦予管理員最高權(quán)限，對用戶注冊、登錄進(jìn)行控制等。系統(tǒng)的數(shù)據(jù)庫設(shè)計(jì)如圖2所示，為了保證數(shù)據(jù)安全，數(shù)據(jù)庫有嚴(yán)格的權(quán)限管理和數(shù)據(jù)加密策略。

圖 1 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的結(jié)構(gòu)

圖 2 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的數(shù)據(jù)庫設(shè)計(jì)

2 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的功能設(shè)計(jì)

結(jié)合三峽庫區(qū)事故型水環(huán)境污染評估系統(tǒng)的特點(diǎn)，本系統(tǒng)主要由5個功能模塊構(gòu)成，分別是基本操作模塊、數(shù)據(jù)輸出模塊、風(fēng)險評估模塊、信息查詢模塊和系統(tǒng)管理模塊，系統(tǒng)功能結(jié)構(gòu)見圖3。ArcGIS API for Silverlight組件式的開發(fā)，實(shí)現(xiàn)了地圖基本操作、信息查詢和數(shù)據(jù)輸出功能；ArcSDE結(jié)合SQL Server數(shù)據(jù)庫管理，實(shí)現(xiàn)了風(fēng)險評估和系統(tǒng)管理模塊等功能。

圖 3 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的功能結(jié)構(gòu)

2.1 基本操作

基本操作主要是通過ArcGIS API for Silverlight接口實(shí)現(xiàn)的，主要包括地圖切換、地圖縮放、地圖漫游、地圖面積測量、距離測量、全圖顯示以及鷹眼等功能。

2.2 數(shù)據(jù)輸出

數(shù)據(jù)輸出主要以圖片形式導(dǎo)出為主，提供PDF、PNG、JPG、BMP等4種常用的格式。

2.3 風(fēng)險評估

風(fēng)險評估是系統(tǒng)的核心模塊，主要包括風(fēng)險源分級評估、敏感目標(biāo)分級評估、風(fēng)險分區(qū)評估3個功能子模塊。

2.3.1 風(fēng)險源分級評估 系統(tǒng)的風(fēng)險源包括4種，分別為化工廠、污水處理廠、?；反a頭、水上加油站(船)，對這4種風(fēng)險源進(jìn)行分級評估時，先調(diào)用數(shù)據(jù)庫的風(fēng)險源數(shù)據(jù)，計(jì)算出風(fēng)險源的風(fēng)險值R1并保存到數(shù)據(jù)庫，然后對計(jì)算出的風(fēng)險值R1進(jìn)行分級，可以將風(fēng)險源劃分為3個不同的級別：特大風(fēng)險源、重大風(fēng)險源和一般風(fēng)險源，分別以紅、藍(lán)、綠不同顏色的圓形符號在地圖上顯示出來，以便區(qū)別。

2.3.2 敏感目標(biāo)分級評估 確定敏感目標(biāo)時，以容易受到風(fēng)險源的污染威脅為標(biāo)準(zhǔn)，敏感目標(biāo)包括多種類型，其中飲用水源地與人類生產(chǎn)、生活最為密切，故系統(tǒng)中只對飲用水源地進(jìn)行分級評估。具體過程為：先調(diào)用數(shù)據(jù)庫的敏感目標(biāo)數(shù)據(jù)，計(jì)算出敏感目標(biāo)風(fēng)險值R2并保存到數(shù)據(jù)庫，然后根據(jù)計(jì)算出的敏感目標(biāo)風(fēng)險值R2進(jìn)行分級，可以將敏感目標(biāo)劃分為3個不同的級別：特大敏感目標(biāo)、重大敏感目標(biāo)和一般敏感目標(biāo)，分別以紅、藍(lán)、綠不同顏色的五角星符號在地圖上顯示出來，以便區(qū)別。

2.3.3 風(fēng)險分區(qū)評估 根據(jù)風(fēng)險源風(fēng)險值R1、敏感目標(biāo)風(fēng)險值R2的不同，采用分別基于風(fēng)險源的風(fēng)險分區(qū)以及基于風(fēng)險源和敏感目標(biāo)耦合后的風(fēng)險分區(qū)。在對風(fēng)險源和敏感目標(biāo)分級的基礎(chǔ)上，根據(jù)不同的R1和R2值，計(jì)算各種不同的風(fēng)險分區(qū)結(jié)果，然后根據(jù)分區(qū)標(biāo)準(zhǔn)將不同風(fēng)險區(qū)劃分為3類：高風(fēng)險區(qū)、中風(fēng)險區(qū)和低風(fēng)險區(qū)，分別以紅、藍(lán)、綠不同顏色在地圖上顯示出來,以便區(qū)別。

2.4 信息查詢

2.4.1 屬性查圖 1)全部查詢。顯示數(shù)據(jù)庫中所有的風(fēng)險源和敏感目標(biāo)的屬性信息、空間信息。

2)邏輯查詢。根據(jù)用戶設(shè)置的多個關(guān)鍵字段和邏輯謂詞，提取出符合條件的風(fēng)險源和敏感目標(biāo)。

3)范圍查詢。用戶設(shè)置查詢范圍及查詢中心點(diǎn)，搜索距中心點(diǎn)一定范圍內(nèi)的風(fēng)險源和敏感目標(biāo)。

4)直接查詢。即模糊查詢，根據(jù)用戶輸入的關(guān)鍵字段，即可搜索出包含關(guān)鍵字段的風(fēng)險源和敏感目標(biāo)。

2.4.2 圖查屬性 主要包括三角形查詢、多邊形查詢、矩形查詢及圓查詢，均是在地圖上設(shè)置好空間范圍，再查詢范圍內(nèi)包含的風(fēng)險源和敏感目標(biāo)的屬性信息。

2.5 系統(tǒng)管理

系統(tǒng)管理主要是賦予管理員權(quán)限，可以對空間數(shù)據(jù)和屬性數(shù)據(jù)進(jìn)行編輯、查詢、統(tǒng)計(jì)等操作。屬性數(shù)據(jù)存放在SQL Server數(shù)據(jù)庫中。

3 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的關(guān)鍵技術(shù)研究

3.1 風(fēng)險源分級方法

不同風(fēng)險源的污染物差異很大，其中化工廠的污染物以鹽酸、乙酸、氫氧化鈉等為主，水上加油站(船)的污染物以汽油、柴油為主，?；反a頭的污染物以汽油、柴油、瀝青、潤滑油為主，污水處理廠的污染物以生活污水為主。參照《危險化學(xué)品重大危險源辨識標(biāo)準(zhǔn)GB18218-2009》，可知不同危險化學(xué)品有不同的臨界量，利用臨界量值來判斷風(fēng)險級別的大小，采用單元內(nèi)經(jīng)校正系數(shù)校正后的各種危險化學(xué)品的實(shí)際存在量與其臨界量的比值(L)對指標(biāo)分級[17]，計(jì)算公式如下：

L=?1q1/Q1+?2q2/Q2+…+?nqn/Qn。

(1)

式中：q1，q2，…，qn表示每種危險化學(xué)品的實(shí)際存在量，單位為t；Q1，Q2，…，Qn表示與各危險化學(xué)品相對應(yīng)的臨界量，單位為t；?1，?2，…，?n表示與各危險化學(xué)品相對應(yīng)的校正系數(shù)。

除了考慮污染物的臨界量以及校正系數(shù)之外，還應(yīng)考慮其他影響因素。根據(jù)風(fēng)險源的性質(zhì)，對飲用水源地的污染風(fēng)險評價，參照《生活飲用水衛(wèi)生標(biāo)準(zhǔn)GB5749-2006》，可知不同污染物質(zhì)有不同的允許限值B。此外，風(fēng)險源突發(fā)風(fēng)險事故發(fā)生概率P以及化學(xué)物品的濃度C[18]也是主要影響因素，其中P的計(jì)算公式如下：

P=P平均+Ks+Ka+Kp+Km。

(2)

式中：P平均為行業(yè)平均事故概率；Ks為風(fēng)險源類型影響系數(shù)，Ks=Ps/P平均，Ps為不同類型風(fēng)險源事故概率，有固定源及移動源；Ka為受體影響系數(shù)，Ka=Pa/P平均，Pa為不同受體事故概率；Kp為風(fēng)險事故場所影響系數(shù)，Kp=Pp/P平均，Pp為不同事故場所事故概率；Km為管理水平影響系數(shù)，Km=Pm/P平均，Pm為不同管理水平事故概率。

結(jié)合式(1)、(2)及提到的影響因素，可計(jì)算一個單位風(fēng)險源的污染風(fēng)險值R1，計(jì)算公式見式(3)：

(3)

風(fēng)險源等級劃分是在風(fēng)險值R1的基礎(chǔ)上通過統(tǒng)計(jì)計(jì)算得來的，分別確定各因素的平均值(μ)和標(biāo)準(zhǔn)偏差(σ)[19-20]。設(shè)定風(fēng)險值R1≥μ+1.15σ時為特大風(fēng)險源，在μ-0.05σ≤R1<μ+1.15σ時為重大風(fēng)險源，在R1<μ-0.05σ時為一般風(fēng)險源。風(fēng)險源定性等級及相應(yīng)的風(fēng)險值范圍如表1所示。

表 1 風(fēng)險源定性等級與相應(yīng)的風(fēng)險值范圍

3.2 敏感目標(biāo)分級方法

本研究中，敏感目標(biāo)只考慮飲用水源地，參考《飲用水源地保護(hù)區(qū)劃分技術(shù)規(guī)范》(HJ/T338-2007)劃定，上游河段1 000 m和下游100 m以中泓線為界的同側(cè)水域?yàn)橐患壉Ｗo(hù)區(qū)，二級保護(hù)區(qū)的上游邊界為一級保護(hù)區(qū)上游邊界向上游延伸不得小于2 000 m，二級保護(hù)區(qū)的下游邊界距一級保護(hù)區(qū)下游邊界的距離不得小于200 m，這里飲用水源風(fēng)險評估范圍以風(fēng)險源處于一級保護(hù)區(qū)和二級保護(hù)區(qū)為準(zhǔn)。

飲用水源的風(fēng)險評估除了考慮評估范圍內(nèi)不同污染物的毒性(毒性大小以臨界量Q和允許限值B來判斷)之外，以《國家突發(fā)環(huán)境事件應(yīng)急預(yù)案》為依據(jù)，還應(yīng)考慮污染物擴(kuò)散時間T、暴露人口M、跨界影響系數(shù)K(不跨界取1，跨縣界取2，跨省界取 2.5，跨多個界需要疊加跨界系數(shù))等其他因素。飲用水源地風(fēng)險值R2的計(jì)算公式如下：

(4)

同理，敏感目標(biāo)等級劃分以R2≥μ+1.15σ時為特大敏感目標(biāo)，在μ-0.05σ≤R2<μ+1.15σ時為重大敏感目標(biāo)，在R2<μ-0.05σ時為一般敏感目標(biāo)。敏感目標(biāo)定性等級及相應(yīng)的風(fēng)險值范圍見表2。

表 2 敏感目標(biāo)定性等級與相應(yīng)的風(fēng)險值范圍

3.3 風(fēng)險分區(qū)方法

在本研究范圍內(nèi)，風(fēng)險分區(qū)主要是以長江主干流以及重要的支流為主，以10 km為一個單位區(qū)段，值得注意的是，河道范圍越小則分析精度越高，比如按每5 km為一個單位區(qū)段計(jì)算，精度比按每10 km為一個單位區(qū)段要高。

在確定了風(fēng)險源風(fēng)險值R1、敏感目標(biāo)風(fēng)險值R2的基礎(chǔ)上，利用ArcGIS API提供的Dijkstra最短路徑算法，將不同級別的風(fēng)險源和敏感目標(biāo)歸并到河流上不同區(qū)段內(nèi)，然后統(tǒng)計(jì)出各個不同河流區(qū)段內(nèi)不同級別風(fēng)險源或敏感目標(biāo)的個數(shù)，分別累加求和(∑R)。最后通過ArcGIS API提供的自然斷點(diǎn)法(Natural Breaks)分類出不同的拐點(diǎn)值來定義分區(qū)標(biāo)準(zhǔn)，以分區(qū)標(biāo)準(zhǔn)為依據(jù)，將風(fēng)險區(qū)劃分為3種類型，分別為高風(fēng)險區(qū)、中風(fēng)險區(qū)和低風(fēng)險區(qū)。

4 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)的應(yīng)用實(shí)例

本研究以三峽庫區(qū)地理數(shù)據(jù)、風(fēng)險源和敏感目標(biāo)等為試驗(yàn)數(shù)據(jù)，構(gòu)建了三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)(圖4)。該系統(tǒng)對研究區(qū)域內(nèi)的風(fēng)險源及敏感目標(biāo)的分級情況進(jìn)行了實(shí)例化(圖5、圖6)；同時亦對基于風(fēng)險源以及基于風(fēng)險源對敏感目標(biāo)的影響整合后的風(fēng)險分區(qū)情況進(jìn)行了分析及可視化(圖7和圖8)，從而可以較為直觀地反映研究區(qū)域內(nèi)的風(fēng)險分布情況。

圖 4 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)主界面

從圖7、圖8可以看出，2種風(fēng)險分區(qū)方法得到的結(jié)果差異明顯。以局部放大圖為例，圖7主城區(qū)中，由于重慶江北區(qū)、南岸區(qū)是風(fēng)險源的密集區(qū)，而且特大風(fēng)險源、重大風(fēng)險源比較多，所以計(jì)算結(jié)果顯示該地區(qū)為高風(fēng)險區(qū)；而渝北區(qū)嘉陵江支流入口風(fēng)險源數(shù)量不多，但是有特大風(fēng)險源而且風(fēng)險值特大，另外大渡口區(qū)風(fēng)險源等級大多是一般風(fēng)險源以及重大風(fēng)險源，所以計(jì)算結(jié)果顯示該地區(qū)為中風(fēng)險區(qū)。圖8中，大渡口區(qū)以及渝北區(qū)嘉陵江支流入口的敏感目標(biāo)比較多，而且大多是特大敏感目標(biāo)、重大敏感目標(biāo)，該地區(qū)附近以及上游的風(fēng)險源也比較多，計(jì)算結(jié)果顯示該地區(qū)為高風(fēng)險區(qū)；南岸區(qū)的敏感目標(biāo)雖然數(shù)量比較少，但是有重大敏感目標(biāo)且附近以及上游的風(fēng)險源比較多，所以計(jì)算結(jié)果顯示該地區(qū)為中風(fēng)險區(qū)。而江北區(qū)敏感目標(biāo)數(shù)量比較少，并且上游的風(fēng)險源距離江北區(qū)太遠(yuǎn)影響較小，所以計(jì)算結(jié)果顯示該地區(qū)為低風(fēng)險區(qū)。

圖 5 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)中基于風(fēng)險源的風(fēng)險分級

圖 7 三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)中基于風(fēng)險源的風(fēng)險分區(qū)

5 小 結(jié)

本研究在對系統(tǒng)總體設(shè)計(jì)、功能設(shè)計(jì)和關(guān)鍵技術(shù)進(jìn)行分析的基礎(chǔ)上，建立了三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)，該系統(tǒng)實(shí)現(xiàn)了基本的GIS功能及可視化的風(fēng)險評估效果功能，可反映研究區(qū)域內(nèi)的風(fēng)險分布情況，從而可為政府部門在三峽庫區(qū)水環(huán)境污染風(fēng)險評估方面提供較為科學(xué)的決策支持，也可為三峽庫區(qū)水環(huán)境可持續(xù)發(fā)展提供較為有效的信息參考。由于三峽庫區(qū)事故型水環(huán)境污染風(fēng)險評估系統(tǒng)涉及到海量數(shù)據(jù)，而敏感目標(biāo)分級評估子模塊中關(guān)于整合風(fēng)險源的影響后敏感目標(biāo)的價值用威脅度計(jì)算的算法效率較低，因此在風(fēng)險源、敏感目標(biāo)分級評估及風(fēng)險分區(qū)評估的數(shù)學(xué)模型驗(yàn)證方面，還需要在實(shí)際應(yīng)用中不斷完善與優(yōu)化。此外，對水環(huán)境風(fēng)險污染進(jìn)行預(yù)警，也是下一步的研究重點(diǎn)。

[參考文獻(xiàn)]

[1] 陳蓓青，譚德寶，宋 麗．GIS技術(shù)在突發(fā)性水污染事件應(yīng)急響應(yīng)系統(tǒng)中的應(yīng)用研究 [J].長江科學(xué)院院報，2010，27(1)：29-32．

Chen B Q,Tan D B,Song L.Application research of GIS technology in suddenly-happened water pollution accident emergency response system [J].Journal of Yangtze River Scientific Research Institute,2010, 27(1): 29-32．(in Chinese)

[2] Emerson S D， Nadeau J．A coastal perspective on security [J].Journal of Hazardous Materials,2003,104(1/2/3):1-13．

[3] 張 勇，徐啟新，楊 凱，等．城市水源地突發(fā)性水污染事件研究述評 [J].環(huán)境污染治理技術(shù)與設(shè)備，2006，7(12)：1-4．

Zhang Y，Xu Q X，Yang K，et al．Review of studies on water pollution incidents in areas of urban water supply [J].Techniques and Equipment for Environmental Pollution Control，2006，7(12)：1-4．(in Chinese)

[4] Zografos K G，Vasilakis G M，Giannouli I M．Methodological framework developing decision support system(DSS) for hazardous materials emergency response operations [J].Journal of Hazardous Materials,2000,71:503-521．

[5] 楊愛玲，朱顏明．城市地表飲用水源保護(hù)研究進(jìn)展 [J].地理科學(xué)，2000，20(1)：72-77．

Yang A L，Zhu Y M．The study of urban surface drinking water sources protection [J].Scientia Geographica Sinica，2000，20(1)：72-77．(in Chinese)

[6] 鄭天柱，倪深海，熊遠(yuǎn)明，等．城市水源污染的預(yù)防和生態(tài)治理綜述 [J].水資源保護(hù)，2001，66(4)：16-19．

Zheng T Z，Ni S H，Xiong Y M，et al．Prevention of water source pollution and ecological environment remediation in cities [J].Water Resources Protection，2001，66(4)：16-19．(in Chinese)

[7] 曾維華，程聲通，楊志峰．流域水資源集成管理 [J].中國環(huán)境科學(xué)，2001，21(2)：173-176．

Zeng W H，Cheng S T，Yang Z F．Integrated water resources management on river basin [J].China Environmental Science，2001，21(2)：173-176．(in Chinese)

[8] 石劍榮．水體擴(kuò)散衍生公式在環(huán)境風(fēng)險評價中的應(yīng)用 [J].水科學(xué)進(jìn)展，2005，16(1)：92-102．

Shi J R．Application of extended waterbody diffusion formulate to environmental risk assessment [J].Advances in Water Science，2005，16(1)：92-102．(in Chinese)

[9] 錢家忠，李如忠，汪家權(quán)，等．城市供水水源地水質(zhì)健康風(fēng)險評價 [J].水利學(xué)報，2004，24(8)：90-93．

Qian J Z，Li R Z，Wang J Q，et al．Environmental health risk assessment for urban water supply source [J].Journal of Hydraulic Engineering，2004，24(8)：90-93．(in Chinese)

[10] 曾光明，卓 利，鐘政林，等．突發(fā)性水環(huán)境風(fēng)險評價模型事故泄漏行為的模擬分析 [J].中國環(huán)境科學(xué)，1998，18(5)：403-406．

Zeng G M，Zhuo L，Zhong Z L，et al．Assessment model about the accident water environment risk modeling of the toxicant accident leakage [J].China Environmental Science，1998，18(5)：403-406．(in Chinese)

[11] 曾光明，卓 利，鐘政林，等．水環(huán)境健康風(fēng)險評價模型 [J].水科學(xué)進(jìn)展，1998，9(3)：213-217．

Zeng G M，Zhuo L，Zhong Z L，et al．Assessment models for water environmental health risk analysis [J].Advances in Water Science，1998，9(3)：213-217．(in Chinese)

[12] 張 羽．城市水源地突發(fā)性水污染事件風(fēng)險評價體系及方法的實(shí)證研究 [D].上海：華東師范大學(xué)，2006．

Zhang Y．Empirical studies on risk assessment system of accidental water pollution events in city’s source water area [D].Shanghai：East China Normal University，2006．(in Chinese)

[13] 張 鈞．江河水源地突發(fā)事故預(yù)警體系與模型研究 [D].南京：河海大學(xué)，2007．

Zhang J．Research on the early-warning system and models for abrupt accidents in water source areas of rivers [D].Nanjing：HoHai University，2007．(in Chinese)

[14] 呂連宏，羅 宏，路超君．沿江化工園突發(fā)水污染事故預(yù)警指標(biāo)體系研究 [J].工業(yè)安全與環(huán)保，2011，37(1)：30-32．

Lü L H，Luo H，Lu C J．Research on early warning system for sudden water pollution accident in chemical industry park along the Yangtze river [J].Industrial Safety and Environmental Protection，2011，37(1)：30-32．(in Chinese)

[15] 袁 潔，張子凡．南京市飲用水源地環(huán)境管理地理信息系統(tǒng)的建立 [J].污染防治技術(shù)，2005，18(3)：23-25．

Yuan J，Zhang Z F．Development of geographic information system for environmental management of drinking water resource in Nanjing [J].Pollution Control Technology，2005，18(3)：23-25．(in Chinese)

[16] 張艷軍．基于SOA的三峽庫區(qū)水環(huán)境安全預(yù)警平臺集成 [J].四川環(huán)境，2010，29(1)：48-50．

Zhang Y J．A safety precaution system of water environment for the area of Three Gorges Reservoir based on SOA [J].Sichuan Environment，2010，29(1)：48-50．(in Chinese)

[17] 魏利軍，吳宗之．一種簡單可行的危險化學(xué)品重大危險源分級方法 [J].中國安全科學(xué)學(xué)報，2009,12(12):45-50．

Wei L J，Wu Z Z．A simple and feasible grading method for major dangerous chemical hazard installations [J].China Safety Science Journal，2009,12(12):45-50．(in Chinese)

[18] 馬 越，彭劍峰，宋永會，等．飲用水源地突發(fā)事故環(huán)境風(fēng)險分級方法研究 [J].環(huán)境科學(xué)學(xué)報，2012,32(5):1211-1218．

Ma Y，Peng J F，Song Y H，et al．Study on the classification method of environmental risk for sudden accidents in drinking water sources [J].Acta Scientiae Circumstantiae，2012,32(5):1211-1218．(in Chinese)

[19] Colorado Department of Public Health and Environment Water Quality Control Division.Source water assessment methodology for surface water sources and ground water sources under the direct influence of surface water [R].Colorado:EnecoTech Inc,2004.

[20] 劉 琰，鄭丙輝，萬 峻，等．C市飲用水源風(fēng)險評價實(shí)例分析 [J].環(huán)境科學(xué)研究，2009,22(1):52-59．

Liu Y，Zheng B H，Wan J，et al．Risk assessment for drinking water sources of C City:A case study [J].Research of Environmental Sciences，2009,22(1):52-59．(in Chinese)