文章快速检索     高级检索
  大地测量与地球动力学  2019, Vol. 39 Issue (10): 1081-1085  DOI: 10.14075/j.jgg.2019.10.018

引用本文  

张东华, 李志娟, 刘全明, 等. 大坝变形的奇异谱分析预测[J]. 大地测量与地球动力学, 2019, 39(10): 1081-1085.
ZHANG Donghua, LI Zhijuan, LIU Quanming, et al. Singular Spectrum Analysis for Analyzing and Forecasting the Dam Deformation[J]. Journal of Geodesy and Geodynamics, 2019, 39(10): 1081-1085.

项目来源

国家自然科学基金(51569018, 51509131)。

Foundation support

National Natural Science Foundation of China, No.51569018, 51509131.

通讯作者

刘全明,博士,副教授,主要研究方向为大地测量数据处理与遥感反演,E-mail:61654346@qq.com

Corresponding author

LIU Quanming, PhD, associate professor, majors in geodesy data processing and remote sensing inversion, E-mail:61654346@qq.com.

第一作者简介

张东华,讲师,主要研究方向为大地测量数据处理、GIS算法与开发,E-mail:zdh@imau.edu.cn

About the first author

ZHANG Donghua, lecturer, majors in geodesy data processing and GIS algorithm and program development, E-mail:zdh@imau.edu.cn.

文章历史

收稿日期:2018-10-08
大坝变形的奇异谱分析预测
张东华1     李志娟2     刘全明1     罗艳云1     
1. 内蒙古农业大学水利与土木建筑工程学院,呼和浩特市昭乌达路306号,010018;
2. 内蒙古自治区航空遥感测绘院,呼和浩特市兴安南路42号,010010
摘要:利用奇异谱分析法对大坝变形数据进行分析,提取趋势和周期分量,分析影响因子与各变形分量的相关性。结果表明,大坝变形主要与水位变化和时效因子有关,温度变化对大坝变形的周期成分贡献最大,其次为水位。另外在准确提取信号的基础上, 利用奇异谱分析迭代法对大坝变形进行预测,并与多元回归分析方法和高斯过程模型进行对比,发现其预测精度明显高于后两者。
关键词奇异谱分析大坝变形变形预测

大坝变形受温度、水位和时间效应等多因素影响,很难直接建立较为准确的变形量与影响因素之间的非线性关系模型。目前常用的预测方法主要有多元线性回归、时间序列模型、灰色理论和神经网络等[1-2],但大坝变形具有非平稳性和非线性,传统方法的应用受到一定限制。奇异谱分析(singular spectrum analysis,SSA)从时间序列向空间动力结构出发,以经验正交的方式展开时间序列子空间,不基于正弦波假定前提,不需要先验信息,能够准确识别时间序列中的周期信号[3]。在使用GPS进行大坝变形监测时,根据SSA理论准确提取大坝变形信号,可提高变形监测的精度,准确获取大坝的真实变形信息,进一步分析各影响因子对大坝变形的贡献情况,深入研究大坝变形的内在规律并进行预测[4-5]。考虑到SSA提取信号的优势,在准确分析序列提取信号的基础上,基于SSA迭代法对大坝变形进行预测,并与高斯过程模型、多元回归分析方法进行对比,期望能够更好地预测大坝的变形情况。

1 SSA预测方法 1.1 SSA原理

SSA是根据所观测到的时间序列构造出轨迹矩阵,对矩阵进行分解、重构,提取代表原时间序列不同成分的信号,通过分析时间序列的结构进行预测[3]。首先,对于中心化后的时间序列xi(1≤iN),根据窗口L按照式(1)构造L×(NL+1)的轨迹矩阵X

