樂至縣蟠龍河水庫泄庫水質影響預測及跟蹤評估研究

郭明昆 鄒俊良 周安瀾 楊長軍 湯冰冰

（1.四川省生態環境科學研究院，四川 成都 610064；2.四川省環科院科技咨詢有限責任公司，四川 成都 610064；3.四川省環保科技工程有限責任公司，四川 成都 610041）

1 引言

樂至縣蟠龍河水庫位于瓊江支流蟠龍河上游，1959 年建成，大壩為均質土壩，是一座以防洪、灌溉為主，兼發電、水產等綜合利用的中型水庫。蟠龍河水庫因安全問題被水利部大壩安全管理中心核準為三類壩，為消除安全隱患，樂至縣計劃將蟠龍河水庫騰空后進行除險加固。蟠龍河水庫下游有蟠龍河元壩子、瓊江跑馬灘、瓊江大安3 個國控斷面，泄庫水質必然對下游國控斷面水質產生較大影響。本研究擬開展蟠龍河水庫泄水對下游考核斷面水質影響的預測研究，以期為流域的水環境管理提供技術參考和決策依據。

2 數據來源與方法

2.1 蟠龍河水庫概況

蟠龍河水庫屬中型水庫，距離樂至縣城40 km，水庫總庫容4 080 萬m3，正常蓄水位庫容2 100 萬m3，死庫容100 萬m3，水庫正常蓄水位332.00 m，水庫死水位317.43 m。蟠龍河水庫多年平均年徑流深為230.0 mm，多年平均枯水期（11 月—次年4 月）徑流深為22.0 mm；水庫多年平均流量1.36 m3/s，多年平均枯水期流量為0.14 m3/s。

2.2 水環境現狀數據來源

為全面掌握瓊江流域水環境質量狀況及蟠龍河水庫泄庫水質對考核斷面水質影響，本研究系統收集了蟠龍河水庫、元壩子斷面、跑馬灘斷面、大安斷面的水質監測資料，同時在蟠龍河水庫入庫口、庫心、庫尾設置了4 個加密監測點，其中，國控斷面水質監測數據來自樂至生態環境局提供的四川省監測站上報國家核定數據，監測指標包括水溫、pH、懸浮物、溶解氧、CODCr、BOD5、CODMn、NH3-N、TN、TP；采樣監測點數據來自資陽市樂至生態環境局提供的第三方監測數據，監測項目包括CODCr，CODMn，NH3-N，TP；泄庫跟蹤監測數據來自水質自動監測站，采用單因子分析法。按照各監測斷面水質要求，斷面水質對照GB 3838—2002《地表水環境質量標準》Ⅲ類標準限值進行評價。

2.3 水質預測模型

對污染物擴散輸移預測的有效方法就是通過水質模型進行計算［1］。根據污染物進入水體后沿垂直、縱向、橫向的輸移和擴散情況，水質預測模型又分為零維、一維、二維、三維模型［2-3］，其中，一維模型是污染物濃度緊隨縱向變化，適用于河流長度遠大于寬度和深度的情況。本研究中的蟠龍河、瓊江計算河段長度遠大于河道寬度，屬于Q＜150 m3/s 的中小型河段，因此，可采用一維水質模型模擬蟠龍河水庫下泄后河流水質變化情況。

在穩態或標準態的情況下，一維水質數學模型為［2，4-5］：

3 結果與分析

3.1 各斷面水環境質量狀況評價

3.1.1 蟠龍河水庫及下游國控斷面水環境質量狀況評價

