桥梁水平位移混沌特征识别与神经网络预测研究

2021-01-07栾元重梁耀东翁丽媛刘承旭

大地测量与地球动力学 2021年1期

栾元重梁耀东董岳翁丽媛刘承旭

1 山东科技大学测绘科学与工程学院，青岛市前湾港路579号，266590

跨海大桥是人工与自然有机结合的整体，处在海洋这个复杂的环境中，运营后还需承受车载外力的作用，其各部分之间不断地进行能量、物质交换，通过系统信息熵的变化使桥梁整体走向有序[1-3]。因此，桥梁变形是一种非平衡态有序结构。由于我国跨海大桥建设也是近10 a的事情，运营后的跨海大桥变形特征及预测模型研究也处于起步阶段[4-5]。近几年建好的几座跨海大桥均建立了变形监测站，获得了桥梁变形数据，但在数据挖掘、建立动态非线性高精度的桥梁变形时空预测模型方面还远远不够[6-7]。因此，开展跨海大桥时空非线性变形模型研究是理论与实际结合、理论成果又反馈桥梁安全运营的实用课题[8-9]。

本文对青岛胶州湾跨海大桥开展了变形监测工作，桥梁特征点三维坐标数据采集使用的测量机器人为徕卡30全站仪，该仪器测角精度为0.5″,测距精度为0.1 mm+10-6D。测量方法为：在稳定区域布设的D级GPS控制点上安置测量机器人，对桥墩、桥面、桥索塔的特征点自动进行三维坐标观测。由于桥梁处于海洋环境中，受海浪、飓风、船体撞击等因素的影响，使桥梁变形表现为动态非连续变形，并伴有混沌现象发生[10]。

目前用于混沌时间序列预测的数学模型主要有神经网络模型和支持向量机等，但由于支持向量机存在参数准确选取困难[11]，并且局部优化能力差的缺点，因此本文根据桥梁变形监测数据，首先对桥梁变形系统混沌特征进行识别，并对桥梁水平位移单变量进行高精度RBF神经网络混沌预测的深入研究，这对揭示跨海大桥变形机理、实现桥梁安全运营具有重大意义。

1 桥梁水平位移观测数据混沌特征识别

当桥梁系统受外界影响有扰动时，系统可能发生分岔并可能演变成混沌状态。为了判断桥梁变形是否有混沌状态，选取Verhulst生物模型拟合桥梁变形位移时间序列[1]，即

(1)

式中，x为变形位移值，t为时间，k为参数。其差分方程为：

xn+1=k(1-xn)xn

(2)

本文选取青岛胶州湾跨海大桥W1点3 a内的300期水平位移监测数据进行分析，取W1点后100期数据水平位移观测，其在x方向、y方向上的变形量无明显线性与周期性变化趋势，通过重构相空间寻找其内在规律。

1.1 最大Lyapunov指数

Lyapunov指数表示多次迭代中平均每次迭代引起指数分离的度量。Lyapunov指数λ的表达式为：

(3)

式中，λ为Lyapunov指数，F(x)为x的映射值。λ>0时，运动系统存在混沌状态。

基于最大 Lyapunov 指数的混沌时间序列预测，就是利用已经得到的最大Lyapunov指数(即轨道平均发散率)，采用相空间重构技术，通过跟踪最邻近点的演化进行预测[12-14]。格里波基[15]证明，只要最大 Lyapunov 指数大于零，就可以肯定系统中存在混沌状态。

用小数据量法计算桥梁水平位移时间序列最大Lyapunov指数，x方向变形时间序列最大Lyapunov指数为0.015，证明x方向变形时间序列具有混沌性；同理，y方向变形时间序列最大Lyapunov指数为0.014，y方向变形时间序列也具有混沌性，此结果与Verhulst生物模型拟合的桥梁变形位移时间序列具有混沌状态相吻合。

1.2 时间延迟

C-C法求取时间延迟的方法如下。将时间序列分成t个不相交的子序列，每个子序列可定义为：

S(m,N,r,t)=

(4)

式中，m为嵌入维数，N为点的个数，r为元素间相关半径，t为不相交子序列个数。

令N→，有：

(5)

由于时间序列是有限的，且元素间也可能是相关的，故选择对应值中最大和最小的两个半径r，其差量为：

ΔS(m,t)=max{S(m,rj,t)}-{S(m,rj,t)}

