王瑩 張幼寬 梁修雨 謝顯傳

摘要：基于對沙潁河安徽段污染負荷的調(diào)查以及水質(zhì)和流量監(jiān)測數(shù)據(jù)的分析，采用MIKE11模型建立水流與水質(zhì)模型，采用線性規(guī)劃法，以入淮河干流Ⅲ類水為水質(zhì)目標，計算沙潁河安徽段水環(huán)境容量并分析其動態(tài)特征。結(jié)果表明：沙潁河安徽段CODM。和氨氮年水環(huán)境容量分別為2.39萬、0.76萬t/a：受流量及上游水質(zhì)影響，COD和氨氮水環(huán)境容量月變化較大，8月水環(huán)境容量最大，分別為3 557、1363 t，3月水環(huán)境容量最小，分別為621、181 t;上游水質(zhì)對沙潁河安徽段水環(huán)境容量影響顯著。

關(guān)鍵詞：水環(huán)境容量;水質(zhì)模型;線性規(guī)劃法;沙潁河

中圖分類號：X522

文獻標志碼：A

doi：10. 3969/j .issn. 1000- 13 79.2019.01. 018

水環(huán)境容量是指水體環(huán)境在規(guī)定的環(huán)境目標下所能容納的污染負荷[1]。水環(huán)境容量與水體特征、水質(zhì)目標及污染物特性有關(guān)，同時還與污染物的排放方式及排放的時空分布密切相關(guān)[2]，是容量總量控制的核心內(nèi)容之一[3]。水環(huán)境容量計算的基本思路：通過解析或數(shù)值等方法計算河流的水流過程及污染物遷移轉(zhuǎn)化規(guī)律，建立污染負荷排放與水質(zhì)響應(yīng)的關(guān)系，基于給定的環(huán)境目標（如水質(zhì)目標）計算指定區(qū)域所能容納的最大污染負荷。對于特定的河流，在水質(zhì)目標不變的情況下，由于水文特征、自凈能力和排污口排污方式的時空變異性，因此水環(huán)境容量是隨時間和空間變化的。流域?qū)嶋H的水環(huán)境容量需要針對實際背景條件計算，即計算動態(tài)水環(huán)境容量[4]。水環(huán)境容量的計算方法主要包括公式法、模型試錯法、系統(tǒng)最優(yōu)化法（如線性規(guī)劃法和隨機規(guī)劃法）、概率稀釋模型法和模糊數(shù)學(xué)法[5-6]，其中系統(tǒng)最優(yōu)化法適用范圍廣，對邊界條件和設(shè)計流量的適應(yīng)能力較強，是目前水環(huán)境容量計算的主流方法之一[7]。

本研究以沙潁河流域安徽段為例，采用系統(tǒng)最優(yōu)化法中的線性規(guī)劃法計算水環(huán)境容量。整體思路：根據(jù)排污數(shù)據(jù)和水功能區(qū)劃確定排污口和控制斷面的位置：通過河流水流和水質(zhì)數(shù)值模擬，建立排污口與控制斷面的水質(zhì)響應(yīng)關(guān)系，計算實際本底濃度值和各排污口負荷貢獻度：建立線性最優(yōu)化模型，計算流域水環(huán)境容量并分析其動態(tài)特征。

1 研究區(qū)概況

沙潁河是淮河最大的一級支流，全長624 km，流域面積39 980 km2。沙潁河流域年平均氣溫約15℃，年平均降水量約700 mm/a。沙潁河流域跨越河南和安徽兩省，在安徽境內(nèi)長208 km，流域面積4% km2，占沙潁河流域總面積的10%。沙潁河在安徽境內(nèi)被稱為潁河，其主要支流為汾泉河（見圖1）。雖然沙潁河流域僅占淮河流域總面積的1/7，但其污染負荷排放量占淮河流域的1/3。近年來，隨著工業(yè)發(fā)展和人口增長，沙潁河水質(zhì)日趨惡化，污染嚴重，水質(zhì)劣于V類的水功能區(qū)占總數(shù)的60%[8]。為控制流域水污染，實現(xiàn)沙潁河人淮河水質(zhì)達到Ⅲ類水的目標，需要制定科學(xué)合理的污染控制方案，而河流的納污能力（水環(huán)境容量）是制定該方案的基礎(chǔ)。

