兩階段評估體系篩選水源突發(fā)污染應(yīng)急最優(yōu)技術(shù)方案

曲建華，孟憲林，2*，尤 宏，2（.哈爾濱工業(yè)大學市政環(huán)境工程學院，黑龍江 哈爾濱 50090；2.哈爾濱工業(yè)大學城市水資源與水環(huán)境國家重點實驗室，黑龍江 哈爾濱 50090）

兩階段評估體系篩選水源突發(fā)污染應(yīng)急最優(yōu)技術(shù)方案

曲建華1，孟憲林1，2*，尤 宏1，2（1.哈爾濱工業(yè)大學市政環(huán)境工程學院，黑龍江 哈爾濱 150090；2.哈爾濱工業(yè)大學城市水資源與水環(huán)境國家重點實驗室，黑龍江 哈爾濱 150090）

為預判突發(fā)污染對水源的威脅程度，進而篩選適宜污染情景的最優(yōu)應(yīng)急處置技術(shù)，本研究構(gòu)建了包含10個評價指標的水源突發(fā)污染威脅度判別指標體系，依據(jù)威脅度等級判定，確定與不同威脅度對應(yīng)的應(yīng)急處置技術(shù)評估指標體系的指標權(quán)重.同時，為合理體現(xiàn)群決策過程中專家差異對權(quán)重的影響，提出了基于循環(huán)修正的群組G1指標賦權(quán)法.以Spearman等級相關(guān)系數(shù)作為檢驗標準，應(yīng)用平均值法、Board法、Copeland法3種組合評價方法對單一專家評價結(jié)果進行循環(huán)迭代修正，最終得到評價指標的一致性排序和各位專家的權(quán)重，使各項指標的權(quán)重賦值更為準確.將綜合評估模型應(yīng)用于2012年廣西龍江鎘污染事故中，成功篩選出“調(diào)水稀釋—絮凝沉降—水廠強化混凝”的組合應(yīng)急技術(shù)方案，同實際情況吻合，驗證了評估模型的可行性.

威脅度判定；飲用水源地；循環(huán)修正；群組G1；應(yīng)急處置技術(shù)

水源地因關(guān)系到飲用水水質(zhì)安全問題，故其具有社會敏感性，一旦遭受污染，應(yīng)急專家需在短時間內(nèi)篩選出最優(yōu)的應(yīng)急處置技術(shù)方案對污染進行有效處置［1］.目前，國內(nèi)外針對應(yīng)急技術(shù)篩選方面的研究處于起步階段，相關(guān)研究主要集中于通過AHP［2-4］、FAHP［5-6］、G1［7-8］、TOPSIS［9-10］等方法對單一評價指標體系進行權(quán)重賦值，進而對研究對象進行評估.然而，突發(fā)污染對水源地威脅程度的不同勢必導致評估體系中經(jīng)濟型指標和技術(shù)時效型指標權(quán)重比重的不同，因此，固定一套指標體系權(quán)重進行應(yīng)急技術(shù)評估缺乏科學性.同時，應(yīng)急技術(shù)的篩選與評估是多專家群決策過程，上述方法均未合理體現(xiàn)專家差異對指標權(quán)重的影響，無法確保評估結(jié)果的準確性.

為此，本文構(gòu)建了基于威脅度判定的應(yīng)急處置技術(shù)評估模型.依據(jù)威脅度分級標準，采用基于循環(huán)修正的群組G1法確定了與之對應(yīng)的應(yīng)急技術(shù)評估體系的指標權(quán)重，旨在實現(xiàn)多重指標體系對應(yīng)急處置技術(shù)的聯(lián)合評估.

1 水源突發(fā)污染威脅度等級判定

1.1 威脅度評估指標體系的構(gòu)建

為確定影響突發(fā)污染對水源危害程度的因素，對1985～2013年間近300起水源突發(fā)污染案例進行了調(diào)研，依據(jù)調(diào)研結(jié)果從污染情況、水源地規(guī)模、應(yīng)急管理和備用水源等4方面構(gòu)建了包含10個關(guān)鍵評價指標的水源突發(fā)污染威脅度評估指標體系，如圖1所示.

圖1 水源地突發(fā)污染威脅度評估指標體系Fig.1 Threat evaluation index system for drinking water source pollution

1.2 基于循環(huán)修正的群組G1指標賦權(quán)

1.2．1 單專家的G1主觀賦權(quán) （1）邀請專家根據(jù)威脅度評估指標體系中B1～B4指標的相對重要性程度，確定指標間的排序關(guān)系；（2）根據(jù)專家給出的相鄰指標xi-1與xi重要性程度之比rk的理性賦值，計算第p個指標的G1法權(quán)重wp.計算公式為：

（3）依據(jù)權(quán)重wp的值，計算第p-1，p-2，…，3，2個指標的權(quán)重：

應(yīng)用上述方法獲得了各位專家對B1～B4指標的G1賦權(quán)矩陣：

1.2．2 多專家指標排序的一致性檢驗 為統(tǒng)一多位專家意見、合理體現(xiàn)群體決策思想，以Spearman等級相關(guān)系數(shù)作為檢驗標準［11］，剔除未通過一致性檢驗的專家排序.

設(shè)第j位專家給出的排序為Aj=（a1j，a2j，…，apj），則第j和第---k位專家排序的斯皮爾曼系數(shù)βjk為：

通過計算得到各專家排序的βjk矩陣及針對指標間相對重要程度的判定是群決策過程，需反映多數(shù)人意見；即針對某一排序結(jié)果應(yīng)得到至少半數(shù)專家的認可，而針對個別專家的極端指標排序結(jié)果應(yīng)進行刪除.因此，設(shè)定≥0.50時，專家排序通過檢驗，否則剔除該專家排序.

1.2．3 指標排序的二次排序 為逼近指標權(quán)重的理想值，采用平均值法、Board法、Copealand法對通過Spearman系數(shù)檢驗的專家指標排序進行二次排序.

（1）平均值法：將通過一致性檢驗的各位專家的指標排序結(jié)果轉(zhuǎn)換成分數(shù)：

式中：Rkj為第k個指標（k=1，2，…，n）在第j位專家（j=1，2，…，m）排序下的分數(shù)；rkj為指標k在專家j排序中的排名.其次，計算指標k在多專家評價下的綜合評價值：

按各指標綜合評價值的大小對評價指標重新排序，若兩指標評價值相等，則通過計算不同得分的標準差來確定排序，標準差小者為優(yōu).

（2）Board法：通過一致性檢驗的群組專家中，若認為指標a優(yōu)于指標b的人數(shù)大于認為指標b優(yōu)于指標a的人數(shù)，則記為xaSxb.定義Board矩陣B=｛qab｝n×n，其中：

則指標a的組合得分為：

依據(jù)Ba分值大小對指標重新排序，若2個指標的得分相同，則采用公式（6）確定排序先后.

（3） Copealand法：為在Board法的基礎(chǔ)上區(qū)分“相等”和“劣”，以期得到更好的比較結(jié)果，定義Copealand矩陣C=｛cab｝n×n，其中：

則指標a的組合得分為：

表1 指標（B1～B4）的專家一致性排序Table 1 Consistency in group ordering of B1～B4

對均值法、Board法、Copealand法得到的指標排序結(jié)果進行Spearman一致性檢驗，結(jié)果如表1所示.若檢驗結(jié)果一致，則得到指標一致性排序；否則，需對排序結(jié)果進行2次迭代修正，直到得到一致性排序為止.

1.2．4 專家權(quán)重計算 記專家一致性排序為A0，通過一致性檢驗的m位專家的指標排序分別為A1，A2，…，Am；應(yīng)用公式（3）計算m位專家排序和A0的Spearman等級相關(guān)系數(shù)，記為：β01，β02，…，β0m.則第j位專家的權(quán)重為：

指標k的綜合權(quán)重Wk為：

表2 威脅度評估 體系的指標權(quán)重Table 2 Weights of indexes from threat evaluation system

1.3 指標評分標準的確立

在威脅度判定過程中，依據(jù)各指標對水源地的威脅程度，將評價因子劃分為5個等級，并賦予相應(yīng)分值，見表3.

表3 評價因子的評分標準Table 3 Grading for assessment factors

1.4 威脅度分級標準的確定

依據(jù)表2中指標的綜合權(quán)重（Wk）和表3中指標評分標準，采用線性加權(quán)確定突發(fā)污染威脅度的評判結(jié)果：

