王 剛 齊 珺 潘 濤 宋新山 李 夏 張?jiān)砾i(1.北京市環(huán)境保護(hù)科學(xué)研究院，國家城市環(huán)境污染控制工程技術(shù)研究中心，北京 100037；.東華大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 0160；3.天津大學(xué)環(huán)境科學(xué)與工程學(xué)院，天津 30007；.首都師范大學(xué)資源環(huán)境與旅游學(xué)院，北京 10008)

目前，北京市水污染治理仍處于目標(biāo)總量控制階段。水污染以耗氧污染為主，水質(zhì)目標(biāo)為控制耗氧污染、消除黑臭和改善水質(zhì)，氨氮、COD為優(yōu)先控制污染物，也是國家環(huán)境保護(hù)“十五”、“十一五”和“十二五”計(jì)劃中重點(diǎn)確定的主要污染物總量控制因子[1-4]。“十五”以來，北京市對重點(diǎn)污染物實(shí)施總量控制，加大對城鎮(zhèn)污水處理廠建設(shè)的投資力度，開展河道整治工作，河流水質(zhì)達(dá)標(biāo)率整體呈現(xiàn)上升趨勢。但在2008年北京奧運(yùn)以后，隨著城市化進(jìn)程加快，治理成本和難度增大，水污染減排空間大幅減小，水質(zhì)達(dá)標(biāo)率徘徊不前，水質(zhì)改善進(jìn)入平臺(tái)期[5-6]。因此，從流域治理的角度，制定水污染防治措施方案，達(dá)到水質(zhì)改善效果，是當(dāng)前著重關(guān)心的問題[7-9]。河流水質(zhì)模型可以量化評估不同減排措施方案對于改善水質(zhì)的效益。由丹麥水力研究所(DHI)研發(fā)的MIKE11模型可用于模擬河流、河口、河網(wǎng)系統(tǒng)的水流、水質(zhì)過程，其應(yīng)用簡單、模擬準(zhǔn)確，尤其適用于多閘壩河流[10]，在水質(zhì)模擬預(yù)測、水污染治理效果評估中有著廣泛的應(yīng)用[11-14]。

本研究以2010年為基準(zhǔn)年，綜合考慮《北京市“十二五”時(shí)期水資源保護(hù)及利用規(guī)劃》(以下簡稱“十二五”規(guī)劃)中有關(guān)污水處理廠新建或改造升級、規(guī)?；B(yǎng)殖場糞污治理等措施，《北京市加快污水處理和再生水利用設(shè)施建設(shè)三年行動(dòng)方案(2013—2015)》(以下簡稱“三年行動(dòng)方案”)，人口增長、經(jīng)濟(jì)社會(huì)發(fā)展及農(nóng)業(yè)規(guī)模和布局情況，設(shè)置3套方案，并基于MIKE11模型，評估不同方案對北運(yùn)河流域(北京段)主要污染物的減排效果，以期為流域水污染綜合治理和水環(huán)境管理提供技術(shù)支撐。

1 研究區(qū)概況

北運(yùn)河發(fā)源于燕山南麓，是北京市五大水系中唯一發(fā)源于本市境內(nèi)的水系，在北京市境內(nèi)的流域面積為4 249 km2，約占北京市總面積的26%。北運(yùn)河上游為溫榆河，北運(yùn)河干流(北京段)始于通州區(qū)北關(guān)閘，止于通州區(qū)楊洼閘，沿途接納藺溝、清河、壩河、小中河、通惠河、涼水河等主要支流(見圖1)，北運(yùn)河干流(北京段，含溫榆河干流)全長94 km。北運(yùn)河流域(北京段)聚集了北京市約72%的人口，多年平均水量為8.543億m3。其中，地表徑流量為3.369億m3，地表徑流多被山前水庫攔蓄，經(jīng)社會(huì)水循環(huán)過程后，再以城鎮(zhèn)污水處理廠退水的形式進(jìn)入下游平原河道。北運(yùn)河是北京市主要的排水河道，2010年北京市出境水量為8.29億m3，其中北運(yùn)河為7.51億m3，占90%以上。根據(jù)最新的環(huán)境統(tǒng)計(jì)數(shù)據(jù)測算，北運(yùn)河(北京段)流域2013年的COD、氨氮排放量分別為10.33萬、1.23萬t，分別占北京市總排放的53.3%、57.4%。同時(shí)，平原河段又多受閘壩控制，閘壩阻斷河流水力聯(lián)系，水體流動(dòng)緩慢，產(chǎn)生顯著的水文水環(huán)境效應(yīng)。近10年的監(jiān)測資料表明，北京市境內(nèi)的各水系中，北運(yùn)河流域(北京段)超標(biāo)河流數(shù)量最多，其劣Ⅴ類河流占比超過本流域河流數(shù)量的90%。

