大汶河梯級攔蓄水庫富營養(yǎng)化評估

何濤 夏霆 陳懷民 趙計奎 武欽凱

摘 要：

基于大汶河沿線4個監(jiān)測點2007—2014年水質(zhì)資料，分析溶解氧、化學(xué)需氧量、氨氮等9項指標的年際變化特征，采用灰色模式識別模型評價大汶河水質(zhì)狀況。通過9個梯級水庫2015年豐、平、枯水期水質(zhì)現(xiàn)場調(diào)查，分析氮、磷營養(yǎng)鹽時空變化特征，運用加權(quán)綜合營養(yǎng)狀態(tài)指數(shù)法評價梯級水庫營養(yǎng)狀態(tài)。結(jié)果表明：2007—2014年間，受流域內(nèi)各項減排措施影響，大汶河有機污染得到有效緩解；受攔蓄工程影響，河段內(nèi)TP明顯變優(yōu)，TN總體變差；灰色綜合指數(shù)表明，水質(zhì)狀況整體有所好轉(zhuǎn)，且季節(jié)波動趨于平穩(wěn)。因污染程度不同，梯級水庫富營養(yǎng)化狀態(tài)有明顯差異，部分庫區(qū)豐水期和枯水期已分別爆發(fā)藍藻和硅藻水華，不同水期因下泄流量不同等造成流動條件的差異，對庫區(qū)水質(zhì)有較大影響，水污染程度、工程攔蓄與運行造成的流動條件差異共同對庫區(qū)內(nèi)藻類生長與爆發(fā)產(chǎn)生重要影響。

關(guān)鍵詞：大汶河；水質(zhì)變化；梯級攔蓄工程；富營養(yǎng)化

中圖分類號：X524

文獻標志碼：A 文章編號：16744764（2018）05013308

收稿日期：20171123

基金項目：

泰安市水利局科技項目

作者簡介：

何濤（1993），男，主要從事水污染控制研究，Email： hetao226@njtech.edu.cn。

夏霆（通信作者），男，博士，副教授，Email： xiating@njtech.edu.cn。

Received：20171123

Foundation item：

Science and Technology Project of Taian Water Resource Bureau

Author brief：

He Tao（1993），main research interest：water pollution control，Email：hetao226@njtech.edu.cn.

Xia Ting（corresponding author）， PhD，associate professor，Email：xiating@njtech.edu.cn.

Eutrophication assessment in cascade reservoirs of Dawen River

He Tao1， Xia Ting1， Chen huaimin2， Zhao Jikui1， Wu Qinkai1

（1.College of Urban Construction， Nanjing University of Technology， Nanjing 211816，P.R.China；

2.College of Environmental Science and Engineering， Hohai University， Nanjing 210098，P.R.China）

Abstract：

Based on water quality data of 4 stations in the Dawen River from 2007 to 2014， this paper analyzed the interannual variability of 9 indicators such as dissolved oxygen， chemical oxygen demand， ammonia nitrogen etc.， and evaluated water quality by a Grey Model. By the site investigation of 9 cascade reservoirs during the rich， average and dry period in 2015， the spatiotemporal variation of nitrogen and phosphorus nutrients were analyzed， nutritional level was assessed by TLI.The results show that during the period of 2007 to 2014， the organic pollution have been effectively alleviated by the emission reduction measures in the basin， the concentrations of TP have been meliorated obviously while the concentrations of TN have been aggravated integrally by the construction of cascade reservoirs. The grey comprehensive index shows that the water quality of the Dawen River is improved and the seasonal fluctuations tend to be stable. The state of eutrophication has apparent difference among cascade reservoirs because the contamination vary in degree. Part of the reservoirs have been experienced cyanobacteria and diatom blooms during rainy and dry seasons. The different flow conditions induced by different discharge in corresponding season have great influence on water quality of the reservoirs. ， the difference of flow condition caused by the water pollution level， damming and impoundment have a significant impact on the growth and bloom of algae in reservoirs.

Keywords：

Dawen River； water quality variety； cascade reservoir； eutrophication

