杭州市新取水口水質論證

張舒羽，程杭平

（浙江省水利河口研究院，杭州 310020）

杭州市新取水口水質論證

張舒羽，程杭平

（浙江省水利河口研究院，杭州 310020）

杭州市為抵御咸潮，滿足日益增長的用水需求，計劃將自來水取水口由珊瑚沙上移至里山位置。在規劃期里山取水口的水質能否達到飲用水水質標準是公眾關心的問題。在分析里山河段現狀水質的基礎上，借助Mike21軟件，進行了錢塘江多站位水文及水質的驗證，計算了規劃期里山取水口附近的水質情況，并通過改變上游徑流量、上游來水水質、下游潮汐條件、污染負荷控制情況對里山取水口的水質影響作敏感性分析，從而論證在里山位置設置取水口的合理性。通過計算分析表明，在規劃期的排污條件下，上游枯水流量遭遇下游小潮時，里山河段的水質指標CODMn接近二類水體臨界值，為確保飲用水的質量，需削減城鎮污染負荷的排放。上游徑流量、上游來水水質、下游潮汐條件、污染負荷的改變對里山取水口水質的影響敏感。

水質預測；數學模型；取水口

杭州市自來水供水水源85%取自錢塘江，而每年7－11月錢塘江都要受到大潮的影響，當江水含氯度超過250×10－6時，將不能再作為飲用供水水源［1］。隨著城市的發展，用水量日益增加。目前，作為主要水源的珊瑚沙水庫及其取水口的含氯度指標只有80%的保證率，將不能滿足日益增長的用水需求。為提高供水保證率，計劃將取水口再移至上游12 km處的里山（圖1）。

據調查，錢塘江兩岸均設有排污口，在規劃年（2010年）里山位置的水質能否達到飲用水水質的標準是公眾關心的問題。本文通過收集水文水質資料，借助Mike21水質模型，預測規劃期里山位置的水質情況，并對上游徑流量、上游來水水質、下游潮汐條件、污染負荷的改變等情況作敏感性分析［2］。

圖1 里山取水口平面位置圖Fig.1 The sketch of Lishan W ater Intake

1 里山河段現狀水質

水質監測資料包括為取水工程專門所測的資料和常規測驗站的資料。

從收集的資料來看，里山河段枯水期NH3-N為II－III類，劣于平水期（總體II類）；枯水期總磷總體為II類，優于平水期（II－III類）；CODMn和CODCr在2個水期總體基本相當，為I－II類；平水期總氮指標較高，為V－劣V類。分析原因可能是平水期下泄流量大，水環境稀釋擴散作用優于枯水期，因此平水期NH3-N水質優于枯水期，而平水期總氮和總磷指標較高，可能是徑流量增大，帶入農業面源增加所造成。

2 水質模型介紹

采用丹麥水力學研究所研制的平面二維數值模型Mike21研究工程區域的潮流和污染物對流擴散運動。該模型計算軟件具有算法可靠、計算穩定、界面友好、前后處理方便等突出優點，并在70多個國家得到應用。

潮流模型包括1個連續性方程和2個動量方程?；痉匠虨椋?/p>

式中：ζ為潮位，即水面到某一基準面的距離；p，q分別為x，y方向上的垂線平均流量分量；u，v分別為x，y方向上的垂線平均流速分量；h為水深；g為重力加速度；Ω為柯氏力參數；C′為謝才系數；E為渦動粘性系數；t為時間；C為沿程平均污染物濃度；Dx，Dy分別為x，y方向上的擴散系數；Qs為污水排放流量；Cs為污水排放濃度；As為排放口位置的流體微團面積，在數值計算中取x，y向網格寬度的乘積。

上述方程組的初始條件：

根據上述定解條件，Mike21軟件采用交替方向隱格式（ADI）求解二維淺水潮波方程，方程矩陣采用雙向消除（Double Sweep）算法求解。

3 模型驗證

二維上邊界取在上游窄溪，下邊界選在錢塘江閘口站。模型采用矩形網格，空間步長為50m×100m，時間步長5 s，江道地形采用2004年實測地形。

枯水期驗證采用2004年1月的資料，平水期驗證采用2004年9月的資料?，F狀條件下，模型范圍內污水總量為95萬t／d，CODCr排放量為95 t／d，富陽漁山江段NH3-N排放量為1.73 t／d。計算中，CODCr與CODMn的換算值取2.5。各監測斷面的水質驗證結果見表1。

