南四湖水質(zhì)空間分布特征分析與改善效果評(píng)估

武周虎，張 可，金玲仁，楊連寬，張 建

(1.青島理工大學(xué)環(huán)境與市政工程學(xué)院，山東青島 266033;2.山東省環(huán)境保護(hù)廳流域環(huán)境管理處，山東濟(jì)南 250012;3.山東省環(huán)境監(jiān)測(cè)中心站，山東濟(jì)南 250013;4.濟(jì)寧市環(huán)境保護(hù)監(jiān)測(cè)站，山東濟(jì)寧 272045;5.山東大學(xué)環(huán)境科學(xué)與工程學(xué)院，山東 濟(jì)南 250100)

南水北調(diào)東線工程于2002年開工建設(shè)，計(jì)劃2013年正式通水。南四湖是南水北調(diào)東線工程重要的輸水通道和調(diào)蓄湖泊，要求到2012年底水質(zhì)穩(wěn)定達(dá)到GB3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中的Ⅲ類標(biāo)準(zhǔn)。近年來南四湖的水質(zhì)改善狀況備受關(guān)注[1-5]。為考量南四湖實(shí)施“治、用、?！绷饔蛑挝垠w系的效果，以及進(jìn)行南四湖水質(zhì)改善狀況評(píng)估和重大水環(huán)境問題診斷研究，山東省環(huán)境保護(hù)廳分別于2006年11月、2007年4月、2010年11月和2011年5月組織開展了4次南四湖水質(zhì)空間分布監(jiān)測(cè)工作[6]。從時(shí)間上看，2006年和2010年為汛后監(jiān)測(cè)，2007年和2011年為汛前監(jiān)測(cè)，因此，前2次監(jiān)測(cè)的平均值可代表2006—2007年枯水期(對(duì)應(yīng)東線工程實(shí)施后的調(diào)水期，下同)的水質(zhì)狀況，后2次監(jiān)測(cè)的平均值可代表2010—2011年枯水期的水質(zhì)狀況。

國內(nèi)外對(duì)湖泊水環(huán)境質(zhì)量現(xiàn)狀調(diào)查與評(píng)價(jià)的研究成果較多:周豐等[7]基于多元統(tǒng)計(jì)方法進(jìn)行了河流水質(zhì)空間分析，識(shí)別采樣點(diǎn)的空間相似性與差異性，為水質(zhì)監(jiān)測(cè)網(wǎng)絡(luò)優(yōu)化提供支持;劉碩等[8]對(duì)幾種環(huán)境質(zhì)量綜合指數(shù)評(píng)價(jià)方法進(jìn)行探討，提出宜采用最大值法和算術(shù)平均法相結(jié)合的方法來對(duì)全國環(huán)境質(zhì)量進(jìn)行評(píng)價(jià);范瑜[9]給出了平均綜合污染指數(shù)在地面水污染分級(jí)中的應(yīng)用研究，對(duì)地面水污染程度進(jìn)行分級(jí);鄧建才等[10]分析了南四湖水體氮、磷、Chl-a含量及空間分布特征，給出南四湖水體中葉綠素濃度的分布呈現(xiàn)出較大的差異;Kemka等[11]給出了城區(qū)湖泊富營養(yǎng)化研究成果，得出人類活動(dòng)是城區(qū)湖泊富營養(yǎng)化的主因;武周虎等[12]根據(jù)前2次南四湖水質(zhì)空間分布監(jiān)測(cè)數(shù)據(jù)，采用平均綜合污染指數(shù)法篩選出29個(gè)重污染測(cè)點(diǎn)，分析了南四湖污染物的主要來源。隨著南四湖流域扎實(shí)推進(jìn)污染綜合整治和通水時(shí)間的臨近，南四湖水質(zhì)空間分布狀況和水環(huán)境問題發(fā)生了較大變化，針對(duì)南四湖地形情況復(fù)雜的特點(diǎn)和受關(guān)注的水環(huán)境瓶頸問題，選用Surfer 8.0軟件和平均綜合污染指數(shù)法進(jìn)行南四湖水質(zhì)空間分布特征分析與污染控制效果評(píng)估，對(duì)南水北調(diào)東線南四湖輸水水質(zhì)安全保障綜合控制方案的修訂和治污工程實(shí)施具有重要的現(xiàn)實(shí)意義。