大汶河是黃河下游、山東省境內(nèi)的最大支流，也是泰安市唯一的大型河道，是泰安市的母親河[1]，但流域水資源短缺且分布嚴重不均。近年來，為緩解大汶河流域水資源供需矛盾，充分合理地開發(fā)利用大汶河地表徑流，在干流和兩大支流上規(guī)劃新建和修復(fù)一批攔河樞紐，以保證流域用水需求，減少洪澇與干旱災(zāi)害，提高水資源調(diào)配能力。但河道攔蓄工程被認為具有阻隔、截留和累積等效應(yīng)[24]，對河流生態(tài)環(huán)境產(chǎn)生一系列的影響：因改變河流自然水文模式，破壞河流自然“連續(xù)性”[56]，河流生態(tài)系統(tǒng)因此而向“靜水”生態(tài)系統(tǒng)轉(zhuǎn)變，強水動力下的河流搬運作用逐漸演變成“湖泊”沉積作用[7]，營養(yǎng)物質(zhì)容易富集，并使得攔蓄水庫易出現(xiàn)藻類密度大量增長，甚至水華爆發(fā)等現(xiàn)象。

如何協(xié)調(diào)好河流社會服務(wù)與自然生態(tài)環(huán)境保護之間的關(guān)系，一直是河流保護的關(guān)鍵問題，作為南水北調(diào)東線工程重要蓄水走廊——東平湖的惟一匯入河流，大汶河水質(zhì)狀況受到社會各界廣泛關(guān)注。選擇泰安市境內(nèi)大汶河干流、主要支流牟汶河以及河段上新建和修復(fù)的9座攔河水庫為研究對象，基于歷史水質(zhì)監(jiān)測數(shù)據(jù)以及現(xiàn)場監(jiān)測與調(diào)研，分析大汶河水質(zhì)時空變化和梯級水庫富營養(yǎng)化狀態(tài)差異，探討大汶河水質(zhì)狀況變化的影響因素與攔蓄工程對大汶河梯級水庫營養(yǎng)狀況的影響，并提出相關(guān)建議，以促進大汶河水資源開發(fā)利用與生態(tài)健康保護。

1 材料與方法

1.1 研究區(qū)概況

大汶河（北緯35°59′55.6″～36°7′10.1″、東經(jīng)116°12′17.6″～117°56′29.4″）古稱汶水，發(fā)源于泰萊山區(qū)，匯集泰山山脈以南、蒙山山脈以北諸水，自東向西流經(jīng)萊蕪、新泰、泰安、肥城、寧陽、汶上、東平等縣、市，又經(jīng)東平湖清河門出湖閘入黃河，全長231 km，總落差362 m。泰安市境內(nèi)大汶河流域面積6 457.2 km2，占全流域面積比72.2%。泰安屬水資源短缺地區(qū)，水資源年均占有量僅為311 m3，不足全國平均水平的1/7。大汶河流域降雨量年內(nèi)分布嚴重不均，2006—2015年間，泰安大汶河區(qū)平均降雨量為670.22 mm，汛期（6—9月）平均降雨量為494.92 mm，占全年降水量的73.84%，豐、枯水年交替出現(xiàn)，且近10 a來降水量總體有下降趨勢（圖1）。因汛期暴雨集中，而其余時段干旱少雨，大汶河地區(qū)洪澇、旱災(zāi)交替頻發(fā)。以2012年為例，汛期泰安市因洪澇災(zāi)害造成15.03萬人受災(zāi)，受災(zāi)面積8.199千hm2，絕收面積2.2千hm2，而同年又因干旱受災(zāi)面積51 552 hm2，絕收面積達7 519 hm2[8]。

為調(diào)洪補枯、應(yīng)對水資源短缺與洪澇災(zāi)害問題，2010—2014年間，泰安市在大汶河干流和主要支流上新建或修復(fù)了一批攔蓄工程。為盡量減少對大汶河生態(tài)系統(tǒng)的影響，攔蓄工程多采用橡膠壩等低壩形式。目前，關(guān)于梯級攔蓄河流的水質(zhì)評價多圍繞水電開發(fā)等形成的深水大庫開展[9]，尚未見有針對低壩形式開發(fā)的調(diào)研與評估。為深入分析梯級攔蓄工程對大汶河水質(zhì)狀況的影響，此次評估以唐莊（TZ）、顏張（YZ）、泉林（QL）、顏謝（YX）、汶口1號（WK1）、汶口2號（WK2）、堽城（GC）、琵琶山（PPS）和戴村（DC）9個梯級水庫為研究對象。

