简介这份资源聚焦声发射AE信号的时域特征计算面向材料无损检测、结构健康监测领域的研究者与工程师可用于识别裂纹扩展、塑性变形等典型声发射源。包内为单个MATLAB脚本通过读取声发射数据、滤波与阈值处理计算上升时间、幅度、持续时间、计数和能量等关键参数辅助评估材料内部损伤程度与结构稳定性。压缩包仅含1个m文件大小约2KB轻量易用适合需要快速验证信号处理逻辑或进行二次开发的初学者。目前已有495人学习下载。借助该脚本读者可掌握声发射信号参数化分析的基本流程理解各时域特征在故障诊断中的物理意义并在此基础上扩展出更完整的声发射数据分析方案。1. Acoustic_signal_computing 的核心不是算出波形而是算准“时间”Acoustic_signal_computing 放在声发射AE检测里翻译过来就是“声信号计算”而这个领域最磨人的参数恰恰是上升时间。做过声发射监测的人多半有这种经历一次清晰的突发型裂纹事件第一套系统显示上升时间 12μs换一套系统变成 24μs同一个断铅源反复触发波形肉眼看起来几乎一样结果特征值忽大忽小。这不是设备玄学而是整条采集链路在改写信号。上升时间又恰恰是声发射定位的基础波达时刻怎么标直接决定定位坐标偏出去多远。下面的内容写给做声发射检测、结构健康监测和声源定位的工程师目标是让你搞清楚这个参数从哪来、怎么算、怎么用才不会翻车。2. 声发射信号的上升时间是被测量系统“改写”出来的传感器谐振、放大带宽与采样率2.1 谐振式传感器决定了上升时间的下限声发射信号本身是材料内部应变能释放产生的弹性波。裂纹扩展、断铅、撞击这类突发型源物理过程可以很快快的时候几微秒甚至更短。但声发射系统记录到的波形不是这个原始弹性波而是经过传感器这一层“机械滤波器”之后的输出。绝大多数 AE 传感器是压电陶瓷谐振式换能器厂家给出的谐振频率150kHz、300kHz、500kHz对应着很高的 Q 值通带很窄。窄带系统对阶跃输入的响应输出脉冲上升时间的经验下限可以粗估为 0.35/BWBW 是 -3dB 带宽。一只 150kHz 谐振、Q 值约 5 的传感器有效带宽约 30kHz0.35/30000 ≈ 11.7μs。这意味着就算原始事件在 1μs 内完成你从波形上看到的信号前沿也会被拉伸到接近 12μs。换成 500kHz 谐振、带宽 80kHz 的传感器下限降到 4.4μs。所以在上升时间这件事上传感器谐振特性是最强的一只手它决定了你不可能看到比它更快的沿。一个更反直觉的推论是两个不同谐振频率的传感器装在同一块板上测同一个断铅源测到的上升时间可能差出一倍以上。这不代表物理事件不一样只代表系统对高频能量的响应不同。尤其在做不同批次、不同系统之间的 AE 数据对比时只要传感器没有统一标定上升时间就没有可比性别一上来就拿它当“物理过程快慢”的判据。2.2 前置放大器、采样率与记录长度都在决定“算出来的上升时间”传感器之后是前置放大器常见增益 40dB。前放的带通滤波与输入带宽对上升时间的影响常被人忽略。声发射前放大多是 AC 耦合低截止频率若设在几十 kHz对长衰减的 AE 波形底部有畸变高截止频率若只有几百 kHz带宽进一步收窄上升沿会变缓。同样一只传感器接 300kHz 高截止前放和接 1MHz 高截止前放上升时间可能差开一倍。你调的所谓“上升时间滤除”有时候其实是在滤除前放带宽。采样率是另一道闸门。上升沿如果只剩一两个采样点跨阈位置的量化误差就是一个采样周期上升时间的统计结果会乱跳。做上升时间计算和声发射定位的采集系统我一般建议不低于 5MS/s只做振铃计数和能量统计的巡检1MS/s 勉强够用但别拿那个数据去算上升时间。记录长度会影响起点回溯。上升时间测量需要一段完整包含噪声底、上升沿和峰值位置的波形窗口。只存“上升段”而不存衰减包络的触发式记录10-90% 法还能勉强用但记录长度最好至少覆盖到预计上升时间的 20 倍以上峰值位置才可信。我常用的一套典型采集参数如下参数推荐值对上升时间的影响传感器谐振频率150kHz 通途300kHz 以上做定位决定上升时间下限高频传感器沿更陡前置放大增益40dB 起步按噪声底调整过高则饱和限幅过低则触发抖动低截止/高截止20-50kHz / 1-2MHz带宽越窄输出上升沿越缓采样率不低于 5MS/s上升时间/定位场景上升沿采样点不足时跨阈量化误差大ADC 位数不低于 12bit常用 16bit峰值与噪声底量化精度影响 10-90% 计算记录长度至少为上升时间的 20 倍太短会截峰影响峰值与起点回溯注意表里的推荐值是常见 AE 检测系统的经验区间不同厂家传感器和采集卡有各自规格选型以手册参数为准。2.3 选传感器时怎么权衡上升时间与传播距离选高频传感器还是低频传感器不只是灵敏度问题。谐振频率高上升时间小时间分辨率好对声发射定位有利但高频成分在结构里衰减快传感器离源越远信号可能小到抬不起头。低频 150kHz 系统能听到更远但上升时间被拉伸到 10μs 级通道之间的到时差如果只有几十微秒这个系统偏差就不能忽视。常见的折衷做法是 300kHz 传感器做定位150kHz 做大面积巡检。如果系统里同时混着 150kHz 和 300kHz 传感器它们输出的上升时间天生不在一个量级。做统一阈值触发时请分通道设定触发阈值别用一个全局值否则上升时间这个特征在后端分类里就变成纯噪声。3. 上升时间的软件计算起点判定才是真正决定结果的那只手3.1 三种上升时间定义声发射为什么用“阈值到峰值”先理清定义。电子学里上升时间一般指 10%-90%信号从 10% 幅度上升到 90% 幅度用的时间。这个定义与阈值无关只要峰值幅度标定可靠。声发射领域传统 AE 参数里的上升时间则是在波形越过触发阈值后到第一个峰值的时间它跟阈值设在哪直接挂钩。还有一个工程上常用的“包络上升时间”对波形做 Hilbert 包络再取 10-90%。三种定义各有适用场景定义计算路径受什么影响适合场景阈值→峰值触发样本到峰值样本阈值、噪声、触发比较器老采集系统、在线触发10-90% 法峰值归一化后找 0.1/0.9 交叉点峰值位置、噪声底离线处理、事件分类包络 10-90%Hilbert 包络后取交叉点包络平滑程度起点回溯、抗振荡传统 AE 系统全是硬件阈值触发加计数器所以“阈值到峰值”的定义沿用了几十年。但做声信号计算时我建议两个一起算阈值到峰值负责跟历史数据对齐10-90% 负责作为稳定特征来分类。只保留其中一个遇到阈值设置不合理时没有任何后悔药。3.2 一段可直接改用的上升时间测量代码下面的 Python 函数输入一段波形和采样率分别输出“阈值到峰值上升时间”和“10-90% 上升时间”。噪声底用波形前面一段无事件区间估计阈值取噪声 RMS 的 4 倍。import numpy as np def find_cross_time(sig, level, fs): idx np.argmax(sig level) if idx 0: return 0.0 frac (level - sig[idx - 1]) / (sig[idx] - sig[idx - 1]) return (idx - 1 frac) / fs def measure_rise_time(wave, fs, pre_s0.02, thr_factor4.0): wave np.asarray(wave, dtypenp.float64) # 前 pre_s 秒通常还没有事件用于估计噪声底 noise_seg wave[: int(pre_s * fs)] noise_rms np.std(noise_seg) 1e-12 thr thr_factor * noise_rms # 峰值位置和幅度取绝对值AE 信号极性不定 peak_idx int(np.argmax(np.abs(wave))) peak_amp wave[peak_idx] # 从峰值向前找第一个低于阈值的样本做线性插值 j peak_idx while j 0 and np.abs(wave[j]) thr: j - 1 t_cross (j (thr - np.abs(wave[j])) / (np.abs(wave[j 1]) - np.abs(wave[j]))) / fs rise_thr (peak_idx / fs - t_cross) * 1e6 # 微秒 # 10%-90%在峰值之前的上升段中找两个交叉点 pre np.abs(wave[: peak_idx 1]) a10 0.1 * np.abs(peak_amp) a90 0.9 * np.abs(peak_amp) t10 find_cross_time(pre, a10, fs) t90 find_cross_time(pre, a90, fs) rise_10_90 (t90 - t10) * 1e6 # 微秒 return rise_thr, rise_10_90, thr, peak_idx / fs说明几个关键点。find_cross_time返回信号首次越过指定电平的时间用相邻两个样本做线性插值避免采样点之间的量化误差。measure_rise_time里噪声底取波形前 20ms前提是这段确实没有 AE 事件前触发记录如果设得够长这 20ms 通常会落在事件之前。thr_factor取 4.0 是常用起点信噪比很高时可以降到 2.5环境噪声大时提到 6。阈值设太低起点会回溯到噪声底里上升时间整体偏大设太高又可能漏掉弱事件。如果手里有历史数据建议先用 10-90% 的结果画分布再拿阈值到峰值版本去和老系统对齐。这段代码有一个隐藏限制它假设波形里只有一个突发事件。如果上升阶段不止一个峰值比如两个源先后到达argmax找到的可能就是第二个峰上升时间会被测成两个峰之间的时间。碰上这种叠波需要在峰值后加一个时间窗约束比如只取峰前 200μs 范围并在记录前段检查是否已经有一个峰值。没有这个约束算出来的上升时间和做声发射定位时拿两个波混叠的到时去解方程一样结论全偏。3.3 起点回溯用短时能量和 AIC 代替“第一个过阈点”硬件触发器的原始时间戳是“信号越过阈值”的时刻这个时刻对上升沿内部的小振荡和噪声调制非常敏感。同一个断铅源重复打十次触发点会前后跳上升时间也跟着跳。经验做法是离线重新估计声信号真正的“起振点”。常用手段是短时平均/长时平均STA/LTA或 AIC赤池信息准则分割。思路可以浓缩成一段短代码。先对波形平方算短时能量再从峰值位置往回找“能量突然跌到噪声底”的那个拐点def onset_from_energy(wave, fs, win1e-4): n max(1, int(win * fs)) env np.convolve(wave**2, np.ones(n) / n, modesame) peak_idx int(np.argmax(env)) base np.mean(env[: n * 2]) k peak_idx while k 0 and env[k] base * 2.0: k - 1 return k / fs代码核心是当短时能量持续低于 2 倍噪声能量时认为信号还没到一旦能量开始往上爬就把这个位置作为起点。阈值系数 2 可以改成 3-5噪声越大系数越大。AIC 更精细在峰值附近一个滑动窗口里算每个分割点的 AIC 值最小值对应的样本就是最优起点但计算量大一些适合离线定位预处理。起点回溯在声发射定位里尤其关键。阈值触发的时间戳只适合在线事件记录不适合直接作为定位的时间基准。做定位之前我一般把所有通道的同一次事件拿出来统一用能量或 AIC 重新标一遍到时再去做时差计算。这样上升沿的随机抖动会被平滑掉定位的重复性能上一个台阶。时间基准先抖动几微秒后面的滤波、分类、定位全白干。4. 声发射定位时差估计为什么受上升沿形态支配TDOA 计算与最小二乘落地4.1 到达时刻误差和定位误差是线性关系声发射定位的基本原理是时差定位也叫 TDOA。已知 N 个传感器的坐标 (x_i, y_i)记录到同一事件到达各传感器的时刻 t_i假设波速 v 已知事件源坐标 (x0, y0) 和发震时刻 t0 满足(x_i - x0)² (y_i - y0)² v² (t_i - t0)²这个方程组是非线性的工程上常用最小二乘求解。这里最要紧的问题是t_i 从哪来在线采集系统给的 t_i 往往是“某通道越过阈值”的硬件时间戳而上升时间越长的信号阈值触发点在上升沿上的位置不确定性越大。原因很直观上升沿越缓同一电平附近的振幅变化率越小噪声一叠跨阈点的变动就是几微秒起步。换算一下波速按 3200m/s 算时差误差 10μs 对应定位误差约 3.2cm如果阈值抖动到 40μs定位误差就是 13cm。对一套传感器间距只有 20-30cm 的定位系统来说这种抖动足以把源位置定到传感器阵列外面去。所以在做定位前信号要先经过事件筛选只保留突发型事件上升时间过大或持续时间过长的信号不适合做 TDOA。经验上定位实验里能用的典型事件上升时间在 5-50μs 之间泄漏、摩擦这类连续信号没有明确的上升沿别拿来做定位。具体门限要结合传感器带宽来定不能跨系统照抄一个数值。提示上升时间筛选的门限不要拍脑袋先看你那套传感器的上升时间下限再乘 2-5 倍作为事件级过滤的下限。4.2 平面四传感器 TDOA 定位的最小二乘代码下面这段代码可以直接运行传感器坐标和到时需要按秒传入。代码用 scipy 的 least_squares 做非线性拟合。import numpy as np from scipy.optimize import least_squares def locate_tdoa(sensors, times, v3200.0): sensors np.asarray(sensors, dtypefloat) times np.asarray(times, dtypefloat) n len(times) # 初值传感器几何中心发震时间用最早到时往回推一点 x0 np.mean(sensors[:, 0]) y0 np.mean(sensors[:, 1]) t0 np.min(times) - 0.5e-3 init np.array([x0, y0, t0]) def residual(p): x, y, t p r [] for i in range(n): dist np.hypot(sensors[i, 0] - x, sensors[i, 1] - y) r.append(dist / v - (times[i] - t)) return np.array(r) res least_squares(residual, init, methodlm, max_nfev100) x, y, t0 res.x return (x, y, t0), res.cost调用时把四个传感器的坐标按米传入times 按秒传入。v 取 3200 是多数钢板上最低阶 Lamb 波的经验值严格做定位前应该用一块和被测对象同样材料的板在已知距离处敲击实测波速。初值设定很关键。t0如果给得离谱误差函数容易陷在局部极小。我一般用最早到达时刻减去“传感器间距/波速”的一半作为上界再把传感器几何中心作为位置初值基本都能收敛。另一个容易被忽略的点t0 是发震时刻的拟合结果不是我们要的定位状态。定位的绝对误差取决于波速标定相对误差取决于各通道到时的一致性。判断结果可不可信不能只看 x、y要看res.cost和每个通道的残差。如果某个通道残差总是特别大说明这个通道的到时标记有问题可能是传感器耦合不良也可能是这路信号上升形态和其他通道差太多。波速怎么标定不要只查手册直接在试件或现场结构上选一条已知距离的路径用断铅或脉冲激励实测一次到达时间差。标定要分方向和距离做几次结构里有焊缝、加强筋时波速并不均匀。频散强烈的板中群速度会随传播距离变化单一波速只是工程近似残差能说明误差量级就够了。4.3 通道一致性校准断铅源是最可靠的验证手段通道一致性怎么查在传感器阵列中间找一个已知点用 HB 铅芯在板面做断铅试验。断铅产生的声发射信号频带宽、重复性好是工程上公认的人工声源。每次断铅都当作事件记录下来重算到时统计各通道触发点之间的固定差值。如果某一通道每次都比另一通道晚 3μs那不是巧合是前放延迟或线缆长度不同。把固定偏差记下来后续定位时提前扣除。常见做法要求每根铅芯只断一次断过的部位不再重复在距离传感器最远的位置加做一次检查定位残差与距离的关系。如果远场误差明显变大先查波形幅值是不是过低导致到时抖动再考虑波速频散。板中 Lamb 波频散意味着不同频率成分传播速度不同上升时间越短的宽带信号到时定义越敏感。工程上一般用“主要能量到达时刻”作为到时而不是硬件触发时间能减小频散带来的整体偏移。声发射定位项目里前面采集和参数都正常定位就是不准十有八九是到时的定义和通道一致性没处理好而不是最小二乘算法本身有问题。5. 声发射信号处理避坑上升时间是干扰最容易“伪装”的参数5.1 假“快事件”ESD 放电与电磁干扰现象系统忽然记录一堆“上升时间只有几微秒甚至零点几微秒、峰值接近满量程、持续时间很短”的突发事件而且多通道几乎同时触发。把事件放进定位算法里算出来的源位置不在传感器阵列里甚至落在结构外面。原因这不是声是电磁干扰。ESD 放电的时间极短脉冲宽度在 1ns 以下到 200ns 之间上升时间从 200ps 以下到 10ns。这样的电磁脉冲耦合进屏蔽不良的传感器线缆后前置放大器看到的就是一个超快的电压冲击和压电晶体输出的声信号完全不像。电磁波比板中的弹性波快几十万倍所以多通道几乎同时触发定位自然发散。解决硬件上把线缆换成双层屏蔽同轴屏蔽层在整个链路里保持连续传感器和结构之间做好接触与等电位前放输入加限幅保护和带通滤波把带外高频能量直接滤掉。软件上给事件加两个快速过滤条件上升时间低于系统信号带宽允许的最小值就丢弃多通道到达时间差如果小到不可能由声速解释也丢弃。常见做法是设一个最小到达时间差例如 0.5μs只有大于它的通道对才进入定位计算。5.2 阈值抖动让上升时间忽大忽小现象同一个断铅源连续触发 10 次第 3 次显示上升时间 8μs第 7 次显示 16μs波形肉眼看起来几乎一样。原因触发阈值落在了信号前沿附近而声发射前沿不是干净斜坡是带有小振荡的包络。噪声每次不同跨阈点就会在振荡波峰之间前后跳动跳一个振荡周期上升时间就多出来一个周期。解决阈值不要贴着噪声底设至少高出噪声 RMS 3-5 倍在线触发只做事件发现不做精确定时。存完波形后在软件里用短时能量或 AIC 重新标起点。这是血泪经验硬件阈值时间戳除了筛选事件不要直接拿去做上升时间和定位的最终值。5.3 ADC 饱和导致“限幅上升时间”现象事件幅值很大波形峰值出现平顶此时算上升时间结果反而偏小10-90% 法的 90% 点从限幅值出发算出一个没有物理意义的结果。原因前置放大器增益太高或量程设置太小信号被 ADC 限幅。峰值被硬截断峰值位置从真实峰提前到限幅起始点波形尖部被削平后90% 电平的交叉点也错位。解决正式采集前先用已知源确认峰值不超过满量程的 90%。如果信号饱和降低前放增益再测一遍噪声底。事件峰值接近满量程时在参数表里标记为“饱和”上升时间直接放弃。保存原始波形的好处就在这种时候只要没被硬件完全毁掉还能重新提特征。5.4 定位残差大先别怀疑算法去检查时间戳现象定位解出来了残差有几十微秒坐标在传感器之间乱跳增加传感器数量后残差不降反升。原因各通道的到时标记本身不一致。有的通道信号上升沿陡触发点可靠有的通道信号弱、上升沿缓触发点明显偏晚还有前放延迟和线缆长度不同造成的固定偏差。这些偏差在 TDOA 方程里不表现为随机噪声而是作为系统性偏移进入所有时差残差自然下不去。解决先做断铅一致性校准把每通道的固定时延扣掉再把各通道到时统一到同一物理参考点例如峰值 50% 或能量拐点重新标记。如果仍然乱跳检查波速标定或者缩小通道间距重新布置。翻车案例里十次有八次是到时不齐剩下两次才是波速问题。没有原始波形时这类问题很难追查所以从一开始就要带着波形存。6. 把上升时间变成有效特征三层过滤、断铅验证和我的数据习惯做在线声发射监测时我不推荐把全部触发事件直接送给定位或分类器。我会在事件入库前加三层过滤。第一层是硬件带通排除工频、机械振动和大部分电磁干扰。第二层是事件级参数过滤上升时间低于系统可分辨下限的丢弃高于上限的标记为连续型事件到达时间差小于声速允许窗口的多通道组合直接打上“电磁串扰”标签。第三层是定位残差一致性检验定位后各通道残差超过提前设定的阈值事件标记为低可信用于统计但不出报告。验证这套流程我用断铅试验居多一支 0.5mm HB 铅芯在被测对象表面上距传感器阵列最近边 10cm 以上的位置压住铅芯弯断每次产生的声发射作为校准脉冲。至少做 5 个分散的已知源点记录定位结果与真实坐标的误差。源点编号真实坐标定位坐标距离误差平均上升时间P1(0.00, 0.00)(0.01, -0.02)2.2cm14.1μs误差在 10cm 以内基本可以认为是波速建模与到时刻标定在合理水平如果某个源点误差明显大于其他点回查那点到各传感器的信号幅度和上升时间多半是触发时间戳或信号太弱。最后说一个我的个人习惯只要采集系统允许原始波形永远和参数一起保存。上升时间、振铃计数、幅值这些参数只当索引不当原始证据。因为处理流程一旦发现阈值设错、增益饱和、传感器贴装不良还可以用波形重算所有参数如果只存参数表那只能重新测一遍。这个习惯帮过我太多次也希望帮到你。本文还有配套的精品资源点击获取