img

QQ群聊

img

官方微信

  • CN 62-1112/TF 
  • ISSN 1005-2518 
  • 创刊于1988年
高级检索

黄金科学技术, 2024, 32(1): 109-122 doi: 10.11872/j.issn.1005-2518.2024.01.116

采选技术与矿山管理

基于改进XGBoost算法的深部巷道松动圈智能预测研究

凡兴禹,, 王雪林

核工业井巷建设集团有限公司,浙江 湖州 313000

Research on Intelligent Prediction of EDZ Around Deep Tunnels Based on Improved XGBoost Algorithm

FAN Xingyu,, WANG Xuelin

Nuclear Industry Jingxiang Construction Group Co. ,Ltd. ,Huzhou 313000,Zhejiang,China

收稿日期: 2023-08-11   修回日期: 2023-10-31  

Received: 2023-08-11   Revised: 2023-10-31  

作者简介 About authors

凡兴禹(1995-),男,云南宣威人,高级工程师,从事地下巷道稳定性分析和采矿方法研究Fanxingyu95@gmail.com , E-mail:Fanxingyu95@gmail.com

摘要

深部巷道爆破开挖后由于爆炸冲击和原位应力动态卸载耦合作用,围岩内不可避免地产生松动圈,进而影响结构的稳定性,因此对松动圈厚度进行超前预测显得非常重要。依托多座地下矿山松动圈测试作为研究对象,共获取300组有效数据样本。采用4种主流的超参数优化算法,即遗传算法(GA)、灰狼优化算法(GWO)、粒子群优化算法(PSO)和樽海鞘算法(SSA)对XGBoost算法进行优化,并以此构建4种松动圈预测混合模型。采用R2RMSEMAEMAPE指标对预测模型的性能进行对比分析,并开展松动圈厚度参数的敏感性分析。最后,将最优的PSO-XGBoost模型应用于地下矿山运输巷道进行工程验证。结果表明:在群体规模分别为90、70、60和100时,GA-XGBoost、GWO-XGBoost、PSO-XGBoost和SSA-XGBoost模型取得了最佳的预测表现。其中,PSO-XGBoost模型在训练集和测试集中的相关系数分别为0.9244和0.8787,具有最佳的预测性能。相比基准模型(XGBoost、RF、SVM和LightGBM),优化后模型松动圈的预测精度和性能均得到显著提升。巷道当量直径(TD)和围岩地质强度指标(GSI)对松动圈厚度的影响最为显著,垂直主应力也具有明显的影响。优化后的XGBoost模型在实际工程中的应用结果显示实测值与预测值误差在10%以内,PO-XGBoost具有工程应用价值。

关键词: 松动圈 ; 深部巷道 ; 机器学习 ; 人工智能 ; 地应力 ; 优化XGBoost算法

Abstract

During deep tunnelling using drill-and-blast method,excavation damaged zone (EDZ) is inevitably induced in surrounding rocks due to the coupled impacts of blast loading and dynamic initial stress unloading and thus affect the structure stability.Therefore,it is very important to predict EDZ depth before roadways excavation.Relying on the field measurements of EDZ in several underground mines as the research object,300 data samples were collected.Four mainstream hyperparametric optimization algorithms,i.e.,genetic algorithm (GA),gray wolf optimization algorithm(GWO),particle swarm optimization algorithm(PSO),and salp swarm algorithm (SSA),were used to optimize the XGBoost algorithm and to construct four hybrid models for EDZ prediction.Comparative analysis of predictive model performance was conducted in terms of R2RMSEMAE and MAPE,along with a sensitivity analysis of the influencing parameters.Finally,the optimal PSO-XGBoost model was applied to a transportation roadway in an underground mine for engineering validation.The results show that the GA-XGBoost,GWO-XGBoost,PSO-XGBoost,and SSA-XGBoost models achieve the best predictive performance with swarm sizes of 90,70,60 and 100,respectively.Among them,the PSO-XGBoost model demonstrates the best predictive performance with correlation coefficients of 0.9244 and 0.8787 in the training and testing sets,respectively.Moreover,compared to bench models(XGBoost,RF,SVM and LightGBM),both the prediction accuracy and stability of the optimized models are improved.The tunnel diameter(TD) and rock mass geological strength index(GSI) have the most significant influence on the loosened zone thickness,along with a noticeable impact from the vertical principal stress.The application results of the optimized XGBoost model in practical engineering show that the error between the measured value and the predicted value is within 10%,indicating that the PO-XGBoost is of significance for engineering application.

Keywords: excavation damaged zone(EDZ) ; deep tunnels ; machine learning ; artificial intelligence ; in-situ stress ; optimized XGBoost algorithm

