一种缺资料流域径流量智能预测方法与流程

1.本发明涉及水文预测领域，尤其涉及一种缺资料流域径流量智能预测方法。


背景技术：

2.径流量预测一直是水文学研究的重点之一，可为流域水旱灾害防治、库群安全经济运行、水资源科学分配及社会可持续发展等提供重要决策支撑。其精度依赖于水文模型的构建。水文模型是对复杂水文过程进行概化的科学模型，是径流预报研究的重点及核心之一，按照建模基本原理可其分为具有物理机制的水文模型和数据驱动模型。
3.传统的降雨-径流模型是基于水文的物理过程来估算流量，可以考虑土地利用、坡度、土壤和气候的空间变异性，成功的物理水文建模通常需要准确的输入，如降水和流量数据。但是，高精度的水文和气象站不能部署在偏远地区。此外，由于资金的限制，许多发展中国家水文监测和预报信息的水文和气象观测站非常有限，且蒸发、入渗和地下水流等水文循环的组成部分难以测量。这意味着在缺资料流域的径流量模拟和预测是水文学家面临的最重要和最具挑战性的问题之一。
4.随着观测数据的增多和计算能力的提升，数据驱动模型也被应用于径流量预测。如专利cn112801416a发布了一种基于多维水文信息的lstm流域径流量预测方法。专利cn112818608a公开了一种基于改进的粒子群优化算法与支持向量机的中长期径流预报方法。专利cn111275253b公开了一种融合dnn和误差校正的径流概率性预测方法及系统。但是大多数深度学习模型只能捕获变量之间的相关性，它们无法理解因果关系，也就是说深度学习模型可以很好地拟合观测数据，但预测可能在物理上不一致，轻微的扰动也可能导致较大的变化。专利cn113379029b，实现了物理定律与深度学习模型相结合进行水质预测，但仅限于将物理定律运用到深度学习模型中，对物理机制的运用比较片面，且对于资料匮乏地区缺乏相似性分类，不利于迁移学习。


技术实现要素：

