任曉強,李 娜,管孝艷,3,陶 園,3,彭 舟
(1.中國水利水電科學(xué)研究院水利研究所,北京100048;2.中國灌溉排水發(fā)展中心,北京100054;3.國家節(jié)水灌溉北京工程技術(shù)研究中心,北京100048)
隨著我國人口不斷增長、工業(yè)化進程加快以及社會經(jīng)濟發(fā)展,誘使水環(huán)境問題凸顯,水環(huán)境風險分析越來越被關(guān)注,其可為水環(huán)境風險管理提供參考和依據(jù)[1]。作為“華北明珠”,白洋淀是雄安新區(qū)發(fā)展建設(shè)的重要依托,其在促進流域周邊城市生態(tài)建設(shè)、調(diào)節(jié)生態(tài)氣候、保護物種多樣性和維持華北地區(qū)的生態(tài)平衡等方面發(fā)揮著十分重要的作用[2]。上游入淀河流作為白洋淀的天然水源,其水環(huán)境質(zhì)量直接影響白洋淀的環(huán)境治理和生態(tài)恢復(fù),因此,白洋淀上游入淀河流水環(huán)境風險評價十分必要[3]。
水環(huán)境風險評價是水環(huán)境風險管理和污染防治的基礎(chǔ)。目前,國內(nèi)外研究學(xué)者在水環(huán)境風險評價的理論基礎(chǔ)和實際應(yīng)用方面都取得了一定的進展。20世紀80年代初,美國國家科學(xué)研究院(NAS)和環(huán)保局(USEPA)首先提出了基于危害鑒定、劑量-效應(yīng)評價、暴露評價、風險表征的“四步走”方法用以評價水中污染物給人體健康所帶來的風險[4];Nazar 等[5]較早提出了表征環(huán)境風險的量化指標,為環(huán)境風險從定性評價轉(zhuǎn)為定量評價提供了基礎(chǔ);祝慧娜等[6]綜合考慮水環(huán)境污染的脆弱性和對人體健康的危害性,構(gòu)建了河流水環(huán)境污染風險模糊綜合評價模型,并將其成功運用于湘江水環(huán)境風險評價中;魏娜等[7]將世界自然基金會(WWF)提出的水環(huán)境風險評估指標進行了本土化的調(diào)整,構(gòu)建了適應(yīng)中國國情的流域水環(huán)境風險評估指標體系。經(jīng)過了多年研究發(fā)展,國內(nèi)外已基本形成了比較成熟的區(qū)域水環(huán)境風險評價指標和評價體系,推動了水環(huán)境風險評價在實踐中的應(yīng)用。
然而,目前水環(huán)境風險評價多以流域、水源地為研究對象[8],對河流水環(huán)境風險的研究尚少,且大多注重水質(zhì)、水量或者單一污染物的風險研究[6],因此亟待開展綜合性河流水環(huán)境風險評價研究。本研究以白洋淀上游入淀河流為研究對象,考慮河流上下游河段風險傳遞性,從風險源、河流特性、自然地理及社會發(fā)展條件、水污染治理能力等因素出發(fā),構(gòu)建了白洋淀上游河流對白洋淀水環(huán)境影響的風險評價體系,并利用改進層次分析法(IAHP)和水污染綜合指數(shù)法對白洋淀上游入淀河流進行風險評價,為河流及白洋淀水環(huán)境污染防治和風險管理提供參考。
白洋淀流域位于華北平原中部,地處山西、河北、北京三省交界處,地理位置為東經(jīng)113°39′~116°11′,北緯39°4′~40°4′,屬于海河流域大清河水系,流域面積約為3.11 萬km2[2]。白洋淀位于流域下游地區(qū),是華北平原最大的淡水湖[2]。白洋淀上游共有8條入淀河流,分別為潴龍河、孝義河、唐河、府河、漕河、瀑河、萍河和白溝引河(潴龍河和萍河因常年斷流,不納入研究范圍)(圖1)。白洋淀流域人口稠密,水資源總量相對匱乏,社會正常經(jīng)濟活動的用水量持續(xù)增加,除府河、孝義河、白溝引河、瀑河長期有水外,其余入淀河流為季節(jié)性河流,部分依靠引黃入冀或南水北調(diào)維持水面[3]。白洋淀上游入淀河流水質(zhì)較差,絕大多數(shù)河流監(jiān)測斷面水質(zhì)處于Ⅳ~劣Ⅴ類水平,主要污染物為化學(xué)需氧量(COD)、總氮(TN)、氨氮(NH3-N)等[11]。
圖1 白洋淀上游入淀河流Fig.1 The rivers into Baiyangdian Lake in upper Baiyangdian basin
收集白洋淀流域數(shù)字高程模型(DEM)數(shù)據(jù),并利用ArcGIS中的水文分析模塊劃分白洋淀子流域,同時考慮流域內(nèi)的行政區(qū)劃、水功能區(qū)劃以及《水污染防治行動計劃》(即‘水十條’)的要求,對白洋淀子流域的邊界范圍進行修正[12]。子流域劃分完成后,將白洋淀上游入淀河流,按照縣區(qū)分成若干個河段單元,同時參考白洋淀各子流域的控制范圍,以及各個縣區(qū)的行政邊界,劃分每個河段單元的流域控制范圍,最終生成22 個水環(huán)境風險評價單元,文中以“河流名稱—縣級行政區(qū)名稱”的方式命名評價單元。
水環(huán)境風險評價指標體系的建立是進行水環(huán)境風險評價的前提和基礎(chǔ)。考慮到風險沿河流方向的傳遞性和累積性,把一條河流相鄰兩個河段中的上游河段看作下游河段的潛在風險來源之一,在遵循科學(xué)性、代表性、全面性、可比性、方便性等原則基礎(chǔ)上[13],從風險源、河流特性、自然地理及社會發(fā)展條件、水污染治理能力等因素出發(fā),構(gòu)建了考慮河流上游河段對下游河段影響的水環(huán)境風險評價體系(表1)。
表1 河流水環(huán)境風險評價指標體系Tab.1 River water environment risk assessment index system
(1)風險源子系統(tǒng)。風險源的有效辨識是水環(huán)境風險評價的前提和基礎(chǔ)。風險源子系統(tǒng)中選擇單位耕地面積化肥施用量、畜禽糞便排泄量、農(nóng)村生活污水排放量、農(nóng)村人體糞尿排泄量作為面源污染的評價指標。其中,單位耕地面積化肥施用量反映研究區(qū)域種植業(yè)污染程度,畜禽糞便排泄量反映研究區(qū)域畜禽養(yǎng)殖業(yè)污染程度,農(nóng)村生活污水排放量、農(nóng)村人體糞尿排泄量反映研究區(qū)域農(nóng)村生活源污染程度[14]。選取單位GDP 工業(yè)廢水排放量作為點源污染的評價指標。此外,在河流下游河段水環(huán)境風險評價中,還需考慮來自相鄰上游河段的水環(huán)境影響,因此,選擇相鄰上游河段的風險值及其到相鄰上游河段的距離作為下游河段污染的評價指標。
(2)河流特性子系統(tǒng)。河流特性可以直接反映河流水資源與水環(huán)境現(xiàn)狀,也是水環(huán)境評價系統(tǒng)的重要組成部分。以河流水質(zhì)現(xiàn)狀、河流水功能區(qū)達標率指標直接反映河流水環(huán)境現(xiàn)狀,多年平均年排水量或徑流深×流域面積直接反映河流水資源現(xiàn)狀以及其對污染物濃度及污染物遷移的影響。
(3)自然地理及社會發(fā)展條件子系統(tǒng)。自然地理及社會發(fā)展條件是水環(huán)境風險的間接驅(qū)動要素。在流域的發(fā)展過程中,人類的社會經(jīng)濟活動起主導(dǎo)作用。選取人口數(shù)量、人口自然增長率和人均GDP 作為社會發(fā)展的評價指標,其中人口數(shù)量、人口自然增長率反映研究區(qū)域人口規(guī)模和變化情況,人均GDP 反映研究區(qū)域經(jīng)濟發(fā)展水平的高低[15]。選取污染源到白洋淀的距離指標來反映研究區(qū)域內(nèi)污染物遷移到流域出口-白洋淀的衰減作用。
(4)水污染治理能力子系統(tǒng)。水污染治理能力是河流水環(huán)境污染防治和風險管理的主要驅(qū)動要素。采用工業(yè)廢水排放達標率反映工業(yè)點源污染控制和治理水平,采用城鎮(zhèn)生活污水集中處理率反映研究區(qū)域內(nèi)城鎮(zhèn)生活污水集中收集處理設(shè)施的配套程度,采用城鎮(zhèn)垃圾收集率反映研究區(qū)域內(nèi)對城鎮(zhèn)生活垃圾的處理能力[16],采用環(huán)保投資占GDP比重以反映研究區(qū)域內(nèi)對環(huán)境保護投入的重視程度。
層次分析法(AHP)是一種較為常見的賦權(quán)方法。然而,傳統(tǒng)層次分析法(AHP)采用Saaty[17]提出的1~9 標度法來構(gòu)造判斷矩陣。當評價體系中指標較多時,不僅迭代次數(shù)較多,運算量大,而且判斷矩陣很可能因一致性較差(甚至不具有一致性)而導(dǎo)致計算出的指標權(quán)重可靠性較低。因此,本研究運用已有的改進層次分析法(IAHP)[17,18],采用三標度法(0,1,2),對指標進行兩兩比較,建立比較矩陣,進而確定判斷矩陣。該方法具有自調(diào)節(jié)功能,不需要進行一致性檢驗,與傳統(tǒng)層次分析法(AHP)相比,標度值具有合理性和良好的判斷傳遞性,在比較判斷過程中提高了準確性。該方法計算步驟包括構(gòu)建風險評估指標體系、建立調(diào)查問卷并邀請水環(huán)境領(lǐng)域的專家對各評價指標的相對重要性給予判斷、構(gòu)造判斷矩陣、求解特征向量、層次單排序、層次總排序確定各評價指標的最終權(quán)重。最終的計算結(jié)果見表1。
研究區(qū)域內(nèi)各縣區(qū)的人口數(shù)量、經(jīng)濟生產(chǎn)總值(GDP)、人口自然增長率、耕地面積、化肥施用量(折純量)、畜禽養(yǎng)殖量、工業(yè)廢水排放達標率、城鎮(zhèn)生活污水集中處理率、城鎮(zhèn)垃圾收集率、農(nóng)村人口數(shù)量等數(shù)據(jù)主要來自于各地區(qū)2017年年鑒或經(jīng)濟統(tǒng)計年鑒;水質(zhì)數(shù)據(jù)主要來自于2017年政府網(wǎng)站發(fā)布的部分斷面水環(huán)境質(zhì)量以及實地采樣檢測數(shù)據(jù);節(jié)能環(huán)保支出數(shù)據(jù)主要來自于2017年各縣區(qū)政府公共預(yù)算支出;工業(yè)廢水數(shù)據(jù)主要來自于生態(tài)環(huán)境廳/局給出的污染重點監(jiān)測企業(yè)相關(guān)信息;畜禽糞便排泄系數(shù)、農(nóng)村人體糞尿產(chǎn)生系數(shù)主要來自于文獻[19]、農(nóng)村人均生活污水產(chǎn)生系數(shù)主要來自于文獻[20]。
畜禽糞便排泄量、農(nóng)村生活污水排放量、農(nóng)村人體糞尿排放量等面源污染指標采用輸出系數(shù)模型來計算[21,22]。研究單元到上游河段的距離、污染源到白洋淀的距離等指標利用河流水系圖以及ArcGIS 軟件來實現(xiàn)。上游河段的風險值是采用綜合指數(shù)加權(quán)求和法來實現(xiàn)[7]。
白洋淀上游入淀河流水環(huán)境風險評價指標等級劃分的依據(jù)是:①優(yōu)先采用國家標準或者省、市標準規(guī)定的等級分類值;②相關(guān)研究領(lǐng)域的參考文獻較多時,參照文獻中指標的分級標準,并結(jié)合研究區(qū)域的實際情況來對指標的評價等級進行劃分[23];③當沒有國家、省市標準或者相關(guān)的參考文獻分級數(shù)據(jù)可以參考的情況下,使用韋伯-費希納定律(W-F 定律)來計算指標的分級標準界限值[24]。白洋淀上游入淀河流水環(huán)境風險評價指標分級標準見表2。
表2 河流水環(huán)境風險評價指標分級標準Tab.2 Classification standard of river water environment risk assessment index
白洋淀上游入淀河流水環(huán)境風險綜合評分值的計算采用綜合指數(shù)加權(quán)求和法,即用各評估指標的評分值與它的權(quán)重系數(shù)相乘并求和來反映河流水環(huán)境風險的大小。其數(shù)學(xué)表達式為:
式中:R為水環(huán)境風險綜合評分值;Si為各評價指標的評分值;Wi為各評價指標的權(quán)重值;i為各評價指標;n為評價指標的個數(shù)。
河流水環(huán)境風險綜合評分值的計算結(jié)果越大,水環(huán)境風險程度越大,反之越小。若最終計算結(jié)果屬于(0,1],則風險等級為Ⅰ(無風險或者可接受風險),風險產(chǎn)生概率極低或破壞性極弱;屬于(1,2],風險等級為Ⅱ(低風險),要約束用水行為來防范風險;屬于(2,3],風險等級為Ⅲ(中風險),風險發(fā)生或者潛在存在造成一定損害;屬于(3,4],風險等級為Ⅳ(高風險),風險極易發(fā)生并造成極大破壞;屬于(4,5],風險等級為Ⅴ(極高風險),風險發(fā)生頻繁且造成不易恢復(fù)性破壞[7]。
2017年白洋淀上游入淀河流水環(huán)境風險評價計算結(jié)果如圖2 和表3、4所示。從計算結(jié)果(圖2 和表3)來看,2017年白洋淀上游入淀河流總體風險等級為Ⅲ級(個別河段除外),屬于中風險區(qū)域,即存在潛在風險并會給下游的白洋淀造成一定損害。孝義河的定州段、安國段、博野段及蠡縣段風險等級為Ⅳ級,屬于高風險區(qū)域,該區(qū)域極易發(fā)生水環(huán)境風險并會對下游的白洋淀造成極大影響;唐河—淶源段風險等級為Ⅱ級,是研究區(qū)域內(nèi)唯一的低風險區(qū)域。
表3 河流水環(huán)境風險評估結(jié)果Tab.3 Assessment results of water environment risk of the rivers
圖2 河流水環(huán)境風險分布Fig.2 Distribution of water environment risk of the rivers
對不同河段水環(huán)境風險評價指標得分熱力圖(表4)進行分析可知,影響白洋淀上游入淀河流水環(huán)境風險水平的主要原因在于環(huán)保投資占GDP 比重、單位耕地面積化肥施用量、河流水質(zhì)現(xiàn)狀和河流水功能區(qū)達標率等指標。白洋淀上游入淀河流周邊多數(shù)縣區(qū)對水環(huán)境保護的重視程度稍顯不足,環(huán)保投資占GDP 比重較低,按照國際經(jīng)驗,當?shù)貐^(qū)的環(huán)保投資占GDP 比重小于1%時,難以控制環(huán)境惡化的趨勢。白洋淀上游入淀河流主要分布在大清河淀西平原內(nèi)部,其農(nóng)業(yè)生產(chǎn)條件很好、農(nóng)業(yè)經(jīng)濟比較發(fā)達,是我國的主要產(chǎn)糧區(qū)[25]。絕大多數(shù)風險評估單元控制的流域范圍內(nèi)化肥施用強度達到國際公認化肥施用安全上限(225 kg/hm2)的3~5倍之多,給周邊河流帶來較大的環(huán)境壓力。此外,白洋淀上游入淀河流多數(shù)河段的水質(zhì)現(xiàn)狀較差,河流水功能區(qū)達標率較低,河流水環(huán)境長期面臨較大的風險壓力,特別是對孝義河定州段、安國段、博野段及蠡縣段等水環(huán)境高風險河段。
表4 不同河段水環(huán)境風險評價指標得分熱力圖Tab.4 The height map of water environment risk assessment index scores in different river sections
除上述評價指標影響外,孝義河—曲陽段、唐河—渾源段、唐河—靈丘段和唐河—安新段的主要風險成因還包括這些河段控制范圍內(nèi)人均GDP 指標水平比較低;漕河—易縣段、漕河—徐水段、唐河—靈丘段、唐河—唐縣段和唐河—定州段等河段控制范圍內(nèi)畜禽養(yǎng)殖業(yè)較發(fā)達,畜禽糞便排泄量指標處于較高水平,加之基礎(chǔ)設(shè)施的不完善,畜禽糞便無法得到妥善處理,給周邊的河流水環(huán)境帶來較大風險。
通過上述分析,在白洋淀上游入淀河流生態(tài)環(huán)境的治理與修復(fù)方面,未來應(yīng)進一步加大環(huán)境保護投資,妥善處理經(jīng)濟增長與農(nóng)業(yè)面源污染的關(guān)系,降低化肥施用強度,優(yōu)化畜禽養(yǎng)殖模式,實施畜禽糞便資源化利用,同時進一步加強對河流水體的治理力度,改善河流水質(zhì)條件。
構(gòu)建了考慮河流上下游河段風險傳遞性的白洋淀上游入淀河流水環(huán)境風險評價指標體系,并運用改進層次分析法(IAHP)和綜合指數(shù)加權(quán)求和法計算了各評估單元的水環(huán)境風險綜合評分值。結(jié)果顯示。
(1)2017年,白洋淀上游入淀河流總體處于中風險區(qū)域,孝義河—定州段、安國段、博野段和蠡縣段屬于高風險區(qū)域,極易發(fā)生水環(huán)境風險并會對白洋淀區(qū)域造成極大破壞;唐河—淶源段是研究區(qū)域內(nèi)僅有的低風險區(qū)域。
(2)影響白洋淀上游入淀河流水環(huán)境風險等級的主要原因是環(huán)保投資占GDP比重、單位耕地面積化肥施用量和河流水質(zhì)現(xiàn)狀等指標。
(3)在白洋淀上游入淀河流生態(tài)環(huán)境的治理與修復(fù)方面,應(yīng)加大環(huán)境保護投資,妥善處理經(jīng)濟增長與農(nóng)業(yè)面源污染的關(guān)系,降低化肥施用強度,優(yōu)化畜禽養(yǎng)殖模式,實施畜禽糞便資源化利用,同時加強對河流水體的治理力度,改善河流水質(zhì)條件。□