工程中三維有限長聲屏障的優化設計及應用

阮學云，邵良友，章林柯，魏 玥，李 達，許祥濤

（1.安徽理工大學 機械工程學院，安徽 淮南 232001；2.武漢理工大學 能源與動力工程學院，武漢 430063）

隨著城市和工業化進程的推進，我國環境問題日益凸顯，噪聲已經成為一大焦點。長期接觸比較強烈的噪聲，可以引起病理性改變，還會給人帶來一些其他危害。聲屏障因其經濟有效，常作為主要的一種降噪措施。常規聲屏障的降噪效果往往會受到固定幾何形狀限制，Yu等[1]基于Miura-ori晶胞構建三維屏障結構，通過折疊其形狀改變聲影區和衍射路徑；吳小萍等[2]基于NSGA-Ⅱ算法對高速鐵路聲屏障高度進行優化設計；Liu等[3]提出了二維聲學中的等幾何快速多極邊界元方法以及基于相干靈敏度的聲屏障形狀優化算法；Toledo等[4]基于進化法提出一種通過將薄平面聲屏障理想化為具有零邊界厚度的輪廓來優化薄平面聲屏障的整體形狀和頂部邊緣的方法；Mun 等[5]提出了一種使用模擬退火（SA）算法的全局優化設計所需噪聲的新方法；陳磊磊等[6]采用移動近似算法（MMA）進行二維聲屏障結構形狀的優化分析；徐圣輝等[7]研究了Cadna/A軟件在高速鐵路聲環境影響評價中的應用。目前國內外提出的大多是優化聲屏障形狀，很少涉及到優化聲屏障結構尺寸。

1 有限長聲屏障衰減理論及算法

1.1 聲屏障降噪原理以及二維聲屏障衍射路徑

在聲源和接受點之間加入一個有一定面密度的物體，這樣在聲音的傳播過程中就會產生衰減現象,這個物體就成為“聲屏障”。聲波遇到屏障時將產生反射、透射和衍射三種傳播現象如圖1。一部分從屏障上方、左右兩端繞射到接受點，一部分直接透過聲屏障直接到接受點，另一部分在屏障處反射。其中衍射是聲屏障達到降噪效果的主要手段，是其通過延長傳播路線來增加聲衰減。

圖1 聲波遇到屏障的三種傳播現象

聲波傳播到屏障后方的區域稱之為“聲影區”，聲源經過聲屏障直行的路徑與聲屏障到接受點路徑的夾角為繞射角，一般繞射角越大，聲衰減越多[8]。

S′和R′是聲源和接受點關于地面的投影。其主要四條路徑為：SO′R，SO′BR，SAO′R，SAO′BR，如圖2所示接受點聲壓為四條路徑聲壓的矢量和。二維聲屏障聲壓計算只考慮了聲屏障的高度和形狀因素，三維聲屏障考慮了長度因素，即側面繞射聲。

圖2 二維聲屏障衍射路徑

1.2 三維有限長聲屏障衰減的計算方法

目前工程常用的是由前川純一歸納出的以無量綱量菲涅爾數N為參量的計算方法。

菲涅爾數N定義為：

其中：δ為聲程差，如圖所示δ=A+B-D，λ為聲波波長（m），f為聲波頻率（Hz），c為聲速，一般取340 m/s。

Kurze[9]通過研究實驗結果圖表，利用Kirchhoff積分定理，得出一個較為簡單的公式：

本文基于這個公式進行聲屏障衰減計算。在現實中聲屏障都是有限長。對于三維聲屏障，在計算過程中，需要計算8 條主要的繞射路徑對受聲點聲壓的貢獻[10]，其8條傳播路徑模型建立如圖3。

圖3 三維有限長聲屏障的主要衍射路徑

其中以屏障底邊中點為空間坐標系的原點，且底邊與Y軸重合，垂直于XOY面，聲屏障長度為x2，高度為x3。聲源S高度為h1，距離屏障距離為x1；接受點高度為h2，直線SR長度為D，且與屏障的交點為J點。S為聲源，R為接受點。圖3所示的路線1、2、3、4為聲波從聲屏障上方繞射的最短路徑[11]，路線5、7為聲波從兩端繞射到達接受點，路線6、8為聲波經過地面反射到達接受點。設插入聲屏障之前R點的聲壓級為Lp0，R點坐標（a，b，h2），根據以上信息經計算得出以下各點坐標：

根據空間兩點之間距離公式可得出各路徑長度，具體如下：

8條路徑的聲程差δ由以下公式表示：

菲涅爾數N可表示為：

8條繞射路線的插入損失ΔLd(i)：

因此插入聲屏障之后聲波8 種繞射路徑R點的聲壓級Lpi：

接受點R點某一頻率的聲壓級Lpf為8種路徑通過能量法疊加在一起，具體如下表示：

實際R點的聲壓級Lptot為倍頻帶中心頻率處(63 Hz～8 000 Hz)各個頻率聲壓級的疊加：

因此實際的綜合插入損失ΔLd：

2 結合聲屏障工程案例的優化設計模型

2.1 工程及噪聲源介紹

