[发明专利]适用于高海拔区计算多年平均降雨-融雪径流侵蚀力的方法

权利要求书

1.一种适用于高海拔区计算多年平均降雨-融雪径流侵蚀力的方法，其特征在于：其包括以下步骤：

S1、确定高海拔区需要计算多年平均降雨-融雪径流侵蚀力的计算时间序列，得到侵蚀力计算的起始年与结束年，并在确定计算时间序列后选择出日降雨量大于0.1mm的日降雨，所需降雨基础数据的时间序列为计算时间序列的前一年10月1日至结束年12月31日；

S2、将选择出的日降雨按照降雨发生的月份，以前一年10月份至当年9月份为一个周期作为一个单元年，并将每一个单元年区分为暖季降雨月份和冷季降雨月份，暖季降雨月份为6-9月份，冷季降雨月份为前一年10月份至当年的5月份；

S3、筛选时间序列的所有周期内暖季降雨月份中日降雨量大于12mm的降雨日；

S4、分别按照暖季降雨侵蚀力和冷季融雪径流侵蚀力计算公式，得到时间序列内每一个单元年对应时段的降雨侵蚀力，两者相加得到该一个单元年的年降雨-融雪径流侵蚀力；

R＝R_暖+R_冷

式中：R为一个单元年的年降雨-融雪径流侵蚀力，单位为MJ.mm.hm^-2.h^-1；R_暖为暖季降雨侵蚀力，单位为MJ.mm.hm^-2.h^-1；R_冷为冷季融雪径流侵蚀力，单位为MJ.mm.hm^-2.h^-1；

其中，R_暖的计算公式如下：

式中：P_日为6-9月大于12mm的第j个降雨日的降雨量，单位为mm，k为日降雨大于12mm的第一个降雨日，j为日降雨大于12mm的第j个降雨日；

R_冷的计算公式如下：

R_冷＝0.125P_10-5^1.437

式中：P_10-5为上年10月1日至当年5月31日的总降雨量，单位为mm；

S5、对于不同年份降雨-融雪径流侵蚀力的计算，重复步骤S4，依次得到时间序列内多年的暖季降雨侵蚀力和冷季融雪径流侵蚀力，并依次计算得到时间序列内每个单元年的年降雨-融雪径流侵蚀力；

S6、对步骤S5得到的每个单元年的年降雨-融雪径流侵蚀力求平均值，得到高海拔区时间序列内多个单元年的平均降雨-融雪径流侵蚀力；

式中：为计算时间序列内多个单元年的平均降雨-融雪径流侵蚀力，N为计算时间序列的单元年数，1,2,3,...。

2.根据权利要求1所述的适用于高海拔区计算多年平均降雨-融雪径流侵蚀力的方法，其特征在于：步骤S3中，暖季降雨侵蚀力计算公式的确定方法如下所述：

根据步骤S1中获得的实际测量的暖季期间精细的降雨过程数据计算得到暖季降雨的降雨量和EI₃₀，在excel表格中对暖季降雨的降雨量和EI₃₀散点描出图像曲线，该图像曲线为一条直线，在Matlab软件中采用regress函数作最小二乘法作线性拟合得到暖季侵蚀力关系表达式。

3.根据权利要求1所述的适用于高海拔区计算多年平均降雨-融雪径流侵蚀力的方法，其特征在于：步骤S3中，冷季融雪径流侵蚀力计算公式的确定方法如下所述：

将冷季期间实际测量的降雨量和冷季期间融雪径流侵蚀力散点，在excel表格中对冷季期间降水和融雪径流侵蚀力散点描出图像曲线，该图像曲线为一条直线，在Matlab软件中采用regress函数作最小二乘法作线性拟合得到冷季融雪径流侵蚀力关系表达式。

4.根据权利要求1所述的适用于高海拔区计算多年平均降雨-融雪径流侵蚀力的方法，其特征在于：冷季期间融雪径流侵蚀力的计算包括以下子步骤：

①将河道观测断面的泥沙输移量与泥沙输移比相除得到土壤侵蚀量；

②根据区域坡度坡长因子、植被覆盖因子的空间分布状况，将土壤侵蚀量分别除以坡度坡长因子、植被覆盖因子，得到区域的降雨-融雪径流侵蚀力；

③降雨-融雪径流侵蚀力减去暖季期间的降雨侵蚀力，得到的数值即为冷季期间融雪径流侵蚀力。