安陽河動態納污能力分析

2021-08-11 15:22:44耿文斌趙姍侯俊山閆壽松李曉峰朱文舉

人民黃河 2021年7期

耿文斌趙姍侯俊山閆壽松李曉峰朱文舉

摘要：在不同保證率和不同時間尺度的設計水文條件下，動態計算出安陽河市區段2016年水體納污能力和年內分配，并對計算結果進行了對比分析，結果表明：安陽河市區段整體處于污染狀態，3個水功能區中僅安陽河安陽市排污控制區污染物排放不超標，其余水功能區COD和NH3-N排放量均超過納污能力上限，且NH3-N超標情況更為嚴重;安陽河市區段汛期月均納污能力與非汛期月均納污能力相差較小;安陽河市區段納污能力1—5月整體呈下降趨勢，5—7月呈逐漸上升趨勢，至8月納污能力達到最大，隨后又呈現波動下降趨勢，5月為全年最小;納污能力的大小與時間尺度成負相關關系，不同時間尺度下確定的年納污能力月尺度>水期尺度>年尺度。

關鍵詞：納污能力;保證率;時間尺度;安陽河

中圖分類號：X522;TV212.5 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2021.07.021

引用格式：耿文斌，趙姍，侯俊山，等.安陽河動態納污能力分析[J].人民黃河，2021，43（7）：110-113.

Abstract：Under the designed hydrological conditions of different guarantee rates and different time scales， the paper dynamically calculated the sewage capacity and annual distribution of the urban section of Anyang River in 2016. The calculation results were comparatively analyzed. The results show that the pollutant discharge in Anyang sewage control area is not exceeding the standard， the COD and NH3-N discharge of other water function areas exceed the upper limit of the capacity and the NH3-N is more serious. The monthly average water environmental capacity of the river in flood season is smaller than that in non-flood season， the water environmental capacity of the river in July to October is larger， the water environmental capacity of August is the maximum in the whole year and the water environmental capacity of November to May is gradually reduced， in May， the water environmental capacity is the lowest in the whole year. There is a negative correlation between the size of the capacity and the time scale and the annual capacity is the monthly scale > the water scale > the annual scale.

Key words： water environmental capacity; guarantee rate; time scale; Anyang River

納污能力是指一定條件下水體的最大納污量[1-2]。納污能力是制定污染物限排量的重要依據，能直接反映水體的水質狀況[3-4]。目前我國納污能力大多按照《水域納污能力計算規程》（GB/T 25173—2010）計算，其中設計水文條件選取90%保證率最枯月平均流量或近10 a最枯月平均流量，這種條件下得到的納污能力為一固定值，且設計水文條件較為嚴格，以此作為標準應用于水資源管理和保護必然安全[5]。但實際納污能力與水文過程關系密切，水文過程在自然條件下是動態變化的，納污能力也是一個動態變化的量[6-7]。將一個確定的值作為控制指標難以揭示水體的實際污染負荷，不利于保護和改善水質，也不利于水資源的合理開發利用，因此開展水體動態納污能力評價很有必要[8]。

我國對動態納污能力的研究眾多，但多以不同保證率、不同水文條件來計算水體動態納污能力。李喬臻[9]以95%、90%、75%、50%、20%保證率下的流量為設計流量，選取不同來水水質分析了琶江河段3個水功能區的納污能力;周洋等[10]運用一維穩態水質模型和水環境容量模型，根據水質、水文監測資料和沿岸排污情況，計算了渭河陜西段豐水年、平水年、枯水年COD的水環境容量。這些研究都只分析了不同保證率下年納污能力的動態變化特征，并未結合年內水文特性進行分析，只能確定水功能區全年的污染物限排量，年內的污染物限排量分配仍不清晰，以此為水資源管理的依據會出現全年污染物排放量不超標，但某些月份因污染物排放量過多而造成水功能區水質變差的現象。筆者基于安陽水文站48 a（1969—2016年）長系列資料，在不同保證率條件下，根據不同時間尺度確定設計水文條件，以COD和NH3-N為控制因子，建立河流一維水質模型[11-13]，計算2016年安陽河市區段年、水期和月納污能力，以期為安陽河確定污染物限排總量、年內污染物限排量分配提供數據支撐，為安陽河管理與保護提供科學依據。

1 研究區概況

安陽河是海河流域漳衛南運河水系的第二大支流，發源于林州市清泉寺，自西向東流經林州市、安陽縣、安陽市區、內黃縣（見圖1），在內黃縣李大晁入衛河，全長約164 km，流域面積1 920 km2。上游穿過林州盆地后，在橫水鎮附近進入山區，于彰武水庫以下進入平原。干流上建有彰武水庫和小南海水庫。安陽河的主要補給水源包括巖溶泉水（泉主要有小南海泉和珍珠泉，泉域面積935 km2）和大氣降水，據安陽水文站多年觀測資料統計，最大流量可達2 060 m3/s（1982年8月2日），最小流量為0.29 m3/s（1982年7月8日），多年平均流量為10.73 m3/s[14]。

