CN109488321A - 一种盾构刀盘扭矩确定方法及系统 - Google Patents

一种盾构刀盘扭矩确定方法及系统 Download PDF

Info

Publication number
CN109488321A
CN109488321A CN201910004164.4A CN201910004164A CN109488321A CN 109488321 A CN109488321 A CN 109488321A CN 201910004164 A CN201910004164 A CN 201910004164A CN 109488321 A CN109488321 A CN 109488321A
Authority
CN
China
Prior art keywords
cutter head
unit
head torque
characteristic
shield 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
Application number
CN201910004164.4A
Other languages
English (en)
Other versions
CN109488321B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201910004164.4A priority Critical patent/CN109488321B/zh
Publication of CN109488321A publication Critical patent/CN109488321A/zh
Application granted granted Critical
Publication of CN109488321B publication Critical patent/CN109488321B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18—Special adaptations of signalling or alarm devices
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Primary Health Care (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Animal Husbandry (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明公开一种盾构刀盘扭矩确定方法及系统。方法包括：确定影响盾构刀盘扭矩的参数；根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型；根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。采用本发明的方法或系统能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。

Description

一种盾构刀盘扭矩确定方法及系统
技术领域
本发明涉及盾构掘进领域，特别是涉及一种盾构刀盘扭矩确定方法及系统。
背景技术
盾构是一种自动化程度高、施工速度快的隧道施工装备，主要用于软土地层的隧道挖掘。其刀盘扭矩是重要的性能参数，用于维持刀盘不断旋转切削前方土层。目前对盾构刀盘扭矩进行预测的研究，主要分为基于软土地层的力学分析的理论建模与基于工程数据分析进行预测。其中，在基于工程数据分析进行预测的这类方法中，通过量纲分析进行计算的方法普适性高。然而，已有基于量纲分析的方法考虑的参数众多，从而导致模型结构复杂，不利于软土地层中盾构掘进的实际应用需求。
发明内容
本发明的目的是提供一种盾构刀盘扭矩确定方法及系统，能够筛选出影响盾构刀盘扭矩的主要参数，能够帮助盾构技术管理人员及时调整盾构刀盘扭矩的大小。
为实现上述目的，本发明提供了如下方案：
一种盾构刀盘扭矩确定方法，包括：
确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
可选的，所述根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩，具体包括：
根据公式分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；T为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
可选的，所述根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果，具体包括：
将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法采用公式得到无量纲量的系数；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩T；λ为惩罚权重。
可选的，所述将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型，具体包括：
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
一种盾构刀盘扭矩确定系统，包括：
参数确定模块，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
无量纲参数确定模块，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
系数识别结果确定模块，用于根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
模型建立模块，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
盾构刀盘扭矩计算模块，用于根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
可选的，所述无量纲参数确定模块，具体包括：
无量纲量参数确定单元，用于根据公式 分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
无量纲化刀盘扭矩确定单元，用于根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；T为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
可选的，所述系数识别结果确定模块，具体包括：
系数识别结果确定单元，用于将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩作为因变量，利用LASSO算法采用公式得到无量纲量的系数α1,α2,α3,α4,α5；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩λ为惩罚权重。
可选的，所述模型建立模块，具体包括：
模型建立单元，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供一种盾构刀盘扭矩确定方法，方法包括：确定影响刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。本发明能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。在基于LASSO算法的基础上，本发明不仅融入了基本的力学规律，而且通过LASSO算法剔除了无关的影响参数，从而为盾构技术管理人员及时调整刀盘扭矩大小提供了有效的决策支持。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例盾构刀盘扭矩确定方法流程图；
图2为本发明实施例盾构刀盘扭矩确定系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种盾构刀盘扭矩确定方法及系统，能够筛选出影响盾构刀盘扭矩的主要参数，能够帮助盾构技术管理人员及时调整刀盘扭矩的大小。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例盾构刀盘扭矩确定方法流程图。如图1所示，一种盾构刀盘扭矩确定方法，包括：
步骤101：确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
步骤102：根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
步骤103：根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
步骤104：将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
步骤105：根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
步骤102，具体包括：
根据公式分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
步骤103，具体包括：
将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法采用公式得到无量纲量的系数α1,α2,α3,α4,α5；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩T；λ为惩罚权重。
步骤104，具体包括：
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
在基于套索算法(LASSO算法)的基础上，本发明不仅融入了基本的力学规律，而且通过LASSO算法剔除了无关的影响参数，从而为盾构技术管理人员及时调整刀盘扭矩大小提供了有效的决策支持，能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。
图2为本发明实施例盾构刀盘扭矩确定系统结构图。如图2所示，一种盾构刀盘扭矩确定系统，包括：
参数确定模块201，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
无量纲参数确定模块202，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
系数识别结果确定模块203，用于根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
模型建立模块204，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型；
盾构刀盘扭矩计算模块205，用于根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
所述根据无量纲参数确定模块202，具体包括：
无量纲量参数确定单元，用于根据公式 分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
无量纲化刀盘扭矩确定单元，用于根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；T为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
所述系数识别结果确定模块203，具体包括：
系数识别结果确定单元，用于将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩作为因变量，利用LASSO算法采用公式得到无量纲量的系数α1,α2,α3,α4,α5；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，即无量纲化刀盘扭矩T；λ为惩罚权重。
所述模型建立模块204，具体包括：
模型建立单元，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
具体实施例1：
下面结合实施例对本发明的盾构刀盘扭矩的计算方法做出详细说明。具体步骤如下：
1、列出影响盾构刀盘扭矩T的参数：pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；T为盾构刀盘扭矩，单位：kNm。
2、根据下列公式(1)-(6)，分别计算无量纲量及无量纲化刀盘扭矩大小。其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
3、将按公式(1)-(6)式计算得到的值汇总，并将(π1，π2，π3，π4，π5)这五个无量纲量作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法，得到无量纲量的系数α1～α5的识别计算式如下：
α1＝0,α2＝0,α3＝6.02×10-5,α4＝0,α5＝0。
4、将各个无量纲量系数识别结果α1～α5与五个无量纲量π1～π5分别相乘，得到盾构刀盘扭矩T的计算模型：如下式所示：
T＝6.02*10-5EHD2 (7)
将如下表所示的工程数据：
E D ω p<sub>e</sub> γ H W ν
67650 6.34 1.22 168 40 9.6 160.59 0.03
带入式(7)可得盾构刀盘扭矩T的大小如下：
T＝1571.50(kNm)。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种盾构刀盘扭矩确定方法，其特征在于，包括：
确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
2.根据权利要求1所述的盾构刀盘扭矩确定方法，其特征在于，所述根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩，具体包括：
根据公式分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
3.根据权利要求2所述的盾构刀盘扭矩确定方法，其特征在于，所述根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果，具体包括：
将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩作为因变量，利用LASSO算法采用公式得到无量纲量的系数；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩λ为惩罚权重。
4.根据权利要求3所述的盾构刀盘扭矩确定方法，其特征在于，所述将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型，具体包括：
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
5.一种盾构刀盘扭矩确定系统，其特征在于，包括：
参数确定模块，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
无量纲参数确定模块，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
系数识别结果确定模块，用于根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
模型建立模块，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
盾构刀盘扭矩计算模块，用于根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
6.根据权利要求5所述的盾构刀盘扭矩确定系统，其特征在于，所述无量纲参数确定模块，具体包括：
无量纲量参数确定单元，用于根据公式 分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
无量纲化刀盘扭矩确定单元，用于根据公式确定盾构掘进过程中的无量纲化刀盘扭矩
其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
7.根据权利要求6所述的盾构刀盘扭矩确定系统，其特征在于，所述系数识别结果确定模块，具体包括：
系数识别结果确定单元，用于将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩作为因变量，利用LASSO算法采用公式得到无量纲量的系数α1,α2,α3,α4,α5；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩λ为惩罚权重。
8.根据权利要求7所述的盾构刀盘扭矩确定系统，其特征在于，所述模型建立模块，具体包括：
模型建立单元，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
CN201910004164.4A 2019-01-03 2019-01-03 一种盾构刀盘扭矩确定方法及系统 Active CN109488321B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910004164.4A CN109488321B (zh) 2019-01-03 2019-01-03 一种盾构刀盘扭矩确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910004164.4A CN109488321B (zh) 2019-01-03 2019-01-03 一种盾构刀盘扭矩确定方法及系统

Publications (2)

Publication Number Publication Date
CN109488321A true CN109488321A (zh) 2019-03-19
CN109488321B CN109488321B (zh) 2019-11-29

Family

ID=65713791

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910004164.4A Active CN109488321B (zh) 2019-01-03 2019-01-03 一种盾构刀盘扭矩确定方法及系统

Country Status (1)

Country Link
CN (1) CN109488321B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113312763A (zh) * 2021-05-19 2021-08-27 上海交通大学 一种采用无量纲参数的离合器压力分布获取方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102262712A (zh) * 2011-08-18 2011-11-30 天津大学 上下分层地质条件下盾构刀盘正面载荷的计算方法
CN107992447A (zh) * 2017-12-13 2018-05-04 电子科技大学 一种应用于河流水位预测数据的特征选择分解方法
CN108388689A (zh) * 2018-01-16 2018-08-10 天津大学 一种基于量纲分析的盾构刀盘扭矩计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102262712A (zh) * 2011-08-18 2011-11-30 天津大学 上下分层地质条件下盾构刀盘正面载荷的计算方法
CN107992447A (zh) * 2017-12-13 2018-05-04 电子科技大学 一种应用于河流水位预测数据的特征选择分解方法
CN108388689A (zh) * 2018-01-16 2018-08-10 天津大学 一种基于量纲分析的盾构刀盘扭矩计算方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113312763A (zh) * 2021-05-19 2021-08-27 上海交通大学 一种采用无量纲参数的离合器压力分布获取方法和装置

Also Published As

Publication number Publication date
CN109488321B (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
CN107577862B (zh) 一种tbm在掘岩体状态实时感知系统和方法
CN109685378A (zh) 一种基于数据挖掘的tbm施工围岩可掘性分级方法
CN108089227B (zh) 新的基于三维地震资料的地层孔隙压力预测方法
CN106777520A (zh) 一种埋入式抗滑桩的设计方法
US20150142406A1 (en) Directional drill ahead simulator: directional wellbore prediction using bha and bit models
CN113223158B (zh) 一种基于爆破岩体三维实体模型的炮孔装药量计算方法
CN106526329A (zh) 一种测量地表到深层大地电阻率的方法及装置
CN109488321B (zh) 一种盾构刀盘扭矩确定方法及系统
WO2012021292A1 (en) Reservoir upscaling method with preserved transmissibility
US9951560B2 (en) Axial motion drill bit model
US10012025B2 (en) Lateral motion drill bit model
CN110197013B (zh) 基于Morphing的河床基岩面建模方法
US11112514B2 (en) Systems and methods for computed resource hydrocarbon reservoir simulation and development
CN109736819A (zh) 一种盾构掘进总推力确定方法及系统
CN113158561B (zh) 适用于多种岩体条件的tbm操作参数优化方法及系统
CN116291406A (zh) 一种海上疏松砂岩油藏大孔道识别方法和系统
CN109711079A (zh) 一种tbm掘进总推力确定方法及系统
CN111611705B (zh) 一种盾构掘进地层参数的取值方法
CN109766621A (zh) 一种tbm刀盘扭矩确定方法及系统
CN109138969B (zh) 一种钻进状态变量的预测方法、设备及存储设备
CN113971334A (zh) 基于三维应力场数值模拟预测裂缝有效性的方法及装置
CN113625339B (zh) 正断层地层厚度确定方法及装置
CN113032959B (zh) 基于施工参数的盾构机掘进的安全状态确定方法、装置
CN111911173B (zh) 一种隧道开挖方法
CN118278265B (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