DOI：10.19981/j.CN23-1581/G3.2023.32.028

摘" 要：北方冬春兩季氣候干燥，降水不具備形成徑流條件，農(nóng)田殘留化肥、腐敗有機物等污染物在溝渠內(nèi)聚集。汛期溝渠內(nèi)污染物隨地表徑流匯入河道，對干流河道水質(zhì)造成影響。為研究初汛期污染物匯入對河道水質(zhì)的影響，選擇北方淮河流域重點河流南涑河作為研究對象，針對COD和氨氮，建立MIKE11模型模擬在初汛期隨污染物匯入后南涑河水質(zhì)的變化規(guī)律，探索最佳調(diào)度方案解決初汛期污染物含量過高的問題。研究表明，初汛期南涑河省控斷面老屯橋不滿足水質(zhì)要求。通過模型模擬，截污能夠保證初汛期南涑河省控斷面水質(zhì)滿足水質(zhì)目標。該方法為MIKE11模型的廣泛應用提供基礎，并為北方河道初汛期水質(zhì)污染提供解決思路。

關鍵詞：MIKE11；初汛期；聯(lián)合調(diào)度；水質(zhì)污染；水質(zhì)變化規(guī)律

中圖分類號：X522" " " " 文獻標志碼：A" " " " " 文章編號：2095-2945（2023）32-0112-04

Abstract： The climate in winter and spring in northern China is dry， the precipitation does not have the conditions to form runoff， and the pollutants such as residual chemical fertilizer and rotten organic matter gather in the ditch. During the flood season， the pollutants in the ditch converge into the river with the surface runoff， which has an impact on the water quality of the main river. In order to study the influence of pollutant influx on river water quality in the early flood season， Nansu River， a key river in the Huaihe River Basin in the north， was selected as the research object. Aiming at COD and ammonia nitrogen， a MIKE11 model was established to simulate the change law of Nansu River water quality with the influx of pollutants in the early flood season， and to explore the best dispatching scheme to solve the problem of excessive pollutant concentration in the early flood season. The study shows that Laotun Bridge in the control section of Nansu River does not meet the requirements of water quality in the early flood season. Through the model simulation， sewage interception can ensure that the water quality of the control section of Nansu River in the early flood season can meet the water quality target. This method provides a basis for the wide application of MIKE11 model， and provides a solution to water pollution in the early flood season of northern rivers.

Keywords： MIKE11; early flood season; joint dispatching; water pollution; water quality change law

南涑河屬于淮河流域，種植業(yè)、畜禽養(yǎng)殖業(yè)發(fā)達，但相應污水處理設施建設不完善，農(nóng)村生活污水集中處理率低，會對周圍環(huán)境產(chǎn)生不利影響。北方冬春兩季氣候干燥，降水不具備形成徑流條件，殘留化肥、腐敗有機物等污染物容易在農(nóng)田溝渠內(nèi)聚集。在汛期，農(nóng)田溝渠內(nèi)的污染物隨地表徑流匯入河道，造成河道水質(zhì)惡化[1]。此外，流域內(nèi)城市老舊小區(qū)、城郊、城中村等生活污水收集管網(wǎng)不完善，未被收集的生活污水也會隨降雨進入河道，引起河道水質(zhì)惡化。

目前世界上針對河道水質(zhì)水量的研究相對較多，丹麥水利單位（DHI）研究開發(fā)的MIKE系列軟件是目前驗證最多的水文權威軟件[2]。MIKE模型廣泛應用于汛期對行洪安全的研究，張林[3]利用MIKE11模型對章水流域南康區(qū)河道典型斷面的防洪能力進行了評價。蘇曉天[4]以MIKE FLOOD為平臺，研究不同設計降雨情景下區(qū)域的水文區(qū)域和內(nèi)澇風險。在水環(huán)境方面，田凱達等[5]利用MIKE11構建合肥市十五里河水動力及水質(zhì)模型成功研究該河道水質(zhì)改善效果。張斯思[6]利用MIKE11研究渦河污染物變化規(guī)律，并對點源截污、污水廠升級改造等方案進行了預測，為河流綜合治理工程決策提供借鑒和依據(jù)。董增川等[7]通過搭建水質(zhì)水量調(diào)度模型研究水工建筑物調(diào)度對水質(zhì)的改善效果。

本文以淮河流域典型河道南涑河作為研究對象，利用MIKE11模型研究初汛期南涑河水質(zhì)變化規(guī)律，同時采用水質(zhì)水量聯(lián)合調(diào)度的方式，解決初汛期污染的問題。這對于研究淮河流域初汛期河道水質(zhì)變化規(guī)律及水質(zhì)水量聯(lián)合調(diào)度具有較大研究意義。

1" 區(qū)域概況及研究思路

南涑河位于淮河流域邳蒼郯新片區(qū)，長42 km，流域面積279.1 km2，河道比降約1/1 000～1/2 000。2020年對南涑河進行治理，河道治理長度21.55 km，防洪標準為20年一遇。南涑河流域上游流經(jīng)臨沂工業(yè)園區(qū)，中游流經(jīng)羅莊區(qū)城區(qū)，下游流經(jīng)大面積農(nóng)田，區(qū)域第一產(chǎn)業(yè)、第二產(chǎn)業(yè)比重大，村莊、社區(qū)人口分布較集中，上游小型工業(yè)企業(yè)分布較為密集，下游農(nóng)田種植業(yè)、畜禽養(yǎng)殖業(yè)等產(chǎn)業(yè)較為發(fā)達。南涑河入河污染物主要來自沿河城鎮(zhèn)生活污水、農(nóng)村生活污水、污水處理廠尾水、農(nóng)業(yè)種植面源污染及畜禽養(yǎng)殖廢水。

南涑河流域共有3處污水處理廠，出水均執(zhí)行GB 18918—2002《城鎮(zhèn)污水處理廠污染物排放標準》一級A標準；南涑河流域小型工業(yè)企業(yè)、鄉(xiāng)鎮(zhèn)企業(yè)分布較廣，主要產(chǎn)業(yè)類型包括化工、鑄造、食品、建材和制藥等，工業(yè)廢水主要經(jīng)由污水處理廠處理達標后進入水體；南涑河流域穿越城區(qū)段生活污水基本都納入市政管網(wǎng)經(jīng)污水處理廠處理后排放，管網(wǎng)覆蓋率超過90%，根據(jù)各縣區(qū)統(tǒng)計年鑒估算，南涑河流域城鎮(zhèn)人口按9.18萬計，污水集中處理率按95%計算，入河系數(shù)按0.2計算，南涑河未集中處理城市生活污水入河排污量1.809萬t/a，COD總?cè)牒恿考s14.640 t/a，氨氮總?cè)牒恿考s2.009 t/a；根據(jù)現(xiàn)場調(diào)查，流域約50%村莊尚未建設排水管道，農(nóng)村生活污水排水系數(shù)為0.33，則農(nóng)村生活污水排放量為16.5 L/（人·日）。按照南涑河流域農(nóng)村人口及排污水質(zhì)計算農(nóng)村生活污水產(chǎn)污量，農(nóng)村生活污水集中處理率按35%計，根據(jù)村莊及鄉(xiāng)鎮(zhèn)駐地距離河道遠近，入河系數(shù)按平均0.1取值，入河排放量7.184萬t/a，COD總?cè)牒恿考s14.368 t/a，氨氮總?cè)牒恿考s1.437 t/a。根據(jù)《第一次全國污染源普查-農(nóng)業(yè)污染源 肥料流失系數(shù)手冊》及南涑河流域范圍內(nèi)農(nóng)田種植作物類型、種植面積，根據(jù)肥料流失系數(shù)確定南涑河流域農(nóng)田徑流污染源強為COD 1.797 kg/（畝·年）（1畝約等于667 m2，下同），氨氮0.104 kg/（畝·年），入河量按20%計，流域耕地面積約37萬畝，經(jīng)核算南涑河農(nóng)田種植業(yè)COD總?cè)牒恿考s132.59 t/a，氨氮總?cè)牒恿考s7.67 t/a。南涑河流域主要養(yǎng)殖種類為雞、鴨、鵝、豬和羊，流域內(nèi)規(guī)?；笄蒺B(yǎng)殖場污水處理設施和糞污資源化利用設施配建率基本達到100%，本次污染源核算主要考慮規(guī)模以下畜禽養(yǎng)殖場產(chǎn)生的糞污排放情況。畜禽養(yǎng)殖產(chǎn)污及排污系數(shù)采用《第一次全國污染源普查畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊》相關系數(shù)，本次研究范圍屬于華東區(qū)，采用華東區(qū)相關系數(shù)，COD入河排污量約為10.57·t/a。

南涑河現(xiàn)設有老屯橋斷面1個地表水環(huán)境省控考核斷面，如圖1所示，位于羅莊區(qū)老屯村大橋處，通過分析南涑河2019—2021年歷年水質(zhì)數(shù)據(jù)可知，南涑河流域水質(zhì)滿足地表水Ⅳ類水水質(zhì)要求。但由于南涑河流域污染源較多，初汛期溝渠內(nèi)累積冬春兩季的污染物隨地表徑流匯入南涑河，污染難以得到有效控制等原因，省控斷面在初汛期間仍有不滿足近期水質(zhì)目標Ⅴ類水要求的可能。

2" 流域模型構建

2.1" NAM水文模型構建

本文利用NAM水文模型耦合計算多個旁側(cè)入流模擬流域內(nèi)產(chǎn)匯流過程，基于90 m×90 m的DEM數(shù)據(jù)和南涑河集水區(qū)域，將研究區(qū)域劃分為不同的水文分區(qū)，每個水文分區(qū)的參數(shù)均為本區(qū)域的參數(shù)平均值。根據(jù)收集到的流域內(nèi)羅西、羅莊、角沂和臨沂等8個水文站點的2020年降雨數(shù)據(jù)，利用泰森多邊形法計算南涑河水文分區(qū)內(nèi)平均降雨量，流域蒸發(fā)數(shù)據(jù)采用2020年老屯站蒸發(fā)數(shù)據(jù)。

將每個水文分區(qū)概化成一個NAM模型，水文模型的產(chǎn)匯流參數(shù)是模型的率定參數(shù)，其余均為確定性參數(shù)。

2.2" 水動力學模型

2.2.1" MIKE11 HD簡介

MIKE11 HD水動力計算模型基于一維非恒定流Saint-Venant方程組來模擬河流或河口的水流狀態(tài)[8]，利用Abbott-Ionescu六點隱式有限差分格式求解[9]。

2.2.2" 水動力學模型構建

1）河網(wǎng)概化。概化河網(wǎng)綜合考慮了南涑河及其主要支流包括老龍溝、漁梁溝、降水溝和藕蒲溝等，最終確定概化15條河道，概化河道總長度196.773 km。

2）斷面概化。南涑河、小黑溝、降水溝、藕蒲溝和老涑河采用實測斷面數(shù)據(jù)，每隔100 m設置1個斷面；老龍溝、洪溝、漁梁溝等其他支流利用地形數(shù)據(jù)、相應河道的河寬水深數(shù)據(jù)及河道坡降數(shù)據(jù)，以梯形斷面的方式進行概化。

3）邊界設置。HD模塊將NAM模型得到的徑流結(jié)果通過設置耦合鏈接匯入相應支流河道[10]。模型的外部邊界為上游入流和下游出口的邊界條件，河道上游端點給定一個恒定的流量值作為源頭入流邊界條件，區(qū)間子流域入流同樣也通過NAM水文模型計算，并連接到河道中，模型下游采用實測水位進行邊界設置。

4）水工結(jié)構物設置。MIKE11 SO將水工建筑物分為不涉及調(diào)度規(guī)則的一般水工建筑物和涉及調(diào)度規(guī)則的水工建筑物[11]。通過設置不同的調(diào)度規(guī)則控制水工建筑物的運行方式，實現(xiàn)水質(zhì)水量聯(lián)合調(diào)度。

2.3" 水質(zhì)模型

MIKE11 AD模塊能夠?qū)λw中的可溶性物質(zhì)和懸浮性物質(zhì)對流擴散過程進行模擬，其根據(jù)HD模塊產(chǎn)生的水動力條件，應用對流擴散方程進行計算并通過設定衰減常數(shù)來模擬非保守物質(zhì)[12]。