1.2 數(shù)據(jù)來源與分析

收集與分析2007—2014年大汶河沿線4個水文站（萊蕪、北望、大汶口和戴村壩）奇數(shù)月份的水質(zhì)監(jiān)測數(shù)據(jù)，水質(zhì)指標選?。喝芙庋酰―O）、pH值、氟化物（F-）、化學(xué)需氧量（COD）、高錳酸鹽指數(shù)（IMn）、五日生化需氧量（BOD5）、總氮（TN）、總磷（TP）和氨氮（NH3—N）。

針對大汶河干流和支流牟汶河上的9個梯級水庫，分別于2015年1月（枯水期）、5月（平水期）和7月（豐水期）進行采樣分析，樣品采集、保存和分析按照文獻[10]進行。調(diào)查指標包括：pH值、溶解氧、透明度（SD）、總氮、氨氮、硝態(tài)氮（NO-3—N）、高錳酸鹽指數(shù)、總?cè)芙庑缘═DN）、總磷、總?cè)芙庑粤祝═DP）、葉綠素（Chla）和浮游植物等。

1.3 評價方法

水質(zhì)狀況評價參照文獻[11]運用的水質(zhì)評價灰色模式識別模型：1）基于《地表水環(huán)境質(zhì)量標準》（GB 3838—2002）建立參考序列及比較序列；2）數(shù)據(jù)歸一化處理；3）關(guān)聯(lián)度、關(guān)聯(lián)離散度、隸屬度計算；4）水質(zhì)灰色識別模式綜合指數(shù)（GC）計算。

水庫富營養(yǎng)評估參照中國環(huán)境監(jiān)測總站推薦的加權(quán)綜合營養(yǎng)狀態(tài)指數(shù)法，計算式為

2 結(jié)果與分析

2.1 大汶河水質(zhì)狀況年際變化

大汶河水質(zhì)因子年際變化如路圖2所示。DO是衡量河流自凈能力及能否維持生態(tài)平衡的重要指標。2007—2014年間，大汶河DO平均濃度變化范圍為6.9～14.95 mg/L，平均值為8.99 mg/L。2010年后穩(wěn)定在一個較高水平，超地表水環(huán)境質(zhì)量Ⅰ類標準（GB 3828—2002）。pH值反映水體酸堿強度與氫離子濃度，決定眾多物化生過程。2007—2014年間，大汶河pH值變化范圍為7.18～7.90 mg/L，平均值為7.50 mg/L，呈弱堿性。2010年后下降了約0.5個單位，整體變化趨勢具有明顯的趨同性。F廣泛分布于自然水體，是生物體必須的微量元素之一，過量則會抑制植物的光合作用。2007—2014年間，大汶河F平均濃度變化范圍為036～0.98 mg/L，平均值為0.57 mg/L。處于Ⅰ類水質(zhì)標準，且2010年后變化較為平穩(wěn)，各站點水平相當。

COD、IMn和BOD5綜合反映地表水體受有機污染的程度。2007—2014年間，大汶河COD平均濃度變化范圍為19.45～50.72 mg/L，平均值為2887 mg/L。2007—2011年間整體上在Ⅳ～Ⅴ類標準波動，2012年后穩(wěn)定在Ⅳ類標準。IMn平均濃度變化范圍為3.78～11.85 mg/L，平均值為5.69 mg/L。2007—2011年間整體上呈現(xiàn)顯著下降趨勢，尤其萊蕪與大汶口站，2012年后穩(wěn)定在Ⅳ類標準。BOD5代表水體可生物降解有機物含量，2007—2014年間平均濃度變化范圍為2.27～7.69 mg/L，平均值為5.69 mg/L。2007—2011年間，上游萊蕪站穩(wěn)定在Ⅴ類標準，其余3站呈現(xiàn)顯著上升趨勢。2012年后除下游戴村站處于Ⅲ類標準外，整體上均穩(wěn)定在Ⅳ類標準。

TN、TP和NH3—N是表征水體富營養(yǎng)化狀態(tài)的重要指標。2007—2014年間，大汶河TN平均濃度變化范圍為7.335～16.68 mg/L，平均值為1149 mg/L。超地表水環(huán)境質(zhì)量標準中Ⅴ類標準的限值，達劣Ⅴ類。萊蕪與戴村站變化較為平穩(wěn)且濃度相對較低，北望、大汶口2010年后處于較高水平且存在較為明顯的上升，2014年出現(xiàn)穩(wěn)定趨勢。NH3—N平均濃度變化范圍為0.65～4.54 mg/L，平均值為1.48 mg/L。除大汶口站超Ⅴ類標準限值外，整體在在Ⅳ～Ⅴ類標準波動。TP平均濃度變化范圍為0.02～0.39 mg/L，平均值為0.16 mg/L。除個別年度超Ⅴ類標準外，總體上在Ⅲ～Ⅳ類標準波動。2012年后下游3站點均穩(wěn)定在較低水平（003 mg/L）。

根據(jù)大汶河2007—2014年48個奇數(shù)月DO、pH、F-、COD、IMn、BOD5、TN、TP和NH3—N 9項指標監(jiān)測數(shù)據(jù)，利用灰色模式識別模型對大汶河水質(zhì)狀況進行綜合分析與評價。通過主成份分析（KMO檢驗值=0.675）確定各參評指標權(quán)重，權(quán)向量W=0.039、0.040、0.085、0.090、0.134、0.082、0.147、0160、0.223。計算水質(zhì)灰色識別模式綜合指數(shù)，結(jié)果如圖3所示，基于GC值，2007—2014年間大汶河水質(zhì)轉(zhuǎn)好趨勢明顯，且2012年后季節(jié)性波動變小，介于2～3之間。整體上大汶河水質(zhì)變化趨勢與上游來水（萊蕪站）具有明顯趨同性，上、中游（萊蕪、北望、大汶口站）水平相當，下游（戴村壩站）相對較低。

2.2 梯級水庫營養(yǎng)狀況時空變化

沿河筑壩蓄水形成梯級水庫群，改變原有自然水文、水動力模式，河流自凈能力下降，營養(yǎng)物質(zhì)易于在庫區(qū)富集，富營養(yǎng)化問題嚴峻。調(diào)查期梯級水庫氮、磷營養(yǎng)鹽時空變化如圖4所示。大汶河梯級水庫TN、NH3—N和NO-3—N濃度沿程變化趨勢基本一致，總體呈波浪式變化且沿程降低，在泉林、堽城處到達波峰，其中，泉林處達最大值，分別為1637、2.98和8.24 mg/L。泉林庫區(qū)主要受重污染支流泮汶河匯入影響，后者承納市區(qū)工業(yè)廢水及40多萬城鎮(zhèn)居民生活污水，水質(zhì)超地表水Ⅴ類標準（GB 3838—2002）。季節(jié)變化方面，TN濃度枯水期最高，豐水期最低，表明庫區(qū)TN污染可能主要來自沿途城鎮(zhèn)和生活污水排放等點源污染，豐水期稀釋作用加強；NH3—N濃度除受泮汶河影響較大的泉林庫區(qū)外，均為豐水期最高，說明NH3—N污染主要來源于沿河區(qū)域面源污染，同時，可能經(jīng)有機氮微生物在夏季有氧條件下分解形成；NO-3—N濃度枯水期最高，豐、平水期基本一致，主要由于枯水期NO-3—N污染主要來自外源化肥輸入，豐、平水期則以生活污水和內(nèi)部生物作用為主。TP濃度整體上在泉林庫區(qū)達最大值，約0.59 mg/L。豐水期相對較高，自顏謝下泄后TP明顯偏低，均低于0.01 mg/L，表明梯級攔蓄對TP的攔截與沉積作用明顯。對比上、下游梯級水庫營養(yǎng)鹽變化，并無明顯的累積性效應(yīng)；而泉林等水庫營養(yǎng)鹽的峰值性凸顯變差，應(yīng)主要受污染匯入的影響，如泉林水庫水質(zhì)應(yīng)與污染嚴重的泮汶河支流匯入密切相關(guān)。

