一种水文气象序列样本增扩方法、装置和设备制造方法及图纸

本发明专利技术公开了一种水文气象序列样本增扩方法、装置和设备，对原始水文气象序列进行正态化处理，得到服从正态分布的水文气象序列；消除服从正态分布的水文气象序列的自相关性，得到服从同一分布的正态水文气象序列；采用改进的基于LSTM的Seq2seq模型对服从同一分布的正态水文气象序列的两端进行端点延拓，得到端点延拓后的序列；采用改进的GR

【技术实现步骤摘要】
一种水文气象序列样本增扩方法、装置和设备


[0001]本专利技术涉及水文频率计算领域，具体涉及一种水文气象序列样本增扩方法、装置和设备。

技术介绍

[0002]水文气象频率计算主要利用一定时期内的实测水文气象序列样本计算某种水文变量的出现次数，从而预测出某种频率下该变量对应的设计值，进而在防洪保障、水利工程运行、工程水文设计等方面提供重要的依据。近年来由于水资源短缺、气候变化及人类活动的加剧，各流域的下垫面条件显著改变，使得水文气象序列样本的代表性及一致性不能得到保障，从而影响后续的频率计算及工程水文设计的准确度，甚至影响水利工程的建设及安全。
[0003]研究表明进行合理的样本增扩可以有效提高水文气象序列代表性，目前比较常用的样本增扩方法有蒙特卡洛法、支持向量机和Bootstrap方法。其中，蒙特卡洛法原理简单，但在极小子样的计算中误差比较大；支持向量机方法在处理小子样数据回归问题上，能克服模型对样本容量的敏感问题，优势相对明显，但其准确度受参数的影响较大，而参数的优化过程非常繁琐；Bootstrap法通过重复抽样增加样本数量，无需预设样本分布，在工程实际中被广泛应用，但由于其生成的增扩样本的取值范围局限于原始小样本，使得增扩样本分布一定程度偏离原样本的真实分布，样本代表性有限，故而难以满足实际工程中的水文设计需要。因此，亟需提出一种更具有代表性的水文气象序列样本增扩法，从而提高水文变量设计值精度，为后续工程水文设计奠定基础。

技术实现思路

