引水改善平原感潮河網水質效果評估

許益新，李一平，羅育池，余香英，蔣婧媛，羅 凡

(1.廣東省環境科學研究院，廣東 廣州 510045； 2.河海大學環境學院，江蘇 南京 210098)

平原感潮河網地勢平坦，水流緩慢，同時受潮汐作用和泵閘調度的影響，往復流現象明顯，抗污染能力較弱。隨著經濟的快速發展，平原感潮河網地區水資源高度開發利用，水環境質量日益下降[1]。利用平原感潮河網緊鄰大江大河的優勢，根據潮汐變化進行引水、排水，構建科學有序的水循環體系，促進河道水體流動，提高水體稀釋能力，從而達到改善水環境質量的效果[2]。水環境數學模型是用來模擬和預測污染物混合、遷移、轉化過程的重要工具之一，在定量分析引水改善河網水環境質量效果方面發揮著重要作用[3]。

目前，對于引水改善河網地區水質效果的研究多集中于特定引水流量[4]、引水水質[5]、水系結構優化[6]等，顧爐華等[7]基于太湖水量水質數學模型研究不同引水量和雨型下七浦塘引水對陽澄湖和周邊河網水質的改善效果，結果表明引水量為20～40 m3/s時，河網水質改善效果最佳；趙軒等[8]基于常州市河網水環境數學模型研究上游引水、污染源削減和提高自凈能力對河網水質的改善效果，結果表明保證上游來水水質良好的情況下，上游引水對主干河道水質有較好的改善效果；賈海峰等[9]基于城鎮河網水環境數學模型研究水系結構優化和水利調度對河網水質的改善效果，結果表明水系結構優化方案能夠降低BOD5和NH3-N濃度。目前對復雜的平原感潮河網地區潮汐作用下調水引流工作研究較少，且尚未建立起潮位、引水量和水質改善效果的響應關系。

本文選取典型平原感潮河網地區張家港市為研究對象，根據張家港市緊鄰長江的優勢和水循環體系的建設，開展調水引流改善平原感潮河網水質效果研究，建立“長江潮位-內河引水量-水質改善效果-引水服務面積”的響應關系，以期為平原感潮河網調水引流工作提供科學的指導和建議。

1 研究區域概況

張家港市位于長江三角洲沖積平原，境內河網密布，北鄰長江，屬典型的平原感潮河網地區，境內共有37個水功能區，隨著污染源削減方案的實施和水循環體系的構建，水功能區達標率逐步提高，2011—2016年水功能區達標率從57.4%提升到79.6%，主要超標因子為氨氮。

張家港市目前正在建設、完善中部、東部和西部三大水循環體系，其中中部水系水循環體系較為完善，主要依靠朝東圩港水利樞紐和一干河水利樞紐候潮開啟閘門進行自流引水，引水時內河側閘前控制水位為4.1 m，二干河為排水通道；東部水系主要依靠三干河水利樞紐、四干河節制閘和六干河節制閘候潮開啟閘門進行自流引水，引水時三干河內河側閘前控制水位為3.8 m，四干河和六干河內河側閘前控制水位為3.7 m，五干河和七干河為排水通道；西部水系主要依靠五節橋泵閘工程、十字港閘和護漕港節制閘候潮開啟閘門進行自流引水，引水時內河側閘前控制水位為4.0 m，太字圩港為排水通道。各引水口門距離長江天生港潮位站均在20 km以內，故選取天生港潮位站作為長江潮位代表站[10](圖1)。

2 平原河網非穩態數學模型構建及評價方法

2.1 模型的建立

在張家港市中部水系水動力水質耦合數學模型研究[11]的基礎上，采用一維圣維南方程組構建張家港市平原感潮河網水量模型離散格式為六點中心隱式差分格式，采用“雙掃”法計算求解；采用一維對流擴散方程作為水質模擬基本方程，根據長江三角洲河網和張家港市河網水質現狀、特點及總量控制要求，選取高錳酸鹽指數和氨氮為模擬因子，構建張家港市平原感潮河網水質模型，以評估張家港市三大水循環體系的構建對區域河網水質的改善效果[12]。河網模型共概化368條河流(市級河道、鎮級河道和部分村級河道)，151個閘站(沿江水利樞紐、內河節制閘和部分泵站)；河流斷面數據采用走航式ADCP監測資料進行概化；模型上邊界為各引水通道與長江交匯處，為流量邊界，流量時間序列根據各水文自動監測站點高頻率長序列監測數據分析所得；模型下邊界為各排水通道與長江交匯處，為水位邊界。

圖1 研究區域位置及評價斷面分布

2.2 模型的率定和參數的選取

