单模光纤耦合仿真：光线追迹与物理光学传输的耦合效率差异

简介面向光学设计与光纤耦合仿真学习者的一份实操范例文档聚焦如何在 Zemax 中搭建并优化单模光纤到透镜的耦合系统。范例以康宁 SMF-28e 光纤与 SUSS MicroOptics FC-Q-250 微透镜阵列为对象围绕旁轴高斯光束计算、GBPS 优化操作数、单模光纤耦合效率能量传输与模式匹配等环节展开并附元件参数表与优化后的结果数据讨论孔径光阑设定、高斯切趾、数值孔径多种定义对结果的影响。读者可据此复现束腰尺寸与光纤到透镜间距的调优过程理解系统对称性与 Pick-up solve 的作用。压缩包共 1 个 doc 文档约 504KB篇幅紧凑是一份可直接对照 Zemax 界面操作的完整范例说明。目前已有 3603 人学习适合光通信器件、光纤耦合封装方向的工程师与相关专业学生作为建模练习与参数参考。1. 一条光线都没丢耦合效率却从 0.99 掉到 0.57把透镜间距从 2mm 拉到 20mm 之后发生了什么单模光纤耦合这件事最反直觉的地方不在透镜选型而在你怎么算。同一个系统两根单模光纤对着两组相同的微透镜透镜间距按实验取 2mm 时基于光线的耦合计算给出 0.99物理光学传输给出 0.994谁看都对得上。可只要把间距改到 20mm光线法照样报 0.99物理光学法直接掉到 0.57——接近一半的能量凭空蒸发而追迹图上的光线一条没少。这就是单模光纤耦合的玄学所在它耦合的不是能量到没到而是波前和光纤模式长得像不像。光线只带位置和方向不带位相所以它天然看不见衍射、看不见像差怎么把高斯波前揉歪。搞懂这个差异才能明白为什么一套完整的光纤耦合仿真要分三层来做也才能解释为什么软件里那个数值孔径参数换个定义填结果能差 50%。下面从一个具体的光纤耦合器文件出发把构型约束、三种算法的数学内核和容易出错的地方一层层拆开。2. 对称性先锁死一半自由度Pick-up solve 与浮动光阑在约束什么这套耦合器的物理构型很朴素源光纤端面 → 第一片微透镜 → 2mm 空气间隔 → 第二片微透镜 → 接收光纤端面。用的光纤在 1.31μm 波长下模场直径 9.2μm、纤芯直径 8.3μm、数值孔径 0.14微透镜是熔融石英基片上的透镜阵列单片透镜直径 240μm、节距 250μm、曲率半径 330μm、圆锥常数 0、数值孔径 0.17、基片厚 0.9mm、内部透过率优于 0.99。看上去参数不多但每一个都在影响后面能优化的自由度。真正决定这套模型能不能算出可信结果的不是这些数字本身而是建模时先把哪几件事钉死了。2.1 反演对称为什么倒过来用必须得到同样结果两组光纤和两组透镜在制造公差内是完全一致的这意味着整个系统天然具备空间反演对称和轴对称性——把耦合器掉个头耦合效率应该几乎不变。为了让这个物理事实在模型里成立建模时用了一个细节透镜到像面的距离不是独立变量而是通过 Pick-up solve 跟随物面到透镜的距离。前者一动后者自动同步。这一步看着像偷懒实际是在砍自由度。如果两端距离各自独立优化器完全可能给出一个左端 0.08mm、右端 0.15mm的解数值上效率略高一点但物理上不成立——因为实物是对称的你没法造出一个两端不对称的耦合器。用 Pick-up 把两个距离绑成一个等于把优化空间从二维压到一维同时也把倒过来用这个约束编码进了模型。初始值给的是 0.1mm这个数字不是算出来的是接近实际装配的经验值。原文说得很清楚后面优化时它会变。2mm 的透镜间距同理来自实验数据后续也会被优化。初始值只负责让优化器站在一个合理的起点上不负责正确。提示凡是物理上应该对称的系统都建议用 Pick-up 或等价约束把对称关系写死而不是优化完再检查两端是否一致。前者是约束后者是碰运气。2.2 float by stop size孔径光阑为什么交给透镜口径去定系统孔径光阑设在第一片透镜的后表面类型选根据光阑尺寸浮动float by stop size。这两件事必须一起理解。普通成像系统里光阑通常是一个明确的光圈口径固定入瞳位置和大小由它反推。但在这个耦合器里没有额外的光圈——限光的物理实体就是微透镜本身的 240μm 口径。把光阑放在第一片透镜后表面、并允许光阑尺寸浮动等于告诉软件这个系统的入瞳由透镜实际口径决定光阑口径不是一个待优化参数而是几何结构的结果。为什么放在后表面而不是前表面因为光束在前表面时还是从光纤刚出来的发散光口径很小经过透镜准直后在后表面附近达到最大横向尺寸这里限光效应最真实。数据也印证了这一点表面 3第一片透镜后表面上 1/e² 光束直径约 65μm表面 4第二片透镜前表面约 70μm而这两个面的物理半口径是 120μm——光束直径大约只有口径的一半多一点。也就是说绝大部分能量能通过但边缘确实会被切掉一部分这部分就是渐晕后面算系统效率时会体现出来。如果这里把光阑口径写死成一个不等于 120μm 的数或者把光阑放到别的面上渐晕量就错了而渐晕恰恰是系统效率里唯一与模式无关的损失项。2.3 数值孔径的三种定义0.14 和 0.09 其实不矛盾这是整个案例里最容易出错的一处也是最值得展开的一处。同一个数值孔径在光纤厂商数据表、光学软件的不同模块、不同的教科书中至少对应三种定义边缘光线倾角的正弦光强降到 1/e² 处的半角正弦光强降到峰值 1% 处的半角正弦。光纤厂商数据表给的 0.14 用的是第三种而软件在计算高斯模式耦合时用的是第二种算下来约 0.09。两者不矛盾差的是一个纯几何的比例因子。高斯光束的强度分布是 I(r) ∝ exp(-2r²/w²)令 I/I₀ 0.01解出 r/w √(-ln0.01/2) ≈ 1.517。而 1/e² 对应的正是 r/w 1。所以 1% 强度处的半角约等于 1/e² 半角的 1.517 倍。反过来验算模场半径 w 4.6μm波长 1.31μm高斯光束远场发散半角 θ λ/(πw) 1.31/(π×4.6) ≈ 0.0906 rad小角度下 NA ≈ sinθ ≈ 0.09。再乘 1.517 得 0.1375正好落在 0.14 附近。两个数字就此闭环。import math lam 1.31e-6 # 波长单位 m mfd 9.2e-6 # 模场直径单位 m w0 mfd / 2 # 模场半径 4.6 μm # 高斯光束远场发散半角 theta lam / (math.pi * w0) na_1e2 math.sin(theta) # 1% 峰值强度处的半角1/e^2 半径放大 sqrt(-ln(0.01)/2) 倍 k math.sqrt(-math.log(0.01) / 2) na_1pct math.sin(k * theta) print(f1/e2 NA ≈ {na_1e2:.4f}) # 约 0.0906 print(f1% NA ≈ {na_1pct:.4f}) # 约 0.1375这段代码做的是同一件事的两种表述先用模场直径反推发散角再用强度阈值换算两种 NA 定义。k就是那个 1.517 的换算系数它只跟高斯分布的形状有关跟具体光纤无关。实际建模时如果软件要求填 1/e² 定义的 NA而手头只有数据表的 0.14就必须除这个系数直接填 0.14等于把接收光纤的模式放大了一半重叠积分立刻掉下来。另外模型在孔径上加了高斯切趾来产生高斯分布。这只是一个近似——真实的单模光纤出射场是贝塞尔型而不是纯高斯但在模场直径附近两者吻合得很好。原文也注明后面会做进一步精确计算这正对应后面用物理光学传输时导入真实模式场的做法。3. 旁轴高斯与单模耦合两个积分把光到了和模式对上了拆开三层算法里前两层都还在光线的框架内只是对光束的描述精度不同。旁轴高斯光束只跟踪一个复参数 q单模光纤耦合计算则开始算重叠积分。两者给出的结果差得不多0.117mm 对 0.11mm但它们能回答的问题完全不是一个层级。3.1 从 9.2μm 模场直径到 4.6μm 束腰再到 65μm 准直直径旁轴高斯光束计算分析菜单下的物理光学分支里那个设置很简单束腰位置固定在表面 1和物面重合束腰半径 4.6μm也就是模场直径的一半。然后让这个高斯光束一路传过去。中间发生了什么可以从瑞利长度看出来。z_R πw₀²/λ π×(4.6μm)²/1.31μm ≈ 50.7μm。光纤到透镜的距离初始 0.1mm约等于两个瑞利长度此时光束半径已经涨到 w₀√(1(z/z_R)²) ≈ 4.6×2.24 ≈ 10.3μm。也就是说光在到达透镜时已经不是细小平行光了而是明显发散的。透镜把它准直回来。由曲率半径 330μm 和熔融石英约 1.45 的折射率估算对称双凸结构的焦距落在 0.36mm 上下把 0.09 的发散半角乘上去准直后的束腰半径约 33μm1/e² 直径 66μm——和文件里报出的表面 3 上 65μm 几乎重合。这个自洽性说明模型的光焦度设置是对的。表面 4 的 70μm 比 65μm 大一点这个大一点本身就是衍射的证据。准直后的束腰并不在透镜后表面而是在透镜后约一个焦距处之后光束开始重新发散2mm 的传输距离已经超过准直光束的瑞利长度量级约 2.6mm所以到达第二片透镜时必然比出射时宽一些。旁轴高斯能捕捉到这个变化因为它跟踪的就是复参数随距离的演化。3.2 GBPS 一行优化函数凭什么够用优化前有个明显的失配像面接收光纤端面处的光斑是 5.3μm而不是对称性要求的 4.6μm。这说明焦点没落在像面上耦合必然打折。优化用的操作数是旁轴高斯光束尺寸GBPS它返回指定面上的高斯光束尺寸可以直接作为优化目标。由于系统被 Pick-up 约束成一维变量只有表面 1 的厚度表面 5 的厚度自动跟随目标只有一个让像面上光束尺寸回到 4.6μm。所以优化函数编辑器里就只有一行。优化函数示意一行 GBPS Surf6(像面) Target4.6 Weight1 变量 表面1 厚度Pick-up 联动 表面5 厚度一行优化能收敛靠的不是运气是自由度已经被前一步的对称性约束砍干净了。这也是为什么建模顺序很重要先把物理约束写进模型再谈优化反过来先优化再补约束通常要推翻重来。优化结果光纤到透镜的距离从 0.1mm 变成 0.117mm。束腰回到了 4.6μm焦点落回像面。3.3 S 和 T能量传输与相干重叠积分是两个分母不同的比值旁轴高斯的局限在于它假设光束一路都是理想高斯完全不管孔径截断和像差。单模光纤耦合计算往前迈了一步它拆成两个量分别算。系统效率 S分子是通过入瞳、经过系统考虑渐晕、体吸收、偏振传输之后的能量分母是源光纤辐射的全部能量。它衡量的是能量留下了多少跟模式形状无关。接收端效率 T这是相干重叠积分描述汇聚波前和接收光纤模式的匹配程度。分子是接收光纤复振幅 Fr 与出瞳处波前复振幅 W 的重叠积分模平方分母是两个模的自身内积。因为两个函数都是复数这个积分保留了位相信息——位相差一点实部虚部对不上模平方就掉下来。总耦合效率可以理解为这两个量的乘积光得先到还得长得像。S 掉通常是口径或吸收的问题T 掉通常是像差或离焦的问题分开看才知道该往哪改。3.4 FICL 优化为什么系统效率没变接收效率涨了直接用 FICL 操作数优化它就是把这个乘积作为目标几个循环后光纤到透镜的距离变成 0.11mm比旁轴高斯的 0.117mm 略小。两个数字的差异来自评价标准不同GBPS 只要求像面光斑尺寸等于 4.6μm而 FICL 要求波前和模式在振幅和位相上都匹配。位相匹配的最优点未必是光斑尺寸最小的点。优化后有几个现象值得注意。系统效率几乎没变因为渐晕是由透镜口径和光束尺寸决定的0.007mm 的离焦改变不了通过口径的能量比例。接收端效率提升了因为重新聚焦让波前曲率半径和光纤端面的模式波前对齐了。这说明在这个系统里优化的收益全部来自位相匹配不来自能量收集。至于最大耦合效率能不能再往上推——通过加非球面、加额外表面可以但本例中已经接近这套光学结构的极限。透镜是球面、圆锥常数 0球差是固有限制。4. 物理光学传输位相图上的抛物线和四次曲线到底在报什么错前两层算法共用一个隐含假设光束始终保持高斯形状。当衍射可以忽略、像差轻微时这个假设成立一旦传输距离变长或孔径截断严重假设就崩了。物理光学传输POP把复振幅场真正在面上传播才让衍射效应这四个字有了具体含义。4.1 256×256 采样与 Auto 矩阵采样不足会让结果自己抖POP 的设置里表面 1 上是束腰半径 4.6μm 的高斯模式一路传输到像面在像面上用完全相同的模式算重叠积分。采样率填 256×256但关键动作是点一下 Auto 去设定起始矩阵。这个 Auto 不是可有可无的辅助按钮。物理光学传播把场离散成网格网格范围太小时边缘能量被强行截断等于凭空加了一个硬光阑网格太粗时高频位相信息球差那类四次曲线采样不足算出来的重叠积分会随矩阵大小漂移。Auto 的逻辑是根据系统里最小和最大的特征尺寸——光纤模场、透镜口径、传输距离引起的衍射展宽——自动选一个既能覆盖全部能量、又能分辨位相细节的窗口。手动固定一个数然后改传输距离是翻车的常见起点。顺带说一句这个 256×256 的网格在透镜间距 2mm 时够用间距拉到 20mm光束在中间空间里展宽到 1/e² 宽 0.15mm已经碰上了 0.12mm 的透镜半口径此时如果矩阵范围还按 2mm 情况设置第二个透镜上的截断和衍射都会算错。4.2 位相才是要看的量光强图几乎是骗人的POP 报告里最有信息量的不是辐照度分布而是位相分布。原因很直接辐照度几乎是理想高斯形M² 只有 1.082光看光强你会以为一切正常。位相不一样接收端模式失配的地方光强几乎为零位相却把误差全暴露出来。位相图上能看到两件事。一是抛物线形状它等价于离焦——波前曲率半径和光纤端面的模式波前不匹配对应一个二次位相项。二是叠在抛物线上的四次曲线成分它等价于球差——球面透镜边缘和中心的光焦度不同波前边缘被多弯了一点。再加上透镜边缘处位相曲线的突然跳变那是孔径截断造成的衍射。把这三样分开看才知道效率还能不能提、该改哪个面。如果位相主要是抛物线调距离就能解决如果四次项占主导调距离没用得上非球面。4.3 POPD 操作数与 20mm 间距0.57 是怎么算出来的POP 分析窗口的所有数据都可以被优化函数里的 POPD 操作数取到这比在窗口里读数有用得多因为可以拿去优化。本例中取到的关键项数据项含义数值0总光纤耦合效率0.9941系统效率0.9982接收端效率0.99610导引光束束腰4.57μm23有效光束宽度4.84μm26M²1.082系统效率 0.998 说明透镜口径切掉的能量不到 1%和前面光束直径只有口径一半多的判断一致。然后是那个关键实验把透镜间距从 2mm 改成 20mm。光线法FICL算出 0.99几乎不动——因为光束在两片透镜之间近似准直光线追迹看不到准直光束在自由空间里的衍射。POP 算出 0.57。原因在数据里传输 20mm 后高斯模式的 1/e² 宽度涨到 0.15mm而透镜半口径只有 0.12mm。光束已经比口径宽了第二个透镜的物理光阑结结实实地切掉了一圈能量切掉的部分还会衍射向中心扩散把原本干净的高斯波前污染掉。聚焦到接收光纤上的光束 M² 涨到 2.45形状已经不是高斯。这里有个细节很能说明两类算法的差别把 2D 外形图的 Y 方向放大后可以看到由于球差光线法给出的第二个透镜上的光束口径反而比第一个透镜上小。光线法相信这个结论POP 不相信——因为实际发生的是衍射展宽占主导。把光纤到透镜的距离固定、把透镜间距设为变量、用 POP 优化得到的最优间距是 2.072mm。实验用的 2mm 和它差不到 4%说明当初的实验设计是对的——但这个对是 POP 验证出来的光线法验证不了。用通用绘图工具扫一遍间距对耦合效率的影响还能看到效率随间距变化的曲线以及 M² 随间距的同步变化。4.4 偏振与镀膜14% 的能量去哪了前面的计算都默认没有界面反射和材料吸收。打开偏振计算在单模耦合和 POP 的设置里勾上并在系统设置里把入射光定义成 y 线偏振两边的效率都掉到 86%。值得注意的是掉的只有系统效率能量传输接收端模式耦合几乎不变。原因是菲涅尔反射的大小随角度变化很缓慢在这么小的角度范围内它只是均匀地砍掉一部分能量不改变波前的形状。波前形状不变重叠积分的归一化结果就不变。14% 这个数字可以手算验证熔融石英折射率约 1.45单面正入射反射率 ((1.45-1)/(1.451))² ≈ 3.4%。两片透镜四个面总透过率 ≈ (1-0.034)⁴ ≈ 0.871也就是损失 13% 左右再加上基片体吸收落在 14% 完全合理。n 1.45 # 熔融石英折射率1.31 μm 附近 R ((n - 1) / (n 1)) ** 2 # 单面菲涅尔反射率 surfaces 4 # 两片透镜四个空气-玻璃界面 T_total (1 - R) ** surfaces print(f单面反射率 {R:.4f}, 四面总透过率 {T_total:.4f}, 损失 {1-T_total:.2%})这段的作用是给14%一个可复现的来源而不是接受软件吐出来的数字。参数只有两个折射率和界面数。如果模型里加了胶合面或保护玻璃界面数要跟着改损耗估算才对得上。再给每个玻璃面加一层单层氟化镁增透膜耦合效率回升。这里体现的是设计里的一个取舍单层膜只能把反射率压到一个有限水平且最佳工作波长固定要在这个 1.31μm 单波长系统里榨出最后几个点单层膜够用宽谱系统就得上多层膜。5. 排查清单五个让耦合效率对不上的常见坑这一节记录的是实际建模时最容易反复踩的几处。每一条都按看到什么现象 → 底层原因 → 怎么改来写。5.1 现象耦合效率比预期低一大截怎么优化都上不去原因几乎总是 NA 定义混用。把光纤数据表上的 0.14 直接填进要求 1/e² 定义的位置接收光纤的模式被认为是发散的重叠积分里两个模式的宽度差了一半位相再怎么匹配也救不回来。旁轴高斯那边也会跟着出错——发散角被高估 1.5 倍准直后的光束直径、透镜上的光斑尺寸全部偏移。解决方式是先确认当前输入框要的是哪种定义再用 1.517 这个系数换算数据表值 ÷ 1.517 得到 1/e² 值× 1.517 反过来。更稳妥的做法是用模场直径反算 λ/(πw)两条路算出来的数应该自洽对不上就说明有一处填错了。5.2 现象POP 结果随采样矩阵变化改个网格效率就跳几个点原因是采样窗口没有随系统变化重新设定。窗口太小边缘能量被强行截断等效于多加了一个硬光阑窗口太大网格变粗四次位相项和透镜边缘的陡变采样不足。两种都会让重叠积分失真而且失真方向不一样所以结果看起来忽高忽低。解决方式是在每次改动传输距离或透镜口径后重新点 Auto并且在每个面的报告里检查实际采样点数和光束宽度的比值。经验上光束 1/e² 宽度至少要有十几个采样点透镜口径要完整落在窗口内。传输距离拉长时需要重新评估不能沿用短距离下的设置。5.3 现象光线法效率 0.99换成物理光学变 0.57不知道信哪个原因不是哪个算错了而是两者回答的问题不同。光线法在光束近似准直、衍射可忽略时是准的传输距离一长衍射展宽让光束碰到透镜口径被切掉的能量和随之而来的衍射光线法根本表示不了。解决方式是用传输距离和口径的比值做判据光束在中间空间的 1/e² 宽度接近或超过透镜半口径时必须切到物理光学。本例中 2mm 间距下光束 65~70μm、口径 240μm光线法可信20mm 间距下光束 150μm、口径 240μm光线法失效。做优化时更要注意评价函数用什么算法优化结果就带那个算法的偏差——用 FICL 优化出来的间距拿去 POP 里验可能不是最优。5.4 现象系统效率怎么算都是 1.0渐晕像是没生效原因通常是光阑设置不对。光阑口径被写死成一个大于透镜实际口径的值或者光阑放在了一个光束还没展开的面上软件就认为没有任何光被挡住。这个系统里如果光阑放在第一片透镜前表面那里光束刚从光纤出来只有二十几微米永远碰不到口径渐晕恒为零。解决方式是把光阑设在光束横向尺寸最大的那个面本例是第一片透镜后表面类型选根据光阑尺寸浮动让口径由几何结构决定而不是由某个输入框决定。设置完可以看系统效率和接收端效率的差值——差值为零通常就是渐晕没生效的信号。5.5 现象效率和文献/实测差十几个百分点各模块内部却都自洽原因是偏振和镀膜被关着。默认状态下软件算的是无界面反射、无吸收的理想传输四个玻璃面各 3.4% 的反射被完全忽略加起来正好十几个点。所有模块内部都一致因为它们一致地忽略了同一件事。解决方式是在单模耦合和 POP 两个设置里都勾上偏振计算在系统设置里明确入射偏振态然后跑一遍看系统效率掉多少。如果掉得和手算的菲涅尔损耗吻合说明这部分对了再给玻璃面加增透膜看效率回升多少。注意勾了偏振之后所有优化结果都要重跑一遍最优距离会轻微移动。6. 方法分层与代价一个系统的三个耦合效率该信哪一个回到最开始那个矛盾同一套参数0.99 和 0.57 并存。经过上面几层拆解现在可以给这三种方法划一条清晰的分工线。方法跟踪的量能看见什么看不见什么适用判据旁轴高斯光束复参数 q束腰位置、光斑尺寸、准直度孔径截断、像差、模式失配只需要一阶参数和初始尺寸单模光纤耦合光线 高斯模渐晕、能量传输、位相匹配衍射、非高斯模式光束宽度远小于口径物理光学传输复振幅场衍射、任意模式、像差位相计算代价高、设置敏感光束宽度接近口径或模式复杂这张表背后是一个更根本的取舍你优化用什么算法结果就带那个算法的世界观。用 GBPS 优化得到 0.117mm用 FICL 得到 0.11mm用 POPD 优化透镜间距得到 2.072mm——三个数都对但只有最后一个考虑了衍射。如果系统的实际工作点在衍射显著的区域前两个数就是精致的错误。计算代价也得算进来。旁轴高斯几乎是瞬时的FICL 一次优化几秒钟POP 每次评价都要在多个面之间做一次完整的场传播优化循环跑起来是前者的几十上百倍。所以合理的做法不是一律用最准的而是分层先用旁轴高斯把一阶参数和束腰定下来用 FICL 在光线框架内优化距离和口径最后用 POP 做终验和精调。哪一层发现假设不成立比如光束宽度逼近口径就把那一层的结论降级往上提一层重算。还有一个容易被忽略的取舍这个系统的最大耦合效率已经接近结构极限。透镜是球面、圆锥常数 0球差是固有的位相图上那条四次曲线就是它的签名。想再往上走代价是加非球面或增加表面数——从两片球面透镜变成非球面制造复杂度和成本都上一个台阶。在 1.31μm 单波长、耦合效率已经到 0.99 量级的情况下这笔账通常不划算。真正值得投入的地方往往在装配公差和波长漂移上而不是再挤那零点几个点。我自己的血泪经验是每次换算法验证之前先把 NA 定义、光阑位置、采样矩阵这三样重新核一遍。这三处任何一处错了换再准的算法也只是把一个错误的数算得更精确。从那以后我每次拿到一个新的光纤耦合模型不管它看起来多简单都强制走一遍NA 反算 → 渐晕是否为非零 → Auto 采样点检查这三步再开始谈优化。希望帮到你。本文还有配套的精品资源点击获取