公路隧道空气幕与射流风机通风方法对比及优化

2020-07-07 04:13张博文王海宁张迎宾
有色金属科学与工程 2020年3期
关键词:射流风量风速

张博文, 王海宁, 张迎宾

(1. 江西理工大学,.a 资源与环境工程学院; b.能源与机械工程学院,江西 赣州341000;2. 中国计量大学质量与安全工程学院,杭州310018)

射流风机作为特殊的一种轴流风机[1],因为其布置灵活、资金成本较少等优点,在隧道实际通风项目中,一直受到专家学者的青睐。 近几十年来无论在理论还是模拟方面,国内外学者对射流风机都有着大量的研究。 在理论方面,赫海仙[2]通过平板和风机的实验,说明了风机的安装距离不同和风机距离平板的间距能够很大程度影响平板受到阻力和风机的推力比。崔德振等[3]采用了实地测量和模型实验的研究方法,对比不同汽车模型在压力差作用下形成风量和风速的变化情况, 同时验证了自模区的存在临界雷诺数,为纵向射流风机模型的设置提供了方法。 李晓菲等[4]通过理论分析,发现影响纵向最小间距的影响因素主要分为4 个方面,分别是叶轮直径,射流出口风速,风机轴线之间的横向间距以及风机与拱顶的接近度。邢利英等[5]验证了在射流风机出口处增设收缩管的方法,从而提高隧道射流通风效率,能够很有效地改善隧道通风效果,作为通风优化的基本研究方法,对实际问题进行现场测试以及数值模拟,提出改进方案并得到优化效果的通风优化处理模式。 在数值模拟方面。 周家才等[6]编制了一套纵向通风的程序模拟气流,对车辆阻力、污染物的浓度分布,进行了有效的分析,并且在实际中进行了运用。 韩国的Joon-Hyung 等[7]在三维PANS 和PBA 代理模型的基础上提出了能够优化隧道射流风机性能的程序。 国内WANG F 等[8]模拟探究弯曲公路隧道中射流风机的各个横截面位置对气流流场的影响,他们发现当风机和隧道顶棚的距离变化越来越大的时候射流速度越来越快,当射流风机和隧道壁的距离越来越近或者风机之间的距离越来越近时,风机需要更长的距离来运行。 吴珂等[9]模拟指出合理的洞口距可以避免进口气流和风机出口气流相互的影响,并且对射流通风的效果有良好的促进作用。

虽然射流风机在现今的隧道通风中研究较成熟,但随着研究的不断深入慢慢地也开始出现一些问题。例如, ①射流风机在运行过程中存在能量利用率低、耗电量大、通风效率不高,尤其在大浓度的区域处理尾气不理想, 很多时候达不到所需要的风流调控[10]。②射流风机安装的位置主要是在隧道的顶部,虽然很大程度增加了车辆的行驶空间,但是,在清理维护方面会存在车辆伤害, 对操作人员有着很大的危险性。③在传统射流风机安装过程时,风机安装高度对风机的通风效率有着很大影响,同时风机在运行过程中长时间承受风机动荷载,对风机支架和减震装置的安装也有一定技术的要求等。

据了解,一些井下大断面且长度大的巷道因为需要采掘和运输,大型器械需要不断进出,同时还伴随着炮烟的冲击影响, 巷道内的空气质量状况与长隧道相比有过之而无不及, 而硐室型风流调控技术——空气幕能有效地解决此类问题[11-13]。硐室型空气幕顾名思义安装在侧壁的硐室内的空气幕, 由于安装位置的特殊,不仅不影响车辆运输,而且还可以实现对风流的增阻、引射、隔断的作用[14-15]。 现今,硐室型空气幕在矿井的研究和应用较多, 在隧道通风中应用尚未见报告, 但空气幕的特性决定其具有很广泛的运用领域,因此,论文以瓮福磷矿隧道通风为研究背景,运用FLUENT[16]进行模拟和分析,对比模拟过程中空气幕和射流风机在隧道中风量和风速变化状态,找到较优组合,为今后空气幕引射风流在隧道中的可行性增加了筹码。

1 实际隧道的通风效果分析

1.1 瓮福磷矿穿岩洞隧道概况

