4. 穿透曲线数据预处理

从吸附床在线分析仪或离线取样得到的原始穿透数据,几乎不可能直接拿来画曲线、算传质区长度或拟合模型。原始数据里通常混着三类“杂质”:一是高频噪声,来自检测器电子漂移、泵脉动或取样间隔过密;二是基线偏移,来自零点漂移、载气本底或检测器预热不足;三是时间与单位的不一致,比如进气浓度用 ppm、出口浓度用 mg/m³,或者不同通道的采样时刻错位。数据预处理的任务,就是把这些“杂质”清理掉,让后续的穿透点、传质区、吸附容量计算建立在可信的数据基础上。本节按“去噪—基线校正—时间轴对齐—单位统一—格式标准化”的顺序展开,并说明每一步的判断依据。

一、去噪:先分清噪声与真实信号

去噪不是越平滑越好。穿透曲线的关键信息集中在穿透初期浓度开始抬升的“拐点”和接近进气浓度时的“平台”,过度平滑会把拐点抹圆,导致穿透时间偏晚、传质区被高估。工程上常用的做法是先看原始数据的采样间隔和噪声幅度:若采样间隔为 1 min、噪声峰峰值约为满量程的 1%~2%,可采用 5 点或 7 点滑动平均;若噪声更大或存在明显尖峰(如阀门切换引起的瞬时扰动),应先用 3σ 准则或中值滤波剔除野值,再做滑动平均。

以一组示例数据说明:某活性炭床出口浓度在 0~40 min 内基本在 2~6 ppm 波动,40 min 后开始单调上升。若直接对全段做 11 点滑动平均,40~50 min 的拐点会被拉平约 3~4 min。更合理的做法是分段处理:平台段可用较大窗口(如 9 点)压制噪声,拐点段用较小窗口(如 3 点)保留变化趋势。判断去噪是否合格,可看两点:去噪后曲线仍保持单调上升段的基本形状;残差(原始减平滑)中没有明显的系统性偏差,只有零均值随机波动。

二、基线校正:把“零点”拉回真实位置

基线校正针对的是系统性的零点偏移。常见情形是:吸附床尚未穿透时,出口检测器读数并不为零,而是稳定在某个本底值 Cb(例如 3 ppm),这可能是载气携带的微量杂质或检测器零点漂移。若不校正,穿透点会被误判为“浓度超过某个绝对阈值”的时刻,而不是“浓度相对本底开始显著上升”的时刻。

校正方法取决于基线是否稳定。若穿透前有一段足够长的稳定平台(示例中 0~30 min 读数在 2~4 ppm),可取该段中位数作为 Cb,对所有时刻做 C′(t)=C(t)−Cb,并规定负值截断为 0。若基线存在缓慢漂移(如从 2 ppm 漂到 5 ppm),则应采用线性或分段线性基线拟合,再逐点扣除。校正后应检查:穿透前平台是否贴近零线;穿透后平台是否接近“进气浓度减本底”的预期值。若两者都对不上,说明基线模型选错了,需要回到原始数据重新判断。

穿透曲线预处理四步链路与关键判据 ① 去噪 野值剔除 + 滑动平均 平台段大窗口 / 拐点段小窗口 判据:拐点不被抹平 残差零均值 ② 基线校正 C′(t)=C(t)−C_b 稳定平台取中位数 / 漂移用线性拟合 判据:平台贴零线 负值截断为 0 ③ 时间轴对齐 统一 t=0 起点 线性插值到统一采样网格 判据:进气/出口时刻一致 无重复时间戳 ④ 单位统一 ppm ↔ mg/m³ 格式标准化 判据:C/C₀ 无量纲 字段齐全可追溯 预处理前后对比(示例数据,非实测) 时间 t / min 浓度 原始(含噪、基线偏移) 预处理后(去噪+基线校正) 基线 C_b 已扣除 单位换算示例 C[mg/m³] = C[ppm] × M / 22.4 × (273/(273+T)) 示例:甲苯 M=92,25 ℃,100 ppm ≈ 100×92/24.45 ≈ 376 mg/m³ 统一后以 C/C₀ 无量纲量作图

三、时间轴对齐:让进气与出口在同一时钟上

穿透曲线本质上是“出口浓度随时间变化”,但这个“时间”必须与进气浓度的时间基准一致。实际中常见的问题是:进气浓度由钢瓶标气或发生装置给出,记录时刻为 tin;出口由在线分析仪按自己的采样周期记录,时刻为 tout,两者起点可能相差几十秒到几分钟。若直接把两条序列按行号对齐,穿透时间会产生系统性偏差。

处理方法是:先确定一个公共零点,通常取“吸附床开始通气的时刻”为 t=0;再将各通道的原始时间戳减去各自零点,得到相对时间;最后在统一的时间网格(如每 1 min 一个点)上做线性插值。插值时注意两点:一是不要外推,网格范围应取各通道时间范围的交集;二是若某通道存在重复时间戳(如仪器重发),应先按时间排序并取平均或保留最后一条,避免插值出错。对齐后可用一个简单检查:穿透前出口浓度应平稳,穿透后上升段与进气浓度的阶跃时刻在图上应能对应上。

四、浓度单位统一与数据格式标准化

单位不统一是穿透曲线分析中最容易被忽视、后果却最直接的问题。ppm 是体积分数,mg/m³ 是质量浓度,两者换算依赖分子量和温度、压力。常用换算式为:C[mg/m³]=C[ppm]×M/22.4×(273/(273+T)),其中 M 为吸附质摩尔质量(g/mol),T 为气体温度(℃)。例如甲苯 M=92,25 ℃ 时 100 ppm 约合 376 mg/m³。若进气用 ppm、出口用 mg/m³ 而不换算,C/C₀ 会完全失真。

格式标准化则要求每条记录至少包含:相对时间 t(min)、出口浓度 C(统一单位)、进气浓度 C₀(统一单位)、计算得到的 C/C₀,以及数据来源标记(通道号、采样方式)。建议统一保留 3~4 位有效数字,时间保留 1 位小数,避免因舍入造成平台段抖动。标准化后的数据表应能直接用于绘图和后续计算,且每一列的含义、单位、零点定义都可追溯。

预处理顺序建议固定为“先去噪、再基线校正、后时间对齐、最后单位统一”。若先做单位换算再做基线校正,本底值也会被换算,容易在换算系数上引入额外误差;若先对齐再基线校正,插值会改变本底段的统计特性,使中位数估计不稳。
去噪窗口、基线模型和插值网格一旦确定,应在同一批数据中保持一致,并在报告中记录。不同样品用不同窗口平滑后再比较穿透时间,结论不可比。

完成上述四步后,建议输出一份预处理记录,至少包含:原始文件与通道、采样间隔、去噪方法与窗口、基线取值或拟合方式、公共零点定义、插值网格、单位换算系数及温度假设、最终数据行数。这份记录既是数据可追溯性的要求,也是后续复现和排障的依据。预处理做得扎实,后面的穿透点判读和传质区计算才有意义。