马世帮,李栋,解琪,李宏光,张灯,储隽伟,孙宇楠
(1.西安应用光学研究所,陕西 西安 710065;2.陆装航空局洛阳室,河南 洛阳 471000)
太赫兹时域光谱仪是火炸药、生化战剂等材料进行特征光谱测量的核心仪器,广泛应用于目标识别、危险品成分鉴别、穿透探测等军事应用领域[1-5]。信噪比参数和光谱范围参数是太赫兹时域光谱仪的关键技术指标参数[6-8],目前已具备太赫兹时域光谱仪特征光谱、线性度和透射比参数的计量能力。但由于光谱范围和信噪比参数存在着名称术语和指标定义不统一的问题,导致太赫兹名称术语常常被混用,从而引起仪器性能评价的误解。
太赫兹时域光谱仪光谱范围的定义不明确,导致同一台仪器使用不同方法进行光谱范围测量的结果相差很大。光谱范围的定义方式主要有两种[9],一是在频谱数据中,以超过N倍噪声限的频谱数据对应的上下限频率为光谱范围;二是在频谱数据中,以超过1/M频谱数据最大值的频谱数据对应的上下限频率为光谱范围,M取值为10~80。
信噪比是反映太赫兹光谱仪计量性能的关键指标[10-12],信噪比差的光谱仪很容易将噪声尖峰误判为特征光谱[13-14]。目前缺乏太赫兹信噪比测量和评价的统一方法,各种商用太赫兹光谱仪的测量结果无法准确统一。当前,太赫兹光谱仪存在几种不同的信噪比名词术语。如美国俄克拉荷马州立大学把主要信号段最大值定义为信号,信号前段峰峰值的 1/2 定义为噪声。首都师范大学把主要信号段峰峰值定义为信号,信号前段均方差定义为噪声[15]。不同术语之间的混用造成结果相差几个量级,严重影响太赫兹光谱仪性能判断。因此,有必要统一太赫兹光谱仪信噪比定义方法,使不同系统具有可比性。信号计算分为两种,分别是最大值和峰峰值,对于信号值的评估采用峰峰值更合理。
综上,需要研究太赫兹时域光谱仪光谱范围和信噪比校准技术,统一光谱范围和信噪比指标定义,解决太赫兹时域光谱仪光谱范围和信噪比校准难题,从而实现太赫兹时域光谱仪的量值统一。
太赫兹时域光谱仪按照不同的时间延迟方式,可分为机械延迟式时域光谱仪和异步采样式时域光谱仪。其中,机械延迟式时域光谱仪组成框图如图1 所示。
图1 机械延迟式时域光谱仪组成框图Fig.1 Composition block diagram of mechanical delayed time-domain spectrometer
机械延迟式时域光谱仪主要由飞秒激光器、分束器、时域延迟系统、太赫兹发射模块、准直镜1、会聚镜、样品室、准直镜2、带孔会聚镜、太赫兹探测模块、太赫兹弱信号处理系统和计算机组成。飞秒激光器的出射激光经分束器分为两束,一束作为泵浦光路,另一束作为探测光路。其中,泵浦光谱经时域延迟系统产生时间延迟之后入射至太赫兹发射模块,太赫兹发射模块经准直和汇聚后进入样品室之中,穿过样品,携带被测样品信息的太赫兹波经准直镜2 准直并经带孔反射镜反射之后,入射至太赫兹探测模块。另一束太赫兹探测光路穿过带孔反射镜之后,也入射至太赫兹探测模块。探测光路和携带样品信息的泵浦光路同时入射至太赫兹探测模块,太赫兹弱信号处理系统对太赫兹探测模块的输出信号进行放大滤波,并采集太赫兹弱信号处理系统随不同延迟时间下的信号曲线,即可得到太赫兹时域光谱仪的时域信号。经傅里叶变换之后,即可得到频谱信息。异步采样式时域光谱仪组成框图如图2 所示。
图2 异步采样式时域光谱仪组成框图Fig.2 Composition block diagram of asynchronous sampling time-domain spectrometer
不同于机械延迟式时域光谱仪,异步采样式时域光谱仪利用2 台重频不同的飞秒激光器替换了机械延迟式时域光谱仪的飞秒激光器、分束器和时域延迟系统。通过2 个激光器的重频之差,实现了泵浦光和探测光的时间延迟[15]。
对于定义方式一,由于噪声限不好确定,不同的噪声限确定方式对于光谱范围的影响很大,故排除该方法;对于定义方式二,频谱数据的最大值很容易确定,因此最终选择定义方式二作为光谱范围的测量依据。
太赫兹时域光谱仪光谱范围示意图如图3所示。以超过频谱数据最大值1/M的频谱数据作为光谱范围,其中M取10~80。
图3 太赫兹时域光谱仪光谱范围示意图Fig.3 Schematic diagram of spectral range for THz timedomain spectrometer
太赫兹时域光谱仪的信噪比分为时域信噪比和频域信噪比。时域信噪比可以反映时域光谱仪时域信号的信噪比水平,但不能反映每个频点下的信噪比水平;频域信噪比不仅可以反映不同频率下的信噪比水平,还可以反映不同频点的信噪比。
1)时域信噪比校准原理
太赫兹脉冲时域信号如图4 所示。其中a 段为信号前段,其波动由噪声引起;b 段为主要信号段,集中了THz 信号的绝大部分能量;c 段为信号衰减部分,THz 信号能量很少。由于频谱的可靠分辨率完全取决于c 段所取的长度,因此太赫兹时域脉冲信号的采集时间主要是在c 段。将图4 中的b 段峰峰值定义为信号,a 段均方差定义为噪声。由此可以得到太赫兹时域光谱时域信噪比计算公式:
图4 太赫兹时域光谱信号Fig.4 THz time-domain spectral signal
式中:Smax为时域脉冲信号的最大值,即b 段信号的最大值;Smin为时域脉冲信号的最小值,即b 段信号的最小值;为时域信号噪声的均方根值,即a 段信号的均方根值。
2)频域信噪比校准原理
不放入样品,利用时域光谱仪测量时域信号曲线,对采集的时域信号进行傅立叶变换,得到此时的信号功率谱曲线PSλ;将金属板移入样品室,利用时域光谱仪测量时域信号曲线,对采集的时域信号进行傅立叶变换,得到此时的噪声功率谱曲线PNλ;利用式(2)计算频域信噪比曲线RSNλ:
式中:RSNλ为频域信噪比;PSλ为信号功率谱曲线;PNλ为噪声功率谱曲线。
以日本爱德万公司的异步采样时域光谱仪作为被校对象,型号为TAS-7500TS。对空测时域信号进行处理,得到光谱范围。当光路内相对湿度为5%、且M取值不同时,光谱范围的校准结果如表1 所示。
表1 不同M 值的光谱范围校准结果Table 1 Calibration results of spectral range with different M values
由表1 可以看出,M值越大,光谱范围就越宽。当M值≥50时,光谱范围变化不大,因此建议M取值50。当M取值50时,不同湿度环境下的光谱范围校准结果如表2 所示。
表2 不同湿度环境下的光谱范围校准结果Table 2 Calibration results of spectral range with different humidity
由表2 可以看出,相对湿度越小,光谱范围就越宽。当相对湿度≤50%时,相对湿度对于光谱范围影响不大;当相对湿度>50%时,由于空气中水蒸气对于太赫兹的吸收急剧增加,所以光谱范围急剧缩小。因此在使用太赫兹时域光谱仪时,建议将光路内的相对湿度保持在50%及以下。
以日本爱德万公司的异步采样时域光谱仪作为被校对象,型号为TAS-7500TS。对空测时域信号进行处理,得到时域信噪比。不同湿度环境下的时域信噪比校准结果如表3 所示。
表3 不同湿度环境下的时域信噪比校准结果Table 3 Calibration results of signal-to-noise ratio in time domain with different humidity
由表3 可以看出,相对湿度越小,时域信噪比就越大。当相对湿度≤50%时,相对湿度对时域信噪比的影响不大;当相对湿度>50%时,由于空气中水蒸气对于太赫兹的吸收急剧增加,所以时域信噪比急剧缩小。因此在使用太赫兹时域光谱仪时,建议将光路内的相对湿度保持在50%及以下。
以日本爱德万公司的异步采样时域光谱仪作为被校对象,型号为TAS-7500TS。对空测时域信号进行处理,得到信号功率谱曲线。将金属板放入泵浦光路,对测得的时域信号进行处理,得到噪声功率谱曲线。当光路内的相对湿度为5%时,得到频域信噪比曲线如图5 所示。
图5 频域信噪比曲线Fig.5 Signal-to-noise ratio curve in frequency domain
由表5 可以看出,当频点从0.1 THz 变化至0.2 THz时,该时域光谱仪的频域信噪比迅速增加,之后增速变缓,待增加至1.4 THz 之后,该时域光谱仪的频域信噪比开始缓慢减小。另外,频域信噪比的最大值为70 dB 左右,和时域信噪比的最大值接近。
比较了目前常用的光谱范围和时域信噪比校准方法的优缺点,并优选其中一种方法,设计了太赫兹时域频域信噪比参数校准方法。研究了M取值对光谱范围的影响,M值越大,光谱范围就越宽。当M值≥50时,光谱范围变化不大,因此在光谱范围校准时建议M取值50。同时研究了光路内相对湿度变化对于光谱范围和时域信噪比的影响,相对湿度越小,时域信噪比越大,光谱范围就越宽。当相对湿度≤50%时,相对湿度对于光谱范围和时域信噪比的影响不大;当相对湿度>50%时,因为空气中水蒸气对于太赫兹的吸收急剧增加,光谱范围和时域信噪比急剧缩小。此外,完成了频域信噪比校准实验,得到了频域信噪比曲线,频域信噪比曲线的最大值和时域信噪比接近。在太赫兹时域光谱实验中,光路中的相对湿度必须≤50%,才能使时域光谱仪处在比较良好的工作状态之下。光路内相对湿度变化对光谱范围和时域信噪比的影响研究结论,对于太赫兹时域光谱仪的使用、太赫兹光谱仪相关标准的制定具有重要的指导意义。