基于MLR-SSA-GRU的混凝土坝裂缝开度组合预测模型

1、基于MLR-SSA-GRU的混凝土坝裂缝开度组合预测模型混凝土抗拉强度低，在复杂的外力环境作用下不可避免会产生裂缝1。裂缝会显著影响混凝土坝结构耐久性，严重者会引发渗漏、溶蚀，危害坝体安全，因此，采用混凝土坝裂缝原型监测资料建立具有强解释性和高预测精度的模型，对准确判断裂缝发展及演化趋势、掌握大坝运行状况、及时进行安全示警具有重要意义2。统计模型简单明确，可直观反映各因素对裂缝开度变化的影响，在大坝安全监控领域得到广泛应用。受限于当前监测技术和理论2-4，统计模型在考虑影响因素及裂缝行为的非线性关系上仍较为困难，存在水位温度等多重共线性等问题，难以达到高精度要求。已有研究表明，大坝观测数据残差

2、项中存在混沌现象5。包腾飞等6对残差时间序列进行相空间重构，提出统计混合新模型以提高预测精度；邱莉婷等7基于混沌理论，建立统计回归-Legendre组合模型进行预测，降低误差；严春丽8利用多元线性回归分离影响变量，并结合递归最小二乘(recursive least square,RLS）算法对残差序列进行预测。然而大坝观测数据序列复杂多变，传统方法建模时未能充分考虑数据长短期关联性，易忽略数据内部特征，从而影响结果的可靠性。近些年来人工智能技术快速发展，不依赖于外部特征的深度学习引起了广泛关注。其学习过程不依赖于先验知识，通过特有的学习机制从大量数据中挖掘有效信息，在众多领域如金融、交通流量等

3、领域9，10取得了显著效果11-14。作为一种递归神经网络，门控循环单元（gated recurrent unit，GRU）神经网络可充分利用历史记忆，通过门结构控制与传播信息，在时间序列处理预测中表现优异。裂缝反映了材料内部损伤的不断累积，经统计回归的残差数据本质上是一个伴随复杂变化、噪音以及不确定性的时间序列，包含系统特征及运行规律，应用GRU神经网络挖掘残差序列中长短期序列的潜在联系，可有效弥补传统统计模型回归精度低、忽略残差项等缺陷。此外，为有效消除残差序列中的噪音干扰，引入奇异谱分析（singular spectrum analysis，SSA）对噪声进行分离，提取残差序列的主要成分

4、。本文提出一种MLR-SSA-GRU（multiple linear regression-SSA-GRU）组合模型应用于混凝土坝裂缝开度预测研究，该模型采用统计模型作常规预测，分离水压、温度以及时效影响，分析影响大坝裂缝开度的主要因素。对同样包含大量有效信息的残差项（如各变量非线性影响等），通过SSA分析进行序列重构，并采用GRU神经网络进行训练拟合及预测。通过多种模型融合可以发挥统计模型及深度学习的优势，以达到兼顾强解释性及高精度要求。1 残差重构及预测原理1.1奇异谱分析原理混凝土坝裂缝开度残差序列包含常规统计模型无法解释的成分，影响因素复杂，且具有非平稳性、非线性噪声特点。为分离噪声空

5、间提取主要有效信息，本文对残差序列进行分解重构。SSA是一种高效的时间序列分析方法，可在不依赖于先验知识前提下，基于时间序列本身建立轨迹矩阵，分析结果稳定可靠。具体步骤如下：1）嵌入。对残差序列D=(y1,y2,yN)中心化，其中yN为裂缝开度残差，N为序列长度。选择窗口长度L对序列进行滞后排列得到轨迹矩阵Y：Y=y1y2yLy2y3yL+1yi+1yi+2yi+LyNL+1yNL+2yN(1)2）分解。计算协方差矩阵C=YYT，对其特征值进行分解，得到特征值k及特征向量Vk，其对应的时间主成分aki为aki=YiVk=j=1Lyi+jVk,j,0iNL(2)3）重构。根据时间主成分重构时间序

6、列，得yki=1ij=1iakijVk,j,1Lj=1LakijVk,j,1Ni+1j=iN+LLakijVk,j,1iL1LiNL+1NL+2iN(3)原序列为yi=k=1Lyki(4)根据实际情况，通常只考虑提取原序列的主要成分，即根据奇异值大小选取前q(qL)个重构原序列。重构序列包含原始序列的主要信息，过滤掉部分噪声数据干扰，可显著改善预测模型的表现。1.2门控循环单元神经网络裂缝开度残差序列是单变量长时间跨度序列，其数据在时间维度上具有关联性，GRU神经网络特有内置的门结构可突出重要特征、弱化非重要特征，实现对历史监测数据的有效选择积累，预测表现准确高效。GRU模型的标准结构如图1所

7、示，图1中，为激活函数，=1/1+ex。图1GRU模型结构图Fig.1Structure diagram of GRU model由图1可见，GRU单元由重置门rt、更新门zt构成，其中下标t表示当前时刻，重置门控制前一时刻状态的遗忘程度，更新门决定当前输入对整体的重要性，两者均通过激活函数来决定是否启动。GRU结构通过偏置和相关参数实时调整激活函数的输入，控制门的开闭。具体算法如下：rt=(U(r)xt+W(r)ht1+br)(5)zt=(U(z)xt+W(z)ht1+bz)(6)ht=tanh(U(c)xt+W(c)rtht1+bc)(7)ht=ztht1+(1zt)ht(8)式中：U和W

8、为当前时刻输入xt和上一时刻输出ht1对应的权重矩阵，即相对重要程度，通过梯度下降更新；br、bz、bc为重置门、更新门、候选输出状态的偏置量；tanh为激活函数；当前输入与重置门运算和叠加，形成候选状态输出ht，ht与上一时刻输出ht1通过更新门动态控制形成最终输出ht。GRU模型通过特有机制对既往信息及当前输入筛选记忆，可有效挖掘裂缝开度残差长序列的关联性。2 MLR-SSA-GRU组合模型2.1模型融合为直观反映水位温度等因素对裂缝开度的影响，根据坝工理论2，采用统计回归模型进行了研究，t1时刻观测值为(t1)=H(t1)+T(t1)+(t1)+(t1)(9)式中：(t1)为裂缝开度观测

9、值；H(t1)为水压分量；T(t1)为温度分量；(t1)为时效分量；(t1)为残差项。下一时刻t1+T对应的预测模型可表示为(t1+T)=H(t1+T)+T(t1+T)+(t1+T)(10)这种预测模型仅考虑主要因素（水位、温度、时效）对裂缝开度的影响，机理明确，但未能考虑各变量的非线性影响、裂缝的塑性变形以及测值突变等因素。本文对该预测模型未考虑的残差项进行补充，以提高预测精度。考虑到残差序列复杂多变，采用SSA分离噪声空间，可提取有效信息降低时间序列的复杂程度。本文应用MLR-SSA-GRU组合模型进行混凝土坝裂缝开度预测，分为7个步骤，具体流程见图2。图2MLR-SSA-GRU用于混凝土

10、坝裂缝开度预测流程图Fig.2Flowchart of application of MLR-SSA-GRU model in concrete dam crack opening prediction1）采集裂缝开度监测数据，建立统计回归模型分析各变量影响，并将预测值与实测值对比得到残差序列。2）对残差序列中心化处理，确定窗口长度构建轨迹矩阵。3）确定特征值及特征向量，重构残差序列，根据奇异值大小，选取前q个趋势项及周期项作为有效成分，得到处理后的残差序列。4）划分数据，将残差序列划分为训练集、验证集及测试集。5）构建GRU网络结构：输入当前裂缝开度序列xt以及上一序列xt1计算出裂缝开度数

11、据ht1；xt、ht1经过rt、zt，由函数决定是否激活；通过激活函数tanh得到候选输出状态ht；ht与ht1通过更新门控制组合得到最终输出状态ht。6）采用Adam优化算法更新参数，选取均方误差（mean-sqare error，MSE）作为评价GRU预测模型的依据，对训练好的模型进行测试，得到测试集的裂缝开度预测值。7）在统计回归模型的基础上并入残差序列预测值，构建统计回归-奇异谱分析-门控循环单元神经网络组合预测模型。2.2评价指标为衡量预测模型的精度优劣，选取相关系数R2、均方根误差（root mean square error，RMSE）和平均绝对误差（mean absolute

12、error，MAE）3个指标进行评价：R2=1i=1n(yiyi)2i=1n(yiy)2(11)RMSE=1ni=1n|yiyi|2(12)MAE=1ni=1n|yiyi|(13)式中：yi为实测值；yi为预测值；y为均值。3 工程实例本文以某混凝土重力拱坝为例进行分析，该坝位于皖南长江支流青弋江上游，坝顶高程126.3 m，坝顶宽8 m，自左到右分为28个坝段，对26号坝段105 m高程处测缝选取2022-01-242022-12-31采集到的裂缝开度监测序列，其中前900组用于训练模型，后164组为测试集。3.1统计回归模型的建立混凝土坝裂缝变形主要受温度、水位和时效等因素的影响，结合该坝

13、实际运行状态，建立统计模型：=a0+i=14ai(HiHi0)+i=12(b1i(sin2it2365sin2it0365)+b2i(cos2it2365cos2it0365)+c1(0)+c2(lnln0)(14)式中：a0为常数项；ai、b1i、b2i、c1、c2为回归系数；H0为建模资料系列初始日的水深；H为坝前水深；t2为初始日至观测日天数；t0为始测日到建模资料初始日累计天数；为初始日至观测日天数除以100；0为始测日到建模资料初始日累计天数。根据式（14）对裂缝监测数据训练集进行线性回归分析，裂缝开度拟合效果及各分量影响如图3所示，相应的残差序列如图4所示。图3裂缝开度拟合效果及各

14、分量影响Fig.3Fitting result of crack opening and component influence图4裂缝开度残差序列Fig.4Residual sequence of crack opening由图3可见，温度分量变化趋势与裂缝开度基本一致，约占裂缝开度年变幅的70%，为影响大坝裂缝开度变化的主要因素；水压分量约占裂缝开度年变幅的25%，为次要因素；时效分量约占裂缝开度年变幅的5%，对裂缝开度影响较小。由此可见，温度变化是裂缝开度变化的主要影响因素，是运行监测的重点关注项目，在大坝日常运行管理时需引起重视。3.2残差序列重构及预测混凝土坝裂缝开度残差序列性态复杂

15、，可采用SSA分解提取序列中有效信息。在SSA中，窗口长度的选取会影响周期识别及噪声分离效果15。本文选取不同的窗口长度分别进行残差序列重构，根据重构效果确定最佳取值为30 d，分解结果如图5、6所示。图5残差序列SSA分解周期和趋势分量Fig.5SSA decomposition periodic and trend terms of residual sequence图6SSA分解各分量贡献Fig.6Contribution of each component decomposed by SSA在SSA分解得到的若干分量中，根据奇异值大小选取前5个分量（PC1PC5）作为有效信息，其累计贡

16、献率已达92.36%，视为已涵盖序列主要成分。其中，前2个分量（PC1、PC2）为低频信号，变化幅度较小，主要受到水位温度等周期性荷载作用；后3个分量（PC3PC5）为高频信号，受到不确定荷载及噪声影响。重构后的残差序列如图7所示。图7残差序列及重构值Fig.7Residual sequence and reconstructed value对图7所示的重构序列应用GRU神经网络进行训练拟合，GRU神经元层数过少会导致模型解释能力不足，过多则容易过拟合，层数一般为13层。结合现有研究16，取GRU的层数为2层，采用贝叶斯优化确定模型具体参数，如表1所示，对应MSE为0.000 14，此时训练集

17、预测相关系数R2为0.999 9。表1GRU参数范围及优化值Table 1Parameters range and optimization values of GRU神经元层数参数范围寻优值随机失活层(0.1,0.3)0.25第1层(100,300)200第2层(100,300)247训练好的MLR-SSA-GRU组合模型应用于测试集进行裂缝开度预测，结果如图8所示。图8MLR-SSA-GRU模型预测结果Fig.8The prediction result of MLR-SSA-GRU model由图8可见，本文的MLR-SSA-GRU模型预测值与实测序列吻合良好，且预测值大都落在95%置信

18、区间，表明其预测性能稳定，鲁棒性强。为进一步评估本文所提MLR-SSA-GRU模型应用于混凝土坝裂缝开度预测的性能优劣，引入多层神经网络（multi layer perception,MLP）、随机森林（random forest，RF）以及支持向量回归机（support vector regression，SVR）进行对比分析，评估指标选用R2、RMSE和MAE，如表2所示。表2各模型预测性能比较Table 2Comparison of prediction performance of different models模型R2RMSE/mmMAE/mmMLRSSAGRU0.9970.0120.008MLP07810.0970.077RF0.8490.0810.067SVR0.9380.0520.043除预测精度评估外，引入Diebold-Mariano（DM）假设检验用于评估2个不同预测模型