杜　娟，李懷恩，趙湘璧，李家科

(1.榆林學(xué)院建筑工程系，陜西 榆林　719000；2.西安理工大學(xué)西北水資源與環(huán)境生態(tài)教育部重點實驗室，陜西 西安　710048)

?

基于SWAT模型的渭河流域污染控制效果模擬

杜娟1，李懷恩2，趙湘璧1，李家科2

(1.榆林學(xué)院建筑工程系，陜西 榆林719000；2.西安理工大學(xué)西北水資源與環(huán)境生態(tài)教育部重點實驗室，陜西 西安710048)

摘要：在實測數(shù)據(jù)的基礎(chǔ)上，利用率定驗證后的SWAT模型對渭河流域污染負荷進行模擬計算，并對非點源和點源污染控制措施的模擬效果進行評價。結(jié)果表明：非點源污染綜合控制措施可較好地削減污染物，且削減效果明顯高于單一非點源措施。點源污染控制措施對N、P削減效果較明顯，且削減效果隨著強度的增加而增強。點源和非點源相結(jié)合的污染控制綜合方案對流域N、P的削減效果比單一措施效果明顯，采取綜合措施后2020年渭河流域污染負荷產(chǎn)出均有所下降，且滿足“十二五”規(guī)劃中主要污染物控制目標，說明綜合措施治理效果較好，流域水環(huán)境質(zhì)量得到明顯提高。

關(guān)鍵詞：渭河流域；SWAT模型；非點源污染；點源污染；污染控制措施

按照污染源發(fā)生類型,水環(huán)境污染可分為點源污染和非點源污染2大類。即使點源污染得到全面控制,水質(zhì)達標率也僅為42%～65%[1]。我國是農(nóng)業(yè)大國,人口持續(xù)增加及不合理人類活動的增多加劇了水體非點源污染,導(dǎo)致可用水資源量大幅度減少。目前,我國正處于水污染構(gòu)成的轉(zhuǎn)變初期,在總污染負荷中點源污染貢獻率逐漸下降,非點源污染貢獻率逐步上升。若能在點源得到削減的基礎(chǔ)上,加強對非點源污染的控制治理,必然能夠抑制水質(zhì)惡化,改善生態(tài)環(huán)境。

渭河流域貫穿陜甘寧3省,是我國重要的工業(yè)生產(chǎn)基地和糧棉油產(chǎn)區(qū),更是關(guān)中地區(qū)的重要水源地,直接影響著兩岸各地的經(jīng)濟和生活水平。作為關(guān)中地區(qū)唯一的污廢水排泄通道,流域內(nèi)大部分生活污水、工業(yè)廢水、農(nóng)業(yè)污水等進入渭河,導(dǎo)致渭河干流水質(zhì)惡化為Ⅴ類,甚至劣Ⅴ類,水體基本喪失原有使用功能。該研究旨在利用SWAT模型識別流域內(nèi)污染來源,并通過定量化的模擬分析不同污染控制方案的效果,提出治理污染的有效控制措施,為渭河流域的水環(huán)境治理提供理論依據(jù)。

1材料與方法

1.1研究區(qū)概況

渭河是黃河水系最大的一級支流,全長818 km,流域面積13.48萬km2,涉及陜甘寧3個省(自治區(qū))的10個市(地區(qū)、州、盟)共72個縣(旗、市)。全流域按河段劃分為3段:渭河源頭鳥鼠山至寶雞峽以上為上游段,屬山區(qū)河流:河谷狹窄,川峽相間,水流湍急,全長430 km;寶雞峽至咸陽鐵橋為中游段,屬山區(qū)向平原過渡區(qū)河流:河道寬,多沙洲,水流分散,全長177 km;咸陽鐵橋至潼關(guān)口為下游段,屬典型平原區(qū)河流:水流緩慢,比降小,河道淤積嚴重,全長211 km(圖1)。

圖1　渭河流域水系分布

1.2研究方法

