一种土壤对流域地表径流调节能力的评价方法技术

本发明专利技术公开了一种土壤对流域地表径流调节能力的评价方法，该方法通过对流域土壤进行离散化处理，将流域土壤分成若干个土壤单元，并计算每个土壤单元在场次降水的不同时段的地表产流量，再根据在场次降水的不同时段的地表产流量计算出单个土壤单元中土壤对地表径流调节能力的评价指标，最后得出整个流域尺度土壤对地表径流调节能力进行定量评价得到流域尺度土壤对地表径流调节能力的评价指标，还方法从区域水循环机理层面实现流域尺度上土壤对地表径流调节作用进行定量评价，其评价客观、准确，满足当前生态海绵流域建设、流域旱涝灾害有效应对的现实需求。

【技术实现步骤摘要】
一种土壤对流域地表径流调节能力的评价方法
本专利技术涉及一种土壤对流域地表径流调节能力的评价方法。
技术介绍
当前，我国面临严峻的水资源短缺、流域旱涝灾害频发等多种水问题；因此，流域地表径流过程的“削峰补枯”日益受到国内相关管理机构及科研学者的重视。而流域土壤是由大小不同、形状不规则的土壤颗粒组成，是一个多孔体，具有吸水、蓄水等的“海绵”功能，对流域地表径流过程具有显著的调节能力；这在学术界已经得到共识。但是，由于当前土壤对流域地表径流调节能力尚未形成统一的认识，土壤对径流调节能力没有的可靠的评判标准，导致不同流域甚至同一流域的评价结果难以对比，不能为生态海绵流域及水生态文明建设提供充分的参考数据。
技术实现思路
本专利技术的目的是提供一种土壤对流域地表径流调节能力的评价方法，可以对不同流域的土壤对径流调节能力进行客观的、准确的定量评价，为生态海绵流域及水生态文明建设提供充分的参考数据。为解决上述技术问题，本专利技术提供一种土壤对流域地表径流调节能力的评价方法，包括以下步骤:S1：根据流域土壤质地、土层厚度、土壤孔隙度、土体构型和有机质含量对流域尺度土壤进行空间离散化处理，将流域土壤分成若干个土壤单元；S2：根据土壤的土壤特性参数计算出各土壤单元不同土层的土壤水动力学特性参数；S3：根据各土壤单元不同土层的土壤水动力学特性参数，采用Green-Ampt入渗模型计算土壤入渗率；再根据土壤单元初始含水量计算各土壤单元在场次降水不同时段的地表产流量；S4：根据场次降水过程中累计降雨量与累计地表产流量构建单个土壤单元中土壤对地表径流调节能力的评价指标；S5：根据单个土壤单元中土壤对地表径流调节能力的评价指标，采用加权平均法对流域尺度土壤对地表径流调节能力进行定量评价得到流域尺度土壤对地表径流调节能力的评价指标。进一步地，步骤S2中土壤特性参数为土壤粒径级配数据、土壤有机质含量和土壤盐度；步骤S2中土壤水动力学参数为土壤容重、土壤凋萎系数、田间持水量、饱和导水率和水利传导度。进一步地，步骤S3的根据每个土壤单元的初始土壤含水量、降雨量、实际蒸发量、降雨强度、土壤入渗率和地表产流量计算出每个土壤单元格在场次降水的不同时段的地表产流量的具体计算方法为：当SW0>SWs时：若P-E<0，RI<In，则不产流：R＝0；若P-E<0，RI>In，则超渗产流：R＝(RI-In)·det(t)；若P-E>0，RI<In，则蓄满产流：R＝P-E；若P-E>0，RI>In，则超渗产流+蓄满产流：RT＝P-E；当SW0<SWs时：若P-E<(SW0<SWs)，RI<In，则不产流：R＝0；若P-E<(SW0<SWs)，RI>In，则超渗产流：R＝(RI-In)·det(t)；若P-E>(SW0<SWs)，RI<In，则蓄满产流：R＝P-E-(SW0-SWs)；若P-E>(SW0<SWs)，RI>In，则超渗产流+蓄满产流：RT＝P-E-(SW0-SWs)；其中，SW0为初始土壤含水量，SWs为最终土壤含水量，P为每个土壤单元的降雨量(mm)；E为每个土壤单元实际蒸发量(mm)；RI为降雨强度(mm/h)；In为土壤入渗率(mm/hr)；R为地表产流量(mm)。进一步地，步骤S4的根据场次降水过程中累计降雨量与累计地表产流量构建单个土壤单元中土壤对地表径流调节能力的评价指标的具体计算公式为：式中，η为单个土壤单元中土壤对地表径流调节能力的评价指标，Ri为场次降水过程中的累计降雨量，Pi为场次降水过程中的累计径流量，baseflow为计流量，进一步地，步骤5的根据单个土壤单元中土壤对地表径流调节能力的评价指标，采用加权平均法对流域尺度土壤对地表径流调节能力进行定量评价采用的计算公式为：式中，ηT为流域尺度土壤对地表径流调节能力的评价指标，ηi为单个土壤单元中土壤对地表径流调节能力的评价指标。本专利技术的有益效果为：本申请通过从区域水循环机理层面实现流域尺度上土壤对地表径流调节作用的客观的、准确的定量评价，满足当前生态海绵流域建设、流域旱涝灾害有效应对的现实需求。具体实施方式下面对本专利技术的具体实施方式进行描述，以便于本
的技术人员理解本专利技术，但应该清楚，本专利技术不限于具体实施方式的范围，对本
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本专利技术的精神和范围内，这些变化是显而易见的，一切利用本专利技术构思的专利技术创造均在保护之列。一种土壤对流域地表径流调节能力的评价方法，包括以下步骤：S1：根据流域土壤质地、土层厚度、土体构型、土壤孔隙度和有机质含量对流域尺度土壤进行离散化处理，将流域土壤分成若干个土壤单元。S2：根据土壤的土壤粒径级配数据、土壤有机质含量和土壤盐度等土壤特性参数计算出土壤单元不同土层的土壤容重、土壤凋萎系数、田间持水量、饱和导水率和水利传导度等土壤水动力学特性参数。S3：根据各土壤单元不同土层的土壤水动力学特性参数，采用Green-Ampt入渗模型计算土壤入渗率；再根据土壤单元初始含水量计算各土壤单元在场次降水不同时段的地表产流量；根据每个土壤单元的初始土壤含水量、降雨量、实际蒸发量、降雨强度、土壤入渗率和地表产流量计算出每个土壤单元格在场次降水的不同时段的地表产流量，其具体计算方法如下：当SW0>SWs时，即前期土壤水己蓄满时：若P-E<0，RI<In，则不产流：R＝0；若P-E<0，RI>In，则超渗产流：R＝(RI-In)·det(t)；若P-E>0，RI<In，则蓄满产流：R＝P-E；若P-E>0，RI>In，则超渗产流+蓄满产流：RT＝P-E；当SW0<SWs时，即前期土壤水未蓄满时：若P-E<(SW0<SWs)，RI<In，则不产流：R＝0；若P-E<(SW0<SWs)，RI>In，则超渗产流：R＝(RI-In)·det(t)；若P-E>(SW0<SWs)，RI<In，则蓄满产流：R＝P-E-(SW0-SWs)；若P-E>(SW0<SWs)，RI>In，则超渗产流+蓄满产流：RT＝P-E-(SW0-SWs)；其中，SW0为初始土壤含水量，SWs为最终土壤含水量，P为每个土壤单元的降雨量(mm)；E为每个土壤单元实际蒸发量(mm)；RI为降雨强度(mm/h)；In为土壤入渗率(mm/hr)；R为地表产流量(mm)。S4：根据场次降水过程中累计降雨量与累计径流量构建单个土壤单元中土壤对地表径流调节能力的评价指标，具体采用的计算公式为：式中，η为单个土壤单元中土壤对地表径流调节能力的评价指标，Ri为场次降水过程中的累计降雨量，Pi为场次降水过程中的累计径流量，baseflow为计流量，S5：根据单个土壤单元中土壤对地表产流量调节能力的评价指标，采用加权平均法对流域尺度土壤对地表径流调节能力进行定量评价得到流域尺度土壤对地表径流调节能力的评价指标，具体采用的计算公式为：式中，ηT为流域尺度土壤对地表径流调节能力的评价指标，ηi为单个土壤单元中土壤本文档来自技高网...

