基于SWAT模型的湯河流域面源污染時(shí)空分布研究

崔杰石

（遼寧省阜新水文局，遼寧阜新110000）

?

基于SWAT模型的湯河流域面源污染時(shí)空分布研究

崔杰石

（遼寧省阜新水文局，遼寧阜新110000）

摘要：本文運(yùn)用SWAT模型模擬湯河面源污染時(shí)空分布，并定量分析氣候變化對(duì)湯河面源污染的影響。研究結(jié)果表明：SWTA模型適用于湯河流域面源污染模擬，氣溫變化、降雨量變化都會(huì)引起湯河總氮和總磷的變化。研究成果為湯河面源污染模擬及污染控制提供了參考。

關(guān)鍵詞：SWAT模型；面源污染時(shí)空模擬；氣候變化情景模式；湯河流域

流域面源污染模擬一直是國(guó)內(nèi)外水污染模擬的難點(diǎn)和重點(diǎn)，面源污染由于缺乏面源污染監(jiān)測(cè)信息以及時(shí)空分布復(fù)雜性的特點(diǎn)，面源污染模擬效果一直難以達(dá)到滿意的模擬精度。由于地理信息技術(shù)和計(jì)算機(jī)水平發(fā)展，分布式水文水質(zhì)耦合模型為流域面源污染模擬提供一個(gè)有效的模擬平臺(tái)，其中由美國(guó)學(xué)者研究的分布式水文模型SWAT模型具有模型需求參數(shù)較少、操作簡(jiǎn)單的特點(diǎn)，被國(guó)內(nèi)外許多學(xué)者在不同流域的面源污染模擬進(jìn)行相關(guān)運(yùn)用。在國(guó)內(nèi)面源污染模擬研究也較多，但是還未在湯河流域的面源污染模擬中進(jìn)行過(guò)相關(guān)運(yùn)用，湯河流域內(nèi)有湯河水庫(kù)，該水庫(kù)是遼陽(yáng)市、鞍山市這兩大城市居民主要水源。本文引入在國(guó)內(nèi)運(yùn)用較好的SWAT模型，模擬湯河流域的面源污染時(shí)空分布。此外，考慮到氣候變化對(duì)流域水循環(huán)各個(gè)要素都將產(chǎn)生一定的影響，影響流域水環(huán)境各種污染物的遷移和轉(zhuǎn)化，因此在氣候變化大背景下，設(shè)定4種氣候變化情景模式分析氣候變化下湯河流域面源污染響應(yīng)。

1　研究方法及資料處理

1.1 研究區(qū)域概況

湯河流域?qū)儆谔雍幼蟀兑患?jí)支流，發(fā)源于遼寧省遼陽(yáng)縣吉洞峪滿族鄉(xiāng)粉城墻村，又名湯河?xùn)|支，流經(jīng)遼陽(yáng)縣、遼陽(yáng)市弓長(zhǎng)嶺區(qū)，在遼陽(yáng)縣小屯鎮(zhèn)雙廟村注入太子河。流域面積1422km2，河流長(zhǎng)度87km，河流平均比降2.82‰，多年平均年降水量763.0mm，多年平均年徑流深239.8mm，流域平均寬度為16.3km。

圖1　研究流域水系圖

1.2 研究方法

SWAT模型為水量水質(zhì)耦合的模型，模型可模擬不同形式下總氮、總磷的遷移轉(zhuǎn)化過(guò)程，模型具體模擬原理可見(jiàn)參考文獻(xiàn)［9］，本文運(yùn)用SWAT模型模擬湯河流域2000～2010年面源污染時(shí)空分布，并結(jié)合流域內(nèi)二道河子實(shí)測(cè)2000～2012年水質(zhì)監(jiān)測(cè)數(shù)據(jù)，對(duì)比分析SWAT模型在湯河流域面源污染模擬的適用性，并通過(guò)設(shè)定4種氣候變化情景模式，定量模擬不同氣候變化情景模式對(duì)湯河流域總氮和總磷的影響，并對(duì)比氣候變化前后流域總氮和總磷的變化率。

1.3 模型輸入資料準(zhǔn)備

圖2　研究流域遙感地理數(shù)據(jù)