選擇渭河流域華縣水文站斷面以上流域部分,在現(xiàn)有研究資料(1986—1990年)基礎(chǔ)上,加入1990—2010年數(shù)據(jù),利用SWAT模型對該區(qū)域非點源和點源污染進行模擬,識別區(qū)域內(nèi)主要污染物,包括有機氮(ORGN)、有機磷(ORGP)、氨氮(以NH4+-N計)、礦物磷(MINP)等的污染產(chǎn)出,結(jié)合景觀生態(tài)學(xué)理論與方法,定量模擬、評價土地利用結(jié)構(gòu)的調(diào)整以及農(nóng)業(yè)管理措施的改變等非點源控制措施效果,提出主要非點源污染物的控制管理措施和對策;加入點源負荷后,對點源和非點源負荷產(chǎn)出進行模擬,分析綜合措施的模擬效果。

1.3模型構(gòu)建

在已建立的模型基礎(chǔ)數(shù)據(jù)庫[2]中,添加1990—2010年的水文、氣象和水質(zhì)資料,包括日最低最高氣溫、日風(fēng)速和日蒸發(fā)量等氣象資料,降雨、徑流和泥沙數(shù)據(jù),NH4+-N、NO3--N和NO2--N等水質(zhì)資料以及1986—2009年流域內(nèi)各大排污單位的數(shù)據(jù)(地理位置,廢、污水排放量,污染負荷量等),并將上述圖形文件進行坐標轉(zhuǎn)換,氣象水質(zhì)等資料轉(zhuǎn)化為DBF格式,輸入模型生成新的SWAT模型數(shù)據(jù)集。根據(jù)流域數(shù)字高程模型(DEM)和水系劃分子流域及水文響應(yīng)單元(HRU)進行參數(shù)敏感性分析,加入點源資料后,應(yīng)用實測數(shù)據(jù)對徑流、水質(zhì)等相關(guān)系數(shù)進行率定、驗證,選用相對誤差(ER)、確定性系數(shù)(R2)和Nash-Suttcliffe模擬效率系數(shù)(Ens)對模型的適用性進行評價[3-4](表1和圖2)。

表1模型適用性評價

Table 1Applicability of the SWAT model

項目相對誤差(ER)/%R2Nash-Suttcliffe模擬效率系數(shù)(Ens)率定期驗證期率定期驗證期率定期驗證期徑流-1.27-7.900.930.870.880.84泥沙-16.6625.120.860.780.690.66氨氮-27.55-25.990.910.890.830.89硝態(tài)氮11.5725.630.780.690.590.66亞硝態(tài)氮-20.28-11.130.870.610.800.61

月均徑流模擬值與實測值的相對誤差ER在±10% 內(nèi),Ens≥0.5且R2≥0.8,說明調(diào)參后的模型滿足地表徑流的模擬精度要求;月均泥沙模擬值與實測值的相對誤差ER在±30%內(nèi),Ens≥0.5且R2≥0.7,可認為泥沙部分的精度基本滿足模型模擬要求;月均NH4+-N、NO3--N、NO2--N的模擬值與實測值的相對誤差ER均在±30%以內(nèi),Ens≥0.5 且R2≥0.6,說明模擬效果較好,基本滿足模擬要求。在完成徑流、泥沙和水質(zhì)的率定驗證后,與其相關(guān)的敏感參數(shù)值得到確定(表2)。

2結(jié)果與討論

2.1不同代表年非點源污染模擬結(jié)果

對SWAT模型進行驗證后,將1955—2009年共55 a的長系列實測渭河徑流量作為基礎(chǔ)資料,利用P-Ⅲ型水文頻率分析軟件對徑流量進行頻率分析,確定不同代表年,并模擬相應(yīng)代表年的非點源污染負荷產(chǎn)出,結(jié)果見表3。由表3可知,豐水年、平水年、枯水年的降雨量、徑流量和泥沙量均依次減小,且與當年降雨徑流量有密切關(guān)系。泥沙量與N、P非點源負荷分布有很好的相關(guān)性,說明控制渭河流域的水土流失和減少土壤侵蝕可以減小流域非點源負荷輸出。另外,研究流域人口多,密度大,加之耕地、農(nóng)村用地的化肥施用量不斷增大,有機肥施用量下降等原因,都造成渭河流域非點源污染負荷比例的增加。