圖1 北運(yùn)河流域(北京段)概況Fig.1 The overview of the North Canal River (Beijing section)

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源

本研究以2010年為基準(zhǔn)年，2015年為規(guī)劃水平年。2010年各閘壩日下泄流量、閘壩上下游日水位以及降雨、蒸發(fā)等數(shù)據(jù)均由北京市水文總站提供。各水質(zhì)監(jiān)測斷面數(shù)據(jù)由北京市環(huán)境保護(hù)監(jiān)測中心提供。北京市城區(qū)入河排污口流量、主要污染物濃度及排污類型等信息來自北京市污普信息數(shù)據(jù)庫。構(gòu)建模型所需的北運(yùn)河干流(北京段，含溫榆河干流)及其支流共11個(gè)河流的斷面信息通過實(shí)地監(jiān)測獲取。工業(yè)企業(yè)、農(nóng)業(yè)(種植業(yè)、畜禽養(yǎng)殖業(yè)、水產(chǎn)養(yǎng)殖業(yè))和集中處理設(shè)施的主要污染物排放數(shù)據(jù)來源于《北京市環(huán)境統(tǒng)計(jì)數(shù)據(jù)手冊(2013年)》，生活源直排量由生活污染物的產(chǎn)生量減去污水處理廠的處理量得到，而生活污染物產(chǎn)生量由產(chǎn)污系數(shù)法得到。人口數(shù)據(jù)來源于《北京市人口統(tǒng)計(jì)年鑒(2000—2013年)》，人口增長率根據(jù)2000—2013年北京市城鎮(zhèn)和農(nóng)村人口變化趨勢得到。規(guī)劃水平年污水處理廠升級改造、再生水廠建設(shè)規(guī)劃等信息來自“十二五”規(guī)劃和“三年行動(dòng)方案”。

2.2 MIKE11模型基本原理

分別應(yīng)用MIKE11模型中的水動(dòng)力學(xué)(HD)模塊、平流—擴(kuò)散(AD)模塊構(gòu)建北運(yùn)河干流(北京段)的水動(dòng)力模型和水質(zhì)模型。

MIKE11模型中的HD模塊是基于垂向積分的物質(zhì)守恒方程和動(dòng)量守恒方程，即一維Saint-Venant方程組，其模擬結(jié)果為河道各個(gè)斷面、各個(gè)時(shí)刻的水位和流量等水文要素信息。

(1)

式中：A為過水?dāng)嗝婷娣e，m2；t為時(shí)間，s；Q為流量，m3/s；x為沿河道水流方向的距離，m；q為單寬旁側(cè)入流量，m2/s；g為重力加速度，m/s2；h為水位，m；B為謝才系數(shù)，m1/2/s；R為水力半徑，m。

AD模塊根據(jù)HD模塊的水動(dòng)力條件，應(yīng)用對流擴(kuò)散方程計(jì)算污染物濃度。該模塊的基礎(chǔ)是一維對流—彌散方程，方程的基本假設(shè)為：污染物在斷面上完全混合，遵循一級反應(yīng)動(dòng)力學(xué)，且符合Fick擴(kuò)散定律(擴(kuò)散與濃度梯度呈正比)。一維對流—彌散方程如下：

(2)

式中：c為污染物質(zhì)量濃度，mg/L；D為縱向擴(kuò)散系數(shù)，m2/s；K為線性衰減系數(shù)，s-1；cs為污染物源匯的質(zhì)量濃度，mg/L。