通過水流、水質的數值模擬和驗證表明，本文所建立的數學模型基本能反映錢塘江河道的實際情況，驗證結果合理可靠，率定的參數可用于設計潮型下的水質預測。

表1 水質驗證結果Table 1 Verification result of water quality

4 規劃期水質計算及敏感性分析

4.1 規劃期水質計算

按照國家規范，水環境評價必須要在上游流量85%～90%設計保證率的水文條件下，其水質仍能滿足水質標準。統計上游富春江電站多年徑流量，保證率在90%的徑流流量為300 m3／s，但近兩年，水資源日趨緊張，小于200 m3／s的流量頻頻出現，針對具體情況，本著慎重的態度，認為對于杭州市的重要水源保護區宜采用保證率95%的流量為妥，其值為200 m3／s。

小潮稀釋水體較少，污染物排放口形成的混合區濃度高，影響較大。因此，取小潮汛作為下邊界是最不利的工況條件。

規劃期，富陽至小沙江段南北岸的負荷之比約為8∶1。江北片污水總量為15.7萬t／d，CODCr排放量為9.5 t／d，較現狀污染負荷增幅為25%；江南片污水總量為76萬t／d，CODCr排放量為75.2 t／d，較現狀污染負荷增幅22%。江北片NH3-N排放量為2.35 t／d，江南片NH3-N排放量為4 t／d。

對于污染物CODMn，模型范圍內所有污染源都有其對應的規劃期負荷，因此可以計算絕對值。對于污染物NH3-N，僅收集到富陽漁山江段的規劃期負荷，缺乏其他眾多污染源的NH3-N負荷，因此NH3-N的計算采用“增量法”，集中反映富陽城鎮污水排放對于杭州市飲用水源保護區的“貢獻值”。

表2為上游枯水流量遭遇下游小潮時在里山取水口水源保護區的CODMn和NH3-N增量全潮平均值?？梢?，規劃期保護區上游斷面（取水口上游4 km）的污染物CODMn值為4.00 mg／L，保護區下游斷面（取水口下游3 km）的CODMn值為3.98 mg／L，總體接近Ⅱ類水體的臨界值；NH3-N全潮濃度增量為0.12 mg／L。

表2 規劃期排污量下里山取水口水源保護區的水質情況Table 2 Water quality of Lishan Water Intake at Water Source Protection Area in programm ing pollution load mg／L

4.2 敏感性分析

里山取水口的水質與上游徑流量、上游來水水質、下游潮汐條件、污染負荷控制情況有一定的關系，但具體的敏感性如何，需要方案計算進行確定。本文擬定了7種工況，以CODMn指標為例，分析不同條件下的水質敏感性，具體的工況設置和計算的成果見表3。

4.2.1 富陽南片污染負荷削減對保護區水質的影響

通過方案1，2，3的比較，分析污染負荷削減對取水口水質的影響。方案1為規劃期負荷的工況，和現狀相比，江北片CODcr負荷增幅25%；江南片CODcr負荷增幅22%，在這種工況下，保護區CODMn指標已經接近Ⅱ類水體臨界值。方案2將富陽南片的負荷維持在現狀水平，北片為規劃排污量，在這種工況下，保護區的CODMn平均值有所下降，基本保持在Ⅱ類水體的標準值內。如果進一步削減南片的污染物，即方案3削減江南片污染負荷為現狀的75%，江北片為規劃排污量，則保護區的CODMn平均值將對應下降，但幅度不是很大。

表3 里山取水口不同位置、不同工況下的水質敏感性分析Table 3 Sensitivity analysis on water quality at Lishan W ater Intake under different conditions and locations mg／L

4.2.2 上游來水水質對保護區水質的影響

通過方案3，4的比較，分析上游來水水質對保護區的影響。隨著上游城鎮和經濟的發展，污染物排放量增加，使上游來水的水質變差，其濃度的改變將對取水口水質產生影響。在富陽南片污染負荷為現狀的75%，北片為規劃排污量的負荷下，方案4上游邊界濃度在方案3的基礎上，增加15%?？芍?，盡管富陽南片的污染負荷在現狀基礎上削減了25%，但只要其上游的污染物沒有加以控制，保護區的CODMn平均值仍滿足不了Ⅱ類水體的水質要求。因此，應嚴格控制上游的污染物排放，其總量應保持現狀的水平。