瓮福磷矿穿岩洞隧道是一座双向单洞隧道,隧道入口处于蚂蚁坟地段, 坐标为X:87668.514 m,Y:39795.449 m,Z:1364.24 m;出口处为原厂地段,坐标为X:88394.867 m,Y:39232.036 m,Z:1327.48 m。 隧道全长近919.8 m,纵坡4.0%,隧道路面横坡定为2%。该隧道区域属中山地区,山间沟谷地貌,总体地势东高西低。 隧道所经地段主要为中山山脊地带,海拔高程在1100~1500 m 间,相对高差近400 m。隧道出、入口段均为山脊斜坡地段,地表自然坡度在25~30 ℃左右,同时隧道围岩工程地质条件较差。

该隧道按双向二车道二级公路单洞隧道设计,隧道建筑限界按40 km/h 行车速度确定, 隧道内轮廓采用三心圆形式,单洞建筑限界净宽9.5 m、净高5.0 m,行车道(含路缘带)宽8.0 m,检修道宽0.75 m。 单洞净空面积正常路拱段为56.16 m2, 带仰拱地段为69.56 m2。 同时,在左右两侧检修道及余宽下设置一个尺寸为50 cm×50 cm 的电缆沟, 路面两侧各设置直径为20 cm 的圆形开口型排水边沟。

关于隧道通风方面,考虑到解决隧道内行车尾气排出的通风问题,瓮福穿岩洞隧道在原有的通风设计中,安装了两道射流风机进行纵向通风,即在距离隧道进、 出口近100 m 处分别在隧道顶部安装了2 台射流风机并联通风,其中1 台运行,1 台备用。

1.2 瓮福磷矿穿岩洞隧道存在的通风问题

考虑到隧道封闭性、狭长性、四面不可及性等特点, 通过分析瓮福磷矿穿岩隧道的通风和运营状况,得到主要的通风问题如下:

1)瓮福磷矿穿岩洞隧道是一座双向单洞隧道,其内部的风流大小和方向会随车流方向、车速、交通量的变化而变化,并且隧道内空气流动的大小和方向还受隧道内、外气候条件的影响很大,这导致了隧道内的自然风向表现出极大的不确定性,自然风不能有效排出行车产生的污浊空气。

2) 穿岩洞隧道内原安装的两道射流风机相距较远,如图1 所示,射流风机接力效果减弱,加上自然风压在不同自然条件下对其产生的积极或者消极的影响,使得射流风机引射风流的作用大大削弱。 当射流风机运行时,隧道进、出口附近的通风效果较好,但在中部易出现污风停滞, 污浊空气难以及时排出隧道,导致汽车运行的可见度降低,增加行车危险度;当射流风机停止运行时,隧道内行车产生的污浊空气主要靠自然通风排出,由于自然风压的不稳定,因而隧道内的污浊空气紊乱,且难以及时排出。

图1 原射流风机现有安装示意Fig. 1 Existing installation diagram of original jet fan

3)在两处各安装2 台射流风机, 其中1 台运行,1 台备用,不仅射流风机数量偏少,而且排风的接替能力偏弱。

4)由于隧道全长约920 m,隧道内双向行车速度慢,交通风相互影响,易导致风流紊乱,汽车排放尾气难以排出,尤其是在通车量较大的情况下,隧道内空气质量差,且行车能见度较低,如图2 所示,不仅影响了行车安全和行车效率,而且严重危害到汽车司机和行人的身体健康。

2 实际隧道的需风量计算

2.1 公路隧道通风卫生标准

根据JTGTD70.02.2014《公路隧道通风照明设计规范》相关部分的说明,计算CO、烟雾和异味进行稀释所需的3 种需风量中的最大值为隧道通风的最终需风量标准。 并且隧道通风标准要求如下:

图2 隧道通风现状Fig. 2 Tunnel ventilation status map

1) 双向交通的隧道设计风速不应大于8 m/s;单向隧道设计风机不应大于10 m/s, 特殊情况可取12 m/s;人车混合隧道风速小于7 m/s。

2)人车混合隧道长度小于1000 m,同时CO 设计浓度为150×10-6。

3)采用钠灯光源时,行车速度在40 km/h 时,烟雾设计浓度:0.009 m-1,浓度大于0.012 m-1时,采取交通管制。

4)在隧道进行纵向通风换气时,隧道内换气风速不应低于2.5 m/s。

