基于水量調控的城市湖泊水環境保障措施研究*

孫東旭，丁愛中

(北京師范大學 水科學研究院/教育部地下水污染控制與修復工程研究中心，北京 100875)

水是人類生存發展的基本條件，對于人類生活居住的城市來說，這一條件體現于城市湖泊。一個地區的湖泊存在狀況可以反映一個城市的特征和發展水平。我國有許多城市湖泊，大多數湖泊屬于小型淺水湖泊，可用于防洪航運、娛樂、發電和供水。多數城市湖泊還能作為景觀湖泊提高周圍居民的生活質量，改變周邊的生態系統。

隨著人口增長、工業化和城鎮化進程的推進，人類的生活和生產活動產生了大量的污染物質，城市湖泊極易成為這些污染物的受納體。調查表明，我國大部分城市湖泊、水庫己達到富營養化或超富營養化程度。如杭州西湖、南京玄武湖、江西九江甘棠湖、廣州麓湖、流花湖、武漢墨水湖等。我國大多數湖泊水的總氮(TN)質量濃度超過0.02 mg/L，92%以上水體的總磷(TP)質量濃度超過0.02 mg/L，近50%水體的TP質量濃度為0.2～1.0 mg/L，水質已普遍遭受污染。

針對城市湖泊特點國內外學者對其治理進行了大量的研究。如泰晤士河建立城市污水和廢水的治理系統。美國通過與自然相協調的可持續河流管理理念，建立了完善的濕地生態系統。日本強調用生態工程方法治理河流環境、恢復水質、維護景觀多樣性和生物多樣性。我國對人工生態工程凈化系統和水生植物凈化湖泊方面的研究也取得了巨大進展。

相對于缺水和水質較差的城市，水量調控是城市湖泊水資源維持和改善的重要措施。水量調控治理城市湖泊見效較快，調水時通過改善水動力條件，可以提高水體的自凈能力；水量調控也是湖泊水資源維持的先決條件，不僅可以改善湖泊水質，也可以應對湖泊水體在豐水期和枯水期水位的急劇變化；水量調控還可以對水資源進行合理的分配和調度，避免水資源的浪費。沒有一個合理的水量調控措施做為保障，城市的水資源利用效果將大打折扣。

筆者以北京南海子湖泊為實驗模擬案例，根據水質的監測結果和保護目標，以TP為水質的控制目標，通過建立的水動力學和水質模型，模擬了綜合考慮水質水量的湖泊水調度問題，研究不同水量調控方案對改善城市湖泊水體水質的總體效果和最大收益。為基于水量調控治理城市湖泊水資源狀況這一措施提供依據和建議。

1 北京南海子湖泊概況

南海子郊野公園位于大興區東北部南五環南側、大興新城與亦莊新城之間，全園規劃面積約8 000萬m2。其中，一期工程總占地面積1 524萬m2，湖泊水體面積30萬m2，引水渠總長4.5 km，湖體平均水深1.5 m，總蓄水量約70萬m3。

據預測，湖泊蒸發、滲漏及綠化灌溉需水量約291.67 m3/h，湖泊擬通過引水渠引水入湖416.67 m3/h。由于湖泊尚未建設完成，一期景觀湖仍為目前普遍采用的“盲腸式用水”模式，即湖泊通過水閘等手段獨立于外部水系，不與外部水系一起循環，只有在需要的時候從水源引水進行一定的補水。

目前，南海子湖泊水質保持主要存在以下幾個問題：(1)湖泊水系主要水源為大興區內污水處理廠和再生水廠的再生水，相對公園內景觀環境用水水質具有一定差距；(2)水資源量極度短缺，總供水水量不足，生態換水量和換水頻率不足，導致公園水體流動性差，存在發生水華的潛在風險；(3)湖泊水系為新建，尚未建立自身生態系統，水體自凈能力低，環境容量低；(4)缺乏水體健康監測與評估系統；(5)缺乏水質惡化應急處理系統。

2 模型的建立

在湖泊中選取8個采樣點作為模型建立、模擬、驗證的參考數據，具體位置見圖1，采樣點基本覆蓋了湖泊的整體范圍。

圖1 模型建立的采樣點分布圖

由于南海子湖泊屬于城市淺水湖泊，所以方程需考慮柯氏力、水面風應力及湖底摩擦力，應用Navier-Stokes方程，可得到湖泊的控制方程[1-2]，如式(1)～式(3)所示。

(1)

(2)

(3)

式中：δ為湖面的水位，m；h為水深，m；q為湖泊的流量，m3/h；u為沿x方向的流速分量；v為沿y方向的流速分量；g為重力加速度，m/s；ρ為水體密度，g/m3；f為柯氏力系數；τbx和τby為湖面風應力分量；τwx和τwy為湖底摩擦力分量。

對上述8個點TN、TP的分析結果見圖2。

采樣時間(a) 各點變化

采樣時間(b) 各點平均值變化圖2 南海子湖區采樣點氮磷比隨時間的變化趨勢

根據Stumn研究，氮磷比在小于7時為氮限制，大于7為磷限制。

