CN113156492B - 一种应用于tbm隧道岩爆灾害实时智能预警方法 - Google Patents
一种应用于tbm隧道岩爆灾害实时智能预警方法 Download PDFInfo
- Publication number
- CN113156492B CN113156492B CN202110432239.6A CN202110432239A CN113156492B CN 113156492 B CN113156492 B CN 113156492B CN 202110432239 A CN202110432239 A CN 202110432239A CN 113156492 B CN113156492 B CN 113156492B
- Authority
- CN
- China
- Prior art keywords
- microseismic
- precursor
- parameter
- hyper
- sequence
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000002243 precursor Substances 0.000 claims abstract description 65
- 238000012544 monitoring process Methods 0.000 claims abstract description 57
- 239000011159 matrix material Substances 0.000 claims abstract description 56
- 230000009467 reduction Effects 0.000 claims abstract description 30
- 230000008569 process Effects 0.000 claims abstract description 19
- 239000011435 rock Substances 0.000 claims abstract description 18
- 238000012549 training Methods 0.000 claims abstract description 10
- 230000005641 tunneling Effects 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 31
- 239000000203 mixture Substances 0.000 claims description 28
- 230000010354 integration Effects 0.000 claims description 15
- 238000009412 basement excavation Methods 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 3
- 229910052731 fluorine Inorganic materials 0.000 claims description 3
- 125000001153 fluoro group Chemical group F* 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 239000003795 chemical substances by application Substances 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 6
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 208000013201 Stress fracture Diseases 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种应用于TBM隧道岩爆灾害实时智能预警方法,其步骤:进行微震监测平台搭建;降维获得降维后的微震数据矩阵;进行降维后的微震数据矩阵中的微震事件聚类;组建前兆微震序列数据库;建立分类回归树;训练分类回归树;对自适应增强集成模型中的超参数寻优,本发明提供了全断面隧道掘进机掘进过程岩爆灾害烈度实时预警方法,克服了传统人工方法效率低、报不准、主观性和经验依赖性强等难题,对全断面隧道掘进机安全高效掘进具有重要意义。
Description
技术领域
本发明涉及隧道掘进机施工技术领域,更具体涉及一种应用于TBM隧道岩爆灾害实时智能预警方法。
背景技术
全断面隧道掘进机是一种高效的隧道开挖机械,特别适用于长距离隧道的施工,它已经在隧道工程中得到越来越广泛的应用。岩爆是地下工程开挖过程中一种常见的动力失稳地质灾害,在时间上具有突发性,在空间上具有随机性,且破坏力极强,特别是在全断面隧道掘进机掘进过程中岩爆的发生往往会造成机毁人亡的巨大损失。然而目前国内外尚无完备的针对全断面隧道掘进机掘进过程岩爆灾害实时智能预警方法,主要存在以下不足:
(1)由于全断面隧道掘进机现场施工条件的限制,在线状狭小的可利用空间中难以搭建高效的微震监测平台。
(2)现有的岩爆灾害智能预警模型难以实现实时预警,受岩体参数(单轴抗压强度、岩石脆性指标、弹性应变能等)快速准确获取技术的限制。
(3)依赖于现场工程师的工程经验,岩爆烈度主要还是由人为确定,主观性强且效率低。
针对上述不足,本发明专利提出了一种应用于TBM隧道基于微震监测的岩爆灾害实时智能预警方法。
发明内容
为解决全断面隧道掘进机安全高效掘进的难题,本发明的目的在于提供了一种应用于TBM隧道岩爆灾害实时智能预警方法。它将人工智能与微震监测相结合,随着微震数据量的不断增加,人工智能为数据处理提供了强有力的有段,进而为岩爆灾害实时智能预警提供了坚实的基础。
为实现上述目的,本发明采用以下技术方案:
一种应用于TBM隧道岩爆灾害实时智能预警方法,包括以下步骤:
步骤1、全断面隧道掘进机掘进过程中进行微震监测平台搭建;
步骤2、对微震监测平台监测获得的微震数据库降维获得降维后的微震数据矩阵X*;
步骤3、利用高斯混合模型进行降维后的微震数据矩阵中的微震事件聚类,为降维后的微震数据矩阵中的微震事件分配聚类标签;
步骤4、根据降维后的微震数据矩阵组建前兆微震序列数据库;
步骤5、建立分类回归树,分类回归树包括根节点、内部节点和叶节点,
根节点用于存储步骤4建立的前兆微震序列数据库中所有的前兆微震序列,内部节点表示属性条件,内部节点将上一层内部节点或根节点中包含的前兆微震序列分割为满足不同属性条件的前兆微震序列子集,叶节点包含岩爆发生与否以及岩爆烈度;
步骤6、训练分类回归树,并通过加权投票法得到自适应增强集成模型的输出;
步骤7、利用贝叶斯优化算法对自适应增强集成模型中的超参数寻优。
如上所述的步骤1包括以下步骤:
步骤1.1、选取监测断,距离开挖掌子面由近到远依次布设有第1组监测断面、第2组监测断面和第3组监测断面,
步骤1.2、每组监测断面在拱肩位置左右对称布置2个传感器;
步骤1.3、全断面隧道掘进机向前掘进设定距离后,将原第3组监测断面前移成为新的第1组监测断面,原第1组监测断面和原第2组监测断面分别作为新的第2组监测断面和新的第3组监测断面。
如上所述的步骤2包括以下步骤:
步骤2.1、构建原始微震数据矩阵X=(xij)m×n,m是微震数据库中微震事件的总数量,n是微震事件的震源参数的总数量,xij是第i个微震事件的第j个震源参数的原始微震数据,
步骤2.3:对标准化后的微震数据计算相关系数矩阵R=(rhl)n×n,其中rhl为第h个震源参数和第l个震源参数之间的皮尔逊相关系数,
步骤2.4:计算相关系数矩阵的特征值λ1,λ2,…,λn及其对应的特征向量v1,v2,…,vn,
步骤2.5:将数值大于1且累计贡献率超过80%的特征值对应的特征向量组成矩阵V,特征值的累计贡献率根据公式2计算,
其中,ηg指g个特征值的累计贡献率,
步骤2.6:根据公式3可以得到降维后的微震数据矩阵X*,
X*=X×V 公式3
其中,X*指降维后的微震数据矩阵;X指原始微震数据矩阵;V指步骤2.5中由特征向量组成的矩阵。
如上所述的步骤3包括以下步骤:
步骤3.1、按照公式4、公式5和公式6构造高斯混合模型P(y|θ),
其中,φk(y|θk)指高斯混合模型中的第k个高斯密度函数;αk(αk≥0)指第k个高斯密度函数的权重系数;指第k个高斯密度函数中的特征参数,且y为微震事件,σk和μk均为高斯密度函数中的待定参数,K为高斯混合模型中高斯密度函数的数量,
期望最大化算法的流程如下:
步骤3.2.2、根据公式7计算高斯混合模型中的第k个高斯密度函数φk(y|θk)对降维后的微震数据矩阵中第i个微震事件yi的响应γik,
其中,γik指高斯混合模型中的第k个高斯密度函数对降维后的微震数据矩阵中第i个微震事件的响应,m是微震数据库中微震事件的总数量,yi为第i个微震事件,
步骤3.3、迭代反演得到高斯混合模型中的参数α={α1,α2,...,αK}和后,利用高斯混合模型可以对降维后的微震数据矩阵中的微震事件进行聚类,为降维后的微震数据矩阵中的微震事件分配聚类标签,还获得不同微震事件属于各个聚类标签的概率。
如上所述的步骤4包括以下步骤:
为降维后的微震数据矩阵中的微震事件分配聚类标签后,将每5个微震事件对应的聚类标签按照微震事件捕捉时间的先后顺序排列即构造1个前兆微震序列,进而建立由前兆微震序列组成的前兆微震序列数据库。
如上所述的步骤5中,采用属性条件A的基尼指数作为分割判据,属性条件A的基尼指数Gini_index(D,A)如公式11所示:
其中,D指由待分割节点所包含的前兆微震序列组成的前兆微震序列数据集,待分割节点包括内部节点或者根节点;Dv指由属性条件A分割前兆微震序列数据集得到的前兆微震序列数据子集,V指前兆微震序列数据集D中样本的类数,
基尼值Gini(Dv)定义为:
其中,pk指前兆微震序列数据子集Dv中第k类样本所占的比例。
如上所述的步骤6包括以下步骤:
步骤6.1、初始化前兆微震序列数据库中前兆微震序列的权重,
步骤6.2、根据步骤5利用前兆微震序列数据库建立的分类回归树,
步骤6.3、计算分类回归树的预测误差,
步骤6.4、计算分类回归树在自适应增强集成模型中的权重,
步骤6.5、根据分类回归树预测结果更新前兆微震序列数据库中前兆微震序列的权重,
步骤6.6、如果已建立的分类回归树数量达到预先设定值,则停止训练分类回归树,并通过加权投票法得到自适应增强集成模型的输出;如果已建立的分类回归树数量未达到预先设定值,则重复步骤6.2~6.5。
如上所述的步骤7包括以下步骤:
步骤7.1、初始化自适应增强集成模型中的超参数,建立超参数数据库,
步骤7.2、根据建立的超参数数据库训练如公式13所示高斯过程代理模型:
f(B)~GP(0,κ(xp,xq)) 公式13
其中,κ(xp,xq)是关于超参数的协方差函数矩阵;xp是超参数数据库中的第p个超参数;xq是超参数数据库中的第q个超参数;f(B)是自适应增强集成模型在不同超参数组合下的精度,GP为高斯过程函数,
步骤7.3、通过公式14所示的最大化概率提升获取函数PI(B)确定下一组待评价超参数,并补充到超参数数据库中,
其中,Φ指标准正态分布的累积分布函数;μ(B)和σ(B)分别指根据步骤7.2建立的高斯过程代理模型得到的自适应增强集成模型在不同超参数组合下的精度的平均值和标准差;f(B*)指当前超参数数据库下自适应增强集成模型在不同超参数组合下的精度的最优值,
步骤7.4、若达到设定的最大迭代次数则终止迭代,否则重复步骤7.2~7.3,迭代结束后,将超参数数据库中实现最高自适应增强集成模型精度的一组超参数作为最优超参数。
本发明与现有技术相比,具有以下优点和效果:
(1)根据全断面隧道掘进机现场施工条件,设计和优化了微震监测传感器的布设方式,满足随掘进随监测的要求,同时还能有效地捕捉掘进过程中产生的岩体微破裂信号。
(2)微震监测是一种原位三维空间实时监测技术,利用微震数据开展岩爆预警研究相比传统利用岩性参数开展岩爆预警研究更具有实时性,能及时防控岩爆灾害。
(3)本发明提供自动化智能预警系统,具有高效性,同时可以最大限度减少人为因素的干扰。它由知识驱动而非简单由现场施工人员的工程经验主导,在准确性上可以得到保证。
附图说明
图1为全断面隧道掘进机掘进过程中微震监测平台搭建示意图。
图2为本发明的整体流程示意图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
本发明依次进行微震监测平台搭建、微震数据降维、微震数据聚类、前兆微震序列构造和智能识别。
一种应用于TBM隧道岩爆灾害实时智能预警方法,包括以下步骤:
步骤1、全断面隧道掘进机掘进过程中微震监测平台搭建,详见图1。
步骤1.1:选取监测断面。距离开挖掌子面由近到远依次布设有第1组监测断面、第2组监测断面和第3组监测断面,第1组监测断面距离开挖掌子面80m,这是为了保证工作人员和设备的安全。第2组监测断面和第1组监测断面之间的间距、第3组监测断面和第2组监测断面之间的间距均为35m。
步骤1.2:传感器钻孔安装。传感器选用单轴加速度计,其频率响应范围为50Hz-5kHz,灵敏度为30V/g,钻孔深度约为1.5m,每组监测断面在拱肩位置左右对称布置2个传感器。
步骤1.3:全断面隧道掘进机每掘进35m时,将原第3组监测断面前移成为新的第1组监测断面,原第1组监测断面和原第2组监测断面分别作为新的第2组监测断面和新的第3组监测断面。
步骤2、对微震监测平台监测获得的微震数据库降维获得降维后的微震数据矩阵X*,详见图2。
步骤2.1:构建原始微震数据矩阵X=(xij)m×n,m是微震数据库中微震事件的总数量,n是微震事件的震源参数的总数量,xij是第i个微震事件的第j个震源参数的原始微震数据。
步骤2.3:对标准化后的微震数据计算相关系数矩阵R=(rhl)n×n,其中rhl为第h个震源参数和第l个震源参数之间的皮尔逊相关系数。
步骤2.4:计算相关系数矩阵的特征值λ1,λ2,…,λn及其对应的特征向量v1,v2,…,vn。
步骤2.5:将数值大于1且累计贡献率超过80%的特征值对应的特征向量组成矩阵V。特征值的累计贡献率根据公式2计算。
其中,ηg指g个特征值的累计贡献率。
步骤2.6:根据公式3可以得到降维后的微震数据矩阵X*。
X*=X×V 公式3
其中,X*指降维后的微震数据矩阵;X指原始微震数据矩阵;V指步骤2.5中由特征向量组成的矩阵。
步骤3、利用高斯混合模型进行降维后的微震数据矩阵中的微震事件聚类,为降维后的微震数据矩阵中的微震事件分配聚类标签。
步骤3.1:高斯混合模型是多个高斯密度函数的加权和,构造高斯混合模型P(y|θ),如公式4、公式5和公式6。
其中,φk(y|θk)指高斯混合模型中的第k个高斯密度函数;αk(αk≥0)指第k个高斯密度函数的权重系数;指第k个高斯密度函数中的特征参数,且y为微震事件,σk和μk均为高斯密度函数中的待定参数,K(大写)为高斯混合模型中高斯密度函数的数量。
期望最大化算法的流程如下:
步骤3.2.2、根据公式7计算高斯混合模型中的第k个高斯密度函数φk(y|θk)对降维后的微震数据矩阵中第i个微震事件yi的响应γik。
其中,γik指高斯混合模型中的第k个高斯密度函数对降维后的微震数据矩阵中第i个微震事件的响应,m是微震数据库中微震事件的总数量,yi为第i个微震事件。
步骤3.3、迭代反演得到高斯混合模型中的参数α={α1,α2,...,αK}和后,利用高斯混合模型可以对降维后的微震数据矩阵中的微震事件进行聚类,为降维后的微震数据矩阵中的微震事件分配聚类标签,还可以获得不同微震事件属于各个聚类标签的概率。
步骤4、根据降维后的微震数据矩阵组建前兆微震序列数据库。
为降维后的微震数据矩阵中的微震事件分配聚类标签后,将每5个微震事件对应的聚类标签按照微震事件捕捉时间的先后顺序排列即可构造1个前兆微震序列,进而建立由前兆微震序列组成的前兆微震序列数据库。
步骤5、建立分类回归树。
分类回归树由一系列节点组成,包括根节点、内部节点和叶节点。
根节点用于存储步骤4建立的前兆微震序列数据库中所有的前兆微震序列,在一个分类回归树中只有一个根节点。内部节点表示属性条件,当前层内部节点包含上一层内部节点或根节点中满足当前层内部节点表示的属性条件的前兆微震序列。换句话说,内部节点将上一层内部节点或根节点中包含的前兆微震序列分割为满足不同属性条件的前兆微震序列子集。叶节点包含岩爆发生与否以及岩爆烈度等信息。
基尼指数用于确定哪种属性条件作为分割判据,完成分类回归树的生长过程。基尼指数越小,节点纯度越高,分割效果越好。属性条件A的基尼指数Gini_index(D,A)定义为:
其中,D指由待分割节点所包含的前兆微震序列组成的前兆微震序列数据集,待分割节点包括内部节点或者根节点;Dv指由属性条件A分割前兆微震序列数据集得到的前兆微震序列数据子集,V指前兆微震序列数据集D中样本的类数,每个前兆微震序列都会对应岩爆发生与否以及岩爆烈度的信息,将对应相同信息的前兆微震序列称为同一类样本,在这里分为五类:无岩爆,轻微岩爆,中等岩爆,强岩爆,极强岩爆,即在这里V=5。
基尼值Gini(Dv)定义为:
其中,pk指前兆微震序列数据子集Dv中第k类样本所占的比例,每个前兆微震序列都会对应一个岩爆发生与否以及岩爆烈度的信息,将对应相同信息的前兆微震序列称为同一类样本,岩爆烈度分为五类:无岩爆,轻微岩爆,中等岩爆,强岩爆,极强岩爆。
步骤6、训练分类回归树,并通过加权投票法得到自适应增强集成模型的输出。
步骤6.1:初始化前兆微震序列数据库中前兆微震序列的权重。
步骤6.2:根据步骤5利用前兆微震序列数据库建立分类回归树。
步骤6.3:计算分类回归树的预测误差。
步骤6.4:计算分类回归树在自适应增强集成模型(AdaBoost)中的权重。
步骤6.5:根据分类回归树预测结果更新前兆微震序列数据库中前兆微震序列的权重。
步骤6.6:如果已建立的分类回归树数量达到预先设定值,则停止训练分类回归树,并通过加权投票法得到自适应增强集成模型(AdaBoost)的输出;如果已建立的分类回归树数量未达到预先设定值,则重复步骤6.2~6.5。
步骤7、利用贝叶斯优化算法(Bayes)对自适应增强集成模型(AdaBoost)中的超参数寻优,自适应增强集成模型(AdaBoost)中的超参数包括分类回归树数目、最大树深和学习率。
步骤7.1:初始化自适应增强集成模型(AdaBoost)中的超参数。本实施例中,初始化100组自适应增强集成模型(AdaBoost)中的超参数从而建立超参数数据库。
步骤7.2:根据建立的超参数数据库训练高斯过程代理模型,即公式13。
f(B)~GP(0,κ(xp,xq)) 公式13
其中,κ(xp,xq)是关于超参数的协方差函数矩阵;xp是超参数数据库中的第p个超参数;xq是超参数数据库中的第q个超参数;f(B)是自适应增强集成模型(AdaBoost)在不同超参数组合下的精度。GP为高斯过程函数(Gaussian Process)。
步骤7.3:通过公式14所示最大化概率提升获取函数PI(B)(PI,probabilityimprovement)确定下一组待评价超参数,并补充到超参数数据库中。
其中,Φ指标准正态分布的累积分布函数;μ(B)和σ(B)分别指根据步骤7.2建立的高斯过程代理模型得到的自适应增强集成模型(AdaBoost)在不同超参数组合下的精度的平均值和标准差;f(B*)指当前超参数数据库下自适应增强集成模型(AdaBoost)在不同超参数组合下的精度的最优值。
步骤7.4:若达到设定的最大迭代次数则终止迭代,否则重复步骤7.2~7.3。迭代结束后,将超参数数据库中实现最高自适应增强集成模型(AdaBoost)精度的一组超参数作为最优超参数。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,包括以下步骤:
步骤1、全断面隧道掘进机掘进过程中进行微震监测平台搭建;
步骤2、对微震监测平台监测获得的微震数据库降维获得降维后的微震数据矩阵X*;
步骤3、利用高斯混合模型进行降维后的微震数据矩阵中的微震事件聚类,为降维后的微震数据矩阵中的微震事件分配聚类标签;
步骤4、根据降维后的微震数据矩阵组建前兆微震序列数据库;
步骤5、建立分类回归树,分类回归树包括根节点、内部节点和叶节点,
根节点用于存储步骤4建立的前兆微震序列数据库中所有的前兆微震序列,内部节点表示属性条件,内部节点将上一层内部节点或根节点中包含的前兆微震序列分割为满足不同属性条件的前兆微震序列子集,叶节点包含岩爆发生与否以及岩爆烈度;
步骤6、训练分类回归树,并通过加权投票法得到自适应增强集成模型的输出;
步骤7、利用贝叶斯优化算法对自适应增强集成模型中的超参数寻优。
2.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤1包括以下步骤:
步骤1.1、选取监测断,距离开挖掌子面由近到远依次布设有第1组监测断面、第2组监测断面和第3组监测断面,
步骤1.2、每组监测断面在拱肩位置左右对称布置2个传感器;
步骤1.3、全断面隧道掘进机向前掘进设定距离后,将原第3组监测断面前移成为新的第1组监测断面,原第1组监测断面和原第2组监测断面分别作为新的第2组监测断面和新的第3组监测断面。
3.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤2包括以下步骤:
步骤2.1、构建原始微震数据矩阵X=(xij)m×n,m是微震数据库中微震事件的总数量,n是微震事件的震源参数的总数量,xij是第i个微震事件的第j个震源参数的原始微震数据,
步骤2.3:对标准化后的微震数据计算相关系数矩阵R=(rhl)n×n,其中rhl为第h个震源参数和第l个震源参数之间的皮尔逊相关系数,
步骤2.4:计算相关系数矩阵的特征值λ1,λ2,…,λn及其对应的特征向量v1,v2,…,vn,
步骤2.5:将数值大于1且累计贡献率超过80%的特征值对应的特征向量组成矩阵V,特征值的累计贡献率根据公式2计算,
其中,ηg指g个特征值的累计贡献率,
步骤2.6:根据公式3得到降维后的微震数据矩阵X*,
X*=X×V 公式3
其中,X*指降维后的微震数据矩阵;X指原始微震数据矩阵;V指步骤2.5中由特征向量组成的矩阵。
4.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤3包括以下步骤:
步骤3.1、按照公式4、公式5和公式6构造高斯混合模型P(y|θ),
其中,φk(y|θk)指高斯混合模型中的第k个高斯密度函数;αk指第k个高斯密度函数的权重系数,αk≥0;指第k个高斯密度函数中的特征参数,且y为微震事件,σk和μk均为高斯密度函数中的待定参数,K为高斯混合模型中高斯密度函数的数量,
期望最大化算法的流程如下:
步骤3.2.2、根据公式7计算高斯混合模型中的第k个高斯密度函数φk(y|θk)对降维后的微震数据矩阵中第i个微震事件yi的响应γik,
其中,γik指高斯混合模型中的第k个高斯密度函数对降维后的微震数据矩阵中第i个微震事件的响应,m是微震数据库中微震事件的总数量,yi为第i个微震事件,
5.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤4包括以下步骤:
为降维后的微震数据矩阵中的微震事件分配聚类标签后,将每5个微震事件对应的聚类标签按照微震事件捕捉时间的先后顺序排列即构造1个前兆微震序列,进而建立由前兆微震序列组成的前兆微震序列数据库。
7.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤6包括以下步骤:
步骤6.1、初始化前兆微震序列数据库中前兆微震序列的权重,
步骤6.2、根据步骤5利用前兆微震序列数据库建立的分类回归树,
步骤6.3、计算分类回归树的预测误差,
步骤6.4、计算分类回归树在自适应增强集成模型中的权重,
步骤6.5、根据分类回归树预测结果更新前兆微震序列数据库中前兆微震序列的权重,
步骤6.6、如果已建立的分类回归树数量达到预先设定值,则停止训练分类回归树,并通过加权投票法得到自适应增强集成模型的输出;如果已建立的分类回归树数量未达到预先设定值,则重复步骤6.2~6.5。
8.根据权利要求1所述的一种应用于TBM隧道岩爆灾害实时智能预警方法,其特征在于,所述的步骤7包括以下步骤:
步骤7.1、初始化自适应增强集成模型中的超参数,建立超参数数据库,
步骤7.2、根据建立的超参数数据库训练如公式13所示高斯过程代理模型:
f(B)~GP(0,κ(xp,xq)) 公式13
其中,κ(xp,xq)是关于超参数的协方差函数矩阵;xp是超参数数据库中的第p个超参数;xq是超参数数据库中的第q个超参数;f(B)是自适应增强集成模型在不同超参数组合下的精度,GP为高斯过程函数,
步骤7.3、通过公式14所示的最大化概率提升获取函数PI(B)确定下一组待评价超参数,并补充到超参数数据库中,
其中,Φ指标准正态分布的累积分布函数;μ(B)和σ(B)分别指根据步骤7.2建立的高斯过程代理模型得到的自适应增强集成模型在不同超参数组合下的精度的平均值和标准差;f(B*)指当前超参数数据库下自适应增强集成模型在不同超参数组合下的精度的最优值,
步骤7.4、若达到设定的最大迭代次数则终止迭代,否则重复步骤7.2~7.3,迭代结束后,将超参数数据库中实现最高自适应增强集成模型精度的一组超参数作为最优超参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110432239.6A CN113156492B (zh) | 2021-04-21 | 2021-04-21 | 一种应用于tbm隧道岩爆灾害实时智能预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110432239.6A CN113156492B (zh) | 2021-04-21 | 2021-04-21 | 一种应用于tbm隧道岩爆灾害实时智能预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113156492A CN113156492A (zh) | 2021-07-23 |
CN113156492B true CN113156492B (zh) | 2022-07-26 |
Family
ID=76867864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110432239.6A Active CN113156492B (zh) | 2021-04-21 | 2021-04-21 | 一种应用于tbm隧道岩爆灾害实时智能预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113156492B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113936104A (zh) * | 2021-12-17 | 2022-01-14 | 矿冶科技集团有限公司 | 裂纹识别方法、装置、电子设备及介质 |
CN116297868B (zh) * | 2023-05-22 | 2023-08-08 | 中国矿业大学(北京) | 一种岩爆声发射前兆特征的计算方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697999B (zh) * | 2013-12-30 | 2015-10-14 | 中国科学院武汉岩土力学研究所 | 一种高应力硬岩tbm施工隧道微震波速实时获取方法 |
CN103941282B (zh) * | 2014-05-07 | 2017-02-01 | 中国科学院武汉岩土力学研究所 | 一种深埋大断面隧道下台阶开挖岩爆风险区判别方法 |
CN106054243A (zh) * | 2016-06-03 | 2016-10-26 | 河南理工大学 | 一种基于微震监测的岩爆多指标预测方法 |
CN109343111B (zh) * | 2018-09-30 | 2021-09-07 | 中国科学院武汉岩土力学研究所 | 一种软硬互层地质条件下长距离tbm隧洞岩爆微震监测方法 |
CN110298503B (zh) * | 2019-06-26 | 2022-11-25 | 东北大学 | 基于微震信息和深度卷积神经网络的隧洞岩爆预警方法 |
CN111123355B (zh) * | 2020-01-07 | 2020-12-11 | 山东大学 | 一种基于微震监测数据的岩爆预测方法及系统 |
-
2021
- 2021-04-21 CN CN202110432239.6A patent/CN113156492B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113156492A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112183993B (zh) | 一种tbm掘进控制参数智能预测及优化决策方法 | |
CN113156492B (zh) | 一种应用于tbm隧道岩爆灾害实时智能预警方法 | |
CN112116162B (zh) | 基于ceemdan-qfoa-lstm的输电线覆冰厚度预测方法 | |
CN113610945B (zh) | 一种基于混合神经网络的地应力曲线预测方法 | |
CN112949202A (zh) | 一种基于贝叶斯网络的岩爆概率预测方法 | |
CN113408190B (zh) | 基于Bayes-LSTM模型的公路隧道施工期围岩变形预测方法 | |
CN113537566B (zh) | 一种基于dccso优化深度学习模型的超短期风电功率预测方法 | |
CN114723095A (zh) | 缺失测井曲线预测方法及装置 | |
CN112712025A (zh) | 一种基于长短期记忆神经网络的复杂岩性识别方法 | |
CN115983465A (zh) | 一种基于小样本学习的冲击地压时序预测模型构建方法 | |
CN116384554A (zh) | 机械钻速预测方法、装置、电子设备及计算机存储介质 | |
CN112926267A (zh) | 一种基于掘进参数反演的tbm隧道岩爆等级预测方法及系统 | |
CN115481565A (zh) | 一种基于lstm和蚁群算法的土压平衡盾构机掘进参数预测方法 | |
CN117252236A (zh) | 一种基于dbn-lstm-bwoa的隧道爆破振动峰值预测方法、系统、设备及介质 | |
CN115526096A (zh) | 盾构机施工掘进参数的预测方法、装置及可存储介质 | |
CN117150875A (zh) | 基于深度学习的钻前测井曲线预测方法 | |
CN112149255A (zh) | 基于lstm神经网络的采煤工作面液压支架工作阻力预测方法 | |
Zhang et al. | Performance and Applicability of Recognizing Microseismic Waveforms Using Neural Networks in Tunnels | |
CN115964667A (zh) | 基于深度学习和重采样的河流-湖泊岩相测井识别方法 | |
CN115983333A (zh) | 基于混合神经网络的盾构掘进速度智能预测方法 | |
CN115906926A (zh) | 一种盾构法隧道掘进速度智能控制系统 | |
CN113435055B (zh) | 盾构刀盘扭矩领域自适应迁移预测方法和系统 | |
Zhang et al. | Time Series Prediction on Settlement of Metro Tunnels Adjacent to Deep Foundation Pit by Clustering Monitoring Data | |
CN113361476B (zh) | 一种基于人工智能技术的张衡一号震前异常信号识别方法 | |
Wang et al. | Comparative study on prediction models of EPB shield tunnelling parameters in the water-rich round-gravel formation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |