CN112100841A - 盾构机姿态预测方法、装置、终端设备以及存储介质 - Google Patents

盾构机姿态预测方法、装置、终端设备以及存储介质 Download PDF

Info

Publication number
CN112100841A
CN112100841A CN202010944273.7A CN202010944273A CN112100841A CN 112100841 A CN112100841 A CN 112100841A CN 202010944273 A CN202010944273 A CN 202010944273A CN 112100841 A CN112100841 A CN 112100841A
Authority
CN
China
Prior art keywords
angle data
predicted
data
yaw
pitch angle
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
Application number
CN202010944273.7A
Other languages
English (en)
Other versions
CN112100841B (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 20th Bureau Group Corp
Original Assignee
China Railway 20th Bureau Group Corp
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 20th Bureau Group Corp filed Critical China Railway 20th Bureau Group Corp
Priority to CN202010944273.7A priority Critical patent/CN112100841B/zh
Publication of CN112100841A publication Critical patent/CN112100841A/zh
Application granted granted Critical
Publication of CN112100841B publication Critical patent/CN112100841B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明公开一种盾构机姿态预测方法，所述方法包括以下步骤：获取盾构机的盾构姿态数据，所述盾构姿态数据包括所述盾构机的俯仰角数据、横摆角数据和滚动角数据；将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据；根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果。本发明还公开了一种盾构机姿态预测装置、终端设备以及存储介质。根据预测俯仰角数据、预测横摆角数据和预测滚动角数据获得所述盾构机的姿态预测结果准确率高。

Description

盾构机姿态预测方法、装置、终端设备以及存储介质
技术领域
本发明涉及盾构机施工工程技术领域，特别涉及一种盾构机姿态预测方 法、装置、终端设备以及存储介质。
背景技术
在盾构施工过程中，盾构姿态控制不利容易引起地面沉降、管片错台和 土体超挖等工程问题。
相关技术中公布了一种基于混合深度学习的盾构掘进姿态智能预测方法 及系统(申请号：CN110195592A)，将影响盾构姿态的主要因素作为输入变 量，盾构姿态的响应参数作为输出变量，在采用小波变换进行数据去噪的基 础上，利用CNNs进行特征提取、LSTM算法构建预测模型，以实现对盾构 姿态的预测。
但是，采用上述方法对盾构姿态进行预测，盾构机的预测结果准确率较 低。
发明内容
本发明的主要目的是提供一种盾构机姿态预测方法、装置、终端设备以 及存储介质，旨在解决现有技术中盾构机的预测结果准确率较低的技术问题。
为实现上述目的，本发明提出的一种盾构机姿态预测方法，所述方法包 括以下步骤：
获取盾构机的盾构姿态数据，所述盾构姿态数据包括所述盾构机的俯仰 角数据、横摆角数据和滚动角数据；
将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角 预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预 测横摆角数据和预测滚动角数据；
根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角数据， 获得所述盾构机的姿态预测结果。
可选的，所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分 别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯 仰角数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
获取盾构机的历史俯仰角数据、历史横摆角数据和历史滚动角数据；
利用公式一分别对所述历史俯仰角数据、所述历史横摆角数据和所述历 史滚动角数据进行处理，获得训练俯仰角数据、训练横摆角数据和训练滚动 角数据；
利用训练俯仰角数据、训练横摆角数据和训练滚动角数据分别对第一门 控循环单元神经网络、第二门控循环单元神经网络和第三门控循环单元神经 网络进行训练，获得所述俯仰角预测模型、所述横摆角预测模型和所述滚动 角预测模型。
所述公式一为：
Figure BDA0002674018810000021
其中，当x为历史俯仰角数据时，x′为训练俯仰角数据，xmin为历史俯仰 角数据中的最小值、xmax为历史俯仰角数据中的最大值；当x为历史横摆角数 据时，x′为训练横摆角数据，xmin为历史横摆角数据中的最小值、xmax为历 史横摆角数据中的最大值；当x为历史滚动角数据时，x′为训练滚动角数据， xmin为历史滚动角数据中的最小值、xmax为历史滚动角数据中的最大值。
可选的，所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分 别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯 仰角数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
利用所述俯仰角数据绘制俯仰角箱形图，并利用所述俯仰角箱形图在所 述所述俯仰角数据中筛选有效俯仰角数据；
利用所述横摆角数据绘制横摆角箱形图，并利用所述横摆角箱形图在所 述所述横摆角数据中筛选有效横摆角数据；
利用所述滚动角数据绘制滚动角箱形图，并利用所述滚动角箱形图在所 述所述滚动角数据中筛选有效滚动角数据；
利用线性内插法分别对所述有效俯仰角数据、所述有效横摆角数据和所 述有效滚动角数据进行处理，获得输入俯仰角数据、输入横摆角数据和输入 滚动角数据；
所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯 仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、 预测横摆角数据和预测滚动角数据，获得预测俯仰角数据、预测横摆角数据 和预测滚动角数据的步骤包括：
将所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据分 别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯 仰角数据、预测横摆角数据和预测滚动角数据，获得预测俯仰角数据、预测 横摆角数据和预测滚动角数据。
可选的，所述将所述输入俯仰角数据、所述输入横摆角数据和所述输入 滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中， 获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，获得预测俯仰角 数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
利用经验模态分解法分别对所述输入俯仰角数据、所述输入横摆角数据 和所述输入滚动角数据进行分解，获得具有时间序列的多组俯仰角数据、具 有时间序列的多组横摆角数据和具有时间序列的多组滚动角数据；
所述将所述有效俯仰角数据、所述有效横摆角数据和所述有效滚动角数 据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预 测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤包括：
将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角数据分 别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得多组预 测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据；
所述根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角 数据，获得所述盾构机的姿态预测结果的步骤包括：
根据所述多组预测俯仰角数据、所述多组预测横摆角数据和所述多组预 测滚动角数据，获得所述盾构机的姿态预测结果。
可选的，所述将所述多组俯仰角数据、所述多组横摆角数据和所述多组 滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中， 获得多组预测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据的步 骤包括：
将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角数据分 别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，利用滑动窗 口的方式进行滚动预测，获得多组预测俯仰角数据、多组预测横摆角数据和 多组预测滚动角数据。
可选的，所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据 和所述多组预测滚动角数据，获得所述盾构机的姿态预测结果的步骤包括：
获取所述多组预测俯仰角数据的俯仰角权重、所述多组预测横摆角数据 的横摆角权重和所述多组预测滚动角数据的滚动角权重；
根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所述 多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及所述滚动角权重， 获得结果俯仰角数据、结果横摆角数据和结果滚动角数据；
根据所述结果俯仰角数据、所述结果横摆角数据和所述结果滚动角数据， 获得所述盾构机的姿态预测结果。
可选的，所述根所述根据所述多组预测俯仰角数据、所述多组预测横摆 角数据、所述多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及所 述滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据的 步骤包括：
利用调和平均误差最小化法分别对所述所述俯仰角权重、所述横摆角权 和所述滚动角权重进行处理，获得选定俯仰角权重、选定横摆角权重和选定 滚动角权重；
根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所述 多组预测滚动角数据、所述选定俯仰角权重、所述选定横摆角权以及所述选 定滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据。
此外，为实现上述目的，本发明还提出了一种盾构机姿态预测装置，所 述装置包括：
获取模块，用于获取盾构机的盾构姿态数据，所述盾构姿态数据包括所 述盾构机的俯仰角数据、横摆角数据和滚动角数据；
预测模块，用于将所述俯仰角数据、所述横摆角数据和所述滚动角数据 分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测 俯仰角数据、预测横摆角数据和预测滚动角数据；
结果获得模块，用于根据所述预测俯仰角数据、所述预测横摆角数据和 所述预测滚动角数据，获得所述盾构机的姿态预测结果。
此外，为实现上述目的，本发明还提出了一种终端设备，所述终端设备 包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的盾 构机姿态预测程序，所述盾构机姿态预测程序配置为实现如上述任一项所述 的盾构机姿态预测方法的步骤。
此外，为实现上述目的，本发明还提出了一种存储介质，所述存储介质 上存储有盾构机姿态预测程序，所述盾构机姿态预测程序被处理器执行时实 现如上述任一项所述的盾构机姿态预测方法的步骤。
本发明技术方案采用一种盾构机姿态预测方法，通过获取盾构机的盾构 姿态数据，所述盾构姿态数据包括所述盾构机的俯仰角数据、横摆角数据和 滚动角数据；将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输 入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角 数据、预测横摆角数据和预测滚动角数据；根据所述预测俯仰角数据、所述 预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果。 由于，将盾构姿态数据数据拆分为俯仰角数据、横摆角数据和滚动角数据，并分别利用俯仰角预测模型、横摆角预测模型和滚动角预测模型进行预测， 获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，同时，预测俯仰 角数据、预测横摆角数据和预测滚动角数据相对于现有技术的整体预测结果 准确率高，平稳性好，所以，根据预测俯仰角数据、预测横摆角数据和预测 滚动角数据获得所述盾构机的姿态预测结果准确率高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实 施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面 描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的 附图。
图1是本发明实施例方案涉及的硬件运行环境的盾构机姿态预测设备结 构示意图；
图2为本发明盾构机姿态预测方法第一实施例的流程示意图；
图3为本发明利用箱型图法筛选有效数据原理图；
图4为本发明盾构机姿态预测方法第二实施例步骤S12之前的流程示意图；
图5为本发明盾构机姿态预测装置第一实施例的结构框图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步 说明。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行 清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例， 而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有 作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
随着城市地铁的大规模建设，盾构法施工因其特有的施工工程优势得以 推广。在盾构施工工程过程中，盾构姿态控制不利容易引起地面沉降、管片 错台和土体超挖等工程问题，因此，盾构姿态控制是施工过程中需要重点关 注的关键业务场景之一。针对盾构姿态的控制问题，目前主要聚焦在基于姿 态纠偏的事后反馈控制和基于预测的事前主动控制两个方面。前者由于存在 时间上的滞后性，导致其控制措施的执行不及时，无法针对盾构姿态的实时 演变给出及时有效的反馈；后者能够在姿态发生变化前对姿态的未来趋势给 出及时预测，可辅助现场施工人员和管理人员制定下一步的掘进方案，以保 障施工过程中姿态得以及时控制。基于此，盾构姿态预测在盾构姿态控制方 面具有十分重要的意义。
从目前公开的资料来看，关于盾构姿态预测方法的研究已经公开了若干 发明，目前最为有效的预测方法为基于人工智能的预测方法。
周诚等人公开了申请公布号为CN 110195592 A，名称为“基于混合深度 学习的盾构掘进姿态智能预测方法及系统”的盾构姿态预测方法，主要是将 影响盾构姿态的主要因素作为输入变量，盾构姿态的响应参数作为输出变量， 在采用小波变换进行数据去噪的基础上，采用CNNs进行特征提取，进而依托 LSTM算法构建预测模型进行盾构姿态预测。
尹清锋等人公开了申请公布号为CN 111365015 A，名称为“一种基于 XGBoost的盾构掘进参数特征提取与姿态偏差预测方法”的盾构姿态预测方法， 主要是在基于专家经验和数据分析进行影响盾构姿态的相关因素特征提取的 基础上，采用XGBoost算法构建盾构姿态预测模型以实现相应的偏差预测。
上述关于盾构姿态的智能化预测，有效的实现了盾构姿态的实时、较为 准确的预测，但仍存在以下几点问题：
(1)既有公开专利均是从掘进参数与盾构姿态之间的耦合关系出发构建 预测模型进行分析。这种方式可以实现较为准确的预测，但其受盾构姿态相 关影响参数选择的限制，姿态参数众多，如何选择适宜的参数进行研究分析 不确定性较高，且输入变量维度较高容易导致模型计算开销大，影响其工程 应用。
(2)就盾构姿态预测问题而言，其产生的数据是多种情形的耦合，可能 存在随机项、趋势项等情形，不对这些加以分解的整体式预测，其预测的精 度可能在一定程度上会受到影响。
上述关于盾构姿态的智能化预测，有效的实现了盾构姿态的实时、较为 准确的预测，但仍存在以下几点问题：
(1)既有公开专利均是从掘进参数与盾构姿态之间的耦合关系出发构建 预测模型进行分析。这种方式可以实现较为准确的预测，但其受盾构姿态相 关影响参数选择的限制，姿态参数众多，如何选择适宜的参数进行研究分析 不确定性较高，且输入变量维度较高容易导致模型计算开销大。
(2)就盾构姿态预测问题而言，其产生的数据是多种情形的耦合，可能 存在随机项、趋势项等情形，不对这些加以分解的整体式预测，其预测的精 度可能在一定程度上会受到影响。
针对上述存在的问题，亟待解决的关键技术主要包括：
(1)如果实现从盾构姿态本身出发的单点预测，不考虑耦合回归的多点 预测，而是考虑未来时间盾构姿态变化与过去时刻姿态之间的关系，进而构 建预测模型；
(2)针对获取到的盾构姿态数据，如何将其拆分为多组平稳数据进行预 测研究，以提高预测精度；
(3)采用何种算法进行预测，其预测精度和预测效率更高。
参照图1，图1为本发明实施例方案涉及的硬件运行环境的终端设备结构 示意图。
终端设备可以是移动电话、智能电话、笔记本电脑、数字广播接收器、 个人数字助理(PDA)、平板电脑(PAD)等用户设备(User Equipment，UE)、手持 设备、车载设备、可穿戴设备、计算设备或连接到无线调制解调器的其它处 理设备、移动台(Mobile station，MS)等。设备可能被称为用户终端、便携式 终端、台式终端等。
通常，终端设备包括：至少一个处理器301、存储器302以及存储在所述 存储器上并可在所述处理器上运行的盾构机姿态预测程序，所述盾构机姿态 预测程序配置为实现如前所述的盾构机姿态预测方法的步骤。
处理器301可以包括一个或多个处理核心，比如4核心处理器、8核心处理 器等。处理器301可以采用DSP(Digital Signal Processing，数字信号处理)、 FPGA(Field－Programmable Gate Array，现场可编程门阵列)、 PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实 现。处理器301也可以包括主处理器和协处理器，主处理器是用于对在唤醒状 态下的数据进行处理的处理器，也称CPU(CentralProcessingUnit，中央处理器)； 协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实 施例中，处理器301可以在集成有GPU(Graphics Processing Unit，图像处理器)， GPU用于负责显示屏所需要显示的内容的渲染和绘制。处理器301还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关盾构机 姿态预测方法操作，使得盾构机姿态预测方法模型可以自主训练学习，提高 效率和准确度。
存储器302可以包括一个或多个计算机可读存储介质，该计算机可读存储 介质可以是非暂态的。存储器302还可包括高速随机存取存储器，以及非易失 性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中， 存储器302中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少 一个指令用于被处理器801所执行以实现本申请中方法实施例提供的盾构机 姿态预测方法。
在一些实施例中，终端设备还可选包括有：通信接口303和至少一个外围 设备。处理器301、存储器302和通信接口303之间可以通过总线或信号线相连。 各个外围设备可以通过总线、信号线或电路板与通信接口303相连。具体地， 外围设备包括：射频电路304、显示屏305和电源306中的至少一种。
通信接口303可被用于将I/O(Input/Output，输入/输出)相关的至少一个外 围设备连接到处理器301和存储器302。。在一些实施例中，处理器301、存储 器302和通信接口303被集成在同一芯片或电路板上；在一些其他实施例中， 处理器301、存储器302和通信接口303中的任意一个或两个可以在单独的芯片 或电路板上实现，本实施例对此不加以限定。
射频电路304用于接收和发射RF(Radio Frequency，射频)信号，也称电磁 信号。射频电路304通过电磁信号与通信网络以及其他通信设备进行通信。射 频电路304将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转 换为电信号。可选地，射频电路304包括：天线系统、RF收发器、一个或多个 放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块 卡等等。射频电路304可以通过至少一种无线通信协议来与其它终端进行通信。 该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例 中，射频电路304还可以包括NFC(Near Field Communication，近距离无线通信) 有关的电路，本申请对此不加以限定。
显示屏305用于显示UI(User Interface，用户界面)。该UI可以包括图形、 文本、图标、视频及其它们的任意组合。当显示屏305是触摸显示屏时，显示 屏305还具有采集在显示屏305的表面或表面上方的触摸信号的能力。该触摸 信号可以作为控制信号输入至处理器301进行处理。此时，显示屏305还可以 用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中， 显示屏305可以为一个，电子设备的前面板；在另一些实施例中，显示屏305 可以为至少两个，分别设置在电子设备的不同表面或呈折叠设计；在再一些 实施例中，显示屏305可以是柔性显示屏，设置在电子设备的弯曲表面上或折 叠面上。甚至，显示屏305还可以设置成非矩形的不规则图形，也即异形屏。 显示屏305可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
电源306用于为电子设备中的各个组件进行供电。电源306可以是交流电、 直流电、一次性电池或可充电电池。当电源306包括可充电电池时，该可充电 电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。 本领域技术人员可以理解，图1中示出的结构并不构成对盾构机姿态预测设备 的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同 的部件布置。
此外，本发明实施例还提出一种存储介质，所述存储介质上存储有盾构 机姿态预测程序，所述盾构机姿态预测程序被处理器执行时实现如上文所述 的盾构机姿态预测方法的步骤。因此，这里将不再进行赘述。另外，对采用 相同方法的有益效果描述，也不再进行赘述。对于本申请所涉及的计算机可 读存储介质实施例中未披露的技术细节，请参照本申请方法实施例的描述。 确定为示例，程序指令可被部署为在一个计算设备上执行，或者在位于一个 地点的多个计算设备上执行，又或者，在分布在多个地点且通过通信网络互连的多个计算设备上执行。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程， 是可以通过计算机程序来指令相关的硬件来完成，上述的程序可存储于一计 算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的 流程。其中，上述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(RandomAccessMemory，RAM)等。
基于上述硬件结构，提出本发明盾构机姿态预测方法的实施例。
参照图2，图2为本发明盾构机姿态预测方法第一实施例的流程示意图； 所述方法包括以下步骤：
步骤S12：获取盾构机的盾构姿态数据，所述盾构姿态数据包括所述盾构 机的俯仰角数据、横摆角数据和滚动角数据。
需要说明的是，本发明实施例的执行主体是上文所述的终端设备，盾构 机可以是盾构机和TBM中的任意一种；本发明的盾构姿态数据可以是盾构机 正在运行时获取到的该时刻的盾构姿态数据(即，实时盾构数据)，也可以是 盾构机最近一次运行时候的盾构数据，为避免盾构机因为地形的复杂出现故 障，实时盾构数据取当前时刻的盾构数据为优，本发明不作具体限制。
可以理解的是，获取到的盾构姿态数据是多种用于描述盾构姿态数据的 集合，即，盾构机的俯仰角数据、横摆角数据和滚动角数据数据的集合，需 要将盾构姿态数据进行数据分解，获得盾构机的俯仰角数据、横摆角数据和 滚动角数据数据。
步骤S12：将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输 入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角 数据、预测横摆角数据和预测滚动角数据。
需要说明的是，不同的盾构姿态数据对应的预测模型不同，即，所述俯 仰角数据、所述横摆角数据和所述滚动角数据分别对应俯仰角预测模型、横 摆角预测模型和滚动角预测模型，需要将所述俯仰角数据、所述横摆角数据 和所述滚动角数据分别输入对应的预测模型中进行预测，以获得较准确，较 平稳的预测俯仰角数据、预测横摆角数据和预测滚动角数据。
步骤S13：根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚 动角数据，获得所述盾构机的姿态预测结果。
可以理解的是，将获得的所述预测俯仰角数据、所述预测横摆角数据和 所述预测滚动角数据进行合并，可以获得盾构机的预测盾构姿态数据，并根 据预测盾构姿态数据，获得姿态预测结果。
本实施例技术方案采用一种盾构机姿态预测方法，通过获取盾构机的盾 构姿态数据，所述盾构姿态数据包括所述盾构机的俯仰角数据、横摆角数据 和滚动角数据；将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别 输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰 角数据、预测横摆角数据和预测滚动角数据；根据所述预测俯仰角数据、所 述预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果。 由于，将盾构姿态数据数据拆分为俯仰角数据、横摆角数据和滚动角数据，并分别利用俯仰角预测模型、横摆角预测模型和滚动角预测模型进行预测， 获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，同时，预测俯仰 角数据、预测横摆角数据和预测滚动角数据相对于现有的整体预测结果准确 率高，平稳性好，所以，根据预测俯仰角数据、预测横摆角数据和预测滚动 角数据获得所述盾构机的姿态预测结果准确率高。
进一步的，步骤S12之前，所述方法还包括：利用所述俯仰角数据绘制 俯仰角箱形图，并利用所述俯仰角箱形图在所述所述俯仰角数据中筛选有效 俯仰角数据；利用所述横摆角数据绘制横摆角箱形图，并利用所述横摆角箱 形图在所述所述横摆角数据中筛选有效横摆角角数据；利用所述滚动角数据 绘制滚动角箱形图，并利用所述滚动角箱形图在所述所述滚动角数据中筛选 有效滚动角数据；利用线性内插法分别对所述有效俯仰角数据、所述有效横 摆角数据和所述有效滚动角数据进行处理，获得输入俯仰角数据、输入横摆角数据和输入滚动角数据。
相应的，步骤S12包括：将所述输入俯仰角数据、所述输入横摆角数据 和所述输入滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角 预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，获 得预测俯仰角数据、预测横摆角数据和预测滚动角数据。
需要说明的是，获取到盾构机的姿态数据是原始数据，所以，所述俯仰 角数据、所述横摆角数据和所述滚动角数据也为原始数据，需要在原始数据 中筛选有效数据(剔除原始数据的异常值)，即，利用箱型图法，分别在所述 俯仰角数据、所述横摆角数据和所述滚动角数据中筛选有效俯仰角数据、有 效横摆角数据以及有效滚动角数据；并利用线性内插法分别对所述有效俯仰 角数据、所述有效横摆角数据和所述有效滚动角数据进行缺失值的填充处理， 获得输入俯仰角数据、输入横摆角数据和输入滚动角数据。
参照图3，图3为本发明利用箱型图法筛选有效数据原理图；箱型图法的 各个数据计算方法为：
Figure BDA0002674018810000121
其中，U为上四分位数，L为下四分位数，IQR＝U-L，Upper为上界，Lower 为下界。
可以理解的是，所述俯仰角数据、所述横摆角数据和所述滚动角数据对 应的箱型图不同，需要根据数据信息绘制不同的箱型图，以对所述俯仰角数 据、所述横摆角数据和所述滚动角数据进行有效数据的筛选；上界和下界以 外的数据，即为原始数据中的异常值，需要把异常值在原始数据集中剔除， 获得剔除后的有效数据，即，有效俯仰角数据为所述俯仰角数据中位于俯仰 角箱型图的上下界之间的数据，有效横摆角数据为所述横摆角数据中位于横 摆角箱型图的上下界之间的数据，有效滚动角数据为所述滚动角角数据中位 于滚动角箱型图的上下界之间的数据；获得有效数据后，由于，各个有效数 据中缺失了异常值，需要利用线性内插法对所述有效俯仰角数据、所述有效 横摆角数据和所述有效滚动角数据进行缺失值的填充处理，获得数据信息完 整的输入俯仰角数据、输入横摆角数据和输入滚动角数据。
进一步的，步骤S12之前，所述方法还包括：利用经验模态分解法分别 对所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据进行分 解，获得具有时间序列的多组俯仰角数据、具有时间序列的多组横摆角数据 和具有时间序列的多组滚动角数据。
相应的，步骤S12包括：将所述多组俯仰角数据、所述多组横摆角数据 和所述多组滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角 预测模型中，获得多组预测俯仰角数据、多组预测横摆角数据和多组预测滚 动角数据。
同时，步骤S13包括：根据所述多组预测俯仰角数据、所述多组预测横 摆角数据和所述多组预测滚动角数据，获得所述盾构机的姿态预测结果。
需要说明的是，采用经验模态分解法对所述输入俯仰角数据、所述输入 横摆角数据和所述输入滚动角数据进行分解，获得具有时间序列的多组俯仰 角数据、具有时间序列的多组横摆角数据和具有时间序列的多组滚动角数据。 对于所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据中的 任意一组数据，均进行分解处理，本实施例以所述输入滚动角数据为例进行 讲解。
经验模态分解法对输入滚动角数据进行分解，获得若干本征模函数(IMF) 和一个剩余变量(RES)，即，获得多组滚动角数据，具体的分解步骤为：
1)找到输入滚动角数据y(t)的全部极值点，通过三次样条插值的方式获 取该序列数据的上下包络曲线BLmax(t)和BLmin(t)，并计算其均值AV(t)，即， AV(t)＝BLmax(t)-BLmin(t)；
2)计算输入滚动角数据y(t)与均值AV(t)之差I(t)，即，I(t)＝y(t)-AV(t)；
3)如果I(t)函数的局部极值数据和过零点数目相差不大于1，且所有时点 的上、下包络线取平均均等于零，则说明其为一个本征模函数，记作I1(t)，反 之，令y(t)＝I(t)，返回步骤1)；
4)将得到的本征模函数I(t)从输入滚动角数据y(t)中剔除，可得剩余序列 数据R(t)，即，R(t)＝y(t)-I1(t)；
5)倘若R(t)为单调函数，则将其作为剩余变量RES，否则，返回步骤1)。
基于上述步骤，参照输入滚动角数据分解法，将盾构姿态数据包括的输 入俯仰角数据、输入横摆角数据和输入滚动角数据分解为如下集合：
Figure BDA0002674018810000141
其中，y1(t),y2(t),y3(t)分别为分解之后的多组俯仰角数据的集合、多组横摆 角数据的集合和多组滚动角数据的集合。I11(t),I12(t),…,I1v-1(t),I21(t),I22(t),…,I2v-1(t), I31(t),I32(t),…,I3v-1(t)分别为各姿态参数基于经验模态分解的v-1个本征模函数，RES1,RES2,RES3分别为剩余变量。
然后，将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角 数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得 多组预测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据。例如， 一个I11(t)对应一个预测数据。
进一步的，步骤S12包括：将所述多组俯仰角数据、所述多组横摆角数 据和所述多组滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动 角预测模型中，利用滑动窗口的方式进行滚动预测，获得多组预测俯仰角数 据、多组预测横摆角数据和多组预测滚动角数据。
需要说明的是，采用滑动窗口的方式进行滚动预测；首先设定窗口宽度w， 即，用过去多少时刻的数据来预测当前时刻，对于将所述多组俯仰角数据、 所述多组横摆角数据和所述多组滚动角数据中的任意一种数据，均采用下述 方法进行预测，此处以多组俯仰角数据为例；利用公式2获得多组俯仰角数 据的表达形式，公式2为：
yi(t)＝f(yi(t-1),yi(t-2),…,yi(t-w))，i＝1,2,…v
然后，将yi(t-1),yi(t-2),…,yi(t-w)作为俯仰角预测模型的输入，利用公式3进行预测，获得预测俯仰角数据，公式3为：
Figure RE-GDA0002748827330000161
其中，式中，Xt为t时刻该层的输入变量，Rt和Zt分别为重置门和更新 门，Ht为隐含状态，Ht-1为候选隐含状态，Wxr、Wxz、Wxh∈Rd×h和 Whr、Whz、Whh∈Rh×h为权重参数，br、bz、bh∈R1×h为偏差参数，h为隐藏单元个数， d为输入变量维度，σ和tanh为激活函数，e为矩阵乘法。
可以理解的是，每组俯仰角数据均获得一个预测俯仰角数据，多组俯仰 角数据获得多组预测俯仰角数据，即，所述多组俯仰角数据、所述多组横摆 角数据和所述多组滚动角数据分别对应多组预测俯仰角数据、多组预测横摆 角数据和多组预测滚动角数据。
进一步的，步骤S13包括：获取所述多组预测俯仰角数据的俯仰角权重、 所述多组预测横摆角数据的横摆角权重和所述多组预测滚动角数据的滚动角 权重；根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所 述多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及所述滚动角权 重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据；根据所述结 果俯仰角数据、所述结果横摆角数据和所述结果滚动角数据，获得所述盾构 机的姿态预测结果。
需要说明的是，多组预测横摆角数据的权重是不同的，需要根据多组预 测横摆角各自的横摆角权重，确定结果横摆角数据，即，最终预测横摆角数 据，同理，获得最终预测俯仰角数据和最终预测滚动角数据，利用最终预测 横摆角数据、最终预测俯仰角数据和最终预测滚动角数据，获得盾构机的预 测盾构姿态数据，并根据盾构机的预测盾构姿态数据，获得盾构机的姿态预 测结果；同时，权重在不同的盾构机施工环境可能不同，用户可以根据自己 的需求进行设定，本发明不做限制。
具体应用中，以多组预测俯仰角数据为例，利用公式3获得结果俯仰角 数据，公式4为：
Figure BDA0002674018810000161
其中，y结果俯仰角数据，y1,y2,..,yv分别多组输入俯仰角数据分解后的v-1 个本征模函数和剩余变量的预测结果，
Figure BDA0002674018810000162
为各个预测俯仰角数据的 俯仰角权重。同理，可获得结果横摆角角数据和结果滚动角数据。
进一步的，所述根所述根据所述多组预测俯仰角数据、所述多组预测横 摆角数据、所述多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及 所述滚动角权重，获得结果滚动角数据、结果横摆角数据和结果滚动角数据 的步骤包括：利用调和平均误差最小化法分别对所述所述俯仰角权重、所述 横摆角权和所述滚动角权重进行处理，获得选定俯仰角权重、选定横摆角权 重和选定滚动角权重；根所述根据所述多组预测俯仰角数据、所述多组预测 横摆角数据、所述多组预测滚动角数据、所述选定俯仰角权重、所述选定横 摆角权以及所述选定滚动角权重，获得结果俯仰角数据、结果横摆角数据和 结果滚动角数据。
需要说明的是，获取所述多组预测俯仰角数据的俯仰角权重、所述多组 预测横摆角数据的横摆角权重和所述多组预测滚动角数据的滚动角权重之后， 利用调和平均误差最小化法对获得的权重进行处理，获得选定权重。
具体应用中，以多组预测俯仰角数据为例，利用公式5获得选定俯仰角 权重，公式5为：
min e2
Figure BDA0002674018810000163
s.t.
Figure BDA0002674018810000164
Figure BDA0002674018810000165
其中，e为平均误差，y1,y2,..,yv分别多组输入俯仰角数据分解后的v-1个本 征模函数和剩余变量的预测结果，
Figure BDA0002674018810000166
为各个预测俯仰角数据的选定俯仰 角权重，L表示省略号。同理，利用上述方法获得选定横摆角权重和选定盾冬 权重，以根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、 所述多组预测滚动角数据、所述选定俯仰角权重、所述选定横摆角权以及所 述选定滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数 据。
可以理解的是，将获得的所述结果俯仰角数据、所述结果横摆角数据和 所述结果滚动角数据进行合并，可以获得盾构机的预测盾构姿态数据，并根 据预测盾构姿态数据，获得姿态预测结果。
参照图4，图4为本发明盾构机姿态预测方法第二实施例步骤S12之前的 流程示意图；所述方法包括以下步骤：
步骤S21：获取盾构机的历史俯仰角数据、历史横摆角数据和历史滚动角 数据；
步骤S22：利用公式一分别对所述历史俯仰角数据、所述历史横摆角数据 和所述历史滚动角数据进行处理，获得训练俯仰角数据、训练横摆角数据和 训练滚动角数据；
步骤S23：利用训练俯仰角数据、训练横摆角数据和训练滚动角数据分别 对第一门控循环单元神经网络、第二门控循环单元神经网络和第三门控循环 单元神经网络进行训练，获得所述俯仰角预测模型、所述横摆角预测模型和 所述滚动角预测模型。
所述公式一为：
Figure BDA0002674018810000171
其中，当x为历史俯仰角数据时，x′为训练俯仰角数据，xmin为历史俯仰 角数据中的最小值、xmax为历史俯仰角数据中的最大值；当x为历史横摆角数 据时，x′为训练横摆角数据，xmin为历史横摆角数据中的最小值、xmax为历 史横摆角数据中的最大值；当x为历史滚动角数据时，x′为训练滚动角数据， xmin为历史滚动角数据中的最小值、xmax为历史滚动角数据中的最大值。
需要说明的是，本发明的历史俯仰角数据、历史横摆角数据和历史滚动 角数据可以是连续一段连续时间段内的历史数据，也可以包括多个连续时间 段内各自的历史数据，本发明对连续时间段以及时间段的计时起点均不做限 制，用户可以根据自己的需求进行设定。
另外，本发明的门控循环单元神经网络即为GRU，用户可以根据自己的 需求选用其他类型的初始神经网络模型，并利用训练数据对初始神经网络模 型进行训练，获得预测模型，其中，GRU为一种较优的选择。
本发明为了提高为训练模型的收敛速度，分别采用最大最小值法将所述 历史俯仰角数据、所述历史横摆角数据和所述历史滚动角数据进行归一化至 [0,1]处理，获得获得训练俯仰角数据、训练横摆角数据和训练滚动角数据， 即，采用公式一对所述历史俯仰角数据、所述历史横摆角数据和所述历史滚 动角数据进行处理。
同样，在对历史俯仰角数据、历史横摆角数据和历史滚动角数据进行归 一化至[0,1]处理之前，需要对数据进行有效数据筛选和缺失值填充，筛选和 填充方式参照上述箱型图筛选方法和线性内插填充法，并利用经验模态分解 法对处理后的数据进行分解，获得对应的分解后的训练数据(此处不再赘述)， 并对分解后的数据进行归一化至[0,1]处理。
可以理解的是，利用上述方法训练获得的所述俯仰角预测模型、所述横 摆角预测模型和所述滚动角预测模型对多组俯仰角数据、多组横摆角数据和 多组滚动角数据进行预测时，需要对多组俯仰角数据、多组横摆角数据和多 组滚动角数据进行行归一化至[0,1]处理，然后再利用所述俯仰角预测模型、 所述横摆角预测模型和所述滚动角预测模型对处理后的数据进行预测。
参照图5，图5为本发明盾构机姿态预测装置第一实施例的结构框图，所 述装置包括：
数据获取模块10，用于获取盾构机的盾构姿态数据，所述盾构姿态数据 包括所述盾构机的俯仰角数据、横摆角数据和滚动角数据；
数据获得模块20，用于将所述俯仰角数据、所述横摆角数据和所述滚动 角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获 得预测俯仰角数据、预测横摆角数据和预测滚动角数据；
结果获得模块30，用于根据所述预测俯仰角数据、所述预测横摆角数据 和所述预测滚动角数据，获得所述盾构机的姿态预测结果。
参照具体施工示例，该示例如下：
该施工示例为A项目的示例。完成该实施例所使用的计算机配置和软件 环境为：Intel(R)Core(TM)i5-7200U CPU，16.0GB内存。系统是Windows10 (64位)，程序语言版本为Python3.7.8，集成开发环境为Anaconda包中的 spyder 4.1.4版本，在Keras库中以Tensorflow为后端实现。
步骤1：由于盾构施工在环境较为恶劣的地下进行，因此，原始姿态数据 (即本发明所述的盾构姿态数据：俯仰角数据、横摆角数就和滚动角数据) 可能会因为传感器故障等原因出现一些异常值、缺失值，如果直接用于建模 预测，可能会导致预测精度受到影响。因此，在建模之前，本发明采用箱型 图法进行原始姿态数据的异常值识别，具体的识别方式为：
计算原始姿态数据的上、下四分位数，并基于下式计算该时序数据的上、 下界。
Figure BDA0002674018810000191
其中，U、L分别为上、下四分位数，IQR＝U-L，Upper为上界，Lower 为下界。
将大于上界和小于下界的姿态数据记为该原始姿态数据中的异常值。在 Python语言环境下，基于matplotlib库里pyplot模块里的boxplot()函数实现。
识别出原始姿态数据的异常值之后将其剔除，并采用线性内插的方式完 成缺失值的填充。
步骤2：PLC采集的数据频率为5s/次，由于时间间隔较小时，相应的数 据波动较小。因此，利用过高的采集频率不仅无法提高模型的预测精度，反 而会导致模型计算开销的大幅度增加。因此，本发明将预处理之后的姿态数 据(即本发明的输入数据：输入横摆角数据、输入俯仰角数据和输入滚动角 数据)变换为采集频率为5min/次的姿态数据，共获取数据3000条。
步骤3：对于进行异常值、缺失值预处理之后的姿态数据，为提高数据本 身的规律性，采用经验模态分解的方式对预处理之后的数据进行分解，拆分 成规律性较强的平稳时间序列数据(即，本发明所述的多组横摆角数据、多 组俯仰角数据和多组滚动角数据)，具体的实现方式为：
1)找到预处理之后姿态数据y(t)的全部极值点，通过三次样条插值的方 式获取该序列数据的上下包络曲线BLmax(t)和BLmin(t)，并计算其均值AV(t)，即AV(t)＝BLmax(t)-BLmin(t)；
2)计算盾构姿态时间序列数据y(t)与均值AV(t)之差I(t)，即，I(t)＝y(t)-AV(t)；
3)如果I(t)函数的局部极值数据和过零点数目相差不大于1，且所有时点 的上、下包络线取平均均等于零，则说明其为一个本征模函数，记作I1(t)，反 之，令y(t)＝I(t)，返回步骤1)；
4)将得到的本征模函数I(t)从盾构姿态时间序列数据y(t)中剔除，可得剩 余序列数据R(t)，即R(t)＝y(t)-I1(t)；
5)倘若R(t)为单调函数，则将其作为剩余变量RES，否则，返回步骤1)。
基于上述步骤，可将盾构姿态时间序列数据y(t)分解为如下集合：
y(t)＝{I1(t),I2(t),…,Iv-1(t),RES}
其中，y(t)为预处理之后的盾构姿态时间序列数据，I1(t),I2(t),…,Iv-1(t)分别为基于经验模态分解的v-1个本征模函数，RES为剩余变量。
本发明所提出的经验模态分解在Python语言环境下，基于工具箱pyhht 实现，各个姿态数据分解后的本征模函数个数如表1所示：
表1各姿态参数EMD分解结果
Figure BDA0002674018810000201
步骤4：对于分解后的17个本征模函数和3个剩余变量，在分别构建GRU 预测模型之前，为提高模型的收敛速度，分别采用最大最小值法将数据归一 化至[0,1]，如下式所示：
Figure BDA0002674018810000202
其中，x，x’分别为归一化前、后的样本数据，xmax和xmin为变量x的 最大、最小值。
步骤5：对于归一化之后的17个本征模函数和3个剩余变量，分别建立 基于门控循环单元神经网络的盾构姿态预测模型。全部数据中80％作为训练 集，剩余20％作为测试集。
步骤6：对于划分后的数据集，分别采用滑动窗口的方式进行数据切割， 以得到各预测模型的输入、输出变量。本实例中将窗口宽度设置为10，即用 过去10个时刻的数据对当前时刻的姿态进行预测，如下式所示：
yi(t)＝f(yi(t-1),yi(t-2),…,yi(t-10))
步骤7：当数据划分完成后，即可基于门控循环神经网络进行预测模型的 构建与训练，首先应进行模型初始网络结构的设置。本发明实例中，将初始 网络结构设置为3个GRU层，3个Dropout层(防止过拟合)，1个Dense层。 通过试验法分别确定出姿态预测模型的学习率(lr)、迭代次数(iter)、批处 理大小(epoch)的取值范围，并采用交叉验证网格搜索的方式得到最优超参数 组合。将优化算法设置为Adam算法，GRU层采用Sigmoid函数，Dense层 采用Tanh函数，目标函数为均方误差(MSE)，具体的计算方式为：
Figure BDA0002674018810000211
其中，yi为模型预测值，Yi为工程实测值。
步骤8：采用设置好的初始网络结构和超参数组合，基于训练集进行模型 训练，将训练的标准设定为MSE<7％，当模型训练达到该标准时停止迭代(注： 当模型训练结束，迭代次数达到设定的迭代次数，但尚未达到设定的标准时， 需要增加模型迭代次数)。若模型训练未达到设定标准，需要返回步骤7重新 进行模型网络结构和超参数设定以优化模型训练结果，使其达到MSE<7％的 目标。
步骤9：模型训练完成后，得到的各个模型最优超参数组合如表2所示：
表2GRU预测模型最优超参数组合
Figure BDA0002674018810000212
此时，即可采用测试集进行模型性能测试。将测试集的预测结果和真实 结果进行比对，本实施例中所采用的评价指标为：
Figure RE-GDA0002748827330000242
其中，R2为相关系数，MSE为均方误差，yi为模型预测值，Yi为工程 实测值。均方根误差越小、相关系数越大，说明模型性能越好。
步骤10：将分预测模型的预测结果分别采用最优组合赋权的方式赋予不 同的权重，以得到结果姿态数据(即，结果俯仰角数据、结果横摆角数据和 结果滚动角数据)。具体的实现方式为：
设GRU姿态预测模型得到的预测结果为{y1,y2,..,yv}，采用最优组合赋权得 到结果姿态数据，即：
Figure BDA0002674018810000221
式中，y为各预测模型的结果数据，y1,y2,..,yv分别为v-1个本征模函数和剩 余变量的预测结果，
Figure BDA0002674018810000222
为各个预测模型的预测结果占最终预测结果的权 重。
以使得组合模型的调和平均误差e最小化的思想来实现最优组合赋权， 通过将其转化为求解最优化问题的方式进行计算，即：
min e2
Figure BDA0002674018810000223
s.t.
Figure BDA0002674018810000224
Figure BDA0002674018810000225
i＝1,2,…,v
在Python中编写相应代码通过上述公式计算得，各个模型中分模型赋权 情况如表3所示：
表3各个模型中分模型赋权情况
Figure BDA0002674018810000226
步骤11：为验证本发明公开的准确性和有效性，将其与经典的BP、支持 向量回归(SVR)算法的预测结果进行对比(如表4所示)。结果表明，本发 明公开的基于最优组合赋权的EMD-GRU盾构姿态预测模型的准确性较高， 能够为盾构施工过程中姿态预测提供一种新的思路，进而有助于实现盾构姿 态的主动控制。
表4各模型预测结果对比
Figure BDA0002674018810000227
以上所述仅为本发明的可选实施例，并非因此限制本发明的专利范围， 凡是在本发明的发明构思下，利用本发明说明书及附图内容所作的等效结构 变换，或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范 围内。

