南水北調(diào)中線陜西水源區(qū)污染源排放時(shí)空變化特征

王 蕾，姚志鵬，吳 蕊，關(guān)建玲，羅儀寧，柴楊揚(yáng)，和 瑩

1.陜西省環(huán)境監(jiān)測中心站，陜西 西安 710054 2.中國環(huán)境監(jiān)測總站，國家環(huán)境保護(hù)環(huán)境監(jiān)測質(zhì)量控制重點(diǎn)實(shí)驗(yàn)室，北京 100012 3.中聯(lián)西北工程設(shè)計(jì)研究院有限公司，陜西 西安 710082

南水北調(diào)中線陜西水源區(qū)污染源排放時(shí)空變化特征

王 蕾1，姚志鵬2，吳 蕊3，關(guān)建玲1，羅儀寧1，柴楊揚(yáng)1，和 瑩1

1.陜西省環(huán)境監(jiān)測中心站，陜西 西安 710054 2.中國環(huán)境監(jiān)測總站，國家環(huán)境保護(hù)環(huán)境監(jiān)測質(zhì)量控制重點(diǎn)實(shí)驗(yàn)室，北京 100012 3.中聯(lián)西北工程設(shè)計(jì)研究院有限公司，陜西 西安 710082

基于2008—2012年污染源環(huán)境統(tǒng)計(jì)數(shù)據(jù)，采用Spearman秩相關(guān)系數(shù)法等分析南水北調(diào)中線陜西水源區(qū)污廢水及污染物年排放總量變化趨勢，借助ArcGIS空間分析功能表征水源區(qū)污染源排放的空間分布特征。結(jié)果表明：2008—2012年水源區(qū)污廢水和NH3-N年排放總量呈顯著上升趨勢，COD、As、Pb、Cd、Cr和Hg年排放總量呈拋物線型變化，總體上COD、Pb、Cd年排放總量增加，As、Cr和Hg年排放總量降低。水源區(qū)污廢水及污染物排放量空間差異明顯，污廢水、COD和NH3-N排放涉及流域所有區(qū)縣，其排放量從干流到流域邊緣呈現(xiàn)較明顯的梯度變化，即位于流域中心或地級(jí)市行政中心的區(qū)縣排放量明顯高于位于流域邊緣的區(qū)縣；As、Pb、Cd、Cr和Hg排放量呈現(xiàn)明顯的區(qū)域分布，主要分布在勉縣等8個(gè)區(qū)縣。

南水北調(diào)中線；陜西；水源區(qū)；污染源；時(shí)空變化

流域污染源排放是影響河流水質(zhì)的主導(dǎo)因素。自2007年國家污染源普查工作開展以來，作為防治流域水污染的有效手段，污染源調(diào)查和解析一直備受關(guān)注[1]。陜西省漢丹江流域位于陜西省南部的秦巴山地，承擔(dān)著南水北調(diào)中線工程供水水源丹江口水庫70%以上的供應(yīng)水量，是南水北調(diào)中線的主要水源涵養(yǎng)區(qū)[2]。近年來陜西省漢丹江流域水質(zhì)總體保持良好，漢丹江出陜斷面水質(zhì)保持在國家《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)Ⅱ～Ⅲ類[3]。然而，隨著陜西省漢丹江流域社會(huì)經(jīng)濟(jì)等的快速發(fā)展，污染物直接或間接進(jìn)入河道水體，使得南水北調(diào)中線陜西水源區(qū)水環(huán)境質(zhì)量面臨巨大的壓力[4]。隨著南水北調(diào)中線工程的正式通水，陜西水源區(qū)水環(huán)境狀況及污染源控制已經(jīng)受到社會(huì)各界廣泛關(guān)注。本研究依據(jù)南水北調(diào)中線陜西水源區(qū)環(huán)境統(tǒng)計(jì)數(shù)據(jù)，結(jié)合實(shí)地調(diào)研結(jié)果，分析水源區(qū)污廢水及各污染物排放總量變化趨勢及其空間分布特征，為區(qū)域污染源控制及水資源保障提供理論依據(jù)和數(shù)據(jù)支持。