沙潁河安徽段的主要污染物來自城鎮(zhèn)生活廢水排放和工業(yè)點源污染。根據(jù)2012年安徽省環(huán)境統(tǒng)計數(shù)據(jù)，該段共有排污口約166個，其空間分布及排污負荷見圖2。

為計算水環(huán)境容量，采用重心概化法將干流沿岸的排污口概化為6個排污口。

2 數(shù)值模擬

2.1 水流模擬

河流水流過程是水質(zhì)模擬和水環(huán)境容量計算的基礎(chǔ)，采用數(shù)值方法模擬水流過程。研究僅考慮污染物沿河流流向的一維遷移過程[9]，因此選用一維水動力模型MIKE11的HD模塊模擬河流徑流量。MIKE11水動力計算模型是基于一維非恒定流Saint -Venant方程組來模擬河流或河口的水流狀態(tài).Saint - Venant方程組為

對式（1）采用Abbott -Ionescu六點隱式有限差分格式求解。剖分網(wǎng)格尺寸最大值為5 000 m，并在局部加密，共生成179個計算節(jié)點。時間步長為5 min，總計算時段為2 a（2011-01-01-2012-12-31）。以界首閘和沈丘閘實測日流量為上游流量邊界，以潁上閘實測日河水位為下游水位邊界，初始河水位通過上下游觀測水位線性插值獲得。采用阜陽閘和潁上閘的實測日流量與水位進行模型率定驗證，率定期為2011年，驗證期為2012年。根據(jù)MIKE11模型按過水?dāng)嗝嬖O(shè)置糙率的方法，將斷面分為上、中、下三部分，率定參數(shù)為斷面上、中、下三部分的曼寧糙率n。研究區(qū)河道斷面變化較大，為有效校正模型[10].斷面三部分的初始n值均為0.02[11]?；谟^測流量和水位，通過調(diào)整參數(shù)，獲得最佳擬合參數(shù)。

2012年阜陽閘和潁上閘流量與水位的模擬值與觀測值見圖3、圖4.其中：2012年阜陽閘流量的相對誤差為11.4%（Nash系數(shù)為0.83），水位的相對誤差為5.3%（Nash系數(shù)為0.93）;潁上閘流量的相對誤差為12.3%（Nash系數(shù)為0.75），水位的相對誤差為0.3%（Nash系數(shù)為0.97）。模擬結(jié)果基本符合實際水文情況。通過驗證，獲得斷面上、中、下三部分的n值分別為0.020，0.022和0.025。驗證的水流模型可用于河流COD和氨氮運移模擬以及水環(huán)境容量計算。

2.2 水質(zhì)模擬

污染物在河流中遷移的控制方程為一維對流彌散方程，考慮一級降解，其方程為

式（3）主要參數(shù)是彌敢度a和降解系數(shù)k，筆者參考沙潁河流域已有研究成果[12]，選取合適的彌散度a，見表1。姚瑞珍[13]通過室內(nèi)試驗方法確定了CODM。和氨氮在河流中的初始降解系數(shù)k分別為0.053 6和0.034 9，并基于實測楊湖斷面水質(zhì)，通過敏感度分析及手動調(diào)整參數(shù)，獲得了最佳擬合參數(shù)，見表1。2011年楊湖斷面水質(zhì)模擬值和實測值見圖5，其中COD模擬的相對誤差為28%，氨氮模擬的相對誤差為10%.水質(zhì)模型基本符合實際情況，可用于沙潁河安徽段的水環(huán)境容量計算。