筆者根據(jù)2006—2007年(2次平均，下同)和2010—2011年兩個(gè)枯水期的水質(zhì)狀況，采用Surfer 8.0軟件繪制南四湖水質(zhì)空間分布等值線圖，進(jìn)行南四湖及分湖區(qū)水質(zhì)指標(biāo)的綜合比較分析和平均綜合污染指數(shù)的計(jì)算，開展南四湖水質(zhì)空間分布特征分析、污染控制效果評(píng)估與水環(huán)境問題的變化研究。

1 南四湖概況、采樣點(diǎn)布設(shè)及監(jiān)測(cè)指標(biāo)

南四湖是南陽湖、獨(dú)山湖、昭陽湖和微山湖4個(gè)由西北向東南相連而成的湖泊的總稱，流域面積31700 km2，東、西、北三面承接蘇、魯、豫、皖4省32個(gè)縣(市、區(qū))53條河流的來水。南四湖湖內(nèi)總面積1266 km2，南北狹長126 km，東西寬5～25 km，岸線總長311 km。為提高蓄洪、抗旱、供水、航運(yùn)和養(yǎng)殖等綜合效益，1960年在南四湖湖腰修建了7.36 km長的二級(jí)閘壩水利樞紐工程，將其分成上級(jí)湖和下級(jí)湖。二級(jí)閘壩以北為上級(jí)湖，湖內(nèi)面積606 km2，二級(jí)閘壩以南為下級(jí)湖，湖內(nèi)面積660 km2。南四湖正常蓄水位相應(yīng)的總庫容為20.36億m3，其中上級(jí)湖的庫容為9.24億m3(水位34.2 m)，下級(jí)湖的庫容為11.12億m3(水位33.0 m)。南四湖上級(jí)湖承納27500 km2的來水，占總來水面積的86.8% ，而庫容僅占南四湖總庫容的45.4%。根據(jù)南四湖湖區(qū)萬分之一地形圖(2001年CAD版)分析，南四湖湖盆淺平，北高南低，平均水深約2.0 m，天然情況是上級(jí)湖泄水入下級(jí)湖再入中運(yùn)河，而南水北調(diào)輸水是由下級(jí)湖南端進(jìn)入、經(jīng)二級(jí)壩泵站、由上級(jí)湖北端入梁濟(jì)運(yùn)河。南四湖湖內(nèi)開闊水域面積為781.45 km2，僅占湖內(nèi)總面積的61.7%，其余包括池塘、網(wǎng)圍和網(wǎng)箱等水產(chǎn)養(yǎng)殖區(qū)，總面積達(dá)到27 525 km2(約90%為網(wǎng)圍和池塘)、農(nóng)田15400 km2、陸生植被5530 km2。湖內(nèi)臺(tái)田中化肥和農(nóng)藥施用、水產(chǎn)養(yǎng)殖、航運(yùn)和分散居民生活污水排放等造成的面源污染，對(duì)湖泊水質(zhì)造成直接影響。

南四湖水質(zhì)空間分布監(jiān)測(cè)采用網(wǎng)格布點(diǎn)法在上級(jí)湖和下級(jí)湖均勻布設(shè)監(jiān)測(cè)點(diǎn)位，每個(gè)網(wǎng)格南北長約為3.32 km，東西寬約為 3.03 km，面積約為 10 km2，同時(shí)在主航道和湖區(qū)邊界增設(shè)采樣點(diǎn)。南四湖水質(zhì)空間分布監(jiān)測(cè)采樣點(diǎn)位及分區(qū)參見文獻(xiàn)[12]。

根據(jù)南四湖的歷史監(jiān)測(cè)資料分析，該水體的污染特征主要是結(jié)構(gòu)性污染，主要超標(biāo)因子為COD、TP、TN等[5]，且南四湖已處于不同程度的富營養(yǎng)狀態(tài)。因此，在南四湖水質(zhì)空間分布監(jiān)測(cè)中，確定主要監(jiān)測(cè)項(xiàng)目為水溫、水深、pH、COD、NH3-N、TP、TN、CODMn(2010年和2011年增加)、透明度、電導(dǎo)率和Chl-a。以2011年為例，共完成90個(gè)湖區(qū)點(diǎn)位和35個(gè)入湖河口及3個(gè)出湖河口點(diǎn)位的監(jiān)測(cè)，共報(bào)出有效數(shù)據(jù)2017個(gè)，其中現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)705個(gè)，實(shí)驗(yàn)室監(jiān)測(cè)數(shù)據(jù)1312個(gè)。利用河口數(shù)據(jù)另文分析污染物外源輸入的主要通道。