圖2　華縣水文站率定期(1987—1988年)和驗證期(1995和2009年)污染負荷模擬值與實測值對比

2.2渭河流域非點源污染控制措施效果模擬

在流域建立5種方案25個情景對不同農(nóng)業(yè)管理污染控制措施的效果進行模擬,具體方案見表4。對以上方案進行綜合分析,以確定不同措施對水質(zhì)指標的靈敏度。選用平水年非點源污染負荷量來說明各措施的削減效果(表5)。從表5可知,采取水土保持措施(方案1)后N、P的削減效果較明顯;其中ORGN的削減率最高達68.72%,ORGP達58.71%,NH4+-N達56.59%,MINP達58.32%,TP達58.53%,NO3--N削減效果相對較差。與情景1-6相比,減少表層土壤施肥量占比(情景2-4)對NSP負荷的削減較少,但對MINP和TP的削減較明顯,MINP削減率可達7.30%,TP削減率達5.23%。合理減少施肥量(方案3)對P、NO3--N的削減效果較為明顯,且隨著化肥施用量的減少,削減率隨之增加,MINP達5.25%,TP為4.70%,NO3--N為4.05%。改進灌溉方式并減少灌溉用水量對N、P負荷流失的影響比施肥量更加明顯。農(nóng)田中N流失主要以硝態(tài)N為主,P則以溶解性P和顆粒P為主,這是因為P主要與土壤懸浮顆粒有關(guān),“故灌溉量增加時,P的流失會急劇增大,且多為可溶性磷[5]”。采用現(xiàn)代農(nóng)業(yè)中推廣的灌溉方式和節(jié)水措施(方案4)進行模擬,對NSP負荷的削減量較大,其中ORGN削減率最高達27.19%,NO3--N削減率達14.36%,MINP削減率達65.08%,TP削減率達32.29%。同時采取水保措施、合理減少施肥量和表層土壤施肥量占比并改進灌溉方式(情景5-4)顯示對NSP負荷的削減效果顯著,其中ORGN達29.42%,ORGP達32.90%,NO3--N達20.23%,MINP達29.45%,TP達30.67%,由此可見,綜合措施對NSP各項指標均有很好的削減作用,可以有效控制流域內(nèi)NSP負荷的輸出。

表2模型參數(shù)率定結(jié)果

Table 2Calibration of parameters of the model

變量 模擬過程定義 典型值范圍 參數(shù)最終值 CN2徑流徑流曲線數(shù)0～100;微調(diào)范圍±8%45;微調(diào)值-5.9%CANMX徑流樹冠最大截留量0～1000.50mmSOL-AWC徑流土壤水有效利用量0～1;微調(diào)范圍±8%0.22;微調(diào)值-0.03%EPCO徑流植物吸收補償系數(shù)0～10.85ESCO徑流土壤蒸發(fā)補償系數(shù)0～11.00GWQMN徑流最小基流出流閾值0～5000280.00GW-REVAP徑流土壤再蒸發(fā)系數(shù)0.02～0.200.02REVAPMN徑流土壤再蒸發(fā)閾值0～5000.50SMFMX徑流6月21日雪融系數(shù)0～1010SMFMN徑流12月21日雪融系數(shù)0～106.8USLE-C泥沙土壤植被管理因子0.003～0.450耕地0.100;林地0.038;未利用荒草地0.040Spexp泥沙泥沙輸移指數(shù)系數(shù)1.00～1.501.20PRF泥沙河道峰值流量系數(shù)0～21.00CH-COV泥沙河道覆蓋因子-0.001～11.00CH-EROD泥沙河道可侵蝕因子-0.05～0.600.13CMN水質(zhì)有機氮的礦化速度0.001～0.0030.002NPERCO水質(zhì)氮的下滲系數(shù)0～10.10

表3不同頻率年的計算結(jié)果

Table 3Calculation in different typical hydrologic years

