技术特征：
1.一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，包括以下步骤：s1、在损失函数中加入对违反物理定律的预测的惩罚项，实现对深度学习模型的物理约束；s2、使用基于过程的水质模型模拟生成水质指标的模拟时间序列数据，并进行预处理，生成训练数据集和验证集；s3、使用训练集对物理定律约束的深度学习模型进行预训练，学习机理模型中人类对过程机理的理解，根据模型在验证集上的表现最优化模型的超参数，得到prpgdl预训练模型；s4、使用历史实测时间序列数据对prpgdl预训练模型进行训练，得到prpgdl模型；s5、将当前和当前时刻之前预设时间段内的实测时间序列数据输入到训练完成的prpgdl模型中，输出得到未来水质指标的预测数据。2.根据权利要求1所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，所述prpgdl模型为基于本发明提出的技术体系建立的一类水质预测模型的统称。3.根据权利要求1所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，步骤s1具体为：s11、根据所要预测的水质指标的特征选择合适的物理定律，对选定的物理定律进行概化，计算预测值对应的物理状态偏离真实值对应的物理状态的程度，计算公式概化为：pr(y
m
，y
p
，**arges)＝|s(y
m
，**arges)
‑
s(y
p
，**arges)|
‑
α#(1)其中y
m
为观测值，y
p
为预测值，**arges为除观测值和预测值之外的计算选定物理状态的变量，s(y，**arges)为计算物理状态的函数，α为确定的偏离程度的阈值；s12、将pr(y
m
，y
p
，**arges)函数集成到深度学习模型的损失函数中，作为损失函数惩罚项的一部分，损失函数修改为：loss function＝|pr(y
m
，y
p
，**arges)|
×
loss(y
m
，y
p
)#(2)其中loss(y
m
，y
p
)为深度学习模型未经修改的损失函数，y
m
为观测值，y
p
为预测值。4.根据权利要求3所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，所述物理定律包括物质守恒定律、能量守恒定律以及与水体污染物迁移转化过程有关的物理学或湖沼学的定律或规律中的一种或几种。5.根据权利要求1所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，步骤s1所述的深度学习模型包括rnn、gru、lstm、lstm
‑
seq2seq、transformer以及这些深度学习模型的常见变体中的一种或几种。6.根据权利要求1所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，s2所述的水质模型包括efdc、qual2k、swat以及基于过程的水质模型中的一种或几种。7.根据权利要求1所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，所述s2具体过程为：s21：随机生成选定的水质模型所需要的边界条件、参数以及输入变量的模拟数据；s22：将s21生成的边界条件和参数设置为水质模型的边界条件和参数，将输入变量的
模拟数据输入到水质模型中，输出得到水质模型的模拟时间序列数据；s23：重复步骤s21和s22，直到生成的模拟时间序列数据的数据量达到设定值；s24：输入变量的模拟数据与水质指标的模拟时间序列数据构成预训练数据集，将预训练数据集按照一定比例划分为训练集和验证集，对训练集、验证集进行标准化处理。8.根据权利要求7所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，步骤s21中的输入变量的模拟数据，为随机生成的所选定的水质模型所需要的输入变量的数据，包括气象、水文、地形、地理属性数据以及水质数据中的一种或几种。9.根据权利要求7所述的一种基于物理定律和过程驱动的深度学习模型的水质预测方法，其特征在于，输出的水质指标包括溶解氧、高锰酸盐指数、五日生化需氧量、化学需氧量、叶绿素浓度、氨氮、总氮、总磷、重金属离子的浓度、ph和水温中的一种或几种。
技术总结
本发明提供一种基于物理定律和过程驱动的深度学习模型的水质预测方法，根据物理定律对深度学习模型的损失函数进行修改；使用水质模型生成水质指标的模拟时间序列数据；使用模拟数据对深度学习模型进行训练，得到预训练模型；使用水质指标的历史实测数据对预训练模型进行调整优化，得到物理约束和过程驱动的深度学习模型PRPGDL；最后，基于PRPGDL模型预测未来的水质指标数据。本发明相比水质模型需要更少的边界条件和参数、有更高的预测准确度、速度和灵活性；相比深度学习模型具有更高的准确性和通用性，并且需要更少的实测数据；提供准确度更高、泛化能力和适用性更强、实测数据需求更少的水质预测方法。求更少的水质预测方法。求更少的水质预测方法。
技术研发人员：李强 王永桂
受保护的技术使用者：中国地质大学（武汉）
技术研发日：2021.04.22
技术公布日：2021/9/9