式中：Wk為指標層各指標的綜合權(quán)重；Vk為指標層各指標對應(yīng)的評分值.以威脅度評判結(jié)果的最大值與最小值為界，均分為4個區(qū)間，定義為威脅度的4個級別，如表4所示.

2 基于威脅度判定的應(yīng)急處置技術(shù)評估模型

2.1 構(gòu)建應(yīng)急技術(shù)評估指標體系

通過歷史案例的調(diào)研發(fā)現(xiàn)［15-16］，突發(fā)污染應(yīng)急處置技術(shù)的可行度主要與技術(shù)的時效性、可操作性、處置成本等因素有關(guān)；基于上述考慮，建立了包含6個準則指標和19個評價指標的水源污染應(yīng)急處置技術(shù)評估指標體系，如圖2所示.

2.2 計算指標權(quán)重

當突發(fā)污染威脅度為I/II級時，表明城市供水安全受到極大威脅，此時應(yīng)急技術(shù)評估的核心指標應(yīng)為技術(shù)時效性；III/IV級威脅度表明突發(fā)污染對水源地破壞程度有限，不會對供水安全造成嚴重影響，此時在考慮技術(shù)性能的同時，應(yīng)盡可能減少應(yīng)急處置的成本.

為在不同威脅度下對應(yīng)急處置技術(shù)進行準確評估，本文采用循環(huán)修正的群組G1賦權(quán)法計算得到了2套與威脅度等級對應(yīng)的應(yīng)急技術(shù)評估指標體系的權(quán)重值，如表5所示.

表4 威脅度分級標準Table 4 Criterion of threat classification

圖2 水源污染應(yīng)急處置技術(shù)評估指標體系Fig.2 Technology evaluation index system for water source pollutions disposal

表5 指標（C1～C19）對目標層（A1）的權(quán)重Table 5 Weights of indexes C1～C19

續(xù)表5

2.3 確定指標評分標準

表6 應(yīng)急技術(shù)評估體系事中控制指標評分標準Table 6 Scoring criteria of the mid-control indexes on emergency treatment technology

表7 應(yīng)急技術(shù)評估體系事后處置指標評分標準Table 7 Scoring criteria of the post-disposal indexes on emergency treatment technology

本文構(gòu)建的應(yīng)急技術(shù)評估指標體系綜合考慮了突發(fā)污染的事中控制和事后處置.事中控制指標突出強調(diào)了應(yīng)急技術(shù)性能和處置成本，各指標評分標準如表6所示.

事后處置指標著重考慮了應(yīng)急技術(shù)實施及應(yīng)急廢棄物的處置對周圍生態(tài)環(huán)境、居民社會生活所產(chǎn)生的影響，指標評分標準如表7所示.

3 案例應(yīng)用

本文以廣西龍江鎘污染事故為例對綜合評估模型的可行性進行了驗證.

2012年1月15日，由于企業(yè)違法排污導致約21t含鎘廢水通過地下溶洞進入龍江，對下游飲用水源地柳江造成極大威脅，嚴重危及柳州市的供水安全.實時監(jiān)測數(shù)據(jù)顯示龍江中鎘濃度峰值為0.408mg/L，超出《地表水環(huán)境質(zhì)量標準》III類水體中鎘濃度標準限值約83倍［17-18］.

根據(jù)鎘的化學性質(zhì)及相似的歷史應(yīng)急處置案例，篩選出了5套可行的應(yīng)急處置技術(shù)方案，如表8所示［19-20］.來自水體污染防治領(lǐng)域的20位專業(yè)人士對污染事故進行了詳細調(diào)研，采用基于威脅度判定的應(yīng)急處置技術(shù)評估模型對可行的應(yīng)急技術(shù)方案進行評估.依據(jù)實際污染情況對威脅度評估指標進行打分，判定事故為II級威脅（6.2566），因此選用與II級威脅度對應(yīng)的應(yīng)急技術(shù)評估指標體系的權(quán)重值對各應(yīng)急技術(shù)方案進行評估，綜合各位專家的打分結(jié)果，得出的方案評價值及優(yōu)選順序為：方案1（9.0350）＞方案3（8.7342）＞方案4（8.7081）＞方案2（8.6995）＞方案5（8.5364）.

