-
变形分析经过多年的研究和发展,取得丰硕的成果。由于引起变形的原因十分复杂,一般来说与变形体自身的物理特性、力学性质、地质条件和外界环境等诸多因素有关,使得变形具有不确定性和错综复杂性,从而诞生了多种变形分析的理论和方法。
到目前为此,回归统计分析法仍然是变形分析的主要方法,应用多元回归分析法建立变形与变形因子之间函数关系然后进行定量预报,常用于大坝变形预测[1-2]、边坡变形分析和预测[3-5]等。回归分析在确定影响变形的因子时有时只能推断,并且有的因子不可测,使得回归分析在有限的条件下受到限制[6]。
由于引起变形的原因复杂,许多工程实践中的时间序列数据经常表现出较明显的非线性特征。例如边坡的变形进入加速期,其变形值时间序列多数情况下会表现出较显著的非线性特点,另外确定影响变形的因素也是一个难点。
随着机器学习等新兴学科的发展,人工智能和机器学习理论逐渐引起了测绘科学界的重视。目前,应用于监测数据分析的人工智能机器学习主要代表方法是BP神经网络[7-10]和支持向量机(support vector machine,SVM)[11-13]。这些方法具有分布存储、并行处理、自学习、自组织以及非线性映射等优点,由这些方法与其他技术的结合形成的混合方法和混合系统也有许多研究应用案例[14-15]。将智能机器学习方法引入到测绘科学领域进行学科交叉研究,能够较好地解决变形监测中变形机理模糊、处理非线性系统等复杂问题。
近几年,高斯过程回归(Gaussian processes regression,GPR)方法在国际上成为智能机器学习领域的研究热点[16]。但GPR应用于变形数据分析和预测的研究相对较少。国内已有学者开始应用GPR对地下工程岩体非线性行为进行预测[17-18]、边坡可靠度分析[19]和边坡变形预测[20],这些研究成果主要是利用GPR已有的核函数开展研究,但关于GPR用于变形分析时核函数的选择、全局最优超参数快速求解、超参数动态更新及训练样本集优化等问题还有待进一步研究。
本文主要利用GPR非线性、自适应等优点对变形监测时空序列数据进行分析,研究基于GPR变形智能预测模型。
-
回归分析方法就是用来确定两种或两种以上变量,即因变量y和自变量x之间相互依赖的定量关系的一种统计分析方法,回归模型中的x和y的关系由任意假设的函数f(x)表示,回归问题的一般模型表示为:
$$ \boldsymbol{y} = f\left( x \right) + \boldsymbol{\varepsilon} ,\boldsymbol{\varepsilon} \sim N(0,\boldsymbol{\sigma} ^2_n{\boldsymbol{I}_n}) $$ (1) 其中,y=[y1, y2, …, yn]T是受噪声污染的n×1维的观测矢量,ε是相互独立服从高斯分布的观测噪声向量,σn2为噪声方差,In是单位阵;x=[xi, 1, …, xi, d]T(i=1, 2, …, n),x∈Rd是n×1维的服从高斯分布的随机变量。
GPR经过对观测数据的学习,能间接地精确地描绘出f(x),GPR是基于函数的,由均值函数m(x)和协方差函数k(x, x)唯一确定,因此,GPR通常定义为:
$$ f\left( x \right) \sim GP\left( {m\left( x \right),k\left( {x,x} \right)} \right) $$ (2) 其中,
$$ \left\{ \begin{array}{l} \boldsymbol{m}\left( x \right) = \boldsymbol{E}\left[ {f\left( x \right)} \right]\\ \boldsymbol{k}\left( {\boldsymbol{x},\boldsymbol{x}} \right) = \boldsymbol{E}[\left( {f\left( x \right) - m\left( x \right)} \right){\left( {f\left( x \right) - m\left( x \right)} \right)^{\rm{T}}}] \end{array} \right. $$ 函数f(x)是由多维高斯分布组成的,根据多维高斯分布的性质,由式(1)和式(2)可得观测值y的先验分布为:
$$ \boldsymbol{y} \sim \boldsymbol{N}(0,\boldsymbol{k}\left( {x,x} \right) + \boldsymbol{\sigma} ^2_0{\boldsymbol{I}_n}) $$ (3) 通常作数据预处理时减去均值,使m(x)=0,由高斯分布的性质进一步得到观测值y与输出样本y*(预测值)的联合先验分布:
$$ \left[ {\begin{array}{*{20}{c}} \boldsymbol{y}\\ {{\boldsymbol{y}_*}} \end{array}} \right] \sim N\left( {0,\left[ {\begin{array}{*{20}{c}} {\boldsymbol{k}\left( {x,x} \right) + \boldsymbol{\sigma} ^2_n{\boldsymbol{I}_n}}&{\boldsymbol{k}(x,{x_*})}\\ {\boldsymbol{k}({x_*},x)}&{\boldsymbol{k}({x_*},{x_*})} \end{array}} \right]} \right) $$ (4) 式中,k(x, x)是输入样本x的协方差矩阵,且为n×n的对称正定矩阵,各个元素k(xi, xj)表示xi与xj的相关性;k(x, x*)=k(x*, x)T是输入样本x与待预测输入值x*之间的n×1的协方差矩阵;k(x*, x*)为待预测输入值x*自身的方差。
由x和y组成观测数据集D={(xi, yi)|i=1, 2, …, n)},D称为训练样本集或学习样本集,在已获得的训练集D条件下,y*的后验分布为:
$$ \boldsymbol{p}({y_*}|D,{x_*}) \sim \boldsymbol{N}(m({y_*}),\boldsymbol{k}({y_*},{y_*})) $$ (5) 其中y*的均值和方差分别为:
$$ \left\{ \begin{array}{l} \boldsymbol{m}({y_*}) = \boldsymbol{k}({x_*},x){\left[ {\boldsymbol{k}\left( {x,x} \right) + \boldsymbol{\sigma} ^2_n{\boldsymbol{I}_n}} \right]^{ - 1}}y\\ \boldsymbol{k}\left( {{y_*},{y_*}} \right) = \boldsymbol{k}\left( {{x_*},{x_*}} \right) - \boldsymbol{k}\left( {{x_*},x} \right){[\boldsymbol{k}\left( {x,x} \right) + \\ \boldsymbol{\sigma} ^2_n{\boldsymbol{I}_n}]^{ - 1}}\boldsymbol{k}{\left( {{x_*},x} \right)^{\rm{T}}} \end{array} \right. $$ (6) 式中, m(y*)即为待预测值x*的对应的输出值y*的均值; k(y*, y*)是输出预测值的后验方差,可以用来度量预测结果的不确定性,即可信程度。这是GPR区别于BP和SVM方法的一个独特优点。
-
由式(6)可知,GPR中的协方差函数(核函数)对GPR的预测性能起着决定性的影响,协方差函数是用来衡量训练样本与预测样本之间相似的程度,其中平方指数协方差函数(squared exponential covariance function,SE)因其无穷可微使得GPR非常平滑,是最为常用协方差函数,其形式为:
$$ \boldsymbol{k}({x_i},{x_j}) = \boldsymbol{\sigma} ^2_f\exp\left( { - \frac{1}{{2{\boldsymbol{l}^2}}}{r^2}} \right) + \boldsymbol{\sigma} ^2_n{\delta _{ij}} $$ (7) 式中,r2=$\sum\limits_{k = 1}^d {{{\left( {{x_{i,k}} - {x_{j,k}}} \right)}^2}} $; σf2为信号方差,用于控制局部相关性的程度; l是特征长度尺度参数,表示输出输入的相关性大小,其值越大相关性越小。σf2和l使得协方差函数具有很好的适用性。σn2为噪声方差,δij为克洛内克尔函数,习惯上将参数集θ={l, σf2, σn2}称为超参数。
通过建立训练样本条件概率的负对数作为似然函数求偏导,再采用牛顿法、共轭梯度法等优化方法进行最小化求得超参数的最优解, 有关其他协方差函数和超参数求解在文献[21]中有详细论述。
应用训练样本集求解超参数这一过程称为训练,当训练完成后,即可求得相应的协方差阵,将待预测值x*代入式(6)得到预测值y*及其预测方差。
-
GPR用于变形分析和预测,不是简单的照搬和套用,需要在解决好以下两个基本问题的基础上,再一步研究基于GPR的时间驱动智能预测模型。
-
GPR用于长序列变形数据分析时,理论上选定核函数后,经初始学习获取最优超参数,在后续的预测环节只需要用待预测的输入量更新协方差矩阵即可,这在开始一段时间内是可行的。然而随着监测时间的推移,新的待预测量与超参数变得越来越不相适应,同时GPR的泛化性能也会降低,如果继续使用原有的超参数进行预测,预测误差会逐渐增大,在适当的时候需要借助训练样本集重新更新超参数。
考虑到变形量在时间域上的相关性,用后继的观测数据加入到初始训练样本集重新求解超参数,同时更新协方差矩阵由此来预测下一个时刻的变形量。每次更新时,由更新前已存在的超参数作为本次更新的初值,这样做的目的是为减少迭代次数,加快求解速度。但这样会存在一个问题,训练样本集会越来越大,不利于超参数的快速求解,应该维持一个数量稳定的能够自动更新的训练样本集,更新模式如下:
假设有d个样本对组成初始训练样本集D={(xi, yi), i=1, 2, …, d},待预测值xd+1作为GPR的输入量,yd+1为GPR的输出量,当实测完yd+1后,将(xd+1,yd+1)作为新的样本对加入到原来的训练集中,截去距离新加入的观测值最远的历史观测值(x1,y1),以保持样本个数的恒定,进一步预测下一个观测时刻变形量,以此类推。将这种更新模式称为“递进-截尾式”动态更新模式,图 1为“递进-截尾式”更新模式示意图。
-
应用GPR进行变形预测时,如何确定训练集的数量d是一个关键因素,如果样本数量少达不到训练的效果,预测精度低,但也并不意味着样本数量多预测精度就高,数量太多反而会影响后续超参数的求解效率,需要综合考虑确定样本数量。
样本数量是否合理的评判依据可以从两方面考虑,一是定性分析法,即通过预测结果是否落入GPR置信区间以及置信区间的宽度定性的评判预测结果的好坏;二是定量分析法,定量的评判方法需要计算有关评价指标,方法如下。
选定d个观测数据组成初始训练样本,同时设定k个已测的数据作为检验数据。将k个实测值和d个样本值作为待预测值输入到GPR预测模型中,得到k个实测值对应的预测值和d个训练样本对应的预测值。由k个预测值和实测值求取残差,用式(8)计算出预测中误差σk用于衡量预测精度;由d个样本值和预测值求取残差,用式(9)计算出训练中误差σd用于衡量训练精度。
逐渐增加d的个数并计算每次d改变后的σk和σd,当二者趋于稳定时,d不再增加,作为最终的学习样本数量值。
$$ {\sigma _k} = \sqrt {\frac{1}{k}\sum\limits_{t = d + 1}^{k + d} {v^2_t} } $$ (8) $$ {\sigma _d} = \sqrt {\frac{1}{d}\sum\limits_{t = 1}^{ d} {v^2_t} } $$ (9) 式中,vt=yt-y*t为对应的残差。
-
无论是传统的几何测量还是先进的GNSS及测量机器人系统,这些监测设备仍需人为布置一定数量的监测点,才能发挥其优势,监测点的观测数据是一个时间序列。
“递进-截尾式”动态更新模式是针对单个监测点,以一般变量xi作为输入量,预测出对应的变形量。如果以观测时间ti作为GPR的输入量,时间作为驱动因子一步步预测监测点未来的变形量,本文将其命名为GPR时间驱动智能预测模型(GPR time-driven intelligent prediction model,GPR-TIPM),图 2是GPR-TIPM的预测流程图。
GPR-TIPM用于预测包含为两个过程,首先要确定训练样本数量d并固定样本数量,然后选择某一时刻作为预测的开始时间进行动态预测,假如从当前时间第10天开始预测未来5 d的变形量,则5 d的时间长度作为1个预测周期,当1个观测周期实测完成后,再预测下一个预测周期内的变形,以此类推,图 3是应用该模型的预测过程图。
-
实验数据来自中煤平朔井工二矿边坡(2#边坡)测量机器人监测系统,在监测区域内沿走向方向在各个“台阶”上布设了63个监测棱镜(监测点)构成5条监测线。图 4是边坡外貌和监测棱镜点布置图。
由于矿山边坡所处的环境复杂,监测点的沉降量并非一直表现为连续下沉,监测点有时会出现上升的现象,实验选取了编号为C1-06的监测点连续60 d的时序观测值作为实验数据。之所以选择C1-06点作为实验对象,是因为C1-06可能受局部环境的影响(如施工作业),使其产生了一个明显的加速下沉,变形较其它点更为复杂,用于实验分析说服力更强。
-
GPR的协方差函数对预测性能有很大影响,即使是相同的数据,不同的协方差函数预测结果不完全相同。就2#边坡的监测数据的初步分析结果为长期的下沉趋势伴随着“波动”特征,时间与变形之间是一复杂的非线性关系。由于GPR的基础理论决定了GPR用于预测只能选择一种核函数,一种核函数难以满足变形曲线多样化的变形特征。GPR可以将不同的协方差函数进行组合运算,本文利用组合式协方差函数模型来适应变形过程的多样化特征。
SE无穷可微使得GPR非常平滑,通过设置SE较小的特征长度参数来适应观测数据的“波动”特征,但有时用过渡光滑的假设条件模拟许多物理过程是不现实的,可以考虑使用马特恩(Matérn 32)协方差函数。因此,实验最终选择Matérn 32+ SE作为后续监测点变形预测的核函数模型。
表 1是Matérn 32协方差函数和组合式Matérn 32+SE协方差函计算的训练中误差和平均相对误差(MRE),同时计算了每种模型的计算耗时,显然Matérn 32+SE的效果更好,但计算效率下降。图 5是输出的变形曲线,SE输出的曲线比Matérn 32+SE光滑,但置信区间宽。
表 1 不同核函数误差统计结果
Table 1. The Results of Different Kernel Function
评价指标 协方差函数模型 SE Matérn 32 Matérn 32+SE σd/mm 0.50 0.34 0.3 MRE/% 9.0 5.7 5.0 耗时/s 0.62 0.52 1.9 具体选择哪种协方差函数,需要分析观测数据自身的特点来做出选择,如果观测数据呈现出明显的线性特征,线性协方差函数无疑是理想的选择;如果观测数据表现出很强的周期性特点,显然要选用周期性核函数。
-
选择d为7的初始样本数量进行学习,检验数据集为第13~17天,即k=5,根据图 2中的GPR-TIPM预测流程进行计算,随着d的增加,置信区间宽度有所减小,另外当d的值增加到12时,σk和σd的值基本趋于稳定,最终的学习样本数量取12。图 6是样本数d=7和d=12预测输出结果。
-
实验初始样本选用1~12 d的实测值进行训练,以Matérn 32+SE作为核函数,选用3 d为1个预测周期, 连续预测第13~60天的累积沉降量,图 7是GPR-TIPM输出的未来3 d的时间-沉降曲线。
从图可以看出,距离学习样本集最近的未来第1天的预测值与实测值的偏差最小,第3天最大。未来第1天的预测时间是从第13~58天,有91%的预测值位于置信区间内;未来第2天的预测时间是从第14~59天,有85%的预测值落入置信区间内;未来第3天的预测时间是从第15~60天,有78%的预测值落入置信区间内。
为了客观定量地评价GPR-TIPM的预测性能,实验计算了预测中误差σk和平均相对误差MRE来衡量预测精度,计算结果见表 2。
表 2 预测精度对比
Table 2. The Comparison of Prediction Accuracy
指标 第1天 第2天 第3天 MRE/% 13 16 20 σk/mm 0.9 1.1 1.4 实验结果表明距离学习样本越远,其预测精度越低,如果没有新的实测值更新学习样本,预测误差会越来越大。与其相对的一面,距离学习样本越近,其预测误差总体上越小,但有个别预测周期的第1天的预测预测误差会大一些。预测误差较大的是发生“抖动”较大的位置,如第38天的相对误差为48%。理论上,GPR-TIPM模型可以预测的更远,但预测结果的可靠性程度会降低。
-
本文首先研究了GPR用于变形预测时“递进-截尾式”超参数动态更新模式,进一步给出了训练样本集的选择方法,在此基础上提出GPR-TIPM非线性智能预测模型。GPR是在函数空间内进行贝叶斯推理的,使得GPR-TIPM能够输出光滑预测中心线和置信区间。GPR-TIPM对建模的样本数据要求不苛刻,能够适应于非平稳数据。理论上GPR-TIPM能够进行长远预测,由于输入量与超参数不适应,使得其预测精度会降低。
应用GPR-TIPM对2#边坡的短期内的变形进行预测,并用实测值分析其预测精度,结果表明GPR-TIPM用于变形预测是可行的。GPR-TIPM预测模型仅以时间作为驱动因子,是一维预测模型,GPR-TIPM没有顾及监测数据的自相关性,这是其不足的一面。
Deformation Intelligent Prediction Model Based on Gaussian Process Regressionand Application
-
摘要: 岩体或建构筑物的变形通常具有复杂性和非线性等特性,一般的回归模型难以精确地进行回归预测,应用高斯过程回归理论对变形监测数据呈现出的非线性特征进行时间序列分析。考虑到监测数据的不断更新和累积,以及超参数与样本集的适应性,首先研究了“递进-截尾式”超参数自动更新模式和训练样本集的选择方法;在此基础上构建了以时间作为输入项的高斯过程回归变形智能预测模型(GPR-TIPM);将该模型应用于矿山边坡监测点非线性时间序列分析中,通过分析变形趋势,最终采用Matérn 32和平方指数协方差函数相加的方式进行核函数组合。实验结果表明,采用组合核函数的预测性能较单一核函数有所改善,该方法提高了模型的泛化能力,GPR-TIPM模型在短期内的预测效果较理想。Abstract: The deformation of building structures or rock mass usually has the features of complexity and nonlinearity that the general regression model cannot accurately predict. In this paper, gaussian process regression(GPR) theory is applied in time series analysis of nonlinear deformation monitoring data. Considering the unceasing updates and massive accumulation of monitoring data, the hyper-parameter and the adaptability of sample set, a "progressive~truncation type" hyper-parameter automatic update mode and selection method for training sample set was developed. On this basis, a GPR time-driven deformation intelligent prediction model(GPR-TIPM) was constructed. This model was applied to the nonlinear time series analysis of monitoring points on a mine slope. By analyzing the deformation trend, a composite kernel optimization method including the "Matérn32" and square exponential covariance kernel function is proposed. The experimental results showed that the prediction performance of the combined kernel function is better than that of the single kernel function, and improved the generalization ability of the model The prediction effect of GPR-TIPM model is better in the short term.
-
Key words:
- deformation monitoring /
- gaussian process regression /
- intelligent predict /
- time series /
- mine slope
-
表 1 不同核函数误差统计结果
Table 1. The Results of Different Kernel Function
评价指标 协方差函数模型 SE Matérn 32 Matérn 32+SE σd/mm 0.50 0.34 0.3 MRE/% 9.0 5.7 5.0 耗时/s 0.62 0.52 1.9 表 2 预测精度对比
Table 2. The Comparison of Prediction Accuracy
指标 第1天 第2天 第3天 MRE/% 13 16 20 σk/mm 0.9 1.1 1.4 -
[1] Stojanovic B, Milivojevic M, Ivanovic M. Adaptive System for Dam Behavior Modeling Based on linear Regression and Genetic Algorithms[J]. Advances in Engineering Software, 2013, 65:182-190 doi: 10.1016/j.advengsoft.2013.06.019 [2] 卢骏, 戴吾蛟, 章浙涛.大坝变形变系数回归建模[J].武汉大学学报·信息科学版, 2015, 40(1):139-142 http://ch.whu.edu.cn/CN/abstract/abstract3171.shtml Lu Jun, Dai Wujiao, Zhang Zhetao.Modeling Dam Deformatin Using Varying Goefficient Regression[J]. Geomatics and Information Science of Wuhan University, 2015, 40(1):139-142 http://ch.whu.edu.cn/CN/abstract/abstract3171.shtml [3] Grelle G, Guadagno I F M. Regression Analysis for Seismic Slope Instability Based on a Double Phase Viscoplastic Sliding Model of the Rigid Block[J]. Landslides, 2013, 10(5):583-597 doi: 10.1007/s10346-012-0350-8 [4] 陈晓鹏, 张强勇, 刘大文.边坡变形统计回归分析模型及应用[J].岩石力学与工程学报, 2008, 27(S2):3673-3678 http://www.doc88.com/p-95325619553.html Chen Xiaopeng, Zhang Qiangyoun, Liu Dawen. Deformation Statistical Regressiong Analysis Model of Slope and Its Application[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(S2):3673-3678 http://www.doc88.com/p-95325619553.html [5] Wang Q, Wang C, Xie R. An Improved SCGM(1, m) Model for Multi-point Deformation Analysis[J]. Geosciences Journal, 2014, 18(4):477-484 doi: 10.1007/s12303-014-0012-z [6] 黄声享, 尹晖, 蒋征.变形监测数据处理[M].武汉:武汉大学出版社, 2012 Huang Shengxiang, Yin Hui, Jiang Zheng.Deformation Monitoring Data Processing[M]. Wuhan:Wuhan University Press, 2010 [7] 武雪玲, 任福, 牛瑞卿.多源数据支持下的三峡库区滑坡灾害空间智能预测[J].武汉大学学报·信息科学版, 2013, 38(8):963-968 http://ch.whu.edu.cn/CN/abstract/abstract2728.shtml Wu Xueling, Ren Fu, Niu Ruiqing. Spatial Intelligent Prediction of Landslide Hazard Based on Multi-source Data in Three Corges Reservoir Area[J]. Geomatics and Information Science of Wuhan University, 2013, 38(8):963-968 http://ch.whu.edu.cn/CN/abstract/abstract2728.shtml [8] Li L, Huang G. Procedia Engineering Prediction of Goaf Settlement with Time Sequence of Wavelet Neural Network[J]. Procedia Engineering, 2011, 15:4723-4727 doi: 10.1016/j.proeng.2011.08.884 [9] Suwansawat S, Einstein H H. Artificial Neural Networks for Predicting the Maximum Surface Settlement Caused by EPB Shield Tunneling[J]. Tunnelling and Underground Space Technology, 2006, 21:133-150 doi: 10.1016/j.tust.2005.06.007 [10] 张正禄, 黄全义, 文鸿雁.工程的变形监测分析与预报[M].北京:测绘出版社, 2007 Zhang Zhenglv, Huang Quanyi, Wen Hongyan. Deformation Monitoring Analysis and Prediction for Engineering Constructions[M]. Beijing:Surveying and Mapping Press, 2007 [11] Li P, Tan Z, Yan L, et al. Time Series Prediction of Mining Subsidence Based on a SVM[J]. Mining Science and Technology, China University of Mining & Technology, 2011, 21(4):557-562 https://www.sciencedirect.com/science/article/pii/S1674526411000998 [12] 李德江, 花向红, 李涛.基于支持向量机的建筑物沉降预测模型研究[J].测绘工程, 2009, 18(3):29-31 http://industry.wanfangdata.com.cn/yj/Detail/Periodical?id=Periodical_chgc200903008 Li Dejiang, Hua Xianghong, Li Tao.Researchon Predictive Modeling of Bulding Settlement Based on Suport Vector Mechine[J]. Engineering of Surveying and Mapping, 2009, 18(3):29-31 http://industry.wanfangdata.com.cn/yj/Detail/Periodical?id=Periodical_chgc200903008 [13] Zhang H, Wang Y, Li Y. SVM Model for Estimating the Parameters of the Probability-integral Method of Predicting Mining Subsidence[J]. Mining Science and Technology (China), China University of Mining and Technology, 2009, 19(3):385-388 doi: 10.1016/S1674-5264(09)60072-7 [14] Li S, Zhao H, Ru Z. Deformation Prediction of Tunnel Surrounding Rock Mass Using CPSO-SVM Model[J]. Journal of Central South University, 2012, 19(11):3311-3319 doi: 10.1007/s11771-012-1409-3 [15] 潘平.基于小波神经网络理论的边坡位移预测[J].成都理工大学学报(自然科学版), 2006, 33(2):176-180 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=cdlg200602009&dbname=CJFD&dbcode=CJFQ Pan Ping.Slope Displacement Forecase Based Onwavelet Neural Network[J]. Journal of Chengdu University of Technology:Sci & Technolgy, 2006, 33(2):176-180 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=cdlg200602009&dbname=CJFD&dbcode=CJFQ [16] 何志昆, 刘光斌, 赵曦晶.高斯过程回归方法综述[J].控制与决策, 2013, 28(8):1121-1129 http://www.wenkuxiazai.com/doc/fb36bd03d4d8d15abe234eb2.html He Zhikun, Liu Guangbin, Zhao Xijing. Overview of Gaussian Process Regression[J]. Control and Decision, 2013, 28(8):1121-1129 http://www.wenkuxiazai.com/doc/fb36bd03d4d8d15abe234eb2.html [17] Su G S. Gaussian Process-based Dynamic Response Surface Method for Estimating Slope Failure Probability[J]. Yantu Lixue/Rock and Soil Mechanics, 2014, 35(12):3592-3601 https://www.sciencedirect.com/science/article/pii/S016747301730214X [18] 张研. 地下工程岩体非线性行为预测识别的高斯过程模型与动态智能反馈分析[D]. 南宁: 广西大学, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10593-1013359112.htm Zhang Yan. Gaussina Process Model for Forecasting Andidentifying Nonlinear Behaviorofunderground Engineering Rockmass Anddynamic Intelligent Feedback Analysis[D]. Nanning: Guangxi University, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10593-1013359112.htm [19] 苏国韶, 肖义龙.边坡可靠度分析的高斯过程方法[J].岩土工程学报, 2011, 33(6):916-920 http://www.wenkuxiazai.com/doc/b4caa98a71fe910ef12df8af.html Su Guoshao, Xiao Yilong. Gaussian Process Method for Slope Reliability Analysis[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(6):916-920 http://www.wenkuxiazai.com/doc/b4caa98a71fe910ef12df8af.html [20] 夏自能. 边坡位移非线性时间序列的髙斯过程预测方法[D]. 南宁: 广西大学, 2012 http://cdmd.cnki.com.cn/Article/CDMD-10593-1012495204.htm Xia Zineng. Forecasting Method of Nonlinear Time Series of Slopedeeormat Using Gaussian Process[D]. Nanning: Guangxi University, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10593-1012495204.htm [21] Rasmussen C E, Williams C K. Gaussian Processes for Machine Learning[M]. Cambridge, MA:MIT Press, 2006 -