3. 噪声传播模型：几何发散、大气吸收、地面效应、屏障衰减、ISO 9613

各位好，我是老张。今天咱们聊聊噪声传播模型。说实话，这是整个风电场噪声预测里最核心的部分，也是我这些年花精力最多的地方。你想想看，风机噪声源强测出来了，但声音从轮毂高度传到敏感点，中间经历了什么？衰减了多少？这就要靠传播模型来算。

我个人习惯把传播模型拆成五个部分来理解：几何发散、大气吸收、地面效应、屏障衰减，最后用ISO 9613串起来。咱们一个一个说。

3.1 几何发散：声音的“摊大饼”效应

几何发散是最基本的衰减机制。说白了，就是声音从点源往外传，能量均匀分布在越来越大的球面上，单位面积上的声能自然就少了。

对于点声源，几何发散衰减量用这个公式算：

Adiv = 20 × log10(d / d0) + 11

其中d是距离（米），d0是参考距离（通常取1米）。

举个例子：距离风机200米处，几何发散衰减量大约是：

Adiv = 20 × log10(200 / 1) + 11 = 20 × 2.301 + 11 = 57.02 dB

嗯，这里要注意：这个公式适用于自由场条件。实际项目中，如果地面是硬质反射面（比如水泥地），衰减量会小一些，因为反射声会叠加。我在内蒙古一个项目就遇到过这种情况，风机安装在硬化平台上，实测噪声比预测值高了2-3 dB，后来排查才发现是地面反射搞的鬼。

3.2 大气吸收：空气不是真空

声音在空气中传播，能量会被空气分子吸收，转化成热能。这个衰减跟频率、温度、湿度都有关系。

大气吸收系数α（dB/km）可以用ISO 9613-1查表得到。我一般用这个简化公式估算：

α = 0.005 × f² / (1 + 0.003 × f²) × (T / 293)^(-1.5) × (RH / 50)^(-0.5)

其中f是频率（Hz），T是温度（K），RH是相对湿度（%）。

举个例子，在15°C、70%湿度条件下，不同频率的衰减系数：

频率 (Hz)	63	125	250	500	1000	2000	4000
α (dB/km)	0.1	0.3	0.7	1.5	3.0	6.0	12.0

你看，低频衰减很小，高频衰减很大。这就是为什么远处听到的风机噪声主要是低频“嗡嗡”声——高频成分早就被空气吃掉了。

3.3 地面效应：草地和水泥地差多少？

地面效应是声音在地面和声源之间来回反射产生的干涉现象。说白了，就是直达声和地面反射声叠加，有的频率加强，有的频率抵消。

地面效应衰减量Agr跟地面类型关系很大：

• 硬地面（水泥、沥青、水面）：反射强，地面效应小，Agr ≈ 0 dB
• 软地面（草地、农田、森林）：吸收强，地面效应大，Agr可达3-6 dB
• 混合地面：按比例加权平均

ISO 9613-2里给出了地面因子G的概念：

• G=0：硬地面
• G=0.5：混合地面（比如一半草地一半水泥）
• G=1：软地面

地面效应衰减量计算公式：

Agr = 4.8 - (2 × hm / d) × (17 + 300 / d)

其中hm是声源到接收点的平均高度（米），d是距离（米）。

我记得在河北一个平原项目，敏感点周围全是农田（软地面），地面效应贡献了约5 dB的衰减。但到了冬天，地面结冰变硬，衰减量直接少了2 dB。所以季节性变化也得考虑进去。

3.4 屏障衰减：山体和围墙能挡多少？

屏障衰减是噪声控制里最常用的手段。声波遇到障碍物（山体、围墙、建筑物）时，会发生衍射，一部分能量绕过去，一部分被反射或吸收。

屏障衰减量Abar用菲涅尔数N来算：

N = 2 × δ / λ

其中δ是声程差（直达声路径与绕射声路径之差），λ是波长。

衰减量近似公式：

Abar = 10 × log10(3 + 20 × N)   (当N ≥ 0时)

举个例子：一个3米高的围墙，声源高度10米，接收点高度1.5米，距离50米，对于500 Hz的噪声：

δ ≈ 0.8 米
λ = 340 / 500 = 0.68 米
N = 2 × 0.8 / 0.68 ≈ 2.35
Abar = 10 × log10(3 + 20 × 2.35) = 10 × log10(50) ≈ 17 dB

嗯，这里要注意：屏障衰减对高频效果好，对低频效果差。因为低频波长长，容易绕过去。我做过一个项目，用5米高的围墙挡风机噪声，实测125 Hz只衰减了5 dB，但2000 Hz衰减了20 dB以上。

3.5 ISO 9613：把一切串起来

ISO 9613是目前国际上最常用的户外声传播计算标准。它把上面说的几种衰减机制整合成一个完整的预测模型。

核心公式很简单：

Lp = Lw - Adiv - Aatm - Agr - Abar - Amisc

其中：

• Lp：接收点声压级
• Lw：声源声功率级
• Adiv：几何发散衰减
• Aatm：大气吸收衰减
• Agr：地面效应衰减
• Abar：屏障衰减
• Amisc：其他衰减（如树林、工业区等）

ISO 9613-2里还规定了气象修正：当风速3-8 m/s、温度梯度适中时，用“ favourable conditions ”（有利传播条件）的修正值。说白了，就是顺风时噪声传得更远，逆风时衰减更大。

下面我用SVG画一张图，把整个传播模型的结构展示清楚：

你看，整个模型就是一条减法链。从声源功率级出发，减去各种衰减，最后得到接收点的声压级。每个衰减项都有对应的物理机制和计算公式。

实际项目中，我一般用Python写个脚本，把ISO 9613的整个计算流程自动化。核心代码大概长这样：

def iso9613_predict(Lw, freq, d, T, RH, G, hm, barrier_height=None):
    """
    ISO 9613 噪声传播预测
    Lw: 声功率级 (dB)
    freq: 频率 (Hz)
    d: 距离 (m)
    T: 温度 (°C)
    RH: 相对湿度 (%)
    G: 地面因子 (0-1)
    hm: 平均高度 (m)
    barrier_height: 屏障高度 (m), 可选
    """
    # 几何发散
    Adiv = 20 * np.log10(d) + 11
    
    # 大气吸收 (简化)
    T_k = T + 273.15
    alpha = 0.005 * freq**2 / (1 + 0.003 * freq**2) * (T_k/293)**(-1.5) * (RH/50)**(-0.5)
    Aatm = alpha * d / 1000
    
    # 地面效应
    Agr = 4.8 - (2 * hm / d) * (17 + 300 / d)
    Agr = max(0, Agr * G)  # 按地面因子缩放
    
    # 屏障衰减 (如果有)
    Abar = 0
    if barrier_height:
        # 计算菲涅尔数
        delta = calculate_path_difference(d, barrier_height)
        lam = 340 / freq
        N = 2 * delta / lam
        if N >= 0:
            Abar = 10 * np.log10(3 + 20 * N)
    
    Lp = Lw - Adiv - Aatm - Agr - Abar
    return Lp

好了，关于噪声传播模型就聊到这儿。ISO 9613是个好工具，但记住它也有局限性——比如复杂地形、强风条件下的预测误差会大一些。做项目时，最好结合实测数据做校准，别完全依赖理论计算。