基于启发迁移复形竞争进化算法的水文模型参数率定方法技术

本发明专利技术提供基于启发迁移复形竞争进化算法的水文模型参数率定方法，包括进行数据收集，还包括以下步骤：进行复形划分；提取单纯形；引入单纯形的最优信息，生成单纯形搜索点；所述单纯形搜索点根据单纯形搜索效果对单纯形最优点的引入比例作出调整；替换单纯形最差点为单纯形搜索点，并按原始位置将单纯形样本点重新嵌入复形并按目标函数评估值升序排序复形；将各复形填充回总样本；重复执行上述步骤直至满足收敛标准；输出水文模型最优参数组

【技术实现步骤摘要】
基于启发迁移复形竞争进化算法的水文模型参数率定方法


[0001]本专利技术涉及流域水文模型的
，特别是基于启发迁移复形竞争进化算法的水文模型参数率定方法
。

技术介绍

[0002]洗牌复形进化算法是全局寻优算法，是复形竞争进化算法搭载的外部框架
。
相较于其它全局寻优算法，其特殊的洗牌方式极大地提高了各复形内样本点的空间差异性，具有极强的广度搜索能力
。
[0003]水文模型的误差响应面非常复杂，具有高度非线性
、
不连续不可导的特性，存在大量局部最优解，全局寻优算法的深度搜索方案普遍相对简单，寻优效率和解的精度较差
。
大量实验和实际应用证明，复形进化算法采用的单纯形下山搜索法的深度搜索效果优势明显，在参数率定领域应用较广
。
但其搜索过程中仅仅使用了局部最差点和形心点，并未充分利用局部可用信息
。
[0004]因此，原版的洗牌复形进化算法在复形划分时仍有一定的改进空间，可从样本点之间空间距离的角度划分复形，使复形内各样本点之间的空间距离保持足够的差异性，扩大复形内样本点遍布样本空间的范围，增强算法的广度搜索能力，提高水文模型参数率定的寻优效果
。
[0005]申请公开号为
CN106709181A
的专利技术专利申请公开了一种基于并行编程和模块法的分布式水文模型率定方法，所述方法包括：
1)
基于
MPICH
构建的并行编程环境，耦合开源分布式水文模型源程序与多目标优化算法；
2)
分别针对日径流序列和峰值序列对水文模型进行率定，并建立相应的模型框架；
3)
利用模块法合并上述模型的径流输出，形成重组的模拟径流结果
。
该方法的缺点是仅在算法计算速度上有所提高，并未对算法的机理作出改进，对解的质量没有提升
。
[0006]申请公开号为
CN113255206A
的专利技术专利申请公开了一种基于深度强化学习的水文预报模型参数率定方法，包括：根据流域特性选取合适的水文预报模型，确定模型率定的参数及参数取值范围；建立水文预报模型参数率定的强化学习模型，确定强化学习三要素即状态空间
、
动作空间及奖励函数；应用深度强化学习方法
DQN
，优化水文预报模型的率定参数
。
该方法的缺点是计算效率较低，耗时过长，需要特定的
GPU
硬件支持；易陷入局部最优点，全局性不佳；
AI
训练结果稳定性较差，多次训练结果可靠性有待提升
。

技术实现思路

