基于VirtualLab Fusion的Herriott池建模与仿真全流程解析

做气体吸收测量的人，十有八九绕不开多通池。Herriott池算是其中结构最优雅的一种：两面球面凹镜相对一放，光在里面来回弹几十次，等效光程一下就拉到了几十米甚至上百米，而物理尺寸却只有巴掌大小。最近我在VirtualLab Fusion里完整走了一遍Herriott池的建模与仿真流程，从参数计算、组件搭建到光斑演化分析，过程比预想中曲折，但也确实值得整理出来。这篇文章就是我自己的实操记录，供做光谱检测、激光器腔外增强，或者其他需要长光程小体积方案的朋友参考。

1. 为什么是Herriott池：原理与参数设计

1.1 多通池的基本原理与选型

痕量气体检测的核心矛盾很简单：吸收信号强度正比于光程，但仪器体积和光路稳定性又反过来限制光程。想在几十厘米的空间里获得数十米的有效吸收路径，只能让光“折叠”起来。多通池就是干这件事的器件，常见方案有White池、Herriott池和Chernin池等。White池用三块球面镜实现光路折叠，光斑排布灵活，但结构复杂，装配难度高。Chernin池适合超长光程，但镜面利用率和对准难度都很高。相比之下，Herriott池的结构最为简单：两块曲率半径相同的球面凹面镜面对面放置，一束光从其中一块镜面边缘的小孔射入，在两块镜子之间反复反射，最后从另一侧出射。

Herriott池的优势在于，只要初始入射位置和角度合适，光斑在镜面上会自然排成对称的椭圆图案，出射位置可预测，不需要额外的折返镜组。实际工程里，它被大量用在TDLAS（可调谐半导体激光吸收光谱）、光声光谱和激光腔外增强系统中，主要就是因为它结构皮实、机械公差宽松，而且对装调的容忍度高。我最早接触它，是想把一台体积只有30厘米长的小型光谱仪吸收程做到10米以上，Herriott池几乎是唯一不需要复杂光学平台就能实现的方案。

1.2 关键参数与光斑演化规律

理解Herriott池，核心是搞懂两个参数：镜面曲率半径R和镜间距d。当d接近R时（也就是接近共焦条件），光线在两块镜子之间来回反射时，每次到达镜面的横向位置都会绕着一个中心点旋转一个固定角度。这个旋转角θ的近似关系为：

θ ≈ arccos(1 - d/R)

当d/R = 1时，θ约为90°，也就是说光斑每经过一次镜面反射，就在镜面上转四分之一个圆弧，四次往返后回到初始位置附近。这个旋转角决定了光斑在镜面上的排布方向和间距，也决定了有效反射次数。实际设计时，d/R通常在0.9到1.1之间取值，偏离共焦条件后，光斑轨迹会从规则椭圆变成更稀疏或更密集的图案。

除了d/R，入射光束在镜面上的初始位置和倾斜角同样关键。光斑在两面镜子上各自形成椭圆轨迹，如果入射点偏离镜面中心，椭圆中心跟着偏移；入射角度越大，光斑散开的范围越大，反射次数也越多。但角度过大会导致光斑超出镜面边缘，光学损耗急剧上升。反射次数和等效光程由初始倾斜角、镜面口径和入射点三者的关系共同决定，它们之间没有独立的解析式可以直接套用，这正是需要用仿真软件做参数扫描的原因。

我在设计初期先按经典公式估算了一组参数，然后才在VirtualLab Fusion里建模型验证。下表是我常用的一组初值，适合波长1.55 μm、镜面口径1英寸的常规配置：

对刚开始接触的人，我的建议是不要一上来就追求上百次反射，先按d = R的共焦条件跑通模型，确认光斑分布和出射位置都与预期一致，再去慢慢调整入射角增加反射次数。

2. VirtualLab Fusion中的建模思路与组件选型

2.1 几何追迹还是物理光学追迹

VirtualLab Fusion支持多种追迹引擎，第一步要明确选哪种。Herriott池内的光束来回反射几十次，每一次反射都会改变光束的横向位置和角度，还要考虑衍射、干涉和高斯光束的相位演化，纯几何光学追迹虽然能快速给出光斑轨迹，但无法正确模拟衍射损耗和干涉效应，尤其在入射光束接近衍射极限、镜面口径相对束腰并不大时，几何追迹得到的反射次数往往偏乐观。

我的做法是分两步走：先用几何追迹引擎做光路粗验证，确认光线确实能在两面镜子之间稳定往返，且不跑出镜面范围；再切换到物理光学场追迹引擎，观察高斯光束在多次反射后的束腰变化、光斑形态和出射能量分布。物理光学追迹的结果更接近实验，但计算时间明显更长，参数扫描时如果直接用物理光学，一次扫描可能要跑很久。先用几何追迹把参数空间缩小到一两个候选组合，再做物理光学验证，效率会高很多。