選取TP、TN、CODMn、SD和Chla 5項水質(zhì)因子，運用加權(quán)綜合營養(yǎng)狀態(tài)指數(shù)法評價大汶河梯級庫區(qū)營養(yǎng)狀態(tài)。調(diào)查期大汶河各梯級庫區(qū)綜合營養(yǎng)指數(shù)TLI（∑）沿程變化范圍見圖5，豐、平、枯水期TLI（∑）值均在泉林出現(xiàn)最大峰值（分別為61.54、63.51和68.76），處于中度富營養(yǎng)狀態(tài)，自泉林往下游，各水庫富營養(yǎng)狀況總體有逐級減緩趨勢，且不同水期具有較為明顯的趨同性，豐、平水期戴村水庫已為貧營養(yǎng)狀態(tài)。從富營養(yǎng)化評價結(jié)果來看，主要因污染程度不同，梯級水庫富營養(yǎng)化狀態(tài)差異明顯，泉林水庫污染嚴重，富營養(yǎng)程度嚴重。

3 討論

3.1 大汶河水質(zhì)狀況變化的影響因素

基于2007—2014年水質(zhì)指標主成分分析，主要水質(zhì)污染指標依次為NH3—N、TP、TN和IMn（權(quán)重>0.1）。馬吉讓等[13]調(diào)查大汶河水系30個入河排污口，污染物年入河總量為31 957萬t/a，泰安區(qū)域占比達77.1%，主要污染物COD為19 660 t/a，NH3—N為3 155 t/a，TP為326 t/a，BOD5為2 481 t/a，TN為4 390 t/a，大汶河水體污染物主要是營養(yǎng)鹽和有機物，與本文分析結(jié)果基本一致。近年來，沿河地區(qū)大力推進各項減排措施，尤其萊鋼減排貢獻突出，上游來水水質(zhì)明顯改善，大汶河整體水質(zhì)狀況有所好轉(zhuǎn)。隨著城市化與農(nóng)業(yè)集約化進程加快，沿河生活與農(nóng)業(yè)污染源成為大汶河水體污染的重要來源，2014年生活污水排放量高達 18 966.59萬t，且存在逐年上升趨勢，農(nóng)業(yè)化肥施用量65.70萬t（氮肥26.17萬t，磷肥9.99萬t）。因此，大汶河水環(huán)境治理首要解決沿河污染匯入問題。

調(diào)查期NH3—N和COD濃度整體呈下降趨勢，主要由于近年來泰安市不斷加強對污染排放的監(jiān)管，整治城區(qū)污水管線，提高污水處理能力，截至2014年全市建成污水處理廠12座，大大削減了入河污染，使得河道內(nèi)NH3—N滿足地表水Ⅲ類標準，COD滿足Ⅳ類標準。BOD5濃度2010年后明顯下降，應(yīng)該是周邊污染排放得到有效控制，同時，也可能與梯級水庫建設(shè)后河段內(nèi)新增大面積濕地，水質(zhì)凈化能力提高。受梯級攔蓄工程影響，河段TP濃度2012年后顯著下降，參照其他地區(qū)攔蓄工程的影響分析[14]，這應(yīng)與梯級攔蓄水庫對TP的攔截與沉積作用有關(guān)。與此同時，DO濃度趨于穩(wěn)定化，由于DO受徑流影響較大，這應(yīng)與攔蓄工程調(diào)節(jié)影響使得河道徑流較為穩(wěn)定有關(guān)。

3.2 梯級攔蓄對大汶河富營養(yǎng)化影響

梯級攔蓄影響大汶河氮磷營養(yǎng)鹽的遷移、轉(zhuǎn)化等生物地球化學(xué)行為，個別庫區(qū)富營養(yǎng)化問題嚴峻，但整體上并無明顯的累積效應(yīng)。經(jīng)采樣分析（見圖6），大汶河梯級水庫群調(diào)查期TDN/TN值高達73%，表明庫區(qū)N污染以溶解態(tài)為主，而溶解態(tài)N較顆粒態(tài)受梯級水庫主要的沉積作用影響微弱，因此，TN濃度在庫區(qū)內(nèi)迅速累積，進而導(dǎo)致全河性的嚴重超標。攔蓄工程因加速TP的攔截與沉積而對庫區(qū)與下泄水體TP濃度具有明顯的改善作用。同時，張恩仁等[3]認為發(fā)育出的水庫生態(tài)系統(tǒng)可將上游輸入的13%～42%的溶解態(tài)無機磷固定于浮游生物。大汶河全年TDP/TP值僅為28%，表明大汶河P污染主要以顆粒態(tài)為主，而修建梯級水庫產(chǎn)生阻隔效應(yīng)，水體流動性減小，加速懸浮態(tài)P在庫區(qū)沉積，降低水體中TP濃度。豐水期TP濃度明顯高于平、枯水期，主要由于豐水期庫水交換量大，TP的攔截與沉積作用相對較弱。基于P在庫區(qū)的積蓄，大汶河攔蓄水庫尤其是上游水庫應(yīng)注意防范底泥內(nèi)源P釋放對上覆水質(zhì)的影響。

