[发明专利]无资料地区径流估算方法、装置、计算机设备和存储介质

说明书

本申请涉及一种无资料地区径流估算方法、装置、计算机设备和存储介质。所述方法包括：通过获取驱动数据，所述驱动数据包括预设时间段内流域月尺度降水量、潜在蒸散发量和土壤体积含水率；然后，构建流域月尺度水量平衡变化的目标函数，并确定目标函数的待求解参数；对待求解参数进行初始化，将待求解参数的初始值与驱动数据输入目标函数，得到目标函数的初始函数值；并基于所述待求解参数的初始值和所述初始函数值，对所述待求解参数进行迭代计算；最后，根据所述待求解参数的优化值估算预设时间段内流域月尺度径流值。采用本方法能够在不依赖任何实测径流资料的条件下，通过输入较少的驱动数据即可对无资料地区的径流值进行估算，估算准确率高。

技术领域

本申请涉及水文水资源管理技术领域，特别是涉及一种无资料地区径流估算方法、装置、计算机设备和存储介质。

背景技术

河流径流量是水文循环中的重要变量，也是用于区域水资源评估、生态环境保护和气候变化分析的基础数据。相关的径流监测数据主要来源于水文监测站点，虽然其精度较高，但是全球站点分布十分不均匀且多数站点观测时段较短，受地形及观测成本的制约，许多无资料区缺少可靠的径流数据。

目前的研究通常采用搭建水文模型的方法模拟流域径流量，然而相关方法严重依赖实测径流资料对模型参数进行率定，在无资料地区应用受限。

目前使用相邻流域的实测流量数据率定该流域的模型参数并移植至无资料地区，但是由于不同流域的独特性，模拟准确率不高。

发明内容

基于此，有必要针对上述技术问题，提供一种能够提高无资料地区径流估算准确率的无资料地区径流估算方法、装置、计算机设备和存储介质。

第一方面，提供了一种无资料地区径流估算方法，该方法包括：

获取驱动数据，该驱动数据包括预设时间段内流域月尺度降水量、潜在蒸散发量和土壤体积含水率；

构建流域月尺度水量平衡变化的目标函数，确定该目标函数的待求解参数；

对该待求解参数进行初始化，得到待求解参数的初始值；

将该待求解参数的初始值和该驱动数据输入该目标函数，得到目标函数的初始函数值；

基于该待求解参数的初始值和该初始函数值，对该待求解参数进行迭代计算，得到待求解参数的优化值，将该待求解参数的优化值输入该目标函数得到目标函数值；

在该目标函数值满足预设的迭代终止条件的情况下，根据该待求解参数的优化值估算预设时间段内流域月尺度径流值。

在其中一个实施例中，对该待求解参数进行迭代计算，得到待求解参数的优化值，包括：

计算该目标函数对该待求解参数的梯度函数；

计算该梯度函数在上一轮待求解参数的优化值下的梯度值；

计算本轮该待求解参数在该梯度值下的参数值，得到该待求解参数的优化值。

在其中一个实施例中，计算本轮该待求解参数在该梯度值下的参数值，得到该待求解参数的优化值，包括：

根据该上一轮待求解参数的优化值、该梯度值以及预设的优化精度，将所述上一轮待求解参数的优化值沿着所述梯度值的负方向进行迭代计算，得到该待求解参数的优化值。

在其中一个实施例中，该待求解参数包括多个待求解参数，共同组成参数空间；

计算该目标函数对该待求解参数的梯度函数，包括：

计算该目标函数对该参数空间中每一个待求解参数的梯度函数；

计算该梯度函数在上一轮待求解参数的优化值下的梯度值，包括：