全文HTML
--> --> -->随着社会的进步与发展, 光声光谱仪的需求与日俱增. 目前市面上光声光谱仪用光源大部分直接使用激光器, 少部分使用非相干扩展光源配合其他光学元件进行简单的聚光. 由于不同物质红外吸收峰的差异, 前者不能满足同时检测多种待测物质的要求[5], 而后者多未对光源系统进行过多的杂散光限制, 只在后期处理中进行滤波. 相比而言, 覆盖相同的光谱范围, 非相干扩展光源结构更简单、体积更小、价格更低, 而其产生的杂散光也可以在光学设计阶段通过光场整形来进行消除. 三维非相干扩展光源中, 红外热辐射光源因其输出功率大、光谱范围宽等优点, 是商用多组分、高精度光声光谱仪的首选光源[6].
给定光源和目标光场分布, 逆向求解光学系统结构及光学元件表面参数是非成像光学系统设计的一般过程. 在光声光谱仪用三维扩展光源光场整形系统中, 光源和目标光场两端的要求都极为苛刻. 对于输入端光源来说, 三维非相干扩展光源方向性极差, 向全空间辐射, 能量密度很低, 光学系统需要收集足够多的光能, 供系统进行光场整形; 对于输出端目标光场来说, 要求光线在光声池腔体内形成一个毫米级半径、厘米级长度的柱状分布; 而对于光场整形系统本身, 由于光声光谱仪仪器内安装尺寸的限制, 体积不宜过大, 不能使用过复杂的光路结构和过多的光学元件.
虽然这种极端的设计要求相对较少, 但国内外也有一些设计案例分别针对输入和输出两端的情况进行讨论. 1996年, 美国科学应用国际公司的Ong等[7]完成了一个用于管状光源照明的边缘光线系统的设计, 此系统专为管状光源的边缘光线设计, 在保证最大辐射效率的同时, 改善光通量的分布情况. 1997年, Shatz等[8]通过非旋转对称的形式完成了针对球形和圆柱形非相干光源的光场整形和传输效率的优化. 2007年, Henning[9]发现非相干耦合系统中若光源位于椭圆反射镜焦点位置不能实现最大耦合效率, 文章通过计算椭圆反射镜的高倍放大特性解释了最大耦合效率位置为何偏离焦点. 2009年, Florian等[10]提出了面向任意光强分布的扩展光源使用单反射镜产生任意形状非均匀光强分布的方法, 并给出了相关设计实例. 2014年, Wester等[11]针对LED (light-emitting diode)光源在光学设计中视为点光源而需要加大光学元件尺寸的问题, 提出了两种基于自由曲面的扩展光源设计方法. 2015年, Wu等[12]提出了一种有效的为扩展光源创建指定强度分布的设计方法, 使边缘杂散光线得到控制. 2020年, Hoffman和llan[13]提出了针对面光源或管状光源的分段反射聚光器优化方法.
本文通过前期实验实际测量光源发光特性, 并根据测量结果在物理尺寸和辐射方式上均进行了扩展光源精准建模. 在光声光谱仪用三维扩展光源光场整形系统光学设计的过程中引入扩展光源概念, 对三维非相干扩展光源辐射的光场进行整形, 致力于从光学角度、整个仪器系统的输入端提升光声光谱仪的性能参数.
2.1.红外热辐射光源
红外热辐射光源工作温度较高, 自身发热量较大. 商用红外热辐射光源发光体形状不规则, 多为柱状或螺旋管状, 严重偏离点光源—非成像光学设计的理想情况, 并且有着非典型的辐射通量空间分布, 在光学仿真软件ZEMAX中使用几何形状类似的管光源、圆柱体光源等非导入光源, 模拟辐射分布与实际光源偏差严重. 因其复杂性, 难以单纯的使用以点光源为默认测量条件的近场光源文件进行相关数据获取和模型化描述, 而红外热辐射光源制造商一般不提供说明发光特性的配光曲线等光源文件.以氧化硅、氮化硅、碳化硅等硅的化合物为发光材料的红外热辐射光源功率较大, 是光声光谱仪光源系统理想的辐射源. 以美国Hawkeye公司IR-Si272型红外热辐射光源为例, 其辐射体为柱状氮化硅, 尺寸为2.8(D) mm × 5(L) mm, 长度与直径比较小[14], 较为接近点光源, 由6 V, 5 A直流电源驱动, 典型工作温度1160 ℃, 电光效率为80%, 辐射波长覆盖2—25 μm范围. 为了在光学设计软件中较为全面而准确的建立红外热辐射光源模型, 以IR-SI272红外热辐射光源为例自行开展辐射通量空间分布情况测量. 考虑商用红外热辐射光源柱/管状辐射体的实际情况, 分别以柱/管状辐射体的中心对称轴夹角和柱/管状辐射体纵深长度划分角度和长度微元, 并分别以角度和长度微元为单位进行离散采样, 实验装置分别如图1[15]和图2所示. 红外探测器选用GENTEC-EO公司型号为XLP12-3S-H2-D2的热电堆探测功率计, 响应波段为0.19—20 μm, 探测孔径为12 mm, 探测孔径外罩内径25.5 mm的管状光阑, 用于消除由空气湍流造成的功率波动并遮挡大角度杂散光. 为保证测量结果的相对准确, 实验中应避免内、外部环境条件的变化对红外热辐射光源及红外探测器产生干扰. 实验全程在暗室中进行, 无其他光源、热源, 且控制暗室内温度为26 ℃, 相对湿度为55%.
图 1 红外热辐射光源发光特性测量实验装置图[15] (a)示意图; (b)实物照Figure1. Experiment configuration for measuring the light distribution characteristic of IR thermal radiation light source[15]: (a) Schematic diagram; (b) photograph.
图 2 红外热辐射光源辐射通量长度分布测量实验装置图 (a)示意图; (b)实物照Figure2. Experiment configuration for measuring the radiant flux distribution on the length of the IR thermal radiation light source: (a) Schematic diagram; (b) photograph.
图1中红外热辐射光源辐射区域中心位于转台旋转轴上, 使二者可以进行同步旋转. 如文献[15]所述, 红外探测器与红外热辐射光源同轴等高, 通过调整二者间距保证所探测的辐射通量数值始终位于有效量程范围内. 测量过程中, 红外热辐射光源发出的光线在传播特定距离后被红外探测器接收, 通过旋转转台改变红外热辐射光源柱/管状辐射体的中心对称轴与探测器之间的夹角, 进行均值和归一化处理后可得到各角度下的相对辐射通量, 即光源模型的角度微元[15].
图2实验中选用铝制反光碗作为锥形光管, 入口直径为4 mm, 出口直径为55 mm, 具有85%—90%的反射率[16]. 红外热辐射光源垂直向下固定于集成水冷的灯座中, 灯座固定于精密位移台上, 可带动红外热辐射光源进行上下移动. 红外热辐射光源正下方放置红外探测器, 探测器靶面上安装锥形光管(反光碗), 红外热辐射光源、锥形光管与探测器靶面三者同心放置. 锥形光管入口比柱状光源横截尺寸略大, 出口尺寸与红外探测器靶面尺寸相当, 可有效地收集红外热辐射光源所辐射的能量, 同时遮挡了其他杂散光. 测量过程中通过位移台改变红外热辐射光源伸入锥形光源中的长度, 获取红外热辐射光源柱状辐射体发出的辐射通量在柱状纵深长度上的分布, 进行等距差值和归一化后可得到纵深各单位长度下的相对辐射通量, 即光源模型的长度微元.
在ZEMAX软件中, 径向光源为轴对称光源, 可基于任意强度-角度数据进行样条拟合, 可满足自定义光源发光特性的需求. 光源模型以径向光源为微元, 按照HAWKEYE IR-SI272红外热辐射光源的几何尺寸进行建模, 并以实验测得IR-SI272红外热辐射光源的相对辐射通量-长度微元分布和相对辐射通量-角度微元分布对每一径向光源的光线总能量和能量方向分布进行设置. 综上, 建立了HAWKEYE IR-SI272红外热辐射光源的三维模型, 如图3所示.
图 3 HAWKEYE IR-SI272光源 (a)三维布局图; (b)实物照Figure3. HAWKEYE IR-SI272: (a) 3D layout; (b) photograph.
2
2.2.锥形光管
相比普通反射镜, 锥形光管能收集更多角度范围的光线, 在非相干光源光学系统中不可或缺. 锥形光管的应用使系统对辐射源有更高的能量收集效率, 通过改变锥形光管内壁反射面的面型, 可以控制光场分布, 进行光场整形. 由2.1节可知, 红外热辐射光源向全空间4π立体角辐射能量, 几乎没有方向性. 具有一定方向的辐射源出光方向对准普通光管入口, 便可实现耦合; 方向性较差的辐射源可以利用抛物面、椭球面等二次曲面的特殊光学特性, 将辐射源置于锥形光管内部, 二次曲面焦点上, 提高能量收集效率. 二次曲面可以使光源辐射出的部分光线收敛角度, 作为锥形光管的初始结构. 二次曲面面型由标准矢高方程[17]描述:| Conic系数取值 | 面型 |
| k = 0 | 球面 |
| k < –1 | 双曲面 |
| k = –1 | 抛物面 |
| –1 < k < 0 | 椭球面 |
| k > 0 | 竖椭球面 |
表1二次曲面conic系数-面型对应表
Table1.Correspondence between conic coefficient of quadric surface and surface type.
理想情况下, 光学元件传输能量没有损耗, 即光学扩展量不变. 对于光声光谱仪用红外热辐射光源, 辐射体为圆柱体, 设直径为d, 长为h, 则光源的光学扩展量简化近似为
图 4 带锥形光管光源的光学扩展量分析Figure4. Analytical diagram of the Etendue of light source with tapered light pipe.
由方程的几何关系可得远焦点处光斑半径和光线发散角:

在进行锥形光管设计时, 还需要综合考虑光学元件的材料折射率和光程差, 借助反射光焦公式[19]:
2
2.3.光声池
光声光谱仪用红外光源系统与光声池衔接使用, 共同工作. 光声池的几何尺寸作为红外光源系统光学设计上的目标光场分布, 依据相关尺寸进行光学建模, 光声池光学模型如图5所示. 根据光声光谱仪用红外光源系统光学指标, 对应于入池光功率、出池光功率和侧壁光功率分别在光声池入口、出口和圆柱侧壁设置探测器. 其中光声池入口和出口位置的探测器为圆形, 半径为5.75 mm, 与光声池圆柱腔体截面尺寸相同, 所探测到的功率数值分别为光声池入池和出池光功率. 在光学模型中可直接获取侧壁光功率. 光声池圆柱腔体侧壁的柱面为半径为5.75 mm、长 21 mm的圆柱管, 将其设定为物体探测器, 其功率值即为光声池侧壁光功率.
图 5 光声池三维布局图Figure5. 3D Layout of photo-acoustic cell.
3.1.评价函数
利用评价函数对自动优化过程进行约束, 主要包含系统光学指标优化方向和结构控制两部分. 光场整形系统的光学指标主要有入池光功率和侧壁噪声比. 其中入池光功率为进入光声池的总光功率, 可直接进行调用, 侧壁噪声比可由侧壁光功率与入池光功率做商而得. 结构控制部分主要用于限制光学元件之间的相对位置, 避免在自动优化过程中因位置的调整造成穿模, 违背实际物理状态, 无法进行后续实验装调. 主要包含红外热辐射光源、锥形光管、斩波器的相对位置.2
3.2.优化结果
根据前文计算得到的初始结构参数及编写的评价函数[20], 使用阻尼最小二乘法和正交下降法, 以光学元件相对位置和锥形光管反射面各面型参数为变量进行局部优化. 设计结果如图6所示, 非球面反光碗光学表面参数如表2所列. 优化设计后, 系统入池光功率为1.38 W, 侧壁噪声为17.7%, 光声池截面能量分布情况如图7所示.
图 6 光源光学系统 (a)三维布局图; (b)实体渲染图Figure6. Total reflection optical system: (a) 3D layout; (b) shaded model.
| 半径/mm | 曲率半 径/mm | 圆锥系数 | 二阶非球 面系数 | 四阶非球 面系数 |
| 12.0 | 9.896 | –0.634 | 5.275×10–3 | –3.870×10–3 |
表2非球面反光碗光学表面参数
Table2.Optical surface parameters of spherical reflective bowl.
图 7 光声池截面能量分布Figure7. The cross-section energy distribution of the photo-acoustic cell.
4.1.光学参数测试
根据设计完成相关光学与机械元件加工后, 进行光声光谱仪用光源系统的搭建与调试, 如图8所示. 比照光声池入口和出口的物理尺寸, 加工光阑, 按光声池位置固定前后两光阑, 等效为光声池的入口和出口, 如图9所示. 分别测量前后两光阑后的功率值, 作为光声池入池功率与出池功率. 实验测得光声池入池功率为1.32 W, 出池功率为1.06 W, 侧壁噪声为19.7%, 侧壁噪声问题得到了明显控制. 三维扩展光源光场整形系统光学参数指标实验值相比设计值略有下降, 主要考虑锥形光管光学反射面镀膜无法达到理想的全反射, 标称镀膜后反射率可达95%以上, 参数指标偏差在其范围内.
图 8 实验装置图Figure8. Experimental setup.
图 9 光源光学指标测量Figure9. Measurement of optical index of light source.
使用HAWKEYE IR-SI272光源配套的MC-234原厂聚光器替换光声光谱仪用三维扩展光源光场整形系统中自行设计加工的锥形光管, 进行对比实验. 在保证出池光功率, 即有用功率最高的条件下, 实验测得入池光功率为1.73 W, 出池光功率为0.86 W, 侧壁噪声为50.3%. 图10为采用两种光场整形系统时, 光源在光声池出口端和探测器上形成的光斑结果. 可见, 相较于原厂聚光器, 自行优化设计的锥形光更好的控制了边缘杂散光, 有效地提升了贯穿光声池的有用功率, 侧壁噪声问题得到了明显改善.
图 10 探测器光斑 (a)自设计锥形光管; (b)原厂MC-234聚光器Figure10. Light spot: (a) Self-designed tapered light pipe; (b) original MC-234 condenser.
2
4.2.气体浓度测试
完成光源系统的搭建后, 将其安装到光声光谱仪中. 在25 ℃恒温下对二氧化碳、一氧化碳、甲烷、乙烷、乙烯、乙炔六种典型气体进行不同浓度(0, 2, 5, 10, 50和100 μL/L)的检测, 记录微音器在不同气体不同浓度下所输出的多组电压值. 对光源系统的效果检验, 是通过气体检测下限来判断的. 气体检测下限的公式[21]为| CO2 | CO | CH4 | C2H6 | C2H4 | C2H2 | |
| 标准差 | 0.9866 | 0.5285 | 0.9118 | 0.7888 | 0.7255 | 0.6947 |
| 斜率 | 11.7280 | 0.7824 | 2.1508 | 2.9748 | 0.4488 | 1.8032 |
| 检测下限/(μL·L–1) | 0.2126 | 1.7076 | 1.0717 | 0.6703 | 4.0870 | 0.9739 |
表3光源激励下气体检测参数表
Table3.Gas detection parameters under excitation of light source.