[0007]为了解决上述的技术问题，本专利技术提出的基于启发迁移复形竞争进化算法的水文模型参数率定方法，在复形竞争进化过程中引入局部最优信息进行单纯形下山搜索，并根据搜索效果启发式调整局部最优信息的引入比例，优化单纯形下山搜索方向，增强算法的深度搜索能力
。
[0008]本专利技术的目的是提供基于启发迁移复形竞争进化算法的水文模型参数率定方法，
包括进行数据收集，还包括以下步骤：步骤1：进行复形划分；步骤2：提取单纯形；步骤3：引入所述单纯形的最优信息，生成单纯形搜索点；步骤4：所述单纯形搜索点根据单纯形搜索效果对单纯形最优点的引入比例作出调整；步骤5：替换单纯形最差点为单纯形搜索点，并按原始位置将单纯形样本点重新嵌入复形并按目标函数评估值升序排序复形；步骤6：将各复形填充回总样本，形成新的总样本，完成洗牌复形进化算法的一代进化；步骤7：重复执行步骤1‑
步骤6，直至满足收敛条件；步骤8：输出水文模型最优参数组
。
[0009]优选的是，所述数据包括特定的水文模型
、
水文模型适用的有资料流域所属的降雨资料和蒸发资料
、
实测径流序列资料或人工参数组生成的理想径流序列资料和水文模型需要的流域相关参数
。
[0010]在上述任一方案中优选的是，所述步骤1包括以下子步骤：步骤
11
：在参数可行空间内随机生成总样本；步骤
12
：依据各样本目标函数评估值升序排序；步骤
13
：以复形个数为间隔数间隔提取样本填充复形；步骤
14
：直至所有复形填充完毕，完成复形划分
。
[0011]在上述任一方案中优选的是，总样本数
s
为复形个数
m
与复形内样本数
n
的乘积，在参数可行空间内，通过梅森旋转算法生成0‑1均匀分布随机数，结合各参数维度允许变化上下限之差的绝对值及各参数维度确定的下限值，随机生成总样本
D = {x
i
, i = 1, ..., s}
，其中，
x
i
为第
i
个样本点
。
[0012]在上述任一方案中优选的是，所述步骤
12
包括依据各样本点
x
i
代入水文模型生成的模拟径流序列与实测或理想径流序列一同代入目标函数式计算目标函数评估值
f
i
，并对总样本
D = {x
i
, f
i
, i = 1, ..., s}
依据目标函数评估值升序排序
。
[0013]在上述任一方案中优选的是，所述步骤
13
包括对于复形
A
k
，从排序后的样本点
x
k
以复形个数
m
作为间隔数间隔提取
n
个样本点填充复形
A
k
，其中，
k=1,
…
,m。
[0014]在上述任一方案中优选的是，所述步骤
13
还包括将总样本
D
划分为
m
个复形，
D = {A
k
, k = 1, ..., m}
，复形划分表示为，其中，为总样本，为复形1，为复形2，为复形
m。
[0015]在上述任一方案中优选的是，所述步骤2包括以下子步骤：步骤
21
：根据线性概率分布非重复提取复形样本点填充单纯形；步骤
22
：单纯形内样本点依据目标函数评估值升序排序；步骤
23
：记录单纯形内最优点和最差点及除最差点外其余点的形心点
。
[0016]在上述任一方案中优选的是，所述步骤
21
包括设定单纯形样本数
q
为复形内样本数
n
的一半，根据线性概率分布非重复提取出
n
个复形样本的序列号，依据
n
个复形样本序列号提取复形样本点，填充至单纯形样本空间
X = {x本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
基于启发迁移复形竞争进化算法的水文模型参数率定方法，包括进行数据收集，其特征在于，还包括以下步骤：步骤1：进行复形划分；步骤2：提取单纯形，包括以下子步骤：步骤
21
：根据线性概率分布非重复提取复形样本点填充单纯形，设定单纯形样本数
q
为复形内样本数
n
的一半，根据线性概率分布非重复提取出
n
个复形样本的序列号，依据
n
个复形样本序列号提取复形样本点，填充至单纯形样本空间
X = {x
i
, f
i
, i = 1, ..., q}
；采用所述线性概率分布提取序列号的公式为，其中，
pos
为提取样本的序列号，
npg
为复形成员个数，
genrand_real1
为梅森旋转算法0‑1均匀分布随机数发生器，为向下取整符号；步骤
22
：单纯形内样本点依据目标函数评估值升序排序；步骤
23
：记录单纯形内最优点和最差点及除最差点外其余点的形心点；步骤3：引入所述单纯形的最优信息，生成单纯形搜索点；步骤4：所述单纯形搜索点根据单纯形搜索效果对单纯形最优点的引入比例作出调整；步骤5：替换单纯形最差点为单纯形搜索点，并按原始位置将单纯形样本点重新嵌入复形并按目标函数评估值升序排序复形；步骤6：将各复形填充回总样本，形成新的总样本，完成洗牌复形进化算法的一代进化；步骤7：重复执行步骤1‑
步骤6，直至满足收敛条件，所述收敛条件为总样本中最优点的目标函数评估值满足一定代数下的变化值小于设定阈值；步骤8：输出水文模型最优参数组
。2.
如权利要求1所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，所述步骤1包括以下子步骤：步骤
11
：在参数可行空间内随机生成总样本；步骤
12
：依据各样本目标函数评估值升序排序；步骤
13
：以复形个数为间隔数间隔提取样本填充复形；步骤
14
：直至所有复形填充完毕，完成复形划分
。3.
如权利要求2所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，总样本数
s
为复形个数
m
与复形内样本数
n
的乘积，在参数可行空间内，通过梅森旋转算法生成0‑1均匀分布随机数，结合各参数维度允许变化上下限之差的绝对值及各参数维度确定的下限值，随机生成总样本
D = {x
i
, i = 1, ..., s}
，其中，
x
i
为第
i
个样本点
。4.
如权利要求3所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，所述步骤
12
包括依据各样本点
x
i
代入水文模型生成的模拟径流序列与实测或理想径流序列一同代入目标函数式计算目标函数评估值
f
i
，并对总样本
D = {x
i
, f
i
, i = 1, ..., s}
依据目标函数评估值升序排序
。5.
如权利要求4所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其
特征在于，所述步骤
13
包括对于复形
A
k
，从排序后的样本点
x
k
以复形个数
m
作为间隔数间隔提取
n
个样本点填充复形
A
k
，其中，
k=1,
…
,m。6.
如权利要求5所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，所述步骤
13
还包括将总样本
D
划分为
m
个复形，
D = {A
k
, k = 1, ..., m}
，复形划分表示为，其中，为总样本，为复形1，为复形2，为复形
m。7.
如权利要求6所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，所述步骤3包括单纯形最差点以单纯形形心点为参考点反射或收缩，并以一定比例向单纯形最优点做线性偏移，生成初步单纯形搜索点：反射点
X
ref_ini
和收缩点
X
con_ini
，公式为
X
ref_ini=
X
w
+2(X
ce
－
X
w
)
，其中，
X
w
为单纯形最差点，
X
ce
为单纯形形心点
。8.
如权利要求7所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法，其特征在于，所述步骤3还包括所述初步单纯形搜索点以比例

【专利技术属性】
技术研发人员：阚光远，侯宇，刘荣华，刘希琛，梁珂，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：