CN110110419B - 一种基于多目标学习的tbm掘进参数预测方法 - Google Patents
一种基于多目标学习的tbm掘进参数预测方法 Download PDFInfo
- Publication number
- CN110110419B CN110110419B CN201910347964.6A CN201910347964A CN110110419B CN 110110419 B CN110110419 B CN 110110419B CN 201910347964 A CN201910347964 A CN 201910347964A CN 110110419 B CN110110419 B CN 110110419B
- Authority
- CN
- China
- Prior art keywords
- tbm
- parameters
- target
- tunneling
- model
- 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
- 230000005641 tunneling Effects 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 72
- 230000001174 ascending effect Effects 0.000 claims abstract description 29
- 239000011435 rock Substances 0.000 claims abstract description 29
- 230000009471 action Effects 0.000 claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 13
- 238000011217 control strategy Methods 0.000 claims abstract description 11
- 239000002245 particle Substances 0.000 claims abstract description 11
- 230000036541 health Effects 0.000 claims abstract description 8
- 238000013473 artificial intelligence Methods 0.000 claims abstract description 7
- 238000011156 evaluation Methods 0.000 claims abstract description 7
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 7
- 239000002893 slag Substances 0.000 claims abstract description 6
- 230000009467 reduction Effects 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 33
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 230000002068 genetic effect Effects 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 abstract description 4
- 230000002452 interceptive effect Effects 0.000 abstract 1
- 208000022971 Tuberculous meningitis Diseases 0.000 description 89
- 208000001223 meningeal tuberculosis Diseases 0.000 description 89
- 238000010276 construction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
- E21D9/08—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining with additional boring or cutting means other than the conventional cutting edge of the shield
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
- E21D9/08—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining with additional boring or cutting means other than the conventional cutting edge of the shield
- E21D9/087—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining with additional boring or cutting means other than the conventional cutting edge of the shield with a rotary drilling-head cutting simultaneously the whole cross-section, i.e. full-face machines
-
- 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
-
- 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
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明提出了一种基于多目标学习的TBM掘进参数预测方法,以解决传统控制策略中人为决策的主观性与决策规则单一性的问题。步骤为:通过智能数据采集系统获取TBM掘进参数;对上述TBM掘进参数进行降噪与增强处理,获取每个循环上升段和稳定段数据;提取循环上升段与稳定段的时域特征,构建反映TBM健康状态的评估指标;利用人工智能算法有效表征上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系,建立基于多目标学习的岩机作用模型;通过非线性优化方法对岩机作用模型的参数进行自适应调整;根据岩机作用模型预估的多目标参数信息,TBM司机可综合围岩地质情况和渣片的颗粒分布优化TBM控制策略。本发明可在线估计TBM掘进参数,同时不干扰TBM正常运行。
Description
技术领域
本发明属于隧道大型掘进设备智能控制领域,具体涉及一种基于多目标学习的TBM掘进参数预测方法,可同时预估多个TBM掘进参数,优化TBM掘进过程中控制策略。
背景技术
在硬岩TBM的掘进过程中,司机需根据岩体地质条件对设备掘进参数进行实时调整。然而,传统的控制策略主要根据TBM司机的驾驶经验和掘进参数的波动情况进行人为判断,但这种方式易受人为主观性和决策规则单一性的影响而导致参数调整与实际值产生较大偏差,诱发TBM控制策略的决策失误,严重威胁人机的安全性与可靠性。
此外,在一个完整的TBM掘进循环内,司机通常需根据不同阶段的操作规则逐步改变TBM掘进及控制参数值,其中如何根据上升段数据来预测TBM稳定段掘进参数的变化趋势是现场TBM司机迫切关注的问题,也是目前TBM智能控制领域中的一个研究热点。因而,如何克服传统决策规则的单一性问题、实时精准预测TBM掘进参数,并以此辅助TBM司机优化控制策略对于实现TBM安全高效掘进具有至关重要的现实意义。
发明内容
针对上述现有技术描述的不足,本发明提供一种基于多目标学习的TBM掘进参数预测方法,通过基于多目标学习的人工智能算法有效表征上升段数据与稳定段多目标变量之间的非线性映射关系,从而实时精准预测多个TBM性能和控制参数。除了借助岩机模型所预估的多目标参数信息之外,TBM司机还可以综合围岩地质情况和渣片的颗粒分布来优化掘进过程中控制策略。
为解决上述问题,本发明所采用的技术方案通过一种基于多目标学习的TBM掘进参数预测方法来实现,其具体步骤如下:
S1,通过TBM智能数据采集系统获取TBM掘进参数的原始数据。
所述TBM掘进参数从现场TBM监控平台中智能数据采集系统获取,包括TBM性能参数、控制参数和设定参数。
所述TBM性能参数包括刀盘扭矩T和总推力F;TBM控制参数包括刀盘转速n和掘进速度V;TBM设定参数包括刀盘转速设定值np和掘进速度设定值Vp。
S2,对采集的TBM掘进参数进行降噪与增强处理,获取循环上升段和稳定段数据。
S2.1,对步骤S1所获TBM掘进参数进行但不限于剔除非工作状态下异常数据的初步处理。
S2.2,将初步处理的TBM掘进参数数据进行标准化,并划分循环上升段和稳定段。
S3,根据预处理后TBM掘进参数数据,提取反映TBM健康状态的评估指标。
所述评估指标,包括循环上升段TBM掘进数据的时域特征和稳定段掘进数据的时域特征。
S3.1,根据预处理后TBM掘进参数数据,提取循环上升段数据的时域特征。
所述循环上升段数据的时域特征,包括所有TBM性能参数、控制参数和设定参数在上升段的时域均值指标、均方根指标和方差指标;时域均值指标和均方根指标反映TBM设备在一段运行时间内的总体健康状态,方差指标有效表征TBM掘进参数在一段运行时间内的的波动情况;
S3.2,根据预处理后TBM掘进参数数据,提取稳定段数据的时域特征。
所述稳定段数据的时域特征为所有TBM性能参数和控制参数在稳定段的均值指标。
S4,利用人工智能算法有效表征上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系,建立基于多目标学习的岩机作用模型。
所述稳定段TBM多目标变量包括设备的刀盘扭矩、总推力、刀盘转速、推进速度。
且所述岩机作用模型包括但不限于BP(Back propagation)神经网络、多输出高斯过程回归(Multiple-Output Gaussian progress regression)、多输出支持向量回归(Multiple-Output SupportVector Regression)等方法。
本发明以多输出高斯过程回归(Multiple-Output Gaussian ProgressRegression-MOGP)为例阐述TBM岩机作用模型的构建过程,且循环上升段数据的时域特征为特征指标,稳定段数据为目标变量;其具体实施步骤如下:
S4.1构建单输出高斯过程回归模型;
单输出高斯过程回归模型(Single-Output Gaussian ProgressRegression-SOGP)首先对上升段TBM特征指标与每个稳定段TBM目标变量间的非线性关系分别进行建模,具体步骤为:
S4.1.1,给定样本特征矩阵X和单目标变量y;
X=[x1,x2,…,xi,…,xm];
y=[y1,y2,…,yi,…,ym];
其中,m为样本数量,xi为样本特征矩阵中第i个元素且由l个特征指标所组成的向量;yi为TBM目标变量中第i个元素。
S4.1.2,基于SOGP方法构造单目标变量与特征指标的预测模型;
yi=f(xi)+ε(i=1,2,…,m);
式中,ε表示均值为0且服从独立同分布的高斯函数,即ε~N(0,δ2),δ2表示高斯函数的方差;f(x)~N(0,K(x,x)),K表示特征样本间的协方差。
S4.1.3,重复步骤S4.1.1-S4.1.2,直至对每个目标变量都建立一个SOGP预测模型fd(x),得到D个SOGP预测模型:
其中,D为目标变量的个数,d为目标变量的序号。
S4.2建立多输出高斯过程回归模型;
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz;
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立。
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
S4.2.4,构造隐含函数的协方差矩阵为:
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
S4.2.6,根据步骤S4.2.5并通过卷积计算,构造多输出高斯过程回归模型,即MOGP预测模型,且MOGP预测模型就是岩机作用模型:
P(Y|X,Θ)=N(Y|0,Kf,f,Σ);
式中,Y表示D个目标变量所组成的矩阵;Kf,f表示卷积计算所得的协方差矩阵;Σ是零均值且服从独立同分布的高斯噪声矩阵,表示每个目标变量的输出值扰动;Θ表示MOGP模型的超参数集合。
S4.2.7,在根据步骤S1-S3获取新的上升段TBM特征向量x*之后,可利用步骤S4.2.6中MOGP模型估计多个目标变量的预测分布,具体如下:
S5,通过非线性优化方法对岩机作用模型的参数进行自适应调整,改善TBM掘进参数的预测精度。
且所述非线性优化方法包括但不限于最大似然估计(Maximum likelihoodestimation),粒子群优化(Particle Swarm Optimization),网格搜索(Grid Search),遗传算法(GeneticAlgorithm)等方法。
S6,根据岩机作用模型预估的多目标参数信息,TBM司机可综合围岩地质情况和渣片的颗粒分布优化TBM控制策略。
与现有技术相比较,本发明具有如下有益效果:
第一,本发明通过对TBM掘进数据提取均值、均方根、方差等特征指标,反映设备掘进过程中健康状态,其均值和均方根特征可描述TBM总体运行状态,方差指标可有效表征TBM掘进参数的波动情况;
第二,本发明利用人工智能算法对上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系进行建模,并通过在线优化方法自适应调整模型参数,从而实时预测TBM稳定段掘进参数的变化趋势;
第三,本发明根据岩机模型所预估的多目标参数信息,在综合围岩地质情况和渣片颗粒分布的基础之上,辅助TBM司机优化掘进过程中控制策略,可显著改善传统决策规则的单一性和人为判断的主观性缺陷。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的系统流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供了一种基于多目标学习的TBM掘进参数预测方法,如图1所示,其步骤如下:
S1,通过TBM智能数据采集系统获取TBM掘进参数的原始数据。
所述TBM掘进参数从现场TBM监控平台中智能数据采集系统获取,包括TBM性能参数、控制参数和设定参数。
所述TBM性能参数包括刀盘扭矩T和总推力F;TBM控制参数包括刀盘转速n和掘进速度V;TBM设定参数包括刀盘转速设定值np和掘进速度设定值Vp。
S2,对采集的TBM掘进参数进行降噪与增强处理,获取循环上升段和稳定段数据。
S2.1,对步骤S1所获得的TBM掘进参数进行但不限于剔除非工作状态下异常数据的初步处理。
S2.2,将初步处理的TBM掘进参数数据进行标准化,并划分成循环上升段和稳定段。
S3,根据预处理后TBM掘进参数数据,提取反映TBM健康状态的评估指标。
所述评估指标,包括循环上升段TBM掘进数据的时域特征和稳定段掘进数据的时域特征。
S3.1,根据预处理后TBM掘进参数数据,提取循环上升段数据的时域特征。
所述循环上升段数据的时域特征,包括所有TBM性能参数、控制参数和设定参数在上升段的时域均值指标、均方根指标和方差指标;时域均值指标和均方根指标反映TBM设备在一段运行时间内的总体健康状态,方差指标有效表征TBM掘进参数在一段运行时间内的的波动情况。
S3.2,根据预处理后TBM掘进参数数据,提取稳定段掘进数据的时域特征。
所述稳定段数据的时域特征为所有TBM性能参数和控制参数在稳定段的均值指标。
S4,利用人工智能算法有效表征上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系,建立基于多目标学习的岩机作用模型。
所述稳定段TBM多目标变量包括设备的刀盘扭矩、总推力、刀盘转速、推进速度。
所述岩机作用模型包括但不限于多输出高斯过程回归(Multiple-OutputGaussian ProgressRegression-MOGP),还包括BP(Back propagation)神经网络、多输出支持向量回归(Multiple-Output Support Vector Regression)等方法。
本发明实施例以多输出高斯过程回归(MOGP)为例阐述TBM岩机作用模型的构建过程,其具体实施步骤如下:
S4.1构建单输出高斯过程回归模型;
为了实现上述目的,单输出高斯过程回归模型(Single-Output GaussianProgress Regression-SOGP)首先对上升段TBM特征指标与每个稳定段TBM目标变量间的非线性关系分别进行建模。
本发明实施例以多输出高斯过程回归(Multiple-Output GaussianProgressRegression-MOGP)为例阐述TBM岩机作用模型的构建过程,且循环上升段数据的时域特征为特征指标,稳定段数据的时域特征为目标变量,其具体实施步骤如下:
S4.1构建单输出高斯过程回归模型;
单输出高斯过程回归模型(Single-Output Gaussian ProgressRegression-SOGP)首先对上升段TBM特征指标与每个稳定段TBM目标变量间的非线性关系分别进行建模,具体步骤为:
S4.1.1,给定样本特征矩阵X和单目标变量y;
X=[x1,x2,…,xi,…,xm];
y=[y1,y2,…,yi,…,ym];
其中,m为样本数量,xi为样本特征矩阵中第i个元素且由l个特征指标所组成的向量;yi为TBM目标变量中第i个元素。
S4.1.2,基于SOGP方法构造单目标变量与特征指标的预测模型;
yi=f(xi)+ε(i=1,2,…,m);
式中,ε表示均值为0且服从独立同分布的高斯函数,即ε~N(0,δ2),δ2表示高斯函数的方差;f(x)~N(0,K(x,x)),K表示特征样本间的协方差。
S4.1.3,重复步骤S4.1.1-S4.1.2,直至对每个目标变量都建立一个SOGP预测模型fd(x),得到D个SOGP预测模型:
其中,D为目标变量的个数,d为目标变量的序号。
S4.2建立多输出高斯过程回归模型;
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz。
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立。
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
S4.2.4,构造隐含函数的协方差矩阵为:
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
S4.2.6,根据步骤S4.2.5并通过卷积计算,构造多输出高斯过程回归模型,即MOGP预测模型,且MOGP预测模型就是岩机作用模型:
P(Y|X,Θ)=N(Y|0,Kf,f,Σ);
式中,Y表示D个目标变量所组成的矩阵;Kf,f表示卷积计算所得的协方差矩阵;Σ是零均值且服从独立同分布的高斯噪声矩阵,表示每个目标变量的输出值扰动;Θ表示MOGP模型的超参数集合。
S4.2.7,在根据步骤S1-S3获取新的上升段TBM特征向量x*之后,可利用步骤S4.2.6中MOGP模型估计多个目标变量的预测分布,具体如下:
S5,通过非线性优化方法对岩机作用模型的参数进行自适应调整,改善TBM掘进参数的预测精度。
所述非线性优化方法包括但不限于最大似然估计(Maximum likelihoodestimation),还包括粒子群优化(Particle Swarm Optimization),网格搜索(Grid Search),遗传算法(Genetic Algorithm)等方法。
本发明实施例以最大似然估计为例对MOGP模型的超参数集合Θ进行寻优,其最优估计值标记为Θ'。此时,上述MOGP模型的预测分布函数更新如下:
S6,根据岩机作用模型预估的多目标变量信息,TBM司机可综合围岩地质情况和渣片的颗粒分布优化TBM控制策略。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种基于多目标学习的TBM掘进参数预测方法,其特征在于,步骤如下:
S1,通过TBM智能数据采集系统获取TBM掘进参数的原始数据;
S2,对采集的TBM掘进参数进行降噪与增强处理,获取循环上升段和稳定段数据;
S3,根据预处理后TBM掘进参数数据,提取反映TBM健康状态的评估指标;
所述评估指标,包括循环上升段TBM掘进数据的时域特征和稳定段掘进数据的时域特征;
S4,利用人工智能算法有效表征上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系,建立基于多目标学习的岩机作用模型;
所述稳定段TBM多目标变量包括设备的刀盘扭矩、总推力、刀盘转速、推进速度;
S5,通过非线性优化方法对岩机作用模型的参数进行自适应调整,改善TBM掘进参数的预测精度;
S6,根据岩机作用模型预估的多目标参数信息,TBM司机综合围岩地质情况和渣片的颗粒分布优化TBM控制策略。
2.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:所述TBM掘进参数从现场TBM监控平台中智能数据采集系统获取,包括TBM性能参数、控制参数和设定参数;
所述TBM性能参数包括刀盘扭矩T和总推力F;TBM控制参数包括刀盘转速n和掘进速度V;TBM设定参数包括刀盘转速设定值np和掘进速度设定值Vp。
3.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S2中,具体步骤如下:
S2.1,对步骤S1获得的TBM掘进参数进行但不限于剔除非工作状态下异常数据的初步处理;
S2.2,将初步处理的TBM掘进数据进行标准化,并划分循环上升段和稳定段。
4.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S3中,具体步骤为:
S3.1,根据预处理后TBM掘进参数数据,提取循环上升段数据的时域特征;
所述循环上升段数据的时域特征,包括TBM性能参数、控制参数和设定参数在上升段的时域均值指标、均方根指标和方差指标;时域均值指标和均方根指标反映TBM设备在一段运行时间内的总体健康状态,方差指标有效表征TBM掘进参数在一段运行时间内的波动情况;
S3.2,根据预处理后TBM掘进参数数据,提取稳定段数据的时域特征;
所述稳定段数据的时域特征为TBM性能参数和控制参数在稳定段的均值指标。
5.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S4中,所述基于多目标学习的岩机作用模型利用人工智能算法有效表征上升段TBM特征指标与稳定段多目标变量之间的非线性映射关系;
且所述岩机作用模型采用包括但不限于BP神经网络、多输出高斯过程回归、多输出支持向量回归方法。
6.根据权利要求5所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S4中,以多输出高斯过程回归(Multiple-Output Gaussian progress regression-MOGP)方法来建立岩机作用模型,且循环上升段数据的时域特征为特征指标,稳定段数据的时域特征为目标变量,具体实施步骤如下:
S4.1构建单输出高斯过程回归模型;
单输出高斯过程回归模型(Single-Output Gaussian progress regression-SOGP)首先对上升段TBM特征指标与每个稳定段TBM目标变量间的非线性关系分别进行建模,具体步骤为:
S4.1.1,给定样本特征矩阵X和单目标变量y;
X=[x1,x2,…,xi,…,xm];
y=[y1,y2,…,yi,…,ym];
其中,m为样本数量,xi为样本特征矩阵中第i个元素且由l个特征指标所组成的向量;yi为TBM目标变量中第i个元素;
S4.1.2,基于SOGP方法构造单目标变量与特征指标的预测模型;
yi=f(xi)+ε(i=1,2,…,m);
式中,ε表示均值为0且服从独立同分布的高斯函数,即ε~N(0,δ2),δ2表示高斯函数的方差;f(x)~N(0,K(x,x)),K表示特征样本间的协方差;
S4.1.3,重复步骤S4.1.1-S4.1.2,直至对每个目标变量都建立一个SOGP预测模型fd(x),得到D个SOGP预测模型:
其中,D为目标变量的个数,d为目标变量的序号;
S4.2建立多输出高斯过程回归模型;
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz;
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立;
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
S4.2.4,构造隐含函数的协方差矩阵为:
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
S4.2.6,根据步骤S4.2.5并通过卷积计算,构造多输出高斯过程回归模型,即MOGP预测模型,且MOGP预测模型就是岩机作用模型:
P(Y|X,Θ)=N(Y|0,Kf,f,Σ);
式中,Y表示D个目标变量所组成的矩阵;Kf,f表示卷积计算所得的协方差矩阵;Σ是零均值且服从独立同分布的高斯噪声矩阵,表示每个目标变量的输出值扰动;Θ表示MOGP模型的超参数集合;
S4.2.7,在根据步骤S1-S3获取新的上升段TBM特征向量x*之后,利用步骤S4.2.6中MOGP模型估计多个目标变量的预测分布,具体如下:
7.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S5中,所述非线性优化方法能自适应调整步骤S4中所搭建的岩机作用模型参数,以提高TBM岩机模型的预测精度;
且所述非线性优化方法包括但不限于最大似然估计、粒子群优化、网格搜索、遗传算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910347964.6A CN110110419B (zh) | 2019-04-28 | 2019-04-28 | 一种基于多目标学习的tbm掘进参数预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910347964.6A CN110110419B (zh) | 2019-04-28 | 2019-04-28 | 一种基于多目标学习的tbm掘进参数预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110110419A CN110110419A (zh) | 2019-08-09 |
CN110110419B true CN110110419B (zh) | 2022-11-18 |
Family
ID=67487143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910347964.6A Active CN110110419B (zh) | 2019-04-28 | 2019-04-28 | 一种基于多目标学习的tbm掘进参数预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110110419B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110852423B (zh) * | 2019-11-12 | 2023-01-24 | 中铁工程装备集团有限公司 | 基于迁移学习的隧道掘进机掘进性能与控制参数预测方法 |
CN110990938B (zh) * | 2019-12-28 | 2023-04-18 | 中铁工程装备集团有限公司 | 一种掌子面在掘岩体状态软测量方法 |
CN110895730B (zh) * | 2019-12-31 | 2022-11-29 | 盾构及掘进技术国家重点实验室 | 一种基于lstm算法的tbm掘进参数的预测方法 |
CN111594198B (zh) * | 2020-04-30 | 2021-03-30 | 东北大学 | 一种开敞式tbm施工数据掘进循环的划分方法 |
CN111594201B (zh) * | 2020-05-28 | 2021-08-13 | 中铁工程装备集团有限公司 | 一种tbm关键参数智能控制系统及方法 |
CN111709650B (zh) * | 2020-06-18 | 2023-05-30 | 中铁十一局集团第四工程有限公司 | 一种滨海复杂地层盾构掘进适应性评价方法 |
CN112183993B (zh) * | 2020-09-22 | 2022-06-17 | 中国科学院武汉岩土力学研究所 | 一种tbm掘进控制参数智能预测及优化决策方法 |
CN112182973B (zh) * | 2020-09-30 | 2022-11-04 | 山东大学 | 考虑刀盘振动与岩渣几何信息的tbm操作参数优化方法 |
CN113033097B (zh) * | 2021-03-25 | 2023-09-19 | 东北大学 | 一种开敞式tbm掘进循环的划分、围岩等级、岩性预测方法 |
CN113946899A (zh) * | 2021-11-10 | 2022-01-18 | 广东粤海珠三角供水有限公司 | 一种基于实时地质信息的泥水盾构掘进参数预测方法 |
CN114692273B (zh) * | 2022-03-29 | 2024-09-17 | 中铁工程装备集团有限公司 | 面向tbm施工隧道地质字典建立方法及系统 |
CN114861039B (zh) * | 2022-03-31 | 2023-03-24 | 北京优特捷信息技术有限公司 | 一种搜索引擎的参数配置方法、装置、设备及存储介质 |
CN115034446B (zh) * | 2022-05-25 | 2024-10-15 | 同济大学 | 一种基于tbm运行数据的不良地质条件识别与预测方法 |
CN115618222B (zh) * | 2022-06-21 | 2023-05-05 | 北京交通大学 | 一种隧道掘进响应参数的预测方法 |
CN115982551A (zh) * | 2022-12-13 | 2023-04-18 | 中铁工程装备集团有限公司 | 掘进设备的故障诊断方法及智能化掘进设备 |
CN116842687A (zh) * | 2023-05-10 | 2023-10-03 | 北京交通大学 | Tbm掘进参数的辅助决策方法、装置及电子设备 |
CN118228117B (zh) * | 2024-05-23 | 2024-07-30 | 石家庄铁道大学 | 基于tbm掘进特征和振动特征参数的围岩等级感知方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CL2010000485A1 (es) * | 2009-05-13 | 2010-08-06 | Univ Sydney | Sistemas y metodos para el analisis y la sintesis de los datos a partir de una pluralidad de fuentes de datos de muestra para el aprendizaje y la regresion del proceso gaussiano. |
CN107577862A (zh) * | 2017-08-30 | 2018-01-12 | 中铁工程装备集团有限公司 | 一种tbm在掘岩体状态实时感知系统和方法 |
CN107632523A (zh) * | 2017-09-30 | 2018-01-26 | 中铁工程装备集团有限公司 | 一种硬岩tbm掘进控制参数智能决策方法及系统 |
CN107918279A (zh) * | 2017-11-20 | 2018-04-17 | 上海交通大学 | 一种基于粒子群优化算法pso的tbm减振控制方法 |
-
2019
- 2019-04-28 CN CN201910347964.6A patent/CN110110419B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CL2010000485A1 (es) * | 2009-05-13 | 2010-08-06 | Univ Sydney | Sistemas y metodos para el analisis y la sintesis de los datos a partir de una pluralidad de fuentes de datos de muestra para el aprendizaje y la regresion del proceso gaussiano. |
CN107577862A (zh) * | 2017-08-30 | 2018-01-12 | 中铁工程装备集团有限公司 | 一种tbm在掘岩体状态实时感知系统和方法 |
CN107632523A (zh) * | 2017-09-30 | 2018-01-26 | 中铁工程装备集团有限公司 | 一种硬岩tbm掘进控制参数智能决策方法及系统 |
CN107918279A (zh) * | 2017-11-20 | 2018-04-17 | 上海交通大学 | 一种基于粒子群优化算法pso的tbm减振控制方法 |
Non-Patent Citations (1)
Title |
---|
TBM掘进参数智能控制系统的研究与应用;张娜等;《隧道建设(中英文)》;20181020(第10期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110110419A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110419B (zh) | 一种基于多目标学习的tbm掘进参数预测方法 | |
CN111079342B (zh) | 一种基于在线岩体等级分类的tbm掘进性能预测方法 | |
CN112183993B (zh) | 一种tbm掘进控制参数智能预测及优化决策方法 | |
CN110852423B (zh) | 基于迁移学习的隧道掘进机掘进性能与控制参数预测方法 | |
CN110807557B (zh) | 一种基于bp神经网络的钻速预测方法和基于bp神经网络以及粒子群算法的钻速优化方法 | |
CN108764540B (zh) | 基于并行lstm串联dnn的供水管网压力预测方法 | |
CN107272403A (zh) | 一种基于改进粒子群算法的pid控制器参数整定算法 | |
CN111182564B (zh) | 一种基于lstm神经网络的无线链路质量预测方法 | |
CN111553118B (zh) | 基于强化学习的多维连续型优化变量全局优化方法 | |
CN110751318A (zh) | 一种基于ipso-lstm的超短期电力负荷预测方法 | |
CN111079343B (zh) | 一种基于宽度学习的风电机组有效风速估计方法 | |
CN112364560B (zh) | 矿山凿岩装备作业工时智能预测方法 | |
CN110245390B (zh) | 一种基于rs-bp神经网络的汽车发动机油耗预测方法 | |
CN113485443B (zh) | 基于深度学习的无人机控制方法、存储介质及设备 | |
CN114384931B (zh) | 一种基于策略梯度的无人机多目标最优控制方法和设备 | |
CN108594793A (zh) | 一种改进的rbf飞控系统故障诊断网络训练方法 | |
CN106708044A (zh) | 基于灰色预测混合遗传算法‑pid全垫升气垫船航向控制方法 | |
CN112149896A (zh) | 一种基于注意力机制的机械设备多工况故障预测方法 | |
CN114239396A (zh) | 一种风机齿轮箱状态预测方法及系统 | |
CN102141778B (zh) | 一种受核糖体rna启发的高阶控制器参数优化方法 | |
CN113807005B (zh) | 基于改进fpa-dbn的轴承剩余寿命预测方法 | |
CN115618497A (zh) | 一种基于深度强化学习的翼型优化设计方法 | |
CN115860176A (zh) | 一种基于组合型神经网络的电推船舶负荷预测方法及系统 | |
CN108538179A (zh) | 一种智能机械制造工人生产培训模拟系统 | |
CN118171742B (zh) | 基于残差估计的知识-数据驱动空战目标意图推理方法及系统 |
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 |