《Eng》：Inverse Acoustic Localization of Pile-Driving Noise from Synchronized Multi-Receiver SPL Measurements
编辑推荐：
施工噪声评估通常假定声源位置已知，这限制了其在无法接近或无法精确测量主动设备时的适用性。本研究提出了一种透明的逆声学模型，用于根据同步现场测量联合估计打桩声源的位置和随时间变化的声压级（SPL）。研究人员在哈萨克斯坦阿斯塔纳的一个建筑工地，在五个全球导航卫星系
施工噪声评估通常假定声源位置已知，这限制了其在无法接近或无法精确测量主动设备时的适用性。本研究提出了一种透明的逆声学模型，用于根据同步现场测量联合估计打桩声源的位置和随时间变化的声压级（SPL）。研究人员在哈萨克斯坦阿斯塔纳的一个建筑工地，在五个全球导航卫星系统（GNSS）定义的接收器位置记录了5秒的噪声，每个接收器以0.125秒的间隔产生40个声压级值。该模型应用球面几何扩散，并通过最小化各接收器获得的距离校正声源级之间的差异来识别声源位置。估计声源位于局部坐标[56.6, 113.8] m处，蒙特卡洛敏感性分析表明，在接收器坐标和声压级测量联合不确定度下，95%置信区域约为5–8 m。重建的声源声压级在1 m参考距离处介于81.6至99.9 dB(A)之间，等效声源声压级为90.4 dB(A)。测量的脉冲模式在重建的声源历史中得以保留，主要峰值约以1 s间隔出现。测量与预测接收器声压级比较得出总体平均绝对误差（MAE）为1.3 dB(A)，均方根误差（RMSE）为1.6 dB(A)，反映了测量与重建接收器声压级之间的内部一致性。随后对实际打入桩的位置进行了独立测量，估计位置的真实定位误差为2.3 m，在上述95%置信区域内。样本外交叉验证和搜索区域内的合成扫描进一步表征了该方法的可辨识性及其在搜索区域中心之外的预期敏感性。
**打桩噪声源逆声学定位研究解读**
**研究背景与问题**
建筑施工噪声是城市环境中变化最剧烈、干扰最强的噪声源之一，其位置、持续时间、设备组成和运行状态在项目全生命周期内不断变化。与相对稳定的交通或工业噪声不同，施工噪声通常具有间歇性、临时性、空间移动性，并以冲击、物料搬运、挖掘和重型机械产生的突发峰值为特征。打桩作业尤为特殊，它产生一系列短暂、高能量的脉冲而非稳态声信号，对周边社区造成显著干扰。传统施工噪声评估通常假定声源位置已知，这限制了其在无法接近或无法精确测量主动设备时的适用性。现有施工噪声模型大多在声源位置和设备特性确定后预测声级、暴露量、成本或缓解需求，而成熟的声学定位方法则依赖于麦克风阵列、波形时序、相位差、方向信息或详细的环境数值表示。目前缺乏一种简化的逆场问题求解方法：仅利用多个空间分布声级计同步测量的标量声压级（SPL，Sound Pressure Level），在声源位置未知的情况下，联合估计声源位置及其时变声学强度。
**研究目的与意义**
本研究旨在开发和现场测试一种透明的逆声学数学模型，利用五个GNSS定义接收器位置的同步短期测量，定位未知打桩声源并估算其时变声压级。该模型将打桩视为5 s区间内的固定点声源，允许声源声压级以0.125 s时间步长变化，通过搜索施工场地内使五个接收器距离校正后的声源级估计差异最小的位置来确定声源坐标。该研究扩展了基于接收信号强度（RSS，Received Signal Strength）的定位范式在施工噪声监测中的应用，其创新点不在于传播模型本身（为标准球面扩散关系），而在于利用低成本非阵列声级计数据联合估计时变声源强度与位置，专门应用于脉冲性打桩噪声。研究成果发表在《Eng》期刊。
**关键技术方法**
研究人员在哈萨克斯坦阿斯塔纳某在建工地打桩作业期间进行了现场测量，五个测量位置（R1至R5）围绕施工活动布置，使用UT352声级计（UNI-T，东莞，中国）以A计权、FAST时间响应（每读数0.125 s）、自动量程模式采集数据。由于设备不提供自动数据记录功能，研究人员使用智能手机对每个声级计屏幕进行视频记录，并通过同一视频应用叠加网络同步系统时钟时间戳建立共同时间基准。测量持续5 s，每个位置获得40个声压级值。坐标通过GNSS接收器确定并转换至局部笛卡尔坐标系（WGS84椭球等距圆柱投影）。模型假设球面几何扩散为声压级衰减主要机制，通过网格搜索最小化目标函数确定声源位置，随后估算时变声源声压级。不确定性分析采用蒙特卡洛方法，分别扰动接收器坐标（标准差3 m）和声压级测量值（标准差1.5 dB），并采用与声级相关的测量不确定度模型。
**研究结果**
*3.1 测量声压级概述*
五个接收器在5 s测量期间各获得40个同步声压级记录，所有时间历程呈现相似的循环模式，具有重复的短期升降，与打桩作业的脉冲特性一致。主要峰值出现在约0.625、1.625、2.625、3.625和4.625 s处，表明主脉冲间隔约1 s的准周期冲击序列。R1记录最高声压级值（最大69.7 dB(A)），R3记录最低总体水平（最大58.9 dB(A)）。等效声压级遵循相同模式，表明声源更靠近R1、R5和R2，远离R3和R4。
*3.2 估计声源位置与声源到接收器距离*
目标函数图显示声源估计位于搜索区中心区域，五个接收器位置之间。目标函数最小值在x = 56.6 m和y = 113.8 m处获得。计算声源到接收器距离分别为：R1为42.1 m，R2为48.8 m，R3为98.5 m，R4为70.6 m，R5为47.4 m。该结果与测量声压级分布一致，R3具有最大源-接收器间距，解释了其相对较低的测量值。随后使用相同GNSS接收器独立测量实际打入桩位置（局部坐标[58.2, 115.4] m），估计位置的真实定位误差为2.3 m。
*3.3 估计声源声压级*
确定最优声源位置后，对每个时间步反算声源声压级。估计声源声压级在1 m参考距离处介于81.6 dB(A)至99.9 dB(A)之间，算术平均值为87.2 dB(A)。最高估计声源声压级出现在1.625 s。5 s区间等效声源声压级为90.4 dB(A)。五个接收器单独反算的声源声压级估计值彼此紧密跟踪并共享相同脉冲模式，其均值作为单一随时间变化的声源估计。
*3.4 重建一致性检验结果*
使用估计声源坐标和声源声压级时间序列计算五个接收器位置的预测声压级，预测值与测量值在量级上表现出良好一致性，时间峰值也通过时变声源声压级估计得以重现。整个数据集（200个样本）总体MAE为1.3 dB(A)，RMSE为1.6 dB(A)。各接收器MAE范围为1至1.6 dB(A)，RMSE范围为1.4至2.2 dB(A)。限制为每个接收器主导脉冲峰值的5个样本（25个样本）时，相关系数R2 = 0.939，略高于完整数据集的R2 = 0.913，表明模型与测量数据的一致性不是脉冲内紧密间隔自相关样本的伪影。R3同时具有最大MAE（1.6 dB(A)）和RMSE（2.2 dB(A)），与其距估计声源距离更大一致。
*3.5 敏感性与不确定性分析结果*
在仅接收器坐标不确定度（σ = 3 m）下，估计声源位置标准差在x方向为2.3 m、y方向为2.6 m，95%置信区域半轴为4.5 m和7.3 m，95百分位径向位移为6.1 m。仅声压级测量不确定度（σ = 1.5 dB）下，位置敏感性较低，标准差分别为1.0 m（x）和1.2 m（y），95%置信区域半轴为2.0 m和3.1 m。采用与声级相关的不确定度模型（σ
SPL = σ? + k·(SPL
min ? SPL)，其中σ? = 1.5 dB，SPL
min = 45.1 dB(A)）后，样本间标准差在响亮时刻为0.66 dB至脉冲间安静时刻1.25 dB之间变化。同时施加两种不确定度时，95%置信区域半轴约5.1 m和8.0 m，95百分位径向位移6.6 m。源位置对假设几何扩散指数的敏感性比测量误差敏感性更显著：指数在1.85–2.15范围内变化时位置偏移1.5–1.6 m；指数为1.5或3.0时偏移约6–8 m。
*3.6 搜索区内源位置敏感性*
在相同模型框架内进行合成验证，源位置覆盖整个搜索区域包括边缘和角落，使用现场重建的源级历史作为代表性脉冲模式生成合成接收器数据。定位误差从搜索区中心到边缘和角落增加三到五倍，确认单一中心定位的现场事件不代表整个搜索区域的性能。对现场数据集应用传播延迟校正后，估计源位置偏移约6.5 m，定位误差为8.5 m，大于未校正基线估计的2.3 m。仅使用每个主导脉冲峰值附近3个样本（共15个样本）进行反演时，定位误差为6.3 m，大于使用完整40个样本序列的2.3 m。
*3.7 样本外交叉验证*
采用留一接收器（LORO，Leave-One-Receiver-Out）交叉验证：依次省略五个接收器之一，仅使用剩余四个接收器重新估计源位置和源级历史，并预测被省略接收器的测量声压级。样本外MAE和RMSE分别为3.1 dB(A)和3.8 dB(A)，高于样本内值（1.3 dB(A)和1.7 dB(A)），确认样本内指标低估了真实独立数据上的误差。
*3.8 与基线定位方法比较*
将所提方法与两种标准RSS基线方法比较：加权质心定位和RSS/路径损耗最小二乘定位。所提方法在平均绝对重建误差上优于加权质心，与RSS最小二乘相当，RMSE更低表明大误差更少。相对于独立测量的地面真值位置，所提方法定位误差（2.3 m）约为RSS最小二乘（5.8 m）的2.5倍更低、加权质心（11.3 m）的5倍更低。
**讨论与结论**
讨论部分指出，五个同步分布接收器的声压级测量可用于估计打桩声源的位置和时变声学强度。测量记录呈现在约1 s间隔的五个主导峰值，与先前受控传播实验观察一致。估计等效源声压级90.4 dB(A)反映了打桩的强非平稳特性，单一恒定源值不足以描述该事件。总体MAE为1.3 dB(A)和RMSE为1.6 dB(A)表明简化模型在样本内一致性检验下以低残差重现测量接收器声级，但该值不与其他施工噪声模型的预测误差直接比较。模型的主要局限包括：忽略地面反射、屏蔽、散射、气象影响等传播机制；假设瞬时传播，而实际源-接收器距离差56 m对应约0.16 s传播时间差，超过一个采样间隔；忽略打桩冲击点与地面接收器间的垂直偏移（对最近接收器影响最多约0.9 dB(A)）；排除地面传播振动再辐射机制；假设单一主导声源。相对于传统波束形成、时间反转、有限元逆方案等方法，所提方法的主要优势在于操作简单性，仅需同步标量声压级值、已知接收器坐标和简单目标函数。
研究结论：本研究开发并现场测试了一种透明逆声学模型，利用五个空间分布接收器的同步测量联合估计未知打桩声源的位置和时变声压级。应用于包含每接收器40个测量的5 s施工场地事件，模型将声源定位在[56.6, 113.8] m处，在组合坐标和测量不确定度下95%置信区域约5–8 m，估计等效源声压级为90.4 dB(A)（1 m参考距离）。重建源历史保留了打桩的脉冲和准周期特性，峰值达99.9 dB(A)。测量与预测接收器声级一致性良好，总体MAE为1.3 dB(A)，RMSE为1.6 dB(A)。独立测量的桩体地面真值坐标给出估计位置的真实定位误差2.3 m。这些发现证明无需麦克风阵列波形数据、复杂环境模拟或大量训练数据集，即可从标量声压级测量实现有用的施工噪声源定位和强度估计。该方法为快速现场评估和未来集成至自动化传感器及GIS监测系统提供了实用基础。更广泛的验证应包括其他场地和事件的独立测量源坐标、更长测量时间、不同接收器布置、多个施工场地以及气象条件、地面效应、屏障、反射、接收器网络传播延迟和多源或移动源修正。