1 研究方法

1.1 研究范圍和數(shù)據(jù)來源

本文污染源研究范圍涉及陜西省境內(nèi)漢丹江干支流主要流經(jīng)的陜南3個(gè)城市27個(gè)區(qū)縣，詳見圖1。污染源數(shù)據(jù)來源于2008—2012年陜西省環(huán)境統(tǒng)計(jì)數(shù)據(jù)庫，包括研究區(qū)域內(nèi)的工業(yè)源、農(nóng)業(yè)源、生活源和垃圾處理廠排污情況。

圖1 南水北調(diào)中線陜西水源區(qū)地理位置示意圖Fig.1 The geographical location in Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project

1.2 分析方法

由于陜西省漢丹江流域重金屬污染負(fù)荷較大[5]，本研究綜合考慮南水北調(diào)中線陜西水源區(qū)污染源類別特征，選取污廢水、總量控制污染物(COD、NH3-N)和重金屬污染物(As、Pb、Cd、Cr和Hg)指標(biāo)作為分析對(duì)象。研究從污染源排放總量統(tǒng)計(jì)等方面進(jìn)行歸納和分析，空間尺度涉及地級(jí)市、區(qū)縣兩級(jí)行政單位，同時(shí)運(yùn)用ArcGIS制圖以描述縣級(jí)行政單位排污空間分布特征。由于略陽縣漢江流域面積占全縣面積的29%，按其總量的1/3核算相關(guān)量。污廢水及各污染物排放年際變化趨勢分析采用Spearman秩相關(guān)系數(shù)法，用Daniel趨勢檢驗(yàn)技術(shù)進(jìn)行定量評(píng)價(jià)，其秩相關(guān)系數(shù)計(jì)算公式：

(1)

di=(Xi-Yi)

(2)

式中：Rs為秩相關(guān)系數(shù)，N為年份，di為變量Xi與Yi的差值；Xi為污廢水或某污染物按年排放總量從小到大排列對(duì)應(yīng)的序號(hào)，Yi為按時(shí)間排列的序號(hào)。將秩相關(guān)系數(shù)Rs的絕對(duì)值同臨界值Wp進(jìn)行比較：當(dāng)Rs>Wp，表明變化趨勢有顯著意義；當(dāng)Rs≤Wp，則表明變化趨勢沒有顯著意義；如果Rs是負(fù)值，則表明指標(biāo)在評(píng)價(jià)時(shí)段內(nèi)呈下降趨勢；如果Rs為正值，則表明指標(biāo)在評(píng)價(jià)時(shí)段內(nèi)呈上升趨勢。

2 結(jié)果與討論

2.1 南水北調(diào)中線陜西水源區(qū)污染源排放時(shí)間變化特征

南水北調(diào)中線陜西水源區(qū)2008—2012年污廢水及各污染物年排放總量變化情況如圖2所示，結(jié)合Spearman秩相關(guān)系數(shù)法判斷污廢水及各污染物年際變化趨勢分析，結(jié)果表明：①南水北調(diào)中線陜西水源區(qū)2008—2012年污廢水排放總量總體呈顯著上升趨勢(Rs=0.9)，其中2010年前增加幅度較大，2010年后基本保持同一水平；COD排放總量呈拋物線型變化，以2010年為拐點(diǎn)，之前呈上升趨勢，之后呈下降趨勢；NH3-N年排放量呈顯著上升趨勢(Rs=0.9)。與2008年相比，水源區(qū)2012年污廢水、COD和NH3-N排放總量分別增長了53.22%、29.44%和114.15%。②南水北調(diào)中線陜西水源區(qū)2008—2012年As、Pb、Cd、Cr、Hg年排放總量均為拋物線型變化，Cr和Hg在2009年排放量最大，As和Pb在2010年排放量最大，Cd在2011年排放量最大(圖2)。這是因?yàn)椋?009年陜西鳳翔血鉛事件以后，陜西省政府積極實(shí)施《重金屬污染防治規(guī)劃》相關(guān)措施，加大了區(qū)域重金屬污染防控力度，促使區(qū)域各重金屬排放量在2009年以后先后呈下降趨勢?？傮w上，南水北調(diào)中線陜西水源區(qū)Cd年排放總量呈顯著上升趨勢(Rs=0.9)，Pb年排放量變化趨勢不明顯(Rs=0.1)，Cr和Hg年排放量呈顯著下降趨勢(Rs=-0.9)，As年排放量呈不顯著下降趨勢(Rs=-0.6)；與2008年相比，2012年水源區(qū)Cd和Pb年排放量分別增加了165.25%和0.11%；As、Cr和Hg年排放量分別減少了15.74%、64.99%和19.84%。