安陽河市區段共包括3個二級水功能區，分別為安陽河彰武水庫飲用水源區、安陽河安陽縣蔣村用水區和安陽河安陽市排污控制區，總長約49.1 km，水質目標分別為Ⅱ類、Ⅲ類和Ⅳ類。

2 納污能力計算模型

2.1 計算模型

本次計算單元為彰武水庫大壩至入衛河口，共包含3個水功能區，計算單元示意見圖2。

安陽河總體上具有流速小、河道寬、水淺等特點，主要河道寬深比小，入河污染物在較短河段內混合均勻，污染物橫向變化不大。因此，可采用概化一維水質模型進行納污能力計算，計算公式如下：

M=31.536（QP+Q）（CseKx86.4u-C0e-KL86.4u）（1）

式中：M為水功能區納污能力，t/a;Q為水功能區設計流量，m3/s;QP為廢污水排放量，m3/s;Cs為水功能區水質控制目標質量濃度，mg/L;C0為水功能區初始斷面質量濃度，mg/L;K為污染物綜合衰減系數，1/d;x為概化點距控制斷面的距離，m;L為水功能區長度，km;u為設計流量下斷面平均流速，m/s。

根據研究區排污口分布以及資料掌握情況[15]，將各個水功能區內多個排污口概化為一個集中排污口，且集中排污口位于河段中點處，因此納污能力計算公式改為

M=31.536（QP+Q）（CseKL2×86.4u-C0e-KL86.4u）（2）

2.2 模型參數

2.2.1 設計水文條件的確定

水體納污能力動態變化受多因素共同影響，而設計水文條件是納污能力計算的重要前提，對計算結果影響較大?；诎碴査恼径嗄瓯O測數據，分別確定不同時間尺度和不同保證率下的設計水文條件，并對納污能力的變化特征進行分析。

（1）基于年尺度的設計水文條件的確定。選取系列資料中逐年最枯月平均流量，采用P-Ⅲ型曲線參數估計法繪制最佳經驗頻率曲線，確定50%、75%、90%保證率下的設計流量分別為2.94、1.76、1.03 m3/s，根據河道大斷面狀況分析得到對應的設計流速。

（2）基于水期尺度的設計水文條件的確定。將年內月份分為汛期（6—10月）和非汛期（11月至翌年5月），選取逐年汛期和非汛期內的最枯月平均流量，采用P-Ⅲ型曲線參數估計法，分別繪制汛期和非汛期的最佳經驗頻率曲線，確定75%保證率下的設計流量分別為2.06、1.88 m3/s，90%保證率下的分別為1.16、1.07 m3/s，并根據河道大斷面狀況分析得到對應的設計流速。

（3）基于月尺度的設計水文條件的確定。選取逐年各月平均流量，采用P-Ⅲ型曲線參數估計法分別繪制1—12月的最佳經驗頻率曲線，確定90%保證率下1—12月的設計流量分別為2.13、2.08、1.88、1.52、1.37、1.48、1.98、2.46、2.24、2.33、2.16、2.12 m3/s，再根據河道大斷面狀況分析得到對應的設計流速。

根據經驗頻率曲線分析結果可知，相同保證率下，時間尺度越小，對應的設計流量越大。主要原因是設計水文條件的確定方法不同，年尺度下選擇逐年最枯月平均流量進行計算，水期尺度則分別選擇逐年汛期最枯月平均流量和非汛期最枯月平均流量進行計算，而安陽河屬于常年河，枯水年和平水年沒有明顯的汛期特征，甚至個別年份汛期內月平均流量小于非汛期內月平均流量，因此汛期和非汛期確定的設計流量更大。月尺度分別選擇對應的月平均流量作為系列值進行分析，確定的設計流量更大。

2.2.2 水質濃度及綜合衰減系數的確定

水功能區初始斷面質量濃度C0采用實測數據;水功能區水質目標質量濃度Cs以水質類別為控制指標。根據《全國地表水環境容量核定技術復核要點》中提供的河道水質綜合衰減系數參考值，結合各水功能區優劣狀況和水文特征選取綜合衰減系數。由經驗法，COD的綜合衰減系數計算公式為

K=0.05+0.68u（3）

式中：u為設計流速。

NH3-N的綜合衰減系數計算公式為

K=0.061+0.551u（4）

2.3 數據來源

安陽河水文數據來源于安陽水文站長系列監測數據;污染物排放數據來源于安陽市水利局和安陽水文水資源勘測局共同出版的《安陽市入河排污口監測報告》，水功能區污染物入河量統計見表1。

3 納污能力計算分析

3.1 年尺度納污能力計算分析

安陽河市區段年尺度納污能力計算結果見表2，由表2可以看出，安陽河市區段3個水功能區中安陽河安陽市排污控制區納污能力最大，安陽河彰武水庫飲用水源區納污能力最小;當上游來水水質優于水功能區水質目標，其他因素一定時，納污能力與保證率的大小成負相關關系，即納污能力隨保證率的增大而變小。

3.2 水期尺度納污能力計算分析

安陽河市區段水期尺度納污能力計算結果見表3，由表3可以看出，安陽河市區段汛期月均納污能力與非汛期月均納污能力相差較小。主要原因是安陽河為常年河，在豐水年，徑流量較為集中，主要集中在6—10月，占年徑流量的60%～80%，平水年、枯水年年內徑流分配較為均勻，汛期特征不明顯，特枯水年基本無明顯的汛期特征，且汛期安陽河上游進行灌溉、水庫蓄水等，導致汛期月均流量與非汛期月均流量相差較小。

3.3 月尺度納污能力計算分析

安陽河市區段月尺度納污能力計算結果見表4，安陽河市區段COD和NH3-N年納污能力分別為1 393.4 t和91.5 t。由表4可知，安陽河市區段逐月納污能力呈波動變化趨勢，并未與流量季節變化保持完全一致，說明安陽河納污能力不僅受自然水文條件的影響，還受到人類活動的影響，且人類活動對納污能力的影響較大，如6月為汛期，但受人類活動影響（上游灌溉、水庫蓄水等），其納污能力較5月增大不明顯。安陽河市區段納污能力1—5月整體呈下降趨勢，5—7月呈逐漸上升趨勢，至8月納污能力達到最大，隨后又呈現波動下降趨勢。

3.4 綜合對比分析

由不同尺度（年、水期、月）納污能力的對比分析可知，納污能力與時間尺度成負相關關系，不同時間尺度下確定的年納污能力月尺度>水期尺度>年尺度。

主要原因是，年尺度納污能力計算分析時考慮多年的不利水文條件，以此為全年的水文條件進行納污能力計算，這種計算方法雖然對水環境的改善有促進效果，但計算尺度較大，確定的設計水文條件不具備代表性，水體水質狀況大多數時段優于設計水文條件，確定的年納污能力較小。隨著計算尺度的減小，設計水文條件以水期、月為尺度確定，更接近水體實際情況，計算得到的納污能力更接近水體污染承載力。

由確定的納污能力與年污染物入河量進行對比分析，安陽河市區段整體處于污染狀態，3個水功能區僅安陽河安陽市排污控制區污染物排放不超標，其余水功能區COD和NH3-N排放量均超過納污能力上限，且NH3-N超標更嚴重。主要原因是，水域周圍農業生產及養殖活動較多，廢水未經處理直接排入河道，因此研究區在控制污染物排放時，對NH3-N排放的限制應更嚴格。

4 結語

納污能力與水體水文過程關系密切，是一個動態變化的量，筆者在不同保證率和不同時間尺度下對安陽河市區段動態納污能力進行了計算，分析了不同保證率和不同時間尺度對納污能力的影響。考慮保證率和時間尺度的變化，可以使納污能力的計算結果與水體實際狀況更相符。目前水功能區存在監測時間間隔過長的問題，且以特定的監測點水質分析結果來判斷整個水功能區的水質狀況，存在極大的偶然性，因此建議以污染物排放量與納污能力比值判斷水功能區是否達標，并分別從月、水期和年尺度進行納污能力考核，以避免傳統考核模式存在的偶然性問題。本文僅進行了年、水期和月尺度的納污能力動態分析，今后可考慮更小時間尺度（日、時）的計算，以了解水體的實際狀況，及時對水資源管控做出調整，提高水功能區納污能力。

參考文獻：

[1] 李如忠，汪家權，王超，等.不確定性信息下的河流納污能力計算初探[J].水科學進展，2003，14（4）：359-363.

[2] 趙鑫，黃茁，李青云.我國現行水域納污能力計算方法的思考[J].中國水利，2012（1）：29-32.

[3] 王林，宓辰羲，宓永寧.遼陽市水功能區納污能力計算研究[J].人民黃河，2017，39（12）：80-84.

[4] 張秀菊，楊凱，蔡愛芳，等.不確定性參數對水體納污能力的影響分析[J].中國農村水利水電，2012（1）：13-17.

[5] 中華人民共和國質量監督檢驗檢疫總局，中國國家標準化委員會.水域納污能力計算規程：GB/T 25173—2010[S].北京：中國水利水電出版社，2011：1-4.

[6] 張海歐，閔濤，羅軍剛，等.渭河干流納污能力動態計算模型系統開發[J].人民黃河，2013，35（6）：70-72.

[7] 張曉，羅軍剛，杜金龍，等.水功能區納污能力動態計算主題模式研究及應用[J].水利信息化，2016（2）：24-28.

[8] 史曉新，禹雪中，馬巍.湖泊納污能力動態特征分析及計算[J].中國水利水電科學研究院學報，2008，6（2）：105-110.

[9] 李喬臻.琶江河段動態納污能力分析[J].水文，2009，29（增刊1）：131-133.

[10] 周洋，周孝德，馮民權.渭河陜西段水環境容量研究[J].西安理工大學學報，2011，27（1）：7-11.

[11] 張川，呂謀.河道納污模型研究與應用[J].青島理工大學學報，2015，36（1）：74-79，113.