[发明专利]一种考虑冰川动态融化的干旱区冰川融水流量计算方法

权利要求书

1.一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，包括以下步骤：

获取流域内的所有冰川的冰川面积；

根据冰川体积-面积比例关系，计算冰川体积；

将冰川面积和冰川体积分配到流域的各个子流域内，计算出各个子流域的总冰川覆盖面积和冰川体积；

根据纬度和日照时数求得流域内的太阳辐射量；

获取流域内和流域周边的气象数据，根据泰森多边形方法，将气象数据插值到子流域；

计算融冰当日的冰川表面温度；

计算每个子流域上的基于冰面温度和太阳辐射量的冰川融化量，从而获得干旱区冰川融水流量。

2.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，计算冰川体积的公式如下：

其中V是单条冰川的冰川体积；A是单条冰川的面积。

3.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，太阳辐射量的计算公式如下：

其中，R_s为太阳辐射量，R_a为地球外辐射，N为最大可能日照时数，a_s表示阴天时到达地球表面的地球外辐射的透过系数；b_s表示晴天时到达地球表面的地球外辐射的透过系数，n为常数。

4.根据权利要求3所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，地球外辐射的计算公式如下：

G_sc为太阳常数，d_r为反转日地平均距离，ω_s为日落时角，为纬度，δ为太阳磁偏角。

5.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，融冰当日的冰川表面温度的计算公式如下：

第d_n天的冰川温度为表示为前一天冰川温度和当天气温的函数。

6.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，还包括构建温度辐射因子：

Gmfmx和Gmfmn分别是8月1日和2月1日的冰川融化因子，也就是一年中最大和最小融冰速率出现的时间；dn是一年中的儒略日。

7.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，计算每个子流域上的基于冰面温度和太阳辐射量的冰川融化量，计算公式如下；

其中，GLA_mlt是当天子流域上的冰川融水量；gla_cov是子流域内的冰川覆盖面积比例；T_max是当日最高气温，Gmtmp是冰川融化基础温度，R_s为当天的子流域太阳辐射。

8.根据权利要求1所述的一种考虑冰川动态融化的干旱区冰川融水流量计算方法，其特征在于，还包括根据冰川体积-面积比例关系，更新冰川体积δV和冰川面积δA的相应相对变化，具体如下：

其中，GLA_mlt是当天子流域上的冰川融水量。