另外需要注意VirtualLab Fusion里的追迹设置中有“最大反射次数”或者“终止条件”这类选项。Herriott池动辄几十次反射，默认追迹次数可能不够，必须手动设置足够大的上限，否则光线会在某个反射点被判定为“传播结束”，仿真结果里反射次数肉眼可见地减少。

2.2 组件配置：光源、反射镜与探测面

VirtualLab Fusion里建模Herriott池需要用到这几类组件：光源、反射镜、孔径/挡光板和探测器。光源我习惯用高斯光束模型，波长设为吸收谱线对应值，束腰位置放在入射小孔外侧附近，束腰半径按实际光纤准直器输出估算。反射镜直接用球面反射镜组件，设定曲率半径R、镜面口径和反射率。高反膜层设为99.5%以上，但要留一个出射机制，否则光永远困在池内，探测器收不到信号。实际实验中光通常从镜面上的一个小孔或边缘射出，仿真里可以用一块不透明挡板模拟小孔，或者直接在小范围内放一个探测器统计出射能量。

组件布局上，我会把坐标系原点放在M1（入射镜）的中心。M1位于z = 0处，凹面朝正z方向；M2位于z = d处，凹面朝负z方向，即M1和M2相对而视。光源放在M1后方偏下位置，以一定倾斜角切入M1边缘的小孔。这里有个细节：小孔位置与入射光轴必须仔细对齐，入射光通过小孔后在池内以某个横向位置和角度开始往返，初始横向位置和角度就是决定光斑图案的两个自由度。

为了方便后续扫描，我建议给所有参数命名时带上后缀，比如d_mirror、R_curv1、tilt_x_deg。VirtualLab Fusion的参数管理器支持这种命名方式，参数扫描结果表里一眼就能看出哪个参数影响了反射次数。坐标布局清晰之后再去跑追迹，出问题也容易定位。

3. 实操：从零搭建Herriott池仿真模型

3.1 初始参数定义与光束入射设置

下面给出一个完整的实操流程，以我最近做的模型为例。工作波长选择1550 nm，这是光纤通信和气体检测都很常用的波段；镜面曲率半径R取100 mm，镜间距d先取100 mm，满足共焦条件。具体操作步骤：

1. 新建光学设置，默认单位设为毫米，波长输入1550 nm。
2. 在组件库中添加两个球面反射镜组件，分别命名为M1和M2。
3. 设置M1：曲率半径R = 100 mm，凹面朝向正z，孔径半径12.7 mm，反射率99.5%。
4. 设置M2：R = 100 mm，位于z = 100 mm处，凹面朝向负z方向，其余参数同M1。
5. 添加高斯光源：束腰半径0.2 mm，束腰位置在M1后方10 mm处，波长1550 nm。
6. 在M1镜面上定义入射小孔：x = 2 mm，y = 0，孔径半径0.5 mm。
7. 设置入射束方向：向M2方向倾斜，x方向倾斜角15 mrad，y方向倾斜角5 mrad。

这里解释一下为什么入射角要同时给x和y两个方向分量。x方向的倾斜角决定光斑在水平方向上的铺开范围，y方向倾斜角则让光斑轨迹在垂直方向上也产生位移，两者合成后光斑才能形成完整的椭圆轨迹。如果只有一个方向的倾斜角，光斑只在一条直线上来回移动，反射次数虽然多，但光斑会挤在一起，出射位置和入射位置重叠，实际很难分离，也不利于后续探测器接收。

设置入射小孔时，孔径不能太大，否则M1镜面的有效反射区域变小，光容易从小孔漏出去；也不能太小，否则入射光束被硬切掉一部分，产生额外衍射。0.5 mm半径的圆孔对0.2 mm束腰的高斯光束来说，基本只影响远翼，对追迹结果的干扰很小。

3.2 运行场追迹与结果分析

搭建完成后，先跑一轮几何追迹看光路是否闭合。我用上述参数跑出来的结果：光在M1镜面上形成了44个光斑，M2镜面同样44个，光斑整体呈椭圆排列，椭圆长轴横向分布在x方向。总反射次数约88次（两块镜面各44次），池物理长度100 mm，等效光程约8.8 m。这个结果和用旋转角近似公式算出的数量基本吻合，差距在3%以内。

物理光学追迹的结果要区别看待。场追迹下光斑不再是“几何点”，而是有尺寸有强度分布的亮斑；光斑中心位置与几何追迹一致，但光斑直径会随着反射次数增加发生周期性变化。这个现象来自高斯光束在池内的自再现演化——每次反射相当于一次带有波前曲率半径变化的变换，束腰位置在M1和M2附近周期性地收缩和扩张。如果只关心路径长度，看几何追迹就够；如果关心出射光功率或模式质量，必须看物理光学结果。