PDF (7813KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

凡兴禹, 王雪林. 基于改进XGBoost算法的深部巷道松动圈智能预测研究[J]. 黄金科学技术, 2024, 32(1): 109-122 doi:10.11872/j.issn.1005-2518.2024.01.116

FAN Xingyu, WANG Xuelin. Research on Intelligent Prediction of EDZ Around Deep Tunnels Based on Improved XGBoost Algorithm[J]. Gold Science and Technology, 2024, 32(1): 109-122 doi:10.11872/j.issn.1005-2518.2024.01.116

进入21世纪以来,随着社会经济的高速发展,人类对资源的需求量持续增加。然而,浅部矿产资源日益枯竭,迫使采矿工作转向地下深部。与此同时,为了推动经济发展并克服交通限制,全国各地积极推动深埋隧道的建设,这些深部巷道已成为获取地下资源和连接地区发展的至关重要的通道(谢和平等,2015王勇等,2023余涛等,2023唐登志等,2023)。钻爆法施工具有灵活性强、成本低和开挖能力强等优势,仍是当前最主要的硬岩开挖技术。在使用钻爆法开挖深部地下巷道时,爆破荷载和原岩应力卸载的耦合作用会导致巷道围岩内出现开挖损伤区,即松动圈。松动圈的厚度不仅是评估巷道稳定性的重要指标,而且会影响巷道的施工和后续支护工序。

松动圈厚度的预测方法主要包括经验公式(Perras et al.,2016刘刚等,2021)、原位测量(刘蒙等,2022王旭等,2023)、数值模拟(Sun et al.,2021)和机器学习预测(胡军等,2014朱志洁等,2014)。经验公式法是基于现场测试和弹塑性力学分析结果提出的,如:Wang et al.(2019)修正了Hoek-Brown准则并推导出松动圈的理论解析解;Qiao et al.(2021)基于Drucker-Prager屈服准则推导出松动圈厚度的解析解。原位测量方法主要采用声波测试法和钻孔电视全景摄像法,如:Li et al.(2012)成功地应用数字全景钻孔摄像技术获取了锦屏二号水电站地下隧道的围岩松动圈厚度和分布特征;Wang et al.(2015)在国内地下煤矿中同时采用声波测试和数字全景钻孔摄像技术对地下巷道的松动圈进行测量,并根据测量结果重新优化了支护设计,提高了巷道的长期稳定性。数值模拟方法主要使用有限元软件LS-DYNA和离散元软件FLAC3D,如:Yang et al.(2017)采用LS-DYNA模拟了初始应力条件下圆形巷道在爆破荷载作用下松动圈厚度的演化特征;Hong et al.(2023)基于LS-DYNA研究了松动圈厚度与地应力大小、侧压力系数、巷道形状、装药系数与围岩强度之间的关系。理论分析和数值模拟方法均基于一定的假设条件,其结果可靠性难以保证。现场测量方法虽然可靠性较高,但需要耗费大量人力和成本,且周期较长。

相比之下,机器学习技术在解决非线性问题方面具有巨大优势,已被广泛应用于解决工程问题,包括岩爆预测、地下硐室稳定性分析、边坡和尾矿库稳定性分析等。许多学者也开始使用机器学习方法进行松动圈厚度的预测研究。如:Zhou et al.(2011)基于SVM(支持向量机)算法构建了松动圈厚度预测模型,相比回归模型,该模型具有更高的可靠性和精度;Xie et al.(2019)采用随机森林回归模型预测松动圈厚度,该方法也具有较高的可靠性。此外,BPNN(BP神经网络)、SVR(支持向量回归)和GP(高斯回归)算法也被广泛应用于松动圈厚度的预测研究。现有研究表明,机器学习方法在松动圈预测方面具有成本低、耗时少和预测结果可靠等优势。然而,几乎所有的预测模型均忽略了地应力对松动圈厚度的影响,且考虑的影响参数较少。

XGBoost算法是目前主流的机器学习算法,具有高性能和可扩展性。该算法内部集成了正则化技术,支持交叉验证,从而提高了模型的泛化能力并降低了过拟合风险。本文首先基于XGBoost模型构建了松动圈厚度预测模型,并采用GA、GWO、PSO和SSA对模型进行优化,从而进一步提高松动圈厚度的预测精度和可靠性。其次,依托地下矿山爆破开挖工程获取300组有效的松动圈数据样本。基于R2RMSEMAEMAPE模型性能指标对4种混合模型和4种基准模型进行比较分析,并使用XGBoost算法进行松动圈影响参数敏感性分析。最后,将PSO-XGBoost模型用于武山铜矿运输巷道围岩松动圈预测中,从而验证模型在实际工程中的可靠性。

1 松动圈数据

松动圈常伴随地下巷道开挖完成后产生,其产生的主要原因是炸药爆破产生的冲击荷载和岩体开挖后重分布的应力场的共同作用(图1)。使用SETPLT-02型声波探测仪在玲珑金矿、凡口铅锌矿和武山铜矿开展松动圈原位测试,并对测试地点进行地质调查和钻芯取样后开展系列岩石力学试验。试验矿山地理位置及现场测试过程如图2所示。凡口铅锌矿位于广东省韶关市仁化县境内,是我国大型铅锌矿生产基地。松动圈测试点位于-680 m分段回风巷道内,且所处水平地应力较高,巷道围岩为中等强度灰质砂岩,爆破开挖后监测到巷道出现明显的变形和破坏。

图1

图1   深部巷道松动圈示意图

Fig.1   Schematic diagram of EDZ around deep tunnels


图2

图2   开展松动圈测试的地下矿山

Fig.2   Underground mines conducting EDZ measurements


参考已有研究,基于改进XGBoost模型的智能预测包含地质强度指标、岩体单轴抗压强度、垂直主应力和装药系数等在内的8个松动圈厚度影响因素,输出参数为松动圈厚度。松动圈影响参数和厚度数据如表1所示。

表1   松动圈输入和输出参数

Table 1  Input and output parameters of EDZ

参数类型参数名称符号单位取值范围均值标准偏差P
输入参数地质强度指标GSI-26.43~77.6953.6613.90.049
单轴抗压强度UCSMPa64.91~155.89101.1723.630.000
垂直主应力σvMPa3.75~40.0819.3510.040.000
侧压力系数λ-0.42~2.311.500.280.000
巷道直径TDm2.80~8.354.731.010.015
装药系数PFkg/m31.12~2.001.580.180.000
抵抗线Bm0.50~1.050.760.100.031
孔间距HSm0.66~1.050.860.090.001
输出参数松动圈厚度EDZcm6.7~65.241.413.20.802

新窗口打开| 下载CSV


为了对数据进行可视化分析,将8个输入参数和1个输出参数绘制成小提琴形式,如图3所示。本次样本数据中松动圈厚度集中在30~60 cm之间,整体分布较为均匀。GSI值小提琴图较宽,表明参数值在整个范围内分布很均匀。其他输入参数的小提琴图也较宽,说明数据的分布也很均匀,可靠性较高。需要注意的是,测试地点巷道侧压力系数小于0.5时频次较低,而主要集中在1.5附近。

图3

图3   松动圈输入和输出参数分布

Fig.3   Distribution of EDZ input and output parameters


2 XGBoost及优化算法

2.1 XGBoost算法

XGBoost最初由Chen et al.(2016)提出,是由梯度提升决策树(GBDT)的概念发展而来的一个高度可扩展的机器学习系统。与GBDT相比,XGBoost模型显示出以下优势:(1)有利于惩罚模型的复杂性,避免过拟合;(2)能够以较高的计算速度并行处理海量数据;(3)允许用户定义目标函数,只要其二阶可微;(4)有效处理涉及缺失值的数据集。因此,XGBoost算法已被广泛应用于许多领域,具有优异的性能。

XGBoost算法的计算流程如图4所示。XGBoost中的算法复杂性作为目标函数的一部分被优化,该目标函数的表达式如下:

Obj=i=1nl(yi,y^i)+k=1KΩ(fk)

式中:l为损失函数;Ω为正则项函数。

图4

图4   XGBoost算法智能预测结构示意图

Fig.4   Schematic diagram of intelligent prediction structure by XGBoost algorithm


由于y^i=k=1Kfk(xi)=y^i(t-1)+ft(xi),第t次迭代的目标函数可以用一个更具体的数学模型来描述,表示为

Obj=lyi,y^i(t-1)+ft(xi)+k=1KΩ(fk)

然后,通过3个步骤对目标函数进行优化,目标函数的最终形式如式(3):

Obj=j=1TGjωj+12(Hj+λ)ωj2+γT

式中:Gj=iIjgiHj=iIjhi。总的来说,XGBoost算法通过建立一个变量的二次方程,将寻求最佳目标值转化为一个关于寻找最小值的问题。

2.2 GA算法

GA(遗传算法)是由Holland(1992)首次提出并发展的。该算法的工作原理是维持一个候选解决方案的群体,被称为个体或染色体,代表手头问题的潜在解决方案。群遗传算法模仿自然选择的过程,在几代的时间里反复演化这个群体,它由初始化、评估、选择、繁殖和替换5个关键步骤组成(Yang,2021),其获取最优值的流程如图5所示。

图5

图5   GA算法流程

Fig.5   Flow chart of GA algorithm


通过连续几代,遗传算法倾向于提高种群的整体适配性,向更好的问题解决方案靠拢。通过探索不同的基因组合并对其进行逐渐完善,该算法能够为广泛的复杂问题找到最佳或接近最佳的解决方案。

2.3 GWO算法

GWO(灰狼优化算法)是Mirjalili et al.(2014)提出的一种基于蜂群的新型算法。GWO形成了一个严格的等级结构,其灵感来自于狼的狩猎和社会行为,不同的狼在灰狼群中扮演着不同的角色。每一匹狼都被严酷的社会秩序所支配。GWO算法流程如图6所示。

图6

图6   GWO算法流程

Fig.6   Flow chart of GWO algorithm


2.4 PSO算法

PSO(粒子群优化算法)是众多进化搜索方法中用于解决优化问题的一种强有力的方法,通过模拟鱼群和鸟群来实现(Kennedy,1995)。PSO算法流程如图7所示,在PSO算法中,鸟群代表一组粒子,而食物源代表一个功能目的。通过分享和传递鸟群与食物源之间的距离信息,食物源的位置可以由鸟群来确定。这种合作使整个鸟群能够选择有关食物源位置的最佳信息,并最终聚集在食物周围。在计算过程中,PSO算法中的一些重要参数被更新,通过这些参数可以确定新的位置和速度。

图7

图7   PSO算法流程

Fig.7   Flow chart of PSO algorithm


2.5 SSA算法

与上述2种算法一样,SSA(樽海鞘群算法)属于一组基于蜂群智能的优化算法(Cortés-Caicedo et al.,2022),其灵感来自于樽海鞘群的狩猎和导航行为。在SSA算法中,主要有3个步骤来优化超参数,即樽海鞘群的初始化、领导者的更新和跟随者的更新。在初始化过程中,需要定义樽海鞘群的数量和位置。此外,还需要根据用户要求确定最佳食物来源(适应度函数)。在下一步,优化过程在搜索域中开始,这个搜索域可用一个名为 Si 的矩阵表示:

Si=s11sn1s1msnm

在这个搜索空间中,领导者会首先行动起来,为跟随者提供指引。领导者樽海鞘位置的更新方程为

uji=fj+RCgbj-hbjRB+hbj   (C0.5)fj-RCgbj-hbjRB+hbj   (C<0.5)

式中:gbjhbj 分别为j维度上的搜索上界和下界;fjj维度上的食物位置;RBRC为2个随机参数,数值为[0, 1]。参数RB决定了樽海鞘运动的长度,可以控制运动的方向。参数RC在控制探索和开发方面起着重要作用。

重复上述步骤直至优化达到停止标准。SSA算法的优化过程如图8所示。

图8

图8   SSA算法流程

Fig.8   Flow chart of SSA algorithm


2.6 基于XGBoost的混合松动圈预测模型

选择4种优化算法来调整 XGBoost的超参数,以构建松动圈厚度智能预测混合模型,即GA-XGBoost、GWO-XGBoost、PSO-XGBoost和SSA-XGBoost。虽然4种优化算法的优化思路大不相同,但相应的4种模型的整体架构非常相似,可以概括为以下5个步骤:(1)实验数据预处理;(2)模型参数初始化;(3)种群初始化;(4)迭代循环;(5)输出最优解。以GWO-XGBoost模型为例,松动圈厚度预测流程如图9所示。

图9

图9   GWO-XGBoost模型流程

Fig.9   Flow chart of the GWO-XGBoost model


2.7 模型性能评价指标

所构建的预测模型使用性能指标(PI)进行评估,PI包括均方根误差(RMSE)、平均绝对误差(MAE)、平均绝对百分比误差(MAPE)和标准偏差(R2)。假设RMSE = 0,MAE = 0,MAPE = 0,R2 = 1,则认定预测模型具有十分可靠的预测能力。4种模型性能评价指标的表达式如下:

RMSE=1ni=1n(yi-y^i)2
MAE=1ni=1nyi-y^i
MAPE=100%ni=1nyi-y^iyi
R2=1-i(yi-y^i)2i(yi-y¯i)2

式中:yiy^iy¯i分别为第i个样本的真实值、预测值和真实值的平均值;n为样本数量。

3 结果及分析

为了构建基于XGBoost的深部巷道爆破围岩松动圈厚度预测模型,采用4种元启发式算法对XGBoost进行优化。将地质条件、爆破参数和开挖特征等8个影响因素作为输入,以松动圈厚度作为输出。其中,蜂群数量和迭代次数在优化性能中起着关键作用。因此,对不同蜂群数量下的混合模型进行均方误差(MSE)比较和分析。

3.1 GWO-XGBoost模型

一般来说,随着迭代次数的增加,优化性能趋于稳定,但计算时间也会增加。为了判断优化性能,使用MSE值进行分析。在测试过程中,设置预测性能稳定的迭代次数为500次,对60、70、80、90和100的蜂群规模(粒子数)进行测试,以选择最佳的最优参数。如图10所示,每个蜂群数量(粒子数)的变化过程是不同的,随着迭代次数的增加,MSE值不断下降。每个模型在达到优化终点时,均获得了最低的MSE值。从图10可以看出,当蜂群数量(粒子数)为100时,其MSE最优值大于其他蜂群数量(粒子数)的MSE最优值,且最优值对应的迭代次数远大于其他蜂群数量(粒子数)的迭代次数,表明蜂群数量(粒子数)为100时模型效率低且稳定性差。相反,当蜂群数量(粒子数)为90时,MSE值在最小的迭代次数(121次)时到达稳定值,相比其他蜂群数量的MSE值更小,所以认为蜂群数量(粒子数)为90时GWO-XGBoost模型具有最佳的预测性能。

图10

图10   GWO-XGBoost模型不同粒子数下MSE值随迭代次数变化

Fig.10   Variations of MSE value with iteration number of GWO-XGBoost model under different swarm sizes


为了进一步验证GWO-XGBoost模型在不同蜂群数量下的预测性能,计算训练集和测试集中的MAW、MAPE、RMSER2,结果如表2所示。从表2可以看出,当粒子数为90时,GWO-XGBoost模型获得最佳的预测性能,训练集中MAW、MAPE、RMSER2分别为6.812、6.943、9.522和0.9073,具有较好预测性能的粒子数为60,相关系数在0.9以上。当粒子数为100时,GWO-XGBoost模型的预测性能最差,这与图10中的结论一致,此时,训练集中MAWMAPERMSER2分别为8.212、8.994、11.705和0.8723。综上认为,当粒子数为90时,GWO-XGBoost模型无论是在训练集还是测试集中均具有最好的拟合能力和最佳的预测性能。

表2   GWO-XGBoost模型在不同粒子数下的性能

Table 2  Performance of GWO-XGBoost model under different swarm sizes

粒子数训练集测试集
MAEMAPERMSER2MAEMAPERMSER2
607.0317.32110.1430.90147.5937.90710.9540.8560
707.7438.33111.0430.88028.3628.99711.9260.8361
807.2027.95410.6310.88927.7788.59011.4810.8453
906.8126.9439.5220.90737.3577.49810.2840.8621
1008.2128.99411.7050.87238.8699.71412.6410.8290

新窗口打开| 下载CSV


3.2 GA-XGBoost模型

GA算法可以提升XGBoost模型的参数选择。在使用GA-XGBoost进行松动圈厚度预测之前,选择最有效的GA参数。与GWO-XGBoost模型构建一样,本节将讨论迭代次数和蜂群大小对预测性能的影响。同样设定蜂群数量(粒子数)为60、70、80、90和100,不同蜂群值对应的MSE值和迭代次数曲线如图11所示。

图11

图11   GA-XGBoost模型不同粒子数下MSE值随迭代次数变化

Fig.11   Variations of MSE value with iteration number of GA-XGBoost model under different swarm sizes


图11可以看出,当粒子数为70时,GA-XGBoost模型获得最小的MSE值,且MSE值到达稳定值时对应的迭代次数也较小。如表3所示,当粒子数为70时,训练集中MAW=7.483、MAPE=7.201、RMSE=10.423和R2=0.8734,测试集中MAW=8.802、MAPE=7.777、RMSE=11.257,R2=0.8384。然而,当粒子数为90时,GA-XGBoost模型的综合预测性能最差,这是因为Fitness曲线收敛时对应的MSE值最大,且208次的迭代次数也是最大的,说明该蜂群数量下的模型稳定性和可靠性均处于一个最差的状态。当粒子数为90时,训练集中MAW=8.312、MAPE=8.834、RMSE=12.795和R2=0.8432,测试集中MAW=8.977、MAPE=9.433、RMSE=13.819和R2=0.8098。

表3   GA-XGBoost模型在不同粒子数下的预测性能

Table 3  Performance of GA-XGBoost model under different swarm sizes

粒子数训练集测试集
MAEMAPERMSER2MAEMAPERMSER2
607.5728.13512.1430.85738.1788.78613.1140.8233
707.4837.20110.4230.87348.0827.77711.2570.8384
807.3847.75911.5430.87027.9758.38012.4660.8357
908.3128.73412.7950.84328.9779.43313.8190.8098
1007.9438.42912.4810.85288.5789.10313.4790.8191

新窗口打开| 下载CSV


3.3 PSO-XGBoost模型

图12所示为PSO-XGBoost模型在不同群体规模(粒子数)下的MSE值随迭代次数变化的曲线。从图12可以看出,无论群体规模较小(粒子数为60或70)还是较大(粒子数为100)时,MSE值在初始阶段均迅速下降,然后逐渐趋于稳定。但不难发现,不同群体规模对应的模型收敛MSE值不同,且相应的迭代次数也不同。准确来说,当粒子数为60和70时,MSE值最小,二者之间的差别极小,但是当粒子数为60时模型收敛速度要比粒子数为70时模型收敛速度快。当粒子数为90时,MSE值最大,表明模型稳定性最差,当粒子数为100时,收敛MSE值对应的迭代次数最大,表明模型的收敛速度最慢,性能不佳。

图12

图12   PSO-XGBoost模型不同粒子数下MSE值随迭代次数变化

Fig.12   Variations of MSE value with Iteration number of PSO-XGBoost model under different swarm sizes


为了更加全面地对比不同群体规模下,模型在预测松动圈厚度时的性能,计算各自在训练集和测试集中的性能指标值,结果列于表4。从表4可以看出,当群体规模(粒子数)为60时,PSO-XGBoost模型呈现出最佳的预测性能,此时训练集中MAW=5.231、MAPE=5.513、RMSE=7.143和R2=0.9244,测试集中MAW=5.543、MAPE=5.844、RMSE=77.752和R2=0.8787,各项指标均比粒子数为70时更佳。当粒子数为90时,在训练集和测试集中不仅相关系数是最小的,且MAE、MAPERMSE均最大,表明该群体规模下的模型预测性能最差。

表4   PSO-XGBoost模型在不同粒子数下的预测性能

Table 4  Performance of PSO-XGBoost model under different swarm sizes

粒子数训练集测试集
MAEMAPE/%RMSER2MAEMAPE/%RMSER2
605.2315.5137.1430.92445.5435.8447.5720.8787
705.7415.7347.4550.92126.0816.0787.9020.8752
805.8285.9457.5430.91896.1786.3027.9960.8731
906.3216.7129.3650.90236.7127.1159.9270.8573
1005.9346.2948.0110.91236.2986.6728.4920.8670

新窗口打开| 下载CSV


3.4 SSA-XGBoost模型

作为一种新型的群集智能算法,SSA在各种最优问题上显示出其优越性和可行性。SSA很容易实现和调整。为了便于比较,在基于SSA的模型中也采用相同的超参数值。图13所示为不同蜂群规模下SSA-XGBoost模型的优化过程。由图13可以看出,不同蜂群规模下SSA-XGBoost模型随迭代次数的增加均会逐渐趋于稳定,意味着模型在迭代过程中会逐渐接近最优解,但是收敛时的均方误差大小并不相等,且收敛速度也不一样。当粒子数为100时模型收敛时的MSE值最小,且收敛速度最快,对应的迭代次数为95次。其次收敛性能表现最好对应的蜂群规模(粒子数)为80。很明显,当粒子数为60时,虽然模型收敛速度较快,但MSE值最大,说明模型在该蜂群规模下的稳定性较差。

图13

图13   SSA-XGBoost模型不同粒子数下MSE值随迭代次数变化

Fig.13   Variations of MSE value with iteration number of SSA-XGBoost model under different swarm sizes


表5为不同群体规模下模型的评价指标值。当粒子数为100时,无论是在训练集还是测试集,相关系数均比其他群体规模更大,而MAE、MAPERMSE均更小,表明该模型具有最佳的综合预测性能。与此同时,表5中的数据再次证明当粒子数为60时SSA-XGBoost模型的表现最差,此时在训练集中MAW=8.213、MAPE=8.733、RMSE=12.46、R2=0.8629,测试集中各项评价指标分别为MAW=8.624、MAPE=9.170、RMSE=13.089、R2=0.8208。

表5   SSA-XGBoost模型在不同粒子数下的预测性能

Table 5  Performance of SSA-XGBoost model under different swarm sizes

粒子数训练集测试集
MAEMAPERMSER2MAEMAPERMSER2
608.2138.73312.4660.86298.6249.17013.0890.8208
707.4267.98610.9620.88117.7978.38511.5010.8374
807.1327.32910.5430.88567.4897.69511.0700.8416
908.0128.34511.6070.87438.4138.76212.1870.8313
1007.0877.0129.9140.89337.4417.36310.4100.8496

新窗口打开| 下载CSV


3.5 模型对比分析

基于上述讨论可知,GWO-XGBoost模型、GA-XGBoost模型在、PSO-XGBoost模型和SSA-XGBoost模型在粒子数分别为90,70,60,100时具有最佳的预测性能,因此选择这4种群体规模下的4种模型对60组测试数据进行预测,同时采用未进行优化的XGBoost模型和3种常见的SVM、RF和LightGBM模型对同组数据进行预测,对比分析这8种模型的优劣。图14所示为各种预测模型预测值与实测值的二维分布和线性拟合图。由图14可以看出,PSO-XGBoost模型中的数据点集中在线性拟合曲线附近,95%预测带宽相对其他模型更窄。其他7种松动圈预测模型的数据点离散型均比PSO-XGBoost模型大,尤其是LightGBM和SVM模型,说明PSO-XGBoost模型具有比其他模型更精确和更稳定的预测性能,用于深部巷道爆破开挖后围岩松动圈厚度的预测是最合适的。其次,对比混合模型(GWO-XGBoost、GA-XGBoost、PSO-XGBoost和SSA- XGBoost)与基准模型(XGBoost、SVM、RF、Light-GBM)可知,优化后的模型具有更好的预测性能。

图14

图14   不同模型预测值与真实值对比

Fig.14   Comparison of predicted values and actual values of different models


为了进一步量化分析8种模型的性能,将训练集和测试集的4种性能指标值(MAE、MAPE、RMSER2)绘制成雷达图(图15)。图15(a)所示为不同模型MAE值的分布,很明显PSO-XGBoost模型训练集和测试集的MAE值最小,其次是GWO-XGBoost模型,LightGBM模型具有最大的MAE值,表明松动圈厚度智能预测性能最佳的是PSO-XG-Boost模型,而最差的是LightGBM模型。图15(b)所示为各模型的MAPE值,其分布特征与MAE值一样。PSO-XGBoost模型在训练集和测试集中的MAPE值分别为5.513和5.844,其次是GWO-XGBoost模型,其在下一步训练集和测试集中的MAPE值分别为6.943和7.498。基准模型得到的MAPE值明显比优化模型大,LightGBM模型的MAPE值分别为11.310和12.215,再次证明PSO-XGBoost模型的预测性能。图15(b)和图15(c)所示为各模型的RMSER2。对于RMSE,LightGBM和SVM模型得到的值最大,而PSO-XGBoost和GWO-XGBoost模型得到的值最小。对于R2,情况正好相反,PSO-XGBoost和GWO-XGBoost模型的相关系数是最大的,在训练集中分别为0.9244和0.9073,而LightGBM和SVM模型在训练集中的相关系数分别为0.8243和0.8376。

图15

图15   不同模型的性能评价指标比较

Fig.15   Comparison of performance evaluation indexes of different models


综上所述,无论是在训练集还是测试集中,PSO-XGBoost模型在预测精度、泛化能力和稳定性方面均表现最好,其次是GWO-XGBoost模型。未经超参数调整的基准模型在预测过程中稳定性较差且预测精度不高。因此,在智能预测采用钻爆法开挖的地下巷道的围岩松动圈厚度时使用PSO-XGBoost和GWO-XGBoost模型是可靠的。

3.6 参数敏感性分析

参数敏感性分析是一种检查系统稳定性的方法。在本研究中,网络创建后,对数据集进行敏感性分析。敏感性分析检查与每个模型输入相关的系统倾向和变化。这是通过在期望的范围内改变每个输入的值来完成的。敏感性分析的通用方程式如下:

sk*=sk(ak*)=dfk(ak)dakak=ak*ak*p*k=1,2,,n

式中:sk*k为无量纲的实数组,为正数。大于sk*的值表示相对于ak,指标p的灵敏度更高。通过比较不同的sk*值,可以求得每个输入参数的敏感度强度。

使用XGBoost模型对所有输入参数逐个进行敏感性检测,在某参数进行检测的过程中将其他参数固定。得到各松动圈厚度影响因素的敏感性值并以直方图形式表示,如图16所示。从图16可以看出,深部巷道开挖直径(TD)对超挖厚度影响最大,敏感性参数为0.73;其次是开挖巷道周边围岩的地质强度指标,敏感性参数为0.68。研究表明,通过调整开挖尺寸和强化岩体质量,可以有效减小深部巷道爆破开挖后的松动圈厚度。值得注意的是,作用于开挖巷道上的垂直主应力对松动圈厚度也具有显著影响,敏感性参数为0.56。对于采用光面爆破技术进行开挖的巷道而言,周边孔装药系数和孔间距对松动圈厚度也具有一定强度的影响,因此可以通过优化爆破参数控制围岩损伤区的分布。

图16

图16   松动圈厚度影响因素敏感性分析

Fig.16   Sensitivity analysis of EDZ depth influencing factors


4 工程验证

为了进一步验证PSO-XGBoost松动圈厚度预测模型的可靠性,将该模型应用于武山铜矿地下运输巷道松动圈超前预测中。武山铜矿隶属于江西铜业股份有限公司,位于江西省九江市西瑞昌县白羊镇。矿区包含南北2个矿带,矿床赋存条件差异较大,且地质赋存环境复杂。测试巷道围岩以白色大理岩为主,在整个掘进线路上伴随有蚀变大理岩和砂岩。为了尽量降低爆破开挖过程中爆炸冲击对围岩的损伤,采用光面爆破技术进行掘进。巷道掘进方向与水平最大主应力方向平行,可最大开挖时原位应力卸载对岩体造成进一步损伤。巷道设计宽度为3.4 m,高度为3.0 m。现场使用乳化炸药进行岩体破碎,光面爆破参数根据围岩特性进行调整。

在巷道沿线上选择6个截面进行松动圈厚度预测可靠性验证,不同截面处的岩体特性不同,通过开展室内岩石力学试验和现场节理裂隙调查获得。巷道掘进前,依据岩体参数、爆破参数和巷道开挖直径采用PSO-XGBoost预测模型对松动圈的厚度进行超前预测,待掘进完成后对测试点进行松动圈现场测量,测量结果与预测值列于表6中。由表6可知,6个测试点岩体地质强度指标(GSI)均不相等,完整岩样的单轴抗压强度为87.9~145.8 MPa,软岩的装药系数较硬岩更小,抵抗线和孔间距则更大,说明测试点的选择具有代表性。松动圈厚度测试结果为13.43~35.12 cm,预测值为13.87~38.01 cm,预测误差最大值为9.5%,有效控制在10%以内,表明经过优化后的PSO-XGBoost模型在工程现场具有良好的表现。

表6   松动圈预测值与实测值比较

Table 6  Comparison of EDZ predicted value and mesured value

测试点输入参数输出参数误差/%
GSIUCS/MPaσv/MPaλTD/mPF/(kg·m-3B/mHS/m实测值/cm预测值/cm
Test-171137.412.41.353.11.560.650.5515.3213.879.5
Test-272145.812.41.353.11.500.650.5513.4314.528.1
Test-352101.512.41.353.11.320.800.6722.5623.835.6
Test-45194.812.41.353.11.350.800.6723.7522.266.3
Test-53687.912.41.353.11.181.000.8535.1238.018.2
Test-64090.112.41.353.11.201.000.8532.5434.646.5

新窗口打开| 下载CSV


因此,可以在地下巷道或硐室掘进之前对松动圈厚度进行超前预测,一方面可以根据预测结果对爆破参数进行优化,另一方面也可以根据预测结果制定开挖后的围岩支护方案。

5 结论

(1)群体规模和迭代次数对预测模型的性能具有显著影响。GA-XGBoost、GWO-XGBoost、PSO-XGBoost和SSA-XGBoost模型分别在群体规模(粒子数)为90、70、60和100时取得最佳预测表现。

(2)PSO-XGBoost是4种混合模型中预测性能最佳的模型,其训练集和测试集中的相关系数分别为0.9244和0.8787;优化后的4种模型比XGBoost、RF、SVM和LightGBM基准模型在松动圈厚度预测中具有更好的性能,同时采用GA、GWO、PSO和SSA优化算法对XGBoost优化后能明显提高预测精度和稳定性。

(3)参数敏感性分析结果表明,各输入变量的相对重要性由高到低分别为巷道直径、围岩地质强度指标、单轴抗压强度、垂直主应力、光面孔孔间距、装药系数、光爆层抵抗线和侧压力系数,因此在进行深部巷道围岩松动圈预测和分析时必须要考虑地应力的影响。

http://www.goldsci.ac.cn/article/2024/1005-2518/1005-2518-2024-32-1-109.shtml

参考文献

Chen T QGuestrin C2016.

XGBoost:A Scalable Tree Boosting System

[C]//Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining.San FranciscoAssociation for Computing Machinery: 785-794.

[本文引用: 1]

Cortés-Caicedo BGrisales-Noreña L FMontoya O D2022.

Optimal selection of conductor sizes in three-phase asymmetric distribution networks considering optimal phase-balancing:An application of the salp swarm algorithm

[J].Mathematics,1018):3327.

[本文引用: 1]

Holland J H1992.Adaptation in Natural and Artificial Systems:An Introductory Analysis with Applications to Biology,Control,and Artificial Intelligence[M].CambridgeMIT Press.

[本文引用: 1]

Hong Z XTao MWu C Qet al2023.

The spatial distribution of excavation damaged zone around underground roadways during blasting excavation

[J].Bulletin of Engineering Geology and the Environment,824):155.