2.2 瓮福磷矿穿岩洞隧道需风量计算

依据设计资料可知,穿岩洞隧道的设计行车速度为40 km/h, 则隧道全长稀释CO 的需风量可根据式(1)和式(2)进行计算,而烟雾和异味进行稀释所需的需风量可根据式(3)、式(4)和式(5)进行计算:

式(1)中:Qco为隧道全长CO 排放量,m3/s;qco为CO基准排放量,m3/(辆·km);fa为考虑CO 的车况系数;fd为车密度系数;fh为考虑CO 的海拔高度系数;fm为考虑CO 的车型系数;fiv为考虑CO 的纵坡-车速系数;n为车型类别数;L 为隧道长度,m;Nm为相应车型的设计交通量,辆/h。

式(2)中:Qreq(CO)为隧道全长稀释CO 需风量,m3/s;δ为CO 设计浓度,10-6;P0为标准大气压,kN/m2;P 为隧道设计气压,kN/m2;T0为标准气温,K;T 为隧道夏季的设计温度,K。

式(3)中:QVI为隧道全长烟雾排放量,m3/s;qVI为烟雾基准排放量,m2/(辆·km);fa(VI)为考虑CO 的车况系数;fh(VI)为考虑CO 的海拔高度系数;fm(VI)为考虑CO 的车型系数;fiv(VI)为考虑CO 的纵坡-车速系数;n 为车型类别数。

式(4)中:Qerq(VI)为隧道全长稀释烟雾的需风量,m3/s;K 为烟雾设计浓度,m-1。

式(5)中:Qreq(火)为隧道全长稀释异味的需风量,m3/s;vr为隧道中的平均风速,m/s;Ar为隧道净空断面积,m2。

根据上述隧道通风卫生标准和隧道换气最低风速要求,得到瓮福磷矿穿岩洞隧道的较优平均风速的范围为(2.5~8 m/s),代入式子(5)中得到全长异味所需要的需风量范围为140.4~449.28 m3/s,对比CO、烟雾和异味3 种需风量, 得到当隧道需风量在140.4 ~449.28 m3/s 范围内时, 瓮福磷矿穿岩洞隧道的通风要求达到隧道通风卫生标准。

2.3 瓮福磷矿穿岩洞隧道初始风速确定

为了确定在非风机状态下隧道中的风速, 运用ZRQF-D 系列智能式风速仪[17],采用侧身测量法来对隧道内的风速进行测定,测点主要布置在隧道两侧的进风口100 m 处、出风口100 m 处及隧道中部,共计6 个测点,测点布置示意如图3 所示,测点1、测点2、测点3 的行车方向是测点1→测点3,测点4、测点5、测点6 行车方向是测点4→测点6。 测量结果见表1。

对表1 中瓮福磷矿穿岩洞隧道现场调查的各项数据进行归类和计算,得到在无风机的状态下,隧道内自然风向多为东北风, 风速多在0.3~0.6 m/s 之间变化, 考虑到更好地进行2 种风机组合的模拟对比,将隧道在无风机状态下风速确定为0.5 m/s。

图3 测点布置示意Fig. 3 Schematic diagram of measuring point arrangement

表1 各测点风速、风量测定统计表Table 1 Statistics table of wind speed and wind volume at each measurement point

通过上述瓮福磷矿穿岩洞隧道通风问题的分析,了解到原有隧道通风系统在隧道中的产风量严重不足, 达不到隧道通风设计卫生标准, 为了更好地改善隧道现有的通风状况,同时探究空气幕替代射流风机加强隧道通风的可行性,根据计算得到的瓮福磷矿穿岩洞隧道通风卫生标准需风量的范围(140.4~449.28 m3/s)和平均速度的范围(2.5~8 m/s),提出了2 种新的风机组合优化方式。一种是在原有的射流风机中添加风机台数,即风机组合达到3 组两两并联总数达到6 台,另一种是以空气幕的引射功能为主的3 组两两并联空气幕的组合方式。

3 计算模型建立及边界条件设定

3.1 模型建立及网格划分