因大汶河旱澇不均，攔蓄工程建設(shè)之前，非汛期河道內(nèi)長期斷流，大段河道枯竭，河流生態(tài)健康受損嚴重。攔蓄工程的建設(shè)盡管使得局部水流受阻，但低壩的建設(shè)形式使得河道內(nèi)形成多個首尾相連的淺水水庫，對沿河水資源利用、河道內(nèi)水生生物棲息地修復(fù)與河道健康維護作用巨大。根據(jù)水庫運行水量統(tǒng)計（見圖7），2014年研究區(qū)9個梯級水庫蓄水量達6 680萬m3，改變原有枯水期斷流現(xiàn)象。并使得河道內(nèi)新增河灘型濕地2 885.4 ha，顏張壩下至泉林庫區(qū)段約30 km河段內(nèi)已形成著名的徂徠山汶河自然濕地景觀。現(xiàn)有關(guān)于自然濕地緩沖帶對氮磷去除的研究表明[1516]，濕地對氮、磷的削減作用顯著，生長季節(jié)的濕地緩沖帶對面源TN、TP去除率高達74.1%和84.6%。同時，濕地對河道內(nèi)N、P營養(yǎng)鹽有良好的吸收、利用和轉(zhuǎn)化作用，能有效促進硝化和反硝化，改善水體富營養(yǎng)狀態(tài)。

在外來污染源未得到有效控制條件下，攔蓄工程阻隔使得水體交換受阻，水庫N、P營養(yǎng)鹽蓄積，致使藻類快速生長與爆發(fā)。同步的浮游植物取樣調(diào)研發(fā)現(xiàn)，泉林庫區(qū)夏季豐水期已爆發(fā)嚴重的藍藻（微囊藻）水華，冬季枯水期泉林和堽城庫區(qū)出現(xiàn)了硅藻水華（小環(huán)藻）。

水體Chla含量是表征藻類現(xiàn)存量的重要指標[17]，相關(guān)分析表明，調(diào)查期水庫Chla濃度（見圖8）與綜合營養(yǎng)指數(shù)TLI（∑）顯著相關(guān)（相關(guān)系數(shù)0695**）。進一步分析不同期Chla含量與各水質(zhì)因子的相關(guān)關(guān)系（表1）表明：豐水期Chla濃度與各水質(zhì)因子無明顯相關(guān)性，應(yīng)是夏季藻類生長受光照、溫度和水體滯留狀態(tài)等因素影響較大，且營養(yǎng)鹽本底值較高，營養(yǎng)鹽對藻類生長難以形成限制因素。平、枯水期Chla與TP和NH3N濃度有著極顯著的相關(guān)關(guān)系，應(yīng)是梯級水庫對P的攔截與沉積作用明顯，水體中P含量較低，成為主要限制因素。Redfield等[18]認為浮游植物光合作用中吸收適宜的N/P比為16，當N/P>7時浮游植物生長受磷限制。因平、枯水期河道及水庫內(nèi)TN高，庫內(nèi)TP偏低，水庫富營養(yǎng)化及藻類生長均受P影響大。并且，盡管庫區(qū)TN濃度高，但NH3—N/TN值總體為012，NH3—N作為藻類吸收利用的重要氮源，也成為影響藻類生長的重要因素?？菟贑hla濃度與WT有著顯著的相關(guān)關(guān)系，主要由于枯水期大汶河浮游植物優(yōu)勢門類為硅藻，對光強和溫度的要求相比其他藻類低，能適應(yīng)低溫條件，在營養(yǎng)鹽適宜和緩流條件下以累積生物量的方式快速對水溫作出響應(yīng)[19]。從不同水期Chla分布及其與水質(zhì)因子的關(guān)系來看，大汶河梯級水庫富營養(yǎng)化及藻類爆發(fā)除主要受N、P營養(yǎng)鹽影響外，亦受光照、水溫、水流狀態(tài)等其他理化因素影響較為明顯。

