陳美丹，楊積德，李雋，朱月芳

（1.蘇州市環(huán)境科學(xué)研究所，江蘇蘇州215007；2.澳思環(huán)境科技咨詢有限公司，上海200120；3.蘇州市環(huán)境監(jiān)測(cè)中心，江蘇蘇州215004）

基于源流域水量水質(zhì)模型研究小流域營(yíng)養(yǎng)物管理方案

陳美丹1，楊積德1，李雋2，朱月芳3

（1.蘇州市環(huán)境科學(xué)研究所，江蘇蘇州215007；2.澳思環(huán)境科技咨詢有限公司，上海200120；3.蘇州市環(huán)境監(jiān)測(cè)中心，江蘇蘇州215004）

采用澳大利亞源流域水量水質(zhì)模型模擬東山小流域內(nèi)2001-2010年9種不同土地利用類型產(chǎn)生的降雨徑流及總氮和總磷的污染負(fù)荷，模擬結(jié)果為：流域內(nèi)多年平均徑流量為6150萬m3/a，總氮輸出負(fù)荷為270 t/a，總磷負(fù)荷為22 t/a。模擬結(jié)果表明：水產(chǎn)養(yǎng)殖塘和高地茶果樹是東山地區(qū)主要的營(yíng)養(yǎng)物來源。通過三個(gè)情景方案的模擬，說明當(dāng)?shù)丨h(huán)境管理方案可以有效地削減入湖的營(yíng)養(yǎng)負(fù)荷，其中第二個(gè)情景方案的削減量最大，總氮和總磷負(fù)荷分別削減了18%和25%。

源流域模型；面源；污染負(fù)荷；總氮；總磷

太湖是我國(guó)第三大淡水湖，位于江蘇和浙江兩省交界處，面積達(dá)36900 km2。在過去的20年里，工業(yè)和農(nóng)業(yè)的快速發(fā)展使太湖流域水體富營(yíng)養(yǎng)化成為華南地區(qū)最為嚴(yán)重的環(huán)境問題之一。從20世紀(jì)80年代后期開始，太湖北部梅梁灣就開始頻繁暴發(fā)藍(lán)藻水華，到2006年水華暴發(fā)面積占太湖總面積的一半以上[1～3]；太湖水體水質(zhì)超標(biāo)率也從2001年的65%上升到2006年的95%[4]。近年來，太湖流域工業(yè)點(diǎn)源污染整治力度加強(qiáng)之后，面源污染對(duì)流域水環(huán)境的貢獻(xiàn)及其治理受到越來越多的關(guān)注[5]。許多學(xué)者研究得出總氮的貢獻(xiàn)在34%～52%，總磷的貢獻(xiàn)在17%～54%[6～7]。由此可見，面源污染的影響正在成為繼工業(yè)污染和城市污染之后的最大污染源之一。

研究的目標(biāo)是在太湖流域的一個(gè)示范區(qū)——東山半島構(gòu)建一個(gè)源流域模型（Source IMS），應(yīng)用模型研究流域內(nèi)營(yíng)養(yǎng)物的主要來源；通過模擬不同土地利用類型隨降雨徑流產(chǎn)生的污染負(fù)荷，找出有效改善營(yíng)養(yǎng)物（總氮和總磷）管理的途徑，進(jìn)一步為營(yíng)養(yǎng)物管理提供建議。

1 源流域模型系統(tǒng)（Source IMS）

1.1 模型介紹

在源流域模型系統(tǒng)中，流域水量水質(zhì)模型框架不是一個(gè)模型，而是模型的集合。在這個(gè)集合中，用戶可以選擇對(duì)自己流域內(nèi)的流域水文和污染物生成過程做出最佳描述的徑流、污染物生成、過濾和流內(nèi)模型模塊。污染物生成模型被應(yīng)用于“功能單元”，典型的功能單元是一個(gè)子流域內(nèi)的某種土地利用類型或水文響應(yīng)單元。子流域內(nèi)每個(gè)功能單元的輸出流量和污染物負(fù)荷都在子流域的出口節(jié)點(diǎn)疊加。模型中通過代表河網(wǎng)和流域出口的“節(jié)點(diǎn)”和“鏈接”模擬徑流和污染物的遷移轉(zhuǎn)化過程。除利用源流域模型系統(tǒng)標(biāo)準(zhǔn)配置的模型之外，用戶還可以在模型內(nèi)的任何節(jié)點(diǎn)和鏈接處加入自己的數(shù)學(xué)方程式或表達(dá)式來計(jì)算和記錄所有可能的變量。

1.2 東山概況

東山半島位于太湖的東南部（北緯31°05'，東經(jīng)120°24'），人口約5.3萬人，其中約4.4萬人居住在農(nóng)村，面積為7736 hm2，區(qū)域內(nèi)的土地利用類型面積及所占比例見表1。

1.3 模型構(gòu)建

1.3.1 土地利用類型

