两河口长隧道独头掘进压入式施工通风三维数值模拟

摘要：两河口公路隧道埋深大且无条件修建斜竖井，只能独头掘进3 000余 m，施工通风问题突出。基于CFD理论，对压入式通风气流运动采用三维紊态RNG kε湍流模型进行三维数值模拟，得到了施工期隧道内流场和浓度场随时间在洞内的分布变化规律。研究表明，掌子面附近为回流区，风流结构复杂，回流区以外风流分布逐渐稳定；爆破后的有害气团在隧道内是一个动态“移动”和“扩散”的过程，这个过程中有害气团从爆破掌子面逐渐扩散，并被稀释和排出隧道。同时，在3种工况下，根据CO的进入浓度和允许浓度对工作人员的进洞时间进行了研究。

关键词：长隧道；压入式通风；三维数值模拟；进洞时间

中图分类号：O319.56

文献标志码：A

文章编号：16744764（2014）02003507

Abstract：There is no suitable condition to build any shaft（vertical or oblique shafts） due to the complex geological environment of Lianghekou. Furthermore， the tunnel has to be constructed blind heading. So the problem of ventilation is serious. Based on the theory of CFD，the wind current produced in forced ventilation is simulated by 3D RNG kε Turbulent Model，and then the Timedependent Variation Rules of the flow field and the concentration field during the construction period are shown by numerical simulation of 3D. Results show that the backflow area which is located near the working face has the complex situation of wind current， and the air current distribution becomes stable gradually when it is far away from the backflow area. Meanwhile， the movement of the harmful air produced after blasting can be described as “moving” and “diffusing”. The result of moving makes the harmful air run out of the tunnel from working face. The result of diffusing makes the harmful air deliquated along with the process of moving. In addition， the time that the workers spend on getting into the tunnel has been studied according to the concentrations of CO （in and allowance concentrations respectively）.

Key words：long tunnel；forced ventilation；the 3D numerical simulation；entrance time

随着西部大开发进程的加快，出现了很多特长公路隧道。而目前特长隧道施工通风往往仅根据规范和施工人员的经验进行施工通风[1]，效果很差，一旦措施不合理，对施工人员的身体危害很大。在特长公路隧道施工中，由于独头掘进和隧道长度的关系，如果等炮烟和其他有害气体完全从隧道内排放出去，则需要非常长的时间，会大大影响工程进度；如果选择较短时间进入施工，则对施工人员危害很大。因此，解决特长隧道施工通风问题具有重要意义。

在公路隧道通风问题的研究中，其他国家开展相关研究较早，特别是瑞士、挪威、日本以及奥地利等国家，对公路隧道通风问题进行了大量研究[24]。但由于全真隧道通风的试验费用非常高，所以在这方面取得的成果有限，大部分研究主要集中在计算机数值模拟领域。而且，研究方向主要在隧道建成后的运营通风[58]，对隧道施工阶段的通风问题却少有研究。在中国，相对于研究公路隧道运营通风，隧道施工通风研究非常少，而且研究成果主要集中在公路隧道纵向通风的影响因素、通风方式、有害气体浓度分布规律以及纵向通风计算模型等方面[58]。对仿真试验研究、隧道三维流场数值模拟等方面几乎很少涉及[912]。

目前，随着计算机技术的快速发展，使隧道施工通风的三维数值分析、流场和浓度场动态变化研究成为可能。本文在总结和借鉴已有的公路隧道通风运营、施工通风的研究成果上[1319]，采用大型有限元软件ADINA中的CFD模块对独头掘进的长隧道进行整条隧道有限元数值分析。

1两河口隧道工程概况

两河口水电站交通工程1#公路隧道位于四川省甘孜州雅江县境内，1#公路隧道全长5 855 m，隧道海拔高程3 000 m左右。隧道埋深大，不适宜修建竖井或斜井，做到长隧短打，只能独头掘进3 000余m，通风问题成为制约隧道快速掘进的瓶颈。

两河口1#公路二级公路，隧道内轮廓为三圆心设计，隧道内轮廓净宽度为11.0 m，高7.3 m。隧道施工方法为全断面开挖方法。隧道采用压入式通风，通风管距地面3 m左右，风管直径为1.8 m。爆破产生的主要气体有：一氧化碳、二氧化碳、氮氧化物、二氧化硫和三硝基甲苯（TNT）等[20]。由于有害气体中CO危害较大，在爆破前采用水幕降尘措施，使得爆破后气体中主要剩下CO且浓度稳定，因此，隧道施工通风中以CO浓度大小作为评判通风效果的依据。