4 結(jié)論

1）2007—2014年間，大汶河水體主要污染物為營養(yǎng)鹽和有機物，受沿河地區(qū)各項減排措施影響，NH3—N和COD濃度整體呈下降趨勢，一定程度上緩解了水體有機污染；受梯級攔蓄工程影響，TN濃度逐年提高尤其城市化程度較高區(qū)域，TP濃度2012年后穩(wěn)定在較低水平。水質(zhì)灰色識別模式綜合指數(shù)表明，大汶河水體環(huán)境正向良性方向發(fā)展，水質(zhì)狀況得到明顯改善，且2012年后季節(jié)波動趨于平穩(wěn)。

2）對9個梯級水庫不同水期水質(zhì)調(diào)查表明，除部分庫區(qū)因污染嚴重外，梯級水庫水質(zhì)總體呈現(xiàn)由上游至下游逐步變優(yōu)的趨勢；不同水期因下泄流量不同等造成流動條件的差異，對庫區(qū)水質(zhì)有較大影響；因污染程度不同，梯級水庫富營養(yǎng)化狀態(tài)有明顯差異；水污染程度、工程攔蓄對庫區(qū)內(nèi)藻類生長與爆發(fā)共同產(chǎn)生影響，部分庫區(qū)豐水期和枯水期已分別爆發(fā)藍藻和硅藻水華。

3）梯級攔蓄工程因蓄積水資源、恢復(fù)河灘濕地，改善了河道內(nèi)長期斷流現(xiàn)象，并促進河道生物棲息條件好轉(zhuǎn)。但為減緩水庫富營養(yǎng)化趨勢及防控藻類水華爆發(fā)，應(yīng)嚴控污染源匯入，并優(yōu)化重點庫區(qū)和時段的出、入庫流量，改善水庫流動與營養(yǎng)鹽蓄積條件。

參考文獻：

[1] 趙志鵬，劉慧，周震. 泰安市生態(tài)河道建設(shè)措施[J].山東水利，2010（9）：5152.

ZHAO Z P， LIU H， ZHOU Z. Measures ofecological river construction in Taian [J]. Shandong Water Resources，2010（9）：5152.（in Chinese）

[2] WARD J V，STANFORD J A.Ecological connectivity in alluvial river ecosystems and its disruption by flow regulation [J]. Regulated Rivers Research and Management，1995，11（1）：105119.

[3] 張恩仁，張經(jīng). 三峽水庫對長江N、P營養(yǎng)鹽截留效應(yīng)的模型分析[J].湖泊科學(xué)，2003， 15（1）：4148.

ZHANG E R， ZHANG J. Analysis of the threegorge reservoir impactson the retention of N and P in the Yangtze River[J].Journal of Lake Sciences，2003， 15（1）：4148. （in Chinese）

[4] 鐘華平，劉恒，耿雷華.瀾滄江流域梯級開發(fā)的生態(tài)環(huán)境累積效應(yīng)[J].水利學(xué)報，2007（Sup1）：582586.

ZHONG H P，LIUH，GENG L H. Cumulative effects of Lancang River basin cascade hydropowerdevelopment on ecology and environment [J].Journal of Hydraulic Engineering，2007（Sup1）：582586. （in Chinese）

[5] TUNDISI J G， STRASKRABA M. Theoretical reservoir ecology and its applications [J].Optical Remote Sensing of Air Pollution，1999：95122. （in Chinese）

[6] 申恒倫，蔡慶華，邵美玲，等. 三峽水庫香溪河流域梯級水庫浮游植物群落結(jié)構(gòu)特征[J].湖泊科學(xué)，2012， 24（2）：197205.