數(shù)值模型網(wǎng)格剖分方式及尺寸與水流模型相同，計算節(jié)點位置與水流交叉分布。模擬時間步長與計算時段亦與水流模型設(shè)置相同。以往水環(huán)境容量計算以水質(zhì)目標值代替上游來水水質(zhì)的真實值，忽略了上游輸入的污染負荷對計算河段環(huán)境容量的影響。研究區(qū)上游為豫皖省界斷面，該斷面水質(zhì)的變化對安徽段水環(huán)境容量有不可忽略的影響，因此以上游豫皖省界斷面每周觀測的水質(zhì)為模型上游邊界的濃度進行計算。

為獲取研究區(qū)河流的背景值（考慮上游污染負荷輸入）.以上游豫皖省界斷面實測水質(zhì)為模型上游邊界，假設(shè)區(qū)內(nèi)排污量為0，計算2011-2012年河流污染物濃度。由于研究區(qū)實測水質(zhì)資料缺乏，無法給定符合實際的初始濃度，因此采用數(shù)值模擬方法估算初始條件影響范圍，以此消除初始條件對水質(zhì)模擬的影響。具體方法為，假設(shè)4組COD和氨氮的初始濃度值（C0）分別為1、5、10、15 mg/L。采用上述邊界條件及參數(shù)，計算人淮口濃度時間序列（見圖6）。從圖6可以看出，當(dāng)COD和氨氮模擬時間為50 d時，初始條件的影響可忽略（濃度差別小于0.03 mg/L）。因此選取第51 d的模擬濃度值為初始條件，作為河流背景濃度值的估算值，用于河流氨氮和COD的遷移轉(zhuǎn)化過程模擬。

3 水環(huán)境容量計算

針對水功能區(qū)劃及水質(zhì)目標[8]，劃定3個水環(huán)境容量計算單元及相應(yīng)的控制斷面（見圖2），控制斷面的水質(zhì)目標均為Ⅲ類。以往的水環(huán)境容量研究大都采用單一的水文設(shè)計條件，而根據(jù)不同的水文條件計算出的結(jié)果差別較大，利用最枯流量計算的結(jié)果往往嚴重低于實際環(huán)境容量值[14]。為了解實際動態(tài)水環(huán)境容量動態(tài)特征，從2010-2014年日觀測數(shù)據(jù)中選取2012年作為水文設(shè)計條件，通過線性規(guī)劃模型及水質(zhì)模擬，計算出各排污口允許最大日排污量及河流日環(huán)境容量，累加可得排污口允許最大年排污量和河流月環(huán)境容量。

由于對流彌散方程是線性的，具有可疊加性，因此可以把多個排污口對控制斷面水質(zhì)的綜合影響劃分為單個排污口對其影響的線性疊加[15]。基于疊加原理，李適宇等[16]提出貢獻度的概念，結(jié)合線性規(guī)劃法求解環(huán)境容量。該方法的思路：設(shè)某個排污口（i）單位時間排放單位負荷，而其他排污口無負荷量排放，背景濃度值為0，采用一維水質(zhì)模型計算各水質(zhì)控制斷面（j）的濃度值，該濃度值即為i排污口對各斷面的污染貢獻度（a）;改變排污口，重復(fù)以上步驟，即可求出每個排污口的污染貢獻度：基于排污口的污染貢獻度，考慮河流背景濃度值，以控制斷面水質(zhì)目標為約束，建立線性規(guī)劃模型的約束條件，確定目標函數(shù)，即總排污負荷最大（環(huán)境容量）。具體目標函數(shù)為

4 結(jié)果與討論

4.1 排污口環(huán)境容量

采用MIKE11模型計算出各排污口對各控制斷面COD和氨氮的貢獻度，年均日貢獻度見表2、表3（其中0表示排污口不在此斷面控制范圍內(nèi)）。可以看出，排污口1#、2#、3#對控制斷面I的貢獻度基本相等;僅排污口6#對控制斷面Ⅱ有貢獻;所有排污口對控制斷面Ⅲ均有貢獻，其中排污口5#貢獻度最大，排污口6#貢獻度最小。排污口離控制斷面越近，貢獻度越大，即排污口對控制斷面的影響隨距離的減小而增大。