5.有鉴于此，为了可以更好地实现缺资料地区的数据和径流量预测模型的迁移使用，克服了原有无资料流域径流量预测中数据缺少的问题，本发明提出一种缺资料流域径流量智能预测方法，包括以下步骤：
6.s1、利用水文相似性原则，根据地理环境特点，对流域之间的相似度进行评价；
7.s2、根据流域之间的相似度评价结果，确定与缺资料目标流域相似的流域集，并确定缺资料目标流域过程驱动模型的参数，构建缺资料目标流域的水文过程驱动模型，并对过程驱动模型进行参数率定，确定可用的参数范围；
8.s3、使用s2中率定完成的过程驱动模型根据参数范围随机生成模型参数，输出得到径流量的模拟时间序列数据，共同构成缺资料流域的流域类型-模型参数-径流量数据矩阵，并依据该数据矩阵构建预训练数据集；
9.s4、将径流量的过程知识引入深度学习模型a，得到深度学习模型b，并将预训练数
据集标准化后输入深度学习模型b进行物理引导深度学习模型b的预训练，同时对深度学习模型b的参数进行优化，得到预训练完成后的模型；
10.s5、根据径流量历史实测数据训练s4中预训练完成后的模型，并对优化后的参数进行微调，得到物理引导的深度学习模型c；
11.s6、在s5的模型c中输入当前或历史实测的气候、径流量时间序列数据，输出得到目标时段内的径流量预测数据。
12.本发明提供的技术方案带来的有益效果是：
13.(1)本发明将物理知识植入到深度学习模型中，增加了深度学习网络的物理意义，实现了深度学习模型和水文模型的物理一致性，突出了深度学习和物理机理之间的共生整合：水文模型为数据驱动模型解决了数据短缺的问题，水文模型中的知识赋予了深度学习模型物理感知能力，这种数据和知识的耦合可提高径流量预测的准确度以及迁移鲁棒性，增强跨流域转移能力；
14.(2)相较于单一的深度学习方法，设置了基于径流量的过程知识，考虑了降雨-径流过程的物理机理，解决了深度学习方法用于无资料流域径流量预测中缺乏物理意义、精度不高、可转移性差等问题；相比于过程驱动模型，pgdl模型需要更少的边界条件和实测数据，就可以实现预测精度高、鲁棒可移植性强和推断未观察过程的智能性等优势；
15.(3)本发明针对缺资料地区建立了物理引导深度学习模型，克服了缺资料地区的数据局限性，提高了缺资料地区径流量预报能力，进而可为缺资料地区的洪水预报预警提供一种新思路。
附图说明
16.图1是本发明一种缺资料流域径流量智能预测方法的流程图；
17.图2是本发明实施例一种基于水量平衡的过程驱动efdc-lstm模型对缺资料流域的径流量进行预测的流程框图；
18.图3是本发明实施例模拟各种模型搭建的具体流程图。
具体实施方式
19.为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。
20.本实例选定长江流域的缺资料子流域为研究对象，提出一种缺资料流域径流量智能预测方法，本发明一种缺资料流域径流量智能预测方法的流程参考图1，本发明实施例一种基于水量平衡的过程驱动efdc-lstm模型对缺资料流域的径流量进行预测的流程框图参考图2，本发明实施例模拟各种模型搭建的具体流程图参考图3。
21.本实施例的方法包括以下步骤：
22.s1、利用水文相似性原则，根据地理环境特点，对流域之间的相似度进行评价。
23.具体步骤如下：
24.s11、根据流域的产汇流、气候和下垫面的因素特征，选择产汇流指标、气候指标、下垫面指标作为水文相似性评价指标；并计算所述水文相似性评价指标的相似元值，其计算公式为：
[0025][0026]
其中，当流域a、b的水文相似要素分别为要素ai与要素bi，水文指标值分别为ya与yb时，记水文相似元值为ui(ai,bi)。
[0027]
本实施例中，产汇流指标包括径流系数、流量历时曲线斜率、基流指数等；气候指标包括降水、蒸发、气温、干旱指数等；下垫面指标包括河网密度、土地利用、高程、土壤含水量等。这些指标并不是简单的罗列关系，而是各指标之间相互联系、共同作用形成的有机组合，且指标不是一成不变的，需要根据流域的统计相似、自相似、动态相似、评价目的来选择。
[0028]
s12、使用聚类分析法、层次分析法或者主成分分析法对各水文相似性评价指标进行赋权，得到各指标的权值。
[0029]
s13、根据s11中的相似元值和s12中的权值代入流域水文相似度计算公式，计算流域两两之间的水文相似度，再结合水文相似度评价表，得到流域之间的相似度评价结果，水文相似度计算公式如下：
[0030][0031]
其中，m为流域a中水文要素组成数，n为流域b水文要素组成数，k为流域a、b之间相似水文要素个数，q(ui)为相似元值，βi为水文相似指标的权重，βi∈[0,1]，且∑βi＝1。
[0032]
本实施例选用长江流域缺资料子流域的流域面积、土地利用、河网平均长度、平均地形指数、流量历时曲线斜率作为评价水文相似性的评价指标，利用流域的dem(数字高程模型)并结合arcgis软件中的calmorph模块计算出子流域的河长率、面积率等，结合多年平均降雨数据和土地利用数据计算水文相似元值以及各相似指标的权重，再计算并确定流域的水文相似性结果。
[0033]
s2、根据流域之间的相似度评价结果，确定与目标流域相似的流域集，再根据流域集的过程驱动模型参数确定目标流域的参数，构建缺资料流域的水文过程驱动模型，并对过程驱动模型进行参数率定，确定可用的参数范围。
[0034]
上述的水文过程驱动模型可以为hspf、efdc、wasp、swat、qual2k、mike等数值模拟模型。
[0035]
具体步骤如下：
[0036]
s21、根据s1中的水文相似性评价结果，确定水文资料可以移植到目标流域的流域集，将流域集的过程驱动模型参数移植到目标流域，或使用回归模型建立模型参数与流域属性之间的回归关系，将回归关系移植到目标流域，再计算得到目标流域的模型参数；
[0037]
s22、使用上述模型参数，构建该流域的过程驱动模型，输出得到径流量的模拟数据；
[0038]
s23、使用历史实测水文气候数据对过程驱动模型的参数进行率定，直到径流量模拟值与实测值的一致性最好，确定各参数的阈值范围。
[0039]
本实施例中，通过将与缺资料流域相似的流域集的efdc模型参数进行直接移植或参数平均，确定缺资料流域的模型参数，构建efdc过程驱动模型，模拟得到缺资料流域的径流量数据，并使用平均相对误差、纳什效率系数(nse)、决定系数(r2)对模拟的径流量数据
与实测数据之间的一致性进行评定，完成对模型的湍流参数、粗糙高度等参数的率定验证，确定efdc模型各参数的阈值范围。
[0040]
s3、使用s2中率定完成的过程驱动模型根据参数范围随机生成模型参数，输出得到径流量的模拟时间序列数据，共同构成缺资料流域的流域类型-模型参数-径流量数据矩阵，并依据该数据矩阵构建预训练数据集。
[0041]
s31、使用s23中率定验证后的过程驱动模型，根据各参数范围随机生成模型参数组合，模拟得到径流量数据集；
[0042]
s32、将s31生成的参数和径流量数据构成缺资料流域的流域类型-模型参数-径流量数据矩阵，并按照一定比例划分为训练集、验证集，将其作为数据驱动预训练模型的输入数据集。
[0043]
本实施例中，根据随机生成的模型参数构建efdc模型，模拟生成径流量数据，将其按照一定的比例分为lstm预训练模型的训练集、验证集和测试集。
[0044]
s4、将径流量的过程知识引入深度学习模型a，得到深度学习模型b，并将预训练数据集标准化后输入深度学习模型b进行物理引导深度学习模型b的预训练，同时对深度学习模型b的参数进行优化，得到预训练完成后的模型。其中，深度学习模型包括但不限于lstm、rnn、cnn等深度学习模型。
[0045]
s41、在理论过程中识别出径流量的物理约束，如水量平衡、能量平衡、边界条件、变量之间的单调性等，并简化为等式约束或不等式约束；
[0046]
s42、对于等式约束，一般是守恒方程或物理定律，水量平衡、能量平衡使用等式约束，如公式1，将其引入深度学习模型b中的损失函数，当违反该等式时，在每个时间步上惩罚模型b性能。
[0047]
f(x,y)＝0(3)
[0048]
对于不等式约束，一般是变量之间的单调关系，如公式10。
[0049]
本实施例中，使用relu函数作为深度学习模型b的激活函数，将公式(10)转换为具体损失值，添加到损失函数中。
[0050]
g(x,y)≥0(4)
[0051]
公式(3)、(4)表示目标变量与其他物理变量之间物理关系的一般形式，这些基于物理的方程必须满足与其他损失函数项相同的标准，即连续和可微。
[0052]
综上的物理约束的数学公式表达如下：
[0053][0054]
添加物理约束的损失函数的一般形式如公式：
[0055][0056]
其中，x为与目标变量有联系的其他物理变量；y为目标变量；为模型b对目标变量的预测值；relu函数为深度学习模型b的激活函数；为模型b的经验误差；r(f)为模型b的结构误差；λ为权重超参数；λ
phy
为物理约束损失项的权重超参数。
[0057]
本实施例中，根据水量平衡的物理定律(简化后的公式如公式(7)所示)，以及降雨和径流之间的单调关系修改lstm的损失函数，如公式(8)所示，对违反水量平衡定律的预测结果进行惩罚，实现对lstm模型的物理约束。
[0058][0059][0060][0061][0062]
其中n
t
为输入时间序列长度；y
t
为第t个时间段的径流量；r
t
为第t个时间段的降雨量；δt为时段初至时段末的间隔时间；δw为时段初至时段末的水量变化；e
t
为第t个时间段的蒸散发量；y为径流量观测值；为径流量预测值；f为lstm的模拟函数；δr为降雨量变化量；relu函数为深度学习模型的激活函数；为该模型的经验误差；r(f)为模型的结构误差；λ为权重超参数；λ
phy
为物理约束损失项的权重超参数。
[0063]
s5、根据径流量历史实测数据训练s4中预训练完成后的模型，并对优化后的参数进行微调，得到物理引导的深度学习模型c。模型c具有很强的可推广性和迁移性，可根据步骤s2中提供的不同类型的流域，对步骤s3中的过程驱动模型和步骤s4中的物理引导深度学习预训练模型b的参数进行调整，进而生成适用于不同流域的物理引导的深度学习模型。
[0064]
本实施例中，使用s4中进行物理约束后的lstm模型进行预训练，其输入为随机生成的水文、气候及时间等输入变量的模拟时间序列数据，以及efdc模型模拟的径流量训练集，输出为径流量的预训练模拟结果，通过观察预训练模拟结果与efdc模拟的验证集的rmse、mae、nse，来优化lstm模型的超参数，表现最佳的模型即为pgdl预训练模型，即模型b。
[0065]
s6、在s5的模型c中输入当前或历史实测的气候、径流量时间序列数据，输出得到目标时段内的径流量预测数据。
[0066]
本实施例中，将当前时刻实测的水文、气候等数据输入pgdl预训练模型，对pgdl预训练模型进行训练，调整优化模型的网络连接权重和偏置，得到最终的pgdl模型，并且输出得到未来预设时间的缺资料地区的径流量模拟值。
[0067]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