【技术保护点】
一种土壤对流域地表径流调节能力的评价方法，其特征在于，包括以下步骤S1：根据流域土壤质地、土层厚度、土壤孔隙度、土体构型和有机质含量对流域尺度土壤进行空间离散化处理，将流域土壤分成若干个土壤单元；S2：根据土壤的土壤特性参数计算出各所述土壤单元不同土层的土壤水动力学特性参数；S3：根据各土壤单元不同土层的所述土壤水动力学特性参数，采用Green‑Ampt入渗模型计算土壤入渗率，再根据土壤单元初始含水量计算其在场次降水不同时段的地表产流量；S4：根据场次降水过程中累计降雨量与累计地表产流量构建单个土壤单元中土壤对地表径流调节能力的评价指标；S5：根据所述单个土壤单元中土壤对地表径流调节能力的评价指标，采用加权平均法对流域尺度土壤对地表径流调节能力进行定量评价得到流域尺度土壤对地表径流调节能力的评价指标。

【技术特征摘要】
1.一种土壤对流域地表径流调节能力的评价方法，其特征在于，包括以下步骤S1：根据流域土壤质地、土层厚度、土壤孔隙度、土体构型和有机质含量对流域尺度土壤进行空间离散化处理，将流域土壤分成若干个土壤单元；S2：根据土壤的土壤特性参数计算出各所述土壤单元不同土层的土壤水动力学特性参数；S3：根据各土壤单元不同土层的所述土壤水动力学特性参数，采用Green-Ampt入渗模型计算土壤入渗率，再根据土壤单元初始含水量计算其在场次降水不同时段的地表产流量；S4：根据场次降水过程中累计降雨量与累计地表产流量构建单个土壤单元中土壤对地表径流调节能力的评价指标；S5：根据所述单个土壤单元中土壤对地表径流调节能力的评价指标，采用加权平均法对流域尺度土壤对地表径流调节能力进行定量评价得到流域尺度土壤对地表径流调节能力的评价指标。2.根据权利要求1所述的土壤对流域地表径流调节能力的评价方法，其特征在于，步骤S2中所述土壤特性参数为土壤粒径级配数据、土壤有机质含量和土壤盐度；步骤S2中所述土壤水动力学参数为土壤容重、土壤凋萎系数、田间持水量、饱和导水率和水利传导度。3.根据权利要求1所述的土壤对流域地表径流调节能力的评价方法，其特征在于，步骤S3所述的根据每个土壤单元的初始土壤含水量、降雨量、实际蒸发量、降雨强度、土壤入渗率和地表产流量计算出每个土壤单元格在场次降水的不同时段的地表产流量的具体计算方法为：当SW0>SWs时：若P-E<0，RI<In，则不产流：R＝0；若P-E<0，RI>In，则超渗产流：R＝(RI-In)·det(t)；若P-E>0，RI<In，则蓄满产流：R＝P-E；若P-E>0...

【专利技术属性】
技术研发人员：李淑惠，邢子强，代永江，赵珊，刘姗姗，于志磊，
申请(专利权)人：沈阳市水利建筑勘测设计院，中国水利水电科学研究院，
类型：发明
国别省市：辽宁,21