[发明专利]一种冰川融水径流量的估算方法

权利要求书

1.一种冰川融水径流量的估算方法，其特征在于，包括：

S1、获取冰川地理分布信息并计算得到驱动数据，

所述冰川地理分布信息包括冰川面积、冰川的矢量边界数据，

所述驱动数据的计算方法为：对气象数据进行插值处理、降水雨雪分离处理以及正积温计算后，得到驱动数据，

所述气象数据包括：冰川表面的气温、冰川表面的降水量，所述降水量包括降雪量和降雨量，

所述驱动数据包括冰川表面月降雨量、冰川表面月降雪量、冰川表面正积温、数字高程模型；

S2、模拟冰川消融：将所述冰川表面正积温与潜在融雪正积温进行对比判断，得出冰川的每月融雪量和每月融冰量，再通过度日因子模型的消融原理公式计算冰川表面物质平衡量，即冰川体积变化量，最后通过将冰川体积变化量代入面积-体积转换公式，计算得出冰川面积变化量；

S3、计算冰川融水径流量：通过所述冰川面积变化量结合冰川表面物质平衡量计算冰川融水径流量。

2.如权利要求1所述的方法，其特征在于，所述对气象数据进行插值处理包括以下内容：

S1011、利用气象站经度、纬度、海拔及气象站观测到的气象数据，建立月度气象数据与气象站经度、纬度、海拔之间的插值公式；

S1012、利用数字高程模型将每条冰川切分为等间隔的冰川海拔带，获取冰川海拔带的经度、纬度、海拔，逐月向步骤S1011所得插值公式中代入每条冰川海拔带中心的经度、纬度、海拔，以获取每条冰川海拔带月度的气象数据，每条冰川海拔带月度的气象数据包括每条冰川海拔带月度的降水量、每条冰川海拔带月平均气温。

3.如权利要求2所述的方法，其特征在于，所述降水雨雪分离处理包括以下内容：

S1021、在对气象数据进行插值处理后，得到每条冰川海拔带月度的降水量，通过双临界气温法对每条冰川海拔带月度的降水量进行雨雪分离处理，所述雨雪分离处理的计算公式为：

上式中，S_j,m为第j个冰川海拔带在第m个月的降雪量，单位为mm，P_j,m为第j个冰川海拔带在第m个月的月降水量，单位为mm，T_j,m为月均气温，单位为℃，T_s为降雪临界气温，单位为℃，T_r为降雨临界气温，单位为℃，j为分割的海拔带数量，m为一个消融周期内的月度数量，

S1022、误差校正：

考虑到气象站观测的降水存在误差，对降雨和降雪进行误差校正，所述误差校正的计算公式为：

R_j,m＝B_r×(P_j,m―S_j,m)

S_j,m＝B_s×S_j,m

上式中，S_j,m为第j个冰川海拔带的在第m个月的降雪量，单位为mm，P_j,m为第j个冰川海拔带在第m个月的月降水量，单位为mm，R_j,m为第j个冰川海拔带在第m个月的月降雨量，单位为mm，B_r为降雨的校正系数，B_s为降雪的校正系数，j为分割的海拔带数量，m为一个消融周期内的月度数量。

4.如权利要求2所述的方法，其特征在于，所述正积温计算包括以下内容：

在对气象数据进行插值处理后，得到每条冰川海拔带月平均气温，通过对每条冰川海拔带月平均气温积分获取冰川表面正积温，冰川表面正积温的计算公式为：

上式中，PDD_j,m为第j个冰川海拔带在第m个月的冰川表面正积温，T为计算时段内的平均温度，单位为℃，T_j,m为第m个月内各时间点气温，单位为℃，σ_m为气温分布的标准差，j为分割的海拔带数量，m为一个消融周期内的月度数量。