2.4" 模型率定

2.4.1" 水動力模型率定

河床的糙率系數(shù)是水動力學模型的重要參數(shù)，也是模型率定的重點。根據(jù)研究經(jīng)驗，模型河道糙率初始值設置為0.03 s/m1/3，通過人工手動率定的方法，對河道糙率系數(shù)進行了率定，確定最終取值。模擬過程中，納什系數(shù)為0.84，模擬效果良好。

2.4.2" 水質(zhì)模型率定

選擇南涑河老屯橋斷面為驗證斷面進行水質(zhì)率定，COD誤差均值為10%，氨氮誤差均值為2%，模擬效果良好。

3" 情景分析

3.1" 情景概述

冬春兩季氣候干燥，污染物在溝渠內(nèi)聚集，在初汛期隨降雨徑流匯入河道，造成河道水質(zhì)惡化，為研究初汛期河道水質(zhì)污染現(xiàn)象，選擇2020年6月16日—18日作為研究時段，利用MIKE模型模擬污染物在初汛期的匯入過程。選擇COD、氨氮2項指標反映初汛期水質(zhì)情況。水質(zhì)模擬結(jié)果顯示隨著污染物不斷匯入，河道COD、氨氮質(zhì)量濃度逐漸升高，無法滿足省控斷面水質(zhì)要求。

3.2" 調(diào)度方案

3.2.1" 截污

根據(jù)污染物匯入情況，超標河段主要集中在老龍溝、泄洪渠（小清河）和藕蒲溝（降水溝），其中泄洪渠和藕蒲溝上存在支流閘，出于研究考慮老龍溝上假定存在支流閘，水閘分布如圖 1所示。在滿足行洪安全的前提下，通過控制支流節(jié)制閘，減少污染物匯入南涑河，實現(xiàn)治理截污改善干流水質(zhì)的目的。

3.2.2" 調(diào)水

南涑河上游為涑河，通過調(diào)水稀釋，能夠改善下游水質(zhì)。根據(jù)前期研究成果，調(diào)水量25 m3/s時，水質(zhì)依然不滿足省控斷面水質(zhì)要求，根據(jù)2020年6月18日前涑河（小涑河）周井輔站和鼓樓臺站實測流量值，本次只考慮調(diào)水量5 m3/s的情況。

采用水質(zhì)水量聯(lián)合調(diào)度的方式解決南涑河初汛期水質(zhì)污染的問題，模擬情景為初汛期水質(zhì)變化（支流閘開閘無調(diào)水），調(diào)度方案見表1。

3.3 結(jié)果分析

3.3.1" 模擬結(jié)果

調(diào)度結(jié)果如圖2所示。由圖1可知，調(diào)度方案一、調(diào)度方案五、調(diào)度方案六和調(diào)度方案七的COD與氨氮均滿足省控斷面水質(zhì)要求。調(diào)度方案二、調(diào)度方案三、調(diào)度方案四不滿足水質(zhì)要求，說明泄洪渠、老龍溝支流截污及調(diào)水稀釋效果不明顯，這是因為泄洪渠、老龍溝主要流經(jīng)城區(qū)，不存在面源污染，城市生活污水匯入有限，因此截污效果不好。

3.3.2" 方案分析

由方案一、二、三可知，泄洪渠和老龍溝水質(zhì)污染較輕，對南涑河水質(zhì)影響較小，在關閉閘門的條件下，截污效果不明顯，且由于上游來水量的減少，反而不利于下游水質(zhì)的凈化。藕蒲溝對省控斷面水質(zhì)影響最大，藕蒲溝是南涑河上最大的支流，流域內(nèi)以農(nóng)田為主，大量殘留化肥、農(nóng)藥在溝渠內(nèi)聚集，在初汛期隨地表徑流匯入河道造成污染。在地理位置上，藕蒲溝閘距離老屯橋斷面不足4 km，污染物得不到有效降解稀釋，這也是藕蒲溝對干流水質(zhì)影響最大的一個原因。

