閆雪嫚，盧文喜，宋文博，趙 瑩

(1．吉林大學地下水資源與環(huán)境教育部重點實驗室，長春 130021；2. 吉林大學環(huán)境與資源學院，長春 130021)

0 引 言

水污染問題是當前嚴重的環(huán)境問題之一，準確評價水污染造成的經(jīng)濟損失不僅是進行水環(huán)境規(guī)劃的基礎(chǔ)，更是制定水環(huán)境政策的先決條件[1]。目前，農(nóng)業(yè)非點源污染對地表水環(huán)境的影響日益顯著。調(diào)查資料表明，非點源污染對全球30%到50%的地表水體造成了影響，其中的農(nóng)業(yè)非點源污染是造成地表水水質(zhì)惡化的主要原因[2]。農(nóng)業(yè)非點源污染不僅會危害農(nóng)業(yè)安全，還會帶來巨大的經(jīng)濟損失[3]。因此，運用科學有效的方法對農(nóng)業(yè)非點源污染損失進行定量估算是十分必要的。

伊通河流域農(nóng)牧業(yè)發(fā)達，吉林省乃至全國的主要商品糧生產(chǎn)基地坐落在流域內(nèi)。伊通河兩岸有大面積的耕地，農(nóng)藥和化肥的大量使用在提高糧食產(chǎn)量的同時，也給伊通河帶來了嚴重的污染。隨著點源污染治理工作的不斷加強，非點源污染特別是農(nóng)業(yè)非點源污染已經(jīng)成為伊通河流域面臨的主要問題。而在我國環(huán)境與經(jīng)濟協(xié)調(diào)發(fā)展的背景下，對農(nóng)業(yè)非點源污染進行分析研究時不僅要評價其污染狀況，同時還要分析其經(jīng)濟損失。目前，有很多學者對伊通河流域農(nóng)業(yè)非點源污染進行了模擬研究，但還十分缺少采用機理模型從經(jīng)濟角度對農(nóng)業(yè)非點源污染損失進行分析和評估的研究。

目前，在我國，針對農(nóng)業(yè)非點源污染的研究主要集中在3個方面：調(diào)查評價、模型研究和防治技術(shù)[4]。雖然已有針對流域農(nóng)業(yè)非點源污染損失的相關(guān)研究[5-6]，但在對農(nóng)業(yè)非點源污染損失進行定量估算時，大多數(shù)學者選擇采用輸出系數(shù)法計算農(nóng)業(yè)非點源污染負荷。由于所采用的某些輸出系數(shù)不統(tǒng)一、統(tǒng)計數(shù)據(jù)精度較低以及數(shù)據(jù)不全等問題，采用輸出系數(shù)法計算的農(nóng)業(yè)非點源污染負荷精度較低，直接導致估算的農(nóng)業(yè)非點源污染損失準確度降低。已有研究表明，運用SWAT等機理模型計算農(nóng)業(yè)非點源污染負荷，其精度通常比輸出系數(shù)法高[7]。但目前，基于機理模型的農(nóng)業(yè)非點源污染損失研究還比較少見。

本文以伊通河流域為例，針對伊通河流域的農(nóng)業(yè)非點源污染損失開展研究。與以往不同的是，在對研究區(qū)農(nóng)業(yè)非點源污染損失進行定量估算時，選用SWAT模型計算伊通河流域農(nóng)業(yè)非點源污染負荷。在此基礎(chǔ)上采用影子工程法、機會成本法及恢復費用法分別對泥沙淤積造成的損失、土地廢棄造成的損失及土壤養(yǎng)分流失造成的損失進行估算。以此提高伊通河流域農(nóng)業(yè)非點源污染損失估算的精度，使之更加切合實際，旨在為伊通河流域農(nóng)業(yè)非點源污染的防治和管理提供科學依據(jù)與決策支持。

1 研究區(qū)概況

伊通河流域地處吉林省中部，東北平原腹地，位于東經(jīng)124°32′23″～125°46′42″，北緯43°01′58″～44°50′33″。伊通河是松花江的二級支流，全長共342.5 km，是長春市的母親河，伊通河流域面積為8 840 km2。伊通河流域處于平原地帶，地勢平坦開闊，是吉林省乃至全國重要的糧食生產(chǎn)基地，耕地面積約為6 870 km2，約占總面積的78%；而草地面積及林地面積分別為368、506 km2，所占比例較小。

2 研究方法

2.1 農(nóng)業(yè)非點源污染負荷計算方法