[本文引用: 1]

Hu JunWang KaikaiXia Zhiguo2014.

Support vector machine(SVM)prediction of roadway surrounding rock loose circle thickness optimized by layered fish

[J]. Metal Mine,4311):31-34.

Kennedy JEberhart R1995.

Particle swarm optimization

[C]//Proceedings of the IEEE International Conference on Neural Networks.PerthIEEE.

[本文引用: 1]

Li S JFeng X TLi Zet al2012.

Evolution of fractures in the excavation damaged zone of a deeply buried tunnel during TBM construction

[J].International Journal of Rock Mechanics and Mining Sciences,55125-138.

[本文引用: 1]

Liu GangXiao YongzhuoZhu Junfuet al2021.

Overview on theoretical calculation method of broken rock zone

[J].Journal of China Coal Society,461):46-56.

Liu MengTang HaiMa Yujieet al2022.

Measurement of loose zone of tunnel surrounding rock based on apparent resistivity method

[J].Mineral Engineering Research,374):58-64.

Mirjalili SMirjalili S MLewis A2014.

Grey wolf optimizer

[J]. Advances in Engineering Software,6946-61.

[本文引用: 1]

Perras M ADiederichs M S2016.

Predicting excavation damage zone depths in brittle rocks

[J].Journal of Rock Mechanics and Geotechnical Engineering,81):60-74.

[本文引用: 1]

Qiao SCai ZXu Pet al2021.

Investigation on the scope and influence factors of surrounding rock loose circle of shallow tunnel under bias pressure:A case study

[J]. Arabian Journal of Geosciences,1415):1428.

[本文引用: 1]

Sun QMa FGuo Jet al2021.

Excavation-induced deformation and damage evolution of deep tunnels based on a realistic stress path

[J].Computers and Geotechnics,129103843.

[本文引用: 1]

Tang DengzhiBai GenmingChen Shuanget al2023.

Research on influencing factors and prediction of rockburst in deep-buried tunnel

[J].Sichuan Hydropowder,422):11-17.

Wang H WJiang YXue Set al2015.

Assessment of excavation damaged zone around roadways under dynamic pressure induced by an active mining process

[J].International Journal of Rock Mechanics and Mining Sciences,77265-277.

[本文引用: 1]

Wang RDeng XMeng Yet al2019.

Case study of modified H-B strength criterion in discrimination of surrounding rock loose circle

[J]. KSCE Journal of Civil Engineering,233):1395-1406.

[本文引用: 1]

Wang XuLi Xiaomeng2023.

Research on the present situation and development trend of loose zone test technology for roadway surrounding rock

[J].Coal and Chemical Industry,461):4-7.

