Processing math: 0%

基于改进极限学习机模型的盾构掘进引发地表最大沉降预测

阮永芬, 邱龙, 乔文件, 闫明, 郭宇航. 基于改进极限学习机模型的盾构掘进引发地表最大沉降预测[J]. 水文地质工程地质, 2023, 50(5): 124-133. doi: 10.16030/j.cnki.issn.1000-3665.202210007
引用本文: 阮永芬, 邱龙, 乔文件, 闫明, 郭宇航. 基于改进极限学习机模型的盾构掘进引发地表最大沉降预测[J]. 水文地质工程地质, 2023, 50(5): 124-133. doi: 10.16030/j.cnki.issn.1000-3665.202210007
RUAN Yongfen, QIU Long, QIAO Wenjian, YAN Ming, GUO Yuhang. Prediction of the maximum ground settlement caused by shield tunneling based on the improved limit learning machine model[J]. Hydrogeology & Engineering Geology, 2023, 50(5): 124-133. doi: 10.16030/j.cnki.issn.1000-3665.202210007
Citation: RUAN Yongfen, QIU Long, QIAO Wenjian, YAN Ming, GUO Yuhang. Prediction of the maximum ground settlement caused by shield tunneling based on the improved limit learning machine model[J]. Hydrogeology & Engineering Geology, 2023, 50(5): 124-133. doi: 10.16030/j.cnki.issn.1000-3665.202210007

基于改进极限学习机模型的盾构掘进引发地表最大沉降预测

  • 基金项目: 云南省重点研发计划(社会发展领域)项目(2018BC008)
详细信息
    作者简介: 阮永芬(1964-),女,博士,教授,主要从事工程地质与岩土工程的研究。E-mail:rryy64@163.com
    通讯作者: 邱龙(1998-),男,硕士研究生,主要从事工程地质与岩土工程的研究。E-mail:1034403813@qq.com
  • 中图分类号: TU478

Prediction of the maximum ground settlement caused by shield tunneling based on the improved limit learning machine model

More Information
  • 城市地铁盾构施工引发的地面过大变形会严重影响周边构筑物的正常使用,甚至引发工程事故。针对传统预测方法中的数据维度过大容易导致精度降低、计算复杂等问题,提出了一种基于主成分分析(principal component analysis,PCA)算法和哈里斯鹰优化(Harris Hawks optimization,HHO)算法的极限学习机(extreme learning machine,ELM)预测模型。在地质、几何及盾构参数中初选14个影响因子,利用PCA算法在14维数组中分离和提取5个主成分变量作为模型的输入,利用HHO优化ELM模型的输入层权值和隐含层阈值参数,得到预测模型的最优解。以昆明轨道交通五号线怡心桥站—广福路站隧道区间监测数据进行仿真验证,并将该模型与BP神经网络、RBF、未优化的ELM模型进行对比分析。结果表明:PCA-HHO-ELM预测模型的均方根误差为0.143 5、平均绝对误差为0.026 2、决定系数为0.959 6,相较于其他模型,该模型具有更优的预测性能;与未优化的ELM模型相比,HHO算法能够提高ELM模型的预测精度和泛化能力。PCA-HHO-ELM模型能可靠预测盾构诱发的地表最大沉降,可为类似变形预测提供一种更为可行的新思路。

  • 盾构法因其施工速度快、适应性强、环境污染小等优点,已成为目前城市轨道交通的主要施工方法[13]。但盾构机在掘进过程中会对周围土体产生附加载荷,不可避免地扰动周围土体,进而诱发地表沉降,如果局部地表沉降过大,将严重影响周围结构的正常使用,甚至引发工程事故。所以,研究各类诱发因素与地表沉降间的关系,减小掘进过程带来的地表沉降已然成为隧道施工的重中之重。当前,针对盾构隧道施工引起地面沉降这一岩土工程问题,众多学者进行了深入研究,较有代表性的研究方法有理论法[47]、经验公式法[89]、数值模拟[1014]等。但是严格来讲,上述研究方法均属于传统理论范畴,在参数选取及解决非线性问题等方面均具有很大局限性。

    为解决掘进断面的地质条件、施工过程中产生的各类新参数以及持续时间的影响,近年来国内外学者提出了许多盾构隧道地表沉降预测的机器学习方法。如Zhang等[15]对比了误差反向传播(back propagaton,BP)、径向基函数(radial basis function,RBF)、广义回归神经网络(generalized regression netural network,GRNN)3种神经网络模型对隧道最大地面沉降的预测效果,提出了一种将土层力学特性与几何参数相结合的指标计算方法;赵楠等[16]选用粒子群优化(particle swarm optimization,PSO)算法和遗传算法(genetic algorithm,GA)对支持向量机参数进行优化处理,针对多步算法产生的较大偏差,结合实际工程,验证了融合长短时记忆(long short-term memory,LSTM)神经网络支持向量机模型的优越性;李建生等[17]结合Peck公式与遗传算法的相关优点,提出了一种用融合理论公式预测隧道地面沉降的智能学习方法。随着网络技术的发展,越来越多的新兴优化算法不断提出。Chen等[18]基于3种神经网络,提出了一种定义输入参数物理意义的改进指标,用于量化地质参数,提高了人工神经网络模型的预测精度;乔金丽等[19]采用变步长的方法构建BP网络模型,经广州地铁二号线实测沉降检验,其预测效果较好;Wang等[20]提出了基于小波核函数的平滑相关向量机,预测隧道开挖引起的地表沉降;陈仁鹏等[21]通过PSO算法确定了BP神经网络和随机森林算法的超优参数,提出了一种预测盾构隧道纵向及最大地表沉降的方法;宫思艺等[22]在XGBoost地层动态识别模型的基础上,利用盾构施工参数推断地层变化,提出基于BP-SVR的地面沉降预测融合模型;林荣安等[23]为提高由盾构施工引起的软硬不均地层地表沉降预测的准确性,建立基于粗糙集-支持向量回归的地表沉降预测模型;李洛宾等[24]研究处理盾构掘进的序列化数据时,发现循环神经网络的沉降预测结果优于BP神经网络。上述预测方法均取得了一定的效果,但数据维度过大导致精度低、模型学习速度慢、计算复杂,并且由于影响盾构施工地表沉降因素的复杂性、各类方法的缺陷性等问题,在实际应用中仍有一定的局限性。

    极限学习机(extreme learning machine,ELM)是一种典型的单隐层前馈神经网络(single-hidden layer feed-forward netural network,SLFN)[2526],与BP神经网络等传统学习算法相比,ELM具有泛化性能好、学习速度快、预测精度高等优点[27]。但ELM模型输入层权值和隐含层阈值具有随机性,在不确定条件下可能导致模型预测精度降低。为解决这一问题,本文引入哈里斯鹰优化(Harris Hawks optimization,HHO)算法 [28]对ELM中的权值和阈值参数进行优化,弥补ELM中存在的缺陷。针对上述预测方法中存在的数据维度过大、模型学习速度慢、精度低等问题,本文提出了结合主成分分析(principal component analysis,PCA)算法、HHO算法和ELM模型的盾构隧道地表最大沉降预测方法,并以昆明轨道交通五号线怡心桥站—广福路站隧道区间盾构掘进实测数据为例,通过仿真训练和对比分析验证模型的预测性能。

    昆明轨道交通五号线怡心桥站—广福路站(以下简称“怡—广站”)隧道区间起点里程右DK11+911.523,终点里程右DK13+056.149,全长1 144.626 m。线路出怡心桥站后,沿滇池路西北侧敷设,绕避怡心桥、采莲桥桩基后回归滇池路下方,于滇池路广福路路口进入广福路站。本段线路设置曲线半径分别为350 m及1000 m的平面曲线。区间隧道埋深9.71~21.35 m,底板标高1867.08~1872.54 m,隧道洞径6.20 m,采用盾构法施工。研究数据取自怡—广站DK11+966.221—DK12+966.733区段,间隔5 m设置一个监测点,包含隧道轴线处地面沉降监测点位200个。模型所使用的数据集由输入参数和沉降值组成。ELM模型输入是序列形式,故输入值是由沉降监测点前后各2环(共4环)的数据组成,而ELM模型的输出只包括一环。

    数据集中测量沉降值取掘进过程中记录数据的最大沉降值,模型的输入参数(A1A14)包括:地质参数、几何参数、盾构参数3个部分。数据集数据统计信息如表1

    表 1.  数据集数据统计信息表
    Table 1.  Data set statistics table
    参数 最小值 最大值
    地质参数 地下水位(A1)/m 1.50 1.92
    浮密度(A2)/(g·cm–3 0.54 0.77
    压缩模量(A3)/MPa 3.37 4.42
    泊松比(A4 0.28 0.33
    侧压力系数(A5 0.49 0.59
    几何参数 埋深(A6)/m 18.4 25.3
    盾构参数 土压力(A7)/MPa 1.8 2.7
    总推力(A8)/kN 8 760 16 280
    掘进速度(A9)/(mm·min–1 61 94
    出土量(A10)/ m3 45.8 48.0
    刀盘扭矩(A11)/(N·m) 1 200 2 275
    拼装时间(A12)/h 0.35 0.85
    注浆量(A13)/m3 3.60 4.40
    注浆压力(A14)/bar 3.0 3.8
    地表最大沉降/mm −15.69 10.44
     | Show Table
    DownLoad: CSV

    地质参数:怡—广站区间盾构隧道穿越范围内的土层主要包括粉质黏土、粉砂、黏质粉土、泥炭质土。为综合反映各断面不同土质对沉降的影响,选取地下水位、浮密度、压缩模量、泊松比及侧压力系数作为地质参数。其中,浮密度为隧道上方平均加权浮密度,压缩模量、泊松比及侧压力系数取掌子面地层加权平均值。

    几何参数:盾构掘进过程中对地面沉降产生影响的有隧道埋深、盾构机直径及盾尾间隙等。因盾构掘进中,其直径和盾尾间隙是常数,对模型训练没有影响,故仅考虑隧道埋深。

    盾构参数:盾构掘进过程中,会产生大量包括盾构机姿态、注浆参数及刀盘状态在内的64项相关参数。结合文献[29]和实际施工过程,选取8项盾构参数作为输入:土压力、盾构总推力、掘进速度、出土量、刀盘扭矩、管片拼装时间、注浆量、注浆压力。

    PCA是一种常用的数据降维算法之一,其原理是在高维数组中分离和提取多个相互正交的特征变量,充分反映原始数据的差异性,在尽量减少原始指标信息损失的同时,减少待分析的指标[3031]。PCA数据降维通过重构主成分因子表达式,将原始的N维特征映射到低维的K维特征(K<N)中,得到K个主成分变量。

    设原始N个影响因子构成的N维特征向量 {\boldsymbol{X}} = {({x_1},{x_2}, \cdots ,{x_N})^{\rm{T}}} ,则X的协方差矩阵D(X)为:

    {\boldsymbol{D}}({\boldsymbol{X}}) = {(c_{ij})_{N \times N}} = \left( {\begin{array}{*{20}{c}} {{c_{11}}}&{{c_{12}}}& \cdots &{{c_{1N}}} \\ {{c_{21}}}&{{c_{22}}}& \cdots &{{c_{2N}}} \\ \vdots & \vdots & \vdots & \vdots \\ {{c_{N1}}}&{{c_{N2}}}& \cdots &{{c_{NN}}} \end{array}} \right) (1)

    式中:cij——N维特征构成的协方差因子。

    通过QTD(X)Q变换得到D(X)特征值{\boldsymbol{\lambda}} = \lambda_1,\lambda_2, \cdots , \lambda_N (\lambda_1 \geqslant \lambda_2 \geqslant \cdots \geqslant \lambda_N),最后得到重构后的综合向量Z,其表达式为:

    Z_i = {L_{1i}}{x_1} + {L_{2i}}{x_2} + \cdots + {L_{Ni}}{x_N} (2)

    式中:L——特征值λ的正交特征向量,L=(Lij)。

    为全面反映原始N维特征且不减少指标信息损失,一般取方差最大的前KK<N)个变量作为主成分因子,且定义累计贡献率(即累计方差)大于或等于80%时为筛选出的K个主成分因子,将得到的K个主成分因子通过式(2)计算,即可得到重构后的Z1ZK

    SLFN大多采用梯度下降法迭代,容易导致模型训练慢、计算复杂等问题。与之相比,ELM模型的输入层权值和隐含层阈值具有随机性,模型收敛速度快,训练时只需设置隐含层神经元个数,便可通过迭代获得全局最优解,所以与传统SLFN相比,ELM模型的泛化性更好、学习速度更快。其网络拓扑结构如图1所示。

    图 1.  ELM算法模型结构图
    Figure 1.  Structure diagram of ELM algorithm model

    图1知,ELM模型的输入层、隐含层及输出层之间全连接,分别有KLS个神经元节点。输入层的K个神经元节点,对应2.1节中的K个主成分变量;若模型用于分类问题,则S>1,若模型用于回归预测,则S=1。设Q个训练集输入矩阵X和输出矩阵Y分别为 {\boldsymbol{X}} = {\left( {{x_{i1}},{x_{i2}}, \cdots ,{x_{iK}}} \right)^{\rm{T}}} \in {{\boldsymbol{R}}_K} {\boldsymbol{Y}} = {\left( {{y_{i1}},{y_{i2}}, \cdots ,{y_{iS}}} \right)^{\rm{T}}} \in {{\boldsymbol{R}}_S} ,假定任意小误差激活函数gx),则网络输出矩阵Z

    {Z}_j = \sum\limits_{i = 1}^L {{{{\beta}} _i}g({{\boldsymbol{\omega}} _j} \cdot {{\boldsymbol{X}}_j} + {{{b}}_j})} ,j = 1,2, \cdots ,Q (3)

    式中: {{\boldsymbol{\omega}} _i} = \left( {{\omega _{i1}},{\omega _{i2}}, \cdots ,{\omega _{iK}}} \right) ——输入层权值矩阵;

    {\boldsymbol{\beta}} = \left( {\beta_1,\beta_2, \cdots ,\beta_S} \right) ——输出层权值矩阵;

    {\boldsymbol{b}} = \left( {b_1,b_2, \cdots ,b_L} \right)——隐含层阈值矩阵。

    G为ELM的隐含层输出矩阵,则式(3)可简化为

    {\boldsymbol{G \cdot \beta}} = {{\boldsymbol{Z}}^{\rm{T}}} (4)

    Q较大时,为减少计算量,L一般取小于Q的数,这时ELM的误差就可逼近一个无限小的值εε>0),如式(5)所示。所以,当激活函数任意无限小且可微时,ELM参数就不再需要全部迭代调整,即当ωb随机设定后,ELM就可直接运行且不再修改权值和阈值参数。综上,ELM模型的训练结果可等效为求解方程(6)的最小二乘解。

    \sum\limits_{j = 1}^Q {\left\| {Z_j - Y_j} \right\|} = 0 (5)
    \mathop {\min }\limits_{\boldsymbol{\beta}} \left\| {{\boldsymbol{G \cdot \beta }}- {{\boldsymbol{Z}}^{\rm{T}}}} \right\| (6)

    式(6)的解为:

    \overline {\boldsymbol{\beta}} = {{\boldsymbol{G}}^ + } \cdot {{\boldsymbol{Z}}^{\rm{T}}} (7)

    式中:G+——G的Moore-Penrose广义逆。

    HHO算法是Heidari等[28]基于哈里斯鹰的捕食行为,于2019年提出的一种群算法,该算法需要调节的参数较少,具有较强的全局搜索能力。具体包括探索阶段、转换阶段及开发阶段,具体捕食行为如下。

    在此阶段,哈里斯鹰处于侦察猎物状态,若未找到猎物,需不断按式(8)更新位置寻找猎物。

    X(t + 1) = \left\{ \begin{split} &\left. {X_{\rm{rand}}(t) - r_1} \right|\left. {X_{\rm{rand}}(t) - 2r_2X(t)} \right|&q \geqslant 0.5 \\ &[X_{\rm{rabbit}}(t) - X_{\rm{m}}(t)] - r_3[lb + r_4(ub - lb)]&q \lt 0.5 \\ \end{split} \right. (8)

    式中:r1r2r3r4——[0,1]之间的随机数;

    q——用于判定选择采用的策略;

    ublb——搜索空间的上、下限;

    Xrand(t) ——哈里斯鹰个体的随机位置;

    Xrabbit(t)——猎物当前位置;

    X(t)、X(t+1) ——当前和下一次更新时的个体位置;

    Xm(t)——当前种群所有个体的平均位置,表达式为:

    X_{\rm{m}}(t) = \frac{1}{N}\sum\limits_{i = 1}^N {X_i(t)} (9)

    式中:Xi(t)——种群中第i个个体的位置;

    N——种群数量。

    在此阶段,猎物体能极大消耗,哈里斯鹰根据猎物的逃逸能量(C)调整自身位置,将全局搜索范围变为局部搜索,C按式(10)计算:

    C = 2C_0\left(1 - \frac{t}{{M_{\text{a}}}}\right) (10)

    式中:C0——猎物的初始能量,每次迭代时在[−1,1]中 随机自动更新;

    Ma——最大迭代数;

    t——当前迭代数。

    当|C|≥1时,猎物逃逸能量高,猎物容易逃脱,则返回2.3.1的搜索阶段;当|C|<1时,猎物能量较低,则进入下一阶段。

    在此阶段,哈里斯鹰通过4种策略对猎物进行围捕,其数学模型见式(11)—(17),设猎物的逃脱概率为r,当r<0.5时猎物逃脱,当r≥0.5时逃脱失败。

    (1)软包围策略

    X(t + 1) = \Delta X(t) - \left. E \right|JX_{\rm{rabbit}}(t) - \left. {X(t)} \right| (11)
    \Delta X(t) = X_{\rm{rabbit}}(t) - X(t) (12)

    式中:J——猎物逃逸的跳跃距离,逃逸时在0~2之间随机选择;

    ΔX——哈里斯鹰的最优个体与当前个体的差值;

    E——猎物动态逃逸能量。

    (2)硬包围策略

    X(t + 1) = X_{\rm{rabbit}}(t) - \left. E \right|\left. {\Delta X(t)} \right| (13)

    (3)渐进式快速俯冲软包围策略

    X(t + 1) = \left\{ \begin{gathered} Y\quad f(Y) \lt f(X(t)) \\ Z\quad f(Z) \lt f(X(t)) \\ \end{gathered} \right. (14)
    Y = X_{\rm{rabbit}}(t) - \left. E \right|JX_{\rm{rabbit}}(t) - \left. {X(t)} \right| (15)
    Z = Y + {\boldsymbol{S}} \cdot LF(D) (16)

    式中:f()——适应度函数;

    D——已知种群维度;

    S——D维随机矩阵;

    LF()——莱维(Levy)飞行函数。

    (4)渐进式快速俯冲硬包围策略

    Y = X_{\rm{rabbit}}(t) - \left. E \right|JX_{\rm{rabbit}}(t) - \left. {X_{\rm{m}}(t)} \right| (17)

    由于ELM模型的输入层权值和隐含层阈值都是随机产生的,这些参数在不确定的情况下可能导致ELM模型预测精度降低或者模型产生过拟合现象,因此本文利用HHO算法对ELM模型的权值和阈值参数进行优化,寻求全局最优解,从而使ELM具有更好的学习和预测能力。基于以上原理,本文提出了一种集ELM、PCA、HHO于一体的混合人工智能模型。ELM用于构建盾构隧道地表最大沉降与其影响因素之间的关系;PCA用于数据降维,简化ELM网络拓扑结构;HHO用于优化ELM的权值和阈值参数。PCA-HHO-ELM模型的计算步骤如下,流程图如图2所示。

    图 2.  PCA-HHO-ELM模型流程图
    Figure 2.  Flow chart of PCA-HHO-ELM model

    步骤1:对数据进行预处理,剔除非正常数据,利用PCA算法将数据降维,提取出K个主成分因子。

    步骤2:为使数据具有代表性,将提出的主成分因子及地表最大沉降值随机排列,把数据分为训练集和测试集。

    步骤3:对训练集及测试集归一化处理。

    步骤4:利用HHO算法对ELM中的输入层权值和隐含层阈值参数进行优化,重构ELM模型。

    步骤5:迭代更新适应度值及位置,将训练误差作为HHO适应度函数,判定是否满足终止条件,若满足则退出循环,输出寻优结果;否则,重复步骤5。

    步骤6:将测试集输入优化后的ELM中,输出预测结果。

    表1知,盾构隧道地表最大沉降的影响因素具有14维,属于高维数组,为避免模型学习速度慢、计算复杂、精度低等问题,现利用PCA算法对14维影响参数进行数据降维,设定累计贡献率大于80%的影响参数为主成分因子。运用Matlab编写程序代码,经多次试验分析,最终确定出了5个影响因子,其提取后的变量特征值与累计贡献率如图3所示。

    图 3.  PCA主成分提取变量
    Figure 3.  PCA principal component extraction variables

    图3主成分提取变量结果可知,前5个成分的累计贡献率分别为36.661%、53.114%、66.596%、74.333%、81.510%,到第5成分时累计贡献率达到了81.510%,大于80%。因此,在尽量减少信息损失的同时,前5个成分能够较为全面地反映原始指标体系,即能够将初始的14维变量数据降至5维,并且由5个公因子表征出来。所得到的主成分系数矩阵如表2所示。

    表 2.  主成分系数矩阵
    Table 2.  Principal component coefficient matrix
    影响因子主成分
    12345
    A10.8970.0230.0860.0770.052
    A20.937−0.0770.0430.0410.051
    A3−0.2280.5860.376−0.460−0.163
    A4−0.271−0.259−0.7240.3400.105
    A50.956−0.068−0.130−0.007−0.028
    A60.3760.3350.3650.4920.315
    A70.683−0.093−0.116−0.5200.060
    A8−0.0780.3100.425−0.0530.748
    A9−0.4740.530−0.5820.0190.216
    A10−0.2270.867−0.2470.123−0.128
    A11−0.809−0.112−0.056−0.2560.188
    A12−0.790−0.4640.277−0.0710.080
    A13−0.201−0.6360.0460.0190.247
    A14−0.445−0.0110.6220.384−0.365
     | Show Table
    DownLoad: CSV

    根据表2主成分系数矩阵,赋予A1A14相应权重,得到重构后Z1Z5的表达式,如式(18)—(22)所示。

    \begin{split} Z_1 &= 0.897A_1 + 0.937A_2 - 0.228A_3 - 0.271A_4 + 0.956A_5 \\ &+ 0.376A_6 + 0.683A_7 - 0.078A_8 - 0.474A_9 - 0.227A_{10} \\ &- 0.809A_{11} - 0.790A_{12} - 0.201A_{13} - 0.445A_{14} \end{split} (18)
    \begin{split} Z_2 &= 0.023A_1 - 0.077A_2 + 0.586A_3 - 0.259A_4 - 0.068A_5 \\ & + 0.335A_6 - 0.093A_7 + 0.310A_8 + 0.530A_9 + 0.867A_{10} \\ & - 0.112A_{11} - 0.464A_{12 }- 0.636A_{13} - 0.011A_{14} \\ \end{split} (19)
    \begin{split} Z_3 &= 0.086A_1 + 0.043A_2 + 0.376A_3 - 0.724A_4 - 0.130A_5 \\ & + 0.365A_6 - 0.116A_7 + 0.425A_8 - 0.582A_9 - 0.247A_{10} \\ & - 0.056A_{11} + 0.277A_{12} + 0.046A_{13} + 0.622A_{14} \end{split} (20)
    \begin{split} Z_4 &= 0.077A_1 + 0.041A_2 - 0.460A_3 + 0.340A_4 - 0.007A_5 \\ & + 0.492A_6 - 0.520A_7 - 0.053A_8 + 0.019A_9 + 0.123A_{10} \\ & - 0.256A_{11} - 0.071A_{12} + 0.019A_{13} + 0.384A_{14} \end{split} (21)
    \begin{split} Z_5 &= 0.052A_1 + 0.051A_2 - 0.163A_3 + 0.105A_4 - 0.028A_5 \\ & + 0.315A_6 + 0.060A_7 + 0.748A_8 + 0.216A_9 - 0.128A_{10} \\ & + 0.188A_{11} + 0.080A_{12} + 0.247A_{13} - 0.365A_{14} \end{split} (22)

    表1可知,数据集量纲不同,需对数据进行标准化处理,从而加快极限学习机模型的收敛,并在训练过程中防止梯度爆炸。标准化处理算法将重构的5维主成分变量和地表最大沉降数据进行线性变换使其值域在[−1,1]附近。计算过程如下:

    \tilde x = \frac{{x - \mu }}{\sigma } (23)

    式中:μ、σ——输入数据某一维度的最小值及标准差;

    \tilde x ——标准化处理后的数据。

    为避免连续监测点位之间沉降差异较小,导致训练集沉降范围具有局限性,现将200组经PCA算法降维的数据随机排列,前170组用于模型训练(即训练集),后30组用于模型性能测试(即测试集)。利用HHO算法对ELM的输入层权值和隐含层阈值参数进行寻优,将得到的最优结果作为ELM的最优参数,最大程度提高模型的预测性能。经反复试验,将HHO种群大小设为30,最大迭代次数设为100。

    由前文知,Z1Z5作为模型输入且用于回归预测(S=1),故ELM输入层神经元和输出层神经元分别5和1。选择Sigmoid函数作为激活函数gx),如式(24)所示。由式(24)知gx)无穷可微,满足激活函数定义。根据Kolomogorov定理,SLFN隐含层神经元不小于2K+1,即不小于11个。

    g(x) = \frac{1}{{1 + {{\text{e}}^{ - x}}}} (24)

    为验证模型预测效果,选用平均绝对误差(MAE)、均方根误差(RMSE)和决定系数(R2)3个指标对模型进行评价。RMSE弥补了MSE与目标变量量纲不一致的缺陷,MAE避免了正负误差相抵消,R2用于评价模型与数据的相关性,R2的范围为[0,1],R2值越大,模型拟合效果越好。具体公式如下:

    RMS E = \sqrt {\frac{1}{m}\sum\limits_{i = 1}^m {{{\left[ {h(x_i) - y_i} \right]}^2}} } (25)
    MAE = \left| {\frac{1}{m}\sum\limits_{i = 1}^m {\left[ {h(x_i) - y_i} \right]} } \right| (26)
    {R^2} = 1 - \frac{{\displaystyle\sum\limits_{i = 1}^m {{{\left[ {y_i - h(x_i)} \right]}^2}} }}{{\displaystyle\sum\limits_{i = 1}^m {{{(y_i - \overline y )}^2}} }} (27)

    式中:m——样本总数;

    h(xi)——模型的预测沉降值;

    yi——实测沉降值;

    \overline y ——实测沉降数据的平均值。

    经反复试验,将隐含层节点数设为20,预测效果较好。在PCA-HHO-ELM模型仿真训练后,对30组测试样本进行测试,训练及测试结果如图4所示。

    图 4.  PCA-HHO-ELM模型预测结果
    Figure 4.  Prediction results of PCA-HHO-ELM model

    图4(a)中可看出,训练集数据点基本分布在理想拟合直线附近,训练样本的R2为0.941 9、RMSE为0.154 0、MAE为0.011 8,这说明PCA-HHO-ELM模型的预测值与实测值基本一致,模型训练效果较好。从图4(b)中可看出,测试样本数据点同样集中分布于理想拟合直线附近,其R2为0.959 6、RMSE为0.1435、MAE为0.026 2,这说明PCA-HHO-ELM模型具有很好的学习与预测能力,能够有效预测盾构隧道地表最大沉降。

    为进一步验证PCA-HHO-ELM模型的泛化能力,在同等试验条件下将未优化的ELM模型、BP神经网络模型和RBF模型与PCA-HHO-ELM模型进行多次试验分析对比,对比结果如图5所示。

    图 5.  各模型预测结果对比
    Figure 5.  Comparison of prediction results of various models

    图5预测曲线可知,4种预测模型的预测值与实测值整体上变化趋势一致,均能在一定程度上对盾构隧道地表最大沉降进行可靠预测,但在局部区域,PCA-HHO-ELM模型的预测结果较其他3种算法更接近地表沉降实测值,其预测值与实测值咬合得更加紧密。相较于RBF、ELM模型,BP神经网络模型同样具有较高的可靠性。将4类模型预测结果(表3)带入式(22)—(24),得到各模型的预测性能指标,其结果如表4所示。

    表 3.  4类模型预测结果
    Table 3.  Prediction results of four types of models
    组号123456789101112131415
    实测值/mm2.899.865.034.452.0710.440.823.640.4413.429.694.962.541.931.43
    BP预测值/mm1.62−8.493.61−8.92−0.5612.08−0.973.510.12−13.7810.10−3.742.532.240.70
    RBF预测值/mm5.040.492.122.00−2.862.810.314.480.50−11.035.07−11.392.032.010.50
    ELM预测值/mm−0.24−6.174.210.11−0.9412.900.27−0.51−1.51−11.1311.95−4.052.632.120.26
    PCA-HHO-ELM预测值/mm2.11−7.834.93−5.250.149.38−1.907.34−1.23−12.6612.93−4.723.833.46−0.82
    组号161718192021222324252627282930
    实测值/mm0.016.335.2214.1913.741.144.175.311.610.089.901.2410.842.2412.35
    BP预测值/mm−0.064.854.89−14.13−14.250.494.105.931.73−0.05−9.73−0.01−10.461.30−11.83
    RBF预测值/mm0.503.592.56−11.32−7.450.504.475.070.500.50−9.670.50−8.020.58−7.45
    ELM预测值/mm−0.535.725.79−11.76−15.911.342.905.711.65−0.58−9.780.70−8.43−0.12−14.09
    PCA-HHO-ELM预测值/mm−1.336.506.89−13.27−14.072.164.352.434.00−0.99−9.010.13−10.761.18−12.41
     | Show Table
    DownLoad: CSV
    表 4.  模型性能评价结果
    Table 4.  Performance evaluation results of the model
    评价指标RMSEMAER2最大误差/mm平均误差/mm
    PCA-HHO-ELM模型0.143 50.026 20.959 63.701.29
    BP预测模型0.232 30.042 40.894 34.480.82
    RBF模型0.392 10.071 60.698 610.352.52
    ELM预测模型0.571 80.104 40.359 24.561.48
     | Show Table
    DownLoad: CSV

    表4中可看出,本文提出的PCA-HHO-ELM模型的5项评价指标基本优于其余3种预测模型,仅平均误差略大于BP模型,其RMSEMAE分别为0.1435和0.0262,说明该模型更为接近实际沉降值且误差最小,R2高达0.959 6,说明预测值与实测值拟合效果较好,可靠性较高。4种模型的RMSEMAE尽管都很小,但是ELM模型误差是本文模型的近4倍,经HHO优化后,PCA-HHO-ELM模型的单样本最大误差值和平均误差均小于未优化的ELM模型,决定系数也从0.3592提升至0.9596,其预测效果较好。综上所述,4类预测模型预测效果为:PCA-HHO-ELM模型>BP模型>RBF模型>ELM模型,HHO算法能够提高ELM的预测精度,PCA-HHO-ELM模型具有更强的泛化能力。

    (1)利用PCA算法进行数据降维,在尽量减少信息损失的同时,将原始14维影响参数降至5维,根据主成分系数矩阵重构主成分变量,将得到的主成分变量作为模型输入,简化了ELM模型的网络拓扑结构,从而减少所需优化的权值和阈值个数。

    (2)利用PCA-HHO-ELM模型进行盾构沉降预测,训练集及测试集数据点基本分布在理想拟合直线附近,模型预测效果较好。PCA-HHO-ELM预测模型的各项评价指标基本均优于BP神经网络模型、RBF模型和ELM模型。

    (3)经HHO算法优化后,ELM模型精度显著提高,R2从0.3592提升至0.9596,其单样本最大误差和平均误差也得到相应改善。PCA-HHO-ELM模型在一定程度上提高了盾构沉降预测效果,使其具有更优的泛化性能,可为类似隧道工程提供一种更为可行的新思路。

  • 图 1  ELM算法模型结构图

    Figure 1. 

    图 2  PCA-HHO-ELM模型流程图

    Figure 2. 

    图 3  PCA主成分提取变量

    Figure 3. 

    图 4  PCA-HHO-ELM模型预测结果

    Figure 4. 

    图 5  各模型预测结果对比

    Figure 5. 

    表 1  数据集数据统计信息表

    Table 1.  Data set statistics table

    参数 最小值 最大值
    地质参数 地下水位(A1)/m 1.50 1.92
    浮密度(A2)/(g·cm–3 0.54 0.77
    压缩模量(A3)/MPa 3.37 4.42
    泊松比(A4 0.28 0.33
    侧压力系数(A5 0.49 0.59
    几何参数 埋深(A6)/m 18.4 25.3
    盾构参数 土压力(A7)/MPa 1.8 2.7
    总推力(A8)/kN 8 760 16 280
    掘进速度(A9)/(mm·min–1 61 94
    出土量(A10)/ m3 45.8 48.0
    刀盘扭矩(A11)/(N·m) 1 200 2 275
    拼装时间(A12)/h 0.35 0.85
    注浆量(A13)/m3 3.60 4.40
    注浆压力(A14)/bar 3.0 3.8
    地表最大沉降/mm −15.69 10.44
    下载: 导出CSV

    表 2  主成分系数矩阵

    Table 2.  Principal component coefficient matrix

    影响因子主成分
    12345
    A10.8970.0230.0860.0770.052
    A20.937−0.0770.0430.0410.051
    A3−0.2280.5860.376−0.460−0.163
    A4−0.271−0.259−0.7240.3400.105
    A50.956−0.068−0.130−0.007−0.028
    A60.3760.3350.3650.4920.315
    A70.683−0.093−0.116−0.5200.060
    A8−0.0780.3100.425−0.0530.748
    A9−0.4740.530−0.5820.0190.216
    A10−0.2270.867−0.2470.123−0.128
    A11−0.809−0.112−0.056−0.2560.188
    A12−0.790−0.4640.277−0.0710.080
    A13−0.201−0.6360.0460.0190.247
    A14−0.445−0.0110.6220.384−0.365
    下载: 导出CSV

    表 3  4类模型预测结果

    Table 3.  Prediction results of four types of models

    组号123456789101112131415
    实测值/mm2.899.865.034.452.0710.440.823.640.4413.429.694.962.541.931.43
    BP预测值/mm1.62−8.493.61−8.92−0.5612.08−0.973.510.12−13.7810.10−3.742.532.240.70
    RBF预测值/mm5.040.492.122.00−2.862.810.314.480.50−11.035.07−11.392.032.010.50
    ELM预测值/mm−0.24−6.174.210.11−0.9412.900.27−0.51−1.51−11.1311.95−4.052.632.120.26
    PCA-HHO-ELM预测值/mm2.11−7.834.93−5.250.149.38−1.907.34−1.23−12.6612.93−4.723.833.46−0.82
    组号161718192021222324252627282930
    实测值/mm0.016.335.2214.1913.741.144.175.311.610.089.901.2410.842.2412.35
    BP预测值/mm−0.064.854.89−14.13−14.250.494.105.931.73−0.05−9.73−0.01−10.461.30−11.83
    RBF预测值/mm0.503.592.56−11.32−7.450.504.475.070.500.50−9.670.50−8.020.58−7.45
    ELM预测值/mm−0.535.725.79−11.76−15.911.342.905.711.65−0.58−9.780.70−8.43−0.12−14.09
    PCA-HHO-ELM预测值/mm−1.336.506.89−13.27−14.072.164.352.434.00−0.99−9.010.13−10.761.18−12.41
    下载: 导出CSV

    表 4  模型性能评价结果

    Table 4.  Performance evaluation results of the model

    评价指标RMSEMAER2最大误差/mm平均误差/mm
    PCA-HHO-ELM模型0.143 50.026 20.959 63.701.29
    BP预测模型0.232 30.042 40.894 34.480.82
    RBF模型0.392 10.071 60.698 610.352.52
    ELM预测模型0.571 80.104 40.359 24.561.48
    下载: 导出CSV
  • [1]

    胡群芳,李敏,王飞. 我国城市地铁隧道盾构施工刀具磨损性分区研究[J]. 现代隧道技术,2016,53(2):26 − 34. [HU Qunfang,LI Min,WANG Fei. Research on regionalization of cutting-tool wear conditions in shield tunnels of urban subways in China[J]. Modern Tunnelling Technology,2016,53(2):26 − 34. (in Chinese with English abstract)

    HU Qunfang, LI Min, WANG Fei. Research on regionalization of cutting-tool wear conditions in shield tunnels of urban subways in China[J]. Modern Tunnelling Technology, 2016, 532): 2634. (in Chinese with English abstract)

    [2]

    路德春,马一丁,王国盛. 近接隧道列车运行时地表振动响应数值模拟[J]. 吉林大学学报(地球科学版),2021,51(5):1452 − 1462. [LU Dechun,MA Yiding,WANG Guosheng. Numerical study on ground surface vibration response under train load in multi adjacent tunnels[J]. Journal of Jilin University (Earth Science Edition),2021,51(5):1452 − 1462. (in Chinese with English abstract)

    LU Dechun, MA Yiding, WANG Guosheng. Numerical study on ground surface vibration response under train load in multi adjacent tunnels[J]. Journal of Jilin University (Earth Science Edition), 2021, 515): 14521462. (in Chinese with English abstract)

    [3]

    XU Qianwei,ZHU Hehua,MA Xianfeng,et al. A case history of shield tunnel crossing through group pile foundation of a road bridge with pile underpinning technologies in Shanghai[J]. Tunnelling and Underground Space Technology,2015,45:20 − 33. doi: 10.1016/j.tust.2014.09.002

    [4]

    蒋彪,皮圣,阳军生,等. 长沙地铁典型地层盾构施工地表沉降分析与预测[J]. 地下空间与工程学报,2016,12(1):181 − 187. [JIANG Biao,PI Sheng,YANG Junsheng,et al. Analysis and prediction of ground surface settlements due to EPB shield tunneling of Changsha metro[J]. Chinese Journal of Underground Space and Engineering,2016,12(1):181 − 187. (in Chinese with English abstract)

    JIANG Biao, PI Sheng, YANG Junsheng, et al. Analysis and prediction of ground surface settlements due to EPB shield tunneling of Changsha metro[J]. Chinese Journal of Underground Space and Engineering, 2016, 121): 181187. (in Chinese with English abstract)

    [5]

    PARK K H. Analytical solution for tunnelling-induced ground movement in clays[J]. Tunnelling and Underground Space Technology,2005,20(3):249 − 261. doi: 10.1016/j.tust.2004.08.009

    [6]

    赖文杰,齐昌广,郑金辉,等. 含分数阶的灰色模型及其在地基沉降预测中的应用[J]. 水文地质工程地质,2019,46(3):124 − 128. [LAI Wenjie,QI Changguang,ZHENG Jinhui,et al. Gray model with fractional order and its application to settlement prediction[J]. Hydrogeology & Engineering Geology,2019,46(3):124 − 128. (in Chinese with English abstract)

    LAI Wenjie, QI Changguang, ZHENG Jinhui, et al. Gray model with fractional order and its application to settlement prediction[J]. Hydrogeology & Engineering Geology, 2019, 463): 124128. (in Chinese with English abstract)

    [7]

    郭社军,郝晓龙,舒国明. Peck公式在小半径曲线隧道沉降分析的应用[J]. 公路,2020,65(8):405 − 408. [GUO Shejun,HAO Xiaolong,SHU Guoming. Application of Peck formula in settlement analysis of small radius curved tunnel[J]. Highway,2020,65(8):405 − 408. (in Chinese).

    GUO Shejun, HAO Xiaolong, SHU Guoming. Application of Peck formula in settlement analysis of small radius curved tunnel[J]. Highway, 2020, 658): 405408. (in Chinese).

    [8]

    PECK R B. Deep excavations and tunneling in soft ground[C]//Mexico:Proceeding of 7th ICSMFE,1969:225 − 290.

    [9]

    范珊珊,郭海朋,朱菊艳,等. 线性回归模型在北京平原地面沉降预测中的应用[J]. 中国地质灾害与防治学报,2013,24(1):70 − 74. [FAN Shanshan,GUO Haipeng,ZHU Juyan,et al. Application of linear regression model for land subsidence prediction in Beijing plain[J]. The Chinese Journal of Geological Hazard and Control,2013,24(1):70 − 74. (in Chinese with English abstract)

    FAN Shanshan, GUO Haipeng, ZHU Juyan, et al. Application of linear regression model for land subsidence prediction in Beijing plain[J]. The Chinese Journal of Geological Hazard and Control, 2013, 24(1): 70 − 74. (in Chinese with English abstract)

    [10]

    KIM C Y,BAE G J,HONG S W,et al. Neural network based prediction of ground surface settlements due to tunnelling[J]. Computers and Geotechnics,2001,28(6/7):517 − 547.

    [11]

    DIAS D,KASTNER R. Movements caused by the excavation of tunnels using face pressurized shields:Analysis of monitoring and numerical modeling results[J]. Engineering Geology,2013,152(1):17 − 25. doi: 10.1016/j.enggeo.2012.10.002

    [12]

    于德海,舒娇娇,秦凯凯. 盾构地铁隧道穿越既有铁路桥的沉降分析[J]. 水文地质工程地质,2020,47(2):148 − 152. [YU Dehai,SHU Jiaojiao,QIN Kaikai. An analysis of the settlement of a shield tunnel passing under the operating railway bridge[J]. Hydrogeology & Engineering Geology,2020,47(2):148 − 152. (in Chinese with English abstract)

    YU Dehai, SHU Jiaojiao, QIN Kaikai. An analysis of the settlement of a shield tunnel passing under the operating railway bridge[J]. Hydrogeology & Engineering Geology, 2020, 472): 148152. (in Chinese with English abstract)

    [13]

    潘涛. 软土地区双线区间盾构隧道施工对周边地表以及建筑物沉降的影响[J]. 水文地质工程地质,2022,49(1):101 − 108. [PAN Tao. Influences of double-track shield tunnel construction on settlements of adjacent ground and buildings in a soft soil area[J]. Hydrogeology & Engineering Geology,2022,49(1):101 − 108. (in Chinese with English abstract)

    PAN Tao. Influences of double-track shield tunnel construction on settlements of adjacent ground and buildings in a soft soil area[J]. Hydrogeology & Engineering Geology, 2022, 491): 101108. (in Chinese with English abstract)

    [14]

    徐明祥,黄强兵,王庆兵,等. 西安地裂缝地段浅埋暗挖地铁隧道施工沉降规律[J]. 水文地质工程地质,2020,47(1):161 − 170. [XU Mingxiang,HUANG Qiangbing,WANG Qingbing,et al. Settlement rules of shallow-buried metro tunnel construction in the Xi’an ground fissure section[J]. Hydrogeology & Engineering Geology,2020,47(1):161 − 170. (in Chinese with English abstract)

    XU Mingxiang, HUANG Qiangbing, WANG Qingbing, et al. Settlement rules of shallow-buried metro tunnel construction in the Xi’an ground fissure section[J]. Hydrogeology & Engineering Geology, 2020, 471): 161170. (in Chinese with English abstract)

    [15]

    ZHANG Pin,CHEN Renpeng,WU Huaina. Real-time analysis and regulation of EPB shield steering using Random Forest[J]. Automation in Construction,2019,106:102860. doi: 10.1016/j.autcon.2019.102860

    [16]

    赵楠,李洁. 基于LSTM-SVM的隧道围岩位移预测[J]. 公路,2021,66(6):404 − 407. [ZHAO Nan,LI Jie. Displacement prediction of tunnel surrounding rock based on LSTM-SVM[J]. Highway,2021,66(6):404 − 407. (in Chinese)

    ZHAO Nan, LI Jie. Displacement prediction of tunnel surrounding rock based on LSTM-SVM[J]. Highway, 2021, 66(6): 404 − 407. (in Chinese)

    [17]

    李建生,阳军生,杨铠,等. 改进遗传算法在浅埋隧道施工倾斜地表沉降预测中的应用[J]. 公路工程,2008,33(6):46 − 49. [LI Jiansheng,YANG Junsheng,YANG Kai,et al. Application of improved genetic algorithm in prediction of inclined ground settlement in shallow tunnel construction[J]. Highway Engineering,2008,33(6):46 − 49. (in Chinese)

    LI Jiansheng, YANG Junsheng, YANG Kai, et al. Application of improved genetic algorithm in prediction of inclined ground settlement in shallow tunnel construction[J]. Highway Engineering, 2008, 336): 4649. (in Chinese)

    [18]

    CHEN Renpeng,ZHANG Pin,KANG Xin,et al. Prediction of maximum surface settlement caused by earth pressure balance (EPB) shield tunneling with ANN methods[J]. Soils and Foundations,2019,59(2):284 − 295. doi: 10.1016/j.sandf.2018.11.005

    [19]

    乔金丽,范永利,刘波,等. 基于改进BP网络的盾构隧道开挖地表沉降预测[J]. 地下空间与工程学报,2012,8(2):352 − 357. [QIAO Jinli,FAN Yongli,LIU Bo,et al. Predicting the surface settlement by shield tunneling based on modified BP network[J]. Chinese Journal of Underground Space and Engineering,2012,8(2):352 − 357. (in Chinese with English abstract)

    QIAO Jinli, FAN Yongli, LIU Bo, et al. Predicting the surface settlement by shield tunneling based on modified BP network[J]. Chinese Journal of Underground Space and Engineering, 2012, 82): 352357. (in Chinese with English abstract)

    [20]

    WANG Fan,GOU Biancai,QIN Yawei. Modeling tunneling-induced ground surface settlement development using a wavelet smooth relevance vector machine[J]. Computers and Geotechnics,2013,54:125 − 132. doi: 10.1016/j.compgeo.2013.07.004

    [21]

    陈仁朋,戴田,张品,等. 基于机器学习算法的盾构掘进地表沉降预测方法[J]. 湖南大学学报(自然科学版),2021,48(7):111 − 118. [CHEN Renpeng,DAI Tian,ZHANG Pin,et al. Prediction method of tunneling-induced ground settlement using machine learning algorithms[J]. Journal of Hunan University (Natural Sciences),2021,48(7):111 − 118. (in Chinese with English abstract)

    CHEN Renpeng, DAI Tian, ZHANG Pin, et al. Prediction method of tunneling-induced ground settlement using machine learning algorithms[J]. Journal of Hunan University (Natural Sciences), 2021, 487): 111118. (in Chinese with English abstract)

    [22]

    宫思艺,孔宪光,刘丹,等. 融入复杂地层动态识别的盾构施工地表沉降预测方法研究[J]. 仪器仪表学报,2019,40(6):228 − 236. [GONG Siyi,KONG Xianguang,LIU Dan,et al. An approach for predicting shield construction ground surface settlement of complex stratum using dynamical strata identification[J]. Chinese Journal of Scientific Instrument,2019,40(6):228 − 236. (in Chinese with English abstract)

    GONG Siyi, KONG Xianguang, LIU Dan, et al. An approach for predicting shield construction ground surface settlement of complex stratum using dynamical strata identification[J]. Chinese Journal of Scientific Instrument, 2019, 406): 228236. (in Chinese with English abstract)

    [23]

    林荣安,孙钰丰,戴振华,等. 基于RS-SVR的上软下硬地层盾构施工地表沉降预测[J]. 中国公路学报,2018,31(11):130 − 137. [LIN Rongan,SUN Yufeng,DAI Zhenhua,et al. Predicting for ground surface settlement induced by shield tunneling in upper-soft and lower-hard ground based on RS-SVR[J]. China Journal of Highway and Transport,2018,31(11):130 − 137. (in Chinese with English abstract)

    LIN Rongan, SUN Yufeng, DAI Zhenhua, et al. Predicting for ground surface settlement induced by shield tunneling in upper-soft and lower-hard ground based on RS-SVR[J]. China Journal of Highway and Transport, 2018, 3111): 130137. (in Chinese with English abstract)

    [24]

    李洛宾,龚晓南,甘晓露,等. 基于循环神经网络的盾构隧道引发地面最大沉降预测[J]. 土木工程学报,2020,53(增刊1):13 − 19. [LI Luobin,GONG Xiaonan,GAN Xiaolu,et al. Prediction of maximum ground settlement induced by shield tunneling based on recurrent neural network[J]. China Civil Engineering Journal,2020,53(Sup 1):13 − 19. (in Chinese with English abstract)

    LI Luobin, GONG Xiaonan, GAN Xiaolu, et al. Prediction of maximum ground settlement induced by shield tunneling based on recurrent neural network[J]. China Civil Engineering Journal, 2020, 53(Sup 1): 13 − 19. (in Chinese with English abstract)

    [25]

    LIU Peng,LI Zhengtian,ZHUO Yixin,et al. Design of wind turbine dynamic trip-off risk alarming mechanism for large-scale wind farms[J]. IEEE Transactions on Sustainable Energy,2017,8(4):1668 − 1678. doi: 10.1109/TSTE.2017.2701348

    [26]

    TANG Jiexiong,DENG Chenwei,HUANG Guangbin. Extreme learning machine for multilayer perceptron[J]. IEEE Transactions on Neural Networks and Learning Systems,2016,27(4):809 − 821. doi: 10.1109/TNNLS.2015.2424995

    [27]

    HUANG Guangbin,ZHU Qinyu,SIEW C K. Extreme learning machine:A new learning scheme of feedforward neural networks[C]//Budapest,Hungary:2004 IEEE International Joint Conference on Neural Networks,2005:985 − 990.

    [28]

    HEIDARI A A,MIRJALILI S,FARIS H,et al. Harris Hawks optimization:Algorithm and applications[J]. Future Generation Computer Systems,2019,97:849 − 872. doi: 10.1016/j.future.2019.02.028

    [29]

    BOUAYAD D,EMERIAULT F. Modeling the relationship between ground surface settlements induced by shield tunneling and the operational and geological parameters based on the hybrid PCA/ANFIS method[J]. Tunnelling and Underground Space Technology,2017,68:142 − 152. doi: 10.1016/j.tust.2017.03.011

    [30]

    BOUVEYRON C,LATOUCHE P,MATTEI P A. Exact dimensionality selection for Bayesian PCA[J]. Scandinavian Journal of Statistics,2020,47(1):196 − 211. doi: 10.1111/sjos.12424

    [31]

    苗凤娟,孙同日,陶佰睿,等. PSO与PCA融合优化核极限学习机说话人识别算法仿真[J]. 科学技术与工程,2019,19(21):195 − 199. [MIAO Fengjuan,SUN Tongri,TAO Bairui,et al. Algorithmic research on kernel extreme learning machine for speaker recognition based on PSO and PCA optimization[J]. Science Technology and Engineering,2019,19(21):195 − 199. (in Chinese with English abstract) doi: 10.3969/j.issn.1671-1815.2019.21.029

    MIAO Fengjuan, SUN Tongri, TAO Bairui, et al. Algorithmic research on kernel extreme learning machine for speaker recognition based on PSO and PCA optimization[J]. Science Technology and Engineering, 2019, 1921): 195199. (in Chinese with English abstract) doi: 10.3969/j.issn.1671-1815.2019.21.029

  • 加载中

(5)

(4)

计量
  • 文章访问数:  1414
  • PDF下载数:  118
  • 施引文献:  0
出版历程
收稿日期:  2022-10-08
修回日期:  2023-01-07
刊出日期:  2023-09-15

目录

  • 表 1.  数据集数据统计信息表
    Table 1.  Data set statistics table
    参数 最小值 最大值
    地质参数 地下水位(A1)/m 1.50 1.92
    浮密度(A2)/(g·cm–3 0.54 0.77
    压缩模量(A3)/MPa 3.37 4.42
    泊松比(A4 0.28 0.33
    侧压力系数(A5 0.49 0.59
    几何参数 埋深(A6)/m 18.4 25.3
    盾构参数 土压力(A7)/MPa 1.8 2.7
    总推力(A8)/kN 8 760 16 280
    掘进速度(A9)/(mm·min–1 61 94
    出土量(A10)/ m3 45.8 48.0
    刀盘扭矩(A11)/(N·m) 1 200 2 275
    拼装时间(A12)/h 0.35 0.85
    注浆量(A13)/m3 3.60 4.40
    注浆压力(A14)/bar 3.0 3.8
    地表最大沉降/mm −15.69 10.44
     | Show Table
    DownLoad: CSV
  • 表 2.  主成分系数矩阵
    Table 2.  Principal component coefficient matrix
    影响因子主成分
    12345
    A10.8970.0230.0860.0770.052
    A20.937−0.0770.0430.0410.051
    A3−0.2280.5860.376−0.460−0.163
    A4−0.271−0.259−0.7240.3400.105
    A50.956−0.068−0.130−0.007−0.028
    A60.3760.3350.3650.4920.315
    A70.683−0.093−0.116−0.5200.060
    A8−0.0780.3100.425−0.0530.748
    A9−0.4740.530−0.5820.0190.216
    A10−0.2270.867−0.2470.123−0.128
    A11−0.809−0.112−0.056−0.2560.188
    A12−0.790−0.4640.277−0.0710.080
    A13−0.201−0.6360.0460.0190.247
    A14−0.445−0.0110.6220.384−0.365
     | Show Table
    DownLoad: CSV
  • 表 3.  4类模型预测结果
    Table 3.  Prediction results of four types of models
    组号123456789101112131415
    实测值/mm2.899.865.034.452.0710.440.823.640.4413.429.694.962.541.931.43
    BP预测值/mm1.62−8.493.61−8.92−0.5612.08−0.973.510.12−13.7810.10−3.742.532.240.70
    RBF预测值/mm5.040.492.122.00−2.862.810.314.480.50−11.035.07−11.392.032.010.50
    ELM预测值/mm−0.24−6.174.210.11−0.9412.900.27−0.51−1.51−11.1311.95−4.052.632.120.26
    PCA-HHO-ELM预测值/mm2.11−7.834.93−5.250.149.38−1.907.34−1.23−12.6612.93−4.723.833.46−0.82
    组号161718192021222324252627282930
    实测值/mm0.016.335.2214.1913.741.144.175.311.610.089.901.2410.842.2412.35
    BP预测值/mm−0.064.854.89−14.13−14.250.494.105.931.73−0.05−9.73−0.01−10.461.30−11.83
    RBF预测值/mm0.503.592.56−11.32−7.450.504.475.070.500.50−9.670.50−8.020.58−7.45
    ELM预测值/mm−0.535.725.79−11.76−15.911.342.905.711.65−0.58−9.780.70−8.43−0.12−14.09
    PCA-HHO-ELM预测值/mm−1.336.506.89−13.27−14.072.164.352.434.00−0.99−9.010.13−10.761.18−12.41
     | Show Table
    DownLoad: CSV
  • 表 4.  模型性能评价结果
    Table 4.  Performance evaluation results of the model
    评价指标RMSEMAER2最大误差/mm平均误差/mm
    PCA-HHO-ELM模型0.143 50.026 20.959 63.701.29
    BP预测模型0.232 30.042 40.894 34.480.82
    RBF模型0.392 10.071 60.698 610.352.52
    ELM预测模型0.571 80.104 40.359 24.561.48
     | Show Table
    DownLoad: CSV