土壤侵蚀量估算

1、1 土壤侵蚀量计算模型关于土壤侵蚀量的计算，目前国内外主要采用的是美国的通用土壤流失方程 USLE（Universal Soil- Loss Equation）,作为一个经验统计模型，它是土壤侵蚀研究过程中的一 个伟大的里程碑，在土壤侵蚀研究领域一度占据主导地位，并深刻地影响了世界各地土壤侵蚀模型研究的方向和思路。由于USLE模型形式简单、所用资料广泛、考虑因素全面、因子具有物理意义，因此不仅在美国而且在全世界得到了广泛应用。“通用土壤流失方程式”的形式如下：A R K L S C P1-1式中：A土壤流失量（吨/公顷年）R降雨侵蚀力指标；K土壤可蚀性因子。它是反映土壤吝易遭受侵蚀程度的一个数

2、字。其单位是，在标准条件下，单位侵蚀力所产生的土壤流失量；L坡长因子。当其它条件相同时，实际坡长与标准小区坡长（22.1米）土壤流失量的比值；S坡度因子。当其它条件相同时，实际坡度与标准小区坡度（9%）上土壤流失量的比值；C作物经营因子。为土壤流失量与标准处理地块（经过犁翻而没有遮蔽的休 闲地）上土壤流失量之比值；P土壤保持措施因子，有土壤保持措施地块上的土壤流失量与没有土壤保持 措施小区（顺坡梨耕最陡的坡地）上土壤流失量之比值。通用土壤流失方程的计算结果只适用于多年平均土壤流失量，而不能够代表当地某一 年或某一次降雨所产生的土壤流失量。当方程式右边每个因子值都是已知数时，即地块内的土壤种类、

3、坡长、坡度、作物管理情况、地块内的土壤保持措施以及降雨侵蚀力都已知，且 都被分别赋于一个适当的数值时，它们相乘后，就得出在此特定条件下所预报的年平均土壤流失量。2模型中各参数确定依据降雨侵蚀力指标R值的确定R值的确定有以下三种途径：（1） R值的经典算法：美国学者威斯奇迈尔和史密斯（1985年）利用美国35个土壤保持试验站8250个休闲小区的降雨侵t资料统计得出R指标与降雨动能 E及最大30分钟降雨强度I30的经验关系，计算式如下:1-2R E?%0（2） R值的简易计算：上式在实际应用中，计算降雨动能E需要降雨过程，其计算是件繁杂的事情，故 R值简易计算的关键在于寻求一个通过常规降雨资料就可

4、得到的参数，并建立它与R值的经典算法的关系，省去动能E的计算。美国学者威斯奇迈尔提出了一个直接利用年平均降雨量和多年各月平均降雨量推求年降雨侵蚀力R值的经验计算公式：12P212（1.5 lg 0.8188）R 1.735 10 P1-31式中：R年降雨侵蚀力值;P年平均降雨量（mm ）;P1月平均降雨量（mm）。（3）根据各地的降雨资料获彳#全国降雨侵蚀力R值的等值线图（图1）,若无相关降雨资料，可查等值线图取得。（参考文献：王广忠、焦菊英，中国的土壤侵蚀因子定量评价研究，水土保持通报，）咫尔演1000土壤可蚀因子K值的确定USLE中的K值，是一个由实验确定的定量数值，对于一个具体土壤，它对

5、于“单元” 小区上每单位侵蚀指标的土壤流失率。单位小区的认为定义是：坡长 22.1m,均一坡度9%, 顺坡梨耕的连续休闲地。当这些条件都满足时，L, S, C, P各为，则K=A/EI。由于直接测定K值比较费事，因此美国在进行区域性土壤性质与K值关系研究的基础上，建立了通过土壤基本性质查算 K值的办法，并绘制了用以查算的诺谟图。我国K值的获取有两个方面：一是应用诺谟图；二是实测 K值。同时K值可以通过土壤质地和有机质含量从表1中查找。表1 USLE中的可蚀因子 K值土壤有机质含量（） 土壤质地24细沙壤质细沙土粘壤土沙质粘土粘土地形因子LS的确定L （坡长）和S （坡度）常作为一个复合因子进行

6、综合计算，LS就是特定坡度与坡长上的土壤流失量与标准小区上土壤流失量的比率。美国学者威斯奇迈尔和史密斯得出的适用于大于9%t度的LS关系为：LS （）03（产1-4225.16式中：丫一一坡地坡长（米）；0坡地坡度（度）。作物经营因子C与土壤保持措施因子P的确定C值确定：植被作用系数 C与植被类型、覆盖度有关，地面无任何覆盖时C值等于；有很好的植被或其它保护层时，C值等于。不同植被类型及覆盖率下的C值按表2求取。表2作物经营因子C值表不同植被覆盖的C值地面覆盖度（%）020406080100草地灌木乔灌混交茂密森林P值确定：P为特定水土保持措施的土壤流失与起伏耕地的相应土壤流失之比。措施好的P

