基于SWAT模型的黑水河流域非點源污染控制方案研究

徐凌云

(上海勘測設計研究院有限公司,上海 200335)

隨著工業點源污染控制水平的提高，非點源污染成為了流域水環境污染的主要原因，而流域內的土地利用是非點源污染的主要影響驅動因子。土地利用通過各種直接或間接參與生態環境系統的物質能力轉換與流動，對生態環境質量及其組分如水、大氣等產生巨大的影響。不適當的土地利用方式和農田管理模式會導致土壤侵蝕和過量的N、P隨地表徑流流失[1-2]，從而形成流域內的大面積非點源污染，引發水體水生態環境質量的惡化。非點源污染負荷的防治是提高流域水生態環境質量的重要舉措之一，目前對非點源污染負荷的研究應用最廣泛的是美國農業部(USDA)農業研究局(ARS)開發的SWAT模型[3-4]，該模型在非點源污染負荷的測算與時空分布[5-7]、污染防控或管理措施的效果評估[8-11]、污染源的識別與控制[12-14]、關鍵污染區域與污染時間的識別[15-17]等多方面均得到了廣泛的應用。本文利用SWAT模型建立四川省黑水河流域非點源污染負荷預測模型，模擬流域泥沙及污染物空間分布及不同控制措施情景下流域非點源污染負荷的變化情況，探討流域不同土地利用方式與非點源污染之間的定量關系，為流域土地利用規劃和非點源污染負荷的控制提供實際解決方案。

1 研究區域概況

黑水河是金沙江左岸一級支流，位于四川省涼山彝族自治州境內，介于東經102°20′～102°53′、北緯28°7′～26°48′，發源于昭覺縣瑪果梁子，自北向南流經昭覺、普格、寧南3縣，于寧南縣東南部葫蘆口注入金沙江(圖1)。流域總面積3 591 km2，干流全長173.0 km，天然落差1 931 m，平均比降11.05‰，河口處多年平均流量為80 m3/s，徑流量25.25億m3。目前黑水河流域土地利用方式主要有林地、坡耕地、耕地、建設用地、裸地、水域及水利設施用地等，其中坡耕地和耕地主要沿黑水河、則木河、西洛河及其兩側支溝的河谷地帶分布。根據對黑水河流域的初步調查分析，威脅流域內水文、水質及自然生態系統的主要因素包括地表徑流污染、采石挖沙、農田施肥及水產畜禽養殖等。

2 模型建立及率定驗證

2.1 基礎資料

a)DEM數據。采用流域30 m×30 m的DEM數據。DEM圖中可以看出，黑水河流域北高南低，地貌多以山嶺和河谷相間，河谷到分水嶺高差較大，在700～4 000 m，沿河支溝發育。

b)土地利用。根據黑水河流域1∶250 000的土地利用圖，結合SWAT模型的模擬要求，將土地利用類型重新劃分為5類。其中林地面積最大，占總面積的52.26%，其次是草地和旱地，分別占31.25%和13.81%，水田占比2.47%，裸地占比0.21%，為了研究不同控制措施的效果，忽略占比較小的采石挖沙、居民點、工礦養殖等土地利用類型。

c)土壤數據。土壤數據采用第二次全國土地調查南京土壤所提供的1∶100萬土壤數據，包括土壤厚度、土壤顆粒組成、持水率等基本信息，土壤類型采用收集到的1∶250 000土壤類型空間分布圖。土壤屬性部分在HWSD數據庫里查詢得到，部分采用模型默認值，另一部分利用SPAW軟件計算或者公式法計算得到屬性值。

d)植被覆蓋。SWAT模型使用簡單的植被生長模型來模擬所有的陸地覆蓋類型，模型能夠區分一年生植物和多年生植物，一年生植物從種植日期開始到收獲日期，或直接積累的熱量單元等于植物的潛在熱量單元；多年生植物全年維持其根系系統，在冬季月份中進行冬眠；當日均空氣溫度超過最小即基準溫度時，重新生長。植物生長模型用來評估水分和營養物質從根區的遷移、蒸發以及生物產量。SWAT模型中有關土地利用和植被覆蓋的數據通過文件crop.dat進行存儲和計算。

e)氣象數據庫。SWAT模型需要的氣象數據主要包括流域的日降水量、最高/最低溫度、太陽輻射、風速和相對濕度，這些數據可以是統計數據，也可以通過SWAT模型的天氣模擬程序(天氣發生器)生成，或者是統計和模擬生產數據的結合。本次模擬氣象數據來自會理、西昌、昭覺3個氣象站連續8年的資料，包括平均風速、平均氣溫、日最高氣溫、日最低氣溫、平均相對濕度、太陽輻射以及降水量數據等。