2 隧道通风数值计算原理

采用有限元软件ADINA中的CFD模块对不同工况下整条隧道的通风进行数值模拟。隧道压入式通风为贴壁受限射流，流场中存在射流、回流、滞流等气体流动，因此，采用三维紊态RNG kε湍流模型进行分析。

2.1基本假定

基本假定主要有以下3个方面：1）假定通风气流为三维粘性不可压缩流体；2）气体流动不产生热能耗散，忽略隧道壁的传热性；3）从已有成果来看，爆破后产生的有害气体在稀释和向洞外移动过程，主要取决于初始有害气体量，而与初始分布关系不大，因此，假定通风前初始有害气体均匀分散在出风口和掌子面之间。

2.2数学模型

根据假定，隧道风流模型适用于Navierstokes方程，并采用RNG kε 紊流模型使方程组封闭。隧道内空气湍流流动和浓度扩散控制方程如下。

4数值计算结果及分析讨论

4.1隧道流场分析

根据总结现有研究成果，隧道独头掘进中，风流变化规律在不同掘进距离规律大致相同[19]。因此，模型对隧道独头掘进时，掘进距离为3 130 m时进行重点分析。

在分析中，隧道内流场很快稳定，而且从洞口到掌子面流场稳定所需的时间逐渐变短。图2为洞口断面平均风速随时间的变化规律。洞口断面风速在100 s后就已基本稳定，洞内其余部分的流场稳定时间均小于100 s，在后面的分析中，可不考虑流场随时间的变化，把浓度场随时间的变化规律放在稳定的流场中分析。

根据其他研究成果和本模型计算结果，把隧道流场分为3个区域：回流区、回流影响区、稳定区。

1）回流区回流区位于掌子面到风筒出口附近。根据计算结果，截取掌子面附近Z=3.8 m截面的流场矢量，如图3所示。从图中可见，从掌子面到风筒出口附近的风流形成一个漩涡。风流在射向掌子面后，折射后从隧道另一侧返回，当其返回到风管出口附近时，部分风流受风管出口风流影响又流回掌子面，部分风流流出此回流区域。根据计算，回流区位于风管出口与掌子面之间，距离掌子面约18 m。图4为回流区中心（X=3 112 m）横断面风速分布图，图中G是风速为0的等风速线，此线基本与隧道断面中线重合，其左侧风速为流向掌子面方向，其右侧风速为流向洞口方向。结合图3、4，可知掌子面附近流场是一个以距掌子面18 m处横断面的中轴线为中心，以风管出口附近到掌子面为范围的大漩涡。

2）回流影响区根据分析结果，回流影响区为风筒出口附近到距离掌子面约300 m左右的一个区域。这个区域受回流区流场的影响，风速因断面的位置不同而差异较大，并且规律性不明显。图5为隧道掘进至3 130 m工况下不同断面和同一断面不同高度的风速分布图。从图可见，距离掌子面越近，风速波动越大。当距离掌子面超过300 m，风速逐渐趋于稳定。

3）稳定区从图5可见，距离掌子面300 m，风速区域稳定。同时，风速在断面的分布是中间大，逐渐向隧道侧壁减小。这一规律，与流体力学理论完全一致。图6为X=1 000 m处横断面风速分布等值线图，从距离掌子面约300 m以后各横断面风速分布均基本与此图相同。

4.2CO浓度随时间变化规律

根据前面对流场的研究，可以看出掌子面附近的回流区会阻碍有害气体的排出。同时，风速中间大，周围小的特点对有害气体的排出也有一定的影响。下面仅对3 130 m情况下的规律进行研究。

隧道内浓度场变化规律与流场分布规律密切相关，图7、8分别为隧道掘进3 130 m时通风1 min和5 min时回流区（Z=3.8 m）CO浓度场分布。从图中可见，在距离掌子面18 m左右形成一个大漩涡，其中心CO浓度最高。并且中心处与流场的回流区中心重合。随着时间的推移回流区CO浓度逐渐降低，梯度基本不变，但其回流区中心浓度大并向周围逐渐减小的规律保持不变。