4.2.3 上游徑流量改變對保護區水質的影響

通過方案1，5，6的比較，分析徑流量改變對保護區水質的影響。在方案1富陽城鎮污染負荷規劃工況的基礎上，考慮將上游徑流量進一步削減為150 m3／s（方案5）和抬升為300 m3／s（方案6）2種工況。通過比較發現，徑流量的改變直接影響保護區的CODMn平均值，在徑流量削減為150 m3／s時，保護區CODMn平均值超過Ⅱ類水體的水質要求；而將上游徑流量提高到300 m3／s，則保護區CODMn平均值有較大幅度的下降，可滿足Ⅱ類水體的水質要求。

4.2.4 下游潮汐條件改變對保護區水質的影響

通過方案1，7的比較，分析下游潮汐條件改變對保護區水質的影響。一般情況下，小潮期污染物稀釋擴散條件差，而大潮期流速大，污染物遷移、擴散作用強于小潮，污染物推移上溯的距離會長一些。因此，在規劃工況的基礎上，考慮大潮的邊界條件（方案7）。通過比較發現，大潮工況對于保護區CODMn平均值的改變量是明顯的，它并沒有將下游排放的污染物推移至影響保護區，而是有效地稀釋了保護區的污染物濃度。

5 結 論

（1）為了確保杭州市取水口的安全，必須設置一級保護區來加以保護，并嚴格遵守國家“飲用水水源保護區污染防治管理規定”中一級保護區有關禁止事項，嚴格加以保護。根據具體情況在里山取水口的上游4 km、下游3 km設定為一級保護區是必須而妥當的。

（2）為了確保杭州市飲用水源保護區受到嚴格保護，必須限制富陽城鎮污染負荷總量的排放，富陽北片可控制在規劃水平；富陽南片必須控制在現狀水平；條件許可的話，應調整富陽城鎮的產業結構，進一步削減污染負荷排入量，提高杭州飲用水源保護區的水質。

（3）上游來水的水質、上游的徑流量和下游的潮汐條件對保護區范圍內的水質影響敏感。對于富陽上游的污染物排放，應嚴格控制其總量的排放，要求污染負荷總量排放必須維持在現狀水平。

［1］ 韓曾萃，戴澤衡，李光炳，等.錢塘江河口治理開發［M］.北京：水利電力出版社，2003：61－62.（HAN Zeng-cui，DAI Ze-heng，LIGuang-bing，et al.Regulation and Exploitation of Qiantang Estuary［M］，Beijing：hydraulic and electric publisher，2003：61－62.（in Chinese））

［2］ 黃少燕.重要水源地水質問題及保護對策研究［J］.亞熱帶水土保持，2006，（12）：65－67.（HUANG Shaoyan.Water quality problem and research on protection measures of main water resources［J］.Soil and Water Conservation of Subtropical Region，2006，（12）：

65－67.

（編輯：周曉雁）

Argumentation on W ater Quality of New W ater Intake in Hangzhou City

ZHANG Shu-yu，CHENG Hang-ping

（Zhejiang Institute of Hydraulics and Estuary，Hangzhou 310020，China）

In order to resist saltwater intrusion and meet the demand of increasing water consumption，Hangzhou City government plans tomove thewater intake from Shanhusha to upstream a site called Lishan.Whether thewater quality of new water intake can reach quality standard has become the focus of public concerns.Firstly，this paper analyzed presentwater quality；secondly，field data of hydrology and water quality of several stations in Qiantang River were validated with the 2-D mathematicalmodel Mike21.Then concentration in programming pollution load was calculated.Furthermore，on the basis of analyzing the influence factors on water quality such as upstream runoff amountand upstream water quality，and downstream tide range and pollution load，the rationality ofwater intake at Lishan sitewas evaluated.From these studies，when upstream low flow meets downstream neap tide，the concentration of water quality index CODMnis close to national second-class seawater standard（GB 3097－1997）.In order to insure water quality，pollution load should be reduced.The factors influencing water quality，such as upstream runoff amount and upstream water quality，and downstream tide range and pollution load，significantly effectwater quality.

water quality prediction；mathematicalmodel；water intake

X522

A

1001－5485（2010）06－0014－04

2009-06-12；

2009-12-07

張舒羽（1976-），女，浙江寧波人，工程師，碩士，主要從事河口海岸工程研究，（電話）0571-86437021（電子信箱）zsyjhy2000＠sina.com。