根據張家港市水文自動監測站點實時監測數據和野外水量水質同步監測結果，對河網水動力水質模型進行率定，結果表明：張家港市現狀河道糙率為0.01～0.04，其中城區主要河道為直立式護岸，河道糙率較小，為0.01～0.02，鄉鎮河道兩側雜草樹木較多，河道糙率較大，為0.03～0.04；張家港市現狀高錳酸鹽指數降解系數為0.04～0.08 d-1，氨氮降解系數為0.03～0.08 d-1。各斷面流量模擬相對誤差在25%以內，平均相對誤差為13.1%；各斷面水質模擬相對誤差在30%以內，高錳酸鹽指數平均相對誤差為15.2%，氨氮濃度平均相對誤差為12.0%，具體見表1。

表1 張家港市河網模型水動力水質模擬相對誤差

2.3 調水引流改善水環境質量評價方法

根據張家港市水動力水質數學模型計算結果，評價調水引流改善水環境質量的方法主要有：①分析評價引水期間監測點位水質指標濃度隨時間變化趨勢；②分析評價引水前后監測點位水質指標濃度的改善效果；③分析評價引水前后監測點位水質類別的變化。因此，結合前人研究成果，采用水質濃度改善率、類別變化指數和濃度變化指數3個指標評估平原感潮河網水質改善效果與引水量的響應關系[13]。公式為

(1)

式中：R為水質濃度改善率；C0i為第i種水質指標引水前的濃度；Cai為第i種水質指標引水后的濃度；i為參加評估的水質指標數目。R值反映引水前后水質濃度改善的程度，正值表示水質濃度下降，負值表示水質濃度上升。

(2)

式中：G為類別變化均值指數；G0i為第i種水質指標引水前的水質類別；Gai為第i種水質指標引水后的水質類別。G屬于均值型多因子指數，反映引水前后水質類別變化的等級，數值為正表示水質類別改善，數值為負表示水質類別下降。

(3)

式中P為濃度變化指數。P屬于均值型多因子指數，反映引水前后水質濃度變化的程度，數值為正表示水質濃度下降，數值為負表示水質濃度上升。

3 結果與討論

3.1 平原感潮河網內河引水量與長江潮位響應關系

長江張家港段受到上游大通來水和沿海潮汐的雙重作用，每天有兩個高潮位和低潮位，屬于不規則半日潮，且高潮位具有明顯的半月周期特性，低潮位則沒有明顯的周期特征[14-15](圖2)。當長江潮位高于張家港市內河水位時，開啟沿江水利樞紐使長江水自流引進。半月周期內，隨著長江高潮位的升高，各引水口門引水量呈現線性增長趨勢(圖3)，其中朝東圩港引水量與長江潮位相關性最強，且引水量最大，一干河次之；由于三干河、四干河和六干河引水過程中受到內河防洪水位的限制，未能按照最大引水量進行引水，故引水量與長江潮位相關性較弱，且引水量較少。西部水系的五節橋港、十字港和護漕港引水通道尚未正式實施，故暫不討論。

圖2 天生港潮位站高、低潮位與內河水位示意圖

不同引水口門引水量除了受到長江潮位的影響，也和閘門開啟高度、閘門寬度以及內河防洪控制水位等密切相關[16]，以三干河和四干河為例，二者相距2.6 km，引水過程中，長江潮位基本一致，內河防洪控制水位均為3.8 m，閘門開啟高度均為1.2 m，根據引水量分析結果，三干河與四干河引水量比例基本為1.2∶1，三干河和四干河閘門開啟寬度比例也是1.2∶1，表明閘門開啟寬度是影響三干河和四干河引水量差異的關鍵因素。

圖3 內河引水量與天生港潮位響應關系

3.2 平原感潮河網水質改善效果與引水量響應關系

利用構建的模型，計算各片區水質改善效果與引水量的響應關系。結果表明：隨著引水量的增加，各片區水質濃度改善率、類別變化指數和濃度變化指數均呈現逐漸上升的趨勢，表明增大引水量對提高水質改善效果有顯著作用，與崔廣柏等[2]、顧爐華等[7]、童朝鋒等[18]、逄勇等[19]的研究結果相似。由于中部水系水循環體系最為完善，引水時內河側閘前控制水位為4.1 m，朝東圩港水利樞紐和一干河水利樞紐可按照最大引水能力進行引水，且引水過程中，沿線小型內河節制閘均處于關閉狀態，使得長江清潔水源充分到達中部水系末端，對中部水系沿線氨氮超標水體產生改善效果。因此，引水對中部水系氨氮濃度改善效果最強，當引水量達到1.4×107m3時，中部水系氨氮濃度改善率達到50%以上，而后，隨著引水量的增大，氨氮濃度改善率穩定在50%～56%(表2)。東北部水系距離長江較近，引水路線短，且沿線水體水質較好，污染源分布少，因此，引水對東北部水系高錳酸鹽指數改善效果最強，當引水量達到2.64×106m3時，東北部水系高錳酸鹽指數改善率接近40%，而后，隨著引水量的增大，高錳酸鹽指數改善率穩定在40%～43%(表3)。結果表明隨著引水量的逐漸增大，等量引水對河網水質的改善效率逐漸降低，引水水量不再為提高張家港市河網水質的主要因素。

