基于HS算法优化的EEMD-RNN混凝土坝位移预测模型

范鹏飞; 祝福源; Peng-fei FAN; Fu-yuan ZHU

中国农村水利水电 ›› 2021 ›› (8) : 180-185,194.

水电建设

基于HS算法优化的EEMD-RNN混凝土坝位移预测模型

作者信息 +

Displacement Prediction Model of EEMD-RNN Concrete Dam Based on HS Algorithm Optimization

Author information +

稿件信息 +

摘要

由于混凝土坝位移的监测过程中存在噪声误差，使得智能算法预测位移时容易出现过拟合、局部最小值及收敛速度慢等缺点。为了解决上述存在的问题，首先采用集合经验模态分解法（EEMD）对位移实测值进行分解，得到不同频次的分量IMF和剩余量R，然后通过递归神经网络（RNN）对各分量和剩余项进行训练，得到环境量和时效与IMF及R的映射关系，最后通过模型得到的映射关系可求各分量和剩余项的响应预测值，并进行等权求和可得位移的预测值。在RNN训练过程中引入和声搜索算法（HS）对其进行优化，系统性去噪，以一定的概率限差函数进行扰动以得到最优解，从而优化了RNN的权重及阈值，提高了模型的健壮性。以某混凝土坝为计算案例，结果表明，HS-EEMD-RNN模型的整体拟合和预测的精度较高、过拟合程度小，与EEMD-RNN模型对比可知，HS算法可以明显提高其模型精度，并且对位移过程线的突跳值预测的精度高。

Abstract

Due to the noise error in the monitoring process of concrete dam displacement， the intelligent algorithm is prone to overall fitting， local minimum value and slow convergence speed. In order to solve the above problems， this paper uses the ensemble empirical mode decomposition （EEMD） method to decompose the measured displacement values to obtain the components IMF and residual R with different frequencies. Then， the components and residual terms are trained by recurrent neural network （RNN） to obtain the mapping relationship between environmental variables and time effect and IMF and R. Finally， the predicted response values of each component and residual term can be obtained by the mapping relation obtained from the model， and the predicted value of displacement can be obtained by summation with equal weight. In the process of RNN training， harmony search algorithm （HS） is introduced to optimize the RNN， systematically de-noising， and a certain probability limit function is used to disturb to get the optimal solution. Thus， the weights and thresholds of RNN are optimized and the robustness of the model is improved. Taking a concrete dam as an example， the results show that HS-EEMD-RNN model has higher accuracy of overall fitting and prediction， and less overfitting degree. Compared with EEMD-RNN model， HS algorithm can significantly improve the accuracy of its model， and has high precision in predicting the sudden jump value of displacement hydrograph.

关键词

位移预测 / 混凝土坝 / 集合经验模态分解 / 和声搜索法 / 递归神经网络

Key words

displacement prediction / concrete dam / ensemble empirical mode decomposition （EEMD） / harmony search method （HS） / recurrent neural network （RNN）

基金

国家重点研发计划课题(2018YFC1508603)

国家自然科学基金面上项目(51579085)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

范鹏飞 , 祝福源. 基于HS算法优化的EEMD-RNN混凝土坝位移预测模型[J].中国农村水利水电, 2021(8): 180-185,194

Peng-fei FAN , Fu-yuan ZHU. Displacement Prediction Model of EEMD-RNN Concrete Dam Based on HS Algorithm Optimization[J].China Rural Water and Hydropower, 2021(8): 180-185,194

0 引言

混凝土坝作为一种重要的坝型，在服务经济发展方面起着重要作用。随着现代监测技术的发展，混凝土坝的监测越来越系统化、智能化，监测的数据库也越来越庞大，对于数据的科学分析特别是高坝和特高坝的处理尤为重要。混凝土坝监控系统的监测项目一般包括位移、应力应变、渗流渗压、裂缝发展、坝体温度、钢筋应力应变、环境量等。其中由于混凝土坝相较于其他材料大坝具有较高的刚体结构特性，使得位移成为混凝土坝安全运行管理研究的首要指标^［1］。研究学者们已经对其进行了较为全面的研究，如位移的拟合和预测、影响位移的因素、基于位移的坝体材料参数反演、大坝位移预警指标的拟定等^［2-4］。而位移预测和预判对大坝安全管理决策的作用愈加凸显，并且对位移预测的精度和模型的适用性和可靠性的要求越来越高。因此预测混凝土坝的位移进行和提高预测精度对大坝安全就十分必要。