頻率年降雨量/mm徑流量/(m3·s-1)泥沙量/108t污染物負荷量/t硝態(tài)氮有機氮氨氮亞硝態(tài)氮礦物磷有機磷豐水年631.86300.525.7718088.18130.871269.58856.47100.89903.16平水年573.21232.962.858181.2675.863020.291167.31383.98472.66枯水年469.83193.081.774681.6122.81448.55870.28256.2573.68

2.3渭河流域點源和非點源污染綜合措施模擬分析

隨著關(guān)中地區(qū)社會經(jīng)濟的快速發(fā)展,渭河沿岸的工業(yè)廢水和生活污水逐年增多。通過加大治理力度,大多數(shù)排污單位增加了污染治理設(shè)施,但仍有部分企業(yè)不能達標排放。加之工業(yè)廢水的達標排放標準遠高于地表水排放標準,所以即使工業(yè)廢水全部處理達標排放,對自然河道而言仍是污染濃度很高的廢水。2007年,渭河流域城鎮(zhèn)污水處理率仍不足40%,還有60%未經(jīng)處理的城鎮(zhèn)生活污水直接或間接排入渭河。由此可見,工業(yè)廢水和城鎮(zhèn)生活污水對渭河流域的水污染造成嚴重影響。

在SWAT模型模擬流域NSP的基礎(chǔ)上,加入點源數(shù)據(jù),經(jīng)率定、驗證后,對流域進行模擬。在流域水污染控制中主要表現(xiàn)為點源污染物進入河道的削減量,參考《渭河流域綜合治理五年規(guī)劃》和《渭河污染防治三年行動方案(2012—2014年)》,結(jié)合目前渭河流域水污染治理工作擬定及實施后對點源負荷的削減量,設(shè)定點源模擬情景,同時采取點源和面源污染的綜合控制措施進行模擬(表6)。根據(jù)所設(shè)定的控制措施對渭河流域2009年污染狀況進行模擬計算,結(jié)果見表7。

選取模擬結(jié)果中的NH4+-N、NO3--N、NO2--N、MINP、TN和TP作為研究對象,對平水年的模擬結(jié)果進行分析。方案1中,點源污染控制措施的實施對N、P的削減效果較明顯,且隨著強度加大,削減效果也明顯增強。情景7中,MINP削減率達29.22%,TP達35.79%，NH4+-N達41.22%，TN達46.05%。方案2模擬結(jié)果顯示:綜合措施對流域N、P均有較好的削減效果,情景4中MINP削減較明顯,削減率達50.77%，TP達63.17%，NH4+-N達70.84%，TN達71.60%。

可見,綜合措施對N、P的削減效果最好。但在方案設(shè)定時,點源的削減強度是按照比例設(shè)定的,考慮到實踐操作時存在的問題,認為點源削減方案中情景6的實踐操作性較大。在方案2中,情景4的削減效果最明顯,強度最大,但實施難度也大,其中灌溉用水量僅為原始值的10%,水土保持措施強度加大到0.25,化肥施肥量僅為原始值的15%。

表4渭河流域內(nèi)情景模擬方案

Table 4Scenario simulation of the Weihe River Valley