2.3 主控監(jiān)測斷面及水質(zhì)目標(biāo)選擇

主控監(jiān)測斷面選在北運(yùn)河干流(北京段)的榆林莊，該斷面有完整的水位、流量和水質(zhì)監(jiān)測資料，方便模型的水動(dòng)力參數(shù)率定與水質(zhì)驗(yàn)證。同時(shí)，榆林莊斷面也是北京市出境水質(zhì)考核的國家級斷面，具有典型性。

根據(jù)2013年監(jiān)測資料，榆林莊斷面COD、氨氮的年均質(zhì)量濃度分別為49.4、16.2 mg/L，高于《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)中的Ⅴ類標(biāo)準(zhǔn)限值(40、2 mg/L)，尤其是氨氮的削減、控制難度很大，在近期內(nèi)很難實(shí)現(xiàn)總量控制的目標(biāo)。綜合“十二五”規(guī)劃要求和流域?qū)嶋H情況，將“十二五”末(2015年底)榆林莊斷面水質(zhì)控制目標(biāo)設(shè)定為COD低于40 mg/L、氨氮低于5 mg/L。

3 結(jié)果分析

3.1 模型構(gòu)建

基于MIKE11構(gòu)建了北運(yùn)河干流(北京段)的水動(dòng)力和水質(zhì)模型，模型以上游流量和下游水位為外邊界條件，同時(shí)綜合考慮了10座閘壩、12條溝渠及支流匯入和68個(gè)污染源直接匯入等內(nèi)邊界條件，將農(nóng)業(yè)源概化為隨降雨徑流匯入，農(nóng)村生活污染源概化為沿河線源匯入。

3.1.1 水動(dòng)力邊界條件

基于2010年1月1日至12月31日的日流量、水位資料，河道上游采用沙河閘下泄日流量邊界條件，下游采用楊洼閘日水位邊界，各支流采用流量邊界條件。其中，清河采用羊坊閘下泄流量，通惠河采用高碑店湖下泄流量，運(yùn)潮減河采用北關(guān)分洪閘流量，涼水河采用張家灣閘下泄流量，玉帶河采用碧水污水處理廠流量，其他支流根據(jù)實(shí)測流量估算。另外，根據(jù)通州雨量站和楊洼雨量站2010年的降雨數(shù)據(jù)，通州蒸發(fā)站的蒸發(fā)數(shù)據(jù)，以及北關(guān)攔河閘和楊洼閘的下泄流量，利用MIKE11自帶的NAM模塊計(jì)算北運(yùn)河干流(北京段)和溫榆河干流的降雨徑流量。

3.1.2 水質(zhì)邊界條件

北運(yùn)河干流(北京段)上游有監(jiān)測數(shù)據(jù)的支流采用2010年COD、氨氮月平均濃度作為水質(zhì)邊界條件。其中，清河采用沙子營監(jiān)測數(shù)據(jù)，壩河采用沙窩監(jiān)測數(shù)據(jù)，小中河采用北關(guān)小中河橋監(jiān)測數(shù)據(jù)，通惠河采用通惠橋監(jiān)測數(shù)據(jù)，涼水河采用張家灣監(jiān)測數(shù)據(jù)，藺溝采用后藺溝橋監(jiān)測數(shù)據(jù)，其他無監(jiān)測數(shù)據(jù)的支流及排水溝采用實(shí)測的COD、氨氮濃度作為水質(zhì)邊界條件。北運(yùn)河干流(北京段)下游的楊洼閘采用開邊界。根據(jù)污染源的空間分布，按照就近原則通過較近的支流匯入干流，污染源的年均濃度則是根據(jù)2010年污染源動(dòng)態(tài)更新數(shù)據(jù)，由COD、氨氮的年排放總量除以排水量計(jì)算。

3.2 水動(dòng)力、水質(zhì)模擬結(jié)果驗(yàn)證

將模擬時(shí)段設(shè)置為2010年1月1日至12月31日，時(shí)間步長設(shè)為2 min。首先將MIKE11設(shè)置成自動(dòng)修正參數(shù)模式，利用NAM模塊計(jì)算降雨徑流量，然后進(jìn)行水動(dòng)力試算。將最后一天試算的水位狀況作為水位初始條件,采用水動(dòng)力熱啟動(dòng)選項(xiàng),再次計(jì)算,得到榆林莊斷面的水位變化情況,模擬的水位與榆林莊斷面實(shí)測水位基本吻合(見圖2)。在水動(dòng)力計(jì)算的基礎(chǔ)上，模擬水體中COD和氨氮的變化規(guī)律，COD、氨氮初始質(zhì)量濃度設(shè)定為研究區(qū)全部水質(zhì)監(jiān)測斷面2010年12月的均值，分別為54.3、14.5 mg/L。經(jīng)計(jì)算，得到榆林莊斷面COD、氨氮逐日濃度變化曲線，如圖3所示。分別將模擬值與實(shí)測值進(jìn)行對比，結(jié)果表明：起初幾個(gè)月由于模型初始邊界濃度設(shè)置等問題，模型需要一定的預(yù)熱期，導(dǎo)致模擬效果較差；其余月份的模擬值與實(shí)測值較接近，模擬值與實(shí)測值的誤差基本小于15%，且模擬值與實(shí)測值的變化趨勢吻合較好，氨氮、COD的模擬值與實(shí)測值相比，全年誤差分別為-10.5%、-3.0%。總體模擬效果較好，可以進(jìn)行情景模擬分析。

圖2 2010年榆林莊斷面模擬水位與實(shí)測水位對比Fig.2 Comparison of the simulated and observed water levels of Yulinzhuang section in 2010

圖3 2010年榆林莊斷面COD、氨氮實(shí)測值與模擬值對比Fig.3 Comparison of the simulated and observed COD and ammonia nitrogen of Yulinzhuang section in 2010

3.3 減排方案的情景模擬分析

3.3.1 減排方案設(shè)置

(1) 方案1：污水處理廠提標(biāo)改造

假設(shè)農(nóng)業(yè)結(jié)構(gòu)不作調(diào)整，2015年農(nóng)業(yè)源排放規(guī)模維持在2010年水平；工業(yè)源排放占比很小，不作為減排重點(diǎn)，考慮到GDP增長的實(shí)際情況和“十二五”規(guī)劃中北京市工業(yè)發(fā)展定位，工業(yè)源排放維持在2013年水平。新建和升級改造污水處理廠按“十二五”規(guī)劃和“三年行動(dòng)方案”實(shí)施，其出水執(zhí)行《城鎮(zhèn)污水處理廠水污染物排放標(biāo)準(zhǔn)》(DB 11/890—2012)。城鎮(zhèn)人口和農(nóng)村人口以2010年人口為基準(zhǔn)，分別按2.67%、1.76%的增長率預(yù)測2015年人口規(guī)模。其中，城鎮(zhèn)人口產(chǎn)生的COD和氨氮大部分被污水處理廠收集、削減后，通過污水處理廠退水的形式，以相對較低的濃度進(jìn)入河、湖等水體；而農(nóng)村地區(qū)污水處理設(shè)施比較落后，農(nóng)村人口產(chǎn)生的COD、氨氮按直排統(tǒng)計(jì)。

方案1主要通過增加污水處理廠的處理能力和提高污水廠退水口的水質(zhì)標(biāo)準(zhǔn)來減少城鎮(zhèn)生活源的污染物排放。該方案下北運(yùn)河流域(北京段)2015年工業(yè)源氨氮、COD總排放量比2010年有所增加，主要是經(jīng)濟(jì)增長帶來了排放增量。由于污水處理能力及生活污水處理量大大增加，而同時(shí)出水濃度限值標(biāo)準(zhǔn)提高，最終通過集中處理設(shè)施排放的COD高于2010年的污水處理廠排放量，而氨氮排放下降明顯。同時(shí)，污水集中處理率增加使直排量同步減少，總體上生活源排放大幅下降。氨氮、COD總排放量分別減少至1.12萬、11.62萬t，比2010年(氨氮1.96萬t、COD 15.77萬t)分別降低42.9%、26.3%，削減氨氮0.84萬t，COD 4.15萬t。

