1.一种无资料地区径流估算方法,其特征在于,所述方法包括:获取驱动数据,所述驱动数据包括预设时间段内流域月尺度降水量、潜在蒸散发量和土壤体积含水率;根据由水量输入以及水量输出通量计算的土壤水储量变化与由研究流域月尺度水量状态量的变化之间的均方根误差来构建流域月尺度水量平衡变化的目标函数,将所述目标函数中的未知参数作为待求解参数;对所述待求解参数进行初始化,得到待求解参数的初始值;将所述待求解参数的初始值和所述驱动数据输入所述目标函数,得到目标函数的初始函数值;基于所述待求解参数的初始值和所述初始函数值,计算所述目标函数对所述待求解参数的梯度函数;计算所述梯度函数在上一轮待求解参数的优化值下的梯度值;根据所述上一轮待求解参数的优化值、所述梯度值以及预设的优化精度,将所述上一轮待求解参数的优化值沿着所述梯度值的负方向进行迭代计算,得到所述待求解参数的优化值;将所述待求解参数的优化值输入所述目标函数得到目标函数值;在所述目标函数值满足预设的迭代终止条件的情况下,根据所述待求解参数的优化值估算预设时间段内流域月尺度径流值。
2.根据权利要求1所述的方法,其特征在于,所述待求解参数包括多个待求解参数,共同组成参数空间;所述计算所述目标函数对所述待求解参数的梯度函数,包括:计算所述目标函数对所述参数空间中每一个待求解参数的梯度函数;所述计算所述梯度函数在上一轮待求解参数的优化值下的梯度值,包括:计算所述每一个待求解参数的梯度函数在对应的上一轮待求解参数的优化值下的梯度值;所述根据所述上一轮待求解参数的优化值、所述梯度值以及优化精度,将所述上一轮待求解参数的优化值沿着所述梯度值的负方向进行迭代计算,得到所述待求解参数的优化值,包括:根据所述上一轮每一个待求解参数的优化值、对应的所述梯度值以及优化精度,将所述上一轮每一个待求解参数的优化值沿着所述梯度值的负方向进行迭代计算,得到对应的优化后的参数值,所述每一个优化后的参数值共同组成优化的参数空间;所述将所述待求解参数的优化值输入所述目标函数得到目标函数值,包括:将所述优化的参数空间输入所述目标函数得到目标函数值。
3.根据权利要求1所述的方法,其特征在于,所述构建流域月尺度水量平衡变化的目标函数,包括:根据水量平衡模型中水量的输入通量与输出通量的差值除以研究流域土层的厚度得到第一土壤水储量变化参数;根据研究流域水量的状态量计算得到的土壤水储量变化作为第二土壤水储量变化参数;根据所述第一土壤水储量变化参数、所述第二土壤水储量变化参数的均方根误差构建所述流域月尺度水量平衡变化的目标函数。
4.根据权利要求3所述的方法,其特征在于,所述水量平衡模型的构建过程包括:根据预设时间段内水量输入通量函数、水量输出通量函数以及水量的状态量之间的关系构建所述水量平衡模型。
5.根据权利要求4所述的方法,其特征在于,所述待求解参数包括第一待求解参数、第二待求解参数、第三待求解参数、第四待求解参数、第五待求解参数和第六待求解参数;所述水量输出通量函数的构建过程,包括:根据实际蒸散发量函数、土壤深层渗漏函数以及流域径流量函数联合得到所述水量输出通量函数;所述实际蒸散发量函数包括流域月尺度潜在蒸散发、土壤体积含水率、所述第一待求解参数以及第二待求解参数之间的函数关系;其中,第一待求解参数表示估算的所述潜在蒸散发的系统偏差,第二待求解参数用于确定所述实际蒸散发量函数的曲线形状;所述土壤深层渗漏函数包括土壤体积含水率、所述第三待求解参数以及第四待求解参数之间的函数关系;其中,第三待求解参数表示深层渗漏量的最大值,第四待求解参数表示土壤含水层水量流失的速率;所述流域径流量函数包括所述水量输入通量函数、土壤体积含水率以及所述第五待求解参数之间的函数关系;其中,第五待求解参数表示所述水量输入通量函数中降水产流的速率;所述方法还包括:根据所述水量平衡模型确定第六待求解参数;其中,第六待求解参数表示所述流域研究土层的厚度。
6.一种无资料地区径流估算装置,其特征在于,所述装置包括:获取模块,用于获取驱动数据,所述驱动数据包括预设时间段内流域月尺度降水量、潜在蒸散发量和土壤体积含水率;构建模块,用于根据由水量输入以及水量输出通量计算的土壤水储量变化与由研究流域月尺度水量状态量的变化之间的均方根误差来构建目标函数,将目标函数中的未知参数作为待求解参数;初始化模块,用于对所述待求解参数进行初始化,得到待求解参数的初始值;以及将所述待求解参数的初始值和所述驱动数据输入所述目标函数,得到目标函数的初始函数值;优化模块,用于基于所述待求解参数的初始值和所述初始函数值,计算所述目标函数对所述待求解参数的梯度函数;计算所述梯度函数在上一轮待求解参数的优化值下的梯度值;根据所述上一轮待求解参数的优化值、所述梯度值以及预设的优化精度,将所述上一轮待求解参数的优化值沿着所述梯度值的负方向进行迭代计算,得到所述待求解参数的优化值,将所述待求解参数的优化值输入所述目标函数得到目标函数值;径流计算模块,用于在所述目标函数值满足预设的迭代终止条件的情况下,根据所述待求解参数的优化值估算预设时间段内流域月尺度径流值。
7.根据权利要求6所述的装置,其特征在于,所述待求解参数包括多个待求解参数,共同组成参数空间;优化模块还用于计算所述目标函数对所述参数空间中每一个待求解参数 的梯度函数;以及计算所述每一个待求解参数的梯度函数在对应的上一轮待求解参数的优化值下的梯度值;以及根据所述上一轮每一个待求解参数的优化值、对应的该梯度值以及优化精度,将所述上一轮每一个待求解参数的优化值沿着该梯度值的负方向进行迭代计算,得到对应的优化后的参数值,所述每一个优化后的参数值共同组成优化的参数空间。
8.根据权利要求6所述的装置,其特征在于,所述构建模块还用于根据水量平衡模型中水量的输入通量与输出通量的差值除以研究流域土层的厚度得到第一土壤水储量变化参数;以及根据研究流域水量的状态量计算得到的土壤水储量变化作为第二土壤水储量变化参数;以及根据所述第一土壤水储量变化参数、所述第二土壤水储量变化参数的均方根误差构建所述流域月尺度水量平衡变化的目标函数。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。