根据瓮福磷矿穿岩洞隧道实际长度, 确定模型隧道总长度为920 m,高度为H6=5 m,宽度为W6=9.5 m,隧道断面的面积为69.56 m2。 空气幕模型硐室高H3=H4=2 m, 长度为8 m, 宽2 m, 空气幕风机直径D=2R3=1.3 m,空气幕长为2 m,空气幕出风口与隧道成40°。射流风机模型D6=D7=1.1 m,风机长度L=2 m,风机间的横向间距L10=3.6 m。 应用GEOMETRY 建立的引射风流空气幕和射流风机几何模型如图4 所示。

运用ICEM 对模型进行网格划分[18],为了确保精准,在网格划分期间风机附近网格进行加密,隧道出入口网格设置较为疏松,2 种风机模型都采用结构网格,且质量都达到模拟要求,空气幕模型网格总数量为892264,射流风机模型网格总数量为689654。如图5 所示。

图4 射流风机和空气幕在隧道中的几何模型Fig. 4 Geometric model of jet fan and air curtain in tunnel

3.2 边界条件设定

对隧道模型进行FLUENT 模拟计算时, 主要是以RNG k-ε 的湍流模型为主, 离散项都采用二阶迎风格式,压力速度耦合选用标准SIMPLE。

具体边界条件设置如下:隧道进出口选用压力边界条件,根据实际通风效果评测结果表1,得到隧道在无风机状态下的风速大概为0.5 m/s (空气幕模型也一样)。因此,通过换算得到隧道入口边界条件设置压力为9.5 Pa,出口与外界大气相通设置为0 Pa。 其次, 风机内部压力面, 设置为边界条件FAN。 风机WALL 采用默认壁面光滑, 隧道壁面粗糙度设置为0.0025[19]。

射流风机压力面参数设置,射流风机是一种通过推力和流量来判断通风性能的特殊轴流风机,为了更好地进行对比,需要进行全压换算,根据轴流风机的全压计算公式及:

图5 射流风机和空气幕在隧道中的计算模型网格Fig. 5 Grid diagram of calculation model of jet fan and air curtain in tunnel

式(6)中:P 为全压,Pa;Pd为动压,Pa;Pst为静压,Pa。

获得射流风机30 kW 的动压Pd=0.5dv2=603 Pa,d=1.2 表示空气的密度,kg/m3。关于静压Pst,射流风机的静压主要是用来克服进口处的阻力和消声器的阻力[20],而射流风机静压计算可以通过消声器产品的规格和轮毂比来确定,但射流风机的静压很小,很多时候静压和大气接触后就消失,相比较于以动压为主的射流风机,静压很多时候影响意义不大。因此,为了模拟计算过程的方便,简化风机设计模型的复杂性,文中以风机静压无限小为前提条件, 设计射流风机模型前后无管道, 风机壁面光滑, 使得此时射流风机P=Pd, 射流风机达到最大的推力并且性能达到较优状态。

所以在模拟中可以理解为

式(7)中: P 为全压,Pa;d=1.2 为空气的密度,kg/m3;A为射流风机的出口面积,m2。

空气幕压力面参数设置, 空气幕中的风机为传统的轴流风机,由于模拟需要,根据全压特性曲线,对4 种功率的风机进行拟合,得到全压和风速方程式如表2 所列, 而各功率压力面所对应的压力参数就是全压风速拟合方程[21]。

表2 各功率的风压拟合方程Table 2 Wind pressure fitting equation for each power

4 数值模拟计算结果与分析

4.1 总功率相同下射流风机组合和空气幕组合模拟研究

2 种风机组合的单台风机功率选用30 kW,总功率维持在K空=K射=180 kW 相同情况下进行对比,如图6 表示单个风机功率为30 kW 的2 种风机组合在各自轴线高度上的隧道中线处沿隧道长度方向的速度变化曲线。 图7 表示2 种风机组在隧道中产生的平均风量。

图6 风机组在轴线高度上的速度变化曲线Fig. 6 The speed change curve of the wind turbine at the height of the axis

图7 同功率风机组合隧道中的平均风量Fig. 7 Average air volume of fan combinations of the same power in the tunnel

