基于WEPP模型在烏裕爾河流域的應(yīng)用

■程雪 陳圣波 周歡

（吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院 吉林長(zhǎng)春130026）

基于WEPP模型在烏裕爾河流域的應(yīng)用

■程雪 陳圣波 周歡

（吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院 吉林長(zhǎng)春130026）

本文以小興安嶺西南坡的烏裕爾河流域的氣候數(shù)據(jù)資料為基礎(chǔ)，利用WEPP模型模擬分析坡度因子與土壤侵蝕量的關(guān)系。結(jié)果表明在一定的坡度范圍內(nèi)，隨著坡度的增加，土壤侵蝕量與坡度呈冪函數(shù)遞增關(guān)系，原因主要有坡度的陡緩決定了水力阻力的大小和坡面承雨面積的改變等。表明用WEPP模型來(lái)模擬烏裕爾河流域坡面與土壤侵蝕量的關(guān)系是可行的。對(duì)今后建立適用于我國(guó)東北地區(qū)侵蝕環(huán)境的侵蝕預(yù)報(bào)模型具有指導(dǎo)意義。

坡度WEPP模型土壤侵蝕烏裕爾河流域

1　引言

地面坡度是地形因素中對(duì)坡面土壤侵蝕起重要作用的因子。本文研究WEPP模型（坡面版）在小興安嶺西南坡的烏裕爾河流域不同坡度條件下的土壤侵蝕量，為今后研究烏裕爾河流域的土壤侵蝕預(yù)報(bào)提供技術(shù)支持。

圖1　研究區(qū)概況Figure 1 general situation of the Study area

土壤侵蝕預(yù)報(bào)模型是國(guó)土資源調(diào)查、水土保持效益評(píng)價(jià)和水土保持規(guī)劃設(shè)計(jì)的重要手段，受到世界各國(guó)的普遍重視。新一代水蝕預(yù)報(bào)模型——WEPP(Water Erosion Prediction Project)[1]是美國(guó)農(nóng)業(yè)部農(nóng)業(yè)研究局 (USDA-ARS)、土壤保持局 (USDA-SCS)、林業(yè)局(USDA-FS)、普渡大學(xué)(Purdue University)等共同開(kāi)發(fā)的，其目的是替代通用土壤流失方程—USLE(Universal Soil Loss Equation)[2]，推出更能滿足實(shí)際生產(chǎn)要求的土壤侵蝕預(yù)報(bào)模型。WEPP以物理過(guò)程為基礎(chǔ)的原理，同傳統(tǒng)的經(jīng)驗(yàn)?zāi)Ｐ拖啾雀欣谠诿绹?guó)以外其他地區(qū)的推廣和應(yīng)用。目前，我國(guó)學(xué)者對(duì)WEPP模型在我國(guó)紫色土丘陵區(qū)和黃土丘陵溝壑區(qū)的適用性進(jìn)行了研究，發(fā)現(xiàn)在一定坡度條件下WEPP模型對(duì)侵蝕過(guò)程模擬較為合理[3]，同時(shí)WEPP模型的模擬結(jié)果優(yōu)于USLE模型[4]。但WEPP模型在東北流域的研究還相對(duì)較少。

本文基于克山氣象站1982—2013年30年的氣象資料數(shù)據(jù)，利用WEPP模型模擬在不同坡度條件下的多年平均下的土壤侵蝕量，研究并分析小興安嶺西南坡烏裕爾河流域的坡度與土壤侵蝕量的關(guān)系。從而為建立適用于我國(guó)侵蝕環(huán)境的侵蝕預(yù)報(bào)模型提供科學(xué)依據(jù)。

2　研究區(qū)概況與數(shù)據(jù)收集

2.1研究區(qū)概況

以小興安嶺向東北平原過(guò)渡的烏裕爾河流域?yàn)檠芯繀^(qū)，地貌上從平原向山地過(guò)渡，位于東北平原向小興安嶺的過(guò)渡地帶。烏裕爾河位于黑龍江省西部齊齊哈爾市境內(nèi)，為省內(nèi)最大的內(nèi)陸河。流域面積23110平方公里。

2.2數(shù)據(jù)收集

研究區(qū)氣候資料來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)，年限為1982—2013年。地形數(shù)據(jù)采用烏裕爾河地區(qū)的ASTER GDEM數(shù)據(jù)，土地利用類型數(shù)據(jù)來(lái)源于全球土地覆蓋GlobCover2009數(shù)據(jù)。

3　 WEPP模擬方法

