搞了快十年的气体检测我一直在用催化燃烧传感器做可燃气体报警。直到有一次在天然气场站做维护一台刚标定不到一个月的催化传感器被硅烷蒸汽弄“哑”了——读数直接归零而现场实际上已经漏了半小时的甲烷。那次之后我才真正动了换技术路线的心思。激光可燃气体传感器核心就是把TDLAS可调谐半导体激光吸收光谱技术从实验室搬到了STM32这种几十块钱的MCU上用一束近红外激光代替传统催化珠不吃毒、不用氧气、响应还快。这篇文章从我自己的踩坑经历出发从TDLAS原理讲到STM32实战最后附上已经跑通的底层代码思路。如果你是做气体检测终端、地下管廊监控或者燃气智能硬件的照着这条线一步步来能省不少冤枉路。1. 从“误报到漏报”为什么我们最终选了激光方案1.1 催化传感器最无解的问题在这里大家用惯了催化燃烧催化珠传感器它便宜、原理简单遇到可燃气体时会在催化珠表面无焰燃烧引起电阻变化输出和浓度近似成正比。但它的三个硬伤在实际现场里很难绕开中毒不可逆含硅蒸汽、硫化物、卤化物一旦覆盖催化剂表面传感器灵敏度就会永久衰减只能换新必须有氧环境甲烷浓度如果超过LEL下限的某个范围燃烧机理失效读数反而会掉下来形成“高浓度还显示安全”的假象零点漂移长时间通电后桥路零点漂移明显每隔两三个月就得用标准气重新校准。在燃气场站和地下管廊场景里硅胶灌封剂、防腐漆挥发物都可能要了催化传感器的命。我之前那台设备就是吃了这个亏。1.2 激光传感器解决的是什么激光方案做的是物理光学测量根本不依赖被测气体的化学反应。一束波长可微调的近红外激光穿过被测气室特定分子的特征吸收峰会让光强发生衰减通过衰减程度倒推出浓度。激光器自身就是“选择性过滤器”甲烷只在1.653μm附近有强吸收线别的气体几乎不干扰所以不存在中毒这一说。用STM32做激光传感器的核心工作是把原本一套体积不小的光谱分析系统压缩成“激光管探测器MCU算法”四个模块。激光器温度通过TEC控制在0.01℃量级电流用DAC做高频正弦调制光电二极管出来的微弱信号再经过跨阻放大、数字锁相放大解调最后换算出浓度。整个过程不需要开机预热半小时上电几秒就能稳定出数而且本安性能极好因为激光功率通常在几毫瓦以下低于点燃阈值。1.3 这套方案适合谁如果你做的产品会长期暴露在恶劣工业环境里或者需要快速响应、频繁维护成本高的场景TDLAS是真香。比如燃气管网阀门井、地下综合管廊、沼气发电站、加气站这类地方传统传感器换得勤综合下来激光方案虽然单机物料贵一些但两三年不用换传感头省下的维护人力成本非常可观。如果只是做家用燃气报警器、对成本极度敏感那就先别上激光等供应链再成熟一些。2. 拆开TDLAS光、吸收峰、谐波检测到底怎么回事2.1 比尔朗伯定律是所有计算的地基TDLAS的一切都建立在比尔朗伯定律上[I_t I_0 \times e^{-A} I_0 \times e^{-\alpha \cdot L \cdot C}]其中(I_0) 是入射光强(I_t) 是透过被测气体后的光强(\alpha) 是气体吸收系数由分子吸收谱线决定(L) 是光与气体作用的有效光程(C) 是气体浓度也就是我们要算的目标。公式看着简单实际工程里最难的是把入射光强变化、温度对谱线强度的影响、光电探测器响应率的非线性统统剥离开。如果直接测(I_t)来算浓度激光管老化、镜片污染都会被误判成浓度漂移所以后面才会引入调制和谐波检测。2.2 为什么要选1.653μm这条甲烷吸收线甲烷分子在近红外有很多泛频振动吸收带最常用的两条路线是1.653μm和2.3μm附近的组合频带。2.3μm处吸收强度更高能测到更低的浓度但相应的激光器和探测器成本贵得多。1.653μm处在光通信C波段附近DFB激光器供应链特别成熟很多可调谐激光器模块都是给光通讯做的价格被压得很低所以消费级、工业中低端方案首选1.653μm。这条线上吸收系数虽然只有几十cm⁻¹·atm⁻¹但只要光程拉到20cm甚至更长检测下限做到几百ppm没压力对LEL甲烷爆炸下限5%体积比监测来说完全够用。2.3 直接吸收法和波长调制法的差别最直觉的做法是扫描激光波长扫过吸收峰测一条完整的吸收曲线然后扣基线、算峰高这叫直接吸收法。它的问题在于受到低频噪声、温度漂移、光源波动影响大现场做不好稳定。实际项目里我推荐用波长调制光谱WMS在激光工作点上叠加一个高频正弦波让激光波长在吸收峰附近来回摆动。由于吸收峰是非线性的输出信号会产生谐波分量。其中二次谐波2f信号的峰值与浓度近似成线性而且能巧妙避开1/f噪声和一阶线性干扰。配合STM32做数字锁相放大从含噪信号里把2f分量提出来比直接采样吸收曲线稳得多。2.4 二次谐波“卷积”直觉别被“锁相放大”四个字吓到。打个比方你在一场嘈杂的演唱会里想听清朋友喊你的名字你的耳朵不会去分析所有声音的频率而是专门等那个熟悉的音节出现。数字锁相放大就是做类似的事——用一个和激光调制同频的参考信号把探测器输出乘起来噪声被“平摊”到高频上只有同频的2f分量会变成缓变直流信号。这个东西在STM32上实现就是几行乘加运算后面会细说。3. 硬件选型清单从激光管到STM32的完整链路3.1 核心物料对比我实际搭过几版最稳的组合如下模块器件选型关键参数备注激光器1.653μm DFB蝶形激光器输出功率2~10mW线宽2MHz带TEC和内置热敏电阻激光驱动恒流源调制叠加偏置电流30~80mA调制频率5~20kHz必须低噪声纹波1μA光电探测器PIN光电二极管跨阻放大器带宽1MHz增益10^5~10^6 V/AInGaAs材料响应波段匹配MCUSTM32F407/F41112bit ADC、12bit DAC、高级定时器负责调制波形、锁相解算气室多次反射池或直通式光程20cm~1m现场用耐腐蚀铝合金即可DFB激光器选型时重点关注三个参数中心波长、边模抑制比、TEC控制精度。中心波长必须和甲烷目标吸收线对齐温漂一般控制在±0.3nm/K左右所以TEC必须做闭环否则吸收信号会直接漂出峰外。3.2 STM32型号为什么选F4系列TDLAS解调对实时性有一定要求但不至于上H7那么夸张。STM32F407主频168MHz带FPU做数字锁相和低通滤波够用。关键是它有两个12bit DAC和两个12bit ADCDAC用来生成叠加了正弦波的调制信号ADC同步采集光电探测器输出和参考光信号两路能同时采样省掉很多同步时序上的麻烦。如果你要做双激光器参比方案建议直接上STM32H743它有多路DAC/ADC和更大的RAM可以同时处理两个通道。但起步阶段F407足够逻辑更简单资料也多。3.3 激光驱动电路要注意的细节我第一版电路图用了普通运放做恒流源结果激光器电流噪声大光谱线宽被展宽吸收信号毛刺严重。后来换成“基准源R-R运放VGS控制MOS管”的经典压控恒流源并在激光器供电端并联一个100μF钽电容再加一个0.1μF/C0G电容电流纹波才压下来。另一个容易被忽略的是激光器防静电。蝶形激光器引脚内部有静电敏感结构尤其是靠近激光二极管那两根引脚焊接时必须用接地烙铁最好戴防静电手环。我做第一版的时候用普通烙铁直接焊废了两个管子烧了快一千块这个教训希望大家记住。3.4 探测器电路怎么搭InGaAs PIN光电二极管输出的电流信号极其微弱几百nA到几μA级所以跨阻放大器是第一级的关键。我习惯用AD8605或者OPA2376这类低偏置电流的精密运放跨阻电阻根据最大光电流和ADC量程来定。要注意的是跨阻放大器带宽不能太高太高会引入更多高频噪声也不能太低否则会把调制信号滤掉。一般设计带宽取调制频率的5~10倍比较合适。4. 光路与气室信噪比的第一道分水岭4.1 短光程直通式与多次反射式的取舍TDLAS里最小可测浓度和光程L直接挂钩。直通式气室光程短一般只有2~10cm适合微量级报警器但甲烷检测下限会比较差只能在0~5%VOL的高量程段工作。想要测到更低的ppm级泄漏就得靠多次反射气室让光线折返几十趟等效光程拉到0.5m甚至几米。自己做气室时我建议先别碰White池和Herriott池那种高难度结构调镜片能调到怀疑人生。最简单的是用两片镀金平面反射镜和一对柱面透镜做一个“Z”型几次反射结构。精度要求不高的话一个直通式气室加一个高增益探测器也能出效果只是检测下限能到0.1%VOL就很不错了。4.2 准直和耦合的实操经验激光器出来的光束发散角虽然不大但经过数厘米后光斑会散开。要在激光器前加一个非球面准直镜让光束变成平行光进入气室。出射端再聚焦到光电探测器上探测器感光面直径通常只有0.5mm左右光斑对不准信号就会掉一大截。我调试时习惯用一个可见光红激光笔先沿着光路走一遍把光路对准后再切换红外激光器。因为人眼看不见1.65μm红外光直接用红外光对光路非常痛苦。红激光和红外光同轴校准以后再微调探测器位置让ADC读数最大比盲调快得多。4.3 镜片污染是长期稳定性的宿敌工业现场气室里难免有灰尘、油污镜片一脏入射光强(I_0)直接衰减如果没有参比通道浓度读数会一路往上飘。我的做法是在光路里加一个分束片分出一路参考光不经过气室直接到另一个探测器。浓度计算时把信号光和参考光做比值镜片脏了只会同时影响两个通道比值基本不变。STM32有两个ADC通道做这件事很方便。5. 电路与信号链模拟端的每一步都在决定下限5.1 调制信号的产生DAC还是模拟电路传统TDLAS设备用DDS和模拟乘法器成本高。STM32里直接用DAC就能生成正弦调制信号再叠加一个直流偏置送到激光恒流源去控制激光器。DAC输出速率设为200kHz以上每周期采几十个点波形就已经足够平滑。但要注意DAC输出的阶梯波会引入高频谐波虽然这些谐波离2f信号很远但为了防止混叠还是要在DAC后加一个二阶低通滤波器截止频率设在调制频率的3倍左右。5.2 增益分配这样做才不会爆量程探测器输出经过跨阻放大后还要再经过一级可调增益放大才能进入ADC。增益分配的原则是保证最高浓度时信号不要削顶最低浓度时信号又不要太靠近ADC噪声底。我用STMF407内部的12bit ADC参考电压2.5V。按满量程信号不超过2.4V、本底噪声不低于20mV来算跨阻和二级放大总增益控制在10^5量级。现场调试时用两个电位器分粗调、细调比程序里做数字增益控制直观得多。5.3 噪声从哪来怎么压开关电源纹波给模拟运放供电的LDO要独立出来别让MCU和模拟电路共用一片LDO地环路跨阻放大器的输入输出尽量短模拟地和数字地在ADC处单点连接温度漂移激光器底座的TEC必须压稳电流采样电阻用低温漂的金属箔电阻温度系数选±5ppm/℃左右背向反射DFB激光器对光路背向反射敏感需要在激光器前加一个光隔离器否则会出现频率抖动。我当时排过一版PCB模拟部分和数字部分没做区分导致激光驱动电流抖动直接串进探测器调制信号上全是毛刺。后来重新布局把激光器驱动放在板子的最左端运放和探测器放在最右端MCU搁中间各走各的地噪声一下降了一个数量级。6. STM32固件怎么落从外设配置到数字锁相解调6.1 外设分配策略固件架构不复杂但外设分配要想清楚。我用的分配如下外设用途关键配置DAC1输出直流偏置正弦调制信号12bit输出缓冲开启ADC1采集探测器信号12bit3.3V参考1MspsADC2采集参考光信号与ADC1同步注入触发TIM2产生DAC触发、ADC触发、中断定时PWM模式输出频率和调制频率匹配UART1调试输出浓度值和原始信号115200bps调制频率我选了8kHzSTM32的定时器触发DAC输出正弦表每个正弦周期40个点因此DAC采样率是320kHz。ADC也由同一个定时器触发不过和DAC触发相位错开半个周期等探测器信号稳定后再采样这样能减少采样到建立过程中的毛刺。6.2 正弦调制的实现方法在STM32里生成调制正弦波最直接的做法是查正弦表。建一个40点的正弦数组DAC中断里每来一次就输出下一个点并把这个点同时叠加到直流偏置上。直流偏置对应激光器的中心工作电流正弦幅度对应波长扫描深度这两个参数在标定时调好存到Flash。为避免DAC输出缓冲带来的建立时间过长问题输出驱动电流建议小一点然后在DAC输出端接一个高输入阻抗的运放缓冲器。代码里基本就是查表、赋值、更新索引三步没有复杂计算。6.3 数字锁相放大器的代码思路这一节是子程序中最关键的部分。数字锁相的核心公式很简单[ I \frac{2}{N} \sum_{n0}^{N-1} x[n] \cos(4\pi f_{mod} t_n) ]其中(x[n])是ADC采集的探测器输出(f_{mod})是调制频率。因为我们要提取二次谐波所以参考信号频率是2倍调制频率。STM32每采集完一个完整正弦周期就做一次乘加运算输出一个2f幅度值。下面是我调试中实际用的一段简化代码注意这里只展示核心算法完整工程还需结合中断和外设初始化// 数字锁相放大提取二次谐波幅度 #define TAB_LEN 40 static float ref_sin[TAB_LEN], ref_cos[TAB_LEN]; static float adc_buf[TAB_LEN]; volatile float second_harmonic 0.0f; void lockIn_Init(float mod_freq, float sample_freq) { int i; float dw 2.0f * M_PI * 2.0f * mod_freq / sample_freq; // 2f归一化角频率 for (i 0; i TAB_LEN; i) { ref_sin[i] sinf(dw * i); ref_cos[i] cosf(dw * i); } } void lockIn_process(void) { int i; float sum_sin 0.0f, sum_cos 0.0f; for (i 0; i TAB_LEN; i) { sum_sin adc_buf[i] * ref_sin[i]; sum_cos adc_buf[i] * ref_cos[i]; } // 幅度取平方和开方也可以直接用幅值减少一次开方 second_harmonic 2.0f * (sqrtf(sum_sin * sum_sin sum_cos * sum_cos)) / TAB_LEN; }这里用积化和差玩了一个小技巧ADC信号分别和同相、正交的2f参考信号相乘得到I/Q两个分量幅度就是它们的矢量和。这样做的好处是对调制信号的绝对相位不敏感即使参考信号和实际相位移了90°幅度结果也基本不变省掉了相位对齐的麻烦。6.4 浓度换算与温度补偿算出2f幅度后还不能直接换算浓度。因为我们用的是比值测量把信号光通道的2f幅度除以参考光通道的直流光强得到一个归一化吸收信号。这个归一化信号和浓度之间也不是纯线性尤其是在浓度较高时由于吸收饱和曲线会微微弯曲。我的做法是标定时取5到7个浓度点多项式拟合浓度曲线浓度较高段做二次修正。温度影响分段处理甲烷吸收线强随温度有规律变化在-20℃到60℃范围内近似指数我实验室里在恒温箱里标了一组修正系数查表插值就行。压力影响更直接在压力稳定环境中影响不大但如果对准管道内部测量压力变化必须引入压力传感器做修正。7. 标定与实测比原理更麻烦的是这些坑7.1 纯氮零点为什么会飘新板子装好后先通纯氮理论上浓度应该是0但2f幅度并不完全为0。原因主要有三个残余光学串扰光源本身也有微量谐波、电路板上某级运放的偶次谐波失真、以及反射镜表面散射引起的本底。处理办法是软件里做一个“零气本底扣除”标定前记录纯氮下的2f幅度存到Flash里实际测量时扣掉这个本底。不过本底也有温漂最好在程序里每过一段时间自动执行一次零点校准或者通过气室电磁阀自动切入清洁空气。7.2 标气实操流程我标定甲烷时用1000ppm、5000ppm、1%VOL、2.5%VOL四瓶标气加上纯氮一共五个点。通气的顺序是从低浓度到高浓度每换一个浓度要让气室充分置换。直通式气室置换快几十秒就稳定多次反射气室内部容积大要等三五分钟。然后把每个浓度点下的2f幅度值记录下来。经验公式是[ C a_0 a_1 \cdot U_{2f} a_2 \cdot U_{2f}^2 ]用最小二乘拟合出系数(a_0, a_1, a_2)写入Flash。标定完后回测一遍误差如果超过满量程的1%就得检查气密性——我遇到过回测误差很大的情况最后发现是标气接头处的橡胶密封圈吸附了甲烷吹气半天都不干净。7.3 温度和压力交叉影响的修正一台传感器如果只在一个环境温度下标定拿出去分分钟翻车。我在恒温箱里做了-20℃、0℃、25℃、45℃、60℃五个温度点的浓度曲线发现同一浓度下的2f幅度在低温偏高、高温偏低。后来在程序里加了一个温度校正函数用NTC测气室温度查表做线性插值修正。压力影响比较复杂因为压力会同时改变吸收线型和气体密度。在井下管道场景压力变化大只补偿密度是不够的。我的建议是如果压力范围超过±5kPa就单独接一个气压传感器用理想气体状态方程近似做浓度补偿。7.4 镜片污染与斜率下降的识别现场运行一段时间后镜片污染会导致参考光通道幅值下降这时浓度读数虽然还能靠比值法稳住但如果污染非常严重信号通道的信噪比也会劣化。我在固件里加了一个自检逻辑如果参考通道光强低于标定时阈值的70%就上报“镜片脏污”警告提醒维护人员清洗气室。这个功能很实用别等到读数明显不准了再处理。8. 从能测到好用几个连我自己都想重做的细节8.1 双激光器参比方案是正道但起步别直接上单激光器结构简单但容易受到背景气体影响。比如气室里同时存在甲烷和水蒸气而水蒸气在1.653μm附近也有弱吸收单通道会把水蒸气的信号也当甲烷算进去。后续项目里我给其中一个通道加了一路波长为1.65μm附近但不落在甲烷吸收峰上的激光器做参比把背景吸收扣除掉效果立竿见影。起步阶段先用单激光器跑通全链路别一上来就上双激光器否则定位问题时分不清是光学问题还是电学问题痛苦加倍。8.2 低功耗改造的可能性激光传感器一般功耗比催化传感器高一些因为TEC一直在工作。如果你要做电池供电的无线检测终端可以这样省功耗让TEC在测量前预热1秒测量完成后立刻断电整个测量周期只持续20ms然后MCU进入停止模式。这样虽然平均功耗低了但每次测量的重复性会变差因为激光器没完全恒温。更好的方案是选一个低功耗小TEC配合高精度NTC让温控环在低电流下维持。实测下来我的样机平均功耗可以从1.2W降到0.3W代价是每分钟只能测一次。8.3 对外接口与现场部署要注意的事样机做出来以后别急着填漂亮外壳。我建议先预留几个常见的接口4-20mA工业电流环输出这是现场最通用的模拟信号RS485/Modbus方便接入PLC和SCADA一个干接点继电器报警输出用于联动电磁阀或声光报警器预留UART日志口调试时直接看波形和中间变量。安装时务必注意气室进排气口方向不能把气室完全密封也不能让风直接对着窗口吹。最好的做法是加一个防雨帽和一个微小的气体扩散膜既能保证气体交换又能防尘防直视光。8.4 最后一点个人体会从一块裸板到能稳定测出空气中微量甲烷整个链路里最花时间的不是代码也不是算式而是反复排查“信号为什么变小了”的过程。很多时候不是算法出了问题而是镜片有一层看不见的雾、激光器温度没稳住、或者某条地线多走了几毫米。做激光气体传感器一定要有一个心理准备前期的挫折感很强但一旦光路和电学链路稳了后面的路会越走越顺利。建议你从最小系统开始一步一步来把每一次波形变化的原因都记录下来这些记录比任何论文都有用。