由圖2的結果可以看出，南海子湖氮磷比大于7，為磷限制湖泊，所以水質模型選取TP作為模擬對象，如式(4)、式(5)所示。

(4)

F(C)=Sp-KTP×C×H

(5)

式中：M、N分別為橫向和縱向的單寬流量，m2/s；Ex為橫向擴散系數；Ey為縱向擴散系數，m2/s；C為水體TP質量濃度，mg/L；S為源匯項，g/(m2·s)；F(C)為生化項；Sp為底泥磷釋放速率，g/(m2·d)；KTP為綜合沉降系數，d-1。

模型求解采用笛卡爾坐標下的淺水方程，網格大小為100 m×100 m，如式6、式7所示[3-5]。

(6)

(7)

式中：h為水深，m；u和v為流速在x和y方向的分量；g為重力加速度，m/s；Δx和Δy為x和y方向的特征長度，m；Δt為時間間距，s。

模型的初始條件為2012年1月11日8個監測點的平均值，模型的南邊界條件是監測點S4和S5的監測數據。

模型各項參數的確定分別為式(8)～式(11)[6]。

f=2×ω×sinφ

(8)

式中：f為柯氏力系數；ω為地球自轉角速度；φ為湖泊所處緯度。

(9)

τωx=CD·ρa·w·wx

(10)

τωy=CD·ρa·w·wy

(11)

式中：τωx和τωy為湖面風應力分量；CD為風應力系數；ρa為空氣密度，kg/m3；w為離湖面10 m處風速，m/s；wx和wy分別為x，y方向的風速，m/s。

根據采樣點水質監測數據和公園建設資料對模型進行了驗證，結果反應了南海子公園湖泊的水體指標的變化特征，所以模型可用于模擬計算湖泊水質和水量的變化結果。

3 水量調控策略模擬分析

3.1 南海子湖泊水質水量現狀

根據南海子湖泊的工程建設條件，假設南海子湖泊一次性引水達到蓄水量為70萬m3，此時水位可達到29.5 m，湖心水深2.5 m，湖面面積為30.81 m2，分析其蓄水量的逐月變化情況。經分析，南海子湖三月份蓄水量較少，可以在三月中旬對湖泊進行引水調控，保證其在春夏兩季水量豐沛。同時通過計算得到湖泊引入水須達到10萬m3，才可保證湖泊在維持水量的基礎上，湖泊污染物濃度進一步下降。

采樣結果表明南海子湖夏季磷濃度較高，其余月份較低。全年平均值為0.143 mg/L，各采樣點的水質優劣不同。表1為2012年～2013年各個采樣點TP的平均監測值。由表1數據可以看出湖東的磷質量濃度低于湖西。

表1 2012年～2013年各采樣點TP質量濃度平均值

3.2 水量調控方案設計

為了改善城市湖泊水質環境，降低污染物在湖泊中的總體濃度，筆者提出3種調控策略，選擇在枯水季3月為湖泊進行引水和入湖徑流作為模型的邊界條件，以TP為狀態變量，進行模擬計算。

(1) 湖西入水口引水：從2012年3月11日到2012年4月11日，對湖西入水口進行入水，進水量為0.4 m3/s，引入水經湖東入水口流入湖西，從湖西出水口排出。

(2) 湖東入水口引水：從2012年3月11日到2012年4月11日，對湖東入水口進行引水，進水量為0.4 m3/s，流入湖東后進入湖西，同時從湖東出水口流出。

(3) 對湖西入水口和湖東入水口分別進行引水：從2012年3月11日到2012年4月11日，對湖西入水口進行引水，進水量為0.2 m3/s，在2012年3月11日到2012年4月11日對湖東進行引水，進水量為0.2 m3/s，同時通過湖西的自流引水改善湖泊的整體水動力和水質條件。

3.3 方案模擬結果分析

3.3.1 不同調控方案模擬分析

通過模型模擬發現，第1種方案容易造成湖西污染物流入湖東，湖東水質條件得不到改善，甚至在采樣點S7和S8處水質條件進一步變差，使湖西的污染遷移到湖東區域。第2種方案由于受南海子湖泊水動力條件限制，湖東的引入水不能進入湖西，無法改善湖西的水質環境。造成這2種方案無法有效地治理湖泊污染的主要原因是受到湖泊的水動力條件限制，同時湖西和湖東的污染狀況并不相同，污染物容易通過湖泊的水循環發生遷移，容易使水質較好的區域面臨水質惡化的風險。

在第3種引水方案下，模型模擬得到各采樣點的水質指標TP含量的變化過程，見圖3～圖8。

圖3 2012年采樣點S1處TP含量變化

圖4 2012年采樣點S3處TP含量變化

由圖3和圖4可以看出對湖西入水口處進行引水方案可以較快地改善湖西水質，有效地降低 湖西污染物向湖東遷移的風險。

圖5 2012年采樣點S4處TP含量變化

圖6 2012年采樣點S5處TP含量變化

圖7 2012年采樣點S6處TP含量變化

圖8 2012年采樣點S7處TP含量變化

由圖5～圖7可以看出通過湖東入水口引水和湖西的自流引水，南海子公園湖東各區域水質得到有效改善。通過圖8可以看出，采樣點S7處TP質量濃度并未發生大幅度的升高變化，說明湖西區域污染物和湖東其它區域污染并未在采樣點S7區域滯留，大幅度降低了污染物遷移到湖泊其它區域的風險。

3.3.2 南海子湖東不同引水流量模擬分析

對于水量調控方案三，模擬設計3種不同的引水規模，分別為0.2、0.4、0.6 m3/s，對這3種流量進行模擬計算，連續引水時間為20 d，對各個采樣點區域的TN和TP以及COD平均濃度進行數學統計，計算引水后相對于引水前的各個水質條件的改善效果，如表2～表4所示。

由表2～表4可見，如果按照0.2 m3/s引水流量對湖東入水口進行引水，采樣點S7處的水質容易發生進一步惡化，其主要原因是在南海子湖泊水動力條件的作用下，湖東東北部區域容易使污染物滯留，其它湖泊區域污染物發生遷移，影響湖泊個別區域的水質改善效果，甚至使水質相對較好的區域水質變差。而湖西0.1 m3/s的出水流量不足以在引水時間內使污染物都排出進入湖東。因此0.2 m3/s的引水流量方案不能完全有效地改善南海子湖泊各個區域的水質狀況，引水規模略顯不足。

表2 3種流量分別對各采樣點TN影響效果統計表

表3 3種流量分別對各采樣點TP影響效果統計表

表4 3種流量分別對各采樣點COD影響效果統計表

4 結 論

水量調控這一措施可以快速有效地改善湖泊的水質條件，使城市湖泊水體中污染物濃度大幅度降低，但是采用水量調控治理湖泊還需要考慮多種不同因素。

(1) 針對城市湖泊水動力條件，需要設計一個合理的引水調度方案，對湖泊水質進行補水和換水，如果引水調控方案不合理，不僅不能有效地降低城市湖泊的污染物濃度，而且使城市湖泊中的主要污染物滯留在水體中。

(2) 引水規模對湖泊的治理效果有著極大影響。隨著引水流量的增大，湖泊的主要水質指標如COD、TN和TP的改善結果不斷增強，同時湖泊的改善增值率則不斷下降。因此，在采用水量調控這一措施時，筆者不僅要考慮湖泊的治理改善效果，同時還要考慮治理的投資運行成本，隨著引水規模的不斷增大，投資成本也隨著不斷增加，但是改善效果卻在不斷降低，這會造成投資效益比的明顯下降。因此根據不同城市湖泊的工程數據，需要選取最為合適的引水流量，同時需要綜合考慮湖泊的環境、經濟和技術的相關影響因素。

(3) 雖然合理的水量調控措施可以快速有效地治理城市湖泊污染，但是同時還存在許多風險與局限性。首先，城市湖泊的水動力條件決定了再引水調控過程中，污染物容易發生遷移，使污染物進入湖泊中水質相對較好的區域，造成了水質相對較好的區域的水質惡化，所以選取合適的引水流量和引水周期，使湖泊的主要污染物排出湖泊水體是引水方案的重中之重。其次，引水調控方法治理城市湖泊的效果容易在幾年后發生反彈，所以在引水調控湖泊水體之后，還要控制入湖污染物進入水體。

在引水調控過程中，改善了城市湖泊的水動力條件，通過控制入湖污染，能大幅度改善湖泊水質，避免在非引水期內，湖泊各區域水質產生不同程度的反彈。在控制污染物進入水體、引水調控、生態修復等多種治理措施的綜合作用下，湖泊水質才能逐年得到有效改善。

引水調控方案只適用于相對較小的城市湖泊，對于相對較大的城市湖泊引水換水治理措施不可行，這不僅造成了水資源的浪費，還容易使水體水質進一步惡化。因此在治理城市湖泊污染的過程中，關鍵在于如何選取合理的治理手段，以最小的經濟投入使水體達到最理想的狀態。

參 考 文 獻：

[1] Jin Kang-Ren,Hamrick John H,Todd T s Dale.Application of three dimensional hydrodynamic model for lake okeechobee[J].Journal of Hydraulic Engineering,2000,126(10):758-771.

[2] 馮民權，鄭邦民，周孝德.河流及水庫流場與水質的數值模擬[M].北京：科學出版社，2007.

[3] 劉昌明，陳志愷．中國水資源現狀評價和供需發展趨勢分析[M]．北京：中國水利水電出版社，2001.

[4] 張荔，王曉昌．小流域水量水質綜合模型模擬研究[J]．干旱區資源與環境，2008，22(2)：10-13.

[5] 游進軍，薛小妮，牛存穩．水量水質聯合調控思路與研究進展[J]．水利水電技術，2010，41(11)：7-9.

[6] 董增川，卞戈亞，王船海，等．基于數值模擬的區域水量水質聯合調度研究[J]．水科學進展，2009，20(2)：184-189．