(6)

(7)

(8)

(9)

图1 C-C法计算x方向变形时间序列的时间延迟Fig.1 Time delay of the deformation time series of x direction calculated by C-C method

1.3 嵌入维数m

若{x1,x2,…,xi,…}为实测的一组等间隔变形数据，取{x1,x2,…,xi,…}中部分数据构成一个m维子空间，空间中N个点表达为：

(10)

令

(11)

式中，θ(x)指Heaviside函数，N为点的个数。

吸引子的关联维数计算公式为：

(12)

分别求出不同嵌入维数m所对应的dC值。取dC值第一个到达平稳时的m估值作为嵌入维数，此时dC值可作为吸引子维数的估值，则m=2[d]+1。这就是嵌入维数估计的G-P算法。

x方向时间延迟τ=1，设定维数m=2，3，…，12，分别绘制x方向变形时间序列饱和关联维数法(简称G-P法)的lnC(r)-lnr曲线(图2)，图中曲线自上而下对应于m=2，3，…，12。当维数增加到7时，lnC(r)与lnr的比值不再随m的增大发生变化，因此x方向嵌入维数为7。同理，求得y方向变形时间嵌入维数为7。

图2 x方向变形时间序列lnC(r)-lnr图Fig.2 lnC(r)-lnr diagram of deformation time series of x direction

1.4 二维重构轨迹分析

重构相空间的目的是保持原先变形系统的几何结构，并与其具有等价的动力系统。相空间重构最重要的参数为时间延迟τ和嵌入维数m。延迟滞时重构的变形动力系统相空间简化模型[17-18]为：

X(t)=(x(t),x(t-τ),…,x(t-(m-1)τ))

(13)

式中，m为嵌入维数，τ为时间延迟，t=1,2,…,n。

x方向和y方向最大Lyapunov指数之和为二维最大Lyapunov指数，即λ=λx+λy=0.029，故二维数据同样具有混沌特性。按时间延迟τ=1、m=7等参数对W1点水平位移时间序列进行二维相空间重构(图3)。

图3 水平位移二维相空间重构Fig.3 Two-dimensional phase space reconstruction of horizontal displacement

2 桥梁水平位移混沌时间序列预测

2.1 基于RBF神经网络的混沌预测

X[i]=x(t+(i-1)τ)

(14)

(15)

通过相空间重构得到输入层单元数为m、隐层单元数为k、输出单元数为1的网络拓扑结构m×k×1。

取W1点300期水平位移观测数据，前250期数据用来训练，后50期数据用来预测对比分析。对W1点x方向变形数据进行归一化处理，并将前250组数据作为训练样本，RBF神经网络误差曲线如图4所示。由图可见，在网络训练过程中，只需6步就到达了目标标准差值0.001。

图4 RBF神经网络误差曲线Fig.4 RBF neural network error curve

对后输入的50个数据进行相空间重构，组成一个N-(m-1)τ样本输入，其中τ=1，m=7，输入样本点为44个，第t个输入样本为：

X(t)={x(t),x(t+1),x(t+2),

x(t+3),x(t+4),x(t+5),x(t+6)}

(16)

将44个输入样本输入RBF神经网络预测模型，计算得到RBF神经网络混沌预测值。

2.2 RBF神经网络混沌预测与最大Lyapunov指数混沌时间序列预测的比较

将RBF神经网络混沌预测与基于最大Lyapunov指数的混沌时间序列预测的误差进行比较，分析两种预测模型的精度。

对于x方向的变形，对后50期数据进行预测，第251～259期数据预测结果如表1所示，实测值与基于RBF神经网络的预测值和基于最大Lyapunov指数的混沌时间序列预测值对比如图5(a)所示。

由表1绝对差和RMSE对比分析可知，基于RBF神经网络混沌预测模型的预测值比基于最大Lyapunov指数混沌时间序列预测模型的预测值要更接近于实测值且精度明显提高。由图5(b)可见，随着观测期数的增加，预测值逐渐偏离实测值，残差呈现逐渐变大的趋势。这证明本文基于RBF神经网络所建立的混沌预测模型在短期预测的情况下，比基于最大Lyapunov指数的混沌时间序列预测模型得到的预测结果更接近于真实值。这从一个方面证明混沌预测对时间序列只能进行短期预测。