CN110110419B - 一种基于多目标学习的tbm掘进参数预测方法 - Google Patents

一种基于多目标学习的tbm掘进参数预测方法 Download PDF

Info

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
Application number
CN201910347964.6A
Other languages
English (en)
Other versions
CN110110419A (zh
Inventor
郑赢豪
李建斌
荆留杰
李鹏宇
鞠翔宇
武颖莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway Engineering Equipment Group Co Ltd CREG
Original Assignee
China Railway Engineering Equipment Group Co Ltd CREG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Railway Engineering Equipment Group Co Ltd CREG filed Critical China Railway Engineering Equipment Group Co Ltd CREG
Priority to CN201910347964.6A priority Critical patent/CN110110419B/zh
Publication of CN110110419A publication Critical patent/CN110110419A/zh
Application granted granted Critical
Publication of CN110110419B publication Critical patent/CN110110419B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/06Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
    • E21D9/08Making 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/06Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
    • E21D9/08Making 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/087Making 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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司机还可以综合围岩地质情况和渣片的颗粒分布来优化掘进过程中控制策略。
为解决上述问题,本发明所采用的技术方案通过一种基于多目标学习的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预测模型:
Figure BDA0002043001280000041
其中,D为目标变量的个数,d为目标变量的序号。
S4.2建立多输出高斯过程回归模型;
要想对D个SOGP预测模型
Figure BDA0002043001280000042
用MOGP模型来近似表示,就需要利用光滑核函数
Figure BDA0002043001280000043
和隐含函数u(x)的卷积运算来实现D个SOGP预测模型
Figure BDA0002043001280000044
的近似表示,具体步骤为:
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz;
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
Figure BDA0002043001280000045
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立。
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
Figure BDA0002043001280000051
式中,Sd,q表示由第d个SOGP预测模型fd(x)和隐函数uq(x)计算得到的方差系数,
Figure BDA0002043001280000052
表示第d个目标变量的计算精度矩阵。
S4.2.4,构造隐含函数的协方差矩阵为:
Figure BDA0002043001280000053
式中,
Figure BDA0002043001280000054
表示隐含函数的计算精度矩阵。
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
Figure BDA0002043001280000055
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模型估计多个目标变量的预测分布,具体如下:
Figure BDA0002043001280000061
式中,Kf,f表示特征样本与自身进行卷积计算所得的协方差矩阵;
Figure BDA0002043001280000062
表示新输入样本与特征样本进行卷积计算所得的协方差矩阵。
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预测模型:
Figure BDA0002043001280000101
其中,D为目标变量的个数,d为目标变量的序号。
S4.2建立多输出高斯过程回归模型;
要想对D个SOGP预测模型
Figure BDA0002043001280000102
用MOGP模型来近似表示,就需要利用光滑核函数
Figure BDA0002043001280000103
和隐含函数u(x)的卷积运算来实现D个SOGP预测模型
Figure BDA0002043001280000104
的近似表示,具体步骤为:
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz。
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
Figure BDA0002043001280000105
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立。
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
Figure BDA0002043001280000111
式中,Sd,q表示由第d个SOGP预测模型fd(x)和隐函数uq(x)计算得到的方差系数,
Figure BDA0002043001280000112
表示第d个目标变量的计算精度矩阵。
S4.2.4,构造隐含函数的协方差矩阵为:
Figure BDA0002043001280000113
式中,
Figure BDA0002043001280000114
表示隐含函数的计算精度矩阵。
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
Figure BDA0002043001280000115
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模型估计多个目标变量的预测分布,具体如下:
Figure BDA0002043001280000116
式中,Kf,f表示特征样本与自身进行卷积计算所得的协方差矩阵;
Figure BDA0002043001280000121
表示新输入样本与特征样本进行卷积计算所得的协方差矩阵。
S5,通过非线性优化方法对岩机作用模型的参数进行自适应调整,改善TBM掘进参数的预测精度。
所述非线性优化方法包括但不限于最大似然估计(Maximum likelihoodestimation),还包括粒子群优化(Particle Swarm Optimization),网格搜索(Grid Search),遗传算法(Genetic Algorithm)等方法。
本发明实施例以最大似然估计为例对MOGP模型的超参数集合Θ进行寻优,其最优估计值标记为Θ'。此时,上述MOGP模型的预测分布函数更新如下:
Figure BDA0002043001280000122
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预测模型:
Figure FDA0002043001270000031
其中,D为目标变量的个数,d为目标变量的序号;
S4.2建立多输出高斯过程回归模型;
要想对D个SOGP预测模型
Figure FDA0002043001270000032
用MOGP模型来近似表示,就需要利用光滑核函数
Figure FDA0002043001270000033
和隐含函数u(x)的卷积运算来实现D个SOGP预测模型
Figure FDA0002043001270000034
的近似表示,具体步骤为:
S4.2.1,对单目标变量的SOGP预测模型fd(x)用光滑核函数和隐含函数的卷积运算来近似表示,具体如下:
fd(x)=∫xGd(x-z)u(z)dz;
S4.2.2,为了提高SOGP预测模型的近似精度,将光滑核函数和隐含函数进一步线性分为若干子函数:
Figure FDA0002043001270000035
式中,Q为子函数个数,q为子函数的序号,Gd,q(x)为光滑核函数的子函数、uq(x)为隐含函数的子函数,且Gd,q(x)和uq(x)拥有相同的协方差函数,且相互独立;
S4.2.3,将光滑核函数和隐含函数选择相同的高斯形式,则光滑核函数的表达式如下:
Figure FDA0002043001270000041
式中,Sd,q表示由第d个SOGP预测模型fd(x)和隐函数uq(x)计算得到的方差系数,
Figure FDA0002043001270000042
表示第d个目标变量的计算精度矩阵;
S4.2.4,构造隐含函数的协方差矩阵为:
Figure FDA0002043001270000043
式中,
Figure FDA0002043001270000044
表示隐含函数的计算精度矩阵;
S4.2.5,根据步骤S4.2.2-S4.2.4计算任意两个目标变量之间的协方差,公式为:
Figure FDA0002043001270000045
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模型估计多个目标变量的预测分布,具体如下:
Figure FDA0002043001270000046
式中,Kf,f表示特征样本与自身进行卷积计算所得的协方差矩阵;
Figure FDA0002043001270000051
表示新输入样本与特征样本进行卷积计算所得的协方差矩阵。
7.根据权利要求1所述的基于多目标学习的TBM掘进参数预测方法,其特征在于:在步骤S5中,所述非线性优化方法能自适应调整步骤S4中所搭建的岩机作用模型参数,以提高TBM岩机模型的预测精度;
且所述非线性优化方法包括但不限于最大似然估计、粒子群优化、网格搜索、遗传算法。
CN201910347964.6A 2019-04-28 2019-04-28 一种基于多目标学习的tbm掘进参数预测方法 Active CN110110419B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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减振控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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