CN111079342A - 一种基于在线岩体等级分类的tbm掘进性能预测方法 - Google Patents
一种基于在线岩体等级分类的tbm掘进性能预测方法 Download PDFInfo
- Publication number
- CN111079342A CN111079342A CN201911198511.8A CN201911198511A CN111079342A CN 111079342 A CN111079342 A CN 111079342A CN 201911198511 A CN201911198511 A CN 201911198511A CN 111079342 A CN111079342 A CN 111079342A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- tbm
- state parameters
- tunneling
- tbm tunneling
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 136
- 230000005641 tunneling Effects 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 32
- 239000011159 matrix material Substances 0.000 claims description 14
- 230000001174 ascending effect Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 238000012549 training Methods 0.000 claims description 10
- 230000035515 penetration Effects 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000013473 artificial intelligence Methods 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 238000000611 regression analysis Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000013527 convolutional neural network Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 230000007787 long-term memory Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000008447 perception Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000006403 short-term memory Effects 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 4
- 238000010276 construction Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 7
- 230000001141 propulsive effect Effects 0.000 description 5
- 230000000630 rising effect Effects 0.000 description 4
- 206010057175 Mass conditions Diseases 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000010219 correlation analysis Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000012843 least square support vector machine Methods 0.000 description 2
- 208000022971 Tuberculous meningitis Diseases 0.000 description 1
- 238000012098 association analyses Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 208000001223 meningeal tuberculosis Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Biophysics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于在线岩体等级分类的TBM掘进性能预测方法,包括如下步骤:采集TBM掘进状态参数及在掘掌子面的岩体状态参数;对获取的TBM掘进状态参数进行数据预处理;获得由TBM掘进状态参数预估岩体状态参数的岩体感知模型;采用聚类分析法对岩体状态参数进行岩体等级分类;构建与岩体等级相匹配的TBM性能预测模型;采集新的TBM掘进状态参数,执行步骤S2并根据步骤S3中的岩体感知模型计算新的TBM掘进状态参数所对应的预估岩体状态参数,执行步骤S4确定预估岩体状态参数所对应的岩体等级,进而与TBM性能预测模型相匹配。本发明可随岩体等级变化来自适应匹配TBM性能预测模型,准确预测当前TBM掘进状态。
Description
技术领域
本发明属于隧道工程TBM施工技术领域,具体涉及一种基于在线岩体等级分类的TBM掘进性能预测方法。
背景技术
全断面隧道掘进机(Tunnel Boring Machine,简称TBM),因其具有掘进速度快、安全性高、环保、经济、劳动强度低等优点,目前已成为长大隧道建设的首选施工方法。随着TBM在我国的快速发展和广泛应用,TBM的安全高效掘进与智能控制变得越来越重要,它对于保障TBM掘进效率、避免重大经济损失和工期延误具有至关重要的现实意义。
在实际的隧道开挖中,岩体质量等级在整个掘进过程呈现复杂多变的特点。受岩体等级动态多变的影响,TBM掘进性能的变化趋势在不同围岩条件下呈现不确定性。然而,传统TBM性能预测方法较少考虑TBM掘进参数受岩体变化的影响,这给TBM性能的精准预测带来重大挑战。因而,如何更好地理解岩机相互作用并为硬岩TBMs的掘进参数预测开发一个精度更高、鲁棒性更好的性能预测模型是当前研究的一个热点。
发明内容
针对传统TBM性能预测方法在岩体等级动态多变条件下精度低、鲁棒性差的问题,本发明提出了一种基于在线岩体等级分类的TBM掘进性能预测方法,本发明可根据在线辨识的岩体等级与掘进状态参数来自适应匹配TBM性能预测模型,从而方便司机对岩体条件快速做出响应,并及时调整设备控制参数。
为解决以上技术问题,本发明所采用的技术方案如下:
一种基于在线岩体等级分类的TBM掘进性能预测方法,包括如下步骤:
S1,采集TBM掘进状态参数及在掘掌子面的岩体状态参数,构建岩机相互作用数据库;
S2,对获取的TBM掘进状态参数进行数据预处理,以去除其中的背景噪声,提升数据质量;
S3,对预处理后的TBM掘进状态参数与岩体状态参数进行回归分析,构建拟合方程,获得由TBM掘进状态参数预估岩体状态参数的岩体感知模型;
S4,采用聚类分析法对岩体状态参数进行围岩等级分类,获取不同岩体等级下的岩机数据;
S5,利用人工智能算法对不同岩体等级下岩体状态参数与TBM掘进状态参数之间的复杂非线性关系进行建模,构建与岩体等级相匹配的TBM性能预测模型;
S6,采集新的TBM掘进状态参数,执行步骤S2并根据步骤S3中的岩体感知模型计算新的TBM掘进状态参数所对应的预估岩体状态参数,执行步骤S4确定预估岩体状态参数所对应的岩体等级,进而与所述TBM性能预测模型相匹配。
在步骤S1中,所述TBM掘进状态参数包括刀盘转速n、推进速度v、刀盘贯入度p、刀盘扭矩T和总推力F;所述岩体状态参数包括岩体节理Jv和单轴抗压强度UCS。
在步骤S2中,所述数据预处理包括如下步骤:
a,剔除TBM掘进状态参数中的停机段数据;
b,利用箱线法对剔除停机段数据后的TBM掘进状态参数进行异常值剔除;
c,利用最大类间差方法自适应选取掘进循环阈值,实现循环上升段及稳定段数据自动划分;
d,对划分后的上升段和稳定段数据进行降噪、增强处理,所述降噪方法包括但不限于小波变换(Wavelet Transform)或经验模态分解(Empirical Mode Decomposition)等算法。
在步骤S3中,所述回归分析的方法包括但不限于最小二乘分步回归、BP神经网络或支持向量回归。
在步骤S3中,所述岩体感知模型的公式为:
Jvpre=23.79-0.078×Fn+1.27×p-61.6×T/104+149.4×T2/108-0.43×n;
式中,Jvpre表示预估岩体节理,Fn表示单刀正应力,p表示刀盘贯入度,T表示刀盘扭矩,n表示刀盘转速;
UCSpre=39.5+0.27×Fn-0.71×p+55×T/104-201.75×T2/108-0.43×n;
式中,UCSpre表示预估单轴抗压强度。
在步骤S4中,所述聚类分析法包括但不限于k-means、层次聚类、高斯混合模型或基于密度的聚类算法。
在步骤S5中,所述人工智能算法包括但不限于最小二乘支持向量机、DS证据理论融合、卷积神经网络(Convolutional Neural Networks)或长短时记忆神经网络(LongShort Term Memory Network)等算法。
所述TBM性能预测模型的公式为:
式中,i=1,2,…,m,αi表示第i个训练样本的拉格朗日乘子,b表示偏置,k(x,xi)表示岩机特征参数矩阵x和第i个训练样本xi通过核函数映射所得计算结果,y(x)表示TBM性能预测模型所需预测的总推进力F或者刀盘扭矩T。
所述TBM性能预测模型的公式y(x)根据以下步骤可得:
目标函数可表述为:
式中,i=1,2,…,m,ξ为误差量,γ为惩罚系数且满足γ>0;
引入Lagrange乘数对目标函数进行最大似然估计,化简如下:
式中,α为拉格朗日乘子;
分别对w,b,ξ,α进行求导,可得:
根据上述四个条件得到关于α和b的线性方程组:
根据Mercer条件,核函数矩阵Ω可表示为:
解上述线性方程组,可得TBM性能预测模型:
式中,k(x,xi)表示岩机特征参数矩阵x和第i个训练样本xi通过核函数映射所得计算结果,常用核函数包括高斯核函数、正弦核函数以及Sigmoid核函数等,y(x)表示TBM性能预测模型所需预测的总推进力F或者刀盘扭矩T。
本发明的有益效果:
本发明集岩机数据获取、数据预处理、岩体感知、岩体等级分类及掘进状态预测于一体,可根据在线辨识的岩体等级与掘进运行状态来自适应匹配TBM性能预测模型,从而为主司机对岩体条件快速做出响应、及时调整设备控制参数提供建议与支持,具有实时性和准确性高、鲁棒性强的显著优点,且不影响TBM的正常掘进,从而大大提高设备施工效率、减少因工期延误所导致的经济损失。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种基于在线岩体等级分类的TBM掘进性能预测方法,如图1所示,包括如下步骤:
S1:岩机数据获取:采集反映设备运行状态的TBM掘进状态参数及在掘掌子面岩体状态参数,构建岩机相互作用数据库。
所述TBM掘进状态参数包括刀盘转速n、推进速度v、刀盘贯入度p、刀盘扭矩T和总推力F;所述岩体状态参数包括表示岩体完整性的岩体节理Jv和表示岩体强度特性的单轴抗压强度UCS指标。
S2:数据预处理:对获取的TBM掘进状态参数进行异常值剔除、循环划分、降噪等预处理,以去除其中的背景噪声、降低误报警概率,提升数据质量。
所述数据预处理包括但不限于如下步骤:
a.按照如下规则剔除TBM掘进状态参数中的停机段数据,提取施工段数据;
P=f(n)f(T)f(F)f(v);
函数f(x)定义如下:
b.利用箱线法对剔除停机段数据后的TBM掘进状态参数进行异常值剔除;
下面以刀盘扭矩为例阐述异常值剔除的实施过程:
1)将施工段刀盘扭矩数据按照从小到大的顺序排列,全部数据被分成相等的四部分,其中每部分包括25%的数据,处在各分位点的数值就是四分位数;
2)选取第一、第三四分位数Q1、Q3,Q1和Q3分别等于总体全部数据由小到大排列后第25%、75%的数值;
3)计算四分位距IQR=Q3-Q1,并按照如下原则获取箱线法上、下限阈值;
下限阈值:L1=Q1-IQR*1.5;
上限阈值:U1=Q3+IQR*1.5;
4)通过判断施工段刀盘扭矩数据是否超出[L1,U1]范围来剔除异常值,如果超出范围,则剔除对应刀盘扭矩数据,否则,视为正常数据。
c.利用最大类间差方法自适应选取掘进循环阈值,实现循环上升段、稳定段数据的自动划分;
本实施例主要通过对每个循环总推进力的阈值进行自适应计算,以自动划分循环上升段与稳定段数据,具体实施过程如下;
1)初始化循环总推进力中上升段与稳定段间阈值g0,并统计循环中总推进力小于g0的数据个数记作N0,总推进力大于g0的数据个数记作N1;
2)将循环总推进力中上升段数据占整个循环数据的比例记为w0,并计算循环上升段总推进力的平均值μ0,若总推进力F<g0,则视为上升段;将循环总推进力中稳定段数据占整个循环数据的比例记为w1,并计算循环稳定段总推进力的平均值μ1,若总推进力F>=g0,则视为稳定段;
3)根据循环总推进力中上升段、稳定段数据所占比例w0、w1及其平均值μ0、μ1,可计算整个循环总推进力的均值μ及类间方差S:
w0=N0/(N0+N1)
w1=N1/(N0+N1)
w0+w1=1
μ=w0*μ0+w1*μ1
S=w0(μ0-μ)2+w1(μ1-μ)2
通过进一步化简,可得到等价公式:
S=w0 w1(μ0-μ1)2
4)将整个循环总推进力的均值μ作为新的循环总推进力上升段与稳定段间阈值g1,再按照步骤1)~步骤3)计算新的循环总推进力的均值μˊ及类间方差S′,通过对比类间方差S、S′,如果S′>S,则将循环总推进力的均值μˊ作为循环总推进力上升段与稳定段间最终阈值,此时循环上升段、稳定段数据划分结束;否则,继续执行步骤1)~步骤3),直至循环划分结束。
d.对划分后上升段、稳定段数据进行降噪、增强处理,降噪方法包括但不限于小波变换或经验模态分解算法。
S3:岩体感知模型构建:对预处理后的TBM掘进状态参数与岩体状态参数进行回归分析,构建拟合方程,实现由TBM掘进状态参数估计岩体状态参数的岩体感知模型。
所述回归分析的方法包括但不限于最小二乘分步回归、BP神经网络或支持向量回归等算法。
本实施例以最小二乘分步回归为例阐述岩体感知模型的构建过程,具体实施步骤如下:
1)参数关联分析:通过对岩体状态参数与TBM掘进状态参数进行关联分析,可得:岩体节理Jv、单轴抗压强度UCS均与单刀正应力Fn、刀盘贯入度p、刀盘转速n均呈线性关系,岩体节理Jv、单轴抗压强度UCS均与刀盘扭矩T呈非线性的二阶多项式关系。
2)岩体感知模型构建:根据参数关联分析结果,本发明实施例所构建岩体状态参数X与TBM掘进状态参数矩阵A之间的回归模型如下:
X=A*B;
其中,
式中,岩体状态参数回归系数向量B=[b1,b2,…,b6],掘进状态参数矩阵其中Jv为岩体节理,单位为Num/m3,UCS为岩体单轴抗压强度,单位为MPa,Fn为单刀正应力,Ff表示TBM掘进中所受摩擦阻力,单位为kN,N表示刀盘滚刀的数量,p为贯入度,单位为mm/r,n为刀盘转速,单位为r/min,T为刀盘扭矩,单位为kN·m,m为样本数量。
由于上述回归方程是一个超定方程组,利用最小二乘法来确定回归系数向量B,以使得损失函数H最小:
H=argmin||A*B-X||;
当H最小时,B=(AT*A)-1*AT*X,从而求得最小均方误差意义下的回归系数。
当上述回归系数向量B经过最小二乘法确定后,在掘岩体状态参数与TBM掘进状态参数之间的岩体感知模型构建如下:
Jvpre=23.79-0.078×Fn+1.27×p-61.6×T/104+149.4×T2/108-0.43×n;
UCSpre=39.5+0.27×Fn-0.71×p+55×T/104-201.75×T2/108-0.43×n。
S4:岩体等级分类:采用聚类分析法对岩体状态参数进行岩体等级分类,获取不同岩体等级下的岩机数据。
所述聚类分析法包括但不限于k-means、层次聚类、高斯混合模型或基于密度的聚类算法等。
本实施例以k-means为例阐述岩体等级分类的实施过程,具体实施步骤如下:
1)初始化分类:初始化K个聚类中心,以欧氏距离来表征样本间相似性的度量,计算所有样本到各个初始聚类中心Ck(k=1,2,…,K)的距离J(Ck),按照距离最小的判别准则将样本Xi=[Jvi,UCSi]划分为某类Ck;
2)更新聚类中心:计算每个类别{C1,C2,…,CK}中隶属样本的均值{u1,u2,…,uK}作为该类别的新聚类中心C'={C'1,C'2,…,C'K},而后重新获取所有样本到其所在类别聚类中心的距离平方和J(C');
式中,Xi表示第i个岩体状态参数样本,且Xi=[Jvi,UCSi],uk表示第k类聚类中心Ck中隶属样本的均值,且k=1,2,…,K;
3)收敛规则:判断前后两次聚类分析时聚类中心C'及J(C')值是否变化,如果不再变化,则聚类结束,否则重复执行步骤2更新聚类中心,直至聚类结束。
4)根据步骤1)至步骤3),本实施例利用k-means聚类方法对TBM掘进状态参数所预估岩体状态参数进行岩体等级分类,并将岩体划分为5个聚类类别{C1,C2,C3,C4,C5},即5种岩体等级{Ⅰ,Ⅱ,Ⅲ,Ⅳ,Ⅴ}。根据每个聚类类别的隶属样本,分别提取对应岩体等级下的上升段和稳定段岩机数据。
S5:TBM性能预测模型构建:利用人工智能算法对不同岩体等级下岩体状态参数与TBM掘进状态参数之间的复杂非线性关系进行建模,构建与岩体等级相匹配的TBM性能预测模型。
所述性能预测模型包括但不限于最小二乘支持向量机、DS证据理论融合、卷积神经网络、长短时记忆神经网络等算法。
本实施例以最小二乘支持向量机为例阐述TBM性能预测模型的构建过程,具体实施步骤如下:
a.构建TBM性能预测模型的训练样本集;
根据不同岩体等级下的上升段和稳定段岩机数据构建TBM性能预测模型的训练样本集其中x表示第k种聚类类别Ck(k=1,2,3,4,5)下岩机特征参数矩阵,即对应岩体的具体岩体等级,且x=[UCS,Jv,n,v],m表示训练样本的数量,y表示TBM掘进状态参数中的稳定段刀盘扭矩T或者总推力F。
b.建立TBM性能预测模型;
为了有效表征不同岩体等级下岩机数据与TBM掘进状态参数之间的复杂非线性关系,应用基于最小二乘支持向量回归-LSSVR的人工智能方法,建模过程如下:
LSSVR目标函数可表述为:
式中,i=1,2,…,m,ξ为误差量,γ为惩罚系数且满足γ>0。
引入Lagrange乘数对LSSVR目标函数进行最大似然估计,化简如下:
式中,α为拉格朗日乘子。
分别对w,b,ξ,α进行求导,可得:
根据上述四个条件得到关于α和b的线性方程组:
根据Mercer条件,核函数矩阵Ω可表示为:
解上述线性方程组,可得基于LSSVR的TBM性能预测模型:
式中,k(x,xi)表示岩机特征参数矩阵x和第i个训练样本xi通过核函数映射所得计算结果,常用核函数包括高斯核函数、正弦核函数以及Sigmoid核函数等,y(x)表示TBM性能预测模型所需预测的总推进力F或者刀盘扭矩T。
根据步骤b,分别对5种聚类类别下岩机状态参数与TBM掘进状态参数之间的复杂非线性关系进行建模,可得到5个与岩体等级相匹配的TBM性能预测模型。
S6:TBM性能预测模型验证:采集新的TBM掘进状态参数,执行步骤S2并根据步骤S3中的岩体感知模型计算新的TBM掘进状态参数所对应的预估岩体状态参数,执行步骤S4确定预估岩体状态参数所对应的岩体等级,进而与所述TBM性能预测模型相匹配。
上述系统可随岩体条件变化来自适应匹配TBM性能预测模型,并根据TBM性能预测模型预测结果及时为司机提供设备参数的调整建议,从而快速对当前岩体条件做出响应。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,包括如下步骤:
S1,采集TBM掘进状态参数及在掘掌子面的岩体状态参数;
S2,对获取的TBM掘进状态参数进行数据预处理;
S3,对预处理后的TBM掘进状态参数与岩体状态参数进行回归分析,构建拟合方程,获得由TBM掘进状态参数预估岩体状态参数的岩体感知模型;
S4,采用聚类分析法对岩体状态参数进行岩体等级分类;
S5,利用人工智能算法对不同岩体等级下岩体状态参数与TBM掘进状态参数之间的非线性关系进行建模,构建与岩体等级相匹配的TBM性能预测模型;
S6,采集新的TBM掘进状态参数,执行步骤S2并根据步骤S3中的岩体感知模型计算新的TBM掘进状态参数所对应的预估岩体状态参数,执行步骤S4确定预估岩体状态参数所对应的岩体等级,进而与所述TBM性能预测模型相匹配。
2.根据权利要求1所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S1中,所述TBM掘进状态参数包括刀盘转速n、推进速度v、刀盘贯入度p、刀盘扭矩T和总推力F;所述岩体状态参数包括岩体节理Jv和单轴抗压强度UCS。
3.根据权利要求1或2所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S2中,所述数据预处理包括如下步骤:
a,剔除TBM掘进状态参数中的停机段数据;
b,对剔除停机段数据后的TBM掘进状态参数进行异常值剔除;
c,选取掘进循环阈值,划分出循环上升段及稳定段数据;
d,对划分后的上升段和稳定段数据进行降噪、增强处理。
4.根据权利要求3所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤d中,所述降噪方法包括但不限于小波变换或经验模态分解。
5.根据权利要求4所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S3中,所述回归分析的方法包括但不限于最小二乘分步回归、BP神经网络或支持向量回归。
6.根据权利要求1或5所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S3中,所述岩体感知模型的公式为:
Jvpre=23.79-0.078×Fn+1.27×p-61.6×T/104+149.4×T2/108-0.43×n;
式中,Jvpre表示预估岩体节理,Fn表示单刀正应力,p表示刀盘贯入度,T表示刀盘扭矩,n表示刀盘转速;
UCSpre=39.5+0.27×Fn-0.71×p+55×T/104-201.75×T2/108-0.43×n;
式中,UCSpre表示预估单轴抗压强度。
7.根据权利要求6所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S4中,所述聚类分析法包括但不限于k-means、层次聚类、高斯混合模型或基于密度的聚类算法。
8.根据权利要求1或7所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,在步骤S5中,所述人工智能算法包括但不限于最小二乘支持向量机、DS证据理论融合、卷积神经网络或长短时记忆神经网络。
10.根据权利要求9所述的基于在线岩体等级分类的TBM掘进性能预测方法,其特征在于,所述TBM性能预测模型的公式y(x)根据以下步骤可得:
目标函数可表述为:
式中,i=1,2,…,m,ξ为误差量,γ为惩罚系数且满足γ>0;
引入Lagrange乘数对目标函数进行最大似然估计,化简如下:
式中,α为拉格朗日乘子;
分别对w,b,ξ,α进行求导,可得:
根据上述四个条件得到关于α和b的线性方程组:
根据Mercer条件,核函数矩阵Ω可表示为:
解上述线性方程组,可得TBM性能预测模型:
式中,k(x,xi)表示岩机特征参数矩阵x和第i个训练样本xi通过核函数映射所得计算结果,常用核函数包括高斯核函数、正弦核函数以及Sigmoid核函数,y(x)表示TBM性能预测模型所需预测的总推进力F或者刀盘扭矩T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911198511.8A CN111079342B (zh) | 2019-11-29 | 2019-11-29 | 一种基于在线岩体等级分类的tbm掘进性能预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911198511.8A CN111079342B (zh) | 2019-11-29 | 2019-11-29 | 一种基于在线岩体等级分类的tbm掘进性能预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111079342A true CN111079342A (zh) | 2020-04-28 |
CN111079342B CN111079342B (zh) | 2023-04-07 |
Family
ID=70312392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911198511.8A Active CN111079342B (zh) | 2019-11-29 | 2019-11-29 | 一种基于在线岩体等级分类的tbm掘进性能预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111079342B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111594198A (zh) * | 2020-04-30 | 2020-08-28 | 东北大学 | 一种开敞式tbm施工数据掘进循环的划分方法 |
CN111594197A (zh) * | 2020-04-30 | 2020-08-28 | 东北大学 | 一种基于tbm岩机信息的隧洞断层段超前智能感知方法 |
CN112065421A (zh) * | 2020-10-10 | 2020-12-11 | 中国铁建重工集团股份有限公司 | 一种掘进机刀盘自动定位方法 |
CN112115537A (zh) * | 2020-09-17 | 2020-12-22 | 中铁工程装备集团有限公司 | 一种tbm可掘性分级方法 |
CN112196559A (zh) * | 2020-09-30 | 2021-01-08 | 山东大学 | 基于掘进速度与刀具消耗最优的tbm操作参数优化方法 |
CN112598061A (zh) * | 2020-12-23 | 2021-04-02 | 中铁工程装备集团有限公司 | 一种隧道围岩聚类分级方法 |
CN113158561A (zh) * | 2021-04-14 | 2021-07-23 | 山东大学 | 适用于多种岩体条件的tbm操作参数优化方法及系统 |
CN113420506A (zh) * | 2021-06-30 | 2021-09-21 | 北京交通大学 | 掘进速度的预测模型建立方法、预测方法及装置 |
CN113657515A (zh) * | 2021-08-20 | 2021-11-16 | 盾构及掘进技术国家重点实验室 | 一种基于岩机敏感参量判识和改进fmc模型隧道围岩等级的分类方法 |
CN113779679A (zh) * | 2021-09-14 | 2021-12-10 | 山东高速工程建设集团有限公司 | 一种硬岩盾构或tbm技术参数选型方法及系统 |
CN115142160A (zh) * | 2022-08-22 | 2022-10-04 | 无锡物联网创新中心有限公司 | 一种纱线强力弱环的辨识方法及相关装置 |
CN115438402A (zh) * | 2022-08-22 | 2022-12-06 | 武汉大学 | 基于岩-机互馈感知关系的tbm掘进岩体条件实时预测方法及装置 |
CN115526036A (zh) * | 2022-09-19 | 2022-12-27 | 长安大学 | 一种岩爆规模等级的判别方法及系统 |
CN115618222A (zh) * | 2022-06-21 | 2023-01-17 | 北京交通大学 | 一种隧道掘进响应参数的预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106089222A (zh) * | 2016-06-27 | 2016-11-09 | 中交公局第三工程有限公司 | 一种用于砂岩泥岩单护盾tbm掘进方法 |
CN107201903A (zh) * | 2017-04-07 | 2017-09-26 | 北京工业大学 | Tbm隧道施工的智能化控制方法及系统 |
WO2019042483A2 (zh) * | 2017-08-30 | 2019-03-07 | 中铁工程装备集团有限公司 | 一种tbm在掘岩体状态实时感知系统和方法 |
CN109580916A (zh) * | 2019-01-25 | 2019-04-05 | 黄河勘测规划设计有限公司 | 用于tbm隧洞岩体分级所需参数的快速综合采集系统 |
CN109725129A (zh) * | 2019-01-25 | 2019-05-07 | 黄河勘测规划设计有限公司 | Tbm隧洞岩体分级方法 |
-
2019
- 2019-11-29 CN CN201911198511.8A patent/CN111079342B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106089222A (zh) * | 2016-06-27 | 2016-11-09 | 中交公局第三工程有限公司 | 一种用于砂岩泥岩单护盾tbm掘进方法 |
CN107201903A (zh) * | 2017-04-07 | 2017-09-26 | 北京工业大学 | Tbm隧道施工的智能化控制方法及系统 |
WO2019042483A2 (zh) * | 2017-08-30 | 2019-03-07 | 中铁工程装备集团有限公司 | 一种tbm在掘岩体状态实时感知系统和方法 |
CN109580916A (zh) * | 2019-01-25 | 2019-04-05 | 黄河勘测规划设计有限公司 | 用于tbm隧洞岩体分级所需参数的快速综合采集系统 |
CN109725129A (zh) * | 2019-01-25 | 2019-05-07 | 黄河勘测规划设计有限公司 | Tbm隧洞岩体分级方法 |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111594198A (zh) * | 2020-04-30 | 2020-08-28 | 东北大学 | 一种开敞式tbm施工数据掘进循环的划分方法 |
CN111594197A (zh) * | 2020-04-30 | 2020-08-28 | 东北大学 | 一种基于tbm岩机信息的隧洞断层段超前智能感知方法 |
CN111594198B (zh) * | 2020-04-30 | 2021-03-30 | 东北大学 | 一种开敞式tbm施工数据掘进循环的划分方法 |
CN112115537A (zh) * | 2020-09-17 | 2020-12-22 | 中铁工程装备集团有限公司 | 一种tbm可掘性分级方法 |
CN112115537B (zh) * | 2020-09-17 | 2023-11-17 | 中铁工程装备集团有限公司 | 一种tbm可掘性分级方法 |
CN112196559A (zh) * | 2020-09-30 | 2021-01-08 | 山东大学 | 基于掘进速度与刀具消耗最优的tbm操作参数优化方法 |
CN112065421A (zh) * | 2020-10-10 | 2020-12-11 | 中国铁建重工集团股份有限公司 | 一种掘进机刀盘自动定位方法 |
CN112065421B (zh) * | 2020-10-10 | 2022-04-22 | 中国铁建重工集团股份有限公司 | 一种掘进机刀盘自动定位方法 |
CN112598061A (zh) * | 2020-12-23 | 2021-04-02 | 中铁工程装备集团有限公司 | 一种隧道围岩聚类分级方法 |
CN113158561A (zh) * | 2021-04-14 | 2021-07-23 | 山东大学 | 适用于多种岩体条件的tbm操作参数优化方法及系统 |
CN113420506A (zh) * | 2021-06-30 | 2021-09-21 | 北京交通大学 | 掘进速度的预测模型建立方法、预测方法及装置 |
CN113420506B (zh) * | 2021-06-30 | 2024-08-09 | 北京交通大学 | 掘进速度的预测模型建立方法、预测方法及装置 |
CN113657515A (zh) * | 2021-08-20 | 2021-11-16 | 盾构及掘进技术国家重点实验室 | 一种基于岩机敏感参量判识和改进fmc模型隧道围岩等级的分类方法 |
CN113657515B (zh) * | 2021-08-20 | 2024-06-14 | 盾构及掘进技术国家重点实验室 | 一种基于岩机敏感参量判识和改进fmc模型隧道围岩等级的分类方法 |
CN113779679A (zh) * | 2021-09-14 | 2021-12-10 | 山东高速工程建设集团有限公司 | 一种硬岩盾构或tbm技术参数选型方法及系统 |
CN113779679B (zh) * | 2021-09-14 | 2022-05-24 | 山东高速工程建设集团有限公司 | 一种硬岩盾构或tbm技术参数选型方法及系统 |
CN115618222A (zh) * | 2022-06-21 | 2023-01-17 | 北京交通大学 | 一种隧道掘进响应参数的预测方法 |
CN115438402A (zh) * | 2022-08-22 | 2022-12-06 | 武汉大学 | 基于岩-机互馈感知关系的tbm掘进岩体条件实时预测方法及装置 |
CN115142160B (zh) * | 2022-08-22 | 2023-12-19 | 无锡物联网创新中心有限公司 | 一种纱线强力弱环的辨识方法及相关装置 |
CN115142160A (zh) * | 2022-08-22 | 2022-10-04 | 无锡物联网创新中心有限公司 | 一种纱线强力弱环的辨识方法及相关装置 |
CN115526036A (zh) * | 2022-09-19 | 2022-12-27 | 长安大学 | 一种岩爆规模等级的判别方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111079342B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111079342B (zh) | 一种基于在线岩体等级分类的tbm掘进性能预测方法 | |
CN110110419B (zh) | 一种基于多目标学习的tbm掘进参数预测方法 | |
CN107220734B (zh) | 基于决策树的数控车床车削过程能耗预测系统 | |
CN108803486B (zh) | 基于并联深度学习网络的数控机床热误差预测与补偿方法 | |
CN105321345B (zh) | 一种基于ARIMA模型和kalman滤波的道路交通流预测方法 | |
CN111127364B (zh) | 图像数据增强策略选择方法及人脸识别图像数据增强方法 | |
CN114329826B (zh) | 一种基于生成式对抗网络的平面叶栅定常流动预测方法 | |
CN107909179A (zh) | 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法 | |
CN110561191B (zh) | 基于pca与自编码器的数控机床刀具磨损数据处理方法 | |
CN112859898B (zh) | 一种基于双通道双向神经网络的飞行器轨迹预测方法 | |
CN106067034B (zh) | 一种基于高维矩阵特征根的配电网负荷曲线聚类方法 | |
CN113259331A (zh) | 一种基于增量学习的未知异常流量在线检测方法及系统 | |
CN111917785A (zh) | 一种基于de-gwo-svr的工业互联网安全态势预测方法 | |
CN113420506B (zh) | 掘进速度的预测模型建立方法、预测方法及装置 | |
CN112596031A (zh) | 一种基于知识图谱的目标雷达威胁度评估方法 | |
CN113033899A (zh) | 无人驾驶邻近车辆轨迹预测方法 | |
CN113806889A (zh) | 一种tbm刀盘扭矩实时预测模型的处理方法、装置以及设备 | |
CN112149896A (zh) | 一种基于注意力机制的机械设备多工况故障预测方法 | |
CN114297914A (zh) | 一种面向大电网可靠性评估的深度神经网络结果可信性保障方法 | |
CN112632705B (zh) | 一种基于gmm和高斯过程回归的路感模拟方法 | |
CN116906561B (zh) | 基于短时驾驶风格识别的车辆换挡点优化控制方法及系统 | |
CN117851891A (zh) | 基于数字孪生的cnn-lstm故障检测方法 | |
CN117272139A (zh) | 一种结合特征相关性与时间相关性的盾构机姿态预测方法 | |
CN113722982B (zh) | 一种基于多模型融合的汽车传感器攻击检测与防御方法 | |
CN114943741A (zh) | 一种动态场景下基于目标检测和几何概率的视觉slam方法 |
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 |