CN109736819B - 一种盾构掘进总推力确定方法及系统 - Google Patents
一种盾构掘进总推力确定方法及系统 Download PDFInfo
- Publication number
- CN109736819B CN109736819B CN201910004323.0A CN201910004323A CN109736819B CN 109736819 B CN109736819 B CN 109736819B CN 201910004323 A CN201910004323 A CN 201910004323A CN 109736819 B CN109736819 B CN 109736819B
- Authority
- CN
- China
- Prior art keywords
- unit
- dimensionless
- thrust
- total
- 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.)
- Active
Links
Images
Landscapes
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明公开一种盾构掘进总推力确定方法及系统。方法包括:确定影响盾构掘进总推力的参数;根据所述参数确定盾构掘进过程中的无量纲量和无量纲化总推力;根据所述无量纲量采用LASSO算法,得到所述无量纲量的系数识别结果;将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到掘进总推力计算模型;根据所述掘进总推力计算模型确定盾构掘进总推力。采用本发明的方法或系统能够实现盾构掘进过程中的总推力快速准确的计算。
Description
技术领域
本发明涉及盾构掘进领域,特别是涉及一种盾构掘进总推力确定方法及系统。
背景技术
盾构是一种隧道掘进的专用工程机械,主要用于软土地层的隧道挖掘。其掘进总推力是一种重要的性能参数,以保证装备正常前进,是贯穿装备运行始终的核心参量。目前针对盾构掘进总推力的计算方法,主要分为基于软土地层的力学分析的理论建模与基于工程数据分析进行预测。其中,在基于工程数据分析进行预测的这类方法中,通过量纲分析进行计算的方法普适性高。然而,已有基于量纲分析的方法考虑的参数众多,从而导致模型结构复杂,不利于工程应用。
发明内容
本发明的目的是提供一种盾构掘进总推力确定方法及系统,能够筛选出影响盾构掘进总推力的主要参数,能够帮助盾构技术管理人员及时调整盾构掘进总推力的大小。
为实现上述目的,本发明提供了如下方案:
一种盾构掘进总推力确定方法,包括:
确定影响盾构掘进总推力的参数,所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速;
根据所述参数确定盾构掘进过程中的无量纲量和无量纲化总推力;
根据所述无量纲量采用LASSO算法,得到所述无量纲量的系数识别结果;
将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到掘进总推力计算模型;
根据所述掘进总推力计算模型确定盾构掘进总推力。
可选的,所述根据所述参数确定盾构掘进过程中的无量纲量和无量纲化总推力,具体包括:
其中,π1,π2,π3,π4,π5为五个无量纲量;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;为无量纲化总推力;F为掘进总推力,单位:kN。
可选的,所述根据所述无量纲量采用LASSO算法,得到所述无量纲量的系数识别结果,具体包括:
可选的,所述将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到盾构掘进总推力计算模型,具体包括:
其中,F为盾构掘进总推力,单位:kN;α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为五个无量纲量;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;为无量纲化总推力;F为盾构掘进总推力,单位:kN。
可选的,所述系数识别结果确定模块,具体包括:
可选的,所述模型建立模块,具体包括:
其中,F为盾构掘进总推力,单位:kN;α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为五个无量纲量;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;为无量纲化总推力;F为盾构掘进总推力,单位:kN。
步骤103,具体包括:
步骤104,具体包括:
其中,F为盾构掘进总推力,单位:kN;α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;为无量纲化总推力;F为盾构掘进总推力,单位:kN。
所述系数识别结果确定模块203,具体包括:
所述模型建立模块204,具体包括:
其中,F为盾构掘进总推力,单位:kN;α1,α2,α3,α4,α5为无量纲量的系数;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min。
具体实施例:
下面结合实施例对本发明的盾构总推力的计算方法做出详细说明。具体步骤如下:
1、列出影响盾构掘进总推力F的参数:pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;F为掘进总推力,单位:kN。
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;为无量纲化总推力;F为盾构掘进总推力,单位:kN。
3、将按公式(1)-(6)式计算得到的值汇总,并将(π1,π2,π3,π4,π5)这五个无量纲量作为自变量,无量纲化总推力作为因变量,利用LASSO算法,得到无量纲量的系数α1~α5的识别计算式如下:
α1=0,α2=4.93×10-1,α3=2.22×10-3,α4=0,α5=0.
4、将各个无量纲量系数识别结果α1~α5与五个无量纲量π1~π5分别相乘,得到掘进总推力F的计算模型:如下式所示:
F=4.93*10-1γD3+2.22*10-3EDH (7)
将如下表所示的工程数据:
E | D | ω | p<sub>e</sub> | γ | H | W | ν |
67650 | 6.34 | 1.22 | 168 | 40 | 9.6 | 160.59 | 0.03 |
带入式(7)可得盾构掘进总推力F的大小如下:
F=14166.18(kN)。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种盾构掘进总推力确定方法,其特征在于,包括:
确定影响盾构掘进总推力的参数,所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速;
根据所述参数确定盾构掘进过程中的无量纲量和无量纲化总推力,具体包括:
其中,π1,π2,π3,π4,π5为五个无量纲量;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;为无量纲化总推力;F为盾构掘进总推力,单位:kN;
根据所述无量纲量采用LASSO算法,得到所述无量纲量的系数识别结果,具体包括:
将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到掘进总推力计算模型;
根据所述掘进总推力计算模型确定盾构掘进总推力。
3.一种盾构掘进总推力确定系统,其特征在于,包括:
参数确定模块,用于确定影响盾构掘进总推力的参数,所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速;
无量纲参数确定模块,用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化总推力;
系数识别结果确定模块,用于根据所述无量纲量采用LASSO算法,得到所述无量纲量的系数识别结果;
模型建立模块,用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到掘进总推力计算模型;
盾构总推力计算模块,用于根据所述掘进总推力计算模型确定盾构掘进总推力;
所述无量纲参数确定模块,具体包括:
其中,π1,π2,π3,π4,π5为五个无量纲量;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min;为无量纲化总推力;F为盾构掘进总推力,单位:kN;
所述系数识别结果确定模块,具体包括:
4.根据权利要求3所述的盾构掘进总推力确定系统,其特征在于,所述模型建立模块,具体包括:
模型建立单元,用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘,得到盾构掘进总推力计算模型
其中,F为盾构掘进总推力,单位:kN;α1,α2,α3,α4,α5为无量纲量的系数;pe为土舱压力,单位:kpa;E为土体弹性模量,单位:kpa;γ为土体容重,单位:kN/m3;D为盾构直径,单位:m;H为隧道埋深,单位:m;W为岩土承载力,单位:kpa;v为掘进速度,单位:m/s;ω为刀盘转速,单位:r/min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910004323.0A CN109736819B (zh) | 2019-01-03 | 2019-01-03 | 一种盾构掘进总推力确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910004323.0A CN109736819B (zh) | 2019-01-03 | 2019-01-03 | 一种盾构掘进总推力确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109736819A CN109736819A (zh) | 2019-05-10 |
CN109736819B true CN109736819B (zh) | 2020-02-14 |
Family
ID=66363295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910004323.0A Active CN109736819B (zh) | 2019-01-03 | 2019-01-03 | 一种盾构掘进总推力确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109736819B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111894675B (zh) * | 2020-08-06 | 2022-05-10 | 中铁九局集团有限公司 | 一种地铁隧道塌方推管法抢险推力预测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10246091A (ja) * | 1997-03-05 | 1998-09-14 | Hitachi Constr Mach Co Ltd | シールド工法及びシールド掘削機 |
CN103670383B (zh) * | 2012-09-17 | 2016-08-03 | 中国石油天然气股份有限公司 | 一种识别泥页岩油藏有效储层的方法及设备 |
CN205446412U (zh) * | 2015-12-29 | 2016-08-10 | 瓦戭店轴承集团有限责任公司 | 一种盾构机主轴轴承的保持架 |
CN107992447B (zh) * | 2017-12-13 | 2019-12-17 | 电子科技大学 | 一种应用于河流水位预测数据的特征选择分解方法 |
CN108446418B (zh) * | 2018-01-16 | 2022-03-18 | 天津大学 | 一种基于量纲分析的盾构掘进总推力计算方法 |
CN108763673B (zh) * | 2018-05-16 | 2021-11-23 | 广东省科学院生态环境与土壤研究所 | 基于lasso回归的土地利用变化驱动力筛选方法和装置 |
-
2019
- 2019-01-03 CN CN201910004323.0A patent/CN109736819B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109736819A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108446418B (zh) | 一种基于量纲分析的盾构掘进总推力计算方法 | |
CN109685378A (zh) | 一种基于数据挖掘的tbm施工围岩可掘性分级方法 | |
CN113361043B (zh) | 一种硬岩隧道掘进机刀盘比能的预测方法及系统 | |
CN108388689B (zh) | 一种基于量纲分析的盾构刀盘扭矩计算方法 | |
CN109736819B (zh) | 一种盾构掘进总推力确定方法及系统 | |
CN104453879A (zh) | 钻前压力的预测方法 | |
CN114429231A (zh) | 页岩气井产量预测方法、系统、存储介质以及电子设备 | |
CN115238533A (zh) | 边坡块体稳定性评价方法、系统、设备及可读存储介质 | |
CN118037982A (zh) | 一种岩体地质结构的三维可视化建模方法及其相关设备 | |
CN103477248A (zh) | 用于估算碳氢化合物生产区域中至少一个参数值的计算方法,以在该区域规划和实施作业 | |
CN109711079B (zh) | 一种tbm掘进总推力确定方法及系统 | |
CN112377182A (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
CN109488321B (zh) | 一种盾构刀盘扭矩确定方法及系统 | |
CN109766621B (zh) | 一种tbm刀盘扭矩确定方法及系统 | |
CN116305462A (zh) | 一种考虑爆破与富水扰动影响的岩石Hoek-Brown准则参数确定方法 | |
CN116291406A (zh) | 一种海上疏松砂岩油藏大孔道识别方法和系统 | |
CN116402339A (zh) | 一种盾构隧道施工风险等级评估方法、系统、设备及介质 | |
CN113836617B (zh) | 基于时间序列模型的掘进面前方裂隙类比预测方法及系统 | |
CN112419493B (zh) | 页岩储层三维属性模型建立方法及装置 | |
CN104462792B (zh) | 一种测井数据岩性层数值归约方法 | |
CN115110962B (zh) | 一种不同复合比复合地层盾构掘进参数预测方法及系统 | |
CN108805339B (zh) | 煤矿工作面瓦斯涌出量预测方法、存储介质及电子设备 | |
CN115126474A (zh) | 基于水驱的井间渗流通道表征方法、装置及存储介质 | |
CN115130047A (zh) | 基于化学驱的井间渗流通道表征方法、装置及存储介质 | |
CN116485225B (zh) | 基于钻进参数的施工阶段围岩bq值自动获取方法及系统 |
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 |