研究結(jié)果表明，截污是改善河道水質(zhì)的重要途徑，通過控制節(jié)制閘，減少污染物匯入干流，能明顯提升干流水質(zhì)。截污受支流污染物含量和流量的影響，在本次研究中，藕蒲溝是南涑河主要支流，流域面積大，對干流水質(zhì)影響大，在控制藕蒲溝閘的前提下，能夠有效改善南涑河水質(zhì)。此外，截污需要考慮行洪安全，本次研究藕蒲溝滿足五年一遇除澇設計標準。

汛前開閘泄水不但能夠緩解河道行洪壓力，同時對下游水質(zhì)起到稀釋作用。由方案四可知，自涑河調(diào)水至南涑河能夠令COD含量達到省控斷面水質(zhì)要求，但氨氮含量仍然超標。這是因為涑河本身含有污染物，且調(diào)水量受實際水量的限制，導致稀釋效果不明顯。因此，受水質(zhì)、水量的影響，調(diào)水對水質(zhì)的改善效果有限。

研究結(jié)果表明，截污和調(diào)水是改善河道水質(zhì)的重要途徑，在本次研究中，截污效果更明顯，在保證行洪安全的前提下，關閉藕蒲溝支流閘，能夠減少污染匯入，改善干流水質(zhì)。

4" 結(jié)束語

本文基于MIKE11模型搭建南涑河流域一維水文水動力水質(zhì)耦合模型，可應用于模擬污染物匯入后南涑河水質(zhì)變化規(guī)律。研究發(fā)現(xiàn)，在初汛期污染物隨地表徑流匯入河道，會造成南涑河水質(zhì)超標，為解決這一問題，通過施行截污和調(diào)水2種方式，在保證行洪安全的前提下，能夠?qū)λ|(zhì)起到改善效果，對污染最嚴重的支流實施截污是改善干流河道水質(zhì)的最優(yōu)方式。

參考文獻：

[1] 李海燕，黃宇，車伍，等.非點源污染對北京城區(qū)典型流域水質(zhì)影響的研究[J].給水排水，2007（S1）：60-63.

[2] 姚力瑋.基于MIKE11的嫩江干流水環(huán)境容量模型改進研究[D].北京：華北電力大學（北京），2017.

[3] 張林.基于MIKE11模型的河道防洪能力評價研究[J].黑龍江水利科技，2019，47（12）：70-72.

[4] 蘇曉天.基于MIKE FLOOD的華北典型縣域主城區(qū)內(nèi)澇模擬研究[D].邯鄲：河北工程大學，2023.

[5] 田凱達，劉曉薇，王慧，等.MIKE11模型在合肥市十五里河水質(zhì)改善研究中的應用[J].水文，2019，39（4）：18-23.

[6] 張斯思.基于MIKE11水質(zhì)模型的水環(huán)境容量計算研究——以渦河為例[D].合肥：合肥工業(yè)大學，2023.

[7] 董增川，卞戈亞，王船海，等.基于數(shù)值模擬的區(qū)域水量水質(zhì)聯(lián)合調(diào)度研究[J].水科學進展，2009，20（2）：184-189.

[8] 李昆朋，付東王.基于水動力模型的平原河網(wǎng)引補水研究[J].吉林水利，2019（10）：9-12.

[9] 張海軍，馬永法，葉萍，等.基于數(shù)值模擬的嘉興市區(qū)河道水質(zhì)優(yōu)化調(diào)度方案[J].凈水技術，2020，39（5）：54-68.

[10] 周旭，黃莉，王蘇勝.MIKE11模型在南通平原河網(wǎng)模擬中的應用[J].江蘇水利，2016（1）：52-55，60.

[11] 黃慧慧，張舒，潘靜也，等.水文模型在圩區(qū)水環(huán)境容量計算中的應用研究[J].中國水運（下半月），2018，18（7）：172-174.

[12] 劉開印.洸府河下游中水阻隔消納方案體系研究[D].濟南：山東大學，2023.

基金項目：山東省淮河流域重點平原洼地沿運片邳蒼郯新片區(qū)治理工程科學研究項目（PCTXPQ-KY202002）

第一作者簡介：董坤明（1987-），男，碩士，工程師。研究方向為水利工程規(guī)劃設計。

*通信作者：韓慶港（1997-），男，助理工程師。研究方向為環(huán)境科學等。