2 水質(zhì)空間分布特征分析

2.1 水質(zhì)指標(biāo)分布特征

根據(jù)2006—2007年和2010—2011年兩個(gè)枯水期南四湖90個(gè)點(diǎn)的水質(zhì)監(jiān)測(cè)數(shù)據(jù)，采用Surfer 8.0軟件分別繪制COD、NH3-N、TP、TN和CODMn的空間分布等值線于圖1～5。圖1～5中經(jīng)緯度坐標(biāo)采用十進(jìn)制，各水質(zhì)監(jiān)測(cè)指標(biāo)按GB3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中Ⅲ類(含Ⅰ類和Ⅱ類)、Ⅳ類、Ⅴ類和超Ⅴ類水對(duì)應(yīng)的數(shù)值區(qū)間依次填充為不同灰度，各水質(zhì)指標(biāo)的單位均為mg/L。

圖1 枯水期COD空間分布

由圖1可見，2006—2007年枯水期南四湖COD分布不存在Ⅲ類(達(dá)標(biāo))水，超Ⅴ類水主要分布在小北湖、南陽湖、獨(dú)山湖西側(cè)以及微山湖島北2個(gè)測(cè)點(diǎn)，Ⅳ類水測(cè)點(diǎn)零星分布，其余湖區(qū)幾乎全為Ⅴ類水;2010—2011年枯水期南四湖COD分布在微山湖島南、昭陽湖二級(jí)壩上下游局部出現(xiàn)了Ⅲ類水，超Ⅴ類水僅在南陽湖北端的老運(yùn)河和洸府河入湖河口外混合區(qū)出現(xiàn)了1個(gè)測(cè)點(diǎn)，在該點(diǎn)周圍出現(xiàn)了小范圍的Ⅴ類水，其余湖區(qū)全為Ⅳ類水。由此可見，4年來南四湖COD得到大幅度的改善，COD整體提高了1～2個(gè)水質(zhì)類別。

圖2 枯水期NH3-N空間分布

由圖2可見，2006—2007年枯水期南四湖NH3-N分布只有在南陽湖北端的老運(yùn)河和洸府河入湖河口外混合區(qū)以及微山湖島北新薛河和薛城小沙河入湖河口外混合區(qū)為超Ⅴ類水，在其周圍出現(xiàn)了小范圍的Ⅴ類和Ⅳ類水，其余湖區(qū)全達(dá)到Ⅲ類水的標(biāo)準(zhǔn);到2010—2011年枯水期南四湖NH3-N分布全部達(dá)到Ⅲ類水的標(biāo)準(zhǔn)。

圖3 枯水期TP空間分布

由圖3可見，2006—2007年枯水期南四湖TP在微山湖島南大部分湖區(qū)、微山湖島北西側(cè)和昭陽湖零星分布有Ⅲ類水，南陽湖北段為超Ⅴ類水，南陽湖南段和獨(dú)山湖西側(cè)為Ⅴ類水，其余中段湖區(qū)全為Ⅳ類水，南四湖TP分布北高南低的非均一性相當(dāng)明顯;2010—2011年枯水期南四湖TP分布Ⅲ類水湖區(qū)大大縮減，只有在微山湖島南端有少許分布，超Ⅴ類水湖區(qū)也大為縮減，只分布在南陽湖北端和微山湖島北新薛河入湖河口外混合區(qū)，在其周圍和泗河口外湖段出現(xiàn)了小范圍的Ⅴ類水，其余湖區(qū)全達(dá)到Ⅳ類水。南四湖TP分布變化的特點(diǎn)是達(dá)標(biāo)湖區(qū)和重污染湖區(qū)都縮減了，南四湖TP分布北高南低趨于均一化，最大濃度降低明顯，平均濃度也有所降低。這一現(xiàn)象與2010—2011年枯水期在南四湖南端出湖河流的下泄水量增加有關(guān)。