[0004]针对现有技术中存在的问题，本专利技术提供了一种水文气象序列样本增扩方法、装置和设备，能够使得水文气象序列样本更具有代表性，进而有利于显著提高水文气象频率计算结果的准确性。
[0005]为了解决上述技术问题，本专利技术通过以下技术方案予以实现：
[0006]一种水文气象序列样本增扩方法，包括：
[0007]获取原始水文气象序列；
[0008]判断所述原始水文气象序列是否服从正态分布，若不服从，则进行正态化处理，得到服从正态分布的水文气象序列；
[0009]判断服从正态分布的水文气象序列是否服从同一分布，若不服从，则消除服从正态分布的水文气象序列的自相关性，得到服从同一分布的正态水文气象序列；
[0010]采用改进的基于LSTM的Seq2seq模型对所述服从同一分布的正态水文气象序列的两端进行端点延拓，得到端点延拓后的序列；
[0011]采用改进的GR
‑
CEEMD分解法对所述端点延拓后的序列进行分解，得到多组波形平稳且趋势项单调的新序列；
[0012]采用随机抽样与遗传算法的启发式概念对所述新序列进行随机交叉组合，得到增扩后的水文气象序列样本。
[0013]进一步地，所述采用改进的基于LSTM的Seq2seq模型对所述服从同一分布的正态水文气象序列的两端进行端点延拓，得到端点延拓后的序列，包括：
[0014]设定所述服从同一分布的正态水文气象序列包含n组数据，将其按照时间顺序输入由单层LSTM单元组成的编码器，生成不同时刻的隐藏层状态{h1，h2，
…
，h
n
}，每个时刻生成的隐藏层状态h
t
均与前一时刻的隐藏层状态h
t
‑1和t时刻对应的服从同一分布的正态水文气象序列x
t
有关；
[0015]经过n个时刻的运算，生成最终状态h
n
，存储在中间向量c中；
[0016]将中间向量c输入到Seq2seq模型的编码器中进行不定长时间序列的运算，得到包含序列预测信息的状态值；
[0017]通过Seq2seq的解码器对所述包含序列预测信息的状态值进行解码，得到序列的预测值；
[0018]将所述预测值中距离所述服从同一分布的正态水文气象序列两端最近的极值点添加到所述服从同一分布的正态水文气象序列两端，得到端点延拓后的序列。
[0019]进一步地，所述采用改进的GR
‑
CEEMD分解法对所述端点延拓后的序列进行分解，得到多组波形平稳且趋势项单调的新序列，包括：
[0020]使用基于GR停止准则的EMD分解方法对所述端点延拓后的序列进行m次分解，得到m组第一IMF分量和m组第一单调残余趋势项，每组包含多个第一IMF分量和一个第一单调残余趋势项；
[0021]将各组的第一IMF分量分别叠加，生成叠加后的序列D1(t)；
[0022]对所述叠加后的序列D1(t)进行完全经验模态分解，得到m组第二IMF分量和m组第二单调残余趋势项，每组包含多个第二IMF分量和一个第二单调残余趋势项；
[0023]将m组第一单调残余趋势项与m组第二IMF分量进行组合，得到多组波形平稳且趋势项单调的新序列D2(t)。
[0024]进一步地，所述采用随机抽样与遗传算法的启发式概念对所述新序列进行随机交叉组合，得到增扩后的水文气象序列样本，包括：
[0025]将所述新序列D2(t)视为生物种群中的基因片段，设定所述新序列包含q组IMF1分量、IMF2分量、
…
、IMF
k
分量与q组单调残余趋势项，则相当于所述新序列由q个染色体组成，每个染色体包含k+1个基因；
[0026]利用随机抽样原理，每次从q个染色体中随机抽取两个作为父代，再结合交叉算法，对父代染色体上的k+1个基因进行随机交叉并重组；
[0027]在随机交叉过程中，各组中处于同一位置的分量不进行组合，但与其余分量中的任一分量进行组合，经过多次循环得到增扩后的水文气象序列样本。
[0028]进一步地，采用KS检验法判断所述原始水文气象序列是否服从正态分布，包括：
[0029]将包含n组数据的原始水文气象序列进行排序后记作M1<M2<
…
<M
i
<
…
<M
n
，利用统计量T计算原始水文气象序列是否服从正态分布，计算公式如下：
[0030][0031]T＝max|G(x)
‑
F
n
(x)|，
[0032]式中，F
n
(x)表示原始水文气象序列的经验累积分布函数，G(x)表示正态分布的累积分布函数；
[0033]在预设显著水平α下，当计算得到的T大于此时的临界值T(n，α)，且渐进显著性P值大于α时，原始水文气象序列满足正态分布。
[0034]进一步地，采用贝叶斯检验法判断服从正态分布的水文气象序列是否服从同一分布。
[0035]进一步地，所述得到增扩后的水文气象序列样本后，还包括：
[0036]若增扩后的水文气象序列样本中存在异常值，则去除异常值。
[0037]进一步地，所述消除所述服从正态分布的水文气象序列的自相关性，包括：
[0038]采用去趋势或本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种水文气象序列样本增扩方法，其特征在于，包括：获取原始水文气象序列；判断所述原始水文气象序列是否服从正态分布，若不服从，则进行正态化处理，得到服从正态分布的水文气象序列；判断服从正态分布的水文气象序列是否服从同一分布，若不服从，则消除服从正态分布的水文气象序列的自相关性，得到服从同一分布的正态水文气象序列；采用改进的基于LSTM的Seq2seq模型对所述服从同一分布的正态水文气象序列的两端进行端点延拓，得到端点延拓后的序列；采用改进的GR
‑
CEEMD分解法对所述端点延拓后的序列进行分解，得到多组波形平稳且趋势项单调的新序列；采用随机抽样与遗传算法的启发式概念对所述新序列进行随机交叉组合，得到增扩后的水文气象序列样本。2.根据权利要求1所述的一种水文气象序列样本增扩方法，其特征在于，所述采用改进的基于LSTM的Seq2seq模型对所述服从同一分布的正态水文气象序列的两端进行端点延拓，得到端点延拓后的序列，包括：设定所述服从同一分布的正态水文气象序列包含n组数据，将其按照时间顺序输入由单层LSTM单元组成的编码器，生成不同时刻的隐藏层状态{h1，h2，
…
，h
n
}，每个时刻生成的隐藏层状态h
t
均与前一时刻的隐藏层状态h
t
‑1和t时刻对应的服从同一分布的正态水文气象序列x
t
有关；经过n个时刻的运算，生成最终状态h
n
，存储在中间向量c中；将中间向量c输入到Seq2seq模型的编码器中进行不定长时间序列的运算，得到包含序列预测信息的状态值；通过Seq2seq的解码器对所述包含序列预测信息的状态值进行解码，得到序列的预测值；将所述预测值中距离所述服从同一分布的正态水文气象序列两端最近的极值点添加到所述服从同一分布的正态水文气象序列两端，得到端点延拓后的序列。3.根据权利要求1所述的一种水文气象序列样本增扩方法，其特征在于，所述采用改进的GR
‑
CEEMD分解法对所述端点延拓后的序列进行分解，得到多组波形平稳且趋势项单调的新序列，包括：使用基于GR停止准则的EMD分解方法对所述端点延拓后的序列进行m次分解，得到m组第一IMF分量和m组第一单调残余趋势项，每组包含多个第一IMF分量和一个第一单调残余趋势项；将各组的第一IMF分量分别叠加，生成叠加后的序列D1(t)；对所述叠加后的序列D1(t)进行完全经验模态分解，得到m组第二IMF分量和m组第二单调残余趋势项，每组包含多个第二IMF分量和一个第二单调残余趋势项；将m组第一单调残余趋势项与m组第二IMF分量进行组合，得到多组波形平稳且趋势项单调的新序列D2(t)。4.根据权利要求3所述的一种水文气象序列样本增扩方法，其特征在于，所述采用随机抽样与遗传算法的启发式概念对所述新序列进行随机交叉组合，得到增扩后的水文气象序
列样本，包括：将所述新序列D2(t)视为生物种群中的基因片段，设定所述新序列包含q组IMF1分量、IMF2分量、
...

【专利技术属性】
技术研发人员：兰甜，冶兆霞，任冲锋，李娇娇，李同方，
申请(专利权)人：长安大学，
类型：发明
国别省市：