表8 可行的應(yīng)急處置技術(shù)方案Table 8 The feasible technology schemes

圖3 龍江鎘污染應(yīng)急處置技術(shù)方案概況Fig.3 The sketch map of the emergent disposal plan for the cadmium pollution accident

實際應(yīng)急處置過程為：在龍江河沿岸的拉浪、葉茂、洛東、三岔、糯米灘等5級梯級電站處，投加液堿（將河水的pH值提升至8.1～8.4）和PAC，盡可能地將河流污染團中的溶解性鎘離子沉降至河底，并控制龍江河污染水團的下泄流量；同時，在龍江與融江交匯處，采用調(diào)水稀釋技術(shù)，即通過加大浮石、古頂、大浦等電站的放流量，進一步稀釋污染團的濃度；最后，對柳北、城中等5個水廠進行應(yīng)急改造，采用強化混凝技術(shù)對污染物進行去除，見圖3.通過多重應(yīng)急技術(shù)的聯(lián)合處置，鎘污染團的濃度降至0.0039mg/L，滿足供水要求.案例應(yīng)用結(jié)果表明，本研究構(gòu)建的應(yīng)急技術(shù)綜合評估模型成功篩選出了與實際案例相符的應(yīng)急處置技術(shù)方案.

4 結(jié)論

4.1 提出了突發(fā)水源污染威脅度判定理論，確定了與不同威脅度等級對應(yīng)的應(yīng)急技術(shù)評估指標體系的權(quán)重值；通過多重指標體系對應(yīng)急處置技術(shù)的聯(lián)合評估，確保了評估模型的科學性.

4.2 應(yīng)用基于循環(huán)修正的群組G1法對評估指標進行賦權(quán).通過Spearman等級相關(guān)系數(shù)對未通過一致性檢驗的專家排序的剔除，確保了多專家主觀排序的一致性；通過對指標一致性排序和單專家排序相關(guān)系數(shù)的循環(huán)修正計算得到了單體專家權(quán)重，合理體現(xiàn)了專家差異對權(quán)重的影響. 4.3 通過廣西龍江鎘污染事故實例對評估模型進行了驗證.結(jié)果表明，基于威脅度判定的應(yīng)急技術(shù)綜合評估模型可迅速、準確地篩選出適宜污染情景的最優(yōu)應(yīng)急處置技術(shù)方案，為應(yīng)急決策者提供有效的技術(shù)支撐.

［1］張曉健，陳 超.應(yīng)對突發(fā)性水源污染的城市應(yīng)急供水的進展與展望 ［J］. 城鎮(zhèn)給排水， 2011，37（10）：9-18.

［2］Tsita K G， Pilavachi P A. Evaluation of alternative fuels for the Greek road transport sector using the analytic hierarchy process［J］. Energy Policy， 2012，（48）：677-686.

［3］杜大仲，孟憲林，馬 放.北方某城市河流型飲用水水源地選址方案評價研究 ［J］. 中國環(huán)境科學， 2012，32（2）：359-365.

［4］SHI Shenggang， CAO Jingcan， LI Feng. Construction of a technique plan repository and evaluation system based on AHP group decision-making for emergency treatment and disposal in chemical pollution accidents ［J］. Journal of Hazardous Materials， 2014，（276）：200-206.

［5］Shaverdi M， Heshmati M R， Ramezani I. Application of Fuzzy AHP Approach for Financial Performance Evaluation of Iranian Petrochemical Sector ［J］. Procedia Computer Science， 2014，（31）：995-1004.

［6］Heo E， Kim J， Cho S. Selecting hydrogen production methods using fuzzy analytic hierarchy process with opportunities costs and risks ［J］Internationaln journal of Hydrogen Energy， 2012，（37）： 17655-17662.

［7］王恩旭.基于G1—熵值的智慧旅游城市建設(shè)水平評價模型及實證研究 ［J］. 大連理工大學學報， 2014，35（2）：69-73.

［8］劉仁濤，郭 亮，姜繼平，等.環(huán)境污染應(yīng)急處置技術(shù)的CBRMADM兩步篩選法模型 ［J］. 中國環(huán)境科學， 2015，35（3）：943-952.