由圖4可見，2006—2007年枯水期南四湖TN在南陽湖和微山湖島北、島南的東側(cè)主航道均出現(xiàn)超Ⅴ類水的重污染范圍，在其周圍和獨(dú)山湖西側(cè)出現(xiàn)了小范圍的Ⅴ類水，只有在微山湖西南角出現(xiàn)了小范圍的Ⅲ類水;2010—2011年枯水期南四湖TN水質(zhì)大為改善，除南陽湖北端的老運(yùn)河和洸府河入湖河口外混合區(qū)出現(xiàn)超Ⅴ類和Ⅴ類水外，南陽湖大部、昭陽湖中段二級(jí)壩上下游和薛城小沙河入湖河口外混合區(qū)為Ⅳ類水，其余超過南四湖60%的湖區(qū)達(dá)到Ⅲ類水。由此可見，4年來南四湖TN得到大幅度的改善，TN整體提高了1～2個(gè)水質(zhì)類別。

圖4 枯水期TN空間分布

圖5 2010—2011年枯水期CODMn空間分布

2006—2007年枯水期CODMn未監(jiān)測(cè)，由圖5可見，2010—2011年枯水期南四湖CODMn在南陽湖、獨(dú)山湖和昭陽湖西北局部為Ⅳ類水，昭陽湖和微山湖大部達(dá)到Ⅲ類水的標(biāo)準(zhǔn)，微山湖島北零星存在Ⅳ類水，不存在Ⅴ類和超Ⅴ類水的湖區(qū)。南四湖CODMn整體質(zhì)量較好，但北部水質(zhì)明顯差于南部水體水質(zhì)。

2.2 水質(zhì)達(dá)標(biāo)與重污染范圍分析

表1給出了2006—2007年和2010—2011年兩個(gè)枯水期南四湖水質(zhì)及水深監(jiān)測(cè)結(jié)果的統(tǒng)計(jì)與比較。由表1可見，2010—2011年比2006—2007年枯水期南四湖 COD、NH3-N、TP、TN的下降率在22.9% ～51.1%之間，下降率最大的是TN，最小的是TP，南四湖各水質(zhì)指標(biāo)的平均濃度下降率為38.3%，說明4年來南四湖的水質(zhì)得到大幅度的改善。南四湖水體的透明度下降、電導(dǎo)率上升，2010—2011年枯水期山東發(fā)生了較為罕見的跨越秋季、冬季和春季3個(gè)季節(jié)的持續(xù)干旱，天然徑流水量相對(duì)較小，湖內(nèi)中水比例增加，南四湖流域的24個(gè)污水資源化與截蓄導(dǎo)用工程發(fā)揮了重要作用，才未造成南四湖干涸，且水質(zhì)仍處于好轉(zhuǎn)狀況。

南四湖水質(zhì)空間分布 COD、NH3-N、TP、TN和CODMn達(dá)到Ⅲ類范圍和超Ⅴ類(重污染，下同)范圍所占比例變化比較及位置說明分別列于表2和表3。

由表2可見，2010—2011年比2006—2007年枯水期南四湖COD、NH3-N和TN達(dá)到Ⅲ類的比例得到明顯增加，但TP的達(dá)標(biāo)區(qū)縮減了一半，原因可能與2010—2011年枯水期干旱少雨、天然徑流水量相對(duì)較小、湖內(nèi)中水比例增加以及4年來南四湖底泥中TP含量普遍上升有關(guān)。各水質(zhì)指標(biāo)達(dá)到Ⅲ類范圍的比例平均增大了35.9%，各水質(zhì)指標(biāo)達(dá)到Ⅲ類范圍的比例平均達(dá)到55.1%，即超過一半的湖區(qū)達(dá)到Ⅲ類水標(biāo)準(zhǔn)。監(jiān)測(cè)結(jié)果表明，南四湖底泥中TP平均質(zhì)量比由0.68g/kg增加到1.25 g/kg，平均升高了83.5%，而洸府河口底泥采樣點(diǎn)上升幅度最大，其底泥中TP質(zhì)量比從1.29 g/kg增加到3.71 g/kg，最大升高了187.6%。因此，南四湖入湖河流排入TP和底泥磷富集與釋放是南四湖TP的主要來源。

由表3可見，2010—2011年比2006—2007年枯水期南四湖COD、NH3-N、TP和TN超Ⅴ類的比例得到明顯下降，各水質(zhì)指標(biāo)超Ⅴ類范圍的比例平均縮減了92.0%，各水質(zhì)指標(biāo)超Ⅴ類范圍的比例平均下降到僅占1.3%，說明南四湖重污染湖區(qū)和主要污染物得到比較有效的控制。