圖2 南水北調(diào)中線陜西水源區(qū)污廢水及污染物年排放總量變化Fig.2 Variation of total discharge amounts of wastewater and pollutants in Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project

為了明晰南水北調(diào)中線陜西水源區(qū)污廢水及各污染物的來源情況，進(jìn)一步計(jì)算2008—2012年工業(yè)源、農(nóng)業(yè)源、生活源和垃圾處理源分別對(duì)污廢水及各污染物的排放量占比，結(jié)果顯示：南水北調(diào)中線陜西水源區(qū)污廢水主要來源為生活污水(占55.33%以上)和工業(yè)廢水(占20.77%以上)，COD和NH3-N主要來源為城鎮(zhèn)生活源(占46.89%以上)、農(nóng)業(yè)源(占37.57%以上)和工業(yè)源(占11.50%以上)，重金屬主要來源為工業(yè)源(占85.85%以上)；2008—2012年COD增加量主要來源于生活污水排放，NH3-N增加量主要來源于農(nóng)業(yè)源中畜禽養(yǎng)殖業(yè)排放，重金屬Cd和Pb增加量主要來源于工業(yè)源的鉛鋅冶煉、無機(jī)酸制造等涉重行業(yè)。

2.2 南水北調(diào)中線陜西水源區(qū)污染源排放空間分布特征

2.2.1 污廢水排放總量空間分布特征

南水北調(diào)中線陜西水源區(qū)2012年污廢水排放總量統(tǒng)計(jì)結(jié)果表明：①南水北調(diào)中線陜西水源區(qū)污廢水年排放量空間差異明顯，水源區(qū)所在的漢中市污廢水排放量最大，其次是安康市和商洛市，其污廢水年排放量分別占水源區(qū)排放總量的47.39%、28.24%和24.38%。調(diào)研顯示，漢中市在陜南3個(gè)城市中人口密度最大、人口最多，工業(yè)企業(yè)數(shù)量也最多，生活污水和工業(yè)廢水排放量均較大。②南水北調(diào)中線陜西水源區(qū)27個(gè)區(qū)縣污廢水排放總量分布從漢丹江干流到流域邊緣呈現(xiàn)明顯的梯度變化，即漢中市的勉縣、城固縣、洋縣、漢臺(tái)區(qū)，安康市的漢濱區(qū)和商洛市的商州區(qū)、山陽縣這些接近漢江干流和丹江干流的沿岸城市污廢水排放量處于較高水平，特別是漢江流域地級(jí)市行政中心的漢臺(tái)區(qū)、漢濱區(qū)在流域各區(qū)縣中污廢水排放量最大(圖3)，分別占到流域總排放量的15.97%和12.66%。導(dǎo)致這一現(xiàn)象的可能原因：各地區(qū)水資源分布不均及其城鎮(zhèn)化發(fā)展水平和工業(yè)行業(yè)分布等不同；同時(shí)，各地區(qū)的環(huán)境保護(hù)政策制定和落實(shí)情況不同、企業(yè)排污管理規(guī)范進(jìn)程不一，也會(huì)使污染排放存在差異。

圖3 2012年南水北調(diào)中線陜西水源區(qū)污廢水排放總量空間分布Fig.3 Spatial distribution of total discharge amounts of wastewater in 2012 in Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project

2.2.2 污染物排放總量空間分布特征

2.2.2.1 總量控制污染物排放總量空間分布