模型需要輸入的數(shù)據(jù)包括：流域下墊面信息數(shù)據(jù)（數(shù)字高程數(shù)據(jù)、土地利用類型數(shù)據(jù)、土壤類型數(shù)據(jù)）、水文數(shù)據(jù)（降水、流量、含沙量）、氣象數(shù)據(jù)、水質(zhì)監(jiān)測(cè)數(shù)據(jù)（河道斷面水質(zhì)監(jiān)測(cè)數(shù)據(jù)）。其中：（1）流域數(shù)字高程數(shù)據(jù)下載于中科院地理研究所提供的數(shù)據(jù)共享平臺(tái)上全國(guó)1km× 1km的數(shù)字高程數(shù)據(jù)。結(jié)合流域邊界獲得研究區(qū)域的數(shù)字高程數(shù)據(jù)，見(jiàn)圖2，見(jiàn)參考文獻(xiàn)［10］；（2）土地利用/覆被數(shù)據(jù)。根據(jù)中科院地理所在數(shù)據(jù)共享平臺(tái)上提供了全國(guó)分辨率為1km×1km的土地利用/覆被數(shù)據(jù)［11］，結(jié)合流域邊界運(yùn)用相關(guān)軟件提取研究流域的土地利用/覆被數(shù)據(jù)，結(jié)果見(jiàn)圖2，對(duì)于SWAT模型而言，在土地利用/覆被數(shù)據(jù)，還應(yīng)該對(duì)數(shù)據(jù)進(jìn)行重分類，以滿足SWAT模型的數(shù)據(jù)格式要求，具體分類方法可詳見(jiàn)參考文獻(xiàn)［6］；（3）土壤類型信息數(shù)據(jù)。收集流域內(nèi)1∶100，0000土壤類型數(shù)據(jù)，結(jié)合流域邊界獲取研究流域土壤類型分布圖，見(jiàn)圖2。此外，模型還需要土壤化學(xué)屬性數(shù)據(jù)，對(duì)流域6種不同土壤類型進(jìn)行了采樣，測(cè)定了不同類型土壤的化學(xué)屬性數(shù)據(jù)，結(jié)果見(jiàn)表1。（4）水文數(shù)據(jù)。收集湯河流域內(nèi)二道河子水文站2000～2012年實(shí)測(cè)流量數(shù)據(jù)和流域內(nèi)4個(gè)降雨站點(diǎn)2000～2012年實(shí)測(cè)降水?dāng)?shù)據(jù)；（5）氣象數(shù)據(jù)。收集湯河流域內(nèi)遼陽(yáng)氣象站2000～2012年氣象要素?cái)?shù)據(jù)，氣象要素?cái)?shù)據(jù)主要用來(lái)計(jì)算流域的潛在蒸散發(fā)。（6）水質(zhì)數(shù)據(jù)。收集了湯河流域內(nèi)二道河子水質(zhì)監(jiān)測(cè)點(diǎn)2000～2012年水質(zhì)監(jiān)測(cè)數(shù)據(jù)。（7）農(nóng)業(yè)耕種措施數(shù)據(jù)。結(jié)合當(dāng)?shù)貙?shí)地農(nóng)業(yè)調(diào)查數(shù)據(jù)，建立模型農(nóng)業(yè)耕種措施數(shù)據(jù)庫(kù)，包括耕種時(shí)間、施肥時(shí)間、施肥量數(shù)據(jù)。

表1　土壤化學(xué)屬性數(shù)據(jù)

2　模型運(yùn)用

2.1 總氮和總磷模擬

基模型輸入數(shù)據(jù)，結(jié)合湯河流域內(nèi)二道河子站2000～2012年水質(zhì)監(jiān)測(cè)數(shù)據(jù)，對(duì)SWAT總氮和總氮模擬參數(shù)的率定和驗(yàn)證，結(jié)果見(jiàn)下表和圖。

表2和表3為運(yùn)用SWAT模型模擬的湯河流域2000～2012年總氮和總磷結(jié)果，從表中可以看出，SWAT模型在湯河流域模擬的2000～2012年總氮和總磷相對(duì)誤差均小于20%，模擬的確定性系數(shù)均高于0.65，滿足流域面源污染模擬的精度規(guī)范要求。圖3和圖4為選定的兩個(gè)代表年份實(shí)測(cè)總氮和總磷與模擬的總氮和總磷的過(guò)程對(duì)比圖，從圖中可以看出，模擬的總氮和總磷與實(shí)測(cè)的總氮和總磷具有一定的吻合度，模擬的總氮和總磷的峰值與實(shí)測(cè)總氮和總磷出現(xiàn)的峰值月份吻合。

表2　SWAT模型總氮模擬結(jié)果

表3　SWAT模型總磷模擬參數(shù)率定結(jié)果

2.2 不同氣候模式情景對(duì)湯河流域總氮和總磷的影響分析

為定量分析氣候變化下湯河流域面源污染負(fù)荷響應(yīng)，以SWAT模型為模擬分析平臺(tái)，選定2011年（2000～2012年降水量最大的年份）為代表年份，運(yùn)用SWAT模型模擬四種氣候變化情景模式下（見(jiàn)表4）2011年總氮和總磷的變化量以及變化率，從而定量分析氣候變化對(duì)湯河流域總氮和總磷的影響。

圖3　代表年份總氮模擬成果圖

圖4　代表年份總磷模擬成果圖

表4　設(shè)定的4種氣候變化情景模式

表5　氣溫變化情景方案下的總氮響應(yīng)

表6　氣溫變化情景方案下的總磷響應(yīng)