技术特征：
1.一种缺资料流域径流量智能预测方法，其特征在于，包括以下步骤：s1、利用水文相似性原则，根据地理环境特点，对流域之间的相似度进行评价；s2、根据流域之间的相似度评价结果，确定与缺资料目标流域相似的流域集，并确定缺资料目标流域过程驱动模型的参数，构建缺资料目标流域的水文过程驱动模型，并对过程驱动模型进行参数率定，确定可用的参数范围；s3、使用s2中率定完成的过程驱动模型根据参数范围随机生成模型参数，输出得到径流量的模拟时间序列数据，共同构成缺资料流域的流域类型-模型参数-径流量数据矩阵，并依据该数据矩阵构建预训练数据集；s4、将径流量的过程知识引入深度学习模型a，得到深度学习模型b，并将预训练数据集标准化后输入深度学习模型b进行物理引导深度学习模型b的预训练，同时对深度学习模型b的参数进行优化，得到预训练完成后的深度学习模型；s5、根据径流量历史实测数据训练s4中预训练完成后的深度学习模型，并对优化后的参数进行微调，得到物理引导的深度学习模型c；s6、在s5的模型c中输入当前或历史实测的气候、径流量时间序列数据，输出得到目标时段内的径流量预测数据。2.根据权利要求1所述的一种缺资料流域径流量智能预测方法，其特征在于，步骤s1具体为：s11、根据流域的产汇流、气候和下垫面的因素特征，选择产汇流指标、气候指标、下垫面指标作为水文相似性评价指标；并计算水文相似性评价指标的相似元值，其计算公式为：其中，y
a
为流域a的水文相似要素a
i
的水文指标值，y
b
为流域b的水文相似要素b
i
的水文指标值；s12、使用聚类分析法、层次分析法或者主成分分析法对各水文相似性评价指标进行赋权，得到各指标的权值；s13、根据s11中的相似元值和s12中的权值代入流域水文相似度计算公式，计算流域两两之间的水文相似度，再结合水文相似度评价表，得到流域之间的相似度评价结果，水文相似度计算公式如下：其中，m为流域a中水文要素组成数，n为流域b水文要素组成数，k为流域a、b之间相似水文要素个数，β
i
为水文相似指标的权重，β
i
∈[0,1]，且∑β
i
＝1。3.根据权利要求2所述的一种缺资料流域径流量智能预测方法，其特征在于，所述产汇流指标包括径流系数、流量历时曲线斜率、基流指数；所述气候指标包括降水量、潜在蒸散发量、气温、干旱指数；所述下垫面指标包括河网密度、土地利用、高程、土壤含水量。4.根据权利要求1所述的一种缺资料流域径流量智能预测方法，其特征在于，步骤s2具体为：s21、根据s1中的水文相似性评价结果，确定水文资料可以移植到目标流域的流域集，将流域集的过程驱动模型参数移植到目标流域；或使用回归模型建立模型参数与流域属性
之间的回归关系，将回归关系移植到目标流域，并计算得到目标流域的模型参数；s22、使用目标流域的模型参数，构建目标流域的过程驱动模型，输出得到径流量的模拟数据；s23、使用历史实测水文、气候数据对过程驱动模型的参数进行率定，直到径流量模拟值与实测值的一致性最好，确定各参数的阈值范围。5.根据权利要求1所述的一种缺资料流域径流量智能预测方法，其特征在于，步骤s3具体为：s31、使用率定后的过程驱动模型，根据各参数范围随机生成模型参数组合，模拟得到径流量数据集；s32、将s31生成的参数和径流量数据构成缺资料流域的流域类型-模型参数-径流量数据矩阵，并按照一定比例划分为训练集、验证集，将其作为数据驱动预训练模型的输入数据集。6.根据权利要求1所述的一种缺资料流域径流量智能预测方法，其特征在于，步骤s4具体为：s41、识别径流量的物理约束，包括水量平衡、能量平衡、边界条件、变量之间的单调性，并简化为等式约束或不等式约束；s42、水量平衡、能量平衡使用等式约束，数学公式表达如下：f(x,y)＝0(3)将其引入深度学习模型b中的损失函数，当违反该等式时，在每个时间步上惩罚模型b性能；变量之间的单调关系使用不等式约束，数学公式表达如下：g(x,y)≥0(4)并使用relu函数转换为具体损失值，添加到损失函数中；综上的物理约束的数学公式表达如下：添加物理约束的损失函数的一般形式如公式：其中，x为与目标变量有联系的其他物理变量；y为目标变量；为模型b对目标变量的预测值；relu函数为深度学习模型b的激活函数；为模型b的经验误差；r(f)为模型b的结构误差；λ为权重超参数；λ
phy
为物理约束损失项的权重超参数。

技术总结
本发明公开了一种缺资料流域径流量智能预测方法，包括：根据流域之间的相似度评价结果，确定与目标流域相似的流域集，根据流域集的过程驱动模型参数确定目标流域的参数范围，构建缺资料流域的水文过程驱动模型；基于率定验证后的过程驱动模型，生成径流量的模拟时间序列数据，共同构成缺资料流域的流域类型-模型参数-径流量数据矩阵，依据该数据矩阵构建预训练数据集，用以训练引入径流量过程知识的深度学习模型；用历史实测数据训练预训练完成后的模型，得到最终的物理引导的深度学习模型，并输出目标时段内的径流量指标预测数据。本发明解决了深度学习方法用于无资料流域径流量预测中缺乏物理意义、精度不高、可迁移性差等问题。差等问题。差等问题。


技术研发人员：谢李孜 王永桂 陈卓 李强
受保护的技术使用者：武汉祁联生态科技有限公司
技术研发日：2023.03.17
技术公布日：2023/7/19