(2) 方案2：污水處理廠提標(biāo)改造+提高綜合排放標(biāo)準(zhǔn)

在方案1的基礎(chǔ)上，進(jìn)一步按照《水污染物綜合排放標(biāo)準(zhǔn)》(DB 11/307—2013)提高存在污染物直排的工業(yè)、生活源的排放標(biāo)準(zhǔn)。具體而言，污水未進(jìn)入污水處理廠而存在直排的工業(yè)源以2013年排放為基準(zhǔn)，按DB 11/307—2013進(jìn)行污染物削減，其中氨氮、COD的排放限值分別為1.5、30 mg/L。農(nóng)村人口產(chǎn)生的氨氮、COD按照DB 11/307—2013，經(jīng)過生活污水處理站后，分別按照5、40 mg/L的排放限值進(jìn)入環(huán)境，對于城鎮(zhèn)人口在方案1中未能納入城市污水處理管網(wǎng)系統(tǒng)的直排部分，亦按照此標(biāo)準(zhǔn)執(zhí)行。

方案2在方案1強(qiáng)化污水集中處理措施的基礎(chǔ)上，又依據(jù)DB 11/307—2013對直排污水進(jìn)行有效控制。該方案下北運(yùn)河流域(北京段)氨氮、COD總排放量分別減少至0.56萬、7.89萬t，比2010年分別削減71.4%、50.0%。方案2在方案1的基礎(chǔ)上又進(jìn)一步削減氨氮0.56萬t，削減COD 3.73萬t。

(3) 方案3：污水處理廠提標(biāo)改造+提高綜合排放標(biāo)準(zhǔn)+農(nóng)業(yè)源減排

在方案2的基礎(chǔ)上，進(jìn)一步實(shí)施農(nóng)業(yè)源減排措施。從功能區(qū)優(yōu)化布局的角度，對于《北京城市總體規(guī)劃(2004—2020年)》提出的中心城區(qū)(壩河流域、清河流域、通惠河流域及涼水河流域的部分地區(qū))農(nóng)業(yè)不保留，中心城區(qū)以外區(qū)域種植業(yè)、水產(chǎn)養(yǎng)殖業(yè)保持2013年的規(guī)模不變；六環(huán)路以內(nèi)為禁養(yǎng)區(qū)，六環(huán)路以外規(guī)?；B(yǎng)殖場保持2013年規(guī)模不變，對養(yǎng)殖專業(yè)戶實(shí)行限養(yǎng)措施，不得養(yǎng)殖污染較嚴(yán)重且本地供應(yīng)保障要求不高的生豬和肉牛。

方案3在方案2的基礎(chǔ)上，依據(jù)北京市畜禽養(yǎng)殖業(yè)污染防治規(guī)劃和未來農(nóng)業(yè)結(jié)構(gòu)、布局方面的調(diào)整思路，在北運(yùn)河流域(北京段)分區(qū)、分類實(shí)施農(nóng)業(yè)源減排。在方案3情景下，氨氮、COD總排放量分別降低至0.39萬、5.22萬t，比2010年分別削減80.1%、66.9%。方案3在方案2的基礎(chǔ)上，氨氮、COD分別進(jìn)一步削減0.17萬、2.67萬t。

不同減排措施下，北運(yùn)河流域(北京段)的氨氮、COD相對2010年的削減量及削減比例預(yù)測如表1所示。從表1可以看出，污水處理廠提標(biāo)改造對于氨氮、COD的削減效果最顯著，可削減26.3%的COD和42.8%的氨氮；其次是提高綜合排放標(biāo)準(zhǔn)，能削減23.7%的COD和28.6%的氨氮；再次是農(nóng)業(yè)源減排，可削減16.9%的COD和8.7%的氨氮。

3.3.2 不同減排方案模擬結(jié)果對比分析

根據(jù)近10年以來北京市年降水、徑流資料的排頻計(jì)算結(jié)果，2010年是枯水年之一，可作為水文條件進(jìn)行水質(zhì)模擬分析。保持2010年流量條件不變，基于3種減排方案，模擬預(yù)測榆林莊斷面河流水質(zhì)隨時(shí)間的變化，結(jié)果見圖4。