追迹结束后，我会放三个探测器：M1镜面探测器、M2镜面探测器、出射小孔探测器。前两个统计光斑分布，第三个统计能量透过率。在反射率99.5%、反射次数88次的情况下，理论反射总损耗约(1-0.995)^88，看起来很小，但实际因为光束发散和衍射损耗，出射能量通常比入射能量低一到两个数量级。仿真里如果出射能量出现了“跳崖式”下跌，多半是某个反射点上的光斑已经超出镜面口径，光直接漏掉了，而不是反射膜损耗造成的。

3.3 参数扫描优化反射次数与光斑分离度

Herriott池设计最花时间的就是参数扫描。我用VirtualLab Fusion的Parameter Run功能，分别扫描镜间距d和入射角。先固定R = 100 mm，扫描d从90 mm到110 mm，步长1 mm，观察反射次数变化规律：

扫描入射角时，我固定d = 100 mm，从10 mrad逐步加到30 mrad，每步1 mrad。一个很明显的趋势是：入射角增大时反射次数总体增加，但到某个角度后会出现反射次数“跳变”下降——这是因为光斑椭圆的长轴超出了镜面边缘，大量光斑丢失，追迹算法判定传播结束。这个拐点就是工程上的极限入射角，再往上加角度的意义不大，反而会造成光斑外围拥挤、能量泄漏。

实际调整时没必要追求最大反射次数，我通常把入射角设定在拐点下方10%左右，这样等效光程足够长，又保留了一些机械装配余量。参数扫描结果可以导成表格，直接用于后续的结构设计。

4. 常见问题与排查技巧实录

4.1 光路发散、光斑逃逸的排查

建模过程中最容易遇到的问题是“光路跑飞”——几何追迹里光只反射了几次就从镜面侧面跑了，根本没有形成椭圆图案。我排查这类问题的顺序是：

1. 先查镜间距d是否在稳定区间内。Herriott池要求0 < d < 2R，也就是两个反射镜之间的间距不能超过两倍曲率半径，否则光路在几何上就是发散的。d/R越接近1，稳定性越好；偏离1太多，追迹必然发散。
2. 再查初始入射角。如果倾斜角设置过大，光线经过几次反射后横向位移迅速累计，很快就超出镜面孔径。解决办法是先设一个很小角度（比如5 mrad）跑通模型，确认光路闭合后再逐步增加角度。
3. 确认反射方向设置正确。凹面镜的反射方向是向内还是向外，软件里通常通过面法线方向控制。我自己就踩过这个坑：M2的凹面方向设反了，光线到达M2后直接反射回光源方向，而不是形成往返振荡。

4.2 仿真反射次数与理论估算对不上

几何追迹得到88次反射，理论估算却是90次，这种偏差大概率来自追迹终止条件和入射点偏心位置的设置。VirtualLab Fusion在光斑位置到达孔径边界时会判定反射结束，而理论估算假设镜面无限大，没有边界限制。因此仿真反射次数略低于理论值是正常的，工程上只要偏差在5%以内就可用。

另一种情况是反射次数突然比预期少一半。这时检查入射小孔或探测器的尺寸，它们如果被放在光路上作为挡光板，会把本该继续反射的光提前截断。我在早期版本里曾在出射探测器上设置了大口径接收面，结果它把M2镜面的一部分反射光也挡住了，反射次数直接从44次降到17次。

4.3 实用避坑清单

最后整理一份我反复用到的避坑清单，每一条都是实际调试中踩出来的：

• 反射镜的曲率半径和间距必须作为参数单独管理，不要在模型里写死数字，否则参数扫描时会乱套。
• 入射光束的束腰位置最好放在入射小孔外侧，不要放在池内，否则场追迹时的初始波前曲率会影响反射后的光斑形态。
• 追迹引擎的采样密度要足够。物理光学追迹时，网格过大光斑边缘平滑度差，网格过小运算时间剧增。一般先粗后细，确认结果趋势稳定再加密。
• 反射率不要设置为100%，实用仿真里至少保留0.5%左右的损耗。完全无损耗模型不仅不符合实际，还会让光斑强度分布失真。
• 每次调整参数后先跑几何追迹快速验收，再跑物理光学。几何追迹失败的情况下，物理光学结果再漂亮也没有意义。
• 记录每一组参数对应的反射次数和等效光程，不要只盯着屏幕上的光斑图看。做过几十组扫描之后就明白了，表格比截图可靠得多。

写在后面

这套建模流程我陆续跑了快两周，最深的体会是Herriott池的仿真难点不在软件操作，而在参数间的耦合关系。入射角、镜间距、镜面孔径这三者从来不是独立变量，只调其中一个参数很难得到预期效果。我的个人习惯是先固定d = R做基准模型，把入射角作为主变量，让光斑铺开之后再回头看镜面利用率，最后才回过来调整d/R比值。你在跑模型的时候如果发现光斑分布不理想，也可以试着倒过来，先定入射角，再反推d/R值。 VirtualLab Fusion里参数扫描本身不复杂，真正花时间的其实是判断每一组参数为什么会出现某个结果，这部分只能靠多跑多对比慢慢积累感觉。