南水北調(diào)中線陜西水源區(qū)2012年總量控制類污染物排放總量統(tǒng)計(jì)結(jié)果表明：①南水北調(diào)中線陜西水源區(qū)COD和NH3-N排放量空間差異明顯， COD和NH3-N排放量由大到小依次為漢中市、安康市和商洛市，與污廢水排放量大小排序一致。其中漢中市COD和NH3-N排放量分別占流域排放總量的41.48%和45.49%，安康市COD和NH3-N排放量分別占流域排放總量的32.49%和29.67%，商洛市COD和NH3-N排放量分別占流域排放總量的26.03%和24.84%。②南水北調(diào)中線陜西水源區(qū)27個(gè)區(qū)縣中COD和NH3-N排放總量較大的區(qū)縣一致，COD和NH3-N排放總量空間分布總體情況和污廢水排放總量空間分布類似，其排放總量從干流到流域邊緣均呈現(xiàn)較明顯的梯度變化，其排放總量較高的區(qū)縣分別是漢中市的勉縣、南鄭縣、城固縣、洋縣、西鄉(xiāng)縣，安康市的漢濱區(qū)和旬陽縣，商洛市的丹鳳縣、商南縣和山陽縣(圖4)。漢濱區(qū)和西鄉(xiāng)縣COD和NH3-N排放量在流域區(qū)縣中最大，漢濱區(qū)COD和NH3-N排放量對(duì)流域排放總量貢獻(xiàn)率分別為10.45%和10.23%，西鄉(xiāng)縣COD和NH3-N排放量對(duì)流域排放總量貢獻(xiàn)率分別為8.77%和9.37%。③對(duì)比污廢水排放總量較大的區(qū)縣(圖3)，發(fā)現(xiàn)漢臺(tái)區(qū)和商州區(qū)對(duì)流域污廢水排放總量貢獻(xiàn)較大，而其對(duì)流域COD和NH3-N排放總量貢獻(xiàn)并不大，這說明漢臺(tái)區(qū)和商州區(qū)排放的污廢水中COD和NH3-N濃度較低；而旬陽縣和西鄉(xiāng)縣污廢水排放量貢獻(xiàn)相對(duì)較小，但其COD和NH3-N排放量貢獻(xiàn)卻較大，這說明旬陽縣和西鄉(xiāng)縣排放的廢水中COD和NH3-N濃度較高。調(diào)研顯示，各區(qū)縣污廢水排放量主要來源均為生活污水，但旬陽縣和西鄉(xiāng)縣工業(yè)廢水排放量占比相對(duì)較大，由于牲畜屠宰、化肥制造等工業(yè)企業(yè)較多，COD和NH3-N排放量較大。

2.2.2.2 重金屬污染物排放總量空間分布

南水北調(diào)中線陜西水源區(qū)2012年各重金屬污染物排放總量統(tǒng)計(jì)結(jié)果表明：陜西水源區(qū)各重金屬排放總量呈現(xiàn)明顯的區(qū)域分布。除Hg外，As、Pb、Cd和Cr排放量由大到小依次為漢中市、商洛市和安康市。Hg排放量主要分布在漢中市勉縣(74.73%)和安康市的旬陽縣(15.87%)；As排放量主要分布在漢中市勉縣(67.21%)、商洛市的鎮(zhèn)安縣(20.47%)和商州區(qū)(6.44%)；Pb排放量主要分布在漢中市勉縣(50.92%)、商洛市商州區(qū)(21.06%)和安康市旬陽縣(18.11%)；Cd排放量主要分布在漢中市勉縣(48.28%)和略陽縣(17.03%)、商洛市商州區(qū)(29.65%)；Cr排放量主要分布在商洛市商南縣(49.06%)、漢中市漢臺(tái)區(qū)(16.72%)、西鄉(xiāng)縣(14.15%)和勉縣(12.25%)(圖5)。南水北調(diào)中線陜西水源區(qū)重金屬排放量主要來源為工業(yè)源，因此，重金屬的區(qū)域分布特征與區(qū)域內(nèi)涉重企業(yè)的分布情況及工業(yè)企業(yè)排污情況直接相關(guān)。