基于排污口對各斷面的貢獻度，利用線性規(guī)劃模型計算出2012年各排污口的環(huán)境容量，結(jié)果見圖7。沙潁河安徽段COD環(huán)境容量約2.39萬t/a，氨氮環(huán)境容量約0.76萬t/a。其中：排污口5#的COD和氨氮容量最大，分別為7 892、2 970 t/a，占總水環(huán)境容量的38.90/0和32.9%;排污口3#的COD和氨氮環(huán)境容量最小，分別為948、165 t/a，僅占總水環(huán)境容量的3.9%和2.2%。原因是：排污口5#主要受斷面Ⅲ控制，對斷面I和斷面Ⅱ的貢獻度為0，受到的限制較小;而排污口3#對斷面I和斷面Ⅲ的貢獻度較大，受到的限制較大。排污口4#、5#、6#的環(huán)境容量較大，而排污口1#、2#、3#的環(huán)境容量較小，原因是上游來水（省界斷面）COD和氨氮濃度較大，大部分時期超過Ⅲ類水質(zhì)標準，余留的環(huán)境容量小，而隨水流過程中污染物的降解和稀釋，下游的環(huán)境容量逐漸增大。

4.2 月環(huán)境容量變化分析

2012年沙潁河安徽段月水環(huán)境容量及河流月徑流量見圖8。從圖8可知，COD和氨氮水環(huán)境容量年內(nèi)變幅較大，豐水期（7-10月）環(huán)境容量較大，枯水期（11月至翌年5月）環(huán)境容量較小。COD和氨氮環(huán)境容量均在8月最大，分別為3 557、1363 t;COD和氨氮環(huán)境容量在3月最小，分別為621、181 t;氨氮和COD環(huán)境容量的變化趨勢與徑流量不完全一致，1-7月徑流量逐漸減少，而COD和氨氮環(huán)境容量逐漸增大，直到8月隨徑流量增大而達到頂峰，說明水環(huán)境容量受徑流量、上游來水水質(zhì)等因素的綜合影響。

4.3 省界水質(zhì)對環(huán)境容量的影響

為分析省界斷面水質(zhì)對沙潁河安徽段環(huán)境容量的影響，設(shè)定3種情景討論環(huán)境容量的變化。情景1為省界斷面COD和氨氮濃度為0（即無負荷補給）、背景CODM。和氨氮濃度為0，計算水環(huán)境容量（L0）;情景2為省界斷面水質(zhì)為Ⅲ類水（水質(zhì)目標），計算水環(huán)境容量（L）;情景3為省界斷面COD氨氮濃度為實測值（實際情況），計算COD、氨氮環(huán)境容量（L）。

3種情景月水環(huán)境容量見圖9。在無上游污染負荷輸入的情景1下，環(huán)境容量最大，COD、氨氮年容量分別為36.0萬、9.1萬t：上游污染負荷達標排放情景下，環(huán)境容量較小，河段的環(huán)境容量主要來自CODM。和氨氮在河流中天然降解量.COD、氨氮年容量分別為17.5萬、5.5萬t;情景3，環(huán)境容量最小，COD、氨氮年容量分別為1.96萬、0.62萬t.水環(huán)境容量僅約為情景2的9%。COD和氨氮的環(huán)境容量均為L0>L>L，說明上游水質(zhì)越差，水環(huán)境容量越小。

2012年豫皖省界斷面水質(zhì)情況見圖10.由圖10可知，COD的超標率在泉河省控斷面達到55%，在潁河省控斷面達到29%：氨氮的超標率在泉河省控斷面達到18%.在潁河省控斷面達到46%，水質(zhì)達標率較低。其中，COD超標時間主要為1-6月，氨氮超標時間主要為1-6月（潁河省控斷面）和6-9月（泉河省控斷面）。由此可見，水環(huán)境容量除受河流流量控制外，還受上游水質(zhì)的影響。

5 結(jié)論