方案編號方案內(nèi)容情景設(shè)定說明 1采取水土保持措施,如帶狀種植、等高耕作或梯田等,增加土壤濕度、減少徑流量,以達到保持土壤的目的。改變水土保持因子USLE-P,不同USLE-P值對應(yīng)不同力度的水土保持措施情景1-1(原始狀況)USLE-P=1.0情景1-2(采取措施)USLE-P=0.8情景1-3(采取措施)USLE-P=0.6情景1-4(采取措施)USLE-P=0.5情景1-5(采取措施)USLE-P=0.4情景1-6(采取措施)USLE-P=0.252改進化肥施用方式,如條施、穴施、環(huán)施等,減少耕地中N、P的流失。FRT-LY1值為表層土壤施肥量占總施肥量的比例情景2-1(原始狀況)FRT-LY1=0.4情景2-2(采取措施)FRT-LY1=0.3情景2-3(采取措施)FRT-LY1=0.2情景2-4(采取措施)FRT-LY1=0.13增加農(nóng)家肥、堆肥、漚肥或泥肥的施用量,減少化學(xué)肥料的使用。采用不同的FRT-KG值對應(yīng)每年不同的化肥施用量情景3-1(原始狀況)FRT-KG=原始值情景3-2(采取措施)FRT-KG=原始值的75%情景3-3(采取措施)FRT-KG=原始值的65%情景3-4(采取措施)FRT-KG=原始值的50%情景3-5(采取措施)FRT-KG=04采取改進的灌溉方式以減少灌溉用水量,如噴灌和滴灌等。采用不同的IRR-AMT值對應(yīng)不同的灌溉用水量和灌溉方式情景4-1(原始狀況)IRR-AMT=原始值情景4-2(采取措施)IRR-AMT=原始值的80%情景4-3(采取措施)IRR-AMT=原始值的60%情景4-4(采取措施)IRR-AMT=原始值的50%情景4-5(采取措施)IRR-AMT=原始值的30%情景4-6(采取措施)IRR-AMT=05綜合方案1～4,同時采取水土保持,改變施肥方式、施肥量以及灌溉用水量等措施。同時改變USLE-P值、FRT-LY1值、FRT-KG值和IRR-AMT值,模擬非點源綜合管理措施情景情景5-1(原始狀況)USLE-P=1.0;FRT-LY1=0.4;FRT-KG=原始值;IRR-AMT=原始值情景5-2(采取措施)USLE-P=0.8;FRT-LY1=0.3;FRT-KG=原始值的75%;IRR-AMT=原始值的80%情景5-3(采取措施)USLE-P=0.6;FRT-LY1=0.2;FRT-KG=原始值的65%;IRR-AMT=原始值的60%情景5-4(采取措施)USLE-P=0.5;FRT-LY1=0.1;FRT-KG=原始值的50%;IRR-AMT=原始值的30%

USLE-P為水土保持因子；FRT-LY1為表層土壤施肥量占總施肥量的比例；FRT-KG為化肥施用量；IRR-AMT為灌溉用水量。

雖然近幾年的噴灌、滴灌、井渠雙灌等節(jié)水措施比較成熟,但由于渭河流域面積大，分布廣,故在現(xiàn)有水平下,若要削減90%的灌溉用水量就必須同時采取工程節(jié)水和非工程節(jié)水措施?？刂苹适┯昧靠梢詼p少污染負荷的產(chǎn)出,但在減少化肥施用量的同時要保證農(nóng)業(yè)經(jīng)濟效益,化肥量就不能減少太多。流域水土流失量大,是由于流域包括黃土高原,山地、坡耕地較多,土壤侵蝕嚴重;且支流粗沙區(qū)多,下泄粗泥沙會直接造成下游河道淤積,特別是在黃河洪水對渭河產(chǎn)生頂托倒灌時,北洛河高含沙洪水對下游淤積影響嚴重。故要治理渭河的水土流失需從溝道壩系建設(shè)、坡面治理、減少多沙粗沙區(qū)入渭泥沙量等方面入手?？傮w來看,水土流失治理進展緩慢，淤地壩建設(shè)難度大，邊治理邊破壞現(xiàn)象還相當嚴重,新的人為水土流失還在擴展,使本來就十分脆弱的生態(tài)環(huán)境更加惡化。因此，在近期要將渭河流域的水保措施強度加大到0.25的可操作性較低?？梢?情景4的控制管理措施實際可操作性不高,推薦情景3中的污染控制措施。

2.4預(yù)測渭河流域2020年污染控制措施效果

模擬流域綜合污染控制措施下2020年的污染負荷輸出,結(jié)果見表8。2020年污染負荷輸出較2009年有所下降,其中MINP減少1 516 t,TP減少1 953 t，NH4+-N減少14 141 t，NO2--N減少1 978 t，NO3--N減少4 649 t，TN減少34 122 t。渭河干流ρ(NH4+-N)平均值僅為1.273 mg·L-1,基本上達GB 3838─2012《地表水環(huán)境質(zhì)量標準》中Ⅲ～Ⅳ類水標準,滿足“十二五”規(guī)劃中主要污染物的控制目標。