圖4 2012年南水北調(diào)中線陜西水源區(qū)總量控制污染物排放總量空間分布Fig.4 Spatial distribution of total discharge amounts of total-amount-control pollutants in 2012 in Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project

2.3 南水北調(diào)中線陜西水源區(qū)污染源排放與漢江、丹江水質(zhì)響應(yīng)特征

2.3.1 污染源排放時(shí)間變化與水質(zhì)響應(yīng)分析

統(tǒng)計(jì)分析南水北調(diào)中線陜西水源區(qū)漢江、丹江2008—2012年水體COD、NH3-N和5種重金屬As、Pb、Cd、Cr、Hg的年均值變化情況，結(jié)合第2.1節(jié)水源區(qū)污染源排放時(shí)間變化情況發(fā)現(xiàn)：①漢江、丹江大部分?jǐn)嗝鍯OD和NH3-N年均濃度隨區(qū)域COD和NH3-N年排放量變化，且排放量變化幅度越大，水質(zhì)和排放量的相關(guān)性越高，排放量變化幅度小時(shí)水質(zhì)和排放量的相關(guān)性則較低，這是因?yàn)楹恿魉|(zhì)受區(qū)域污染源排放、污染物遷移轉(zhuǎn)化、河流水量及水體的自凈作用等綜合因素影響，當(dāng)污染物排放量變化幅度較小時(shí)，其他因素影響彰顯，污染物排放量變化影響減小。②水源區(qū)重金屬排放量較大的年度水體中相關(guān)重金屬的檢出率較高，但水體重金屬含量和重金屬排放量沒有明顯的相關(guān)性。

2.3.2 污染源排放空間特征與水質(zhì)響應(yīng)分析

課題組已有研究顯示，依據(jù)我國《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)評(píng)價(jià)，南水北調(diào)中線陜西水源區(qū)2012年漢江、丹江水體5種重金屬As、Pb、Cd、Cr和Hg含量均在Ⅰ類標(biāo)準(zhǔn)限值內(nèi)，其水質(zhì)受重金屬排放的影響很??；COD和NH3-N年均值為Ⅱ類或Ⅲ類，個(gè)別斷面NH3-N月測定值為Ⅵ類，存在偶爾超標(biāo)現(xiàn)象，其水質(zhì)受COD和NH3-N排放的影響相對(duì)較大[6]。根據(jù)監(jiān)測斷面對(duì)應(yīng)的區(qū)域看，陜西境內(nèi)漢江水質(zhì)較差河段所在的區(qū)域位于漢中市的南鄭縣和城固縣，丹江水質(zhì)較差的河段所在的區(qū)域位于商洛市商州區(qū)。對(duì)比第2.2.1節(jié)和第2.2.2節(jié)水源區(qū)污廢水及污染物排放總量空間分布情況，結(jié)果表明：①水源區(qū)水質(zhì)較差的斷面所在區(qū)域與污廢水、總量控制類污染物COD和NH3-N排放量較大的區(qū)域基本一致，這進(jìn)一步說明水源區(qū)污染源排放對(duì)水質(zhì)的影響明顯。②水源區(qū)漢江、丹江不同斷面重金屬含量水平普遍較低，其水質(zhì)并沒有表現(xiàn)出受重金屬排放量的區(qū)域性分布影響。研究顯示，多數(shù)重金屬進(jìn)入河流水體后只有一小部分以自由金屬離子形式溶解在水中，絕大部分會(huì)迅速被懸浮顆粒表面所吸附，并在水動(dòng)力作用的搬運(yùn)過程中逐步沉降下來在河流沉積物中累積[7-8]。正是由于重金屬在河流水體存在的這種“瞬時(shí)性”，使得大多數(shù)河流水體中重金屬的含量很低，即使接近污染物排放口，水體中重金屬的含量也不高[9]。由此可見，目前僅依據(jù)我國監(jiān)測部門實(shí)行的河流水體月監(jiān)測數(shù)據(jù)，在評(píng)價(jià)河流重金屬污染狀況時(shí)具有局限性。對(duì)于河流重金屬長期污染，監(jiān)測評(píng)價(jià)河流沉積物中的重金屬污染水平更具代表性[10]。因此，建議評(píng)價(jià)河流重金屬污染狀況時(shí)同時(shí)監(jiān)測水體重金屬濃度和沉積物重金屬含量，并對(duì)其污染水平進(jìn)行綜合評(píng)價(jià)。