根據資陽市樂至生態環境局提供的2021 年1—9 月的常規監測數據，蟠龍河水庫1—5 月達到Ⅲ類水質考核標準，6—9 月水質為Ⅳ類，主要超標污染物為COD 和CODMn；2021 年1—9 月元壩子斷面達到Ⅲ類水質考核標準的月份僅占44.44%，其中4 月、7—9 月均為Ⅳ類，5 月為Ⅴ類，1—9 月COD 平均值為21.33 mg/L，要實現全年水質達標，元壩子斷面在10—12 月COD 平均濃度需小于16 mg/L，水質達標壓力較大；跑馬灘斷面6 月和8 月水質為Ⅳ類，超標因子為COD 和CODMn，其余月份達到Ⅲ類水質考核標準，全年達標壓力較小；2021 年1—9 月大安斷面有3 個月水質未達標，超標因子為COD，1—9 月COD 平均值為20.13 mg/L，要實現全年水質達標，大安斷面在10—12 月COD 平均濃度需小于19.6 mg/L，仍存在一定達標壓力。2021 年1—9 月蟠龍河水庫及下游國控斷面水環境質量狀況見圖1。

圖1 2021 年1—9 月蟠龍河水庫及下游國控斷面水環境質量狀況

3.1.2 蟠龍河水庫加密監測水質分析與評價

2021 年8 月21—23 日資陽市樂至生態環境局委托第三方檢測公司對入庫壩前、庫心、主壩及泄洪口進行連續3 d 跟蹤監測，根據監測結果，入庫壩前水質優于庫心、主壩及泄洪口等位置，水庫內各點水質NH3-N 和TP 中下層濃度均高于表層水體濃度，而水庫內各點COD 和CODMn未呈現出明顯差異性，但越靠近水庫下游，表層水體水質越優于中下層水體，表面輸入的污染物在水庫發生了沉降或庫內底泥在底層水體向上釋放污染物。蟠龍河水庫加密監測統計結果見表1。

表1 蟠龍河水庫加密監測統計結果mg/L

3.2 水質模型預測及達標影響分析

3.2.1 邊界條件

3.2.1.1 水域概化

本研究模型預測的河流概化主要考慮蟠龍河、跑馬灘水庫及瓊江干流，由于跑馬灘水庫屬于河流型水庫，水庫長12 km、寬120 m，長寬比值為100，在計算污染物擴散時可將其簡化為河段。根據河段劃分原則及旁側入流、污染源和常規監測斷面位置，對蟠龍河、瓊江河道進行簡化，如圖2 所示。

圖2 蟠龍河及瓊江河道概化示意

3.2.1.2 控制因子

根據流域水污染現狀和水污染物總量控制現狀，選擇CODCr，NH3-N 和TP 作為模型預測的控制因子。

3.2.1.3 設計水文條件

參考《瓊江流域水體達標方案》中瓊江流域多年平均流量條件下各計算單元的流量計算成果及河流水文站統計的多年平均流量數據，計算得出流速等參數。多年平均流量下各計算單元水文條件見表2。

表2 多年平均流量下各計算單元設計水文條件

根據樂至縣水務局提供的泄水方案，蟠龍河水庫將按最大下泄量5.6 m3/s 進行泄水，河道粗糙系數取0.027，最大下泄流量下各斷面的河流水文特征見表3。

表3 蟠龍河水庫最大下泄流量下各計算單元設計水文條件

3.2.1.4 水質目標

瓊江流域主要河流水環境功能類別為Ⅲ類，根據“十四五”水環境考核目標，元壩子、跑馬灘、大安3 個國考斷面水質目標為Ⅲ類。

3.2.1.5 不同河段污染物濃度值

（1）在多年平均流量情景下，起始河段蟠龍河水庫C0值取2021 年1—9 月監測水質濃度的算術平均值，COD，NH3-N，TP 濃度為20.11，0.24，0.065 mg/L；元壩子、跑馬灘、大安斷面取2021 年1—9 月的四川省環境監測站上報國家的水質數據的算術平均值，其中，元壩子斷面COD，NH3-N，TP 濃度分別為21.33，0.08，0.09 mg/L，跑馬灘斷面COD，NH3-N，TP濃度分別為17.79，0.05，0.03 mg/L，大安斷面COD，NH3-N，TP 濃度分別為20.13，0.05，0.10 mg/L。（2）在最大下泄流量條件下，起始河段蟠龍河水庫C0值取2021 年8 月21—23 日蟠龍河水庫庫心、主壩、泄洪口等部位水質的算術平均值，COD，NH3-N，TP 濃度分別為24.17，1.55，0.16 mg/L。