在大坝位移的预测方面，传统的方法有统计学，如多元回归法、逐步回归法^［1］；模糊数学预测模型及确定性预测模型等^［5］。近年来随着计算机技术及软件的发展，智能算法在大坝位移预测领域应用也越来越广泛。由单一算法向串联算法发展，而组合智能算法利用各算法的优点并且较好地克服了单一算法的过拟合问题，提高了模型的非线性映射能力、计算效率和预测精度，如徐朗等采用改进的PSO-SVM模型研究了大坝的变形^［6］；郭志扬等构建了卡尔曼滤波的GA-BP大坝变形预测模型^［7］；何启等应用PSO_GM_MC灰色神经网络-加权马尔可夫链组合算法组建了大坝变形监控模型^［8］。而在研究混凝土坝位移预测的问题上，由于混凝土坝位移受非线性变化的水位、温度、时效等因素影响，再加上监测噪声误差的存在，使得位移信号的非平稳性比较明显，无疑增大了对位移的拟合和预测的难度。为解决信号的非平稳性的难题，Huang^［9］于1998年提出了经验模态分解（empirical mode decomposition， EMD）方法，但该算法的信号包络易发生形状畸变并引起端点效应，从而使各个分量的分解不准确。Wu.Z等^［10］在EMD算法的基础上提出了集合经验模态分解（ensemble empirical mode decomposition， EEMD）可以有效地解决上述EMD算法的缺点。EEMD算法应用已较为广泛，如刘淑琼等利用EEMD-RBF神经网络构建了电离层TEC预报模型^［11］；李华等通过EEMD算法实现了轴承故障特征的提取^［12］；谭冬梅等基于EEMD-JADE较好地分离出桥梁挠度监测中的温度效应^［13］。为了提高训练拟合及预测精度，引入递归神经网络（Recurrent Neural Networks，RNN）算法。RNN神经网络是一种结构优化的神经网络，在隐含层中引入了循环结构并且在隐含元单元之间进行内部联系，在训练过程中信息可以在模型中向前或者向后传播，这可以提高算法的训练效率和精度。目前，应用该算法的研究成果颇多，如周暄焯通过RNN及其融合方法研究预测了DNA甲基化的状况^［14］；芦效峰等构建了CNN和RNN组合算法以探索自由文本击键模式持续身份的认证^［15］；范竣翔等基于RNN算法研究了空气污染时空预报模型^［16］。而大坝位移在监测过程中，监测值中包含了测量误差和系统噪声，为了进一步提高训练拟合的精度和算法的健壮性，本文采用和声搜索算法（harmony search algorithm，HS）进行去噪优化，可以有效地提高预测精度，如张磊等利用和声搜索算法可以较好地提取近红外光谱特征变量^［17］；李明等基于和声搜索算法实现了无线传感器网络多重连通覆盖^［18］；王志刚等通过和声搜索算法优化了RGV动态调度模型^［19］。

综上可知，单一算法模型在解决复杂的问题时比较困难，而对于受复杂条件影响的大坝位移具有随机性和非线性，实现较高精度的预测则更加困难。因此，为克服单一算法容易陷入误差函数的局部极小值、收敛速度慢、过度收敛等缺点，本文首先利用EEMD算法进行位移监测数据的平稳化处理，然后利用RNN神经网络较强的非线性映射能力，通过训练样本集学习建立混凝土坝环境量、时效和大坝位移的非线性关系，采用和声搜索算法对其去噪和优化，综合各算法的优势，提高位移预测精度和模型计算效率，构建HS-EEMD-RNN混凝土坝位移预测模型。

1 模型基本理论

1.1 集合经验模态分解（EEMD）

EEMD算法通过在原始时间序列中加入足够多的正态分布白噪声，再将新的时间序列进行EMD分解；然后利用白噪声均值为零的特性对各分解量作算术平均，可得到EEMD分解的IMF分量和剩余项R。分解步骤如下^［13］：