［9］王紅旗，秦 成，陳美陽.地下水水源地污染防治優(yōu)先性研究 ［J］.中國環(huán)境科學， 2011，31（5）：876-880.

［10］Jayakumar D N， Venkatesh P. Glowworm swarm optimization algorithm with topsis for solving multiple objective environmental economic dispatch problem ［J］. Applied Soft Computing， 2014，（23）：375-386.

［11］李浩鑫，邵東國，何思聰，等.基于循環(huán)修正的灌溉用水效率綜合評價方法 ［J］. 農(nóng)業(yè)工程學報， 2014，30（5）：65-72.

［12］周文敏，傅德黔，孫宗光.水中優(yōu)先控制污染物黑名單 ［J］. 中國環(huán)境監(jiān)測， 1990，6（4）：1-4.

［13］張 悅，張曉健，陳 超.城市供水系統(tǒng)應(yīng)急凈水技術(shù)指導手冊［M］. 北京：中國建筑工業(yè)出版社， 2009.

［14］陸艷晨，李翠梅，劉成剛.應(yīng)急狀態(tài)下的城市居民生活供水量與用水量分析 ［J］. 凈水技術(shù)， 2012，31（2）：12-14.

［15］許振成.流域水環(huán)境突發(fā)事件應(yīng)急處置技術(shù)體系研究 ［A］//中國環(huán)境科學學會學術(shù)年會論文集.北京：中國環(huán)境科學出版社，2012：1344-1353.

［16］呂小明.環(huán)境污染事件應(yīng)急處理技術(shù) ［M］. 北京：中國環(huán)境科學出版社， 2012.

［17］張曉健，陳 超，米子龍，等.飲用水應(yīng)急除鎘凈水技術(shù)與廣西龍江河突發(fā)環(huán)境事件應(yīng)急處置 ［J］. 城鎮(zhèn)給水排水， 2013，39（1）：24-32.

［18］張曉健，陳 超，林鵬飛.應(yīng)對水源突發(fā)污染的城市供水應(yīng)急處理技術(shù)研究與應(yīng)用 ［J］. 中國應(yīng)急管理， 2013，（10）：11-17.

［19］柳王榮，虢清偉，楊仁斌，等.混凝法去除鎘及含鎘絮體的穩(wěn)定性研究 ［J］. 水土保持學報， 2012，26（6）：80-84.

［20］卓瓊芳，許振成，虢清偉，等.聚鋁混凝沉淀物含鎘絮體穩(wěn)定性評估 ［J］. 中國環(huán)境科學， 2013，33（S1）：245-248.

A two-stage evaluation system to identify optimum emergency disposal technology schemes of sudden water source pollutions.

QU Jian-hua1， MENG Xian-lin1，2*， YOU Hong1，2（1.School of Municipal and Environmental Engineering，Harbin Institute of Technology， Harbin 150090， China；2.State Key Laboratory of Urban Water Resources and Environment， Harbin Institute of Technology， Harbin 150090， China）. China Environmental Science， 2015，35（10）：3193～3200

In order to predict the threat level of sudden pollution on water source for obtaining the optimum emergency disposal technology scheme， the threat-evaluation index system with 10evaluation indexes is established in this study. The weights of the evaluation indicators are determined in accordance to corresponding threat levels. Besides， Group-G1index weighting method based on circling correction is put forward for reflecting the influence of differences among experts on index weights in group decision making. With Spearman rank correlation coefficient as the test criterion， the consensus ranking of evaluation index and experts' weights are finally obtained through an iterative correction of onefold expert evaluation results with the combined evaluation methods integrating mean value method， Board method and Copeland method， so that the indexes are more accurate. In addition， the feasibility of this synthetic evaluation model has been verified through a sudden cadmium pollution incident that occurred in the Long River in 2012.

threat evaluation；drinking water source；circulation-correction；Group-G1；emergency disposal technology

X522

A

1000-6923（2015）10-3193-08

曲建華（1987-），男，黑龍江密山人，哈爾濱工業(yè)大學博士研究生，主要從事環(huán)境規(guī)劃與管理、環(huán)境工程技術(shù)評價方面的研究.發(fā)表論文5篇.

2015-03-16

國家環(huán)保公益性行業(yè)科研專項項目（201209048）

* 責任作者， 副教授， hgdmxl@163.com