$ \mathit{\boldsymbol{X}} = \left[ {\begin{array}{*{20}{c}} {{x_1}}&{{x_2}}& \cdots &{{x_{i + 1}}}& \cdots &{{x_{N - L + 1}}}\\ {{x_2}}&{{x_3}}& \cdots &{{x_{i + 2}}}& \cdots &{{x_{N - L + 2}}}\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ {{x_L}}&{{x_{L + 1}}}& \cdots &{{x_{i + 1}}}& \cdots &{{x_N}} \end{array}} \right] $ (1)

可得滞后协方差阵为C=XXT,对其进行特征值分解:

$ \mathit{\boldsymbol{C}} = \mathit{\boldsymbol{V \boldsymbol{\varLambda} }}{\mathit{\boldsymbol{V}}^{\rm{T}}} $ (2)

式中,Λ为特征值λk构成的对角阵,λkC的特征值并按降序排列,vk为正交矩阵V的行向量。主成分矩阵A计算公式如下:

$ \mathit{\boldsymbol{A}}{\rm{ = }}\mathit{\boldsymbol{VX}} $ (3)

A的第k行(ak)即为第k个主成分,其中每个数值ak, i为:

$ {a_{k, i}} = \sum\limits_{j = 1}^L {{x_{i + j - 1}}} {v_{k, j}}, 1 \le i \le N - L + 1 $ (4)

其中, vk, jvk的第j个元素。时间序列第k个重构分量可以写为:

$ x_{i}^{k}=\\\left\{ {\begin{array}{*{20}{l}} {\frac{1}{i}\sum\limits_{j = 1}^i {{a_{k, i - j + 1}}} {v_{k, j}}, 1 \le i \le L - 1}\\ {\frac{1}{L}\sum\limits_{j = 1}^L {{a_{k \cdot - j + 1}}} {v_{k, j}}, L \le i \le N - L + 1}\\ {\frac{1}{{N - i + 1}}\sum\limits_{j - i - N + L}^L {{a_{k, i - j + 1}}} {v_{k, j}}, \quad N - L + 2 \le i \le N} \end{array}} \right. $ (5)

SSA的特征值是按降序排列的,因此时间序列可由前几个主成分分量重构,其余则视为噪声,即信号${\hat x}$可以由前q个主成分得到:

$ \hat x = \sum\limits_{k = 1}^q {x_i^k} , i = 1, 2, \cdots , N $ (6)
1.2 SSA预测原理

SSA迭代预测法建立在重构序列基础上,具体步骤如下。

定义矩阵${\mathit{\boldsymbol{X}}_\mathit{M}} = [{{\hat X}_1}, \cdots , {{\hat X}_K}, {Z_{K + 1}}, \cdots {Z_{K + M}}]$,其中,$ [{{\hat X}_1}, \cdots , {{\hat X}_K}]$为基于原有数据序列构建的重构轨迹矩阵,[ZK+1, …, ZK+M]设为未来预测序列构建的相应矩阵向量。因此,只要获得XM便可对其进行对角平均,得到预测序列。定义R=(aL-1, …, a1)T,设I为选取的重构特征值,PiRL(iI)为对应的特征向量, PiPi的前(L-1)行,πiPi的最后一个元素,${v^2} = \sum\nolimits_{i \in I} {\mathit{\pi }_i^2} $

$ \boldsymbol{\mathit{R}}=\frac{1}{1-{{v}^{2}}}\sum\limits_{i\in I}{{{\pi }_{i}}}\underline{{{\boldsymbol{\mathit{P}}}_{i}}} $ (7)

定义线性算子${p_{{\rm{vec}}}}\mathit{\boldsymbol{Y}} = \left( {\begin{array}{*{20}{l}} {\mathit{\boldsymbol{\overline Y}} }\\ {{\mathit{\boldsymbol{R}}^{\rm{T}}}\mathit{\boldsymbol{\overline Y}} } \end{array}} \right)$${\mathit{\boldsymbol{\overline Y}} }$Y的后(L-1)行,写成向量形式即为:

$ {{\boldsymbol{\mathit{Z}}}_{i}}=\left\{ \begin{array}{*{35}{l}} {{{\hat{X}}}_{i}}, i=1, \cdots , K \\ {{p}_{vex}}{{\boldsymbol{\mathit{Z}}}_{i-1}}, i=K+1, \cdots , K+M \\ \end{array} \right. $ (8)

对式(8)进行对角平均可以计算出预测序列,式(8)的等价形式为:

$ {{y}_{i}}=\left\{ \begin{array}{*{35}{l}} {{{\tilde{x}}}_{i}}, i=1, \cdots , N \\ \sum\limits_{j=1}^{L-1}{{{a}_{j}}}{{y}_{i-j}}, i=N+1, \cdots , N+M \\ \end{array} \right. $ (9)

式中,${{\tilde x}_i}$为重构序列的第i个元素。本文利用蒙特卡洛模拟算法[6]估计预测序列的95%可信区间,具体预测原理见文献[3]。

2 大坝变形成因分析

在西南某蓄能水电厂的下库混凝土大坝上布设正倒垂线及4个位移监测点,坝体内安装18处电阻式温度计,于2006-11~2010-08对温度、水位、时间和位移变化量进行1 400次记录。由于各监测点变化规律大致相同,本文选择其中一个点进行分析。对采集数据作粗差探测预处理,并通过平滑处理(平滑半径为7 d)得到200期数据集。

图 1可知,大坝位移变形随时间、温度、水位的变化具有一定的规律性,分析变形量与影响因素的关联性,可以为变形预测打下良好基础。对原始序列利用SSA进行分析,根据先验信息得知大坝变形存在周年、半年和季节性周期,因此窗口长度取1 a(即52周),并给出前30个主成分的W-correlations图(图 2),具体公式见文献[5]。由图 2可知,PC1是趋势项,能量占比最大,为57.29%;PC2-3为年周期分量,能量占比为41.79%;PC4-5为半年周期分量,能量占比为0.39%;PC7-8为季节性周期分量,能量占比为0.13%。前8个主成分能量占比超过总能量的99%,剩余的分量被视为噪声予以忽略,具体的各成分见图 3

图 1 大坝变形数据序列 Fig. 1 The data series of dam deformation

图 2 W-correlations中前30个分量图 Fig. 2 The first 30 PCs of W-correlations

图 3 大坝变形奇异谱分析周期趋势项 Fig. 3 The trend and periodic terms of dam deformation

为分析水位、温度和时效因子对大坝变形的贡献情况,将提取的趋势、周期和余量与影响因子之间进行相关性分析,并利用奇异谱分析提取水位变化中相应的趋势和周期成分(表 1),同时给出变形量、水位和温度变化的周年和半年分量对比情况(图 4图 5)。由表 1图 3可知,大坝变形的每个主成分分量都包含原序列中的不同特征信息,PC1分量变化平缓,主要包含长期变形信息,与水位因子趋势成分相关系数为0.86,与时效因子相关系数为-0.88,即大坝变形长期趋势成分主要与水位的长期趋势变化和时效因子有关,受温度等其他因素的影响较小。大坝周期分量与水位的相应周期分量相关系数为0.38,为正相关;与温度因子相关系数为-0.92,为负相关,即与温度的变化趋势恰好相反,这与混凝土坝体由温度变化引起的坝体变形规律相符合,也与图 4图 5显示的结果一致。由图 4还可以发现,大坝变形与水位及温度的变化存在一定的相位差,这与大坝变形具有滞后效应的规律一致。余量数据关联性较弱,与水位、温度和时效因子的相关性不强,能量占比较小,此处不再考虑。

表 1 各分量与主要影响因素相关性分析 Tab. 1 The results of correlation analysis with the main factors

图 4 变形量、水位和温度变化周年分量对比 Fig. 4 Comparison of annual components of deformation, water level and temperature

图 5 变形量、水位和温度变化半年分量对比 Fig. 5 Comparison of semi-year components of deformation, water level and temperature
3 大坝变形预测

将前190期变形数据{Xi=(ti, Hi, Ti), Yi(i=1, 2, …, 190)}作为建模数据集,191~120期变形数据{Xj=(tj, Hj, Tj), Yj}作为测试数据集。安全监测因果模型的一项重要工作就是将合理有效的影响因子引入模型,以达到准确预测的目的。运行期间大坝变形的主要影响因素是水位、温度和时效。本文使用多元回归分析方法、高斯过程模型(其中高斯过程模型选择的神经网络核函数及均值函数采用MeanZero)和SSA迭代法[7]进行预测建模比较分析,并利用蒙特卡洛模拟算法[6]估计SSA预测序列的95%可信区间,具体预测结果见表 2

表 2 各预测模型预测结果对比 Tab. 2 Comparison of the predicted values for three models

表 2图 6可知,多元回归分析方法预测效果较差,表明针对多因素主导的大坝变形监测,其预测模型未能较好地建立大坝变形的因果关系;高斯过程模型在处理高维数、小样本、非线性等复杂问题上具有很好的适应性,预测效果明显更接近于实测值[8]。与前两种方法相比,SSA能够准确分析并提取大坝变形的趋势和周期信号,并利用迭代法进行预测,预测精度明显提升。

图 6 预测与残差绝对值对比 Fig. 6 The comparison of predictions and absolute errors

为了进一步对比不同预测模型的预测精度,采用均方根误差和平均绝对误差[9]检验预测模型的精度:

$ \text{RMSE}=\sqrt{\frac{1}{N}\sum\limits_{i=1}^{N}{{{\left( {{s}_{i}}-{{{\hat{s}}}_{i}} \right)}^{2}}}} $ (10)
$ \text{MAE}=\frac{1}{N}\sum\limits_{i=1}^{N}{\left| {{s}_{i}}-{{{\hat{s}}}_{i}} \right|} $ (11)

式中,${{\hat s}_i}$为模型预测值,si为大坝变形实测值,N为预测序列长度。

经计算,SSA迭代预测法的RMSE和MAE分别为0.45 mm和0.34 mm,明显优于高斯过程的1.28 mm和0.86 mm及多元回归分析模型的1.75 mm和1.63 mm。由此可知,基于SSA的迭代预测法结果明显小于传统的多元回归分析及高斯过程预测方法,与表 2图 6的结果一致。文中利用蒙特卡洛模拟算法计算了预测序列的可信区间,真实数据基本上处于可信区间内,表明SSA迭代预测法可以有效地预测大坝变形,为大坝安全监测提供保障。

4 结语

利用SSA法分析大坝变形数据并建立大坝变形预测模型,通过实例验证分析可得:1)SSA可以有效分离并提取序列中的趋势和周期成分,发现大坝变形存在年周期、半年周期和季节性周期,并分别分析与相应变形因子之间的相关性。大坝变形的趋势项主要与水位趋势变化和时效因子有关,温度变化对大坝变形的周期成分贡献最大,其次为水位变化。2)考虑到SSA提取信号的优势,在SSA重构序列基础上,利用迭代法对大坝变形进行预测,并与高斯过程和多元回归分析进行对比。结果表明,基于SSA的迭代预测法预测精度均方根误差和平均绝对误差分别为0.45 mm和0.34 mm,明显高于高斯过程模型和多元回归分析方法,且给出了95%可信区间,适合对大坝变形进行短中期预测。

如何根据大坝实时监测数据对其进行动态预测,进一步利用新信息建模,是接下来研究的重点。

参考文献
[1]
周世健, 罗亦泳, 鲁铁定. 局部均值分解(LMD)方法在多尺度变形分析中的应用探讨[J]. 江西科学, 2014, 32(4): 461-466 (Zhou Shijian, Luo Yiyong, Lu Tieding. Discussion for the Local Mend Decomposition(LMD) in Analysis of Multi-Scale Deformation[J]. Jiangxi Science, 2014, 32(4): 461-466) (0)
[2]
王奉伟, 周世健, 周清, 等. 局部均值分解与支持向量回归的大坝变形预测[J]. 测绘科学, 2016, 41(10): 132-135 (Wang Fengwei, Zhou Shijian, Zhou Qing, et al. Forecasting Model of Dam Deformation Based on LMD-SVR Method[J]. Science of Surveying and Mapping, 2016, 41(10): 132-135) (0)
[3]
Golyandina N, Zhigljavsky A. Singular Spectrum Analysis for Time Series[M]. Berlin: Springer, 2013 (0)
[4]
Gourine B, Khelifa S. Analysis of Dam Deformation Using Artificial Neural Networks Methods and Singular Spectrum Analysis[C].Euro-Mediterranean Conference for Environmental Integration, Cham, 2017 (0)
[5]
Khan M A R, Poskitt D S. Forecasting Stochastic Processes Using Singular Spectrum Analysis: Aspects of the Theory and Application[J]. International Journal of Forecasting, 2017, 33(1): 199-213 DOI:10.1016/j.ijforecast.2016.01.003 (0)
[6]
刘岳巍.自适应蒙特卡洛方法和固定宽度置信区间[D].兰州: 兰州大学, 2013 (Liu Yuewei. Adaptive Monte Carlo Method and Fixed Width Confidence Interval[D].Lanzhou: Lanzhou University, 2013) http://cdmd.cnki.com.cn/Article/CDMD-10730-1014135075.htm (0)
[7]
Golyandina N. On the Choice of Parameters in Singular Spectrum Analysis and Related Subspace-Based Methods[J]. Statistics and Its Interface, 2010, 1(5): 403-413 (0)
[8]
徐冲, 刘保国, 刘开云, 等. 基于粒子群-高斯过程回归耦合算法的滑坡位移时序分析预测智能模型[J]. 岩土力学, 2011, 32(6): 1 669-1 675 (Xu Chong, Liu Baoguo, Liu Kaiyun, et al. Intelligent Analysis Model of Landslide Displacement Time Series Based on Coupling PSO-GPR[J]. Rock and Soil Mechanics, 2011, 32(6): 1 669-1 675) (0)
[9]
王奉伟, 周昀琦, 周世健, 等. 大坝变形多尺度分析ELMD-LSSVM预测模型[J]. 辽宁工程技术大学学报:自然科学版, 2016, 35(12): 1 475-1 479 (Wang Fengwei, Zhou Yunqi, Zhou Shijian, et al. Forecasting Model of Dam Deformation Based on ELMD-LSSVM Method[J]. Journal of Liaoning Technical University:Natural Science, 2016, 35(12): 1 475-1 479) (0)
Singular Spectrum Analysis for Analyzing and Forecasting the Dam Deformation
ZHANG Donghua1     LI Zhijuan2     LIU Quanming1     LUO Yanyun1     
1. Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, 306 Zhaowuda Road, Hohhot 010018, China;
2. Institute of Aerial Remote Sensing Surveying and Mapping of Inner Mongolia Autonomous Region, 42 South-Xing'an Road, Hohhot 010010, China
Abstract: We use singular spectrum analysis to extract the corresponding components and compute the correlation coefficients, under the influence of factors. The results show that the trend term of dam deformation mainly relates to water level and aging factor. For seasonal terms, the temperature factor contributes more than water level. Experimental results of dam deformation show that SSA can extract the trends and periodic signal effectively and is useful to forecast dam deformation. Then, recurrent forecasting method of SSA is used for the prediction of dam deformation. Compared with Gaussian process and multiple-regression analysis, the results show that SSA is an effective method with a higher predictive accuracy for analyzing and forecasting dam deformation.
Key words: singular spectrum analysis(SSA); dam deformation; deformation prediction