[发明专利]径流重建方法、装置、计算机设备和存储介质

说明书

本申请涉及一种径流重建方法、装置、计算机设备和存储介质。其中，该方法包括：获取水文站实测的径流时间序列，以及获取多个位置对应的多个帕默尔干旱指数时间序列，其中，多个位置是根据水文站的位置确定的；根据实测的径流时间序列和多个帕默尔干旱指数时间序列，得到水文站对应的替代性指标；将水文站对应的替代性指标以及水文站上游的水文站的重建的径流时间序列输入径流重建模型，得到水文站对应的重建的径流时间序列。可见，该方法中一方面采用与水文站的位置相关联的帕默尔干旱指数时间序列，另一方面考虑到流域内水文站的上下游关系。由于更细致地考虑了径流重建相关数据之间的统计关系，因此提高了重建的径流时间序列的准确性。

技术领域

本申请涉及水资源处理技术领域，特别是涉及一种径流重建方法、装置、计算机设备和存储介质。

背景技术

径流是水循环中的一个重要变量。长时间序列的水文变化和相关气候变化信息是揭示径流长期演变规律和河道径流开发利用的基础。目前全球大多数地区的水文站的监测数据记录较短，最长的器测历史记录也不超过150年。因此，寻找更长的径流变化历史记录亟待解决。

传统技术中，常采用树木年轮资料延长河道径流量时间序列，由此形成了树轮水文学。树轮水文学是一门根据树木年轮指标研究过去长期水文现象的科学。树木年轮已发展成为重要的气候替代性指标，多年来广泛应用于径流重建的过程中。

然而，采用传统技术，由于树轮样点的空间分布不均，会存在重建径流准确性较低的问题。

发明内容

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

一种径流重建方法，所述方法包括：

获取水文站实测的径流时间序列，以及获取多个位置对应的多个帕默尔干旱指数时间序列，其中，所述多个位置是根据所述水文站的位置确定的；

根据所述实测的径流时间序列和所述多个帕默尔干旱指数时间序列，得到所述水文站对应的替代性指标；

将所述水文站对应的替代性指标以及所述水文站上游的水文站的重建的径流时间序列输入径流重建模型，得到所述水文站对应的重建的径流时间序列。

在其中一个实施例中，所述替代性指标包括第一替代性指标；所述根据所述实测的径流时间序列和所述多个帕默尔干旱指数时间序列，得到所述水文站对应的替代性指标，包括：

计算所述实测的径流时间序列和所述多个帕默尔干旱指数时间序列中每个帕默尔干旱指数时间序列之间的相关系数；

从所述多个帕默尔干旱指数时间序列中，选取在所述水文站的控制流域范围内，所述相关系数大于或等于系数阈值的帕默尔干旱指数时间序列作为第一替代性指标。

在其中一个实施例中，所述替代性指标还包括第二替代性指标；所述根据所述实测的径流时间序列和所述多个帕默尔干旱指数时间序列，得到所述水文站对应的替代性指标，还包括：

对所述多个帕默尔干旱指数时间序列进行主成分分析，得到多个主成分和每个所述主成分的方差贡献率；

根据所述每个主成分的方差贡献率，从所述多个主成分中，选取累计方差贡献率大于或等于贡献率阈值的主成分作为第二替代性指标，其中，所述累计方差贡献率是对所述每个主成分的方差贡献率进行累加得到的。

在其中一个实施例中，所述根据所述每个主成分的方差贡献率，从所述多个主成分中，选取累计方差贡献率大于或等于贡献率阈值的主成分作为第二替代性指标，包括：

对所述每个主成分的方差贡献率进行碎石检验，确定所述多个主成分中，累计方差贡献率大于或等于贡献率阈值的主成分并作为第二替代性指标。