3 污染控制效果評(píng)估

3.1 總體污染控制效果評(píng)估

為更全面的評(píng)估4年來南四湖水質(zhì)的改善效果，選用平均綜合污染指數(shù)法[12]進(jìn)行分析。

表1 南四湖水質(zhì)及水深監(jiān)測(cè)結(jié)果統(tǒng)計(jì)與比較

表2 南四湖各水質(zhì)指標(biāo)達(dá)到Ⅲ類比例比較及位置說明

表3 南四湖各水質(zhì)指標(biāo)超Ⅴ類比例比較及位置說明

i指標(biāo)j個(gè)測(cè)點(diǎn)的單項(xiàng)污染指數(shù)

第j個(gè)測(cè)點(diǎn)的平均綜合污染指數(shù)

全湖i指標(biāo)的單項(xiàng)污染分擔(dān)率

式中:i=1，2，…，5，依次表示主要水質(zhì)指標(biāo) COD、NH3-N、TP、TN、CODMn;j為南四湖采樣點(diǎn)編號(hào)(j=1，2，3，…，m);ρij為第 j個(gè)測(cè)點(diǎn) i水質(zhì)指標(biāo)的實(shí)測(cè)值;ρsi為i水質(zhì)指標(biāo)的評(píng)價(jià)標(biāo)準(zhǔn)限值(按Ⅲ類標(biāo)準(zhǔn)依次為20mg/L、1.0mg/L、0.05mg/L、1mg/L、6mg/L，下同)，n為參與評(píng)價(jià)的水質(zhì)指標(biāo)數(shù)。

根據(jù)平均綜合污染指數(shù)的大小可將水體分為達(dá)標(biāo)、基本達(dá)標(biāo)、污染和重污染4類，其對(duì)應(yīng)關(guān)系參見表4，表4還給出了2006—2007年和2010—2011年枯水期南四湖按平均綜合污染指數(shù)的評(píng)價(jià)結(jié)果。圖6給出了兩個(gè)枯水期4項(xiàng)主要污染物分擔(dān)率的比較結(jié)果。

表4 南四湖按平均綜合污染指數(shù)的評(píng)價(jià)結(jié)果

由表4可見，南四湖2010—2011年比2006—2007年枯水期污染和重污染湖區(qū)的比例大幅度減小，特別是重污染湖區(qū)下降到2.2%，達(dá)標(biāo)和基本達(dá)標(biāo)湖區(qū)的比例大幅度提高，達(dá)標(biāo)和基本達(dá)標(biāo)湖區(qū)的比例達(dá)到63.3%。按平均綜合污染指數(shù)逆排序篩選的重污染湖區(qū)(2個(gè)測(cè)點(diǎn))均位于南陽湖北端的老運(yùn)河和洸府河入湖河口外混合區(qū)。2010—2011年枯水期南四湖水質(zhì)平均綜合污染指數(shù)1.03，比2006—2007年枯水期該指數(shù)1.66下降了38.1%。

由圖6可見，南四湖2010—2011年比2006—2007年枯水期TN的分擔(dān)率明顯下降，COD和NH3-N的分擔(dān)率略有下降，而TP的分擔(dān)率明顯上升，TP占到4項(xiàng)主要污染物的41.2%，說明南四湖主要污染物控制不平衡，TN的控制效果相對(duì)較好，而TP的控制效果最差。2010—2011年枯水期南四湖主要污染物按分擔(dān)率排序由大到小依次為TP、COD、TN、CODMn、NH3-N，南四湖的 ρ(TN)/ρ(TP)由2006—2007年枯水期的18.5下降到2010—2011年枯水期的 11.7。

圖6 兩個(gè)枯水期4項(xiàng)主要污染物分擔(dān)率比較

3.2 分湖區(qū)污染控制效果評(píng)估

圖7為2006—2007年和2010—2011年兩個(gè)枯水期南四湖分區(qū)計(jì)算的平均綜合污染指數(shù)比較。

圖7 南四湖分區(qū)平均綜合污染指數(shù)比較

