阮学云,魏 玥,高 艳,李 达,王 相
(1.安徽理工大学 机械工程学院,安徽 淮南 232001;2.中国科学院 声学研究所噪声与振动重点实验室,北京 100190)
廊道作为电厂燃煤、水泥厂石灰石等输送的主要通道,往往距离地面一定高度且布置于企业厂界外。由于其工作时间长、噪声范围广、中低频段传播距离远,对附近居民生活影响较大,是主要的投诉对象。廊道噪声以托辊[1]、带式输送机[2]等噪声源为主,治理方案主要是针对廊道整体进行声屏障包围[3-4]或整体更换低噪声托辊。根据经验确定控制方案,无科学依据。目前Cadna/A 等噪声预测工程软件模拟悬浮声屏障只考虑上方绕射路径,对廊道底部增加声屏障等各组合方案无法进行准确的噪声预测,影响了模型预测计算准确性,从而无法得到科学的噪声控制方案。
国外对声屏障降噪的研究自20 世纪60 年代开始,国内对声屏障降噪的研究起步相对较晚。Mackawa[3]研究出二维声屏障的插入损失计算公式;Zaplaic 等[4]提出一种计算模型避免对运送水泥原料的传送带附近噪声评估时的背景噪声干扰;徐圣辉等[5]利用边界元对近轨不同型式上挑檐进行对比研究;王金瑞等[6]运用ISO9613-2 标准对两种类型的道路声屏障消声效果进行计算与对比;郭萍等[7]在二维声屏障声压计算方法的基础上给出了三维有限长声屏障声压简化算法。以上研究均针对廊道上挑檐运用有限元或边界元算法[8-10],计算复杂、耗时长,不适合工程应用。基于Cadna/A 等工程预测软件对声屏障效果的预测多针对道路等交通噪声[11-13],较少有针对悬浮声源下挑檐绕射声进行相关研究。
针对以上情况,基于Cadna/A 预测软件为建模工具,通过声学元素几何变化方式,实现廊道各类型悬浮声屏障噪声控制方案计算,为科学选择噪声控制方案提供理论依据。
Cadna/A 软件依据ISO9613-2 标准对环境噪声进行模拟与预测,并获得国家认证[14]。根据廊道声源特性,以点声源为例建立测试模型,验证目前噪声预测模型针对悬浮屏障下挑檐的计算结果准确性。其中:声源声功率级80 dB,声源高度10 m,声屏障宽度均为3 m,在垂直声屏障竖面建立高20 m 的垂直声场,观察单侧声屏障与下挑檐声屏障的声场分布,如图1 所示。
图1 单侧及下挑檐声屏障声场分布图
测试结果表明,设置下挑檐模型与单侧声屏障声影区声场分布一致,并经过单点计算,计算声源水平距离20 m、30 m、50 m 及高1.5 m 处噪声值,整理结果如表1 所示。
表1 Cadna/A 软件计算受声点噪声值 (单位:dB/A)
由图1 和表1 可知,Cadna/A 软件计算中,虽可以设置下挑檐型式,但因其软件计算模块对下方绕射声忽略考虑,声影区受声点的数值均没有发生变化,其下挑檐实际不起作用。因此,Cadna/A 预测软件对下挑檐绕射声计算不精确。
考虑预测软件可计算上挑檐绕射声,为将下挑檐计算转化为上挑檐计算,建模时可通过保证声源与受声点声程差不变的情况下转换绕射声路径。根据几何关系,可转换为上绕射进行等效计算,示意图如图2 所示。
图2 下挑檐路径转换示意图
以声源与受声点的水平连线设置为X 轴,垂直于地面的单侧声屏障为Z 轴建立坐标系,声源S(sx,sy,sz),受声点R(rx,ry,rz),声屏障下挑檐绕射点a(ax,0,az),保持S 和R 的位置不变,在XOZ坐标面经过几何关系转换成上挑檐后,保证几何变换后声程差z 不变:
式中:b 为 声源与受声点之间的距离平行于屏障上边界的分量,m;
dSa为声源至下挑檐绕射路程,m;
daR为受声点至下挑檐绕射路程,m;
dSR为声源至受声点的路程,m。
声屏障绕射路径之和:
经过路径的距离计算得出:
根据公式(3)及S,R 和a 点坐标计算出新的绕射点a0(ax0,0,az0)。
根据ISO9613-2 标准[10],屏障绕射声衰减计算公式为:
式中:C2=40 不考虑地面反射影响;
对于单绕射:C3=1 ;
λ为标称频带中心频率的声波波长;
Kmet为气象影响修正因子:
式中:dss为声源到第一绕射边的距离;
dsr为第二绕射边到接收点的距离;
d 为声源到受声点的距离;
z 为声程差;
增设声屏障后受声点声压级:
LW为声源产生的声功率级;
DC为指向性校正;
A 为声源至接收点倍频带衰减;
其中:
Ad为几何发散引起的衰减;
Aa为大气吸收引起的衰减;
Ag为地面效应引起的衰减;
Ab为加屏障引起的衰减;
为其他多方面效应引起的衰减;根据上述式(1)-(7)计算出声波经过下挑檐后受声点处声压级,并利用几何关系通过Cadna/A进行等效计算,计算结果如表2。
表2 下挑檐软件计算与理论计算 (单位:dB/A)
根据上述理论计算与软件计算结果对比,误差均在0.5 dB 以内,因此该等效建模方法可实现Cadna/A 软件中下挑檐绕射声的预测计算。
根据前节提供的方法,对安徽某电厂廊道进行噪声计算,并根据常见的噪声控制方案进行噪声预测分析,给出科学治理方案。该廊道沿途经过一处居民区,如图3 所示,受影响居民较多,为了改善该居民区声环境质量,通过新的几何建模方法,基于Cadna/A 建立廊道周边区域的噪声预测模型,分别计算廊道声屏障的上、下绕射路径,预测声屏障及低噪声托辊对居民区的降噪效果,给出合理可行的降噪措施。
图3 廊道卫星地图
廊道的主要噪声来源皮带输送机呈圆柱形,距地面高度11 m,横截面直径约1 m;居民区距声源水平距离约100 m,在此例中将廊道噪声源等效为线声源。廊道全程2 200 m,居民区受声点高度1.5 m。廊道声源频谱采用廊道近距离实测频谱。
以廊道西边居民区为研究对象,使用Cadna/A软件对该区域进行预测计算,治理前噪声网格示意图如图4 所示。
图4 廊道及居民区治理前噪声网格地图
根据噪声网格地图整理出各敏感点噪声值如表3。
表3 治理前各位置噪声值 (单位:dB/A)
治理前选择居民区距离廊道最近的两处敏感点进行实测,误差均在1 dB 之内,表明了该预测模型的可靠性。
根据治理前预测结果,显示小区内部噪声值基本达标,小区东北角区域噪声值最大可达到59.5(dB/A),根据声环境质量执行《声环境质量标准》(GB12348-2008)2 类标准[15],夜间要求限值50(dB/A),该处敏感点超标。
根据治理前Cadna/A 软件预测结果显示,对敏感点进行治理,建立预测模型,分析三种治理方案的效果及经济性。
方案1:上下挑檐声屏障
设置声屏障宽度、上挑檐声屏障宽度、下挑檐声屏障宽度,均取3 m。根据上文方法选择式(2)和(3)计算出下挑檐几何变化后的声屏障高度,建立模型如图5 所示。
图5 上下挑檐声屏障几何模型建立
在选择设置声屏障区域利用Cadna/A 分别计算出无声屏障、带上挑檐声屏障、下挑檐声屏障在敏感点处声压级Lp1、Lp2、Lp3,根据式(6)和(7)可知:
其中:
Ab1为上挑檐引起的衰减;
Ab2为下挑檐引起的衰减;
根据式(8)-(11)得出上下挑檐声屏障在敏感点处的总声压级 Lp为:
方案2:更换低噪声托辊
低噪声托辊直径较普通托辊更大,且托辊支架厚度更厚,低噪声托辊的寿命是普通托辊的两倍。根据低噪声托辊厂家提供参数,低噪声托辊较普通托辊降低10-15 dB 左右,降噪量按照12 dB 进行设置,更换低噪声托辊。几何示意图如图6 所示。
图6 更换低噪声托辊二维模型
方案3:低噪声托辊配合上下挑檐声屏障
为进一步探究低噪声托辊及上下挑檐声屏障治理效果,运用上文计算方法,对更换低噪声托辊区域增设上下挑檐声屏障,计算其治理效果。几何示意图如图7 所示。
图7 声屏障搭配低噪声托辊二维模型
根据以上三种治理方案进行计算,运用Cadna/A 进行计算,整理结果如表4 所示。
表4 三种方案治理效果 (单位:dB/A)
根据计算结果可知三种方案均可在一定范围内使敏感点达到治理标准。
根据市场调查静音低噪声托辊的价格大约300 元,每100 m 更换低噪声托辊的价格大约20万元,安装费大约1 万元,共计21 万元每百米;增加声屏障每100 米材料费约45 万元,安装费约1万元,共计46 万元每百米。结合治理达标工程整理出预计达标成本,如表5 所示。
表5 廊道噪声治理预计达标成本 (单位:万元)
根据达标成本结果显示,更换低噪声托辊治理廊道噪声,在达标的情况下,其成本最低,噪声预测地图如图8 所示,可作为推荐方案。
图8 更换低噪声托混后噪声网格地图
针对安徽某电厂廊道噪声对周围居民区影响较大,应用上述方法对该居民区进行噪声治理,实施现场如图9 所示。
图9 更换低噪声托辊实物图
根据上述方案分析,更换1 000 m 静音托辊后对敏感点进行噪声数据采集,测试结果如下表6。
表6 治理后各位置噪声值 (单位:dB/A)
根据表6 可以看出,治理效果达到预期值,说明输送廊道噪声治理方案选择合理。
利用几何关系将声屏障下挑檐绕射声进行几何转化,通过新的几何建模方法,实现了Cadna/A准确计算各类型屏障绕射声计算。应用该方法对某电厂廊道建立噪声预测模型,对各方案计算结果进行技术与经济性对比分析,并给出了相关结论,为廊道噪声方案提供了科学依据,其结论也可为同类型输送廊道噪声方案提供参考。