基于河流流量及水質(zhì)觀測數(shù)據(jù)，利用MIKE11水流及水質(zhì)模型，結(jié)合線性規(guī)劃法，計算沙潁河流域安徽段水環(huán)境容量，并分析水環(huán)境容量影響因素，結(jié)果表明：MIKE11模型對水文和水質(zhì)的模擬基本符合實際水文情況，模擬結(jié)果與實測數(shù)據(jù)基本一致;沙潁河安徽段COD№、氨氮年環(huán)境容量總和約為2.39萬、0.76萬t/a，COD和氨氮水環(huán)境容量年內(nèi)變幅較大，8月容量最大，分別為3 557、1 363 t，3月容量最小，分別為621、181 t：上游水質(zhì)對沙潁河安徽段水環(huán)境容量影響較大，實際水環(huán)境容量僅為上游水質(zhì)達標（Ⅲ類水）前提下理論水環(huán)境容量的9%。

參考文獻：

[1] 張玉清，張?zhí)N華，張景霞，河流功能區(qū)水污染物容量總量控制的原理和方法[J].環(huán)境科學(xué)，1998，19（增刊1）：23-33.

[2] 張永良，水環(huán)境容量基本概念的發(fā)展[J].環(huán)境科學(xué)研究，1992，5（3）：59-61.

[3]孟偉，張楠，張遠，等，流域水質(zhì)目標管理技術(shù)研究（I）：控制單元的總量控制技術(shù)[J].環(huán)境科學(xué)研究，2007，20 （4）：1—8.

[4]姜欣，許士國，練建軍，等，北方河流動態(tài)水環(huán)境容量分析與計算[J].生態(tài)與農(nóng)村環(huán)境學(xué)報，2013，29（4）：409-414.

[5]鄭孝宇，褚君達，朱維斌，河網(wǎng)非穩(wěn)態(tài)水環(huán)境容量研究[J].水科學(xué)進展，1997，8（1）：25-31.

[6]周孝德，郭瑾瓏，程文，等，水環(huán)境容量計算方法研究[J].西安理工大學(xué)學(xué)報，1999，15（3）：1-6.

[7] 董飛，劉曉波，彭文啟，等，地表水水環(huán)境容量計算方法回顧與展望[J].水科學(xué)進展，2014，25（3）：451-463.

[8]郝守寧，潁河流域水功能區(qū)面源污染控制研究[D].鄭州：華北水利水電大學(xué)，2014：1-2.

[9]逄勇，陸桂華，水環(huán)境容量計算理論及應(yīng)用[M].北京：科學(xué)出版社，2010：4.

[10] 中國環(huán)境規(guī)劃院，全國水環(huán)境容量核定技術(shù)指南[R].北京：中國環(huán)境規(guī)劃院，2003：46-76.

[11] 陳煉鋼，施勇，錢新，等，閘控河網(wǎng)水文一水動力一水質(zhì)耦合數(shù)學(xué)模型：Ⅱ應(yīng)用[J].水科學(xué)進展，2014，25（6）：856-863.

[12] 劉玉年，施勇，程緒水，等，準河中游水量水質(zhì)聯(lián)合調(diào)度模型研究[J].水科學(xué)進展，2009，20（2）：177-183.

[13] 姚瑞珍，人工控制的半封閉河流水環(huán)境容量測算研究：以老府河為實例[D].武漢：華中科技大學(xué)，2006：34-43.

[14] 李紅亮，王樹峰，不同設(shè)計水文條件下河北省水功能區(qū)納污能力研究[J].南水北調(diào)與水利科技，2010，8（3）：68-70.

[15] LIS Y， TOHRU M. Optimal Allocation of Waste Loadsin aRiver with Probabilistic Tributary Flow Under TransverseMixing[ J]. Water Environment Research， 1999， 71（2）： 156-162.

[16] 李適宇，李耀初，陳炳祿，等，分區(qū)達標控制法求解海域環(huán)境容量[J].環(huán)境科學(xué)，1999，20（4）：97-100.