由图6 可知,2 种风机组合在距离入口50 m 处,因受风机组合的影响, 速度开始不断增加, 当到达150 m 时速度达到平稳且最大, 在入口200 m 处,由于距风机口的距离变远, 风速受风机影响开始减弱,风速也随着距离的增加开始减小,一直减小到距入口300 m 的位置,此时风速保持在平均风速水平线上振动,并且保持一定距离的平稳,直到距入口350 m 风速开始进入新一轮的相同的变化,类似的变化一直持续到末端风机组合的影响结束。 此外,通过2 种风机组合的对比,还发现空气幕风机组合的最高风速和最低风速差约为0.6 m/s, 而射流风机组最大风速和最小风速差约为0.7 m/s。 这说明2 种风机组合速度变化状况大致相同,都是先增加后减小,从最大与最小风速差状况下看,空气幕组合风速变化较为平稳,且速度曲线整体都要高于6.18 m/s,而射流风机组合风速较为逊色,整体速度曲线高于5.95 m/s。其次,从图7 可看出,无风机状态下隧道平均风量为19.71 m3/s,空气幕组合产生的隧道平均风量为244.63 m3/s, 射流风机组合产生的平均风量为234.94 m3/s。 综上图6、图7 得到,2 种风机组合对隧道空气质量都有着显著的提高,并且在同功率状态下,空气幕组合产生的风量要高于射流风机组合,因此如果把射流风机组合的模拟结果当成是一个标准,空气幕组合可以在相对于较低的总功率下完成这个通风要求,从而达到节能减排的目的。

4.2 不同总功率下的风机组合模拟研究

为了验证上述的假设,选用7 种最为接近180 kW射流风机组合通风效果,且低于180 kW 总功率的空气幕组合与180 kW 射流风机组合进行对比,如表3所列,K1,K2,…,K6表示风机型号,图8 表示空气幕各风机功率编号平面图。

表3 7 种风机组合的功率分配Table 3 Power distribution of 7 fan combinations

图9 表示8 种不同功率风机组在各自轴线高度上的隧道中线处沿隧道长度方向的速度变化曲线,图10 表示各个功率在隧道中产生的平均风量。

图8 空气幕各风机功率编号平面示意Fig. 8 Air curtain power number plan of each fan

图9 不同功率风机组在轴线高度上的速度变化曲线Fig. 9 Speed variation curve of different power fan combinations on axis height

图10 各个功率在隧道中产生的平均风量Fig. 10 Average air volume generated by each power in the tunnel

由图9 可知, 其速度变化过程和图6 所述一样,观察图9 得到速度曲线中最为接近180 kW 射流风机组合的是P3和P6空气幕风机功率组合,且P3相对于P6隧道风速比较稳定, 变化的频率也与射流风机组合重合度较高。 而P5和P7的速度曲线由于达不到射流风机组合的要求,因此不满足。 至于其他功率组合,虽然速度曲线都要大于射流风机组合,但考虑最佳的节能原则, 其他几种功率组合都要逊色于P3和P6。 除此之外, 还发现2 种较为特殊功率组合P3和P5, 即P5的风机组合总功率明显要高于P3的风机组合,但P3速度曲线明显要高于P5,这说明合理的功率分配能够在减小总功率的条件下,有效地提高通风效率。 从图10 可以看出,P3功率组合产生风量为234.93 m3/s 能够节约34.5 kW 功率,P6功率组合产生风量为236.55 m3/s 能够节约26.5 kW 功率, 虽然2 种功率组合产生风量近似相等, 但P3总功率明显要小于P6,因此P3功率组合的节能减排优于P6功率组合。 综上所述,除了P7和P5两种风机功率组合,其他风机组合速度曲线都能够高于或者近似等于总功率P=180 kW 射流风机组合的风量和风速要求,而在这些风机组合中P3=145.5 kW 功率组合是满足要求前提下,总功率最低的风机组合。

4.3 总功率相同下不同风机组合的模拟研究

为了进一步进行功率组合的优化,选用上述所提到的P5和P3功率组合进行研究。 由上可知,合理的功率分配能够在减小总功率的条件下,有效地提高通风效率,因此2 种功率组合在确定自身总功率不变的情况下, 分别选3 种不同的功率分配方式进行对比,如表4 和表5 所示。

表4 P3 功率下3 种不同风机功率风机组合的分配Table 4 Allocation of three different fan power fan combinations at P3 power

表5 P5 功率下3 种不同风机功率风机组合的分配Table 5 Allocation of three different fan power fan combinations at P5 power