3.2.1.6 降解系數

根據《全國水環境容量核定技術指南》和同類型項目相關研究［6-8］，本研究根據流速確定的降解系數取值見表4。

表4 降解系數相關研究成果

3.2.2 不同情景下模型預測結果

3.2.2.1 多年平均流量條件下外源污染輸入負荷核算

選取多年平均流量條件下的各段水文參數，起始河段蟠龍河水庫C0值取2021 年1—9 月監測水質濃度的算術平均值，采用河流一維模型計算蟠龍河及瓊江干流不同范圍河段水質，對計算值與斷面實測濃度的差值考慮為施加到河段內的污染負荷。根據計算，元壩子斷面河段施加的CODCr，NH3-N，TP 污染物分別為236.031，0.000，2.793 kg/d，跑馬灘斷面河段施加的CODCr，NH3-N，TP 污染物分別為10 421.52，23.21，16.27 kg/d，大安斷面河段施加的CODCr，NH3-N，TP 污染物分別為4 355.61，6.32，53.94 kg/d。多年平均流量下各計算單元主要污染物濃度預測值及外源輸入負荷見表5。

3.2.2.2 蟠龍河水庫最大泄庫流量條件下污染物濃度預測

選取蟠龍河水庫最大下泄流量條件下的各段水文參數，起始河段蟠龍河水庫C0采用2021 年8 月21—23 日蟠龍河水庫庫心、主壩、泄洪口等部位水質的算術平均值，采用河流一維模型計算該條件下蟠龍河及瓊江干流不同范圍河段水質，考慮各河段施加污染負荷后，元壩子斷面CODCr，NH3-N，TP 的預測濃度分別為24.12，1.52，0.16 mg/L，跑馬灘斷面CODCr，NH3-N，TP 的預測濃度分別為18.29，0.68，0.09 mg/L，大安斷面CODCr，NH3-N，TP 的預測濃度分別為19.26，0.61，0.13 mg/L。根據預測結果，蟠龍河水庫泄庫期間，元壩子斷面平均水質為地表水Ⅳ類，跑馬灘斷面、大安斷面平均水質能夠達到地表水Ⅲ類。蟠龍河水庫最大泄庫流量下各單元主要污染物濃度預測值見表6。

3.2.2.3 跑馬灘及大安斷面達標情況下水庫泄庫水質控制濃度

選取蟠龍河水庫最大下泄流量條件下的各段水文參數，采用河流一維模型計算該條件下蟠龍河及瓊江干流不同范圍河段水質，經計算，當蟠龍河水庫泄庫尾水日平均濃度CODCr超過28.67 mg/L，或NH3-N超過2.32 mg/L，或TP 超過0.40 mg/L 時，預測跑馬灘斷面水質不能達標，見表7。

表7 蟠龍河水庫最大泄庫流量下跑馬灘斷面不達標初始臨界濃度控制值mg/L

當蟠龍河水庫泄庫尾水日平均濃度CODCr超過26.45 mg/L，或NH3-N 超過2.56 mg/L，或TP超過0.32 mg/L 時，預測大安斷面水質超標，見表8。

表8 蟠龍河水庫最大泄庫流量下大安斷面不達標初始臨界濃度控制值mg/L

根據表1 實測水質統計，蟠龍河水庫泄庫期間，元壩子斷面不能達到Ⅲ類水的概率為100%，跑馬灘、大安斷面不能達到Ⅲ類水的概率低于30%。

3.2.2.4 蟠龍河水庫泄庫輸出污染物通量對各斷面水質貢獻