f)水文水質數據。水文水質數據主要通過實地監測、資料收集或購買獲得，用于模型校準驗證的水文、泥沙數據為寧南水文站2008—2013年的實測月徑流、年均輸沙量資料，水質數據采用黑水河河口2008—2013年豐、平、枯水期的實測水質。不同控制措施的研究基于2016年土地利用類型，模型輸入條件采用2016年各子流域的污染源月均日負荷、各入流點的日流量以及流域內主要水庫的相關數據。

2.2 模型建立

2.2.1子流域劃分原則

①以黑水河一級支流作為1個集水區；②規模較大的一級支流(如西洛河、則木河等)，以黑水河二級支流作為1個集水區；③為避免產生過小的集水區或明顯錯誤的破碎集水區，支流匯口局部區域并入干流集水區。根據以上劃分原則，黑水河流域共劃分25個流域單元(圖1)。

圖1 黑水河流域范圍及子流域劃分

2.2.2水文響應單元劃分

本次研究采用每個子流域多個HRUs的劃分方法。根據土地利用和土壤類型的組合和分布特征，確定土地利用和土壤類型的面積閾值，土地利用的閾值為10%，土壤類型的閾值為20%。對土地利用和土壤類型圖進行疊加，在整個子流域內生成138個水文響應單元。

2.2.3農業管理模擬

根據研究區域的農產品種植結構和農業施肥的調查資料，選用流域內主要耕作方式：旱地為冬小麥、油菜/夏玉米、馬鈴薯輪作，水田模擬為水稻。具體的施肥量和種植季節見表1。四川農作物農藥使用強度按15.7 kg/hm2計算。

表1 區域主要作物種植情況及農藥化肥使用量

2.3 模型率定和驗證

2.3.1模型的參數設置

將2008—2010年設置為率定期、2011—2013年設置為驗證期；參數在初始賦值的基礎上通過自動和手動調參以適應流域的實際情況，模型參數取值見表2。

表2 模擬率定參數選取及取值

2.3.2模型率定結果

在SWAT模型中，可以選擇相關系數R2、相對誤差Re和Nash-Suttcliffe系數Ens來評價模型的適用性[18]。相關系數R2=1表示非常吻合，當R2<1時，其值越小反映出數據吻合程度越低。Nash-Suttcliffe系數Ens=1表示非常吻合，若Ens為負值，說明模型模擬平均值比直接使用實際平均值的可信度更低。一般Ens>0.5，R2大于0.6，Re小于15%時表示精度可以滿足模擬要求。率定及驗證情況見表3。

表3 寧南水文站月徑流、泥沙、總氮、總磷率定驗證結果

由表3可知，率定期的徑流、泥沙負荷和氮磷負荷模擬結果除了Re略高外，其余指標均在合理范圍內，驗證期的模擬效果較好，3個指標均在合理范圍內，總體來說，率定期和驗證期的模擬結果相差不大，均可以接受。本研究采用該模型模擬的模擬結果較為可信。

3 結果與討論

3.1 流域產水量及污染負荷空間分布

利用校準驗證后的黑水河流域模型，模擬計算基于2016年土地利用類型的各子流域的產水量及污染負荷(圖2)。

模擬結果顯示，黑水河流域產水量較小的區域主要集中在黑水河流域東北部地區，包括普格縣東部和昭覺縣南部地區，黑水河流域其他地區產水量均較大；黑水河流域產沙強度與產水量大小分布基本一致，同時林地分布面積較大的區域產沙強度小，靠近河流且耕地分布較廣的區域產沙強度較大；從黑水河流域污染負荷空間分布來看，有機氮和有機磷污染負荷強度空間分布特征基本相同，有機氮和有機磷污染負荷強度最大的區域主要集中在流域西北部的和流域西南部。

a)月均產流量

3.2 不同土地利用類型的非點源污染負荷

以月為時間步長，計算黑水河流域基于2016年土地利用類型的非點源污染負荷總量和單位面積負荷量，結果見表4、5。

3.2.1流域總污染負荷

由表4可知，泥沙和氮磷非點源污染負荷的流失與土地利用類型關系密切，旱地產生的非點源污染負荷量最大，各項污染負荷占比均超過50%，其次為草地，再次為林地，最后為水田，由于水田面積占比較小，其污染負荷的貢獻量也最小。綜上可以看出，不同土地利用類型中非點源污染負荷比例不一，主要受人類農業生產活動的影響。

表4 不同土地利用類型非點源污染負荷

表5 不同土地利用類型單位面積非點源污染負荷

3.2.2單位面積污染負荷

不同土地利用類型下，流域單位面積泥沙負荷從高到低依次為旱地、草地、水田和林地；單位面積有機氮和有機磷負荷從高到低依次為旱地、水田、草地和林地。由此可知，在黑水河流域旱地、水田是非點源污染負荷的主要產出方式，林地具有水源涵養功能，在水土保持和污染物截留等方面發揮著巨大的作用。

3.3 不同控制措施及效果分析

