土壤侵蚀的估算方法4830.pdf

《土壤侵蚀的估算方法4830.pdf》由会员分享，可在线阅读，更多相关《土壤侵蚀的估算方法4830.pdf（11页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、土 壤 侵 蚀 的 估 算 方 法 数 据 处 理 流 程 作者：*时间：2011 年 10 月 11 日 北京天合数维科技有限公司 目 录（CONTENT）一、所需数据与参数.错误!未定义书签。1、所需数据.错误!未定义书签。2、所需中间参数.错误!未定义书签。、水土保持因子 P.错误!未定义书签。、地标覆盖因子 C.错误!未定义书签。、地形因子 LS.错误!未定义书签。、土壤可视性因子 K.错误!未定义书签。、降水侵蚀因子 R.错误!未定义书签。3、所需参数.错误!未定义书签。、潜在土壤侵蚀量 Ap.错误!未定义书签。、现实土壤侵蚀量 Ar.错误!未定义书签。、土壤保持量 Ac.错误!未定

2、义书签。4、指标结果参数.错误!未定义书签。、保护土壤肥力的经济效益Ef.错误!未定义书签。、减少土地废弃的经济效益Es.错误!未定义书签。、减轻泥沙淤积的经济效益En.错误!未定义书签。二、处理流程.错误!未定义书签。1、DEM 数据的处理.错误!未定义书签。、坡长 L.错误!未定义书签。、百分比坡度 a.错误!未定义书签。、地形因子 LS.错误!未定义书签。2、气象数据.错误!未定义书签。、月降雨量 Pi 的计算.错误!未定义书签。、土壤侵蚀力指标 R.错误!未定义书签。3、土壤类型数据.错误!未定义书签。4、遥感影像数据.错误!未定义书签。5、土壤理性化数据.错误!未定义书签。三、所需参

3、数的计算.错误!未定义书签。四、指标结果参数计算.错误!未定义书签。一、所需数据与参数 在计算的过程中，总共涉及到的数据有地形数据、遥感影像数据、气象数据、土壤类型数据、土壤理性化数据以及统计数据，涉及到的中间参数有水土保持因子 P，地标覆盖因子 C，地形因子 LS，土壤可视性因子 K，降水侵蚀因子 R，所需要的参数有潜在土壤侵蚀量 Ap，现实土壤侵蚀量 Ar，土壤保持量 Ac，指标结果参数有保护土壤肥力的经济效益 Ef，减少土地废弃的经济效益 Es，减轻泥沙淤积的经济效益 En。1、所需数据 在进行土壤侵蚀的估算过程中，需要以下数据：A、地形数据；B、遥感影像数据；C、气象数据，主要是降雨量

4、数据；D、土壤类型数据；E、土壤理性化数据；F、统计数据。2、所需中间参数 在数据处理的过程中，所涉及到的中间参数与计算公式如下。、水土保持因子 P 按照游松财的方法，水田的 P 值取，其他土地利用方式基本没有采取水土保持措施，因此取值为。、地标覆盖因子 C 地表覆盖因子是根据地面植被覆盖状况不同而反映植被对土壤侵蚀影响的因素，与土地利用类型、覆盖度密切相关。C 值的估算采用如下公式：1,00.65080.3436lg,078.3%0,78.3%cfCccfCfCcfC 式中，fC 表示植被覆盖度，由遥感影像通过 NDVI 计算。、地形因子 LS 通过数字高程模型（DEM），计算获得坡长和坡度

5、，然后根据黄炎和等建立的方程式，获得 LS 的空间分布特征。6.035.008.0aLLS 式中，L，坡长（m）；a，百分比坡度。、土壤可视性因子K 采用陈明华等人建立的土壤可蚀性K 值的计算公式：K=10-3(160.80-2.31X 1+0.38X 2+2.26X 3+1.31X 4+14.67X 5)式中:K 可蚀性K 值(美国习用单位);X 1、X 2、X 3、X 4、X 5 分别表示细砾(1 3mm)%、细砂(0.05 0.25mm)%、粗粉粒(0.01 0.05 mm)%、细粉粒(0.005 0.01mm)%、有机质(10g/kg)。、降水侵蚀因子R 采用周伏建和黄炎和等人（199