Claims (10)

1.一种盾构机姿态预测方法，其特征在于，所述方法包括以下步骤：
获取盾构机的盾构姿态数据，所述盾构姿态数据包括所述盾构机的俯仰角数据、横摆角数据和滚动角数据；
将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据；
根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果。
2.如权利要求1所述的盾构机姿态预测方法，其特征在于，所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
获取盾构机的历史俯仰角数据、历史横摆角数据和历史滚动角数据；
利用公式一分别对所述历史俯仰角数据、所述历史横摆角数据和所述历史滚动角数据进行处理，获得训练俯仰角数据、训练横摆角数据和训练滚动角数据；
利用训练俯仰角数据、训练横摆角数据和训练滚动角数据分别对第一门控循环单元神经网络、第二门控循环单元神经网络和第三门控循环单元神经网络进行训练，获得所述俯仰角预测模型、所述横摆角预测模型和所述滚动角预测模型。
所述公式一为：
Figure FDA0002674018800000011
其中，当x为历史俯仰角数据时，x′为训练俯仰角数据，xmin为历史俯仰角数据中的最小值、xmax为历史俯仰角数据中的最大值；当x为历史横摆角数据时，x′为训练横摆角数据，xmin为历史横摆角数据中的最小值、xmax为历史横摆角数据中的最大值；当x为历史滚动角数据时，x′为训练滚动角数据，xmin为历史滚动角数据中的最小值、xmax为历史滚动角数据中的最大值。
3.如权利要求2所述的盾构机姿态预测方法，其特征在于，所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
利用所述俯仰角数据绘制俯仰角箱形图，并利用所述俯仰角箱形图在所述所述俯仰角数据中筛选有效俯仰角数据；
利用所述横摆角数据绘制横摆角箱形图，并利用所述横摆角箱形图在所述所述横摆角数据中筛选有效横摆角数据；
利用所述滚动角数据绘制滚动角箱形图，并利用所述滚动角箱形图在所述所述滚动角数据中筛选有效滚动角数据；
利用线性内插法分别对所述有效俯仰角数据、所述有效横摆角数据和所述有效滚动角数据进行处理，获得输入俯仰角数据、输入横摆角数据和输入滚动角数据；
所述将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤包括：
将所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据。
4.如权利要求3所述的盾构机姿态预测方法，其特征在于，所述将所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤之前，所述方法还包括：
利用经验模态分解法分别对所述输入俯仰角数据、所述输入横摆角数据和所述输入滚动角数据进行分解，获得具有时间序列的多组俯仰角数据、具有时间序列的多组横摆角数据和具有时间序列的多组滚动角数据；
所述将所述有效俯仰角数据、所述有效横摆角数据和所述有效滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据的步骤包括：
将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得多组预测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据；
所述根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果的步骤包括：
根据所述多组预测俯仰角数据、所述多组预测横摆角数据和所述多组预测滚动角数据，获得所述盾构机的姿态预测结果。
5.如权利要求4所述的盾构机姿态预测方法，其特征在于，所述将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得多组预测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据的步骤包括：
将所述多组俯仰角数据、所述多组横摆角数据和所述多组滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，利用滑动窗口的方式进行滚动预测，获得多组预测俯仰角数据、多组预测横摆角数据和多组预测滚动角数据。
6.如权利要求5所述的盾构机姿态预测方法，其特征在于，所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据和所述多组预测滚动角数据，获得所述盾构机的姿态预测结果的步骤包括：
获取所述多组预测俯仰角数据的俯仰角权重、所述多组预测横摆角数据的横摆角权重和所述多组预测滚动角数据的滚动角权重；
根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所述多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及所述滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据；
根据所述结果俯仰角数据、所述结果横摆角数据和所述结果滚动角数据，获得所述盾构机的姿态预测结果。
7.如权利要求6所述的盾构机姿态预测方法，其特征在于，所述根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所述多组预测滚动角数据、所述俯仰角权重、所述横摆角权以及所述滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据的步骤包括：
利用调和平均误差最小化法分别对所述所述俯仰角权重、所述横摆角权和所述滚动角权重进行处理，获得选定俯仰角权重、选定横摆角权重和选定滚动角权重；
根所述根据所述多组预测俯仰角数据、所述多组预测横摆角数据、所述多组预测滚动角数据、所述选定俯仰角权重、所述选定横摆角权以及所述选定滚动角权重，获得结果俯仰角数据、结果横摆角数据和结果滚动角数据。
8.一种盾构机姿态预测装置，其特征在于，所述装置包括：
获取模块，用于获取盾构机的盾构姿态数据，所述盾构姿态数据包括所述盾构机的俯仰角数据、横摆角数据和滚动角数据；
预测模块，用于将所述俯仰角数据、所述横摆角数据和所述滚动角数据分别输入俯仰角预测模型、横摆角预测模型和滚动角预测模型中，获得预测俯仰角数据、预测横摆角数据和预测滚动角数据；
结果获得模块，用于根据所述预测俯仰角数据、所述预测横摆角数据和所述预测滚动角数据，获得所述盾构机的姿态预测结果。
9.一种终端设备，其特征在于，所述终端设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的盾构机姿态预测程序，所述盾构机姿态预测程序配置为实现如权利要求1至7中任一项所述的盾构机姿态预测方法的步骤。
10.一种存储介质，其特征在于，所述存储介质上存储有盾构机姿态预测程序，所述盾构机姿态预测程序被处理器执行时实现如权利要求1至7中任一项所述的盾构机姿态预测方法的步骤。
CN202010944273.7A 2020-09-09 2020-09-09 盾构机姿态预测方法、装置、终端设备以及存储介质 Active CN112100841B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010944273.7A CN112100841B (zh) 2020-09-09 2020-09-09 盾构机姿态预测方法、装置、终端设备以及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010944273.7A CN112100841B (zh) 2020-09-09 2020-09-09 盾构机姿态预测方法、装置、终端设备以及存储介质