圖5 2012年南水北調(diào)中線陜西水源區(qū)重金屬污染物排放總量空間分布Fig.5 Spatial distribution of total discharge amounts of heavy metal pollutants in 2012 in Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project

3 結(jié)論

1)南水北調(diào)中線陜西水源區(qū)2008—2012年污廢水和NH3-N年排放總量呈顯著上升趨勢，COD、As、Pb、Cd、Cr和Hg年排放總量呈拋物線型變化，總體上COD、Pb、Cd年排放總量增加，As、Cr和Hg年排放總量降低；COD增加量主要來源于生活污水排放，NH3-N增加量主要來源于農(nóng)業(yè)源中畜禽養(yǎng)殖業(yè)排放，重金屬Cd和Pb增加量主要來源于工業(yè)源的鉛鋅冶煉、無機(jī)酸制造等涉重行業(yè)。

2)南水北調(diào)中線陜西水源區(qū)2012年污廢水及污染物排放總量空間分布差異明顯：污廢水、COD和NH3-N排放涉及流域所有區(qū)縣，其排放量從干流到流域邊緣均呈現(xiàn)較明顯的梯度變化，即位于流域中心或地級(jí)市行政中心的區(qū)縣排放量明顯高于其他流域邊緣區(qū)縣；As、Pb、Cd、Cr和Hg呈現(xiàn)明顯的區(qū)域分布，主要分布在漢中市漢臺(tái)區(qū)、略陽縣、西鄉(xiāng)縣和勉縣，安康市旬陽縣，商洛市商州區(qū)、鎮(zhèn)安縣和商南縣。

3)水源區(qū)污染源排放與漢江、丹江水質(zhì)響應(yīng)分析結(jié)果顯示：南水北調(diào)中線陜西水源區(qū)河流水質(zhì)較差的斷面與污廢水、總量控制類污染物排放量較大的時(shí)段或區(qū)域基本一致，表現(xiàn)出污染源排放對(duì)水源區(qū)水質(zhì)有明顯影響；水源區(qū)河流水體重金屬含量水平受重金屬排放量的影響很小，這是由于多數(shù)重金屬進(jìn)入河流水體后易于在河流沉積物中累積，使得僅監(jiān)測河流水體重金屬在綜合評(píng)價(jià)河流重金屬污染狀況時(shí)具有局限性。

[1] 曾祉祥,張洪,單保慶，等.漢江中下游流域工業(yè)污染源解析[J].長江流域資源與環(huán)境, 2014, 23(2): 252-259.

ZENG Zhexiang, ZHANG Hong, SHAN Baoqing, et al. Analysis of industrial pollution sources of the middle and lower reaches of Hanjiang River basin [J]. Resources and Environment in the Yangtze Basin, 2014, 23(2): 252-259.

[2] 賈忠華,趙恩輝.南水北調(diào)中線陜西水源區(qū)水沙變化特征分析[J].西北大學(xué)學(xué)報(bào):自然科學(xué)版,2010, 40(1):155-157, 161.

JIA Zhonghua, ZHAO Enhui. Analysis on changing characteristics of runoff and sediment in the water source area of the Middle Route of South-to-North Water Diversion Project[J]. Journal of Northwest University: Natural Science Edition, 2010, 40(1):155-157, 161.

[3] 張春玲, 周曉強(qiáng). 陜西省漢、丹江流域水資源質(zhì)量近年變化分析與保護(hù)對(duì)策研究[J]. 陜西水利, 2011 (5): 21-26.

ZHANG Chunling，ZHOU Xiaoqiang. Variation analysis and countermeasures studies on the water resources quality in recent years of Han river and Dan river in Shaanxi region [J].Shaanxi Water Resources,2011 (5): 21-26.