SHEN H L， CAI Q H， SHAO M L， et al. Community structure characteristics of phytoplankton along a cascade of reservoirs in Xiangxi River basin， Three Gorges Reservior [J]. Journal of Lake Sciences，2012， 24（2）：197205. （in Chinese）

[7] HART D D，POFF N L. A special section on dam removal and river restoration[J].Bioscience，2002，52（8）：653655.

[8] 泰安市大汶河志編纂委員會. 大汶河志[M].北京：方志出版社，2016：9399.

Compilation Committee of Taian Dawen River Records. Dawen River records [M].Beijing：Fangzhi Press，2016：9399. （in Chinese）

[9] 邱光勝，胡圣，葉丹，等.三峽庫區(qū)支流富營養(yǎng)化及水華現(xiàn)狀研究[J].長江流域資源與環(huán)境，2011， 20（3）：311316.

QIU G S， HU S， YE D， et al. Investigation on the present situation of eutrophication and water bloom in the branches of three gorges reservoir[J].Resources and Environment in the Yangtze Basin，2011， 20（3）：311316. （in Chinese）

[10] 國家環(huán)境保護總局. 水和廢水監(jiān)測分析方法[M].北京：中國環(huán)境科學(xué)出版社，2002：649653.

The State Environmental Protection Administration. Water and wastewater monitoring and analysis method [M].Beijing： China Environmental Science Press， 2002：649653. （in Chinese）

[11] 田偉東，賈克力，史小紅，等. 2005—2014年烏梁素海湖泊水質(zhì)變化特征[J].湖泊科學(xué)，2016（6）：12261234.

TIAN W D，JIA K L，SHI X H， et al.Water quality variation in Lake Wuliangsuhai， 20052014 [J]. Journal of Lake Sciences，2016，28（6）：12261234. （in Chinese）

[12] 金相燦，劉樹坤，章宗涉，等. 中國湖泊環(huán)境[M].北京：海洋出版社，1995：150197.

JIN X C， LIU S K， ZHANG Z S， et al. China lake environment [M].Beijing：Ocean Press， 1995：150197. （in Chinese）

[13] 馬吉讓. 大汶河流域排污口設(shè)置現(xiàn)狀調(diào)查評價與對策建議[C]//全國河湖治理與水生態(tài)文明發(fā)展論壇，2015.

MA J R. The investigation and evaluation of the status quo of sewage outlet in Dawen River basin and its countermeasures [C]//Nationalrivers and Lakesregulation and Water Ecological Civilization Development Forum，2015. （in Chinese）

[14] 黎慧卉，劉叢強，汪福順，等. 貓?zhí)恿饔蛱菁壦畮炝椎南募咀兓卣鱗J].長江流域資源與環(huán)境，2009， 18（4）：368372.

LI H H，LIU C Q，WANG F S， et al. Change characteristics of phosphorus in cascade reservoirs on Maotiao River system in summer [J].Resources and Environment in The Yangtze Basin，2009， 18（4）：368372. （in Chinese）

[15] 王磊，章光新. 濕地緩沖帶對氮磷營養(yǎng)元素的去除研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2006，25：649652.

WANG L， ZHANG G X. Removal of nitrogen and phosphorus in riparian buffer zone of wetland [J].Journal of AgroEnvironment Science，2006，25：649652. （in Chinese）

[16] PETERJOHN W T， CORRELL D L. Nutrient dynamics in an agricultural watershed：Observations on the role of a riparian forest [J].Ecology， 1984， 65（5）：14661475.

[17] 國家環(huán)境保護局科技標準司. 中國湖泊富營養(yǎng)化及其防治研究[M].北京：中國環(huán)境科學(xué)出版社，2001.

Department of Science，Technology and Standards. Study on lake eutrophication and its countermeasure in China [M].Beijing：China Environmental Science Press，2001. （in Chinese）

[18] REDFIELD A C，KETCHUM B H. The influence of organisms on the composition of seawater [J].Sea，1963，40（6）：640644.

[19] 鄧建明，秦伯強. 全球變暖對淡水湖泊浮游植物影響研究進展[J].湖泊科學(xué)，2015，27（1）：110.

DENG J M， QIN B Q. A review on studies of effects of climate change on phytoplankton in freshwater systems [J].Journal of Lake Sciences，2015，27（1）：110. （in Chinese）

（編輯 王秀玲）