太湖錫澄地區TP和CODMn水環境容量的影響因素

王思遠, 許有鵬, 王丹青, 高 斌, 王 強, 溫海燕

(1.南京大學 地理與海洋科學學院, 南京 210023; 2.無錫市城市防洪工程管理處, 江蘇 無錫 214000)

隨著社會經濟的快速發展和城市化進程的不斷推進，太湖流域生態環境急劇惡化，尤其是水域污染和富營養化日趨嚴重，總氮、總磷、生化需氧量(BOD)、化學需氧量(COD)等多項污染物指標超標，水質的急劇惡化逐漸成為經濟發展的制約因素[1-2]。為此各地政府開展了大量的河流污染治理工作，雖然局部水域污染狀況得到一定程度改善，但是水環境惡化趨勢并未得到有效遏制[3]。水環境容量是水體中污染物總量控制的重要依據，尤其是對于高度城市化平原河網地區，其影響因素復雜多變，因此充分利用污染物在水體中的遷移、轉化等機理，研究和控制水環境容量成為了國內外學者的熱點話題[4]。

水環境容量一般指水體在滿足特定功能條件下所能容納污染物的最大負荷量，其大小與水體特征、水質目標及污染物特性有關，因此，又稱水域負荷量或納污能力[5]。水環境容量的計算及其在各功能區間的分配，是水污染總量控制的基礎和核心[6]。目前，國內外對河流水環境容量計算和研究主要是通過公式法[7]、模型試錯法[8]、概率稀釋模型[9]和盲數理論[10]等方法。如李昆等人對豐水期不同程度富營養化的洪湖全湖，采用GIS 空間插值技術，分析豐水期各污染物的變化趨勢[11]；為改善山區水環境狀況，陳玥等人對臨海市COD和NH3—N的水環境容量進行計算分析，并提出污染物排放總量控制方案[12]；馮利忠等采用MIKE三維水動力和水質耦合模型，對“引黃入呼”取水口水環境容量進行動態分析[13]。但是，這些研究多是基于實測流量作為設計水文條件進行計算，并且大多研究未考慮引排水對水環境容量的影響。對于太湖平原河網地區，各個河流流速緩慢，流向和流態呈現出隨機性變化，河流流向不定，相較于流量控制，基于水位作為控制條件，計算設計流量，更適于太湖平原河網水環境容量的計算。

太湖流域平原河網區是我國經濟最發達和最活躍的地區之一，但是城市化、工業化的快速推進也使得區域經濟發展和水環境保護的矛盾十分突出。在此背景下，以流域典型區域錫澄地區為例，從多角度探討分析平原河網區水環境容量的影響因素對區域經濟可持續發展具有重要意義。本研究選取CODMn和TP兩個典型污染因子開展一系列研究工作，旨在從年內四季、典型年、引水量三方面揭示錫澄地區骨干河流的典型因子的水環境容量變化規律及影響因素，以期為城市化下平原河網地區河網水質治理，河道水環境容量的控制與分配提供科學的參考依據。

1 材料與方法

1.1 研究區概況

研究區位于太湖流域平原河網地區，覆蓋錫澄地區的北塘區、南長區、惠山區、崇安區、濱湖區、錫山區、江陰市7個行政區，總面積2 625.31 km2，其中水域面積占40%左右，擁有各級河道2 000多條，以具有引排作用的縱向河道和調節作用的橫向河道為框架，構成了研究區骨干河網(圖1)。該地區氣候屬亞熱帶濕潤季風氣候，四季分明，降水豐沛，根據近30 a的統計資料，年均降水量達1 121.7 mm。區域經濟較為發達，人民生活水平較高，據2017年統計資料，人均生產總值達到16.07萬元，但水質狀況不容樂觀，在2014年，部分水質站點的污染物的濃度已低于Ⅴ類水質，選擇易導致河流富營養化的磷元素以及表示水體中有機物污染的高錳酸鹽指數，開展對錫澄地區骨干河道水環境容量的研究。

圖1 研究區區位圖

1.2 數據來源

本研究涉及資料包括太湖流域錫澄地區的水系、水文、氣象資料以及水質數據。水系資料包括基于矢量地形圖獲取的無錫市2010s 水系圖，以及根據無錫市河湖普查獲取的8條骨干河道的斷面與河長信息。水文、氣象資料包括水位、雨量、引水量數據，水位數據包含研究區內的白芍站、無錫站、甘露站、陳墅站、青旸站5個水位站2005—2014年的水位資料，雨量數據是通過收集具有代表性的雨量站點的2005—2014年雨量資料獲得的，引排水數據主要根據《長江流域水文資料》2010—2014年年鑒綜合整理分析獲得，進而確定研究區骨干河流的引水量。水質資料為無錫市主要河道上的8個水質監測斷面2005—2014年的月實測數據，根據水質監測點的水質數據，選取對河流水質影響較為突出的常規監測指標總磷(TP)和高錳酸鹽指數(CODMn)分析其水環境容量。

1.3 研究方法

由于本研究基于整條河流作為控制單元，是從較為宏觀的角度對研究區的8條河流的水環境容量進行研究。完全混合計算法是指污染物進入到河流水體后，在斷面上污染物完全均勻混合，較好的符合本次研究的需求。在采用此方法的基礎上，考慮河道寬度等河道特征以及不同污染物的降解系數，探討該區域TP和CODMn水環境容量及其影響因素。主要從以下3個方面展開：(1) 考慮研究區功能區劃分的復雜性，確定8條骨干河流作為研究對象；(2) 以近十年最枯月平均水位(90%保證率)作為設計水位，計算2010—2014年的水環境容量，分別從四季、引水量兩方面考慮水環境容量的影響因素。(3) 根據近十年雨量觀測數據，確定研究區枯、平和豐3種年型，計算對應年型的水環境容量。

考慮平原河網分布特征，參考研究區的水環境功能區劃，將天然河道進行概化，確定研究區河網。將水體水質看作完全均勻混合類型，基于水位控制，完成對研究區河網水環境容量數學計算模型的構建。在2020年水質目標控制下，基于設計水文條件，利用研究區河網水環境容量數學模型，對研究區域最小空間單元(河段)和最小時間單元(月)的水環境容量進行計算，公式如下：

(1)

式中:αji為不均勻系數，范圍為(0,1]；wji為控制單元水環境容量(t/月)。

wji=Q0i(Csji-C0ji)+K(Sji·li)Csji

(2)

式中：Q0i為i河段的設計流量(m3/s)；Csji為不同控制單元的水質標準下的水質濃度(mg/L)；C0ji為不同河斷面的每個月的水質濃度(mg/L)；K為污染物綜合降解系數(d-1)；Sji為河道斷面面積(m3)；li為河長(km)。

1.4 水環境容量數學模型參數率定

(1) 近十年最枯月的平均流量。通過謝才—曼寧公式求得流速分布：

(3)

式中：vi為平均流速(m/s)；n為河道糙率；Hi為水深(m)；Ji為水面比降。

(2) 水質目標。依據江蘇省地表水功能區劃所規定的水質目標，計算研究區骨干河道的水環境容量。確定8條河流的水質目標等級，進而得到每種污染物的水質目標濃度。

(3) 污染物的降解系數。污染物的降解系數是計算水環境容量的基礎參數。本研究根據江蘇省太湖流域污染物降解系數的研究成果，確定研究區主要河段的TP降解系數的變化范圍為0.05～0.1，CODMn的降解系數的變化范圍為0.1～0.29[15]。各河流的污染物降解系數是不同的，但相差甚微，故取其在各河段均值，以從宏觀角度估算河流水環境容量。

(4) 不均勻系數。考慮到水資源時空分布的不均勻性，對每條河流求得的水環境容量進行修訂。影響研究區水環境容量計算結果的主要因素是河寬、流量以及水深。根據太湖流域河流的調查結果表明，河寬基本上可以表征河流的水文特征，依據河寬和不均勻系數的關系，確定不同河道的不均勻系數[4]。

2 結果與分析

2.1 季節性差異分析