3.1技術(shù)路線

WEPP模型的運(yùn)行需要建立氣候、土壤、坡面和作物管理4個(gè)數(shù)據(jù)庫(kù)。WEPP模型的氣候數(shù)據(jù)可以由用戶直接輸入或由模型生成。由于從中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)下載的烏裕爾河流域氣象數(shù)據(jù)的氣象參數(shù)只包括WEEP模型需要的日降雨量、日最高溫和日最低溫，所以采用CLIGEN模型來(lái)模擬生成氣候數(shù)據(jù)庫(kù)。根據(jù)已知?dú)庀笳镜臍庀髷?shù)據(jù)，以年降水量的±5%為范圍選擇年降水量比較接近的美國(guó)氣象站，克山站年降水量的選擇范圍為549.2235～607.0365mm，選擇所有年降水量在上述范圍內(nèi)的美國(guó)氣象站，以研究站與參證站間月降水決定系數(shù)（R2）最大為標(biāo)準(zhǔn)選擇各站的參證站[5]，克山站的參證站為Minnesota的GUNFLINT LODGE MN氣象站，研究站與參證站月降水量間的決定系數(shù)為0.9095。

圖2　技術(shù)流程圖Figure 2technical flowchart

3.2WEPP模型坡面版基本原理

WEPP模型坡面版是WEPP模型中最簡(jiǎn)單、最基本的模型版本，將坡面侵蝕分為細(xì)溝侵蝕和細(xì)溝間侵蝕。WEPP模型[6]利用穩(wěn)態(tài)泥沙連續(xù)方程來(lái)描述泥沙的運(yùn)動(dòng)：

式（1）中，x：某點(diǎn)沿下坡方向的距離（m），G：輸沙量（kg/s·m），Dr：細(xì)溝侵蝕速率（kg/s·m2），Di：細(xì)溝間泥沙輸移到細(xì)溝的速率（kg/s·m2）

①水流剪切力（τf）大于臨界土壤剪切力（τc），并且輸沙量（G）小于泥沙輸移能力（Tc）時(shí)，細(xì)溝內(nèi)以搬運(yùn)過(guò)程為主：

式中，Dc：細(xì)溝水流的剝離能力（kg/s·m2），Tc：細(xì)溝間泥沙輸移能力（kg/s·m）。Kr：細(xì)溝可蝕性參數(shù)（s/m），τf：水流剪切壓力(pa)，τc：臨界剪切力(pa)。

②輸沙量（G）大于泥沙輸移能力（T）時(shí)，以沉積過(guò)程為主：

式（4）中，Vf：有效沉積速率（m/s），q：?jiǎn)螌捤髁髁浚╩2/s），β：雨滴擾動(dòng)系數(shù)。

3.3WEPP模型結(jié)構(gòu)

運(yùn)行WEPP模型（坡面版）需要建立氣候、土壤、坡度坡長(zhǎng)和作物管理4個(gè)數(shù)據(jù)庫(kù)（見(jiàn)表1），每個(gè)數(shù)據(jù)庫(kù)都涉及多個(gè)參數(shù)[7]。

表1　WEPP模型（坡面板）輸入?yún)?shù)Table 1 Input parameters of WEPP model(slope panel)

4　結(jié)果分析

4.1土壤侵蝕量隨坡度變化的關(guān)系

利用建立的氣候、土壤、坡面和作物管理4個(gè)數(shù)據(jù)庫(kù)，利用WEPP模型分別模擬坡度為2°、3°、4°、5°、6°、7°、8°、9°、10°、11°、12°十一個(gè)坡度條件下的多年平均土壤侵蝕量。

圖3　土壤侵蝕量與坡度的變化關(guān)系Figure 3therelationship of soil erosion and slope changes

根據(jù)WEPP模型輸入不同坡度參數(shù)，可以得到坡度與侵蝕量的關(guān)系，土壤侵蝕量與坡度的關(guān)系指數(shù)為1.015，從土壤侵蝕量的變化趨勢(shì)來(lái)看，模型模擬的土壤侵蝕量值均隨坡度的增加而呈增加趨勢(shì)。在一定的坡度范圍內(nèi)，輸入不同坡度參數(shù)得到的結(jié)果是隨著坡度的增加，土壤侵蝕量也隨之增加，但是坡度每增加一度，土壤侵蝕量增加的幅度有所不同(如圖4)。

圖4　侵蝕增加量隨坡度的變化關(guān)系Figure 4 erosion increases along with the change of slope

