CN114943125A - 基于XGBoost优化算法的隧道围岩参数智能反演分析方法 - Google Patents
基于XGBoost优化算法的隧道围岩参数智能反演分析方法 Download PDFInfo
- Publication number
- CN114943125A CN114943125A CN202210686819.2A CN202210686819A CN114943125A CN 114943125 A CN114943125 A CN 114943125A CN 202210686819 A CN202210686819 A CN 202210686819A CN 114943125 A CN114943125 A CN 114943125A
- Authority
- CN
- China
- Prior art keywords
- parameter
- value
- sample
- algorithm
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/01—Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Medical Informatics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于隧道工程稳定性分析领域,具体公开一种基于XGBoost优化算法的隧道围岩参数智能反演分析方法,包括以下步骤:S1:建立数值模拟计算模型,获得围岩参数反演样本库;S2:对样本库中的位移、参数数据进行相关性和敏感性分析,对各参数作为待反演参数的可行性进行评价,以及位移特征组合的合理性进行量化评价;S3:利用CART算法确定待反演参数,进行位移特征组合筛选;S4:以CART算法为基学习器,建立隧道围岩参数智能反演分析的XGBoost集成算法模型;S5:应用控制单一变量法和贝叶斯优化法,对XGBoost算法进行优化;S6:将位移特征输入到XGBoost集成算法模型,得到预测的围岩参数值。本发明通过超参数寻优实现了XGBoost集成算法模型的优化,模型的稳定性和预测准确率高。
Description
技术领域
本发明属于隧道稳定性分析领域,具体为一种基于XGBoost优化算法的隧道围岩参数智能反演分析方法。
背景技术
采用数值模拟对隧道进行稳定性分析是一种比较常用的手段,可以较容易获取所需的隧道围岩应力、位移和稳定性等结果。但是,数值模拟计算模型计算结果的准确性,很大程度上受制于围岩参数的输入值。在实际的隧道工程勘察和施工中,受制于工程造价、施工工期等多方面因素的影响,隧道围岩参数难以通过室内试验和原位试验而及时准确地获取。
利用隧道围岩位移值来反演围岩参数的隧道围岩反演分析方法,在20世纪70年代被提出。经过几十年的不断发展和实践,大量的研究结论证明了基于位移值进行围岩参数反演分析的可行性和有效性。
随着数学和计算机科学的发展,利用机器学习算法来处理复杂的数值型变量数据,并探究变量之间的映射关系的方法,开始应用于生产生活的各个领域。从90年代末至今,利用数值模拟分析和各种机器学习算法进行隧道围岩参数反演分析,成为隧道稳定性分析领域的主流方法。
通常,利用数值模拟和机器学习算法进行参数反演的技术流程为:
1.建立拟分析隧道段的三维隧道开挖支护数值模拟计算模型;
2.根据需要确定力学参数的设计组合(例如利用正交试验方案对弹性模量、泊松比、容重、粘聚力、内摩擦角等参数进行组合设计),将参数组合代入数值模拟计算模型计算并提取多组位移特征值,建立围岩参数-位移的样本库;
3.确定待反演参数并进行位移特征筛选;
4.基于机器学习算法,建立围岩位移特征与待反演参数之间的映射关系。将样本库数据划分为训练集和测试集,利用训练集数据训练机器学习算法模型,利用测试集数据检验模型的准确率和可用性。
5.将现场实测的围岩位移代入训练好的模型,即可得到预测的围岩参数值。
从上述流程可知,该项技术工作的关键问题,在于所采用的机器学习算法。在确保样本库可靠性的前提下,如何提高参数反演模型的预测准确率以及模型的稳定性,是该项技术工作需要关注的重点。
对于利用位移值反演隧道围岩参数的研究,前人进行了大量的研究工作,采用了许多不同的机器学习算法,例如人工神经网络、免疫高斯过程算法、差异进化算法、支持向量机算法等。但是,前人研究中普遍使用的是单一算法,单一算法具有预测准确率有限,无法充分发挥算法的预测准确率等问题。虽然已有研究成果同样使用XGBoost集成算法进行类似分析,但是,其未对XGBoost集成算法进行超参数优化,模型的稳定性和预测准确率低。并且,现有研究成果未对围岩各参数作为待反演参数的可行性和位移特征组合的合理性进行评价。
发明内容
本发明的目的在于:针对目前利用隧道位移值进行围岩参数反演分析的研究中,所存在的预测准确率低等问题,提出利用XGBoost集成算法提高模型预测的准确率;在基于机器学习算法进行参数反演时,并不是所有的围岩参数都适合作为待反演参数,因此本发明对各参数作为待反演参数的可行性进行了评价;针对位移特征值高度相关导致反演模型稳定性差等问题,本发明进行位移特征筛选;针对集成算法本身的优化问题,本发明应用控制单一变量法和贝叶斯优化法对XGBoost集成算法的超参数进行优化,提高模型的稳定性和预测能力。
本发明一种基于XGBoost优化算法的隧道围岩参数智能反演分析方法,包括以下步骤:
S1:建立隧道开挖支护的三维有限元数值模拟计算模型,获得围岩参数反演的样本库;
S2:对样本库中的位移、参数数据进行特征相关性和敏感性分析,对各参数作为待反演参数的可行性进行评价,对位移特征组合的合理性进行评价;
S3:利用CART算法进行参数反演分析,进行位移特征筛选;
S4:以CART算法为基学习器,建立隧道围岩参数智能反演分析的XGBoost集成算法模型;
S5:应用控制单一变量法和贝叶斯优化法进行超参数优化,得到优化的隧道围岩参数智能反演分析的XGBoost集成算法模型;
S6:将位移特征输入到训练好的围岩参数智能反演分析XGBoost模型,得到预测的围岩参数值。
作为本发明优选的实施方式,在S1中,对隧道所在岩层的围岩参数进行正交试验设计,获得多组围岩参数组合;将多组围岩参数输入数值模拟计算模型,并提取围岩变形稳定阶段的位移值,获得围岩参数反演的样本库。
作为本发明优选的实施方式,在S2中,特征相关性分析采用皮尔逊相关性系数进行计算,两个变量X和Y之间的皮尔逊相关性系数为两个变量之间协方差和标准差的商,皮尔逊相关性系数的绝对值越接近1,则表示两个特征之间具有越强的线性关系,皮尔逊相关性系数公式为:
作为本发明优选的实施方式,在S2中,取隧道所在岩层参数均值、参数均值的90%和参数均值的110%,分别代入数值模拟计算模型,计算相应的位移,通过比较不同参数变化相同比例时的位移变化量,判断围岩的参数敏感性。
作为本发明优先的实施方式,在S3中,CART算法是一种二叉决策树模型,其具体计算过程如下:
(1)选择切分变量j,切分变量为各个位移特征值;
(2)根据切分原则,选择切分点s;
(3)针对任一切分变量的任一切分点si,将训练样本x切分为2个样本子集R1i和R2i;
其中,R1i={xj|xj≤si},j=1~n;R2i={xj|xj>si},j=1~n;i为切分点的编号;
(4)取c1和c2分别为R1i和R2i中样本目标值y的平均值:
(5)根据平方误差最小化原则,应用筛选最优的切分变量和切分点:
(6)根据筛选的切分变量j和切分点s,将训练样本划分为R1和R2两个样本子集;
(7)在样本子集R1和R2中,重复步骤(1)~(6),将子集R1和R2进一步划分为更小的子集;
(8)重复步骤(1)~(7),对子集进一步划分,直到满足结束条件(如达到树的最大深度,叶子节点对应的样本数达到最小的样本数等);
(9)最终,通过上述划分,将训练样本输入空间划分为R1、R2、R3、...Rm等m个叶子节点,生产决策树:
作为本发明优选的实施方式,在S4中,以CART算法为基学习器,构建XGBoost集成算法模型如下:
(1)对CART算法进行迭代,每次迭代是对上一次得到CART决策树的残差进行拟合,以使残差最小化;残差计算公式为:
r=y-ft-1(x)
其中,r为残差,y为样本输出的真实值,t为迭代次数,x为样本的输入变量,ft-1(x)为第t-1次迭代时CART决策树模型的预测值;
(2)进行第0次迭代:
其中,f0(x)为第0次迭代的预测值,L为目标函数,yi为样本输出变量的真实值,c为目标函数的参数,本专利中取其为样本输出变量的平均值;
(3)进行第1次迭代:以样本输出变量真实值与样本输出变量平均值c的差值,作为样本的输出值y1i;根据指定的目标函数,应用CART算法,构造第1个CART决策树模型,获得样本输出变量的预测值f1(x);
(4)进行第2次迭代:以第1次迭代样本原始输出值y1i与第1次迭代样本输出变量预测值f1(x)的差值,作为样本的输出值y2i;根据指定的目标函数,应用CART算法,构造第2个CART决策树模型,获得样本输出变量的预测值f2(x);
(5)按照与第2次迭代同样的方法,构造第3,4,...,M个CART决策树模型,获得样本输出变量的预测值f3(x),f4(x),...,fM(x),直到残差满足要求;
(6)根据上述M个CART决策树模型,按照以下公式计算预测样本的输出值:
作为本发明优选的实施方式,在S4中,XGBoost算法的目标函数是在常规的平方误差损失函数的基础上,加上一个正则项而得到的:
其中,T为叶子节点的数量;γ和λ分别是叶子节点和叶子节点权重的正则化系数;为叶子节点j下各样本的一阶导数之和;为损失函数对样本输出变量预测值的一阶导数;为叶子节点下j各样本的二阶导数之和;为损失函数对样本输出变量预测值的二阶导数;Ij={i|q(xi)=j}为叶子节点j下的样本集;是前t-1轮学习器累加以后给出的第t-1轮的预测值;yi是样本输出变量的原始值;i为样本的编号。
作为本发明优选的实施方式,在S5中,通过控制单一变量法缩小超参数的取值范围;在控制单一变量法缩小后的取值范围内,通过贝叶斯优化法进一步确定超参数的最优取值。
控制单一变量法的具体过程如下:
(1)分别针对每一个超参数,对超参数取值进行等值离散化,获取一定数量的等值点;
(2)将离散化后的超参数取值代入XGBoost集成算法模型,获取每一个模型的可决系数和泛化误差;
(3)根据可决系数和泛化误差的变化规律,获取每一个超参数的最优取值范围。
在缩小后的超参数范围内,利用贝叶斯优化法对超参数进行寻优的过程为:(1)建立超参数样本为x=[x1,x2,x3,...,xn],XGBoost集成算法模型的目标函数取值为y=[y1,y2,y3,...,yn],y服从多维正态分布;
(2)建立高斯过程的核函数:
(3)根据核函数,建立核向量:
(4)建立y的似然函数:
其中,μ为x的均值;
通过梯度下降法对似然函数求极小值,计算核函数的参数,得到最优的核矩阵;
(5)根据最优核矩阵,建立y的后验概率与先验概率和输入变量之间的高斯过程:
其中y*为后一步迭代得到的输出变量取值,K*为后一步迭代的核向量;
上述高斯过程,建立了模型的目标函数值与超参数之间的函数关系;
(6)以上述高斯过程为概率代理模型,建立POI为习得函数:
其中,f(X)为X的目标函数值,即上述推导得到的高斯过程;f(X+)为目前最优的目标函数值,即f(X)的后验分布;μ(x)和σ(x)为高斯过程所得目标函数的均值和方差;ξ为trade-off系数,控制超参数的优化方向偏向于“开发”或者“搜索”;基于上述POI习得函数,应用蒙特卡洛法进行上述最优超参数的搜寻;
(7)将根据习得函数搜寻的最优超参数带入概率代理模型,建立新的高斯过程;
(8)以新的高斯过程为概率代理模型,建立新的POI,应用蒙特卡洛法进行新一轮的超参数寻优;
(9)重复上述操作,直到满足迭代停止条件,获取最终的最优超参数组合。
解释:POI为Probability of improvement的首字母简称。
综上所述,本发明突出的技术优势是:
1)进行特征相关性分析和参数敏感性分析为确定待反演参数提供支撑,因为并不是所有的围岩参数都可以利用位移进行反演。
2)进行位移特征筛选可以提高模型的稳定性,尽量避免位移特征之间高度相关对模型产生的不利影响,同时也考虑预测准确率来选择最优位移特征组合。
3)CART算法作为非线性算法,可以很好的规避线性模型算法位移特征高度相关对模型的不利影响,CART算法较线性算法在参数反演应用中具有更高的预测准确率,有较好的适用性。
4)基于CART算法建立隧道围岩参数反演智能分析的XGBoost集成算法模型,可以充分发挥集成算法的迭代优化作用,XGBoost集成算法重点关注了CART算法每轮迭代分错的节点,并且调整节点的权重值,可以有效地提高基学习器CART算法的预测准确率。
5)综合利用控制单一变量法和贝叶斯优化法对XGBoost集成算法模型进行超参数优化,可以显著提高模型的预测准确性和稳定性,提高模型的拟合效果。与传统的采用单一方法相比,例如网格搜索、随机搜索和贝叶斯优化等,综合优化方法更加高效和准确。
附图说明
图1为本发明基于XGBoost优化算法的隧道围岩参数智能反演分析方法的流程图;
图2为本发明实施例中隧道开挖的三维模型岩层划分示意图;
图3为本发明实施例中隧道开挖细部网格尺寸示意图;
图4为本发明实施例中隧道围岩开挖以后的应力图;
图5为本发明实施例中隧道围岩开挖以后的应变云图;
图6为本发明实施例中位移提取点位示意图;
图7为本发明基于XGBoost优化算法的隧道围岩参数智能反演分析方法中的CART算法流程图;
图8为本发明基于XGBoost优化算法的隧道围岩参数智能反演分析方法中的XGBoost集成算法流程图;
图9为本发明实施例中超参数n_estimators取不同值时其可决系数得分曲线图;
图10为本发明实施例中超参数n_estimators取不同值时其泛化误差得分曲线;
图11为本发明实施例中泊松比的预测值与真实值的柱状对比图。
具体实施方式
下面结合附图,对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例基于XGBoost优化算法的隧道围岩参数智能反演分析方法,如附图1所示,包括:S1:利用有限元软件,建立隧道开挖三维有限元的数值模拟计算模型,获得围岩参数反演的样本库;S2:对样本库中的位移、参数数据进行特征相关性和敏感性分析,对各参数作为待反演参数的可行性进行评价,对位移特征组合的合理性进行评价;S3:利用CART算法进行参数反演分析,进行位移特征筛选;S4:以CART算法为基学习器,建立隧道围岩参数智能反演分析的XGBoost集成算法模型;S5:应用控制单一变量法和贝叶斯优化法进行超参数优化,得到优化的隧道围岩参数智能反演分析的XGBoost集成算法模型;S6:将位移特征输入到训练好的围岩参数智能反演分析XGBoost模型,得到预测的围岩参数值。
以下通过具体的实例对本发明方案作更加详细的说明:
第一步,首先确定以某隧道工程K4+173到K4+205标段的洞身主线分离式右线隧洞为建模对象,研究K4+173截面开挖支护稳定后的围岩位移和参数。该段隧道埋深154米,隧道围岩强度较高,隧道围岩等级为三级,采用钻爆法进行施工,右线隧道已经推进到足够远的位置。左右线隧道之间隔较远,符合主线分离式隧道对间距的要求,相互间的影响可以忽略不计。采用全断面法开挖、无支护的施工方式,对隧道开挖支护施工阶段进行数值模拟。
样本库的位移数据,根据具体开挖截面的大小以及支护情况进行提取,一般包括拱顶沉降、拱底隆起、水平收敛等多组位移特征;样本库中的参数是岩土体采用摩尔-库伦本构时的输入参数,包括弹性模量、泊松比、容重、粘聚力、内摩擦角等五个参数。根据实测位移时序曲线判断围岩进入稳定阶段所需要的时间,即:统一确定提取的稳定位移值,为隧道开挖以后第32天的稳定位移值。隧道每天进行一个施工步,每一施工步开挖1米。因此,隧道开挖三维模型在深度方向的扩展距离是32米。在宽度和高度方向上,需满足扩展范围约为开挖尺寸的五倍,在本例中,所建隧道开挖三维模型水平向的长度为199.5米,垂直向高度为218.8米。三维模型网格划分采用六面体网格,将隧道开挖部分和所在岩层的网格尺寸控制在0.6米,临近岩层的网格尺寸控制在0.8米,其余部分岩层的网格尺寸控制在1.2米。隧道开挖数值模拟计算模型由有限元计算软件建模获得,可选用的有限元计算软件包括但不限于Midas GTS NX、FLAC 3D、ABAQUS等。数值模拟计算模型的岩土体的本构关系采用摩尔-库伦本构模型,所建的三维模型的岩层划分示意图如附图2所示,隧道开挖细部网格尺寸示意图如附图3所示。
由前期地勘报告可知,该段隧道的岩层共有六层,从上到下分别是强风化砂岩、中风化砂岩、黑色页岩、灰色砂岩、黑色页岩、灰色砂岩。建模时,前五层岩层的围岩参数取用地勘报告中该岩层参数的平均值。第六层岩层灰色砂岩为隧道开挖所在的岩层,该岩层的参数按照正交设计进行输入。根据勘察报告和实际情况确定该岩层参数的最大取值和最小取值,以确定各参数的取值范围。各岩层参数取值如表1所示。
表1
第二步:对隧道开挖所在岩层的灰色砂岩的五个参数进行正交试验设计。根据第一步确定的隧道开挖所在岩层各参数的取值范围进行正交设计,本例采用SPSS软件,对每个参数取6个取值水平,共五个参数。因此正交设计共生成49组围岩参数组合,如表2所示为展示的部分五因素六水平正交试验参数组合表。
表2
确定位移特征提取点的位置。隧道开挖所在岩层围岩的各参数选取参数的均值或中位数,代入数值模拟计算模型计算,观察隧道围岩开挖以后的应力和应变云图,应力图如图4所示,应变云图如图5所示,从图中选择应力-应变结果较大的区域合理布设位移监测点。位移监测点的位移如图6所示。
将正交试验设计的49组围岩参数组合,依次代入数值模拟计算模型进行计算,并提取相应提取点位的位移值,构成参数反演的参数-位移样本库,部分参数-位移数据统计在表3中。
表3
第三步:对数据库中的参数-位移数据进行特征相关性分析。特征相关性分析采用皮尔逊相关性系数进行计算,两个变量X和Y之间的皮尔逊相关性系数为两个变量之间协方差和标准差的商,皮尔逊相关性系数的绝对值越接近1,则表示两个特征之间具有越强的线性关系,皮尔逊相关性系数公式为:
各变量间的皮尔逊相关性系数得分如下表4所示。
表4
第四步:进行参数敏感性分析,首先计算隧道开挖所在岩层(第六层:灰色砂岩)地勘报告所提供的各参数的平均值,然后分别计算各参数的0.9倍和1.1倍。利用控制单一变量法,每次将一种参数的0.9倍、1倍、1.1倍和其他参数的1倍进行组合,每次得到3种参数组合形式,将参数组合代入隧道开挖数值模拟模型进行计算,提取位移特征值,然后用1.1倍参数产生的位移特征值减去0.9倍参数产生的位移值,所得的差取绝对值,用差的绝对值除以1倍参数产生的位移特征值,则得到各参数变化相同比例时各位移特征的变化百分数,该百分数量化了参数变化时产生位移的大小。以容重和弹模为例,参数敏感性的计算过程如表5所示,参数敏感性的计算结果如表6所示。
表5
表6
第五步:根据变量相关性和参数敏感性的计算结果进行分析:由变量相关性的计算结果可知,弹性模量、泊松比和位移特征之间的相关性得分较高,这两个参数与位移特征之间的相关性程度较高;容重、粘聚力、内摩擦角和位移特征之间的相关性得分较低,这三个参数与位移之间的相关性程度较低;由参数敏感性分析的计算结果可知,弹性模量、泊松比变化时,产生的位移变化较大,容重变化时,产生的位移变化很小,粘聚力、内摩擦角变化时,位移不产生变化。综合参数敏感性和变量相关性的分析结果,初步确定待反演参数为弹性模量和泊松比。
从变量相关性计算的分析结果发现,自变量之间即位移特征之间存在较高的线性相关性,其中收敛二和收敛三之间的相关性系数得分为0.96,拱顶降和拱底隆起之间的相关性系数为0.98,自变量之间高度相关,采用线性算法时,会出现模型的可靠度差、置信度低等问题。因此,决定采用非线性算法进行参数反演分析。
第六步:利用CART算法确定待反演参数,进行位移特征组合筛选。CART算法流程如附图7所示。CART算法作为非线性算法,可以很好地避免线性算法模型位移特征之间线性相关和多重共线性对模型的不利影响。CART回归是递归生成二叉树的过程,分裂的过程使用平方误差最小化的准则,其具体计算过程如下:
(1)选择切分变量j,切分变量为各个位移特征值;
(2)根据切分原则,选择切分点s;
(3)针对任一切分变量的任一切分点si,将训练样本x切分为2个样本子集R1i和R2i;
其中,R1i={xj|xj≤si},j=1~n;R2i={xj|xj>si},j=1~n;i为切分点的编号;
(4)取c1和c2分别为R1i和R2i中样本目标值y的平均值:
(5)根据平方误差最小化原则,应用筛选最优的切分变量和切分点:
(6)根据筛选的切分变量j和切分点s,将训练样本划分为R1和R2两个样本子集;
(7)在样本子集R1和R2中,重复步骤(1)~(6),将子集R1和R2进一步划分为更小的子集;
(8)重复步骤(1)~(7),对子集进一步划分,直到满足结束条件(如达到树的最大深度,叶子节点对应的样本数达到最小的样本数等);
(9)最终,通过上述划分,将训练样本输入空间划分为R1、R2、R3、...Rm等m个叶子节点,生产决策树:
表7所示为CART算法参数反演结果统计表。其中,用可决系数来表征参数反演的准确率。可决系数r2的计算公式如下:
表7
如表7所示,CART算法反演弹性模量的准确率较低,因此使用CART算法反演弹性模量的可行性较差,使用CART算法反演泊松比的准确率较高,建立反演泊松比的预测模型的可行性较大。
为了提高CART算法模型预测的准确率,采用XGBoost集成算法对单一算法进行集成。但XGBoost集成算法要求单一算法的预测准确率高于0.5。因此进一步确定待反演参数为泊松比。
为了进一步确定最优的位移特征组合,对位移特征进行筛选。具体操作是对拱顶沉降、拱低隆起、水平收敛一、水平收敛二、水平收敛三五个位移特征组合进行枚举,共19种位移参数组合,将这些参数组合反演泊松比的计算结果统计如表8所示。
表8
由各位移特征组合反演泊松比的计算结果可知,特征组合选为全部位移特征时的预测准确率最高,则此位移特征组合为最优位移特征组合。
第七步:以CART算法作为基学习器,建立基于XGBoost集成算法的隧道围岩参数反演模型。XGBoost集成算法流程如附图8所示。XGBoost算法是一种集成算法,基于Boosting思想将基学习器算法进行集成。XGBoost每次迭代都会调整基学习器CART算法叶子节点的权重,重点关注上一次迭代分错的节点,并修改该节点的的权重值,XGBoost最终的预测结果为所有节点的权重和。
构建XGBoost集成算法模型的具体实施过程如下:
(1)对CART算法进行迭代,每次迭代是对上一次得到CART决策树的残差进行拟合,以使残差最小化;残差计算公式为:
r=y-ft-1(x)
其中,r为残差,y为样本输出的真实值,t为迭代次数,x为样本的输入变量,ft-1(x)为第t-1次迭代时CART决策树模型的预测值;
(2)进行第0次迭代:
其中,f0(x)为第0次迭代的预测值,L为目标函数,yi为样本输出变量的真实值,c为目标函数的参数,本专利中取其为样本输出变量的平均值;
(3)进行第1次迭代:以样本输出变量真实值与样本输出变量平均值c的差值,作为样本的输出值y1i;根据指定的目标函数,应用CART算法,构造第1个CART决策树模型,获得样本输出变量的预测值f1(x);
(4)进行第2次迭代:以第1次迭代样本原始输出值y1i与第1次迭代样本输出变量预测值f1(x)的差值,作为样本的输出值y2i;根据指定的目标函数,应用CART算法,构造第2个CART决策树模型,获得样本输出变量的预测值f2(x);
(5)按照与第2次迭代同样的方法,构造第3,4,...,M个CART决策树模型,获得样本输出变量的预测值f3(x),f4(x),...,fM(x),直到残差满足要求;
(6)根据上述M个CART决策树模型,按照以下公式计算预测样本的输出值:
XGBoost算法的目标函数是在常规的平方误差损失函数的基础上,加上一个正则项而得到的:
其中,T为叶子节点的数量;γ和λ分别是叶子节点和叶子节点权重的正则化系数;为叶子节点j下各样本的一阶导数之和;为损失函数对样本输出变量预测值的一阶导数;为叶子节点下j各样本的二阶导数之和;为损失函数对样本输出变量预测值的二阶导数;Ij={i|q(xi)=j}为叶子节点j下的样本集;是前t-1轮学习器累加以后给出的第t-1轮的预测值;yi是样本输出变量的原始值;i为样本的编号。
表9为XGBoost集成算法反演泊松比效果评价统计表。
表9
XGBoost集成算法超参数众多,在面对具体问题时,具体问题的数据特点显然是和XGBoost的默认的超参数组合不是最匹配的。XGBoost集成算法的超参数和具体问题的数据特点不匹配时,会影响XGBoost集成算法的稳定性和预测准确率。因此,为了获得稳定性更好、预测准确率更高的XGBoost集成算法模型,有必要对XGBoost集成算法进行超参数优化。因此接下来对XGBoost集成算法进行优化,以期提升模型的预测准确率。
第八步:对XGBoost集成算法进行超参数优化,针对参数反演的数据特征,从XGBoost集成算法中筛选出可以进行调优的超参数,可进行调优的超参数见表10。
表10
首先利用控制单一变量法缩小超参数的择优范围,其具体实施过程如下:
(1)分别针对每一个超参数,对超参数取值进行等值离散化,获取一定数量的等值点;
(2)将离散化后的超参数取值代入XGBoost集成算法模型,获取每一个模型的可决系数和泛化误差;
(3)根据可决系数和泛化误差的变化规律,获取每一个超参数的最优取值范围。
分别绘制表10中超参数不同取值时的可决系数得分曲线和泛化误差得分曲线。如图9所示为超参数n_estimators取不同值时的可决系数得分曲线,如图10所示为超参数n_estimators取不同值时的泛化误差得分曲线。
根据超参数n_estimators在取不同值时可决系数和泛化误差的得分曲线,当n_estimators取值在15左右时,其可决系数得分达到峰值左右,其泛化误差得分也达到最低值附近,随着该超参数继续增大,其可决系数和泛化误差都保持不变。因此,根据可决系数和泛化误差的得分曲线,可以缩小该超参数的优化搜索范围为[5,20]。当该超参数的取值大于20时,不仅不会增加模型的准确率减小误差,同时还会造成过拟合,并且增加模型的复杂度,降低模型的运行效率。
通过控制单一变量法确定的XGBoost集成算法的超参数取值范围见表11。
表11
利用贝叶斯优化法对超参数组合进行寻优。贝叶斯优化法以贝叶斯定理为基础,由概率代理函数和采集函数两大部分组成,为了找到目标函数的极值,每次采集函数对数据进行一次采集后都会更新概率代理函数,采集函数根据概率代理函数的反馈,在有可能出现极值的区域进行采集。每一次采集后再次更新概率代理函数,如此进行循环迭代。贝叶斯优化算法充分利用了历史计算数据,少走了很多弯路,其计算量比其他算法少很多,是一种更聪明的优化方法。
利用贝叶斯优化法对超参数进行寻优,其具体过程为:
(1)建立超参数样本为x=[x1,x2,x3,...,xn],XGBoost模型的目标函数取值为y=[y1,y2,y3,...,yn],y服从多维正态分布;
(2)建立高斯过程的核函数:
(3)根据核函数,建立核向量:
(4)建立y的似然函数:
其中,μ为x的均值;
通过梯度下降法对似然函数求极小值,计算核函数的参数,得到最优的核矩阵;
(5)根据最优核矩阵,建立y的后验概率与先验概率和输入变量之间的高斯过程:
其中y*为后一步迭代得到的输出变量取值,K*为后一步迭代的核向量;
上述高斯过程,建立了模型的目标函数值与超参数之间的函数关系;
(6)以上述高斯过程为概率代理模型,建立POI(Probability of improvement)为习得函数:
其中,f(X)为X的目标函数值,即上述推导得到的高斯过程;f(X+)为目前最优的目标函数值,即f(X)的后验分布;μ(x)和σ(x)为高斯过程所得目标函数的均值和方差;ξ为trade-off系数,控制超参数的优化方向偏向于“开发”或者“搜索”;基于上述POI习得函数,应用蒙特卡洛法进行上述最优超参数的搜寻;
(7)将根据习得函数搜寻的最优超参数带入概率代理模型,建立新的高斯过程;
(8)以新的高斯过程为概率代理模型,建立新的POI,应用蒙特卡洛法进行新一轮的超参数寻优;
(9)重复上述操作,直到满足迭代停止条件,获取最终的最优超参数组合。
基于Python,从bayes_opt模块调用BayesianOptimization贝叶斯优化算法,经过30轮迭代寻优计算以后输出的超参数组合见表12。
表12
第九步:将表12输出的超参数组合对XGBoost集成算法的超参数进行定义。然后将数据库中80%的数据集作为训练集数据代入定义超参数的XGBoost集成算法进行训练,利用剩下的20%的数据作为测试集数据对模型进行测试。泊松比的参数反演分析评价指标见表13。
表13
部分预测值与真实值的柱状对比图如图11所示。
为了验证预测参数结果的实用性,将部分预测参数带入数值模拟计算模型,将预测参数产生的位移与真实参数产生的位移进行对比,误差小于10%则证明预测参数的可用性。预测参数与真实参数所产生的位移的比较见表14。可见,所有算例的误差都小于10%,说明基于优化XGBoost集成算法的隧道围岩参数反演模型具有实用价值。
表14
Claims (8)
1.基于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,包括以下步骤:
S1:建立隧道开挖支护的三维有限元数值模拟计算模型,获得围岩参数反演的样本库;
S2:对样本库中的位移、参数数据进行特征相关性和敏感性分析,对各参数作为待反演参数的可行性进行评价,对位移特征组合的合理性进行评价;
S3:利用CART算法进行参数反演分析,进行位移特征筛选;
S4:以CART算法为基学习器,建立隧道围岩参数智能反演分析的XGBoost集成算法模型;
S5:应用控制单一变量法和贝叶斯优化法进行超参数优化,得到优化的隧道围岩参数智能反演分析的XGBoost集成算法模型;
S6:将位移特征输入到训练好的围岩参数智能反演分析XGBoost模型,得到预测的围岩参数值。
2.根据权利要求1所述的基于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,在S1中,对隧道所在岩层的围岩参数进行正交试验设计,获得多组围岩参数组合;将多组围岩参数输入隧道开挖支护的数值模拟计算模型,并提取围岩变形稳定阶段的位移值,获得围岩参数反演的样本库。
4.根据权利要求1所述的基于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,在S2中,取隧道所在岩层参数均值、参数均值的90%和参数均值的110%,分别代入数值模拟计算模型,计算相应的位移,通过比较不同参数变化相同比例时的位移变化量,判断围岩的参数敏感性。
5.根据权利要求1所述的基于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,在S3中,CART算法是一种二叉决策树模型,具体计算过程如下:
(1)选择切分变量j,切分变量为各个位移特征值;
(2)根据切分原则,选择切分点s;
(3)针对任一切分变量的任一切分点si,将训练样本x切分为2个样本子集R1i和R2i;
其中,R1i={xj|xj≤si},j=1~n;R2i={xj|xj>si},j=1~n;i为切分点的编号;
(4)取c1和c2分别为R1i和R2i中样本目标值y的平均值:
(5)根据平方误差最小化原则,应用筛选最优的切分变量和切分点:
(6)根据筛选的切分变量j和切分点s,将训练样本划分为R1和R2两个样本子集;
(7)在样本子集R1和R2中,重复步骤(1)~(6),将子集R1和R2进一步划分为更小的子集;
(8)重复步骤(1)~(7),对子集进一步划分,直到满足结束条件;
(9)最终,通过上述划分,将训练样本输入空间划分为R1、R2、R3、...Rm等m个叶子节点,生产决策树:
6.根据权利要求1所述的基于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,在S4中,以CART算法为基学习器,构建XGBoost集成算法模型如下:
(1)对CART算法进行迭代,每次迭代是对上一次得到CART决策树的残差进行拟合,以使残差最小化;残差计算公式为:
r=y-ft-1(x)
其中,r为残差,y为样本输出的真实值,t为迭代次数,x为样本的输入变量,ft-1(x)为第t-1次迭代时CART决策树模型的预测值;
(2)进行第0次迭代:
其中,f0(x)为第0次迭代的预测值,L为目标函数,yi为样本输出变量的真实值,c为目标函数的参数,取样本输出变量的平均值;
(3)进行第1次迭代:以样本输出变量真实值与样本输出变量平均值c的差值,作为样本的输出值y1i;根据指定的目标函数,应用CART算法,构造第1个CART决策树模型,获得样本输出变量的预测值f1(x);
(4)进行第2次迭代:以第1次迭代样本原始输出值y1i与第1次迭代样本输出变量预测值f1(x)的差值,作为样本的输出值y2i;根据指定的目标函数,应用CART算法,构造第2个CART决策树模型,获得样本输出变量的预测值f2(x);
(5)按照与第2次迭代同样的方法,构造第3,4,...,M个CART决策树模型,获得样本输出变量的预测值f3(x),f4(x),...,fM(x),直到残差满足要求;
(6)根据上述M个CART决策树模型,按照以下公式计算预测样本的输出值:
8.根据权利要求1所述的于XGBoost优化算法的隧道围岩参数智能反演分析方法,其特征在于,在S5中,通过控制单一变量法缩小超参数的取值范围;在控制单一变量法缩小后的取值范围内,通过贝叶斯优化法进一步确定超参数的最优取值;利用贝叶斯优化法对超参数进行寻优的过程为:
(1)建立超参数样本为x=[x1,x2,x3,...,xn],XGBoost集成算法模型的目标函数取值为y=[y1,y2,y3,...,yn],y服从多维正态分布;
(2)建立高斯过程的核函数:
(3)根据核函数,建立核向量:
(4)建立y的似然函数:
其中,μ为x的均值;
通过梯度下降法对似然函数求极小值,计算核函数的参数,得到最优的核矩阵;
(5)根据最优核矩阵,建立y的后验概率与先验概率和输入变量之间的高斯过程:
其中y*为后一步迭代得到的输出变量取值,K*为后一步迭代的核向量;
上述高斯过程,建立了XGBoost集成算法模型的目标函数值与超参数之间的函数关系;
(6)以上述高斯过程为概率代理模型,建立POI为习得函数:
其中,f(X)为X的目标函数值,即上述推导得到的高斯过程;f(X+)为目前最优的目标函数值,即f(X)的后验分布;μ(x)和σ(x)为高斯过程所得目标函数的均值和方差;ξ为trade-off系数,控制超参数的优化方向偏向于“开发”或者“搜索”;基于上述POI习得函数,应用蒙特卡洛法进行上述最优超参数的搜寻;
(7)将根据习得函数搜寻的最优超参数带入概率代理模型,建立新的高斯过程;
(8)以新的高斯过程为概率代理模型,建立新的POI,应用蒙特卡洛法进行新一轮的超参数寻优;
(9)重复上述操作,直到满足迭代停止条件,获取最终的最优超参数组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210686819.2A CN114943125A (zh) | 2022-06-16 | 2022-06-16 | 基于XGBoost优化算法的隧道围岩参数智能反演分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210686819.2A CN114943125A (zh) | 2022-06-16 | 2022-06-16 | 基于XGBoost优化算法的隧道围岩参数智能反演分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114943125A true CN114943125A (zh) | 2022-08-26 |
Family
ID=82910830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210686819.2A Pending CN114943125A (zh) | 2022-06-16 | 2022-06-16 | 基于XGBoost优化算法的隧道围岩参数智能反演分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114943125A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115563831A (zh) * | 2022-10-20 | 2023-01-03 | 北京云庐科技有限公司 | 隧道地层力学参数获取方法、装置、电子设备及存储介质 |
GB2619567A (en) * | 2022-06-09 | 2023-12-13 | Ocado Innovation Ltd | 3-D lattice optimisation |
CN117711140A (zh) * | 2024-02-06 | 2024-03-15 | 山东大学 | 基于多源数据融合的隧洞突水灾害定时预警方法及系统 |
-
2022
- 2022-06-16 CN CN202210686819.2A patent/CN114943125A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2619567A (en) * | 2022-06-09 | 2023-12-13 | Ocado Innovation Ltd | 3-D lattice optimisation |
GB2619567B (en) * | 2022-06-09 | 2024-07-17 | Ocado Innovation Ltd | 3-D lattice optimisation |
CN115563831A (zh) * | 2022-10-20 | 2023-01-03 | 北京云庐科技有限公司 | 隧道地层力学参数获取方法、装置、电子设备及存储介质 |
CN117711140A (zh) * | 2024-02-06 | 2024-03-15 | 山东大学 | 基于多源数据融合的隧洞突水灾害定时预警方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114943125A (zh) | 基于XGBoost优化算法的隧道围岩参数智能反演分析方法 | |
Zhang et al. | Comprehensive review of machine learning in geotechnical reliability analysis: Algorithms, applications and further challenges | |
CN110837669B (zh) | 基于多源异构数据融合的滑坡不确定模型动态构建方法 | |
CN113610945B (zh) | 一种基于混合神经网络的地应力曲线预测方法 | |
Li et al. | Probabilistic back analysis for improved reliability of geotechnical predictions considering parameters uncertainty, model bias, and observation error | |
Fu et al. | A blocking Markov chain Monte Carlo method for inverse stochastic hydrogeological modeling | |
Zhou et al. | Machine learning models to predict the tunnel wall convergence | |
Laalam et al. | Statistical analysis of the petrophysical properties of the Bakken Petroleum System | |
Wang et al. | Deep learning based closed-loop optimization of geothermal reservoir production | |
Shi et al. | Prediction of mechanical behavior of rocks with strong strain-softening effects by a deep-learning approach | |
Tian et al. | Auxiliary Bayesian updating of embankment settlement based on finite element model and response surface method | |
Deng et al. | A hybrid machine learning optimization algorithm for multivariable pore pressure prediction | |
Razak et al. | Embedding physical flow functions into deep learning predictive models for improved production forecasting | |
NO20200978A1 (en) | Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter | |
Obrzud et al. | Optimization framework for calibration of constitutive models enhanced by neural networks | |
Zeng et al. | Local coupled Markov chain model for simulating varied stratigraphy | |
Xie et al. | A back analysis scheme for refined soil stratification based on integrating borehole and CPT data | |
Millán et al. | Application of artificial neural networks for predicting the bearing capacity of the tip of a pile embedded in a rock mass | |
Shahzadi et al. | Deep neural network-based inverse analysis with application to a rockfill dam | |
García Benítez et al. | Neural networks for defining spatial variation of rock properties in sparsely instrumented media | |
Lee et al. | Quick prediction of tunnel displacements using Artificial Neural Network and field measurement results | |
Silva et al. | Generative network-based reduced-order model for prediction, data assimilation and uncertainty quantification | |
Faramarzi | Intelligent computational solutions for constitutive modelling of materials in finite element analysis | |
JIANG et al. | Inversion iterative correction method for estimating shear strength of rock and soil mass in slope engineering | |
Zhang et al. | BPNN-GA-based intelligence information approach to characterize deep fractured strata |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |