简介SL773-2018《生产建设项目土壤流失量测算导则》是水利行业现行标准面向水土保持方案编制人员、水利工程设计者及环境评价从业者用于规范生产建设项目土壤流失量的测算方法与流程。资源包内含1个PDF文件大小约27.85MB完整收录标准正文及附录涵盖土壤流失类型划分、扰动单元与计算单元确定、水力与风力作用下的测算公式以及工程开挖面、堆积体等典型场景的取值表。附录还提供全国各县级行政单元多年平均逐月和年降雨侵蚀力因子、土壤可蚀性因子及风蚀率参考值便于直接查用。目前已有4192人学习下载适合需要对照标准开展水土保持监测、方案编制与预测估算的技术人员可快速掌握测算流程与参数选取依据。1. 从一份被翻烂的 PDF 说起SL 773-2018 到底解决什么问题做水土保持方案的人大概都有过这种经历项目区在山区开挖面、堆积体、扰动地表混在一起专家评审会上被问一句“你这土壤流失量怎么算出来的”如果只能答“按经验估的”基本就凉了。SL 773-2018《生产建设项目土壤流失量测算导则》就是冲着这个场景来的——它把生产建设项目扰动后的土壤流失从“拍脑袋”拉回到“按公式、按参数、按流程”的轨道上。这份标准 2018 年 10 月 23 日发布2019 年 1 月 23 日实施替代了 SL 41-2011 等旧标准中关于土壤流失量测算的内容。它适用于水力和风力作用下生产建设项目土壤流失量的事前预测、事中监测和事后计算核心是把扰动地表分成一般扰动地表、工程开挖面、工程堆积体三大类再按上方有无来水、植被破坏型还是地表翻扰型细分每一类对应一套计算公式。适合谁水土保持方案编制人员、监测单位技术员、水利水电和线性工程的环境工程师以及需要复核土壤流失量的评审专家。如果你手上只有一份 PDF却不知道怎么把里面的公式落到具体项目上这篇笔记就是按“能照着算”的路子拆的。2. 测算流程拆解从扰动单元到计算单元的落地路径2.1 扰动单元划分先搞清楚“哪些地方算一类”标准第 5 章把流程拆成五步划分扰动单元、确定典型扰动单元、现场查勘、划分计算单元、土壤流失量计算。第一步最容易翻车因为很多人直接按施工分区来划结果把不同降雨区、不同土壤质地、不同侵蚀外营力的地块混在一起。扰动单元的划分原则是扰动方式相同、扰动强度相仿、土壤类型和质地相近、气象条件相似、空间上相连续。具体操作时标准给了几条硬杠杠空间不连续的区域必须分开年降雨量按 800mm、800~400mm、400~200mm、≤200mm 四个区间分砂土、壤土、黏土等不同质地分开水力作用和风力作用分开不同下垫面类型一般扰动地表、开挖面、堆积体分开。我一般会先拿项目总平面布置图和施工进度表在 GIS 里叠一下降雨等值线和土壤类型图把明显跨区的块先切出来。这一步偷懒后面计算单元划分和参数取值都会跟着乱。2.2 典型扰动单元抽样10% 和 20 个的底线标准要求按不同类型下的大、中、小规模分别抽样抽样数量不少于同一类型下同等规模扰动单元数量的 10%涉及的各种类型和规模至少选 1 个抽样总数不少于 20 个。实际操作中如果项目扰动单元总数不到 20 个那就全部作为典型扰动单元处理。这里有个常见误解以为抽样就是随便选几个。其实抽样要覆盖“类型 × 规模”的组合比如“水力作用 一般扰动地表 植被破坏型 大规模”至少得有一个代表。如果项目里某种组合确实不存在那就不用选但要在报告里说明。2.3 现场查勘参数怎么测、测什么现场查勘要测的参数包括长度、宽度、坡度植被类型、郁闭度、覆盖度、砾石盖度物质组成和性状上方汇水面积气象资料。标准对测量工具也给了建议长度宽度用皮尺或测距仪坡度用罗盘或坡度测定仪植被郁闭度和覆盖度用照相法或目估法。土壤质地和机械组成有实验条件就现场取样回实验室测没条件就按附录 B 的现场干测法确定。这里提醒一句附录 B 的表 B.1 和表 B.2 是现场鉴别用的别拿实验室数据去套现场干测的等级两者口径不一样。气象参数获取优先级是典型扰动单元附近气象设备监测数据 最近气象站数据 所在县域气象数据。如果连县域数据都拿不到才用附录 C 的县级行政区参考值。2.4 计算单元划分五个“不同”必须分开计算单元是在典型扰动单元内部进一步细分划分依据是坡度不同相差大于 5° 的分开坡长不同相差大于 5m 的分开土壤类型、质地、含水量、有机质含量、砾石含量不同的分开植被类型、郁闭度、覆盖度、砾石盖度不同的分开水土保持措施不同的分开有阻风设施且类型一致的合在一起否则分开上方有来水冲刷和无来水冲刷的分开不同县域或不同气象站控制范围的分开。划分完计算单元后如果现场查勘有漏测的参数必须补测。最后统计每个计算单元的长度、宽度、坡度、面积、土壤参数、植被参数、水保措施、上方汇水面积、降水、风速风向等信息形成一张计算单元要素表。2.5 土壤流失量计算分类型选公式、分时段累加计算时先根据计算单元所属扰动类型按表 1 选对应公式。时间尺度可以按日、月、年、多年平均优先选用过程资料详细的公式。同一计算单元的水力侵蚀和风力侵蚀要分别计算不能混在一起。统计逻辑是计算单元流失量之和 典型扰动单元流失量典型扰动单元流失量之和 项目总流失量扰动单元 ≤20 个时。如果扰动单元 20 个先算同一类型下同等规模典型扰动单元的平均值再乘以该类型和规模的扰动单元数量各乘积求和。3. 水力作用下三大下垫面测算公式、参数与代码实现3.1 植被破坏型一般扰动地表R、K、L、S、B、E、T 七因子植被破坏型一般扰动地表的公式是My R × K × Ly × Sy × B × E × T × A其中 My 是土壤流失量tR 是降雨侵蚀力因子K 是土壤可蚀性因子Ly 是坡长因子Sy 是坡度因子B 是植被覆盖因子E 是工程措施因子T 是耕作措施因子A 是计算单元水平投影面积hm²。降雨侵蚀力因子 R 按资料条件分四种算法# 日降雨资料逐日计算后累加 def R_daily(P_r, P_d12, P_y12): P_r: 第r日降雨量mm需 12mm否则按0计 P_d12: 日降雨量12mm的日均雨量mm P_y12: 日降雨量12mm的年均雨量mm if P_r 12: return 0.0 alpha 21.586 * (P_d12 ** -7.182) beta 0.836 17.144 * (P_d12 ** -1) 24.455 * (P_y12 ** -1) return alpha * (P_r ** beta) # 月降雨资料 def R_monthly(P_m): P_m: 第m月降雨量mm return 0.183 * (P_m ** 1.996) # 年降雨资料 def R_annual(P_n): P_n: 年降雨量mm return 0.053 * (P_n ** 1.655) # 多年平均降雨资料 def R_multi_year(P_d): P_d: 多年平均降雨量mm return 0.067 * (P_d ** 1.627)参数说明日降雨公式里 α 和 β 是统计系数依赖日均雨量和年均雨量这两个值要用项目区长期降雨序列算不能随便拿一年数据凑。月、年、多年平均公式是简化式精度依次降低但资料要求也依次降低。土壤可蚀性因子 K 优先用标准小区实测资料算要求土壤类型相同、实测序列不少于 3 年。没有实测资料时能测粒径就用公式9算不能测就查附录 C。坡长因子 Ly 和坡度因子 Sy 的计算import math def L_y(lam, theta): lam: 计算单元水平投影坡长m100m按100m计 theta: 坡度度 lam min(lam, 100.0) theta_rad math.radians(theta) lam_x lam / math.cos(theta_rad) if theta 1: m 0.2 elif theta 3: m 0.3 elif theta 5: m 0.4 else: m 0.5 return (lam / 20.0) ** m def S_y(theta): theta: 坡度度35度按35度计0度取0 if theta 0: return 0.0 theta min(theta, 35.0) theta_rad math.radians(theta) return -1.5 17.0 / (1.0 math.exp(2.3 - 6.1 * math.sin(theta_rad)))注意用 Office 或 WPS 算三角函数时角度要先转弧度这是血泪经验直接写角度进去结果会差一个数量级。植被覆盖因子 B 按草地、灌木林地、乔木林地分别查表 4 和表 5灌草混合按灌木林地对待。工程措施因子 E 和耕作措施因子 T 查表 6 和表 7、表 8。3.2 地表翻扰型一般扰动地表和植被破坏型的差异地表翻扰型的公式结构和植被破坏型一样区别在于扰动程度不同植被覆盖因子 B 的取值逻辑有变化。标准里地表翻扰型的定义是“人为活动导致地表土壤翻动原有植被覆盖明显减少或裸露维持原有整体地形的扰动地表”。实际计算时B 值要根据翻扰后的植被恢复情况取不能直接套用翻扰前的覆盖度。我一般会按施工进度分时段翻扰刚完成时 B 取接近 1 的值随着植被恢复逐步降低。如果项目有分季度或分月的植被调查数据就按实测覆盖度查表没有的话按经验在翻扰后第一年取 0.6~0.8第二年取 0.3~0.5具体看恢复措施。3.3 工程开挖面上方无来水与有来水的分界工程开挖面的标准下垫面是宽 1m、水平坡长 5m、坡度 50° 的全土质裸露面。上方无来水时公式里不考虑径流冲刷力因子上方有来水时要加上上方汇水带来的冲刷侵蚀。判断“上方有无来水”的关键是开挖面上缘是否达到或越过分水岭顶部是否有截排水沟等坡面径流拦截措施。如果上缘未达到分水岭且顶部无拦截措施就是有来水。开挖面的土质因子反映下垫面质地和土体密度对外营力剥蚀和搬动的影响坡度因子和坡长因子的系数查表 9 到表 11上方无来水和表 12 到表 14上方有来水。这些系数表是按土石质类型和坡度坡长范围给的取值时先确定土石质类型再按实际坡度和坡长查对应系数。3.4 工程堆积体土石质因子与堆放方式工程堆积体的标准下垫面是宽 1m、水平坡长 5m、坡度 25° 的全土质裸露面。堆积体的土石质因子反映物质被外营力剥蚀和搬运的敏感程度和开挖面的土质因子不是一回事。上方无来水的堆积体公式里用表 9 到表 11 的系数上方有来水的用表 12 到表 14。另外风力作用下堆积体还有堆放方式因子查表 17。实际项目里堆积体往往是多级堆放的每级坡度和坡长不同要按计算单元分别算再累加。如果堆积体顶部有较大平台且受上方来水侵蚀就按上方有来水处理。4. 风力作用与新增流失量容易被忽略的测算分支4.1 风力作用下一般扰动地表起动风速与粗糙干扰因子风力作用下的测算分一般扰动地表和工程堆积体两类。一般扰动地表的公式涉及风蚀可蚀性因子、起动风速、粗糙干扰因子等。起动风速是地表颗粒物开始运动时的临界风速粗糙干扰因子表征植被、砾石等地表覆盖对风蚀的削弱作用。表 15 给了典型下垫面的 Z0 取值空气动力学粗糙度表 16 给了不同土壤类型的风蚀可蚀性因子。实际计算时先确定下垫面类型查 Z0再确定土壤类型查风蚀可蚀性因子然后结合风速资料算风蚀量。碾压地表的风蚀量可以参照一般扰动地表计算但要注意碾压后地表粗糙度变了Z0 取值要相应调整。4.2 风力作用下工程堆积体堆放方式因子的影响堆积体的风蚀测算除了土石质因子还要考虑堆放方式因子表 17。不同堆放方式如锥形堆、条形堆、平台堆对风的扰动响应不同因子取值也不一样。我一般会按堆积体的实际形态选堆放方式如果是排土场那种大面积平台按平台堆取如果是临时转运的锥形堆按锥形堆取。如果堆积体表面有覆盖或压实措施粗糙干扰因子要相应调低。4.3 新增土壤流失量估算在建项目与预测项目的区别标准第 10 章把应用规定分成三块在建项目测算规定、预测项目测算规定、新增土壤流失量估算规定。在建项目的土壤流失量测算是事中监测要用实际观测资料或现场查勘数据不能直接用预测阶段的经验值。预测项目是事前预测参数取值可以按附录 C 的参考值结合现场调查确定。新增土壤流失量是项目扰动后的流失量减去扰动前原地貌的流失量扰动前流失量按原地貌类型查相关标准或实测。这里有个容易混的地方新增流失量不是“扰动后流失量”而是“扰动后减去扰动前”。如果项目区原地貌本身流失量就大新增量可能不大但总流失量仍然很高报告里要分开说清楚。5. 避坑与排查参数取值和流程执行中的五个常见翻车点5.1 降雨侵蚀力因子 R 用错时间尺度现象同一个项目不同人算出来的 R 值差好几倍。原因有人用年降雨公式算日尺度有人用日降雨公式算年尺度还有人把多年平均公式套到单年数据上。解决先明确测算时段再选对应公式。日尺度必须用日降雨公式逐日算再累加月尺度用月公式年尺度用年公式多年平均用多年平均公式。不同尺度不能混用混用就是玄学。5.2 坡度因子计算时角度没转弧度现象Sy 算出来是负数或者大得离谱。原因Python 的 math.sin 接受弧度Excel 的 SIN 也接受弧度但很多人直接填角度。解决计算前统一用 math.radians() 或 Excel 的 RADIANS() 转一下。标准里专门提醒了这一点但每年还是有人翻车。5.3 计算单元划分时坡度差 5° 的边界处理现象两个相邻地块坡度差刚好 5°不知道要不要分开。原因标准写的是“相差大于 5° 的划分为不同计算单元”大于 5° 才分等于 5° 不分。解决严格按“大于”执行差 5.0° 不分差 5.1° 就分。实际项目中如果坡度变化连续建议按 5° 间隔分段避免边界争议。5.4 上方来水判断只看地形不看措施现象开挖面上缘没到分水岭但顶部修了截水沟还是按有来水算结果偏大。原因标准里“上方无来水”的条件是“上缘已达到或越过分水岭或在顶部有截排水沟等坡面径流拦截措施”两个条件满足一个就算无来水。解决先看地形再看措施。有截排水沟且能有效拦截的按无来水算没有措施且上缘未到分水岭的按有来水算。5.5 扰动单元数量 20 时的统计逻辑搞反现象项目总流失量算出来比所有典型单元之和还大。原因把“同一类型下同等规模典型单元平均值 × 该类型规模单元数量”理解成了“所有单元平均值 × 总数量”。解决先按类型和规模分组每组算典型单元平均值再乘以该组单元数量最后各组求和。不是所有单元一起平均。6. 从 PDF 到计算表一套可复用的参数取值与复核习惯把 SL 773-2018 用熟之后我养成了一个习惯每接一个新项目先建一张“计算单元要素表”列包括单元编号、扰动类型、下垫面类型、面积、坡度、坡长、土壤质地、植被覆盖度、水保措施、上方汇水面积、降雨资料年限、风速资料年限。这张表填完后面选公式、查参数、算流失量就是流水线作业。参数取值上我一般按这个优先级项目区实测 邻近气象站 县域统计 附录 C 参考值。附录 C 的县级行政区降雨侵蚀力因子和土壤可蚀性因子是多年平均逐月和年值用的时候要注意如果项目测算时段不是多年平均不能直接拿附录 C 的值当某一年或某一月的值用得按实际降雨资料算。复核环节我会重点检查三件事一是 R 值的时间尺度是否和测算时段一致二是坡度因子计算时角度有没有转弧度三是上方来水判断有没有漏掉截排水措施。这三项查完基本能排除 80% 的低级错误。还有一个容易忽略的点标准里说“同一计算单元水力侵蚀和风力侵蚀产生的土壤流失量应分别计算”但没说不可以加总。实际报告里如果项目区既有水力侵蚀又有风力侵蚀要分别列出两项再给合计值不能只给一个总数。最后说一个我自己的教训早期做方案时觉得附录 C 的参考值方便直接拿来用结果评审专家问“你这个 R 值是多年平均还是测算年的”答不上来。从那以后我每次用附录 C 的值都会在报告里注明“因缺乏实测降雨资料采用附录 C 县级行政区多年平均降雨侵蚀力因子”并说明测算时段与多年平均的差异。这个习惯帮我省了不少解释成本。希望帮到你。本文还有配套的精品资源点击获取