在方案1下，榆林莊斷面2015年COD、氨氮的年均值分別為32.3、6.6 mg/L，與2010年的年均值(44.4、16.6 mg/L)相比，分別下降27.2%、60.6%。其中，COD達(dá)到水質(zhì)控制目標(biāo)(低于40 mg/L)，且全年達(dá)標(biāo)天數(shù)為299 d，超標(biāo)時(shí)段主要集中在前汛期6、12月。可見，方案1對于COD減排貢獻(xiàn)很大，但是氨氮離水質(zhì)控制目標(biāo)(低于5 mg/L)還有一定差距，全年達(dá)標(biāo)天數(shù)僅有67 d，在6月5日達(dá)到最大質(zhì)量濃度(11.2 mg/L)。因此，需要在方案1的基礎(chǔ)上配套其他減排措施。

表1 不同減排措施對北運(yùn)河流域(北京段)COD和氨氮的減排貢獻(xiàn)Table 1 Forecast of the contribution of COD and ammonia nitrogen emission reduction in North Canal River (Beijing section) under different reduction measures

圖4 不同減排方案下榆林莊斷面COD、氨氮的變化Fig.4 COD and ammonia nitrogen changes of Yulinzhuang section under different emission reduction programs

在方案2下，氨氮年均值降至4.4 mg/L，與2010年相比下降73.5%，達(dá)到水質(zhì)控制目標(biāo)(低于5 mg/L)，全年達(dá)標(biāo)天數(shù)為267 d。COD與氨氮協(xié)同減排效應(yīng)明顯，COD年均值降低至25.4 mg/L，全年達(dá)標(biāo)天數(shù)增至347 d，超標(biāo)時(shí)段主要集中在前汛期6月上、中旬。可見，加強(qiáng)對污染物直排的控制對氨氮的減排效果相當(dāng)顯著。

在方案3下，COD年均值降低至22.6 mg/L，全年達(dá)標(biāo)天數(shù)為357 d。氨氮年均值降低至4.0 mg/L，全年達(dá)標(biāo)天數(shù)為288 d?？梢?，農(nóng)業(yè)源減排能進(jìn)一步削減COD和氨氮。

3.3.3 減排策略展望

僅僅實(shí)施污水處理廠提標(biāo)改造這一措施難以實(shí)現(xiàn)榆林莊斷面COD、氨氮同時(shí)達(dá)到水質(zhì)控制目標(biāo)。同時(shí)實(shí)施污水處理廠提標(biāo)改造和提高綜合排放標(biāo)準(zhǔn)兩項(xiàng)措施理論上可以實(shí)現(xiàn)COD和氨氮同時(shí)達(dá)標(biāo)，但需要各項(xiàng)減排措施按設(shè)定目標(biāo)落實(shí)到位。目前，城鎮(zhèn)污水處理廠、再生水廠等集中處理設(shè)施已按照北京市“十二五”規(guī)劃和“三年行動(dòng)方案”部署逐步完成提標(biāo)改造，而農(nóng)村生活源等分散面源污染的治理相對滯后，導(dǎo)致DB 11/307—2013的執(zhí)行受阻。模型本身存在一定的不確定性，實(shí)際減排效果可能達(dá)不到預(yù)期目標(biāo)。因此，需要同步實(shí)施農(nóng)業(yè)源減排措施，尤其要加強(qiáng)對畜禽養(yǎng)殖業(yè)污染的集中治理。

需要指出的是，同時(shí)實(shí)施污水處理廠提標(biāo)改造、提高綜合排放標(biāo)準(zhǔn)以及農(nóng)業(yè)源減排3項(xiàng)措施，前汛期6、7月的氨氮濃度仍然可能超標(biāo)，因此上述措施的實(shí)施并不能有效降低前汛期的污染物峰值濃度。這主要是因?yàn)楸边\(yùn)河作為北京市主要排水河道，清河、壩河、通惠河、涼水河等較大支流均流經(jīng)城區(qū)，前汛期的暴雨徑流對河流水質(zhì)影響較大，但減排措施側(cè)重于農(nóng)業(yè)源、生活源、工業(yè)源等，而對城市暴雨徑流污染的影響考慮不夠。因此，還應(yīng)該結(jié)合雨水凈化、雨污分流、雨洪資源利用等措施降低北京市面源污染物的入河量。

