简介面向声发射AE信号分析与材料状态监测人员提供一段MATLAB计算脚本用于从声发射数据中提取上升时间、幅度、持续时间、计数与能量等时域特征为裂纹扩展、塑性变形等内部过程的识别与结构健康评估提供数据基础。压缩包内仅含1个m文件整体约2KB结构轻量紧凑便于直接嵌入已有的AE数据分析流程。目前已有495人浏览/学习适合声发射入门学习者、无损检测从业者及材料损伤研究工程师参考。借助该脚本读者可以理解AE信号从读取、滤波、阈值设定到特征计算的一整套处理思路减少自行搭建代码的时间并能基于自身实验数据调整参数进一步分析材料动态响应与损伤程度为后续特征融合与损伤趋势分析提供支撑。1. 声发射信号里的上升时间为什么这个参数决定定位精度做声发射AE检测的工程师几乎都被波形里的上升时间折磨过。上升时间Rise Time是从声发射信号首次越过阈值到达到峰值幅值之间的时间差这个数看着不大却直接决定声发射源的定位误差和模式识别结果。标题里的 Acoustic_signal_computing本质上就是围绕上升时间这类瞬态参数做声信号计算的一套流程。本文从算法实现、参数整定到定位应用把上升时间在声发射信号处理里的角色讲透。不管你是刚接手声发射系统的新手还是想把定位精度再往上提一档的熟手照着这套思路做都能减少返工。声发射信号不是正弦波而是一个带有快速起跳、然后衰减振荡的瞬态脉冲。它的上升时间通常在几微秒到几十微秒之间取决于源的性质和传播路径。可就是这几十微秒在时差定位里一偏差位置就差出几毫米甚至几厘米。别小看这个参数后面每一步计算都跟它挂钩。2. 声发射信号的上升时间计算从原始波形到参数的四个步骤2.1 先定义上升时间起跳点和峰值点怎么找声发射信号处理领域对上升时间有两种常见定义。一种是“阈值穿越法”信号幅值第一次超过预设阈值到达到最大幅值的时间另一种是“前沿百分比法”从峰值幅值的10%上升到90%的时间。声发射检测仪器里内置的上升时间参数绝大多数用的是第一种因为它和硬件触发逻辑一致。而做声信号计算时我建议先用第一种算一个初值再结合波形前沿做人工校验因为阈值法受噪声影响很大容易把起跳点找偏。起跳点的本质是信号从背景噪声中“露头”的那个时刻。峰值点不能简单选全局最大值因为声发射波形到达峰值后会快速衰减但后续可能是反射波或振铃的次峰更高。正确做法是在起跳点之后开一个固定窗口在窗口内找最大值这个窗口宽度要覆盖一个主波峰但窗口晚了就会把二次波也算进去。2.2 用 Python 实现一个最小可用的上升时间提取脚本我习惯先把波形从数据采集卡导出成 CSV然后用 Python 做计算。下面这个函数是我在处理声发射数据时常用的一段核心代码它把阈值穿越和峰值查找放在一起输入波形和时间轴输出上升时间和对应的索引。import numpy as np def calc_rise_time(wave, time, threshold0.15, peak_win200): 计算声发射信号的上升时间 wave : 一维幅值数组单位 V 或 mV time : 时间数组单位 s threshold : 触发阈值比例相对于波形最大幅值 peak_win : 峰值搜索窗口长度单位“采样点” # 去除直流偏置取前200个点的均值作为基线 baseline np.mean(wave[:200]) wave wave - baseline # 用绝对幅值最大值作为参考 amp_max np.max(np.abs(wave)) if amp_max 0: return 0.0, 0, 0 # 找第一次越过阈值的索引 thresh_val threshold * amp_max cross_indices np.where(wave thresh_val)[0] if len(cross_indices) 0: return 0.0, 0, 0 trigger_idx cross_indices[0] # 在 trigger_idx 后的一段窗口内找峰值 seg_start trigger_idx seg_end min(trigger_idx peak_win, len(wave)) if seg_end seg_start: return 0.0, trigger_idx, trigger_idx peak_offset np.argmax(wave[seg_start:seg_end]) peak_idx seg_start peak_offset # 上升时间 峰值时刻 - 阈值穿越时刻 rise_time time[peak_idx] - time[trigger_idx] return rise_time, trigger_idx, peak_idx这段代码的逻辑分三步。第一步去直流因为采集卡输出常有几十毫伏的偏置不去干净的话阈值穿越点会提前或延后。第二步用最大幅值的比例做阈值而不是固定绝对电压这样对不同衰减程度的信号都能自适应。第三步是峰值搜索这里我限制了一个peak_win窗口防止把后续反射波当主峰。参数方面threshold0.15是经验值信噪比高时可以放到 0.1信噪比低时我会抬到 0.25避免噪声尖峰触发。peak_win200需要根据采样率调整如果采样率是 10 MHz200 个点对应 20 微秒通常够一个完整主波如果是 1 MHz200 点就太长要改成 50 左右。所以这个脚本换采样率时必须同步改参数这是我最早翻车的地方。2.3 采样率与分辨率为什么上升时间计算要卡在微秒级声发射信号的上升时间本身可能就是几个微秒如果采样率只有 500 kS/s一个采样周期就是 2 微秒整个上升时间也就几十个点稍微有点噪声就测不准。我实际操作下来计算上升时间的采样率至少要 2 MHz最好 5 MHz 以上。这里不是玄学而是数字信号的分辨率问题峰值时刻的误差可以控制在半个采样周期内但阈值穿越点的误差受波形斜率和噪声影响更大。给一张我常用参数表方便新手在做 Acoustic_signal_computing 时直接参考。这里的前提是采集卡是 12 位以上、输入带宽够高。应用场景采样率建议阈值比例建议峰值窗口建议金属试件断铅信号5 MS/s0.15100 点复合材料损伤监测10 MS/s0.2150 点高速撞击/ESD 脉冲对比50 MS/s 以上0.150 点低频传感器谐振型2 MS/s0.2580 点注意峰值窗口的单位是点不是时间所以换采样率时要换算。比如 10 MS/s 下 100 点是 10 微秒而 2 MS/s 下 100 点是 50 微秒差别很大。我做数据处理的时候会先把采样率写进配置文件所有窗口参数都以时间为单位再在代码里换算成点数这样不容易踩坑。3. 上升时间与声发射定位时差法里怎样用它校正波达时刻3.1 声发射定位为什么对上升时间敏感声发射定位最常用的方法是时差法TDOA布置多个传感器测量同一个源信号到达不同传感器的时刻差再按波速反算源位置。关键就在“到达时刻”怎么定义。多数采集硬件是把信号超过硬件触发电平的那一瞬间记为到达时刻。可是波形上升时间越长超过固定电平的时刻相对真实源起跳时刻就越晚也就是说到达时刻里包含了一段由上升时间决定的时间偏移。更麻烦的是同一源信号传播到不同距离的传感器时幅值不同、上升时间也不同。近处传感器信号幅值大、上升快触发时刻接近真实到达远处传感器信号幅值小、上升缓触发时刻滞后明显。这样即使波速和传感器位置都精确定位结果也会系统性偏向近处传感器。声发射定位的热搜词里经常有人问“为什么同一组数据不同通道算出来的定位差很多”十有八九就是没做上升时间校正。3.2 用上升时间对到达时刻做补偿一个简单校正模型我常用的校正思路是将触发时刻往回退一段退的量等于上升时间的一半。原理是假设信号起跳后近似线性上升到阈值那么真实起跳点大约在阈值穿越时刻之前半个上升时间的位置。这个模型虽然粗糙但工程上很实用尤其是对于谐振型传感器主频成分起主导作用时这个线性近似误差不大。校正后的到达时刻可以写成def corrected_arrival_time(trigger_time, rise_time, correction_ratio0.5): 用上升时间校正到达时刻 trigger_time : 硬件触发的绝对时间单位 s rise_time : 该通道计算的上升时间单位 s correction_ratio : 回退比例默认 0.5 # 回退半个上升时间使触发时刻更接近真实起跳 t0 trigger_time - correction_ratio * rise_time # 对所有通道校准时保持 correction_ratio 一致 return t0这个校正的物理意义很直白阈值点不一定是波形起点它只是幅值走到阈值那一步的时间。把触发时刻回退相当于把每个传感器的“主观触发时刻”统一对齐到波形起跳点附近。correction_ratio一般在 0.3 到 0.7 之间具体取决于传感器类型和滤波设置。谐振型传感器我会用 0.5宽带传感器用 0.3因为宽带传感器的前沿更陡。做完这步再做 TDOA 定位结果通常会有明显改善。我曾经处理过一组碳纤维板的断铅信号未校正前定位散点半径约 15 mm校正后缩小到 6 mm 以内。但要注意如果各通道的上升时间都是人工读出来的那误差也会被放大所以校正的前提是上升时间算得准。3.3 声发射定位中上升时间的经验阈值设置在声发射检测标准里上升时间没有统一的门限它是一个随传感器类型和试件材料变化的量。我做钢结构、铝合金和复合材料时会分别积累自己的数据库。比如铝合金中模拟源断铅的上升时间多在 520 微秒复合材料中由于衰减大远端信号上升时间可能超过 50 微秒。如果定位时遇到上升时间异常大的通道我不会盲目校正而是先排除这个通道——它可能处于信号盲区或传感器耦合不良。实际定位时我给每个通道设置一个上升时间的合理范围超出范围的通道直接剔除。范围怎么定先采集 20 次已知位置的断铅信号统计每个通道上升时间的均值和标准差然后以均值±3倍标准差作为有效范围。这样能有效避免把松动传感器或噪声误触发的通道算进定位方程组里。这也是声发射信号计算里“先清洗再计算”的原则。4. 声发射信号计算中的常见坑上升时间被高估或低估的排查记录4.1 坑一阈值设太低把噪声当成起跳点现象计算出的上升时间只有几百纳秒明显小于同类信号。看波形图起跳点其实在噪声尖峰上。我遇到过一组数据基线噪声有 20 mV信号幅值 200 mV阈值设为最大幅值的 5%结果噪声尖峰先越过阈值把起跳点提前了。原因阈值比例取太低了没有避开噪声的尖峰幅度。解决先看原始波形的噪声峰峰值把阈值设置为噪声最大值的 23 倍同时用信号最大幅值的比例计算。比如噪声 20 mV、信号最大幅值 200 mV那阈值至少设为 40 mV对应比例 0.2。我会写一个预处理函数先自动估算噪声标准差再动态确定阈值下限。4.2 坑二滤波器把上升沿磨平了高估了上升时间现象同一信号滤波前上升时间 8 微秒滤波后变成 20 微秒。波形的前沿明显变缓峰值位置也往后移了。原因很多数字滤波器比如低通 FIR 滤波器会引入群延迟尤其对高频成分丰富的声发射前沿延迟更明显。解决滤波器的主要目的是去掉低频振荡和高频电磁干扰但不要在那里用截止频率过低的低通滤波器。我一般用带通滤波下限 1 kHz 去掉直流和低频晃动上限 1 MHz 以上保留前沿。如果要做上升时间分析最好在滤波前的原始波形上算否则就要记录滤波器的群延迟并把时间轴往回移。这个坑不亲眼见过很难信滤波后的波形确实更干净但参数已经变了。4.3 坑三峰值点找错把振铃的次峰当成了主峰现象上升时间算出来异常大比如 200 微秒而实际主峰只有 30 微秒。波形上能看到真正的第一个峰值后面跟着更大的反射波。原因峰值搜索窗口开得太宽搜索到后续波包。解决把窗口宽度缩短。窗口宽度应该覆盖“一个主波”的长度而不是整串波列。对于谐振型传感器主波长度通常是一个振荡周期的 35 倍比如 150 kHz 谐振频率振荡周期约 6.7 微秒窗口取 2030 微秒就够。我用 Acoustic_signal_computing 时会把峰值窗口参数单独列出来每次处理前先抽一个通道波形人工看一下再定窗口宽度。4.4 坑四采样率不足上升时间只有三四个采样点现象上升时间计算结果跳变很大同样条件下两次信号一个 3 微秒一个 8 微秒重复性差。原因采样率 500 kS/s每 2 微秒一个点上升沿总共才 3 个点不同的相位差就会造成很大的时间测量误差。解决给采集卡设置足够高的采样率做声发射时至少 2 MS/s最好 5 MS/s 以上。如果硬件已经固定那就只能减小阈值比例让穿越点更靠近前沿陡峭区同时多次平均。另外在计算时可以对波形做插值比如用三次样条插值把时间分辨率提高到 0.2 微秒但不能把插值当作真实信息消除误差。这个坑属于“硬件不够软件来凑”但凑出来的结果要有不确定性标注。4.5 坑五不同通道触发电平不一致导致定位偏差现象四个通道的上升时间测出来都正常但定位结果偏差很大。检查发现各通道的硬件触发阈值分别设成了 30 mV、40 mV、50 mV、35 mV。原因多通道声发射采集系统里每个通道的触发电平是独立可调的如果校准不到位同一个信号在不同通道会有不同的触发延迟。解决把所有通道的触发电平统一设置为同一个电压值并在采集系统自检时用同一个信号发生器输出脉冲逐一检查各通道的触发时刻偏差。通常要求通道间触发时刻偏差小于一个采样周期。这个坑在工程现场很容易出现因为安装人员会为了减少误触发单独调低某个通道结果整个定位系统的时间基准就乱了。我也写过校验脚本在布置传感器后先做一次模拟源激励把所有通道的触发延迟偏差存下来后续定位时自动补偿。5. 用 Acoustic_signal_computing 批处理声发射文件的脚本设计与参数调优5.1 脚本结构从单文件到批量文件夹实际工程中不会只处理一个波形文件而是一测就是几百个文件每个文件包含多通道信号。Acoustic_signal_computing 的思路就是把单文件计算函数封装成批处理管道。下面是一个典型的批处理脚本骨架它扫描输入目录下所有 CSV 文件对每个文件的第一个通道计算上升时间并把结果汇总到 DataFrame 里。import pandas as pd import numpy as np from pathlib import Path # 已有的单通道上升时间计算函数 def calc_rise_time(wave, time, threshold0.15, peak_win200): # 省略具体实现见第2章 pass def process_file(file_path, channel_cols, time_col, threshold0.15, peak_win200): 处理单个CSV文件返回该文件各通道的上升时间 channel_cols : 通道幅值列名列表如 [ch1,ch2,ch3,ch4] time_col : 时间列名 df pd.read_csv(file_path) t df[time_col].values results {file: Path(file_path).name} for ch in channel_cols: wave df[ch].values rt, idx_trigger, idx_peak calc_rise_time( wave, t, thresholdthreshold, peak_winpeak_win ) results[frise_time_{ch}] rt * 1e6 # 转为微秒 return results # 批量处理 input_dir Path(./ae_waveforms) output_rows [] for csv_file in sorted(input_dir.glob(*.csv)): row process_file( csv_file, channel_cols[ch1, ch2, ch3, ch4], time_coltime, threshold0.15, peak_win200 ) output_rows.append(row) summary pd.DataFrame(output_rows) summary.to_csv(rise_time_summary.csv, indexFalse) print(summary.head())这段脚本的核心是process_file。它把读数据、算上升时间、转单位三件事放在一个函数里方便后续替换滤波器或加诊断输出。参数channel_cols和time_col需要根据实际导出的 CSV 列名修改。我习惯在calc_rise_time内部同时返回触发索引和峰值索引这样批处理时可以把索引存下来用于人工回溯。批处理最大的价值是能快速发现异常通道。比如summary里某个通道的上升时间一直比其他通道大一个数量级那就要回去看那个传感器的耦合状态。调优时不要只改阈值还要配合通道间一致性分析。5.2 必调的三个参数阈值、滤波上限、峰值窗口批处理脚本里的参数不能靠拍脑袋。第一个是阈值比例它决定起跳点位置。信噪比高的数据用 0.1信噪比低用 0.2 以上。这个参数会直接影响上升时间绝对值所有通道必须用同一个值不然定位校正没意义。第二个是滤波上限。声发射信号的前沿对应高频成分如果采集系统本身有抗混叠滤波器上限会固定但在软件处理时我会用一个高阶低通滤波器把超过传感器谐振频段的高频噪声滤掉。这个上限不要设置得太低。比如传感器谐振频率 150 kHz我会把软件滤波上限放到 300 kHz只滤掉明显超出传感器响应范围的成分。如果上限设到 150 kHz等于把前边缘削掉上升时间增大。第三个是峰值窗口。前面说过窗口以时间计更稳。我一般默认 30 微秒在 5 MS/s 采样率下就是 150 个点。如果信号波形主波较长比如复合材料低频振铃我会把窗口加大到 80 微秒。这个参数调优的方法很简单取 10 个典型波形分别用不同窗口宽度计算上升时间找到结果稳定且不与后续波包重叠的窗口区间。5.3 怎么验证计算结果用波形标定和人工抽查算完批处理的上升时间不能直接拿去定位必须做验证。我常用的验证方法有两种。第一种是“波形标定”把计算出的触发索引和峰值索引画到原始波形上看这两个点是否分别位于波形起跳处和真实峰值处。我会在脚本里输出一个 PNG 缩略图每 50 个文件抽一个肉眼扫一遍。第二种是“重复性验证”用信号发生器输出固定上升时间的脉冲经过功率放大器激励声发射传感器然后采集计算。如果计算出的上升时间和信号发生器设定的上升时间差在 10% 以内说明算法和参数是可靠的。这个验证过程也是 Acoustic_signal_computing 里最耗时的地方但省不了。我统计过批处理计算一把过的情况很少总有几个文件因为噪声异常导致上升时间离谱。人工抽查的“后悔药”就是保留触发索引和峰值索引发现问题后能快速定位到具体波形而不是重新读数据。这也是我在上面脚本里特意让calc_rise_time返回索引的原因。6. 把上升时间测准的一个技巧用 ESD 脉冲做系统标定最后分享一个我常用的校验技巧。静电放电ESD的放电时间极短脉冲宽度在 1 ns 以下到 200 ns 之间上升时间从 200 ps 以下到 10 ns正好可以作为极快前沿的激励源。做声发射系统标定时我会用一个小型 ESD 发生器在传感器表面附近放电让采集系统记录这个近乎理想的阶跃冲击响应。为什么用 ESD 而不是断铅因为断铅的上升时间本身有离散性而 ESD 的上升时间可以稳定到纳秒级能把系统的带宽和触发延迟暴露出来。具体做法把传感器贴在金属板上ESD 发生器的放电尖端距离传感器 35 mm触发放电后采集波形。分析这个波形时计算出的上升时间不再是声发射源的上升时间而是“传感器 采集卡 滤波器”这个链路对快前沿的响应时间。如果这个响应时间大于 5 微秒说明传感器谐振特性太强用来计算真实声发射源上升时间时必须做反卷积修正否则所有上升时间都是传感器的固有响应不是源信息。这个技巧特别适合评估一套声发射系统能不能用于上升时间敏感的分析。我遇到过一套便宜采集卡标称采样率 5 MS/s但用 ESD 脉冲实测发现上升时间被拉长到 15 微秒后来检查发现是前置模拟滤波器带宽只有 30 kHz。没有这一步标定我还在傻傻地测断铅的“上升时间”其实量的全是滤波器。现在每换一套硬件我都会先打一道 ESD 脉冲花 10 分钟量链路上升时间心里就有底了。希望这个习惯对你也管用。本文还有配套的精品资源点击获取