本次研究采用SWAT模型對伊通河流域農(nóng)業(yè)非點源負荷進行計算。SWAT模型的發(fā)展是以SWRRB(Simulator for Water Resources in Rural Basins)模型為基礎(chǔ)的[8]。它具有很強的物理基礎(chǔ)，采用以日為單位時間進行連續(xù)計算，是一種基于GIS平臺的分布式流域水文模型。SWAT模型由美國農(nóng)業(yè)部農(nóng)業(yè)研究局所研發(fā)，在流域水平衡、河流流量預測和非點源污染控制評價等諸多方面都得到了廣泛的應(yīng)用[9,10]。

SWAT模型所需數(shù)據(jù)由地形、土壤、土地利用、水文和氣象數(shù)據(jù)組成。其中氣象數(shù)據(jù)時段為2008-2012年，包括長春、農(nóng)安、雙陽、煙筒山、四平5個氣象站的逐日氣象數(shù)據(jù)(降水量、最高和最低氣溫、相對濕度、太陽輻射、風速)，數(shù)據(jù)來源及說明見表1。

表1 數(shù)據(jù)來源及說明Tab.1 Source of data and instructions

在相關(guān)數(shù)據(jù)準備完善的基礎(chǔ)上，按照順序輸入數(shù)字高程模型、土地利用類型圖、土壤類型圖等初步建立農(nóng)業(yè)非點源污染模型。然后采用SUFI-2算法進行參數(shù)的率定及模型的驗證，最終建立伊通河流域農(nóng)業(yè)非點源污染模擬模型。本文選用確定性系數(shù)(R2)以及模型效率系數(shù)Nash-Sutcliffe[11]對模型的模擬效果進行評價。

R2的計算公式如下：

(1)

Ens的計算公式如下：

(2)

理論上，R2和Ens的值越接近1，說明模型模擬效果越好[12]。當R2>0.6、Ens>0.5時，可說明模型模擬效果較好、結(jié)果可靠。

在建立伊通河流域農(nóng)業(yè)非點源污染SWAT模型的基礎(chǔ)上，分析水文響應(yīng)單元輸出文件、子流域輸出文件以及主河道或河段輸出文件等多個輸出文件，對伊通河流域農(nóng)業(yè)非點源污染負荷進行計算。

2.2 農(nóng)業(yè)非點源污染損失估算方法

運用不同的環(huán)境經(jīng)濟學分析方法對不同類型的農(nóng)業(yè)非點源污染損失進行計算。

采用影子工程法[13]計算由泥沙淤積所引起的環(huán)境損失，其中泥沙淤積是由土壤侵蝕造成的：

Ds=ZsPs/ρs

(3)

式中：Ds是由泥沙淤積所引起的環(huán)境損失，元；Zs是研究區(qū)每年總的泥沙流失量，t；Ps是為攔截每單位體積泥沙所建工程的投資費用，目前其庫容造價是120 元/m3；ρs是泥沙容重，t/m3，經(jīng)驗取值1.28。

采用機會成本法計算因土地廢棄而失去的經(jīng)濟損失價值。這里的土地廢棄包括兩種情況：①全部耕層喪失生態(tài)價值而廢棄；②因耕層變薄更加貧瘠而廢棄，因此土壤流失總體積按平均耕層厚度折算成土地面積來計算[14]：

Dd=ZB/(0.3×104ρ)

(4)

式中：Dd是土壤廢棄的經(jīng)濟損失價值，元/a；Z是研究區(qū)每年土壤侵蝕總量，t/a；B是土地損失的機會成本，目前造價為1.5萬元/hm2；ρ是土壤容重，t/m3，經(jīng)驗取值1.10。

采用恢復費用法計算土壤養(yǎng)分和有機質(zhì)流失所造成的環(huán)境損失，該方法假設(shè)人們會為了保護自己避免危險而支付錢幣，因此可以用恢復或防護資源不受污染所需支出的費用，作為環(huán)境資源破壞帶來的經(jīng)濟損失的最低估計[13]，按式(5)計算：

Si=TiKiC

(5)

式中：Si是第i種養(yǎng)分流失所損失的價值，元；i是N或P元素；Ti是農(nóng)業(yè)非點源污染土壤第i種養(yǎng)分流失總量，t；Ki是第i種養(yǎng)分折算為磷酸二銨的系數(shù)；C是磷酸二銨肥料的價格，元。根據(jù)中國化肥網(wǎng)(http:∥www.fert.cn/)得知目前磷酸二銨的市場價格為3 200元/t ，N、P折算成磷酸二銨的系數(shù)分別是132/14，132/31。

3 結(jié)果與分析

3.1 農(nóng)業(yè)非點源污染負荷計算結(jié)果