江西某水泥廠，擁有2條5 000 t/d、1條3 500 t/d自動化水泥生產線。由于水泥生產工藝復雜，設備較大、功率較高，所以在運行生產過程中會產生較大的噪聲，因此該水泥廠的廠界噪聲值達到65 dB（A）以上，遠超過50 dB，不符合國家標準GB3096－2008《城市區域環境噪聲標準》和GB12348－2008《工業企業廠界環境噪聲排放標準》的二類噪聲排放要求。在本節將會利用其中一個聲源設備使用上節中提到的算法進行分析，優化聲屏障結構尺寸，并使用Cadna/A 聲學預測軟件模擬驗證。立磨如圖4，作為水泥廠生產線的第一道工序，其噪聲輻射較大。圖5所示是噪聲測量現場，使用聲級計測量立磨倍頻帶中心頻率處頻譜，頻譜圖如圖6所示。插入聲屏障前接受點倍頻帶中心頻率處的聲壓級如表1。

圖6 立磨倍頻帶頻譜圖

表1 插入聲屏障前接受點各頻率的聲壓級

圖4 立磨現場圖

圖5 頻譜測量現場圖

插入聲屏障之前接受點的聲壓級Lp0為70 dB。因其對東邊廠界敏感點影響較大，需要重點治理。聲源主要位于底部的電機、油泵、中部的磨輥、頂部的落料及轉動電機。聲源與接受點直線距離為100 m，遠大于聲源幾何尺寸的兩倍，因此將聲源等效為一個點聲源[12]，聲源高度為8 m。

針對以上情況決定采用聲屏障將其圍擋，從而進行降噪處理。其示意圖如圖7所示。

圖7 立磨聲屏障安裝示意圖

2.2 確定目標函數及設計變量

本文以經濟性為目標函數[13]。其中包括聲屏障成本，鋼結構和基礎費用。聲屏障成本主要是由其面積決定，鋼結構和基礎費用由長度決定，兩者比例為1：α。

因此有限長的聲屏障優化設計的目標函數為：

其中：α為聲屏障成本與鋼結構和基礎費用的比例，本文α取0.75，a、b分別表示聲屏障的長和高，m為實際中每平方的單價，根據實際工程報價，m取1 000，單位為元。

考慮到聲源與聲屏障的距離影響到接受點的聲壓級，本文將以下幾個參數設為設計變量[14]：聲源與聲屏障的距離x1，聲屏障長度x2，聲屏障的高度x3，故設計變量為X=[x1,x2,x3]T。示意圖如圖8。

圖8 設計變量示意圖

根據上述分析，目標函數可簡化為：

2.3 約束條件

設計變量約束：聲屏障類型、建筑界限、聲源高度等。立磨附近還有其他設備噪聲，為了不影響結果，在立磨設備檢修時，我們測得其他設備噪聲，將其作為背景噪聲。由于施工建筑界限，聲屏障與聲源距離要大于1.5 m小于3.8 m。聲源高度為8 m，故聲屏障的高度應大于8 m。

聲學約束：接受點聲壓級滿足《聲環境質量標準》（GB12348－2008）標準，即接受點聲壓級小于等于50 dB。由下式表示：

其中：Lp0表示為聲屏障插入前接受點聲壓級，LpR表示為聲屏障插入后接受點聲壓級，?Lp表示為插入損失。本文已知Lp0=70 dB，因此插入損失?Lp=Lp0-50=20 dB。

本文將其表示為：

2.4 有限長聲屏障優化設計數學模型的一般形式

根據前述，建立如下的聲屏障優化設計數學模型[15]：

2.5 利用MATLAB優化工具箱求解

已知條件整理如下：

根據MATLAB 非線性不等式約束fmincon[16]函數優化求解，主程序如下：

A=[-1 0 0;1 0 0;0-1 0;0 0-1];

b=[-1.5;3.8;-20;-8];

vlb=[0;0;0;];

x0=[3,38,14];

[x,fval,exitflag,output,lambda,]=

fmincon('funn',x0,A,b,[],[],vlb,[],'non')

最終的最優解：

X*= (2.525,41.739,20.869)T

取整后得：

X*= (2.5,42,21)T

2.6 Cadna/A模型驗證

為了驗證優化設計后與優化前相比聲屏障降噪的效果，本小節利用Cadna/A 噪聲預測軟件進行模擬分析，優化前的位置以及聲屏障的尺寸參數是根據水泥廠現場工作人員經驗用噪聲預測軟件試湊法獲得。將優化前后聲屏障的位置以及結構尺寸參數輸入進去，獲得優化前后的聲場分布圖，如圖9、圖10所示。

圖9 聲屏障優化前聲場分布圖

2.7 現場實驗驗證

為了進一步驗證聲屏障的優化結果，現場利用聲級計測量插入優化后聲屏障接受點倍頻帶中心頻率處的聲壓級，如表2所示，其聲壓級為48.9。

表2 插入聲屏障后接受點各頻率的聲壓級

2.8 結果分析

圖9表示優化前聲場分布圖，其接受點聲壓級為48.0 dB，圖10 表示優化后聲場分布圖，其接受點聲壓級為48.2 dB，在現場測得插入優化后聲屏障接受點的聲壓級為48.9 dB，滿足預期效果。表3給出了優化前和優化后各參數（x1,x2,x3）以及所需成本的對比，優化后建造的聲屏障所需費用比優化前節省了29.4%。

表3 屏障優化前后插入損失成本對比

圖10 聲屏障優化后聲場分布圖

3 結語

基于有限長聲屏障衰減的工程算法，以聲屏障綜合成本為目標函數，在滿足一定的插入損失情況下，優化的設計方法相比傳統使用噪聲預測軟件試湊法更加科學，可以降低聲屏障的成本。結果表明，該優化方法今后在聲屏障的實際工程中具有較大的應用前景。