從表5和表6中可以看出，在氣候變化情景模式I下（氣溫升高+2℃，降水不變化），2011年湯河流域總氮和總磷含量在氣候變化后，含量分別減少3.15萬(wàn)t和0.024萬(wàn)t，變化率為-11.0%和-10.2%，可見(jiàn)氣溫升高，流域總氮含量將減少一定的比例，這主要是因?yàn)榱饔驓鉁厣?，流域潛在蒸散發(fā)相應(yīng)增加，在降水不變化的前提下，流域地表徑流相應(yīng)減少，而作為總氮和總磷的運(yùn)移載體，總氮和總磷的含量相應(yīng)減少。而在氣候變化情景模式II下（氣溫降低-2℃，降水不變化），流域總氮和總磷都增加一定的比例，這主要是因?yàn)闅鉁亟档?，使得流域潛在蒸散發(fā)增加，使得流域地表徑流增加，作為面源污染運(yùn)移載體，隨著地表徑流流失的總氮和總磷相應(yīng)減少。而在氣候變化情景模式III下（降水+5%，氣溫不變化），相比于氣候變化前，湯河流域2011年總氮和總磷分別增加28.1% 和26.4%，總氮和總磷均增加的原因是因?yàn)榻邓吭黾?，流域地表徑流相?yīng)增加，使得隨著地表徑流流失的總氮和總磷含量相應(yīng)增加。而在氣候變化情景模式IV下，降水減少5%，氣溫不變化，流域總氮和總磷相應(yīng)減少一定的比例，這同樣因?yàn)榻邓疁p少，使得流域地表徑流深相應(yīng)減少，隨地表徑流流域的面源污染負(fù)荷也相應(yīng)減少，從表中還可以看出，降水量變化引起的總氮和總磷的變化率均要大于氣溫變化引起的總氮和總磷的變化率，可見(jiàn)，降水變化引起的面源污染負(fù)荷影響要大于氣溫變化的影響。

3　結(jié)語(yǔ)

本文結(jié)合SWAT模型模擬湯河流域面源污染時(shí)空分布，并基于SWAT模型定量分析4種氣候變化情景模式對(duì)流域面源污染負(fù)荷的影響，研究取得以下結(jié)論。

（1）SWAT模型適用于湯河流域面源污染模擬，模型模擬的各年份總氮和總磷年相對(duì)誤差小于15%，月尺度過(guò)程確定性系數(shù)達(dá)到0.65以上，符合流域非點(diǎn)源污染模擬精度要求，可用來(lái)模擬分析不同情景下流域的面源污染負(fù)荷時(shí)空分布。

（2）降水量變化將直接影響湯河流域面源污染負(fù)荷含量，且影響率大于氣溫變化引起的湯河流域面源負(fù)荷變化率。

參考文獻(xiàn)

［1］張永勇，王中根，于磊，夏軍，陳向東.SWAT水質(zhì)模塊的擴(kuò)展及其在海河流域典型區(qū)的應(yīng)用［J］.資源科學(xué)，2009（01）：94-100.

［2］楊國(guó)敏，王曉鳳，成杰民，陳雪蘭.GIS技術(shù)在東昌湖水質(zhì)評(píng)價(jià)中的應(yīng)用［J］.環(huán)境監(jiān)測(cè)管理與技術(shù)，2013（02）：57-61.

［3］李穎，王康，周祖昊.基于SWAT模型的東北水稻灌區(qū)水文及面源污染過(guò)程模擬［J］.農(nóng)業(yè)工程學(xué)報(bào)，2014（07）：42-53.

［4］李丹，薛聯(lián)青，郝振純.基于SWAT模型的流域面源污染模擬影響分析［J］.環(huán)境污染與防治，2008（03）：4-7.

［5］周浩.大小凌河地區(qū)地下水資源可利用量預(yù)測(cè)模擬研究［J］.水利技術(shù)監(jiān)督，2015（02）：16-20.

［6］孫小利，王建躍.美國(guó)烏龜溪水庫(kù)降低水藻生長(zhǎng)控制元素磷負(fù)荷的措施［J］.水利規(guī)劃與設(shè)計(jì)，2011（04）：24-26.

［7］趙洋甬，馬靜軍，楊樹(shù)生，肖國(guó)起，蔣蕾蕾.寧波市亭下水庫(kù)總氮超標(biāo)情況成因分析［J］.水利技術(shù)監(jiān)督，2010（05）：13-15.

［8］孫本發(fā).SWAT模型在巢湖典型監(jiān)測(cè)點(diǎn)農(nóng)業(yè)面源污染模擬的評(píng)價(jià)與驗(yàn)證研究［D］.安徽農(nóng)業(yè)大學(xué)，2013.

［9］陳軍寧.基于SWAT模型的鄱陽(yáng)湖區(qū)非點(diǎn)源污染模擬研究［D］.南昌大學(xué)，2013.

［10］DEM數(shù)據(jù)：http：//srtm.csi.cgiar.org/SELEGTION/inputGoord.asp.

［11］土地利用數(shù)據(jù)：http：//westdc.westgis.ac.cn/data/1cad1a63-ca8d-431a-b2b2-45d9916d860d/.

作者簡(jiǎn)介：崔杰石（1966年—），男，高級(jí)工程師。

收稿日期：2015-10-27

DOI：10.3969 /j.issn.1672-2469.2016.02.002

中圖分類號(hào)：X522

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1672-2469（2016）02-0004-03