本文根據(jù)研究區(qū)農(nóng)安站和伊通站的實測月徑流資料以及農(nóng)安站的月輸沙量數(shù)據(jù)進行了參數(shù)的率定及模型的驗證。其中，率定期為2008-2010年，驗證期為2011-2012年。模擬模型參數(shù)的率定及模型的驗證結(jié)果如表2所示。其中徑流率定R2均大于0.90，Ens均大于0.85；驗證期R2均大于0.85，Ens均大于0.80。泥沙模擬率定期R2為0.77，Ens為0.83；驗證期R2為0.76，Ens為0.75。總體上模型模擬結(jié)果符合SWAT模型模擬要求，可以應(yīng)用于伊通河流域研究區(qū)農(nóng)業(yè)非點源污染負荷的模擬計算。

表2 伊通河流域農(nóng)業(yè)非點源污染模型模擬效果Tab.2 The simulation effect of agricultural non-point source pollution model in Yitong River Watershed

在此基礎(chǔ)上，分析模型的輸出文件，計算伊通河流域農(nóng)業(yè)非點源污染負荷，其結(jié)果如表3所示；伊通河流域土壤侵蝕總量的估算結(jié)果如表4所示。由表3可知，在2008-2012這5年內(nèi)，2010年的農(nóng)業(yè)非點源泥沙年負荷、總氮年負荷以及總磷年負荷最大，分別為3 417 200 t、12 776.40 t以及2 614.60 t，2011年最?。黄浯笮『湍杲涤炅坑泻艽蟮年P(guān)系。此外，由表4可以看出，伊通河流域的土壤侵蝕量達到了4 202 877.16 t/a。其中，耕地引起的土壤侵蝕量最大，達到了4 056 071.48 t/a，占總量的96.51%；而水田引起的土壤侵蝕量最小，僅占總量的0.34%。

表3 研究區(qū)2008-2012年農(nóng)業(yè)非點源污染年負荷計算結(jié)果Tab.3 Agricultural non-point source pollution load of the study area in 2008-2012

表4 研究區(qū)土壤侵蝕總量計算結(jié)果Tab.4 The total amount of soil erosion of the study area

3.2 農(nóng)業(yè)非點源污染損失估算結(jié)果

根據(jù)伊通河流域農(nóng)業(yè)非點源污染負荷的計算結(jié)果以及伊通河流域土壤侵蝕總量的計算結(jié)果，結(jié)合式(3)～(5)對伊通河流域農(nóng)業(yè)非點源污染損失進行估算，估算的結(jié)果如表5所示?？梢钥闯觯捎诠浪愕难芯繀^(qū)農(nóng)業(yè)非點源年負荷的不同，估算所得的研究區(qū)農(nóng)業(yè)非點源污染年損失也不同，年降雨量越大的年份損失越嚴重。在2008-2012年這5年內(nèi)，研究區(qū)年均農(nóng)業(yè)非點源污染損失約為52 845.42萬元。其中，2010年的農(nóng)業(yè)非點源污染損失最大，為76 057.47萬元；2012次之，為57 047.08萬元；2011年最小，為40 243.19萬元?？傮w上，伊通河流域內(nèi)農(nóng)業(yè)非點源污染價值損失大小順序為：總氮價值損失>泥沙淤積價值損失>總磷污染價值損失>土地廢棄價值損失，表明在伊通河流域農(nóng)業(yè)非點源污染引起的損失中，以氮養(yǎng)分流失和泥沙淤積為主。以2009年為例進行分析，總氮污染損失達到了20 421.53 萬元，約占總損失的49.42%；泥沙淤積價值損失次之，為16 932.19萬元，約占總損失的40.98%；磷污染價值損失較小，為2151.38萬元，約占總損失的5.29%；土地廢棄價值損失最小，為1 910.40萬元，僅占總損失的4.31%。由此可見，伊通河流域農(nóng)業(yè)非點源污染控制的重點應(yīng)該放在對總氮及泥沙的治理上。

由于沒有考慮到2008-2012年這5年內(nèi)研究區(qū)內(nèi)各土地利用類型面積的變化，估算的農(nóng)業(yè)非點源污染損失可能與實際有偏差。

表5 研究區(qū)農(nóng)業(yè)非點源污染損失計算結(jié)果 萬元

4 結(jié)論與建議

4.1 結(jié) 論

采用SWAT模型對伊通河流域農(nóng)業(yè)非點源污染負荷進行模擬計算，在此基礎(chǔ)之上，運用環(huán)境經(jīng)濟學的方法定量估算研究區(qū)農(nóng)業(yè)非點源污染損失，得出的主要結(jié)論如下。