從圖4中可以看出，坡度在0°～7°范圍內(nèi)，土壤侵蝕量隨坡度變陡而增加，至12°時(shí)隨坡度的增加，其增加率減緩。

4.2影響因素分析

坡度是影響坡面徑流侵蝕能力的動(dòng)力因子，坡度的陡緩決定了坡面徑流動(dòng)力比降的大小[8]。坡面流受降雨的影響，直接表現(xiàn)在流速的變化上，在一定的坡度范圍內(nèi)，隨著坡度的增加，根據(jù)水力學(xué)原理[9]，水流在重力方向上的分力增加，流速呈增加趨勢(shì)。

土壤侵蝕過(guò)程是區(qū)域地理環(huán)境各因素相互作用的過(guò)程，侵蝕狀況取決于侵蝕動(dòng)力和抗侵蝕阻力的對(duì)比關(guān)系。在一定的水深條件下，坡度越陡，其水力比降就越大，陡坡水流對(duì)地表就具有更大的切應(yīng)力，也就具有更大的分散能力。同時(shí)坡度越大，土壤顆粒在順坡方向的分力也愈大，土粒穩(wěn)定性降低，更容易遭受侵蝕。

坡面土壤侵蝕率除了水力阻力隨坡度變化的規(guī)律之外，還與坡面承雨面積有關(guān)[10]，因?yàn)閷?duì)于一定坡長(zhǎng)的坡面，坡度越大，其承雨面積就越小，形成的徑流量就相對(duì)減少，即在一定坡度上符合坡度越陡，一定坡長(zhǎng)上的坡面徑流量就越少。

5　結(jié)論

通過(guò)1982-2013年的烏裕爾河地面氣候數(shù)據(jù)，利用WEPP模型在不同坡度條件下模擬坡度與土壤侵蝕量的關(guān)系是可行的，得出了在一定的坡度范圍內(nèi)，隨著坡度的增加，土壤侵蝕量與坡度呈冪函數(shù)遞增關(guān)系。為模型模擬其它地區(qū)的應(yīng)用提供一個(gè)參考。

WEPP模型是“新一代水蝕預(yù)報(bào)技術(shù)”開(kāi)發(fā)的計(jì)算機(jī)土壤侵蝕預(yù)報(bào)模型，是一個(gè)過(guò)程模型，可應(yīng)用于更廣泛的條件下，不僅可以預(yù)測(cè)產(chǎn)沙量、輸沙量、侵蝕強(qiáng)度等，而且可以估算土壤侵蝕時(shí)空分布即全坡面或坡面任一點(diǎn)的凈土壤流失量及其隨時(shí)間的變化，對(duì)侵蝕過(guò)程的每一個(gè)階段都可以進(jìn)行模擬和預(yù)測(cè)[11]。還可以對(duì)影響因子和侵蝕過(guò)程進(jìn)行分析和描述，只是模型預(yù)測(cè)的準(zhǔn)確性是以大量的參數(shù)為前提的，這就要求我們做好基礎(chǔ)工作，收集大量的基礎(chǔ)數(shù)據(jù)資料，為充分利用WEPP帶給我們的好處，建立更加合理有效的土壤侵蝕預(yù)測(cè)模型。

[1]張玉斌，鄭粉莉，賈媛媛.WEPP模型概述 [J].水土保持研究，2004，11(4)：146-149.

ZHANG Yu-bin，ZHENG Fen-bin，JIA Yuan-yuan.WEPP Model and Its Application[J]. Research of Soil and Water Conservation，2004，11(4)：146-149.

[2]鄭粉莉，劉峰，楊勤科，等.土壤侵蝕預(yù)報(bào)模型研究進(jìn)展 [J].水土保持通報(bào)，2001，21(6)：16-18,32.

ZHENG Fen-fen，LIU Feng，YANG Qin-ke，et al.Review ofResearch Progress in SoilErosion PredictionM odel[J].Bulletin of Soil and Water Conservation，2001，11(4)：146-149.

[3]繆馳遠(yuǎn)，何丙輝，陳曉燕，等.WEPP模型中的CLIGEN與BPCDG應(yīng)用對(duì)比研究[J].中國(guó)農(nóng)學(xué)通報(bào)，2004，20(6)：321-324.

LIAO Chi-yuan，HE Bing-hui，CHEN Xiao-yan，et al.Study on the Application and Contrast of CLIGEN and BPCDG in WEPP Model Application[J].Chinese Agricultural Science Bulletin，2004，20(6)：321-324.