根據我國在氣候上劃分四季的標準，一般以3—5月為春季、6—8月為夏季、9至11月為秋季、12月到次年2月為冬季，以此標準劃分研究區的四季。根據2010—2014年TP和CODMn的水環境容量計算結果，得到兩種污染物四季水環境容量年均值，分別用WTP和WCODMn表示，單位為t。

TP水環境容量的季節性差異較為明顯，呈現出秋季較高，冬季較低的特征。春季TP水環境容量范圍為-1.66～3.42 t，均值為0.37 t；夏季的變化范圍為-1.33～3.38 t，均值為0.82 t；秋季的變化范圍為0.21～3.43 t，均值最大為1.11 t；冬季的變化范圍為-2.7～3.33 t，均值最小為0.24 t，并且TP四季水環境容量顯示為負值的河流百分比分別為37.5%(春)，37.5%(夏)，0%(秋)，50%(冬)，這充分表明大部分河流在秋季對磷的納污能力較強，冬季有一半左右的河流水環境容量為負值，表明TP排放量已超出水體容納能力，水環境質量較差。每年5—10月為太湖流域平原河網區的汛期，根據近十年的雨量資料，研究區汛期、非汛期的年平均降雨量分別為740.73 mm，352.35 mm，降雨量對水位有一定的貢獻，而夏、秋季節正處于汛期，長歷時雨量的積累有助于河道水位上漲，使河流水量有增大趨勢，進而間接影響了TP的水環境容量的大小，因此降雨量的差異性對水環境容量具有一定的影響。另外，農業生產中水田與旱地主要污染來源于種植過程中化肥和農藥的施用，造成磷元素的污染較為嚴重。有研究表明氮磷肥料的使用和水產養殖的農業方式對環境的磷污染貢獻很大[16-17]。水域面積較大的太湖平原河網地區，種植水稻、水產養殖就是主要的農業方式，餌料是磷污染產生的原因之一，非點源污染主要以降雨為載體遷移到河流中，影響磷的濃度，而夏、秋季正是水產養殖業進行的重要時期，這就導致磷濃度有一定季節性差異，但結果顯示秋季磷的水環境容量相對較大，水體中磷濃度相對較低，這可能是因為降雨對水環境容量的影響程度可能要大于水產養殖的影響(表1)。

表1 無錫骨干河流TP水環境容量的季節性分布t

CODMn水環境容量的季節性差異不顯著，根據四季CODMn水環境容量的變化范圍，春、夏、秋、冬的變化范圍分別為17.49～151.16 t，17.72～163.75 t，19.60～155.15 t，23.19～165.01 t，呈現出春季偏低，夏、秋季節偏高的特征。在太陽輻射強度相對較高的季節(夏、秋季節)，水體中有機物更易發生降解，使其濃度降低[18]，并且夏、秋季節降水較多，水位相對升高，水體體積也會隨之變大，使得相應水體的水環境容量數值增大。在空間上，CODMn水環境容量存在一定的季節性特征，京杭運河四季的水環境容量均處于中低值范圍，而伯瀆港的水環境容量四季均高于85 t，即處于中高值偏上。京杭運河是江浙地區的水利重要樞紐，流經鎮江、常州、無錫等經濟較為發達的城市，京杭運河的水質受到的工、農業生產以及日常生活污水的影響要遠高于伯瀆港這樣的小型河道[19]，導致京杭運河的有機物的濃度也相對較高，故京杭運河的CODMn水環境容量要低于伯瀆港(表2)。

表2 無錫骨干河流CODMn水環境容量的季節性分布 t

2.2 引水量對水環境容量影響分析

在2007年太湖流域藍藻爆發后，“引江濟太”調水試驗工程開始啟動，以期改善太湖平原河網的水質。定波閘位于錫澄運河與長江交接處，在“引江濟太”工程中發揮著重要作用，2010—2014年，年平均引水量為1.05×107m3，最大年引水量為1.64×107m3。以定波閘為例，分析引水量對錫澄運河水環境容量的影響。由于降雨量的不同，年內各月引水量存在一定的差異，為了更加準確的呈現引水量與水環境容量的關系，對2010—2014年引水量與水環境容量的五年月均值進行研究，兩種污染物的水環境容量與引水量的對應關系如圖2所示。