6、7,2000,2002）根据实测数据提出的R 值计算式，该值考虑了月降水量。121(1.55270.1792)iRPi 式中，R 为降雨侵蚀力指标（），Pi 为月降雨量（mm）。月降水量Pi：根据气象站点观测的降水点位及数据，进行克里金插值获得。3、所需参数 在计算的过程中主要涉及到潜在土壤侵蚀量 Ap，现实土壤侵蚀量 Ar 以及土壤保持量 Ac 三个参数，可以通过以上 5 个中间参数来获得，其计算公式如下。、潜在土壤侵蚀量 Ap 潜在土壤侵蚀量不考虑地表覆盖类型和土地管理因素，即 c=1，p=1，此时，USLE 的形式为：pARKLS 式中，Ap 为潜在土壤侵蚀量（吨/平方千米），R 为降雨

7、侵蚀力指标，K 为土壤可侵蚀因子，LS 为坡长坡度因子。、现实土壤侵蚀量 Ar 现实土壤侵蚀量主要考虑了地表覆盖类型和土地管理因素，其计算式为：rARKLSCP 式中，Ar 为现实土壤侵蚀量（吨/平方千米），R 为降雨侵蚀力指标，K 为土壤可侵蚀因子，LS 为坡长坡度因子，C 为地表覆盖因子，P 为土壤保持措施因子。、土壤保持量 Ac 由上两式可得到土壤保持量：Ac=ApAr 式中，Ap为潜在土壤侵蚀量（吨/平方千米），Ar为现实土壤侵蚀量（吨/平方千米），Ac为土壤保持量（吨/平方千米）。4、指标结果参数 指标结果参数总共涉及到的数据有保护土壤肥力的经济效益Ef、减少土地废弃的经济效益Es和

8、减轻泥沙淤积的经济效益En三个参数，其计算公式与所需参数值见下。、保护土壤肥力的经济效益Ef 土壤侵蚀使大量土壤营养物质流失，主要表现在N、P、K 等营养物质的流失。()./10000(,)fceociiiEADCR CA C PiN P K 式中，Ef为保护土壤肥力的经济效益（元/a）；Ac，土壤保持量（t/hm2a）；De，土壤密度（t/m3）；Co,土壤有机质含量（kg/t）；R，薪柴转化为土壤有机质的比例，一般为 2；C，薪柴的平均价格（元/kg）；Ci，土壤中 N、P、K 的纯含量；Pi，N、P、K 的价格。注：N、P、K 的价格：参照国家林业局 2008 年 4 月 28 日发布的

9、中华人民共和国林业行业标准森林生态系统服务功能评估规范，氯化钾价格为 2200 元/吨，磷酸二铵价格为 2400 元/吨。元素名称 价格（元/吨）N P K 表 1、NPK 价格表、减少土地废弃的经济效益Es 根据土壤保持量和土壤表土平均厚度（）来推算因土壤侵蚀而造成的废弃土地面积，再根据机会成本法计算因土地废弃而失去的年经济价值。Es=AcPB10000 式中，Es，减少土地废弃的经济效益（元/a）；Ac，土壤保持量（t/a）；P，土壤的容重（t/m3）；B，湿地年均收益（元/hm2）。注：土壤容量取吨/立方米，湿地年均收益根据崔丽娟（2004 年）研究推算，湿地年均收益为元/平方米。、减轻

10、泥沙淤积的经济效益En 按照我国主要流域的泥沙运动规律，全国土壤侵蚀流失的泥沙有24%淤积于水库、江河、湖泊，这部分泥沙直接造成了水库江河、湖泊蓄水量的下降，在一定程度上增加了干旱、洪涝灾害发生的机会，因此可根据蓄水成本计算损失价值。En=Acp24%C 式中，En，减轻泥沙淤积的经济效益（元/a）；Ac，土壤保持量（t/a）；P，土壤的容重（t/m3）；C，水库工程费用（元/m3）。注：土壤容量取吨/立方米，核算过程中采用的水库清淤工程费用为元/立方米，此数据来源于硕士论文。二、处理流程 总体处理流程如下图所示：地形数据DEM遥感影像数据气象数据土壤类型数据土壤理性化数据统计数据坡长百分比坡

11、度地形因子（LS）归一化植被指数植被覆盖度（FC）地表覆盖因子（C）月降雨量Pi土壤侵蚀指标（R）水土保持措施因子（P）土壤可蚀性因子（K）NPK价格PiNPK纯含量Ci土壤容量P土壤密度De土壤有机质含量 Co湿地年均收入 B水库工程费用 C潜在土壤侵蚀量（Ap）现实土壤侵蚀量（Ar）土壤保持量（Ac）减少土地废弃的经济效益保护土壤肥力的经济效益减轻泥沙淤积的经济效益土 壤 侵 蚀 的 估 算方 法 图 1、数据处理流程图 1、DEM 数据的处理 运用 DEM 数据得到的最后结果是求出了地形因子 LS，在求 LS 的过程中，需要坡长 L 和百分比坡度 a，其求解主要在 Arcmap 中完成。

12、、坡长 L 坡长的求解主要参考汤国安老师的 ArcGIS 地理信息系统空间分析实验教程 中的第十一章中的水文分析来完成，Arcmap 中建立模型如下图所示：图 2、无洼地 DEM 生成模型 图 3、生成坡长 L 模型、百分比坡度 a 百分比坡度 a 是用 slope 来实现的。图 4、百分比坡度的生成、地形因子 LS 地形因子 LS 计算的模型为：图 5、地形因子 LS 的计算模型 2、气象数据 根据气象站点观测的降水点位及数据，进行克里金插值获得月降雨量数据 Pi，再根据 Pi 计算土壤侵蚀力指标 R。、月降雨量 Pi 的计算 、土壤侵蚀力指标 R 土壤侵蚀力指标 P 是在 Arcmap 下

13、 Spatial Analyst|Raster Caculator中计算而来。图 6、栅格计算 3、土壤类型数据 土壤可蚀性因子 K 是通过土壤类型数据直接对其进行赋值。4、遥感影像数据 通过遥感影像数据获得地表覆盖因子 C，首先求 NDVI，NDVI 的求解在 ENVI 中实现，主菜单中 transformndvi，算归一化植被指数。接着根据 NDVI 数据来求解植被覆盖度 FC，其计算公式为：FC=(NDVI-NDVI_MIN)/(NDVI_MAX-NDVI_MIN)式中：NDVI_MIN 为 NDVI 的最小值，NDVI_MAX 为 NDVI 的最大值。由此，我们可以根据公式 1,00.

14、65080.3436lg,078.3%0,78.3%cfCccfCfCcfC 来计算 C 值，在 Arcmap 中 Spatial Analyst|Raster Caculator实现。5、土壤理性化数据 按照游松财的方法，水田的 P 值取，其他土地利用方式基本没有采取水 土保持措施，因此取值为。三、所需参数的计算 所需要的参数有潜在土壤侵蚀量 Ap，现实土壤侵蚀量 Ar，土壤保持量 Ac，在进行了上述的计算以后，其计算比较简单，根据公式，在 Spatial Analyst|Raster Caculator中输入对应公示即可实现。最后对 Ac 进行汇总统计，在 Spatial Analyst|Zonal statistics 中进行。四、指标结果参数计算 这一步只是简单的数值计算，在此不再详述，但是在计算的过程中注意单位的统一。注：在进行统计Ac总和的时候，要乘以图像分辨率。还有单位一定要统一！