Wang YongWu AixiangYang Junet al2023.

Progress and prospective of the mining key technology for deep metal mines

[J].Chinese Journal of Engineering,458):1281-1292.

Xie HepingGao FengJu Yang2015.

Research and development of rock mechanics in deep ground engineering

[J].Chinese Journal of Rock Mechanics and Engineering,3411):2161-2178.

Xie QPeng K2019.

Space-time distribution laws of tunnel excavation damaged zones (EDZs) in deep mines and EDZ prediction modeling by random forest regression

[J].Advances in Civil Engineering,(7):6505984.

[本文引用: 1]

Yang J HYao CJiang Qet al2017.

2D numerical analysis of rock damage induced by dynamic in-situ stress redistribution and blast loading in underground blasting excavation

[J].Tunnelling and Underground Space Technology,70221-232.

[本文引用: 1]

Yang X S2021.Genetic Algorithms[M].New YorkAcademic Press.

[本文引用: 1]

Yu TaoZhu NingboYao Zhiganget al2023.

Stability analysis and control of tunnel excavation in deep buried horizontal stratum

[J].Journal of Safety Science and Technology,194):93-99.

Zhou JLi X B2011.

Evaluating the thickness of broken rock zone for deep roadways using nonlinear SVMs and multiple linear regression model

