[發(fā)明專(zhuān)利]基于GIS和CSLE的區(qū)域土壤侵蝕評(píng)價(jià)

在審

• 申請(qǐng)?zhí)枺?/td>: 202111148690.1
• 申請(qǐng)日: 2021-09-28
• 公開(kāi)（公告）號(hào): CN113706357A
• 公開(kāi)（公告）日: 2021-11-26
• 發(fā)明（設(shè)計(jì)）人: 武鵬;郭偉玲;管禮松;蔣婷婷;潘政;王艷麗;宋洪炳;沈?qū)毝?/a>
• 申請(qǐng)（專(zhuān)利權(quán)）人: 安徽理工大學(xué)
• 主分類(lèi)號(hào): G06Q50/26
• 分類(lèi)號(hào): G06Q50/26;G06F17/11;G06T17/05
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

1.一種基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，其具體過(guò)程為：

步驟1，將區(qū)域分成若干個(gè)柵格單元，計(jì)算每個(gè)單元的土壤侵蝕模數(shù)。

步驟2，將收集到的數(shù)據(jù)，用來(lái)提取土壤侵蝕各個(gè)因子值。

步驟3，將計(jì)算的各個(gè)因子?xùn)鸥駡D層相乘得到土壤侵蝕圖。

步驟4，用土壤侵蝕強(qiáng)度進(jìn)行區(qū)域土壤侵蝕評(píng)價(jià)。

2.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟2中利用日降雨數(shù)據(jù)計(jì)算降雨侵蝕力R值：

式中:M表示某半月時(shí)段的侵蝕力[(MJ·mm)/(h·hm²)]；k為該半月時(shí)段內(nèi)的時(shí)間(d)；P為半月時(shí)段內(nèi)第j天的侵蝕性日降雨量，閾值的選擇與我國(guó)侵蝕性降水標(biāo)準(zhǔn)一致，日雨量≥12mm；α和β為模型待定參數(shù)，與研究區(qū)降水特征有關(guān)，計(jì)算公式為：

β＝0.8363+18.177/P_d12+24.455/P_y12

α＝21.586·β^(-7.1891)

式中：P_d12為日降雨量≥12mm的日平均降雨量；P_y12為日雨量≥12mm的多年平均降雨量。

3.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟2中計(jì)算土壤可蝕性K值：

K＝{0.2+0.3exp[-0.0256SAN(1-SIL/100)]}*(SIL/(CLA+SIL))^0.3*[1.0-0.25C/(C+exp(3.72-2.95C))]*[1.0-0.7SN1/(SN1+exp(-5.51+22.9SN1))]

式中，CLA為粘粒、SAN為砂礫、SIL為粉粒；SN1＝1–SAN/100。

4.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟2中計(jì)算坡度坡長(zhǎng)因子LS：

式中：θ為坡度(°)；λ為坡長(zhǎng)(m)；m為可變的坡長(zhǎng)指數(shù)。

5.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟2中水土保持措施因子BET計(jì)算：

園地、林地和草地B因子計(jì)算公式：

式中，WR_i為第i個(gè)半月降雨侵蝕力占全年侵蝕力比例，取值范圍為0-1；SLR_i為第i個(gè)半月園地、林地和草地的土壤流失比例，無(wú)量綱，取值范圍為0-1，計(jì)算公式為：

茶園和灌木林地SLR_i計(jì)算公式：

果園、其他園地、有林地和其他林地SLR_i計(jì)算公式：

SLR_i＝0.44468×e(-3.20096×GD)-0.04099×e(FVC-FVC×GD)+0.025

草地SLR_i計(jì)算公式：

式中，F(xiàn)VC基于NDVI計(jì)算的植被覆蓋度，取值范圍為0-1；GD為喬木林的林下蓋度，取值范圍為0-1，包括除喬木林冠層以外的所以植被(灌木、草本和枯落物)構(gòu)成的林下蓋度，按實(shí)地調(diào)查或經(jīng)驗(yàn)取值。

工程措施因子E有兩種方法計(jì)算：一是根據(jù)公式進(jìn)行計(jì)算；二是根據(jù)不同土地利用進(jìn)行賦值。

方法1：

式中，St為梯田的面積，Sd為淤地壩的控制面積，S為土地的總面積，α、β分別代表梯田和淤地壩的減沙系數(shù)。

方法2：根據(jù)解譯獲取的土壤侵蝕地塊屬性表的“工程措施類(lèi)型或代碼”字段值，查水土保持工程措施因子賦值表，獲取工程措施因子值。

耕地措施因T:根據(jù)解譯獲得的耕作措施輪作區(qū)代碼，查表獲取耕作措施因子值生成T因子?xùn)鸥駡D層。

6.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟3中將計(jì)算的各個(gè)因子?xùn)鸥駡D層相乘得到土壤侵蝕圖。

7.根據(jù)權(quán)力要求1所述的基于GIS與CSLE的區(qū)域土壤侵蝕評(píng)價(jià)方法，其特征在于，在所述步驟4中將土壤侵蝕強(qiáng)度分為微度、輕度、中度、強(qiáng)烈、極強(qiáng)烈和劇烈，用以進(jìn)行區(qū)域土壤侵蝕強(qiáng)度評(píng)價(jià)。