（1）首先把等长度高斯白噪声N次加入原始时间序列，则可得到若干新的时间序列：

y_{i} (t) = y (t) + n_{i} (t)

（1）

式中：

y_{i} (t)

为第i次加白噪声后的时间序列；

y (t)

为原时间序列；

n_{i} (t)

为正态分布白噪声信号。

（2）然后通过EMD算法分解添加白噪声的新时间序列，可得到本征模态函数

{w_{i}}_{j} (t)

及一个剩余项

r_{i} (t)

：

y_{i} (t) = \sum_{i = 1}^{J} {w_{i}}_{j} (t) + r_{i} (t)

（2）

式中：

{w_{i}}_{j} (t)

为第i次添加白噪声后分解所得到的第j个本征模态函数分量；

r_{i} (t)

为剩余项。。

（3）再将N组本征分量及剩余项

r (t)

求平均得到最终的IMF分量及剩余量

R (t)

。

I M F_{j} (t) = \frac{1}{N} \sum_{i = 1}^{N} w_{i j} (t)

（3）

R (t) = \frac{1}{N} \sum_{i = 1}^{N} r_{i} (t)

（4）

式中：

I M F_{j} (t)

为EEMD最终分解分量；

R (t)

为最终剩余量。

（4）最后，原始时间序列信号可分解为：

y (t) = \sum_{j = 1}^{J} I M F_{j} (t) + R (t)

（5）

1.2 递归神经网络（RNN）

RNN神经网络的优点在于一次性激活不同的神经单元，即当前时刻的隐含层节点能够反应前一时刻隐含层节点的输出内容，同时作为输出层的输入以及前一个时间点的隐含层的节点输入，并在处理该时刻的信息之后又导向下一时间点的隐含层节点^［14］。RNN较之传统的神经网络，在训练拟合和预测位移的过程中，可以提高收敛速度和训练精度，其在时间维度上展开后结构如图1。

图1 RNN在时间维度上展开的结构

Fig.1 The structure of RNN in time dimension

Full size|PPT slide

由图1可知，RNN在时间维度是一个不断往后递归的过程，之前时刻的结果对下一个时刻的训练有着直接的影响，但因为在每个时间步上的权值一样，所以在相同时刻不同单元会共享同一组权值。对于长度为T的序列x，RNN的输入层大小为I，隐含层大小为H，输出层大小为K，可得3个权重矩阵U、V和W的维度^［15］：