图10～12以及图9分别为隧道掘进500、1 000、2 000、3 130 m时洞内CO浓度在不同时间的分布。图中进入浓度取为100 mg/m3，安全浓度取为30 mg/m3。这4副图的规律基本相同，与5.2中分析的掘进3 130 m时的基本规律相符，仅因条件不同具体数值不同。这4副图通风5和9 min的曲线都有双峰值，且第一个峰值影响范围均约离掌子面40 m内，这是由于4个模型均把风筒出口设在离掌子面35 m处，但由于前1 500和1 500 m以后掌子面供风量不同，掘进2 000 m和掘进3 130 m情况下，CO浓度第1个峰值的下降速度和第2个峰值的移动速度均明显大于掘进500和 1 000 m情况下的CO浓度相应峰值的变化速度。图10中可看出回流效应的有效影响仅为十几分钟，而该条件下通风时间需二十多分钟，掘进越深其相对影响时间就越短。所以，在较长隧道的通风散烟的数值模拟中对进洞时间的计算可不考虑回流效应对有害气体滞留作用的影响。

从表1中数据可以看出，通风距离越长CO峰值移动速度就越接近于风速。根据该模型的计算原理，当不存在风流的紊动CO的移动速度就是风速。掌子面附近的回流效应对CO的扩散有一定的影响，通风距离越长其相对影响就越小。所以，当通风距离较长时可近似认为CO峰值移动速度等于风速；当通风距离较短时需考虑掌子面附近的回流紊动区对通风时间的影响。

另外，图9～12中曲线以及表中数据反映出，图12基本可与图10前2 000 m重合，风速也相等。图11与图10也有相同的关系。因此，在精度要求不是很高的情况下，工程人员可在3 130 m条件下的数值模拟图中，分别在图9中截取不同掘进距离的CO浓度，并进行推断。

4.4进洞时间及安全时间建议

进洞时间即为洞内CO浓度降至进入浓度以下所需的时间，安全时间为洞内CO浓度完全降至允许浓度之下所需的时间，回流区影响时间为回流区CO浓度峰值减小到30 mg/m3以下所用时间。

根据前文所述分析方法、边界条件以及初始条件等，分别分析不同通风能力，在不同掘进距离条件下工作人员的进洞时间。具体的工况如下：1）工况1，掘进前1 500 m，供风量为3 000 m3/min，掘进1 500 m以后，供风量为4 500 m3/min；2）工况2，掘进前1 500 m，供风量为2 400 m3/min，掘进1 500 m以后，供风量为3 600 m3/min；3）工况3，掘进前1 500 m，供风量为3 600 m3/min，掘进1 500 m以后，供风量为5 400 m3/min。计算结果如表2所示。

在表2中，给出了不同通风能力，分别掘进500、1 000、2 000、3 130 m各情况下进洞时间、安全时间和回流区影响时间。施工人员可在计算进洞时间后进入工作面作业，此时，已过回流效应的影响时间，对施工人员影响较小。同时，从表2中的时间差一项可以看出当隧道掘进距离越长时，掌握进洞时间对隧道快速施工的意义就更为明显。

在建模条件下，回流区影响时间只与掌子面供风量有关，掌子面供风量越大其影响时间就越短，回流区影响时间在较长的通风情况下与进洞时间相比相对越短，所以在较长隧道的通风中可不考虑此效应。

5结论

根据两河口1#公路隧道的实际情况和施工方法，建立了隧道掘进500 、1 000、2 000和3 130 m等4种情况下隧道施工通风的模型，并对隧道流场分布和CO浓度变化规律进行了研究。通过研究掘进3 130 m工况下隧道内流场及CO浓度场的变化，得出独头掘进隧道施工通风系统流场分布和浓度场变化的一般规律。对4种情况下CO浓度场随时间变化规律进行对比，并对不同工况下施工人员合理的进洞时间进行分析。主要结论如下：

1）隧道独头掘进时，采用压入式进行施工通风时，在距离掌子面18 m左右形成回流紊动区，范围为距离掌子面40 m之内；在距离掌子面40～300 m范围为回流影响区；超过300 m，为流场稳定区。同时，有害气体团在逐渐向洞口移动的过程中发生扩散，中心浓度逐渐减小。

2）掌子面附近回流区对CO气体排离掌子面有一定的阻滞作用，回流区中心形成一个CO浓度的局部峰值，CO浓度也以一个与回流区流场相似的漩涡状从中心向周围逐渐降低。在相同的隧道轮廓、相同的风管直径以及相同的风管出口到掌子面距离的情况下，回流区对CO的阻滞时间仅与掌子面供风量有关，供风量越大，阻滞时间越短。在较长隧道的施工通风中回流区影响时间相对较短，可不考虑此效应对进洞时间的影响。