Publications (2)

Publication Number Publication Date
CN112100841A true CN112100841A (zh) 2020-12-18
CN112100841B CN112100841B (zh) 2024-04-19

Family

ID=73752261

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010944273.7A Active CN112100841B (zh) 2020-09-09 2020-09-09 盾构机姿态预测方法、装置、终端设备以及存储介质

Country Status (1)

Country Link
CN (1) CN112100841B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113298220A (zh) * 2021-05-31 2021-08-24 中铁十六局集团北京轨道交通工程建设有限公司 一种基于神经网络优化的盾构机掘进速度预测方法
CN113344256A (zh) * 2021-05-21 2021-09-03 上海隧道工程有限公司 盾构姿态多自由度运动特性预测和控制性能评估系统及方法
CN113837454A (zh) * 2021-09-09 2021-12-24 武汉大学 一种船舶三自由度的混合神经网络模型预测方法及系统
CN114056704A (zh) * 2021-12-03 2022-02-18 惠州市德赛电池有限公司 一种送料纠偏方法、装置及存储介质
CN115450642A (zh) * 2022-11-08 2022-12-09 西南交通大学 一种基于大数据的盾构姿态控制方法及系统
CN116226702A (zh) * 2022-09-09 2023-06-06 武汉中数医疗科技有限公司 一种基于生物电阻抗的甲状腺采样数据识别方法
CN117493837A (zh) * 2024-01-03 2024-02-02 中铁南方投资集团有限公司 一种基于机器学习的盾构掘进机姿态分项预测方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2112713C1 (ru) * 1996-05-16 1998-06-10 Ракетно-космическая корпорация "Энергия" им.С.П.Королева Способ управления космическим аппаратом с помощью реактивных исполнительных органов при поддержании заданной ориентации
JP2006072516A (ja) * 2004-08-31 2006-03-16 Aichi Micro Intelligent Corp 入力支援装置及び、この入力支援装置を組み込んだ携帯電話
CN106126896A (zh) * 2016-06-20 2016-11-16 中国地质大学（武汉） 基于经验模态分解和深度学习的混合模型风速预测方法及系统
WO2017200566A1 (en) * 2016-05-17 2017-11-23 Raytheon Company Gyroscopic attitude control system
CN108664122A (zh) * 2018-04-04 2018-10-16 歌尔股份有限公司 一种姿态预测方法和装置
CN108868805A (zh) * 2018-06-08 2018-11-23 西安电子科技大学 基于统计分析与XGboost结合的盾构纠偏方法
CN109918822A (zh) * 2019-03-15 2019-06-21 中铁工程服务有限公司 基于随机森林回归的盾构掘进参数偏差计算方法及系统
CN110096827A (zh) * 2019-05-09 2019-08-06 中铁工程服务有限公司 一种基于深度神经网络的盾构机参数优化方法
CN110195592A (zh) * 2019-04-30 2019-09-03 华中科技大学 基于混合深度学习的盾构掘进位姿智能预测方法及系统
CN110942210A (zh) * 2019-12-11 2020-03-31 中国铁建重工集团股份有限公司 一种盾构tbm姿态偏移预警方法
CN111365015A (zh) * 2020-03-05 2020-07-03 中建交通建设集团有限公司 一种基于XGBoost的盾构掘进参数特征提取与姿态偏差预测方法
US20200273448A1 (en) * 2019-02-27 2020-08-27 Comcast Cable Communications, Llc Voice Command Detection And Prediction
CN111636891A (zh) * 2020-06-08 2020-09-08 中铁高新工业股份有限公司 盾构姿态实时预测系统和预测模型的构建方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2112713C1 (ru) * 1996-05-16 1998-06-10 Ракетно-космическая корпорация "Энергия" им.С.П.Королева Способ управления космическим аппаратом с помощью реактивных исполнительных органов при поддержании заданной ориентации
JP2006072516A (ja) * 2004-08-31 2006-03-16 Aichi Micro Intelligent Corp 入力支援装置及び、この入力支援装置を組み込んだ携帯電話
WO2017200566A1 (en) * 2016-05-17 2017-11-23 Raytheon Company Gyroscopic attitude control system
CN106126896A (zh) * 2016-06-20 2016-11-16 中国地质大学（武汉） 基于经验模态分解和深度学习的混合模型风速预测方法及系统
CN108664122A (zh) * 2018-04-04 2018-10-16 歌尔股份有限公司 一种姿态预测方法和装置
CN108868805A (zh) * 2018-06-08 2018-11-23 西安电子科技大学 基于统计分析与XGboost结合的盾构纠偏方法
US20200273448A1 (en) * 2019-02-27 2020-08-27 Comcast Cable Communications, Llc Voice Command Detection And Prediction
CN109918822A (zh) * 2019-03-15 2019-06-21 中铁工程服务有限公司 基于随机森林回归的盾构掘进参数偏差计算方法及系统
CN110195592A (zh) * 2019-04-30 2019-09-03 华中科技大学 基于混合深度学习的盾构掘进位姿智能预测方法及系统
CN110096827A (zh) * 2019-05-09 2019-08-06 中铁工程服务有限公司 一种基于深度神经网络的盾构机参数优化方法
CN110942210A (zh) * 2019-12-11 2020-03-31 中国铁建重工集团股份有限公司 一种盾构tbm姿态偏移预警方法
CN111365015A (zh) * 2020-03-05 2020-07-03 中建交通建设集团有限公司 一种基于XGBoost的盾构掘进参数特征提取与姿态偏差预测方法
CN111636891A (zh) * 2020-06-08 2020-09-08 中铁高新工业股份有限公司 盾构姿态实时预测系统和预测模型的构建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HUANG B G: "Online Prediction of Ship Roll Motion Based on a Coarse and Fine Tuning Fixed Grid Wavelet Network", 《OCEAN ENGINEERING》, vol. 160, 31 December 2018 (2018-12-31), pages 425 - 437 *
赵润茂;胡炼;罗锡文;唐灵茂;周浩;杜攀;贺静;: "基于时序的水田平地机俯仰角预测建模与试验", 农业工程学报, no. 11, pages 42 - 47 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113344256A (zh) * 2021-05-21 2021-09-03 上海隧道工程有限公司 盾构姿态多自由度运动特性预测和控制性能评估系统及方法
CN113344256B (zh) * 2021-05-21 2024-03-22 上海隧道工程有限公司 盾构姿态多自由度运动特性预测和控制性能评估系统及方法
CN113298220B (zh) * 2021-05-31 2023-08-04 中铁十六局集团北京轨道交通工程建设有限公司 一种基于神经网络优化的盾构机掘进速度预测方法
CN113298220A (zh) * 2021-05-31 2021-08-24 中铁十六局集团北京轨道交通工程建设有限公司 一种基于神经网络优化的盾构机掘进速度预测方法
CN113837454A (zh) * 2021-09-09 2021-12-24 武汉大学 一种船舶三自由度的混合神经网络模型预测方法及系统
CN113837454B (zh) * 2021-09-09 2024-04-23 武汉大学 一种船舶三自由度的混合神经网络模型预测方法及系统
CN114056704A (zh) * 2021-12-03 2022-02-18 惠州市德赛电池有限公司 一种送料纠偏方法、装置及存储介质
CN116226702A (zh) * 2022-09-09 2023-06-06 武汉中数医疗科技有限公司 一种基于生物电阻抗的甲状腺采样数据识别方法
CN116226702B (zh) * 2022-09-09 2024-04-26 武汉中数医疗科技有限公司 一种基于生物电阻抗的甲状腺采样数据识别方法
CN115450642B (zh) * 2022-11-08 2023-02-17 西南交通大学 一种基于大数据的盾构姿态控制方法及系统
CN115450642A (zh) * 2022-11-08 2022-12-09 西南交通大学 一种基于大数据的盾构姿态控制方法及系统
CN117493837A (zh) * 2024-01-03 2024-02-02 中铁南方投资集团有限公司 一种基于机器学习的盾构掘进机姿态分项预测方法
CN117493837B (zh) * 2024-01-03 2024-03-19 中铁南方投资集团有限公司 一种基于机器学习的盾构掘进机姿态分项预测方法