[4繆馳遠(yuǎn)，何丙輝，陳曉燕.水蝕模型USLE與WEPP在紫色土水蝕預(yù)測(cè)中的應(yīng)用對(duì)比研究 [J].農(nóng)業(yè)工程學(xué)報(bào).2005，1(1)：13-16.

LIAO Chi-yuan，HE Bing-hui，CHEN Xiao-yan.Application and comparison of the erosion model USLE and WEPP in predicting erosion of purple soil[J].Transactions of the CSAE，2005，1(1)：13-16.

[5]項(xiàng)甜甜.大興安嶺東坡侵蝕-沉積過(guò)程及其指示意義研究：碩士學(xué)位論文 [D].長(zhǎng)春：吉林大學(xué)，2014.

XIANG Tian-tian.Studying on Process of Erosion-deposition on Eastern Slope of Great Khingan and ItsIndicative Significance：master's degree thesis[D].Changchun：Jilin University，2014.

[6]劉益軍，王昭艷.水力侵蝕預(yù)報(bào)模型 WEPP氣象數(shù)據(jù)的輸入方法 [J].水土保持科技情報(bào)，2003，(4)：4-5.

LIU Yi-jun，WANG Zhao-yan.The input method of Water erosion prediction model WEPP meteorological data[J].Scientific and Technical Information of Soil and Water Conservation，2003， (4)：4-5.

[7]王建勛.WEPP模型（坡面版）在黃土高原丘陵溝壑區(qū)的適用性評(píng)價(jià)：碩士學(xué)位論文 [D].西安：西北農(nóng)林科技大學(xué)，2007.

WANG Jian-xun.Assessment on WEPP Model(Hillslope Version)Application to Hilly-Gully Region on the Loess Plateau：master's degree thesis[D].Xian：Northwest A&F University，2007.

[8]蔣定生.地面坡度對(duì)降雨入滲影響的模擬試驗(yàn) [J].水土保持通報(bào)，1984，5(4)：83-84.

JIANG Ding-sheng.A simulation test with the impact of the ground slope to rainfall infiltration [J].Bulletin of Soil and Water Conservation，1984，5(4)：83-84.

[9]芮孝芳.水文學(xué)原理 [M].北京：中國(guó)水利水電出版社，1997.

RUI Xiao-fang.Principles of Hydrology[M].Beijing：China Water&Power Press，1997.

[10]曲遠(yuǎn)馨，孫國(guó)勝，許中杰，等.青海賽多浦崗日地區(qū)雀莫錯(cuò)組粒度特征及沉積環(huán)境[J].世界地 質(zhì)，2015，34(2)：437-443.

QU Yuan-xin，SUN Guo-sheng，XU Zhong-jie，et al.Characteristic of granularity of Quemocuo Formation in Saiduopugangri region of Qinghai and its sedimentary environment[J]. Global Geology，2015，34(2)：437-443.

[11]李明飛，葉松青，楊言辰，等.黑龍江老柞山金礦床地質(zhì)、地球化學(xué)特征及其成礦構(gòu)造背景[J].世界地質(zhì)，2014，33(3)：543-553.

LI Ming-fei，YE Qing-song，YANG Yan-chen，et al.Geological and geochemical characteristics of Laozuoshan gold depositin Heilongjiang and its metallotectonic setting[J].Global Geology，2014，33(3)：543-553.

A research of the relationship between slope factor and soil erosion based on the WEPP model——A Case Study of WuyuerRiver

CHENG Xue，CHEN Shen-bo，ZHOU Huan

College of GeoExploration Science and Technology,JilinUniversity, Changchun130026,China

The study used WEPP model to simulation and analysis the relationship between slope factor and soil erosion based on the climate data of WuyuerRiver in the southwest of Lesser Khingan Mountains.Theresult showed that in a certain slope range,with the increase of slope, slope and soil erosion increased in power function relationship.The reasons were that hydraulic friction and change of the acceptant area were decided by steep and slow of the slope.The result also showed that the relationship simulated by WEPP model was feasible.It has significance on establishing erosion forecasting model that suitable for China's northeast erosion environment area in the future.

slope;WEPP model;soil erosion;WuyuerRiver Region

P[文獻(xiàn)碼]B

1000-405X（2016）-2-402-3

程雪（1990～），女，碩士研究生，研究方向?yàn)榇髿膺b感。

陳圣波（1967～），男，教授，博士生導(dǎo)師，研究方向?yàn)榈刭|(zhì)遙感、月球探測(cè)等。