(1)在伊通河流域農(nóng)業(yè)非點源污染模擬中，模型模擬結(jié)果符合SWAT模型模擬要求，可以應(yīng)用于伊通河流域研究區(qū)農(nóng)業(yè)非點源污染負荷的模擬計算。

(2)模型模擬結(jié)果表明，2010年內(nèi)農(nóng)業(yè)非點源泥沙年負荷、總氮年負荷以及總磷年負荷均為最大，2012年次之，2011年最小；其大小和年降雨量有很大的關(guān)系。伊通河流域的土壤侵蝕總量為4 202 877.16 t/a。其中，耕地引起的土壤侵蝕量最大，占總量的96.51%；而水田引起的土壤侵蝕量最小。

(3)根據(jù)估算結(jié)果可知，年降雨量越大的年份損失越嚴重。研究區(qū)年均農(nóng)業(yè)非點源污染損失折算成人民幣約為52 845.42 萬元?？傮w上，伊通河流域內(nèi)農(nóng)業(yè)非點源污染價值損失大小順序為：總氮價值損失>泥沙淤積價值損失>總磷污染價值損失>土地廢棄價值損失。

4.2 建 議

(1)在今后的伊通河流域農(nóng)業(yè)非點源污染治理工作中應(yīng)該加強對土壤侵蝕的控制與管理。

(2)伊通河流域農(nóng)業(yè)非點源污染控制的重點應(yīng)該放在對總氮及泥沙的治理上，特別是關(guān)于總氮的治理應(yīng)為重中之重。此外，雖然流域內(nèi)農(nóng)業(yè)非點源磷污染價值損失占總損失的比例較小，但危害大，也應(yīng)該引起足夠的重視。

[1] 李 昕，董德明，沈萬斌，等.綠色國民經(jīng)濟核算基本問題研究[J]. 地理科學， 2007， 27(2):163-166.

[2] Dennis L C, Peter J V, Keith I. Modeling non-point source pollution in vadose zone with GIS [J]. Environmental Science and Technology, 1997，35(8):2 157-2 175.

[3] 查恩爽.伊通河流域農(nóng)業(yè)非點源污染模擬及最佳管理措施的應(yīng)用[D].長春：吉林大學，2011：3-7.

[4] 劉春光，石 靜. 非點源污染研究進展及趨勢分析[J]. 農(nóng)業(yè)環(huán)境科學學報，2008，27(4)：1 306-1 311.

[5] 段雪梅.平原河網(wǎng)區(qū)農(nóng)業(yè)非點源污染負荷及經(jīng)濟損失估算研究[D].江蘇揚州：揚州大學，2013:89-101.

[6] 鮑秋萍.農(nóng)業(yè)非點源污染氮、磷負荷估算及經(jīng)濟損失評估——以福建省2010年農(nóng)業(yè)非點源污染為例[J].福建工程學院學報，2012，10(6)：621-624.

[7] 劉 莊，晁建穎，張 麗，等.中國非點源污染負荷計算研究現(xiàn)狀與存在問題[J]. 水科學進展，2015，26(5)：432-442.

[8] Lai G Y，Wu D Y，Zhong Y X，et al.Progress in development and applications of SWAT model．Journal of Hohai University: Natural Sciences，2012，40( 3) : 243-251．

[9] 李穎，王康，周祖昊. 基于SWAT模型的東北水稻灌區(qū)水文及面源污染過程模擬[J].農(nóng)業(yè)工程學報，2014，30(7)：42-53.

[10] Cho J，Lowrance R R，Bosch D D，et al.Effect of watershed subdivision and filter width on SWAT simulation of a coastal plain watershed[J].Journal of the American Water Resources Association，2010，46(3)：586-602.

[11] 項 瑛，張余慶，程 婷. 基于SWAT模型的里下河平原區(qū)水文模擬[J].資源科學，2015，37(6)：1 181-1 189.

[12] Moon J，June Y，Lee T, et al.Determining the effective width of riparian buffers in Korean watershed using the SWAT model[J].Environmental Enginering and Management Journal，2013，12(11)：2 249-2 260.

[13] 楊志新，鄭大瑋，李永貴.北京市土壤侵蝕經(jīng)濟損失分析及價值估算[J].水土保持學報，2004，18(3):175 -178.

[14] 王曉燕，張雅帆，歐 洋.北京密云水庫上游太師屯鎮(zhèn)非點源污染損失估算[J].生態(tài)與農(nóng)村環(huán)境學報，2009，25(4)：37-41.

[15] 姚 建.環(huán)境經(jīng)濟學[M]. 成都:西南財經(jīng)大學出版社，2001：25 -60.