依據(jù)不同的土地類型、土地覆蓋、地貌和危害類型，子流域被分成多種功能單元。模型中包括了9種面源土地利用種類：水產(chǎn)養(yǎng)殖、蔬菜種植、低地果樹、高地森林、工業(yè)區(qū)、城鎮(zhèn)區(qū)、農(nóng)村地區(qū)、水面、污水處理廠。

1.3.2 子流域和排水網(wǎng)絡(luò)配置

模型將研究區(qū)域劃分為50個(gè)子流域，由于區(qū)域河道網(wǎng)絡(luò)錯(cuò)綜復(fù)雜，因此將模擬輸出歸納為5個(gè)虛擬節(jié)點(diǎn)（圖1）。節(jié)點(diǎn)1的輸出代表從高地山區(qū)流出半島北側(cè)的降雨徑流和污染物負(fù)荷的集合；節(jié)點(diǎn)3和4代表所有水產(chǎn)養(yǎng)殖塘形成的徑流；節(jié)點(diǎn)2和5代表其余從小河港和運(yùn)河進(jìn)入太湖的徑流。

表1 研究區(qū)域土地利用類型面積及比例

1.3.3 模型運(yùn)行期

模型運(yùn)行時(shí)間為2001-2010年。

1.3.4 降雨徑流生成

本研究項(xiàng)目使用了SIMHYD降雨徑流模型，該模型在濕潤(rùn)溫和氣候條件下模擬效果較好。

1.3.5 小河港、運(yùn)河徑流及泵水過程

徑流經(jīng)過小河港與內(nèi)部交錯(cuò)接連的人工運(yùn)河進(jìn)入太湖。當(dāng)太湖水位漲至3.45 m或更高時(shí)水閘關(guān)閉以防止湖水倒灌，同時(shí)徑流由4個(gè)泵站打入太湖。

1.3.6 氮和磷的生成

對(duì)于污染物質(zhì)的生成，選擇了EMC/DWC模型（事件平均濃度/旱季平均濃度）。該模型的表達(dá)式為：

日負(fù)荷=SF×DWC+QF×EMC

其中：SF——給定日的徑流漫流比例；

QF——給定日的徑流快流比例；

DWC——干旱天氣期間測(cè)得的平均污染物濃度（mg/L）；

EMC——暴雨天氣測(cè)得的流量權(quán)重平均污染物濃度（mg/L）。

2 模型校核與驗(yàn)證

2.1 降雨徑流模型校驗(yàn)

參照當(dāng)?shù)氐奈墨I(xiàn)數(shù)據(jù)對(duì)降雨徑流模型系數(shù)進(jìn)行了校核，所有的模型徑流系數(shù)與文獻(xiàn)值的差異在5%以內(nèi)（表2）。模型中結(jié)合實(shí)際情況，充分考慮了閘泵對(duì)徑流的影響。

表2 研究區(qū)域內(nèi)不同土地利用類型的降雨徑流比例

2.2 營(yíng)養(yǎng)物生成模型

模型輸入?yún)?shù)的總氮和總磷生成率見表3。

表3 分配到各土地利用類型上的基礎(chǔ)EMC/DWC值mg·L-1

3 模擬結(jié)果

在2001-2010年10年運(yùn)行期間，東山多年平均徑流量為6150萬m3/a，總氮輸出負(fù)荷的多年平均值為270 t/a，總磷負(fù)荷多年平均值為22 t/a。具體模型模擬結(jié)果見表4和圖4。

表4 2001-2010年模型模擬結(jié)果

不同土地利用類型徑流量、總氮入湖量、總磷入湖量見圖2。產(chǎn)生年均徑流量最大的是水產(chǎn)養(yǎng)殖塘，約占66%，其次為高地茶果樹。高地茶果樹和水產(chǎn)養(yǎng)殖塘的總氮輸出相近，各占30%～40%。對(duì)于總磷，模型顯示大部分總磷負(fù)荷來自水產(chǎn)養(yǎng)殖塘。

4 規(guī)劃方案及其結(jié)果

4.1 規(guī)劃方案設(shè)計(jì)

根據(jù)當(dāng)?shù)貙?shí)際情況作了3個(gè)規(guī)劃方案，見表5。

4.1.1 改進(jìn)點(diǎn)源管理

模型顯示：將低密度城區(qū)剩余的20%化糞池升級(jí)為二級(jí)處理設(shè)施，同時(shí)將東山污水處理廠納入城市干管，可使進(jìn)入太湖的年均總氮負(fù)荷降低9%，總磷降低13%。此污染削減中的大部分（90%）來自于去除污水處理廠的環(huán)境效應(yīng)。

4.1.2 改進(jìn)面源管理

在蔬菜種植區(qū)和高地茶果樹種植區(qū)采用改進(jìn)的農(nóng)業(yè)管理實(shí)踐，將使這部分總氮和總磷的年均負(fù)荷分別減少9%和2%；而在水產(chǎn)養(yǎng)殖塘區(qū)采用改進(jìn)的管理實(shí)踐，將使這部分總氮和總磷年均負(fù)荷分別減少9%和23%。