圖2 錫澄運河水環境容量與引水量的關系

通過觀察分析，對于太湖平原河網區典型河道錫澄運河，引水量越大，污染物的水環境容量也越大，就TP而言，由于部分河流中磷的濃度較高，超出水質目標，故水環境容量表現為負值，表明磷的排放量大于河流對磷的容納能力。圖2中顯示錫澄運河的TP水環境容量為負值的月份，對應的該月份的引水量也相對較少，但TP的水環境容量呈現出正值的月份，引水量也較大。圖2中錫澄運河的CODMn水環境容量均是正值，引水量較大的月份，CODMn水環境容量也相對較大。以上現象表明閘壩引水較大時，在一定程度上可以提高污染物的水環境容量，間接說明“引江濟太”工程的啟動對于改善平原河網的水質有一定作用[20]。引水調控措施不僅可以增加水體的稀釋能力，還能促進水體流動和交換，進而增強水體的自凈能力，從而增大水體對污染物的容納能力，近年實踐表明其是一種行之有效的改善水環境的措施[21]，張聃等人也曾就引配水對京杭運河水質的影響進行研究，得出引配水對河流水質具有一定的改善作用[22]。

2.3 不同典型年下水環境容量差異分析

根據2005—2014年雨量資料，選擇2010年為枯水年、2012年為平水年、2013年為豐水年，利用3種年型的水質資料，分別計算TP和CODMn典型年的水環境容量，結果見圖3。

通過分析八條河流兩種污染物典型年的水環境容量，得出TP和CODMn的水環境容量在豐水年均高于枯水年。對于TP而言，3種典型年均有河流的水環境容量為負值的情況，以枯水年較為顯著，這些河流包括北興塘河、古運河、錫澄運河等，且主要分布于無錫市市區，由此說明磷含量的超標不僅與農業污染有關，與生活污水和工業廢水的大量排放也有很大關系。CODMn代表河流中有機物的污染狀況，對于不同的典型年，CODMn的水環境容量呈現出相同的規律，九里河和京杭運河均較小，北興塘河、伯瀆港和古運河均較大。這種現象可能與河流所處的位置與功能區劃有關。

圖3 錫澄地區主要河流各典型年TP和CODMn水環境容量分布

3 結 論

(1) TP四季水環境容量顯示為負值的河流百分比分別為37.5%(春)，37.5%(夏)，0%(秋)，50%(冬)，這充分表明大部分河流在秋季對磷的納污能力較強，夏季次之，春、冬季均有一半左右的河流水環境容量為負值，表明TP排放量已超出水體納污能力，水環境質量較差。CODMn水環境容量的季節性差異不太明顯，相對而言，夏、秋季節較好，且無錫市夏季水體對其的容納能力最強。

(2) 兩種污染物的水環境容量與引水量的對應關系較為明顯。污染物的水環境容量大的月份，引水量也較大。且引水量較大的月份大都處于非汛期，降雨量較少的月份，水體中污染物濃度較高，需要靠引水，增加河流的水量，污染物濃度被稀釋，以達到降低污染物濃度的目的，改善了平原河網的水質，間接提高了水體對污染物的容納能力。

(3) 通過分析八條河流兩種污染物典型年的水環境容量，總的規律是豐水年的水環境容量較大，水體所能容納的污染物的負荷量較多。研究的大部分河流位于市中心位置，其水功能區劃大部分是工業用水，表明所研究主要河流的主導功能是用于工業發展，這與市區工業化的快速發展是分不開的，進而導致的水體污染也較為嚴重，相關部門需采取一定的措施，對工廠排放的廢水廢氣處理達標再排放，從源頭上防止污染趨勢進一步惡化。

(4) 基于本研究已有數據，僅研究了總磷和高錳酸鹽指數的水環境容量，對于濃度已經嚴重超標的總氮以及氨氮的水環境容量，其影響機制還有待進行較為深入的研究分析工作。