[4] 王蕾,關(guān)建玲,丁強(qiáng),等. 南水北調(diào)中線陜西水源區(qū)污染源變化特征[J].干旱區(qū)資源與環(huán)境,2015,29 (10): 80-85.

WANG Lei, GUAN Jianling, DING Qiang, et al. The variation characteristics of water pollution source of Shaanxi water source area of the Middle Route of South-to-North Water Diversion Project[J]. Journal of Arid Land Resources and Environment, 2015, 29 (10): 80-85.

[5] 汪新生,郭琦.陜西省重金屬污染特征分析[J].中國環(huán)境監(jiān)測,2011,27(4):22-27.

WANG Xinsheng, GUO Qi. Shaanxi heavy metal pollution characters analysis [J]. Environmental Monitoring in China, 2011, 27(4):22-27.

[6] 王蕾, 關(guān)建玲,姚志鵬,等. 漢丹江(陜西段)水質(zhì)變化特征分析[J].中國環(huán)境監(jiān)測,2015,31(5):73-77.

WANG Lei, GUAN Jianling, YAO Zhipeng, et al. The variation characteristics of water quality in the Han River and Dan River of Shaanxi Region[J]. Environmental Monitoring in China, 2015, 31(5):73-77.

[7] REZA R, SINGH G. Assessment of heavy metal contamination and its indexing approach for river water [J]. International Journal of Environment Science and Technology, 2010, 7(4):785-792.

[8] LI S Y, ZHANG Q F. Risk assessment and seasonal variations of dissolved trace elements and heavy metals in the Upper Han River, China [J]. Journal of Hazardous Materials, 2010, 18(1/3):1 051-1 058.

[10] DAVUTLUOGLU O I, SECKIN G, ERSU C B, et al. Heavy metal content and distribution in surface sediments of the Seyhan River, Turkey [J]. Journal of Environmental Management, 2011, 92(9): 2 250-2 259.

The Spatiotemporal Variation Characteristics of Water Pollution Sources in Shaanxi Water Source Area of the Middle Route of South-to-North Water Diversion Project

WANG Lei1, YAO Zhipeng2, WU Rui3, GUAN Jianling1, LUO Yining1, CHAI Yangyang1, HE Ying1

1.Shaanxi Environmental Monitoring Center, Xi’an 710054, China 2.State Environmental Protection Key Laboratory of Quality Control in Environmental Monitoring, China National Environmental Monitoring Centre, Beijing 100012, China 3.China United Northwest Institute for Engineering Design & Research Co., Ltd, Xi’an 710082, China

Based on the environment statistic data from 2008 to 2012, the water pollution sources variation characteristics were analyzed using the Spearman’s rank correlation coefficient method and the emission spatial distribution characteristics of water pollution sources were described using the spatial analysis function of ArcGIS in Shaanxi water source area of the middle route of South-to-North Water Diversion Project. The results indicated that the annual total discharge amounts of waste water and NH3-N in the region generally increased during 2008-2012, the discharge amounts of COD, As, Pb, Cd, Cr and Hg all showed as a parabolic variation, COD, Pb and Cd generally increased while As, Cr and Hg decreased. Discharge amounts of waste water and pollutants in the region showed obvious differences in spatial distribution, discharge of waste water, COD and NH3-N involved in all the counties, presented a significantly gradient variation, their discharge amounts of all the counties located in the center of river basin or in the administration center of cities were obviously higher than the counties in their marginal inferior. Discharge amounts of As, Pb, Cd, Cr and Hg displayed an apparently regional distribution, and most emission of heavy metals concentrated in eight countries including Mian County.

the Middle Route of South-to-North Water Diversion Project; Shaanxi; water source area; pollution sources; spatiotemporal variation

2015-11-18;

2016-06-24

國家科技支撐計(jì)劃項(xiàng)目(2011BAC12B01-01,2011BAC12B02-01)

王 蕾(1983-)，女，陜西商洛人，博士，高級(jí)工程師。

關(guān)建玲

X824

A

1002-6002(2017)01- 0061- 07

10.19316/j.issn.1002-6002.2017.01.10