根據擬定的蟠龍河水庫泄水方案，計劃從12 月開始泄水，放水歷時57 d，不考慮庫底清淤等工程額外帶來的污染物量，計算得出蟠龍河水庫污染物的通量為：CODCr507.57 t，NH3-N 32.55 t，TP 3.36 t，蟠龍河水庫泄庫污染物輸送到元壩子斷面的總量為：CODCr496.20 t，NH3-N 32.00 t，TP 3.34 t；跑馬灘斷面的總量為：CODCr351.09 t，NH3-N 25.12 t，TP 2.81 t；輸送到大安斷面的總量為：CODCr298.66 t，NH3-N 22.25 t，TP 2.70 t。由于各情景下河流斷面主要為COD，TP 偶爾超標，故計算通量時主要考慮COD 和TP 帶來的影響，蟠龍河水庫泄水帶來的污染物在各個斷面污染物總量的比例分別為：元壩子斷面為70.99%～79.95%，跑馬灘斷面為29.96%～67.00%，大安斷面為24.33%～38.04%。由此可知，蟠龍河水庫泄庫對下游水質的影響隨距離的延長而減小，離水庫最近的元壩子斷面受影響最大。蟠龍河水庫泄庫輸出污染物通量對各斷面水質貢獻情況見表9。

表9 蟠龍河水庫泄庫輸出污染物通量對各斷面水質貢獻情況

3.2.3 水庫泄庫水質影響跟蹤評估研究

2022 年3 月7 日至4 月17 日蟠龍河水庫實施了泄庫，本研究開展了泄庫期間元壩子、跑馬灘和大安斷面的水質跟蹤評估，采用泄庫期間元壩子、跑馬灘和大安斷面的在線監測數據，由于在線監測站沒有COD 監測數據，所以采用同為有機污染程度表征的CODMn進行評估，對CODMn，NH3-N，TP 等主要指標作趨勢圖，見圖3。元壩子斷面未能達到Ⅲ類水的概率為92.86%，其主要超標因子為CODMn和NH3-N；跑馬灘斷面和大安斷面未能達到Ⅲ類水的概率分別為28.57%，11.90%，主要超標因子均為CODMn。跟蹤研究表明，泄庫期間實際水質變化與預測概率相符。

圖3 蟠龍河水庫泄庫期間國控斷面水質實際情況

4 結論與建議

（1）在蟠龍河水庫最大下泄流量的條件下，根據模型預測元壩子斷面不能達到Ⅲ類水質的比例為100%，跑馬灘和大安斷面平均水質能夠達到地表水Ⅲ類；當蟠龍河水庫泄庫尾水日平均濃度CODCr超過28.67 mg/L，或NH3-N 超過2.32 mg/L，或TP 超過0.40 mg/L 時，預測跑馬灘斷面水質不能達標；當蟠龍河水庫泄庫尾水日平均濃度CODCr超過26.45 mg/L，或NH3-N 超過2.56 mg/L，或TP 超過0.32 mg/L 時，預測大安斷面水質超標，預測跑馬灘斷面和大安斷面不達標的概率低于30%。

（2）預測蟠龍河水庫泄庫污染物輸送至元壩子斷面的總量為：CODCr496.20 t，NH3-N 32.00 t，TP 3.34 t；跑馬灘斷面的總量為：CODCr351.09 t，NH3-N 25.12 t，TP 2.81 t；輸送到大安斷面的總量為：CODCr298.66 t，NH3-N 22.25 t，TP 2.70 t。因蟠龍河水庫泄水帶來的污染物在各個斷面污染物總量的比例分別為：元壩子斷面為70.99%～79.95%，跑馬灘斷面為29.96%～67.00%，大安斷面為24.33%～38.04%。由此可知，蟠龍河水庫泄庫對下游水質的影響隨距離的延長而減小，離水庫最近的元壩子斷面受影響最大。

（3）對蟠龍河水庫泄庫期間3 個國控斷面水質的跟蹤研究發現，下游元壩子斷面未能達到Ⅲ類水的概率為92.86%，跑馬灘斷面未能達到Ⅲ類水的概率為28.57%，大安斷面未能達到Ⅲ類水的概率為11.90%，泄庫期間實際水質變化與預測結果相符。

（4）建議根據環境監測考核要求在蟠龍河水庫上游設置臨時替代斷面，作為施工期間及蓄水恢復期元壩子斷面的替代考核斷面；建議在施工期間及蓄水恢復期通過加強上游預警監測及瓊江流域水庫協同放水調度來降低跑馬灘斷面、大安斷面水質超標的概率；建議加強泄庫后蓄水恢復區水生態環境影響研究。