\{\begin{matrix} U \in R^{I \times H} \\ W \in R^{H \times H} \\ V \in R^{H \times K} \end{matrix}

（6）

定义

x^{t}

为序列中第t项输入；

a^{t}

为第t项隐含层的输入；

b^{t}

为对

a^{t}

进行非线性激活，即神经网络的输出。因此，

a^{t}

是由输入层

x^{t}

与上一层隐含层的输出

b^{t - 1}

共同决定：

a_{h}^{t} = \sum_{i} w_{i h} x_{i}^{t} + \sum_{h^{'}} w_{h^{'} h} b_{h^{'}}^{t - 1}

（7）

b_{h}^{t} = f (a_{h}^{t})

（8）

大坝位移监测序列由t=1开始，所以b ₀ = 0，然后把隐含层的信息传导输出层，则输出层的输出结果：

a_{k}^{t} = \sum_{h} w_{h k} b_{h}^{t}

（9）

y_{k}^{t} = \frac{e^{a_{k}^{t}}}{\sum_{j} e^{a_{j}^{t}}}

（10）

在t时刻，普通神经网络的残差为

δ_{k}^{t} = y_{k}^{t} - z_{k}^{t}

，而RNN由于前向传导时隐含层需要接受上一个时刻隐含层的信号，而且在反向传导时还需要接受下一个时刻隐含层的反馈，则RNN输出层的残差项：

δ_{h}^{t} = f^{'} (a_{h}^{t}) (\sum_{k} δ_{k}^{t} w_{h k} + \sum_{h^{'}} δ_{h^{'}}^{t + 1} w_{h h^{'}})

（11）

当位移序列长度为T时，则残差

δ^{T + 1}

均为0，同时整个网络只有一套参数U、V、W，则时刻t的倒数：

U : \frac{\partial O}{\partial w_{i h}} = \frac{\partial O}{\partial a_{h}^{t}} \frac{\partial a_{h}^{t}}{\partial w_{i h}} = δ_{h}^{t} x_{i}^{t}

（12）

V : \frac{\partial O}{\partial w_{h k}} = \frac{\partial O}{\partial a_{k}^{t}} \frac{\partial a_{k}^{t}}{\partial w_{h k}} = δ_{k}^{t} b_{h}^{t}

（13）

W : \frac{\partial O}{\partial w_{h^{'} h}} = \frac{\partial O}{\partial a_{h}^{t}} \frac{\partial a_{h}^{t}}{\partial w_{h^{'} h}} = δ_{k}^{t} b_{h^{'}}^{t}

（14）

则可以写成统一的形式（假定对输入层有

x_{i}^{t} = a_{i}^{t} = b_{i}^{t}

）：

\frac{\partial O}{\partial w_{h i j}} = \frac{\partial O}{\partial a_{h}^{t}} \frac{\partial a_{j}^{t}}{\partial w_{i j}} = δ_{j}^{t} b_{i}^{t}

（15）

最终，对时刻t=1，2，…，T，由RNN神经网络的递归性，将其求和可得RNN网络关于权重参数的导数：

\frac{\partial O}{\partial w_{i j}} = \frac{\partial O}{\partial a_{h}^{t}} \frac{\partial a_{j}^{t}}{\partial w_{i j}} = \sum_{t} δ_{j}^{t} b_{i}^{t}

（16）

1.3 和声搜索算法

对于EEMD分解得到的位移分量，其中高频分量主要由监测误差噪声引起的，尤其是IMF1分量，所以对其系统性去噪以提高模型的精度就很有必要。和声搜索算法的基本思想是模仿乐队演奏的和声调谐过程的元启发式优化算法^［20］，其优化RNN神经网络的实现过程如下^［21］：

步骤一：初始化相关算法参数。

和声搜索算法参数有：和声记忆库容量HMS、和声记忆库取值概率HMCR、微调扰动概率PAR、微调步长bw、迭代次数BT。HMS决定了HS算法的全局搜索能力，但过大的HMS值会影响收敛到最优解的速度，HMCR值大有利于算法的局部搜索，算法参数的具体数值可由仿真试验确定^［22］。

步骤二：定义目标函数并初始化和声记忆库。

定义目标函数：期望输出（实测值）与网络输出（拟合值）的平均相对误差F，其中F值越小表明优化程度越好。

F = \frac{1}{H} \sum_{i = 1}^{H} \frac{|y^{'} (i) - y (i)|}{|y (i)|}

（17）

式中：H为训练集样本总数；

y^{'} (i)

、

y (i)

分别为第i个样本的实测值与拟合值。

构造HMS个随机初始和声

x^{1}

、

x^{2}

、…、

x^{H M S}

添加和声记忆库，记忆库中的和声x_i 对应RNN神经网络的一组权值和阈值。和声记忆库为：

H M = [\begin{array}{l} x^{1} \\ x^{2} \\ ⋮ \\ x^{H M S} \end{array} |\begin{array}{l} f (x^{1}) \\ f (x^{2}) \\ ⋮ \\ f (x^{H M S}) \end{array} \begin{array}{l}  \end{array}] = [\begin{array}{l} x_{1}^{1} x_{2}^{1} \dots x_{N}^{1} \\ x_{1}^{2} x_{2}^{2} \dots x_{N}^{2} \\ ⋮ ⋮ ⋱ ⋮ \\ x_{1}^{H M S} x_{N}^{H M S} \dots x_{N}^{H M S} \end{array} |\begin{array}{l} f (x^{1}) \\ f (x^{2}) \\ ⋮ \\ f (x^{H M S}) \end{array} \begin{array}{l}  \end{array}]

（18）

式中：

x_{j}^{i}

为RNN网络中第i组和声的第j个权值或阈值；f（xⁱ ）为第i组目标函数值，即第i组网络输出与期望输出之间的平均相对误差。

步骤三：新和声的创作。

创作新和声，即构造新的解向量。根据学习和声记忆库、音调微调，随机产生新和声的音调

x_{j}^{'}

（j= 1，2，…，N）：

x_{j}^{'} = \{\begin{array}{l} x_{j}^{'} \in (x_{j}^{1}, x_{j}^{2}, \dots, x_{j}^{H M S}) r_{1} < H M C R \\ x^{L} + r_{1} (x^{U} - X^{L}) r_{1} \geq H M C R \end{array}

（19）

式中：r ₁为［0，1］内均匀分布的随机数；

x_{j}^{'} \in (x_{j}^{1}, x_{j}^{2}, \dots, x_{j}^{H M S})

为在第j列HMS个分量中随机选定的分量；

x^{U}

、

X^{L}

为第j列分量的取值上下限；式（19）为对任意新音调

x_{j}^{'}

，会以HMCR的概率从记忆库中选取，以1-HMCR概率从其对应的取值区间内随机取值^［23］。如果

x_{j}^{'}

是和声记忆库中的解分量，则对其微调为：

x_{j}^{'} = \{\begin{array}{l} x_{j}^{'} + r_{3} b w r_{2} < P A R \\ x_{j}^{'} r_{2} \geq P A R \end{array}

（20）

式中：

r_{2}

、

r_{3}

分别为［0，1］和［-1，1］内均匀分布的随机数。

步骤四：更新记忆库。

如果新和声解向量对应的目标函数值F优于记忆库内的最差解，则用新的和声解替换最差解，更新记忆库。

步骤五：结束循环。

循环步骤三、四，直至满足终止准则或达到最大迭代次数。

2 构建模型及精度评价

2.1 组合模型建模步骤

综合以上算法的优点，基于HS优化的EEMD-RNN组合模型的详细步骤如图2，建模具体流程如下：

图2 基于HS优化的EEMD-RNN组合模型混凝土大坝预测流程图

Fig.2 Flow chart of concrete dam prediction based on HS optimized EEMD-RNN combined model

Full size|PPT slide

（1）确定大坝位移的研究测点，通过EEMD对测值进行分解，得到若干组IMF分量和剩余项R。

（2）应用RNN神经网络算法对IMF分量和剩余项R进行训练拟合，输入为环境量和时效变量，输出为大坝位移。在训练过程中通过和声搜索算法HS优化RNN的阈值和权值，直到得到最优解。

（3）通过训练最优的拟合映射关系，计算预测集所对应的IMF分量和剩余项R。

（4）将各个分量位移IMF和剩余项R的预测值等权求和即为最终的预测结果。

2.2 模型精度评定

为了评价模型拟合及预测的精度，利用平均绝对误差（MAE）、均方根误差（RMSE）和平均绝对百分比误差（MAPE）进行模型的精度评定，分别为：

M A E = \frac{1}{n} \sum_{t = 1}^{n} |Y_{t} - {\overset{\land}{Y}}_{t}|

（21）

R M S E = \sqrt[]{{\frac{1}{n} \sum_{t = 1}^{n} (Y_{t} - {\overset{\land}{Y}}_{t})}^{2}}

（22）

M A P E = \frac{1}{n} \sum_{t = 1}^{n} |\frac{Y_{t} - {\overset{\land}{Y}}_{t}}{Y_{t}}|

（23）

式中：

Y_{t}

为混凝土坝位移实测值；

{\overset{\land}{Y}}_{t}

为模型计算量；n为监测数据总个数；t为监测数据对应的时间。

3 实例计算

某重力坝位于中国东南省份，坝顶高程634.40 m，最大坝高72.4 m，坝体上游面垂直。大坝共分为9个坝段，坝顶长206 m，坝顶宽7.5 m，其中4号和5号坝段是溢流坝段，其余为非溢流坝段。非溢流坝段下游面坝坡比为1∶0.72。其中该坝采用引张线监测坝体的水平位移，共布置1条，固定端布置在坝右01+107.025处，导向端布置在坝左0+93.50处，全长200.75 m，每个坝段布置一个测点，共9个测点，测点名称EX1~EX9，测点布置图如图3。

图3 引张线测点布置图

Fig.3 Layout of measuring points of tension line

Full size|PPT slide

由图3可知，该坝引张线系统的EX5测点位于大坝的中间坝段，监测大坝的顺河向水平位移，其位移变化具有典型性，所以选取测点EX5为算例研究对象，且EX5测点可靠，监测频次为1 次/d，采用自动化监测和人工校测，监测数据连续可靠。

通过基于HS优化的EEMD-RNN混凝土坝位移预测模型对该坝EX5测点所监测的位移进行计算分析。位移监测资料选取2016年6月2日-2017年8月20日，共445组，其中2016年6月2日-2017年6月30日的前395组作为模型训练组，2017年7月1日-2017年8月20日后50组为模型预测组。

EEMD分解参数有两个，高斯白噪声标准差Nstd一般为0.01~0.4，添加噪声次数NE一般为50或100，为尽可能得到较好的分解结果，综合Nstd取0.01，NE取100。模型的训练集EEMD分解结果如图4，其中7组IMF分量及1组剩余量R，IMF1分量、IMF2分量和IMF3分量为高频分量，其余IMF分量为低频。对于高频分量，主要由监测误差噪声引起的，尤其是IMF1分量。所以系统性地去噪以提高模型精度就很必要。

图4 大坝EX5测点位移EEMD分解结果

Fig.4 EEMD decomposition results of displacement of EX5 measuring point of dam

Full size|PPT slide

为了降低测值中的噪声误差对训练结果的干扰，通过和声搜索算法优化RNN训练过程，其中和声记忆库大小HMS取100，搜索范围为［-1，1］，和声记忆库保留概率HMCR取0.85，音调微调概率PAR取0.2，微调步长bw取0.2，最大迭代次数NT取100。在HS优化RNN训练结束后，对各分量IMF和剩余量R的拟合值进行等权求和得到位移的拟合值，拟合结果如图5所示。并且与EEMD-RNN模型的拟合结果进行对比，拟合的精度见表1所示。由图5和表1可知，HS算法可以有效提高训练拟合精度，HS-EEMD-RNN模型的拟合精度明显高于EEMD-RNN模型。

图5 两种模型的拟合结果

Fig.5 Fitting results of the two models

Full size|PPT slide

表1 两种模型的拟合精度对比

Tab.1 Comparison of fitting accuracy of two models

模型	R	MAE/mm	RMSE/mm	MAPE
EEMD-RNN	0.986 5	0.197	0.207	0.107
HS-EEMD-RNN	0.991 8	0.074	0.086	0.041

由HS-EEMD-RNN训练好的网络进行各个分量位移IMF和剩余项R的预测值，等权求和即为最终的大坝位移预测结果。同时计算了EEMD-RNN的位移预测值，位移预测结果如图6。利用均方根误差（MAE）、平均绝对百分比误差（RMSE）、相对误差和绝对误差（MAPE），3种指标对两种模型的预测情况评价，具体对比结果如表2。由图6和表2的对比结果表明，HS-EEMD-RNN模型的预测精度要优于EEMD-RNN模型；而且预测精度高，能够满足工程要求。

图6 EEMD-RNN模型和HS-EEMD-RNN模型对总位移的预测

Fig.6 Prediction of total displacement by EEMD-RNN model and HS-EEMD-RNN model

Full size|PPT slide

表2 EEMD-RNN和HS-EEMD-RNN模型预测值误差对比

Tab.2 Comparison of prediction error between EEMD-RNN model and HS-EEMD-RNN model

模型

平均绝对误差/mm

均方根误差

/mm

平均绝对

百分比误差

EEMD-RNN

0.221

0.250

0.123

HS-EEMD-RNN

0.088

0.101

0.048

其中位移变化过程中的突跳点一般是大坝安全监控工作关注的重点，而以往模型的过拟合及预测误差大的情况在突跳点位置发生的较多，所以选取EX5位移过程线突跳程度最大的10个测值与模型的预测值进行对比，以更好地评估模型的预测精度，位移突跳点预测值对比结果如表3。由表3可知，HS-EEMD-RNN的平均绝对百分比误差为4.93%，比EEMD-RNN的14.89%明显小得多。

表3 位移突跳点预测值对比

Tab.3 Comparison of predicted values of displacement jump points

日期	位移实测值/ mm	位移预测值/mm		绝对百分比误差/%
日期	位移实测值/ mm	HS-EEMD-RNN	EEMD-RNN	HS-EEMD-RNN	EEMD-RNN
2017-07-07	-1.353	-1.348 28	-1.408 28	0.34	4.10
2017-07-09	-1.280	-1.385 30	-1.585 30	8.20	23.82
2017-07-16	-1.381	-1.406 81	-1.466 81	1.87	6.22
2017-07-24	-2.516	-2.610 65	-2.310 65	3.76	8.17
2017-07-28	-3.090	-3.054 94	-2.614 94	1.12	15.36
2017-07-30	-2.907	-2.801 51	-2.651 51	3.63	8.79
2017-08-03	-0.941	-1.098 20	-1.338 20	16.71	42.22
2017-08-10	-3.176	-3.105 14	-2.985 14	2.24	6.02
2017-08-12	-2.439	-2.583 02	-2.683 02	5.91	10.01
2017-08-19	-1.868	-1.970 15	-2.320 15	5.48	24.22
平均值	-	-	-	4.93	14.89

4 结论

本文以提高大坝位移的预测精度为切入点，构建了以HS优化的EEMD-RNN的位移预测模型，主要研究成果如下。

（1）EEMD算法可以有效地处理由于非线性随机变量导致的非平稳位移监测序列；RNN神经网络能够较好地实现位移样本的训练拟合及预测；HS算法可以有效地进行系统性去噪，并能明显提高模型的健壮性及计算精度；HS-EEMD-RNN模型综合上述算法优点能够实现对混凝土坝位移的预测。

（2）通过算例可知，HS-EEMD-RNN模型的误差小，精度高，且对突跳值有较高的预测精度，是一种可靠的混凝土坝位移预测模型。该模型可以拓展到其他结构位移的预测及监控，对工程管理决策提供新的途径。

（3）HS-EEMD-RNN的位移预测模型主要适用于训练集为大样本的研究，为了提高该模型的适应性可以进一步研究训练集为小样本的预测精度及稳定性。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	吴中如. 水工建筑物安全监控理论及其应用［M］.北京：高等教育出版社，2003：2-4. 本文引用 [2]

2	周健，仲静文，郝利朋，等.基于马尔科夫链残差修正的EMD-MELM混凝土坝位移预测模型［J］.水利水电技术，2020，51（6）：40-48. 本文引用 [1]

3	彭圣军. 混凝土坝安全监控模型数值优化及变位预警指标研究［D］.南昌：南昌大学，2014.

4	顾冲时，吴中如.混凝土坝空间位移场的确定性模型反演分析法［J］.河海大学学报，1994（2）：82-87. 本文引用 [1]

5	吴中如. 混凝土坝安全监控的确定性模型及混合模型［J］. 水利学报， 1989（5）：64-70. 本文引用 [1]

6	徐朗，蔡德所.改进的PSO-SVM模型在大坝变形预测中的应用［J］.中国农村水利水电，2018（3）：120-123，128. 本文引用 [1]

7	郭志扬，王建，黄庆.基于卡尔曼滤波的GA-BP模型在大坝变形预测中的应用［J］.中国农村水利水电，2016（12）：113-116. 本文引用 [1]

8	何启，戴波.基于灰色神经网络-加权马尔可夫链的大坝变形监控模型及预报研究［J］.中国农村水利水电，2016（10）：146-150，155. 本文引用 [1]

9	HUANG N E， ZHENG Shen. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc.R.Soc.Lond，1998（454）：903-995. 本文引用 [1]

10	WU Z， HUANG N E. Ensemble empirical mode decomposition： a noise-assisted data analysis method［J］. Advances in Adaptive Data Analysis， 2009，1（1）：1-41. 本文引用 [1]

11	刘淑琼，高鑫，李长春.EEMD-RBF神经网络的电离层TEC预报模型［J］.测绘工程，2020，29（3）：15-19. 本文引用 [1]

12	李华，刘韬，伍星，等.EEMD和优化的频带熵应用于轴承故障特征提取［J］.振动工程学报，2020，33（2）：414-423. 本文引用 [1]

13	谭冬梅，刘晓飞，姚欢，等.基于EEMD-JADE的桥梁挠度监测中温度效应的分离［J］.土木与环境工程学报（中英文），2020，42（3）：90-99. 本文引用 [2]

14	周暄焯. 基于RNN及其融合方法的DNA甲基化预测模型研究［D］.成都：电子科技大学，2020. 本文引用 [2]

15	芦效峰，张胜飞，伊胜伟. 基于CNN和RNN的自由文本击键模式持续身份认证［J］.清华大学学报（自然科学版）， 2018，58（12）： 1 072-1 078. 本文引用 [2]

16	范竣翔，李琦，朱亚杰，等.基于RNN的空气污染时空预报模型研究［J］.测绘科学，2017，42（7）：76-83，120. 本文引用 [1]

17	张磊，丁香乾，宫会丽，等.改进和声搜索算法的近红外光谱特征变量选择［J］.光谱学与光谱分析，2020，40（6）：1 869-1 875. 本文引用 [1]

18	李明，胡江平.基于和声搜索算法的无线传感器网络多重连通覆盖［J］.传感技术学报，2020，33（2）：272-278. 本文引用 [1]

19	王志刚，郭钰卓，周启昕，等.基于和声搜索算法的RGV动态调度优化模型［J］.海南大学学报（自然科学版），2019，37（2）：112-118. 本文引用 [1]

20	YADAY P， KUMAR R， PANDA S K， et al. An intelligent tuned harmony search algorithm for optimisation［J］. Information Sciences， 2012，196：47-72. 本文引用 [1]

21	SUN P M， BAO T F， GU C S， et al. Parameter sensitivity and inversion analysis of a concrete faced rock-fill dam based on HS-BPNN algorithm［J］. Sci China Tech Sci. 2016，59（9）：1 442-1 451. 本文引用 [1]

22	ZHANG L， XU Y， XU G， et al. A catfish effect inspired harmony search algorithm for optimization［J］. International Journal of Nonlinear Sciences & Numerical Simulation， 2013，14（6）：413-422. 本文引用 [1]

DAS S， MUKHOPADHYAY

， ROY

， et al. Exploratory power of the harmony search algorithm： analysis and improvements for global numerical optimization［J］. IEEE Transactions on Systems Man & Cybernetics Part B Cybernetics A Publication of the IEEE Systems Man & Cybernetics Society， 2011，41（1）：89-106.

本文引用 [1]

PDF(1724 KB)

2371

访问

引用

详细情况

段落导航

收稿日期	出版日期
2020-10-29	2021-08-15
发布日期
2021-09-24

选择文件类型/文献管理软件名称

选择包含的内容

摘要

Abstract

关键词

Key words

基金

引用本文

0 引言

1 模型基本理论

1.1 集合经验模态分解（EEMD）

1.2 递归神经网络（RNN）

图1 RNN在时间维度上展开的结构

1.3 和声搜索算法

2 构建模型及精度评价

2.1 组合模型建模步骤

图2 基于HS优化的EEMD-RNN组合模型混凝土大坝预测流程图

2.2 模型精度评定

3 实例计算

图3 引张线测点布置图

图4 大坝EX5测点位移EEMD分解结果

图5 两种模型的拟合结果

表1 两种模型的拟合精度对比

图6 EEMD-RNN模型和HS-EEMD-RNN模型对总位移的预测

表2 EEMD-RNN和HS-EEMD-RNN模型预测值误差对比

表3 位移突跳点预测值对比

4 结论

{{custom_sec.title}}

{{custom_sec.title}}

参考文献

选择文件类型/文献管理软件名称

选择包含的内容

摘要

Abstract

关键词

Key words

基金

引用本文

0 引 言

1 模型基本理论

1.1 集合经验模态分解（EEMD）

1.2 递归神经网络（RNN）

图1 RNN在时间维度上展开的结构

1.3 和声搜索算法

2 构建模型及精度评价

2.1 组合模型建模步骤

图2 基于HS优化的EEMD-RNN组合模型混凝土大坝预测流程图

2.2 模型精度评定

3 实例计算

图3 引张线测点布置图

图4 大坝EX5测点位移EEMD分解结果

图5 两种模型的拟合结果

表1 两种模型的拟合精度对比

图6 EEMD-RNN模型和HS-EEMD-RNN模型对总位移的预测

表2 EEMD-RNN和HS-EEMD-RNN模型预测值误差对比

表3 位移突跳点预测值对比

4 结 论

{{custom_sec.title}}

{{custom_sec.title}}

参考文献

0 引言

4 结论