[J].Procedia Engineering,26972-981.

[本文引用: 1]

Zhu ZhijieZhang HongweiChen Ying2014.

Prediction model of loosening zones around roadway based on MPSO-SVM

[J].Computer Engineering and Applications,5012):1-5.

胡军王凯凯夏治国2014.

分层鱼群优化支持向量机预测巷道围岩松动圈厚度

[J].金属矿山,4311):31-34.

[本文引用: 1]

刘刚肖勇卓朱俊福2021.

围岩松动圈理论计算方法的评述与展望

[J].煤炭学报,461):46-56.

[本文引用: 1]

刘蒙唐海马谕杰2022.

基于视电阻率法测试技术的隧道围岩松动圈测定

[J].矿业工程研究,374):58-64.

[本文引用: 1]

唐登志白根铭陈爽2023.

深埋隧道岩爆影响因素及其预测研究

[J].四川水力发电,422):11-17.

[本文引用: 1]

王旭李小萌2023.

巷道围岩松动圈测试技术现状及发展趋势研究

[J].煤炭与化工,461):4-7.

[本文引用: 1]

王勇吴爱祥杨军2023.

深部金属矿开采关键理论技术进展与展望

[J].工程科学学报,458):1281-1292.

[本文引用: 1]

谢和平高峰鞠杨2015.

深部岩体力学研究与探索

[J].岩石力学与工程学报,3411):2161-2178.

[本文引用: 1]

余涛朱宁波姚志刚2023.

深埋水平岩层隧道开挖稳定性分析及控制

[J].中国安全生产科学技术,194):93-99.

[本文引用: 1]

朱志洁张宏伟陈蓥2014.

基于MPSO-SVM巷道围岩松动圈预测研究

[J].计算机工程与应用,5012):1-5.

[本文引用: 1]

/