表5各措施模擬非點源污染負荷的削減率

Table 5Pollution reducing rates of different control measures controlling simulated non-point source pollution load%

表6渭河流域內(nèi)綜合情景模擬方案

Table 6Program for comprehensive scenario simulation of the Weihe River Valley

方案情景說明10未采取任何點源控制措施1工業(yè)點源削減20%,生活點源不變2工業(yè)點源削減40%,生活點源不變3工業(yè)點源削減60%,生活點源不變4工業(yè)點源削減20%,同時生活點源削減20%5工業(yè)點源削減40%,同時生活點源削減40%6工業(yè)點源削減60%,同時生活點源削減60%7工業(yè)點源削減70%,同時生活點源削減70%21工業(yè)點源削減20%,生活點源削減20%;采取水土保持措施因子為USLE-P=0.8;化肥施用量減少25%,改變化肥施用方式,設(shè)定表層土壤施肥量占總施肥量的比例為FRT-LY1=0.3;改進灌溉方式,減少灌溉用水量為原始值的80%2工業(yè)點源削減40%,生活點源削減40%;采取水土保持措施因子為USLE-P=0.6;化肥施用量減少35%,改變化肥施用方式,設(shè)定表層土壤施肥量占總施肥量的比例為FRT-LY1=0.2;改進灌溉方式,減少灌溉用水量為原始值的60%3工業(yè)點源削減60%,生活點源削減60%;采取水土保持措施因子為USLE-P=0.5;化肥施用量減少50%,改變化肥施用方式,設(shè)定表層土壤施肥量占總施肥量的比例為FRT-LY1=0.1;改進灌溉方式,減少灌溉用水量為原始值的50%4工業(yè)點源削減70%,生活點源削減70%;采取水土保持措施因子為USLE-P=0.25;化肥施用量減少為原始值的15%,改變化肥施用方式,設(shè)定表層土壤施肥量占總施肥量的比例為FRT-LY1=0.1;改進灌溉方式,減少灌溉用水量為原始值的10%

表7模擬方案在2009年的污染變化

Table 7Change in pollution simulated for 2009

污染物方案1情景0情景1情景2情景3情景4情景5負荷輸出/t負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%礦物磷2852.002724.35-4.482497.43-12.432396.73-15.962639.26-7.462445.38-14.26總磷4097.003872.42-5.483473.16-15.233295.99-19.553722.69-9.143381.58-17.46氨氮25037.1423456.43-6.3120646.27-17.5419399.26-22.5222402.62-10.5220001.69-20.11亞硝態(tài)氮3108.002960.89-4.732699.35-13.152583.30-16.882862.81-7.892639.36-15.08硝態(tài)氮7975.617591.27-4.826908.00-13.396604.80-17.197335.04-8.036751.27-15.35總氮58996.6554835.54-7.0547438.00-19.5944155.35-25.1652061.46-11.7645741.19-22.47污染物方案1方案2情景6情景7情景1情景2情景3情景4變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t變化率/%負荷變化量/t礦物磷2284.97-19.882018.62-29.222443.82-14.312185.29-23.381923.02-32.571404.05-50.77總磷3099.35-24.352630.72-35.793361.29-17.962836.08-30.782292.37-44.051508.75-63.17氨氮18015.26-28.0514716.84-4.2219945.59-20.3416335.70-34.7512598.81-49.687301.80-70.84亞硝態(tài)氮2454.49-21.032147.51-30.902638.65-15.102307.19-25.771629.04-47.591062.80-65.80硝態(tài)氮6268.28-21.415466.29-31.466514.50-18.325813.64-27.114053.60-49.182896.75-63.68總氮40512.06-31.3331829.20-46.0546103.52-21.8539910.75-32.3529622.97-49.7916752.53-71.60

表8SWAT模型對2009和2020年渭河流域污染負荷的估算結(jié)果

Table 8Estimated pollution loads in the Weihe River Valley in 2009 and 2020 with the SWAT modelt

3結(jié)論

同時采取水保措施、合理減少施肥量和表層土壤施肥量占總施肥量的比例，并且改進灌溉方式和減少灌溉用水量,在流域中實行非點源綜合措施,對NSP負荷的削減效果顯著,且削減效果比任意一項單獨非點源控制措施明顯,在改善流域水環(huán)境質(zhì)量時可起到較大作用。綜合管理措施對流域N、P的削減效果比單一采取點源或非點源控制措施的效果明顯,且隨著措施強度的加大,污染負荷的輸出減少。

根據(jù)現(xiàn)有資料采取綜合控制措施預(yù)測2020年渭河流域污染負荷的產(chǎn)出情況:流域內(nèi)污染負荷均有所減少,且干流NH4+-N的平均濃度滿足“十二五”規(guī)劃中主要污染物總量的控制目標。這說明在2009年的基礎(chǔ)上,采取點源和非點源的綜合控制措施,對流域污染的治理可取得較好效果,可以有效防止渭河流域污染負荷的產(chǎn)出,提高水環(huán)境質(zhì)量。

參考文獻：

[1]宮瑩,阮曉紅,胡曉東.我國城市地表水環(huán)境非點源污染的研究進展[J].中國給水排水,2003,19(3):21-23.

[2]李家科,劉健,秦耀民,等.基于SWAT模型的渭河流域非點源氮污染分布式模擬[J].西安理工大學(xué)學(xué)報,2008,24(3):278-285.

[3]朱利,張萬昌.基于徑流模擬的漢江上游區(qū)水資源對氣候變化響應(yīng)的研究[J].資源科學(xué),2005,27(2):16-22.

[4]張雪松,郝芳華,楊志峰,等.基于SWAT模型的中尺度流域產(chǎn)流產(chǎn)沙模擬研究[J].水土保持研究,2003,10(4):38-42.

[6]杜娟,李懷恩,李家科.基于實測資料的輸出系數(shù)分析與陜西灃河流域非點源負荷來源探討[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2013,32(4):827-837.

(責任編輯： 陳昕)

收稿日期：2015-07-13

基金項目：國家水體污染控制與治理科技重大專項(2009ZX07212-002);榆林學(xué)院高層次人才啟動基金(14GK31)

中圖分類號：X52

文獻標志碼：A

文章編號：1673-4831(2016)04-0563-07

DOI：10.11934/j.issn.1673-4831.2016.04.008

作者簡介：杜娟(1987—),女,陜西榆林人,碩士生,研究方向為水資源保護。E-mail: dujuan3418@163.com

Simulation of Pollution Control Effect in the Weihe River Valley With SWAT Model.

DU Juan1, LI Huai-en2, ZHAO Xiang-bi1, LI Jia-ke2

(1.Department of Architectural Engineering, Yulin University, Yulin 719000, China; 2.Key Laboratory of Northwest Water Resource and Environment Ecology of Ministry of Education, Xi′an University of Technology, Xi′an 710048, China)

Abstract:Based on measured data, analog calculation of pollution load in the Weihe River Valley was done using the calibrated and validated SWAT model, and then based on the simulation, effects of the measures for control of point source (PSP) and non-point source pollution (NSP) in the valley were evaluated. Results show that comprehensive NSP control measures are better than single-factor NSP control measures in control effect. PSP control measures are quite good at reducing nitrogen and phosphorus pollution, and the effect increases with rising intensity of the measures. Comprehensive pollution control programs for control of PSP and NSP simultaneously are much better than any single-factor control measures in mitigating nitrogen and phosphorus pollution in the valley. Once the comprehensive pollution control programs are launched, it can be expected that the pollution load in the valley will have been reduced by 2020 to a certain extent, that will meet the pollution control targets specified in the “Twelfth Five-Year Plan” for major pollutants, which indicates that the comprehensive pollution control measures are very effective and the water environment quality in the valley will be remarkably improved.

Key words:Weihe River Valley; SWAT model; non-point source pollution; point source pollution; pollution control measure