本发明公开一种盾构刀盘扭矩确定方法及系统。方法包括：确定影响盾构刀盘扭矩的参数；根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型；根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。采用本发明的方法或系统能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。
1.一种盾构刀盘扭矩确定方法，其特征在于，包括：确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩，具体包括：根据公式 分别确定盾构掘进过程中的
根据公式 确定盾构掘进过程中的无量纲化刀盘扭矩其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min； 为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kN·m；
根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果，具体包括：将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩 作为因变量，利用LASSO算法采用公式 得到无量纲量的系数；
其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩 λ为惩罚权重；
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型，具体包括：
将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
其中，T为盾构刀盘扭矩，单位：kN·m；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；
2.一种盾构刀盘扭矩确定系统，其特征在于，包括：参数确定模块，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
无量纲参数确定模块，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
系数识别结果确定模块，用于根据所述无量纲量采用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为盾构刀盘扭矩，单位：kN·m；
系数识别结果确定单元，用于将所述无量纲量π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为盾构刀盘扭矩，单位：kN·m；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压3
力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m ；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
一种盾构刀盘扭矩确定方法及系统
技术领域
[0001] 本发明涉及盾构掘进领域，特别是涉及一种盾构刀盘扭矩确定方法及系统。
背景技术
[0002] 盾构是一种自动化程度高、施工速度快的隧道施工装备，主要用于软土地层的隧道挖掘。其刀盘扭矩是重要的性能参数，用于维持刀盘不断旋转切削前方土层。目前对盾构刀盘扭矩进行预测的研究，主要分为基于软土地层的力学分析的理论建模与基于工程数据分析进行预测。其中，在基于工程数据分析进行预测的这类方法中，通过量纲分析进行计算的方法普适性高。然而，已有基于量纲分析的方法考虑的参数众多，从而导致模型结构复杂，不利于软土地层中盾构掘进的实际应用需求。
发明内容
[0003] 本发明的目的是提供一种盾构刀盘扭矩确定方法及系统，能够筛选出影响盾构刀盘扭矩的主要参数，能够帮助盾构技术管理人员及时调整盾构刀盘扭矩的大小。
[0004] 为实现上述目的，本发明提供了如下方案：
[0005] 一种盾构刀盘扭矩确定方法，包括：
[0006] 确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
[0007] 根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
[0008] 根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
[0009] 将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
[0010] 根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
[0011] 可选的，所述根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩，具体包括：
[0012] 根据公式 分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
[0013] 根据公式 确定盾构掘进过程中的无量纲化刀盘扭矩
[0014] 其中，π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。
[0015] 可选的，所述根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果，具体包括：
[0016] 将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法采用公式 得到无量纲量的系数；
[0017] 其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩T；λ为惩罚权重。
[0018] 可选的，所述将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型，具体包括：
[0019] 将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
[0020] 其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
[0021] 一种盾构刀盘扭矩确定系统，包括：
[0022] 参数确定模块，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
[0023] 无量纲参数确定模块，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
[0024] 系数识别结果确定模块，用于根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
[0025] 模型建立模块，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
[0026] 盾构刀盘扭矩计算模块，用于根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
[0027] 可选的，所述无量纲参数确定模块，具体包括：
[0028] 无量纲量参数确定单元，用于根据公式分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
[0029] 无量纲化刀盘扭矩确定单元，用于根据公式 确定盾构掘进过程中的无量纲化刀盘扭矩
[0030] 其中，π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。
[0031] 可选的，所述系数识别结果确定模块，具体包括：
[0032] 系数识别结果确定单元，用于将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩 作为因变量，利用LASSO算法采用公式 得到无量纲量的系数α1,α2,α3,α4,α5；
[0033] 其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩 λ为惩罚权重。
[0034] 可选的，所述模型建立模块，具体包括：
[0035] 模型建立单元，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
[0036] 其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
[0037] 根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供一种盾构刀盘扭矩确定方法，方法包括：确定影响刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。本发明能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。在基于LASSO算法的基础上，本发明不仅融入了基本的力学规律，而且通过LASSO算法剔除了无关的影响参数，从而为盾构技术管理人员及时调整刀盘扭矩大小提供了有效的决策支持。
附图说明
[0038] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0039] 图1为本发明实施例盾构刀盘扭矩确定方法流程图；
[0040] 图2为本发明实施例盾构刀盘扭矩确定系统结构图。
具体实施方式
[0041] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0042] 本发明的目的是提供一种盾构刀盘扭矩确定方法及系统，能够筛选出影响盾构刀盘扭矩的主要参数，能够帮助盾构技术管理人员及时调整刀盘扭矩的大小。
[0043] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0044] 图1为本发明实施例盾构刀盘扭矩确定方法流程图。如图1所示，一种盾构刀盘扭矩确定方法，包括：
[0045] 步骤101：确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
[0046] 步骤102：根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
[0047] 步骤103：根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
[0048] 步骤104：将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到刀盘扭矩计算模型；
[0049] 步骤105：根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
[0050] 步骤102，具体包括：
[0051] 根据公式 分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
[0052] 根据公式 确定盾构掘进过程中的无量纲化刀盘扭矩
[0053] 其中，π1，π2，π3，π4，π5为五个无量纲量；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min； 为无量纲化刀盘扭矩；T为盾构刀盘扭矩，单位：kNm。
[0054] 步骤103，具体包括：
[0055] 将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法采用公式 得到无量纲量的系数α1,α2,α3,α4,α5；
[0056] 其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，y为无量纲化刀盘扭矩T；λ为惩罚权重。
[0057] 步骤104，具体包括：
[0058] 将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
[0059] 其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
[0060] 在基于套索算法(LASSO算法)的基础上，本发明不仅融入了基本的力学规律，而且通过LASSO算法剔除了无关的影响参数，从而为盾构技术管理人员及时调整刀盘扭矩大小提供了有效的决策支持，能够实现盾构装备掘进过程中的刀盘扭矩快速准确的计算。
[0061] 图2为本发明实施例盾构刀盘扭矩确定系统结构图。如图2所示，一种盾构刀盘扭矩确定系统，包括：
[0062] 参数确定模块201，用于确定影响盾构刀盘扭矩的参数，所述参数包括土体弹性模量、土体容量、隧道埋深、岩土承载力、掘进速度、土舱压力、盾构直径和刀盘转速；
[0063] 无量纲参数确定模块202，用于根据所述参数确定盾构掘进过程中的无量纲量和无量纲化刀盘扭矩；
[0064] 系数识别结果确定模块203，用于根据所述无量纲量采用LASSO算法，得到所述无量纲量的系数识别结果；
[0065] 模型建立模块204，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型；
[0066] 盾构刀盘扭矩计算模块205，用于根据所述刀盘扭矩计算模型确定盾构刀盘扭矩。
[0067] 所述根据无量纲参数确定模块202，具体包括：
[0068] 无量纲量参数确定单元，用于根据公式分别确定盾构掘进过程中的无量纲量π1，π2，π3，π4，π5；
[0069] 无量纲化刀盘扭矩确定单元，用于根据公式 确定盾构掘进过程中的无量纲化刀盘扭矩
[0070] 其中，π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。
[0071] 所述系数识别结果确定模块203，具体包括：
[0072] 系数识别结果确定单元，用于将所述无量纲量π1，π2，π3，π4，π5作为自变量，无量纲化刀盘扭矩 作为因变量，利用LASSO算法采用公式 得到无量纲量的系数α1,α2,α3,α4,α5；
[0073] 其中，n为样本量；X为自变量，包括π1，π2，π3，π4，π5；α为识别系数，包括α1,α2,α3,α4,α5；y为因变量，即无量纲化刀盘扭矩T；λ为惩罚权重。
[0074] 所述模型建立模块204，具体包括：
[0075] 模型建立单元，用于将各个无量纲量系数识别结果与对应的所述无量纲量相乘，得到盾构刀盘扭矩计算模型
[0076] 其中，T为盾构刀盘扭矩，单位：kNm；α1,α2,α3,α4,α5为无量纲量的系数；pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min。
[0077] 具体实施例1：
[0078] 下面结合实施例对本发明的盾构刀盘扭矩的计算方法做出详细说明。具体步骤如下：
[0079] 1、列出影响盾构刀盘扭矩T的参数：pe为土舱压力，单位：kpa；E为土体弹性模量，单位：kpa；γ为土体容重，单位：kN/m3；D为盾构直径，单位：m；H为隧道埋深，单位：m；W为岩土承载力，单位：kpa；v为掘进速度，单位：m/s；ω为刀盘转速，单位：r/min；T为盾构刀盘扭矩，单位：kNm。
[0080] 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。
[0081]
[0082]
[0083]
[0084]
[0085]
[0086]
[0087] 3、将按公式(1)-(6)式计算得到的值汇总，并将(π1，π2，π3，π4，π5)这五个无量纲量作为自变量，无量纲化刀盘扭矩T作为因变量，利用LASSO算法，得到无量纲量的系数α1～α5的识别计算式如下：
[0088] α1＝0,α2＝0,α3＝6.02×10-5,α4＝0,α5＝0。
[0089] 4、将各个无量纲量系数识别结果α1～α5与五个无量纲量π1～π5分别相乘，得到盾构刀盘扭矩T的计算模型：如下式所示：
[0090] T＝6.02*10-5EHD2 (7)
[0091] 将如下表所示的工程数据：
[0092]E D ω pe γ H W ν
67650 6.34 1.22 168 40 9.6 160.59 0.03
[0093] 带入式(7)可得盾构刀盘扭矩T的大小如下：
[0094] T＝1571.50(kNm)。
[0095] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0096] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