將單位面積污染物負荷產出高的子流域作為流域的關鍵源區，根據黑水河流域污染負荷空間分布結果，流域的非點源污染防治的關鍵源區主要在流域西北部和西南部，如子流域11、22、24，這些子流域的耕地分布較廣且靠近河流，需要重點關注。考慮將關鍵源區內單位面積非點源污染負荷輸出最大的旱地改為水田，將零星的裸地改為泥沙輸出負荷最小的林地以及退耕還林、設置緩沖帶等措施，模擬實施不同控制措施后流域非點源污染負荷輸出的變化。

3.3.1水田改旱地

在黑水河流域內，農業種植中旱地占較大比例，旱地是單位面積非點源污染負荷輸出最大的土地利用類型，而水田的非點源污染負荷輸出略小于旱地。因此，將關鍵源區子流域的水田改為旱地，模擬非點源污染負荷，反推通過減少旱地等措施，降低非點源污染負荷。模擬結果顯示(表6)，通過將水田全部改為旱地，流域泥沙負荷增加了7.75%，有機磷污染負荷增加7.40%，有機氮污染負荷增加9.18%。

表6 水田改旱地情景下污染負荷的變化情況

3.3.2裸地改林地

在黑水河流域內，存在部分裸地，裸地是產沙、產污量較大的一種土地利用類型，將關鍵源區內裸地改為水土保持能力和污染物截留能力最好的林地，模擬在此情景下對流域內非點源污染的削減效果。模擬結果顯示，通過將裸地全部改為林地，對泥沙和有機氮、有機磷均有一定的削減效果，其中泥沙污染負荷削減率為0.33%；有機氮污染負荷削減率為0.30%；有機磷污染負荷削減率為0.28%。

表7 裸地改林地情景下污染負荷的削減效果

3.3.3退耕還林

將流域內坡度大于15°的耕地采取退耕還林措施，模擬退耕還林對流域出口處非點源污染的削減作用，模擬結果顯示(表8)，流域出口泥沙及有機污染負荷的輸出量均有一定削減效果，其中泥沙負荷削減率達到33.06%、有機氮污染削減率為18.82%、有機磷污染削減率為13.03%。

表8 退耕還林情景下流域出口污染負荷的削減效果

3.3.4設置緩沖帶

緩沖帶是指利用永久性植被攔截污染物或有害物質的條帶狀、受保護的土地[19]。通常在水體岸邊設置緩沖林帶，保護建立在河湖、溪流和溝谷沿岸的各類林帶或灌木帶，可以有效攔截農業面源污染，減少污染物向水體的排放。SWAT模型可以模擬緩沖林地對泥沙、有機污染物負荷的消減作用。

對關鍵源區子流域分別采取設置10、5 m寬植被緩沖帶措施，模擬不同寬度緩沖帶對流域內非點源污染的削減作用。模擬結果顯示(表9)，5 m植被緩沖帶對泥沙和有機磷、有機氮的削減率分別為49.35%、59.17%、59.16%；10 m植被緩沖帶對泥沙和有機磷、有機氮的削減率分別為60.69%、72.67%、72.67%。

表9 不同寬度緩沖帶情景下污染負荷的削減效果

3.4 討論

結合以上4種不同控制措施情景模擬，以及相應措施對非點源污染負荷削減效果，發現設置植被緩沖帶和退耕還林對非點源污染負荷削減效果比較明顯。其中10 m植被緩沖帶削減效果要比5 m植被緩沖帶的削減效果增加10%～15%，而同樣寬度下，植被緩沖帶對有機磷和有機氮的削減效果相當，退耕還林及旱地改水田情景模擬下，對有機氮的減負效果要好于有機磷，由于流域內裸地面積占比較小，采取裸地改林地后可以一定程度上降低非點源污染負荷，但效果不明顯。若組合以上情景措施，對流域非點源污染負荷的削減效率將會更高。

4 結論

a)建立了適用于黑水河流域SWAT模型，其中校驗期間參數En、R2均在合理范圍內，Re大部分也滿足小于15%的要求；利用校驗后的模型模擬了黑水河流域產水、產沙及產污負荷的空間分布，黑水河流域的非點源污染防治的關鍵源區主要是流域西北部和西南部的子流域，應重點關注。

b)對流域關鍵源區實施4種不同控制措施，模擬發現4種控制措施可以不同程度的削減流域內非點源污染負荷，其中植被緩沖帶、退耕還林的削減效果較好，分別達到了50%和15%左右。研究結果可以直接指導黑水河流域非點源污染負荷控制方案的制定，也可以為其他流域的非點源污染負荷控制研究提供參考。

c)除了本文中提出的改變土地利用方式的控制措施，還可以通過其他工程措施和管理措施來改善流域的水文水質環境，如加強城鎮基礎設施建設、應用農業最佳管理措施等，具體的措施如拆除部分水電站、將坡耕地改造為梯田、科學種植，合理施肥，減少化肥和農藥的使用量等。