[发明专利]一种多弹性对径流变化驱动因素的定量分割方法

权利要求书

1.一种多弹性对径流变化驱动因素的定量分割方法，其特征在于，包括如下步骤：

步骤1，收集流域内至少40年的年径流量资料以及日气象要素资料，其中，气象要素资料包括降水量、最高气温、最低气温、平均相对湿度、平均风速及日照时数；

步骤2，结合最高气温、最低气温、平均相对湿度、平均风速及日照时数，运用Penman-Montieth公式计算流域日潜在蒸发量，进而得到流域多年平均潜在蒸发量；

步骤3，根据流域内的水量平衡方程，结合流域多年平均径流量与多年平均降水量，计算流域多年平均实际蒸发量；

步骤4，采用单参数的Budyko方程，结合多年平均降水量、多年平均潜在蒸发量及多年平均实际蒸发量，计算多年平均的流域特性参数；

步骤5，计算径流对降水量、最高气温、最低气温、平均相对湿度、平均风速、日照时数以及流域特性这七个驱动因素的弹性系数；

步骤6，分析流域内年径流量的变化趋势，判断年径流过程的突变点，将突变点之前作为基准期，突变点之后作为变化期；

步骤7，根据弹性系数、多年平均径流量、各个驱动因素的多年平均值、变化期径流量均值与基准期径流量均值的差值及各个驱动因素变化期均值与基准期均值的差值，计算各个驱动因素对径流变化的贡献率。

2.根据权利要求1所述多弹性对径流变化驱动因素的定量分割方法，其特征在于，步骤2所述日潜在蒸发量的计算公式为：

其中，E_{0_d}为日潜在蒸发量，Δ为饱和水汽压斜率，R_n为净辐射，G为土壤热通量，γ为干湿常数，T为平均温度，U为平均风速，VPD为饱和水汽压差。

3.根据权利要求1所述多弹性对径流变化驱动因素的定量分割方法，其特征在于，步骤3所述多年平均实际蒸发量的计算公式为：

E＝P-Q

其中，E为多年平均实际蒸发量，P为多年平均降水量，Q为多年平均径流量。

4.根据权利要求1所述多弹性对径流变化驱动因素的定量分割方法，其特征在于，步骤4所述多年平均的流域特性参数的计算公式为：

E＝P/[1+(P/E₀)ⁿ]^1/n

其中，E为多年平均实际蒸发量，P为多年平均降水量，E₀为多年平均潜在蒸发量，n为多年平均的流域特性参数。

5.根据权利要求1所述多弹性对径流变化驱动因素的定量分割方法，其特征在于，步骤6所述判断年径流过程的突变点方法为：Pettitt检验方法。

6.根据权利要求1所述多弹性对径流变化驱动因素的定量分割方法，其特征在于，步骤7所述各个驱动因素对径流变化的贡献率计算公式为：

其中，P，TA，TI，U，N，RH，n分别为多年平均降水量、多年平均最高气温、多年平均最低气温、多年平均风速、多年平均日照时数、多年平均相对湿度、多年平均的流域特性参数；η_X为贡献率，ε_X为弹性系数，Q为多年平均径流量，ΔQ为变化期径流量均值与基准期径流量均值的差值，ΔX为各个驱动因素变化期均值与基准期均值的差值。