3）CO气体是隧道流场介质的组成部分，其运动速度应等于流场介质的运动速度，但由于回流区对CO移动的阻滞作用，在通风距离较短时CO平均移动速度明显小于风速，当通风距离较长时CO气体移动速度可近似等于风速。

4）根据规范和具体施工条件，分析了3种工况下，独头隧道掘进500、1 000、2 000、3 130 m施工人员的进洞时间和安全时间。

参考文献：

[1]张恒，杨家松，高辉. 锦屏隧道施工通风研究[J]. 广西水利发电，2008（5）：1113.

[2]Rudin C. Fires in long ruilwuv lunnelsthe ventilation concepts adopted in the alptransit projects [C]//10th International Symposium on the Aerodynamics Ventilation of Vehicle Tunnels， 2000.

[3]Chow W K.Dispersion of carbon monoxide from a vehicular tunnel with the exit located along a hillside [J]. Tunneling and Underground Space Technology， 1989， 2（4）： 231234.

[4]Vardy A. On Semitransverse ventilation systems [C]//7th International Symposium on the Aerodynamics and Ventilation of Vehicle Tunnels， 1991： 629646.

[5]杜翠凤，于明章，李怀宇.单条沿脉巷道通风网路炮烟运动分析及对策[J]. 金属矿山，2000，6（6）：4042.

[6]胡自林，彭立敏，刘胜利.中长公路隧道纵向通风计算模型[J].湘潭矿业学院学报，2003，18（2）：3437.

[7]王永东，夏永旭.公路隧道纵向通风系统局部影响数据模拟[J].西安公路交通大学学报，2001，21（4）：5054.

[8]危宁，李力，王春燕. 隧道施工通风中的有害气体浓度变化分析[J]. 三峡大学学报：自然科学版，2006，28（4）：324327.

[9]高辉，杨家松，陈寿根，等.锦屏辅助交通特长隧道爆破通风数值模拟[J].公路遂道，2008（1）：1013.

[10]刘钊春. 独头掘进隧道施工通风数值模拟[D]. 西安：西安理工大学，2010.

[11]刘钊春，柴军瑞，贾晓梅，等. 压入式通风掘进面有害气体浓度扩散数值模拟[J].岩土力学，2009，12（30）：536539.

[12]陈卫忠，郭小红，曹传林，等. 公路分岔隧道循环风相互影响及其对策研究[J].岩石力学与工程学报，2008，27（6）：11371147.

[13]赵海东.乌鞘岭隧道辅助坑道对运营通风的影响及特长隧道定点防灾研究[D].成都：西南交通大学. 2009.

[14]Patankar S V， Spalding D B. A calculation procedure for heat， mass and momentum transfer in threedimensional parabolic flows [J]. International Journal of Heat and Mass Transfer， 1972， 15（15）： 17871806.

[15]王洪德，林琳，赵轶. 地铁隧道火灾事故通风方式数值模拟[J].辽宁工程技术大学学报：自然科学版，2010，29（2）：177181.

[16]李炎锋，朱滨，孙旋，等. 利用大涡模拟研究地铁区间火灾的烟气扩散[J]. 北京工业大学学报，2007，33（10）：10601065.

[17]Woodburn P J， Britter R E. CFD simulations of a tunnel firePart Ⅰ [J]. Fire Safety Journal， 1996， 26： 3562.

[18]马洪亮. 基于区网藕合模拟的矿井通风系统抗灾能力分析[D]. 北京：中国矿业大学，2008.

[19]Kunsch J P. Simple model for control of fire gases in a ventilated tunnel [J]. Fire Safety Journal， 2002， 37（1）： 6781.

[20]张国枢. 通风安全学[M]. 徐州：中国矿业大学出版社，2000.

[21]赵彬，李先庭，彦启森.入口紊乱参数对室内空气分布的影响研究[J].建筑热能通风空调，2000（1）：14.

[22]高建良，张生华. 压入式局部通风工作面风流分布数值模拟研究[J].中国安全科学学报，2004，14（1）： 9396.

[23]Nakayama S，Uchino K，Masahiro M. 3 Dimensional flow measurement at heading face and application of CFD [J].The Mining and Materials Processing Institute of Japan， 1996， 112（9）： 638644.

[24]TB 10204—2002 铁路隧道施工规范[S].北京：中国铁道出版社，2002.

（编辑王秀玲）

推荐访问： 河口 掘进 数值 通风 隧道

---