表2 張家港市中部水系水質改善效果與引水量響應關系

表3 張家港市東北部水系水質改善效果與引水量響應關系

類別變化指數和濃度變化指數能夠綜合體現引水前后水質類別和水質濃度的變化情況，反映引水對水質的改善效果，結果表明引水對東北部水系類別變化指數提升效果最好，當引水量達到2.64×106m3時，東北部水系整體提升了1.5個水質類別(表3)；引水對中部水系濃度變化指數改善效果最好，當引水量達到1.4×107m3時，中部水系濃度變化指數整體達到0.6以上(表2)。由于東部水系引水時內河側閘前控制水位為3.7～3.8 m，三干河、四干河未能按照最大引水能力進行引水，且東南部水系位于東部水系末端，引水路線長，缺少清潔水源，因此，引水對東南部水系水質改善效果較弱，氨氮濃度改善率均在15%以內，高錳酸鹽指數改善率均在10%以內，類別變化指數均在0.3以內，濃度變化指數均在0.15以內(表4)。西部水系由于沿江節制閘規模較小，引水通道能力有限，引水對西部水系有一定的改善效果，但改善效果較弱，當引水量達到1.84×106m3時，氨氮濃度改善率和高錳酸鹽指數改善率基本達到35%左右，西部水系整體提升了1個水質類別，濃度變化指數整體達到0.5以上(表5)。

表4 張家港市東南部水系水質改善效果與引水量響應關系

表5 張家港市西部水系水質改善效果與引水量響應關系

3.3 平原感潮河網引水服務面積與引水量響應關系

在平原感潮河網水質改善效果與引水量響應關系研究的基礎上，以引水是否使河道水質達到水功能區水質要求為判斷依據，建立引水服務面積與引水量的響應關系。結果表明：引水服務面積隨著引水量的增大呈現線性增長的趨勢，東部水系引水服務面積隨引水量增長幅度最大，西部水系次之，中部水系最小(圖4)。而單位引水量服務面積隨著引水量的增大呈現對數函數下降趨勢，東部水系單位引水量服務面積最大，西部水系次之，中部水系最小。單位引水量服務面積主要與該區域本底水質濃度、污染源分布情況以及引水路線長短密切相關，由于東北部水系水質本底濃度最低，且污染源分布較少，引水路線短，使得其單位引水量服務面積最高，為0.352～0.891 km2/萬m3；西部水系本底水質濃度較高，但其污染源分布較少，且引水路線短，單位引水量服務面積次之，為0.156～0.303 km2/萬m3；中部水系雖然本底水質濃度較低，但其污染源分布較多，且引水路線長，單位引水量服務面積最小，為0.061～0.146 km2/萬m3。這與其他學者的研究成果相似，陳振濤等[5]對杭州市江干區河網的研究表明，在引水水源水質不變的情況下，增大引水水量，單方水起到的水質改善效果逐漸減小，即引水水量起到的效率降低；張秀菊等[17]對南通市通州區新江海河地區的研究表明，隨著引水量的增加，等量引水的水質改善效果逐漸減弱。

圖4 各片區引水服務面積與引水量響應關系

4 結 論

a. 長江張家港段為不規則半日潮，高潮位具有明顯的半月周期特性，隨著高潮位的升高，內河引水量呈現線性增長趨勢，不同引水口門引水量與長江潮位、內河防洪控制水位、閘門開啟高度以及閘門寬度密切相關。

b. 平原感潮河網水質濃度改善率、類別變化指數和濃度變化指數隨著內河引水量的增加均呈現逐漸上升的趨勢，當中部水系引水量達到1.4×107m3時，氨氮濃度改善率達到50%以上，濃度變化指數整體達到0.6以上；當東北部水系引水量達到2.64×106m3時，高錳酸鹽指數改善率接近40%，且整體提升了1.5個水質類別。由于東部水系內河防洪水位偏低，未能按照最大引水量進行引水，引水對東南部水系水質改善效果較弱，建議逐步實施東部水系防洪加固工程，使內河防洪水位逐步提升至4.0～4.1 m。西部水系由于沿江節制閘規模較小，引水通道能力有限，引水對西部水系水質改善效果較弱，建議逐步實施護漕港節制閘擴建工程和太字圩港水系連通工程，逐步提高西部水系引水能力。

c. 平原感潮河網引水服務面積隨著引水量的增大呈現線性增長趨勢，而單位引水量服務面積呈現對數函數下降趨勢，其中，東北部水系單位引水量服務面積最大，為0.352～0.891 km2/萬m3；西部水系次之，為0.156～0.303 km2/萬m3；中部水系最小，為0.061～0.146 km2/萬m3。