观察表4 和表5 得到,P3条件下的第3 种风机组合,由于单个风机功率局限,分配总功率无法达到145.5 kW,所以做了一些相应的调整,改变了单风机功率, 使得总功率达到146 kW, 且满足误差允许范围, 让第3 种风机组合总功率近似和其他两种一样。而P5条件下的第3 种风机组合通过直接选择关闭K6风机,来达到满足总功率150 kW 的目的。

图11 P5 和P3 不同功率组合隧道中的风量变化Fig. 11 Changes in air volume in tunnels with different power fan combinations of P5 and P3

图11(a)和图11(b)表示的是总功率P3=145.5 kW和总功率P5=150 kW 3 种风机功率组合在隧道中产生的风量。 由图11 可知, 在P3和P5的总功率条件下,第1 种风机组合是通过让同功率风机并联的形式组合而成,隧道风量产生分别为P31=234.96 m3/s,P51=232.75 m3/s。 而第2 组风机组合主要是单台同功率风机串联的组合形式, 其中产生的隧道风量为P52=231.61 m3/s,P53=234.25 m3/s。 第3 组风机组合是通过降低风机组合中固定的1 台或者是2 台风机功率来满足P3和P5总功率要求的风机组合形式,产生的隧道风量为P53=223.56 m3/s,P33=232.45 m3/s。 这说明总功率相同,同功率风机两两并联的风机组合产生的风量较多,采用这样的方法,能够一定程度的节约总功率消耗。 图11(c)为P5和P33 种组合隧道中产生风量进行对比,由图11 可知,风量曲线为两条递减的线段。 这表明P5和P3在进行的3 种同功率下不同风机功率组合分配时, 第一种的风机组合形式较为优秀,接着是第2 组,最后是第3 组。与此同时,还观察到曲线前半部分第1 组到第2 组的线段斜率较为的平滑,这表示2 组风量差距不大,P3的第1、 第2 组风机组合,风量差为0.71 m3/s,P5的第1、第2 组的风量差为1.14 m3/s。 而两线段后半部分第2 组到第3 组风机组合的曲线的斜率观察相对较大,并且P3第2 组到第3 组线段的斜率比P5的更低。 这说明通过单一风机功率减小虽然可以达到总功率降低的目的, 但严重的影响风机组合的通风效率。综上分析可知,同功率风机并联的组合形式, 即P3和P5的第1 种风机组合为该功率条件下较优的组合分配方式, 此外总功率在进行分配时,让较多的风机参与功率的分配,这样不仅能够有效的降低功耗, 而且还能极大地提高通风效率。

5 结 论

1)射流风机组合和空气幕风机组合在同个位置以6 台30 kW 风机两两并联进行模拟, 且2 种组合都能够达到瓮福磷矿穿岩洞隧道所需的通风卫生标准:需风量的范围(140.4~449.28 m3/s)和平均风速范围(2.5~8 m/s)。 同时还发现空气幕组合通风效果要优于射流风机的组合, 空气幕组合隧道中风量为244.63 m3/s, 射流风机组合隧道中的234.96 m3/s,两者相差了9.67 m3/s。

2)空气幕组合可以在总功率较低情况下替代射流风机组,通过列举7 种最为接近180 kW 射流风机组合通风效果, 且低于180 kW 总功率的空气幕组合,对比在隧道中运行时产生的风量和风速,得到表3 中的P3=145.5 kW 组合较优,P3组合在近似满足总功率180 kW 的射流风机组风量要求的情况下,功率节约可达到34.5 kW。

3)对较优空气幕组合P3进行组合中风机分配优化,得到在总功率相同的空气幕组合分配条件下,同功率的风机两两并联的设置形式通风效果较优。 另外,在功率分配时,应该尽量让较多的风机参与分配,这样可以有效地降低总功率,并且提高风机通风的效率。

猜你喜欢
射流风量风速
深海逃逸舱射流注水均压过程仿真分析
某乘用车冷却系统进风量仿真及优化
低压天然气泄漏射流扩散特性研究
超超临界660 MW机组二次风量异常下降分析与研究
邯郸市近46年风向风速特征分析
基于最优TS评分和频率匹配的江苏近海风速订正
基于时间相关性的风速威布尔分布优化方法
经坊煤业八采区总需风量计算研究
快速评估风电场50年一遇最大风速的算法
射流齿形喷嘴射流流场与气动声学分析