圖2 徑流及營(yíng)養(yǎng)物入湖量模擬結(jié)果

表5 情景方案的描述和假設(shè)

因此，在蔬菜區(qū)、高地茶果樹區(qū)和水產(chǎn)養(yǎng)殖塘區(qū)實(shí)施行優(yōu)化管理實(shí)踐，將使每年東山進(jìn)入太湖的總氮量減少18%，總磷量減少25%。

4.1.3 大型濕地建設(shè)

在研究區(qū)域南側(cè)建設(shè)大型濕地，覆蓋面積約7 km2。模型顯示，濕地建設(shè)可以削減整個(gè)小流域13%的總氮輸出量，削減16%的總磷輸出量。

4.2 模擬結(jié)果

三個(gè)規(guī)劃方案的模型結(jié)果見圖3。

圖3 規(guī)劃方案模擬結(jié)果

5 結(jié)論

通過模型模擬流域內(nèi)降雨徑流及營(yíng)養(yǎng)物產(chǎn)生遷移，得出以下結(jié)論：

（1）在2001-2010年期間，東山小流域內(nèi)多年平均徑流量為6150萬m3/a，總氮輸出負(fù)荷的多年平均值為270 t/a，總磷負(fù)荷多年平均值為22 t/a。

（2）區(qū)域內(nèi)水產(chǎn)養(yǎng)殖塘和高地茶果樹種植各自貢獻(xiàn)了30%～40%進(jìn)入太湖的總氮負(fù)荷，45%的總磷來自于水產(chǎn)養(yǎng)殖塘。

（3）三個(gè)情景方案的模擬，說明當(dāng)?shù)丨h(huán)境管理方案可以有效地削減入湖的營(yíng)養(yǎng)負(fù)荷，其中第二個(gè)情景方案的削減量最大，總氮和總磷負(fù)荷分別削減了18%和25%。

[1]范成新，季江，陳荷生.太湖富營(yíng)養(yǎng)化現(xiàn)狀、趨勢(shì)及其綜合整治對(duì)策[J].上海環(huán)境科學(xué)，1997，16（8）：4-7，17.

[2]CHEN YUWEI，F(xiàn)AN CHENGXIN，TEUBNER K,et al.Changes of Nutrients and Phytoplankton Chlorophyll：in a Large Shallow Lake，Taihu，China an 8-year Investigation[J].Hydrobiologia，2003，506/509：273-279.

[3]林澤新.太湖流域水環(huán)境變化及緣由分析[J].湖泊科學(xué)，2002，14（2）：111-116.

[4]虞錫君.構(gòu)建太湖流域水生態(tài)補(bǔ)償機(jī)制探討[J].農(nóng)業(yè)經(jīng)濟(jì)問題，2007（9）：56-59.

[5]王蘇民，竇鴻身.中國(guó)湖泊志[M].北京：科學(xué)出版社，1998：3-21.

[6]張維理，武淑霞，冀宏杰，等.中國(guó)農(nóng)業(yè)面源污染形勢(shì)估計(jì)及控制對(duì)策——2l世紀(jì)初期中國(guó)農(nóng)業(yè)面源污染的形勢(shì)估計(jì)[J].中國(guó)農(nóng)業(yè)科學(xué)，2004，37（7）：1008-1017.

[7]金相燦，葉春，顏昌宙，等.太湖重點(diǎn)污染控制區(qū)綜合治理方案研究[J].環(huán)境科學(xué)研究，1999，12（5）：1-5.

The Application of Source IMS Water Quantity and Quality Model into Investigating Nutrient Management Options in Small Watershed Area

Chen Meidan1,Yang Jide1,Li Jun2,Zhu Yuefang3
（1.Suzhou Environment Science Research Institute,Suzhou Jiangsu 215007,China；2.Earth Systems China Environmental&Consulting Co.,Ltd,Shanghai 200120,China；3.Suzhou Environmental Monitoring Center,Suzhou Jiangsu 215004,China）

Using Source IMS water quantity and quality model to simulate rainfall runoff generation and TN and TP nutrients loads in Dongshan Peninsula though 2001-2010.According to the simulation,the average annual runoff was 61500000 m3/a,total nitrogen output load is 270 t/a,total phosphorus load 22 t/a.The aquiculture pond and highland tea fruit is Dongshan area's main nutrient source.Through three scenarios simulation,the local environmental management scheme can effectively reduce the nutrient load into the lake,the second scenarios'reduction is the most,total nitrogen and total phosphorus load respectively by 18%and 25%.

Source IMS；non-point source；pollution load；TN；TP

TV213.4

A

1008-813X（2012）03-0046-04

10.3969/j.issn.1008-813X.2012.03.014

2012-04-12

陳美丹（1981-），女，江蘇常熟人，畢業(yè)于河海大學(xué)環(huán)境工程專業(yè)，工程師，現(xiàn)從事水環(huán)境研究和環(huán)境科研工作。