Also Published As

Publication number Publication date
CN112100841B (zh) 2024-04-19

Similar Documents

Publication Publication Date Title
CN112100841A (zh) 盾构机姿态预测方法、装置、终端设备以及存储介质
Luo et al. Short-term wind speed forecasting via stacked extreme learning machine with generalized correntropy
Liang et al. A novel wind speed prediction strategy based on Bi-LSTM, MOOFADA and transfer learning for centralized control centers
CN112529727A (zh) 基于深度强化学习的微电网储能调度方法及装置及设备
CN109242212A (zh) 一种基于变模式分解与长短记忆网络的风电预测方法
Niu et al. A decomposition-ensemble model with regrouping method and attention-based gated recurrent unit network for energy price prediction
CN113420868A (zh) 一种基于深度强化学习的旅行商问题求解方法及求解系统
Bi et al. Accurate water quality prediction with attention-based bidirectional LSTM and encoder–decoder
Dong et al. A Novel Runoff Prediction Model Based on Support Vector Machine and Gate Recurrent unit with Secondary Mode Decomposition
CN111210051A (zh) 一种用户用电行为预测方法及系统
CN114663710A (zh) 轨迹识别方法、装置、设备及存储介质
CN114969360A (zh) 文本信息处理方法、装置、设备及存储介质
CN114971053A (zh) 低压台区网络线损率在线预测模型的训练方法及装置
CN111931842A (zh) 盾构机导向预测方法、模型训练方法、装置以及设备
CN113516315B (zh) 风力发电功率区间预测方法、设备及介质
Mai Efficient multimedia information mining framework based on deep learning and self-organizing model
CN116307179A (zh) 一种基于vmd-pca-gru的电力负荷预测方法及装置
CN115718536B (zh) 一种调频方法、装置、电子设备及可读存储介质
CN113886688A (zh) 对象的关联关系预测方法、装置、终端设备以及存储介质
CN108965039A (zh) 一种基于流量监控的告警方法以及电子设备
Lv et al. Short-term load forecasting method of IES based on RLA neural network with dual signal processing
Yang et al. A Novel Wireless Propagation Model Based on Bi-LSTM Algorithm
Chen et al. A Hybrid Model for Time Series Prediction Using Adaptive Variational Mode Decomposition
Dong et al. Short-Term Electrical Load Forecasting Based on Whale Optimized Variational Mode Decomposition and DCNN-BiLSTM
CN117939518A (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