7、值小，反之P值大。未采取任何措施的 P值为1。表3为我国不同措施的 P值。表3 土壤保持措施因子P值表中国不同措施P值坡度（0 ）等高带状耕作草田带状间作水平梯田水平沟等高垄作<55 10>103 土壤多年平均侵蚀量估算算例：利用上述计算模型，选取鄂东一安徽段几个典型水毁在灾害点，估算其多年平 均土壤侵蚀量，并对管道安全风险作评价。各水毁点的基本特点和坡体性质参数根据野外管道沿线地质灾害点的现场调查结果，选取灾害点编号为 74、132、137、141、147、149、182及184的水毁点进行土壤多年平均侵蚀量及侵蚀深度估算。各水毁点的基本特征如表4所示，管沟开挖后，坡体基本无任何

8、水工保护措施，坡面上已形成深度不等的冲沟。表5为各水毁灾害点坡体性质的基本参数。表4各水毁点基本特征表室内编号灾害类型附图水毁点基本特征74坡面水毁、 、奉2J塞魂i|LTk - - - t -不五工工科二；党-金电£-：管道向下敷设， 坡体上未修建水 工保护工程，坡 面冲刷成深沟， 目前水土流失严 重。光缆开挖形 成约2m宽的坑， 管道外露，未用土填上。132坡面水毁137坡面水毁141坡面水毁147坡面水毁管道顺斜坡敷设，坡度在15°左右。管沟为砂性碎石土，因管沟开挖致使土体松散，水土流失严重。管沟外侧已形成冲沟。管道上坡，坡体 为风化砂土，由于管沟开挖后，使得坡上土体

9、极 为松散，在坡面流水冲刷下，水 土流失严重，目 前坡面上已形成小冲沟。管道顺坡向上敷设，坡长约150m o上部局部 坡度达30°,整体 坡度为10-15°。 地层为薄层砂页岩，强-全风化， 极易发生冲刷。管道上坡，坡体 为风化砂土，由于管沟开挖后，使得坡上土体极 为松散，在坡面流水冲刷下，水 土流失严重，目 前坡面上已形成冲沟。149坡面水毁 管道上坡，坡体 为风化砂土，由 于管沟开挖后， 使得坡上土体极 为松散，在坡面 流水冲刷下，水 土流失严重。-*1j艮* 二，/：，T 了山 ./ ?*0*5*182坡面水毁, /，3L、zSr.z'-二,”j.二 :/ 髭

10、-r、.A.管道顺坡向下敷 设，由于管沟开 挖后，使得坡上 土体极为松散， 在坡面流水冲刷 下，水土流失严 重，目前坡面上 已形成冲沟。184坡面水毁：*1iJKJL一，1A Tdil管道上方土体松 散，表面水流冲 刷下，水土流失 严重，目前坡面 上已形成冲沟。蓍森 * 一1，:表5各水毁灾害点坡性参数一览表灾害点编号地区土体类型坡长平均坡度相对高度植被覆盖率水土保持措施74大冶市砂性碎石土8015荒草地，小于5%无132池州市砂性碎石土8015灌木，30%无137池州市风化砂土6015灌木，30%无141池州市砂性碎石土15013灌木，30%无147池州市风化砂土12015荒草地，10%无1

11、49池州市风化砂土6020荒草地，5%无182广德县风化砂土6015灌木，10%无184广德县风化砂土4520灌木，5%无各水毁点计算参数的确定土壤多年平均侵蚀量的计算公式为：A R K L S C P （R取多年平均值）依据上述计算模型各参数白确定原则，得到表6所示的各灾害点的计算参数取值表，74号灾害点位于大冶市，其余各灾害点均位于池州市，由于无实际降雨资料可查，故依据前人所作的全国降雨侵蚀力 R值的等值线图进行估算，74号灾害点的年均 R值取，其他灾害点的年均R值取，地形参数LS根据公式1-4确定，C、P参数分别依表2和表3可查的表6各水毁点灾害点计算参数一览表灾害点编号RKLSC&qu

12、ot;CP"P741113211137111411114711149111821118411注：表中C"、P"和C、P分别为管沟开挖前后的作物经营因子和水土保持因子值。表7各个水毁灾害点多年平均侵蚀量估算一览表灾害点编号A"（吨/公顷/年）A（吨/公顷/年）H"（米/年）H （米/年）Y"（年）丫（年）74132137141147149182184注：表中A"、H"和A、H分别为管沟开挖前后的土壤多年平均侵蚀量和侵蚀深度；管道平均埋深取，Y"、Y为土壤侵蚀深度达所需的年数。各水毁点土壤多年平均侵蚀量、侵蚀深度估算及管道安全风险评价各参数确定以后即可估算出各灾害点在管沟开挖前后的土壤多年平均侵蚀量，计算结果如表7所示。由表中计算结果可知：(1)管沟开挖后，土壤多年平均侵蚀量、侵蚀深度显著加大，水土流失现象严重。74号灾害点，土壤多