4 結(jié) 論

(1) 通過污水處理廠提標(biāo)改造，相比2010年，可削減26.3%的COD和42.8%的氨氮，并實(shí)現(xiàn)榆林莊斷面COD和氨氮分別下降27.2%、60.6%，COD達(dá)到水質(zhì)控制目標(biāo)。

(2) 相對COD而言，氨氮減排難度較大。同時(shí)實(shí)施污水處理廠提標(biāo)改造和提高綜合排放標(biāo)準(zhǔn)，榆林莊斷面的氨氮相比2010年下降73.5%，能實(shí)現(xiàn)氨氮達(dá)標(biāo)。

(3) 城市面源污染對前汛期河流水質(zhì)影響較大，在實(shí)施污水處理廠提標(biāo)改造、提高綜合排放標(biāo)準(zhǔn)和農(nóng)業(yè)源減排以外，還應(yīng)加強(qiáng)對城市暴雨徑流污染的管理，以有效降低前汛期水體污染物峰值濃度。

[1] 王亞煒,杜向群,郁達(dá)偉,等.溫榆河氨氮污染控制措施的效果模擬[J].環(huán)境科學(xué)學(xué)報(bào),2013,33(2):479-486.

[2] 郁達(dá)偉,于淼,魏源送,等.1980—2010年溫榆河的水環(huán)境質(zhì)量時(shí)空演變特征[J].環(huán)境科學(xué)學(xué)報(bào),2012,32(11):2803-2813.

[3] 李延,單保慶,陳靜,等.北運(yùn)河流域典型河流水體耗氧污染特征分析[J].安全與環(huán)境學(xué)報(bào),2014,14(5):311-315.

[4] 曾平,趙云,魯欣,等.北京市區(qū)人工湖水質(zhì)時(shí)空變化特征及水污染防治對策探討[J].環(huán)境污染與防治,2014,36(6):111-115.

[5] 范清,孫長虹.北京市2000年—2012年水環(huán)境質(zhì)量及污染規(guī)律研究[J].環(huán)境科學(xué)與管理,2014,39(7):54-57.

[6] 孫長虹,范清,王永剛,等.北京市城市發(fā)展與水環(huán)境演化規(guī)律研究[J].環(huán)境科技,2014,27(4):31-34.

[7] 孫娟,吳悅穎,文宇立,等.北京市“十二五”期間主要水污染物減排對策研究[J].環(huán)境污染與防治,2011,33(1):102-104.

[8] 張文靜,陳巖,劉雅玲,等.綜合考慮減排效益的污染物總量分配方法研究[J].環(huán)境污染與防治,2015,37(3):107-110.

[9] 白輝,趙翠平,陳巖,等.承德市武烈河水環(huán)境容量核算與減排效益評估研究[J].環(huán)境污染與防治,2015,37(3):31-34.

[10] 左其亭,李冬鋒.基于模擬-優(yōu)化的重污染河流閘壩群防污調(diào)控研究[J].水利學(xué)報(bào),2013,44(8):979-986.

[11] 朱茂森.基于MIKE11的遼河流域一維水質(zhì)模型[J].水資源保護(hù),2013,29(3):6-9.

[12] 李梓嘉,董增川,樊孔明,等.MIKE11模型在泗洪縣城城區(qū)河網(wǎng)引水沖污工程中的應(yīng)用[J].水電能源科學(xué),2012,30(8):100-103.

[13] TANG Caihong,YI Yujun,YANG Zhifeng,et al. Water pollution risk simulation and prediction in the main canal of the South-to-North Water Transfer Project[J].Journal of Hydrology,2014,519:2111-2120.

[14] 卿曉霞,張會(huì)波,周健,等.伏牛溪水污染治理效果的數(shù)值模擬研究[J].環(huán)境工程學(xué)報(bào),2015,9(1):65-72.