陈航1，王喆2，芦杰2，李军3，许鸿运3
(1. 中铁隧道局集团有限公司，广东 广州 511458；2. 绍兴市交通工程质量安全监督站，浙江 绍兴 312000；3. 中南大学 土木工程学院，湖南 长沙 410075)
摘 要：铣挖法是一种采用悬臂式掘进机截割隧道掌子面围岩的隧道开挖工法，不同的铣挖步序对围岩稳定性的影响也存在差异。针对此问题，考虑到悬臂式掘进机的设备条件并参考国内隧道工程已经采用的铣挖步序，确定4种铣挖步序方案，使用有限元分析软件abaqus对4种不同的铣挖步序方案进行数值模拟，再运用灰色关联分析法，对4种方案模拟的结果进行综合分析评价。研究结果表明：先从中间竖直开槽，再水平开挖的铣挖步序对围岩扰动最小。
关键词：隧道；铣挖步序；abaqus；灰色关联分析法
随着我国基础设施建设不断推进，我国隧道里程总长度位于世界前列。截至2016年底统计，我国公路隧道总里程达496.63万公里[1]，我国铁路隧道总里程达1.412万公里[2]。隧道作为公路或铁路工程的控制性工程，隧道的稳定性更是重中之重。对于隧道的施工开挖方法，常用的钻孔爆破法对施工周围的环境影响大，施工过程中超、欠挖问题严 重[3]，铣挖法是一种采用悬臂式掘进机截割隧道掌子面围岩的隧道施工工艺，铣挖法具备对围岩的扰动小且超、欠挖问题控制较好，施工过程相对更安全，开挖轮廓圆顺等优点。对隧道施工开挖后地表沉降及围岩的扰动变形控制要求高的工程项目，近年来逐渐会采用铣挖法完成施工[4]，且施工效率尚可，如杭州紫之隧道铣挖施工速度可达3 m/d[5]，或采用铣挖法和其他工法联合施工。目前，随着铣挖法开挖隧道的工程经验不断积累和国内外学者对铣挖法施工工艺的研究，已经在各方面取得了一定的经验和成果。Jagieaao等[6]对单轴抗压强度140 MPa的岩石完成切割试验。王渭明等[7]分析优选了半煤岩巷道快速综掘截割顺序；李国等[8]综合考虑安全和效率2个方面，研究隧道上台阶分左右2次铣挖时，通过对2种工况7种方案数值模拟的结果综合分析并做出了评价，得出2种工况所对应的合理铣挖顺序。王新明[9]提出铣挖与钻爆法联合施工时，上台阶采用铣挖法应遵循由下往上的原则进行开挖。目前的研究并未涉及隧道上台阶采用悬臂式掘进机整体开挖时，不同的铣挖步序对隧道围岩稳定性的影响。悬臂式掘进机铣挖隧道，并没有规定的铣挖掘进步序，理论上来说有无数种铣挖步序，但不同的铣挖步序不仅影响开挖效率，还会影响隧道开挖后围岩的稳定性。本文在现有研究基础上，结合国内隧道工程已采用的铣挖步序，利用有限元分析软件abaqus对不同方案的铣挖步序进行模拟，结合围岩稳定性的判定标准，并引入灰色关联分析法对模拟的结果综合评价，确定有利于隧道围岩稳定的合理铣挖步序，为后续类似工程提供借鉴。
在国内，悬臂式掘进机最开始是应用于煤矿开采，之后才逐渐应用于隧道的开挖。随着悬臂式掘进机功能的不断完善和相应辅助工艺的不断改进，截至目前已采用悬臂式掘进机顺利开挖多条隧道。悬臂式掘进机铣挖隧道施工工法，对围岩扰动小，减小了对施工周围环境的影响。但受限于现有悬臂式掘进机的设备条件，工程实践表明，悬臂式掘进机一般适用于围岩单轴抗压强度40~60 MPa的地层，在60 MPa以上的较硬岩层开挖则显得比较吃力[10]。对于围岩强度较低的隧道工程，为保证工程安全一般选择台阶法进行施工。本文参考国内已采用悬臂式掘进机铣挖隧道工程中的铣挖步序，结合悬臂式掘进机设备条件，总结了4种铣挖路径的方案，采用铣挖法开挖隧道上台阶部分，如图1所示(图中数字表示铣挖步序路径)。第1套方案为切割头从左往右切割，按照从下到上的顺序进行开挖；第2套方案为水平循环切割；第3套方案为切割头从下往上切割，按照从左到右的顺序进行开挖；第4套方案先在隧道中线从下往上竖直开挖中槽，再按照第1套方案的顺序完成开挖。4种方案中铣挖起始位置或铣挖方向的差异必然会引起开挖后隧道围岩受到不同程度的扰动，下文将介绍优选出其中合理步序的方法。
(a) 方案1；(b) 方案2；(c) 方案3；(d) 方案4
图1 铣挖步序示意图
Fig. 1 Excavation sequences schematic diagram
综上所述，本文旨在探讨和研究4种方案中不同的铣挖步序对围岩稳定性的影响，影响隧道围岩稳定性的因素非常复杂，既包括隧道的工程地质和水文地质条件，又与隧道开挖的工法、支护、隧道的断面形式尺寸、隧道埋深等密切相关。本文着重分析悬臂式掘进机铣挖步序的对隧道围岩稳定性影响，其他条件在创建数值模拟计算模型时取相同的参数，即视为一致。分析围岩稳定性的客观判定标准主要有位移值和围岩岩体单元安全度这2个方面。具体包括了拱顶沉降、周边收敛、拱顶单元体安全度、拱腰单元体安全度4个指标。其中拱顶沉降、周边收敛可以通过数值模拟的结果直接查询获得具体的数值。岩体单元安全度则需要根据数值模拟的结果做相应的计算获得。根据文献[11]，对于岩石介质，满足Drucker-Prager准则屈服条件的岩体破坏安全系数Fs可通过式(1)求得，
式中：I1为第一应力不变量；J2为应力偏量第二不变量，σ1，σ2和σ3均为主应力；α，K为与岩石内摩擦角φ和黏聚力c有关的试验常数，当与库伦六边形的外顶点重合时，它们与φ和c的关系为：
综上，总结了隧道围岩稳定性的判定标准以及对应参数的求取方法。
根据前文确定的4种铣挖方案以及与围岩稳定性相关的4个判定标准，要优选出最有利于围岩稳定的铣挖方案，必须要综合分析4种方案各自对应的4个判定标准的参数数值。上述问题属于多因素综合比选问题，灰色关联分析法是一种对各影响因素作用程度进行量化比较的有效工具，本文将采用灰色关联分析法来求解此问题，具体分析评价步骤如下：
将制定的4个铣挖方案视为4个评价对象，4个与围岩稳定性有关的判定标准视为4个评价指标(评价标准)。这样便形成原始数据矩阵Z，
矩阵Z中，Zi(j)表示第i个评价对象中对应的第j的评价指标的参数值，即每一行元素代表着一种方案对应的4个判定标准的参数值。
类冰川地质遗迹分布于黑龙江省中部的广阔区域，地理坐标为东经126°30′～130°30′，北纬45°40′～48°30′之间，目前包括黑龙江省铁力市朗乡花岗岩石林地质公园、宾县二龙·长寿山地质公园、集贤县七星峰地质公园和伊春市汤旺河小兴安岭石林国家地质公园、黑龙江凤凰山国家地质公园。其分布之广、种类之齐全、数量之多、特征之明显和保存之完好，在国内亦属少见(表1)。
第1步，确定参考数列。为了准确直观的进行比较，需要确定参考数列，参考数列是一个理想的比较标准，是以各评价指标的最优值构成的数列，可设Z0={Z0(1)，Z0(2)，Z0(3)，Z0(4)}表示参考数列。上述4个评价指标的最优值分为2种类型。一种为评价指标的参数数值越大则越优的效益型参数，此类型指标包括拱顶单元体安全度和拱腰单元体安全度；另一种为评价指标的参数数值越小则越优的成本型参数，此类型指标包括拱顶沉降和周边收敛。根据上述评价指标最优值的确定标准，便可从原始数据矩阵Z中获得参考数列。至此，将参考数列Z0写入矩阵Z的第一行，构成新的评价矩阵D，见式(7)。
第2步，无量纲化处理，建立新的矩阵E。由于不同的评价指标对应的参数的赋值量级存在差别，需要对矩阵元素进行无量纲化处理，无量纲化处理的方法包括很多种，本文采用的无量纲化处理的函数如式(8)和式(9)。式(8)为越大越优效益型指标无量纲化的处理函数，式(9)为越小越优成本型指标无量纲化的处理函数。
式中：Zi(j)表示指标的实际数值；Ci(j)表示指标无量纲化处理后的数值，Zmin(j)表示第j个指标中包含的所有实际值的最小值，Zmax(j)表示第j个指标中包含的所有实际值的最大值。通过上述方法对矩阵D中的每个元素进行处理，便可以得到矩阵E，见 式(10)。
第3步，关联度计算。将C0={C0(j) j=1,2,3,4}作为参考数列，Ci={Ci(j) j=1,2,3,4}作为比较数列。计算第i个铣挖方案的第j个评价指标Ci(j)与其对应的最优指标C0(j)的关联系数，用xi(j)表示，计算方法如下式(11)，
式中：ρ表示分辨系数，ρÎ[0,1]，一般取0.5[12]，其作用在于提高关联系数之间的差异显著性。环境参数Δmin和Δmax按式(12)和式(13)计算：
所有病例均使用双功能彩色多普勒超声诊断仪（飞利浦iU-22超声诊断仪）进行颈动脉超声检查，由接受过全国脑卒中筛查培训的高年资医师进行检测，以美国放射学年会超声会议公布的诊断标准，颈动脉狭窄程度分为四个等级，分别为正常、轻度狭窄（＜50%）、中度狭窄（50%～69%）、重度狭窄≥70%。
猕猴桃作为新兴水果，市场售价比普通水果更高，而且也是速生丰产的树种，常规管理可达一年栽植、两年见果、三年投产、四年丰产，故可获得较高的经济价值。
根据灰色关联度分析法模型，关联度记为γi，其表达式为：
式中：ω(j)表示各评价指标的权重，一般根据工程实际情况来取值，且所有评价指标的权重之和等于1。这里根据实际情况，拱顶沉降、周边收敛、拱顶岩体单元安全度、拱腰岩体单元安全度的权重分别取ω1=0.3，ω2=0.3，ω3=0.2，ω4=0.2。
太湖流域是一个不可分割的有机整体，上述水问题不是单独、孤立的，许多是经济社会发展过程中出现的新情况、新问题，这些问题相互关联，互为因果，各种水问题之间存在着不可分割的联系。鉴于太湖流域的水问题是流域性综合性水问题，《太湖流域管理条例》的内容涉及防洪，供水，水资源配置、保护和水污染防治等诸多方面，为综合性行政法规，重点将国家水资源管理与保护、水污染防治的法律制度在太湖流域具体化，对水环境实施综合管理与治理，实行比其他流域更加严格的水资源管理、保护和水污染防治制度。
第4步，评选最优方案。根据以上步骤可以计算比较序列和参考序列的相似程度，即得到4种方案与理想方案的接近程度，关联度的数值越大，则表明该方案相对其他方案更有优势，更加合理。
悬臂式掘进机适用于围岩强度较低的隧道工程的开挖，本文参照《铁路隧道设计规范》[13]中Ⅳ级围岩物理力学指标作为数值模拟创建的隧道模型的围岩条件，如表1所示。
表1 围岩物理力学指标
Table 1 Surrounding rock physical and mechanical indicators
围岩级别容重/(kN∙m−3)变形模量/GPa泊松比内摩擦角/(°)黏聚力/MPa Ⅳ201.30.33270.22
采用大型有限元分析软件abaqus对4种铣挖方案进行模拟。本文建立的隧道模型最大净空高度7.2 m，最大跨度5.6 m，采用上下台阶法开挖隧道，并利用悬臂式掘进机开挖上台阶部分。工程中常见的悬臂式掘进机截割头尺寸为976 mm×1 056 mm，截割的面积可近似按1 m2计算，每个循环开挖进尺取0.5 m，因此将模型中的开挖部分切割成1 m×1 m×0.5 m的小块区域。为了保证数值模拟尽可能的准确，数值模拟模型应取分析范围的3~5倍，隧道两侧边界各取3倍隧道最大跨度，上下边界也各取3倍隧道最大净空高度，隧道纵向长度取10 m，其中0~5 m部分上台阶已经开挖完成，5~10 m部分为待开挖部分，模型总尺寸为39.2 m×50.4 m×10 m，分析模型如图2所示，初始应力状态按重力场考虑。
图2 计算模型
Fig. 2 Calculation model
数值模拟结果获得的各评价指标的参数如表2所示。根据数值模拟结果获得的参数，分辨系数ρ取0.5；拱顶沉降、周边收敛、拱顶岩体单元安全度、拱腰岩体单元安全度的权重分别取ω1=0.3，ω2=0.3，ω3=0.2，ω4=0.2，按照前述方法计算求出每个方案的关联度，计算结果如表3所示。
表2 评价指标参数
Table 2 Evaluation index parameter
方案序号拱顶沉降/mm周边收敛/mm拱顶岩体单元安全度拱腰岩体单元安全度 11.930.841.571.11 21.840.781.621.38 31.780.991.311.15 41.670.851.531.44
表3 关联度计算结果
Table 3 Crelation degree calculation result
方案序号关联系数xi(j)(j=1~4)关联度γi 1234 10.333 30.636 40.756 10.333 30.508 8 20.433 31.000 01.000 00.733 40.776 7 30.541 70.333 30.333 30.362 60.401 7 41.000 00.600 00.632 71.000 00.806 5
根据计算结果，4种方案的对应的关联度数值分别为：γ1=0.508 8；γ2=0.776 7；γ3=0.401 7；γ4= 0.806 5。方案对应的关联度数值越大，则表明该方案与理想方案接近程度越大，γ4>γ2>γ1>γ3，其中方案4对应的关联度数值最大，所以方案4为4种铣挖方案中的最优方案，按照方案4的铣挖步序对围岩扰动最小，更有利于隧道围岩稳定。
1) 结合悬臂式掘进机的设备条件和国内已经的隧道工程铣挖经验，确定了4种铣挖步序方案。
偶尔也会有设定的东西能全部用上的例子，那就是《云阳陈寂传》。在写陈寂这个人物时，就设想过他的师门传承，还有神奇的泡茶功夫，经常和朋友聊到。后来给陈统领写番外，居然全部用到，也算是意外之喜吧XD。
2) 通过分析总结了与隧道围岩稳定性有关判定标准以及对应参数的确定方法。
农村人居环境改善直接关系农民群众的切身利益。倪虹介绍，在推进脱贫攻坚和乡村振兴过程中，住建部创造性地提出乡村“美好环境与幸福生活共同缔造”新模式。村里有哪些待解决的问题，由村民自己提出。房屋、厕所改造采用什么方案，由村民和专家商量决定。在共同缔造活动中，村民参与项目决策，通过投工投劳等方式不同程度地参与村庄建设，成果由村民代表参与评价。
3) 采用灰色关联分析方法对4种方案数值模拟的结果进行综合评价，得出先开挖中槽，再水平开挖的铣挖步序有利于围岩稳定。
联系客服