由圖7可見，2010—2011年比2006—2007年枯水期南四湖7個(gè)分湖區(qū)的水質(zhì)平均綜合污染指數(shù)全部處于下降狀態(tài)，下降率在3.1% ～57.1%之間，其中南四湖的主要納污湖區(qū)—南陽湖和薛城小沙河控制單元的影響湖區(qū)—微山湖島北的下降率均達(dá)到50%以上，反映南四湖主要入湖河流的污染物負(fù)荷得到有效控制。在南四湖7個(gè)分湖區(qū)中，2010—2011年枯水期小北湖和南陽湖的水質(zhì)最差，微山湖島南的水質(zhì)相對(duì)較好，總體呈現(xiàn)沿流向由北向南水質(zhì)逐漸變好的走勢(shì)，但在微山湖島北仍出現(xiàn)小的峰值湖區(qū)。這與南四湖入湖河流污染物輸入、地形與水文特征、水力輸移與污染物擴(kuò)散特征、水體富營養(yǎng)化特征以及受污染程度一致。

2010—2011年比2006—2007年枯水期南四湖7個(gè)分湖區(qū)水質(zhì)COD的下降率在34.1% ～55.5%之間，體現(xiàn)了全湖COD整體大幅度降低的趨勢(shì);在南陽湖和微山湖島北NH3-N的下降率分別為61.9%和80.4%，其他5個(gè)湖區(qū)的 NH3-N略有上升，主要是2010—2011年枯水期在采樣期間湖水溫較低、硝化作用較弱所致，但南四湖所有湖區(qū)的NH3-N均可達(dá)到Ⅲ類水標(biāo)準(zhǔn);TN的下降率在14.3% ～72.0%之間，其中南陽湖和微山湖島北的下降率均達(dá)到65%以上，說明南四湖流域工業(yè)和城市生活排水的氮素得到有效控制，特別是微山湖島北經(jīng)新薛河入湖的魯南化肥廠排水得到有效的治理;在南陽湖和獨(dú)山湖TP的下降率分別為55.0%和13.2%，其他5個(gè)湖區(qū)的TP均有所上升，上升率在16.4% ～49.3%之間，其中微山湖島北和微山湖島南的上升率在43.0%以上。一方面是排入南陽湖的主要磷污染河流洸府河等的污染物得到一定控制，另一方面是南四湖入湖河流總體的磷污染物控制相對(duì)較低，在南四湖7個(gè)分湖區(qū)中TP總體呈現(xiàn)與平均綜合污染指數(shù)相同的變化趨勢(shì)。

2010—2011年枯水期南四湖流域干旱少雨，天然徑流水量相對(duì)較小，湖內(nèi)中水比例增加，對(duì)水質(zhì)的影響可以通過南四湖水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率的關(guān)系分析。電導(dǎo)率指標(biāo)常用于推測(cè)水中離子的總濃度或含鹽量，其值越大則反映南四湖水體中來自流域的工業(yè)和城市生活排放的中水量所占比例越高。圖8中分別給出2006—2007年和2010—2011年兩個(gè)枯水期南四湖水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率的關(guān)系。

圖8 南四湖水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率的關(guān)系

由圖8可見，2010—2011年比2006—2007年枯水期南四湖水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率的相關(guān)性更好，南四湖各湖區(qū)電導(dǎo)率總體平均上升了279 μS/cm，說明南四湖內(nèi)中水比例提高，南四湖水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率存在很好的相關(guān)性，因此南四湖水污染物仍然以外源排放為主。值得注意的是2010—2011年比2006—2007年枯水期南四湖的中水比例提高了，而平均綜合污染指數(shù)下降了0.63，水質(zhì)仍然得到大幅度的改善，說明南四湖流域的河流入湖污染負(fù)荷大大降低，在南四湖流域GDP每年以兩位數(shù)增長和天然徑流水量相對(duì)下降的前提下，仍然實(shí)現(xiàn)了水質(zhì)的持續(xù)改善。

4 結(jié)論

a.南四湖水質(zhì)空間分布特征與對(duì)比分析結(jié)果表明，南四湖 COD、NH3-N、TP、TN 和 CODMn的空間分布具有北高南低的非均一性特點(diǎn)。2010—2011年比2006—2007年枯水期南四湖各水質(zhì)指標(biāo)的平均濃度下降率為38.3%，整體提高了1～2個(gè)水質(zhì)類別。

b.按平均綜合污染指數(shù)的評(píng)價(jià)結(jié)果，2010—2011年比2006—2007年枯水期南四湖污染和重污染湖區(qū)的比例大幅度減小，特別是重污染湖區(qū)下降到2.2%，達(dá)標(biāo)和基本達(dá)標(biāo)湖區(qū)的比例達(dá)到63.3%。2010—2011年枯水期南四湖主要污染物按分擔(dān)率由大到小排序依次為 TP、COD、TN、CODMn、NH3-N。

c.按平均綜合污染指數(shù)的評(píng)價(jià)結(jié)果，在南四湖7個(gè)分湖區(qū)中，2010—2011年枯水期小北湖和南陽湖的水質(zhì)最差，微山湖島南的水質(zhì)相對(duì)較好，總體呈現(xiàn)沿流向由北向南水質(zhì)逐漸變好的走勢(shì)，但在微山湖島北仍出現(xiàn)小的峰值湖區(qū)。

d.在南四湖同一次監(jiān)測(cè)中，水質(zhì)平均綜合污染指數(shù)與電導(dǎo)率存在很好的相關(guān)性和北高南低的污染特點(diǎn)，說明南四湖水污染物仍然以外源排放為主。

致謝:南四湖水質(zhì)空間分布監(jiān)測(cè)工作由山東省環(huán)境保護(hù)廳組織，得到美商生化科技公司的技術(shù)支持，監(jiān)測(cè)計(jì)劃由青島理工大學(xué)和山東大學(xué)共同編制，山東省環(huán)境監(jiān)測(cè)中心站、濟(jì)寧市環(huán)境保護(hù)監(jiān)測(cè)站負(fù)責(zé)實(shí)施，在此一并表示感謝!

[1]劉洪林，邢文潔，王素芬，等.南水北調(diào)東線輸水干線南四湖水環(huán)境質(zhì)量研究[J].水文，2003，23(4):49-51.

[2]陳磊，裴海燕，解軍.改善南四湖水質(zhì)的關(guān)鍵問題分析[J].中國給水排水，2007，23(20):6-10.

[3]梁春玲，張祖陸，牛玉生.南四湖水安全評(píng)價(jià)研究[J].安徽農(nóng)業(yè)科學(xué)，2009，37(7):3178-3181，3307.

[4]孫秀敏，王琳，甘泓，等.中水濕地去氮效果研究[J].水資源保護(hù)，2010，26(2):10-12.

[5]武周虎，慕金波，謝剛，等.南四湖及入出湖河流水環(huán)境質(zhì)量變化趨勢(shì)分析[J].環(huán)境科學(xué)研究，2010，23(9):1167-1173.

[6]WU Zhou-hu，ZHANG Jian，ZHU Jie，et al.A monitoring project planning technique of the water quality spatial distribution in Nansi Lake[J].Procedia Environmental Sciences，2011，10:2320-2328.

[7]周豐，郭懷成，黃凱，等.基于多元統(tǒng)計(jì)方法的河流水質(zhì)空間分析[J].水科學(xué)進(jìn)展，2007，18(4):544-551.

[8]劉碩，朱建平，蔣火華.對(duì)幾種環(huán)境質(zhì)量綜合指數(shù)評(píng)價(jià)方法的探討[J].中國環(huán)境監(jiān)測(cè)，1999，15(5):33-37.

[9]范瑜.平均綜合污染指數(shù)在地面水污染分級(jí)中的應(yīng)用[D].上海:中國科學(xué)院上海冶金研究所(今中國科學(xué)院上海微系統(tǒng)與信息技術(shù)研究所)，2000.

[10]鄧建才，劉恩峰，胡維平，等.南四湖水體氮磷及葉綠素空間分布特征[J].灌溉排水學(xué)報(bào)，2010，29(4):119-122.

[11]KEMKA N，NJINé T，ZéBAZé TOGOUET S H，et al.Eutrophication of lakes in urbanized areas:the case of Yaounde Municipal Lake in Cameroon，Central Africa[J].Lakes ＆ Reservoirs:Research and Management，2006(11):47-55.

[12]武周虎，張建，路成剛，等.南四湖水質(zhì)空間分布監(jiān)測(cè)分析與水環(huán)境問題解析[J].長江流域資源與環(huán)境，2011，20(增刊1):137-143.