上周朋友拿了个问题来找我一根 5 米的光纤环想让环内稳定维持 30 mW 的循环光功率输入端到底该给多少毫瓦这种问题在环形激光器、光纤传感和谐振滤波结构里特别常见但它不是简单链路预算能回答的——环内的功率通量由耦合系数、往返损耗和谐振失谐三个因素共同决定通常在谐振点上要比输入功率高一个数量级。我顺手在 OptiSystem 里做了一个环形通量仿真把整套思路完整跑通了一遍。这篇就记录这次仿真的建模顺序、关键参数、结果解读以及几个比较隐蔽的坑。适合正在用 OptiSystem 做环形腔、微环、延迟线一类结构的同学参考尤其是第一次接触“环内通量”这种间接测量量的人。1. 环形通量到底在算哪个物理量1.1 先认清四个关键参数在 OptiSystem 里摆器件之前先把环形腔的基本定义说清楚。这里的“环形通量”不是电磁学里那个磁通量而是环内实际循环传输的光功率通量记作 P_circ。理解它之前必须要定义四个参数耦合功率系数 κ、往返场幅传输系数 a、单程相位失谐 δ以及输入功率 P_in。κ 是 2×2 耦合器里从直通波导耦合进环的功率比例取 0 到 1 之间a 是光在环内跑一圈之后的场幅剩余比例包含光纤衰减、熔接点损耗和耦合器附加损耗δ 是激光频率对应的环路相位偏移谐振点附近 δ 接近 2π 的整数倍P_in 就是光源真正进到耦合器输入口的功率不是激光器面板标称功率。这四个参数不齐全仿真就是空谈。很多人在 OptiSystem 里把器件一顿乱连最后看到输出端有个峰就以为仿真成功了其实环内通量这种量差一个损耗参数结果可能就是 4 倍和 40 倍的差别。1.2 谐振点上的通量放大公式最常见的全通型环形腔all-pass ring里稳态环内通量和输入功率的比值是P_circ / P_in κ a² / [1 a²(1-κ) - 2a√(1-κ) cos δ]当入射波长正好落在谐振点也就是 cos δ 1 时分母取最小值公式可以简化成P_circ / P_in κ a² / [1 - a√(1-κ)]²如果进一步忽略环内损耗也就是 a ≈ 1并且耦合系数不太大那么这个比值近似等于 4/κ。换句话说κ 取 0.1 的时候环内通量大约是输入功率的 40 倍κ 取 0.02 的时候理论上能放大到接近 200 倍。这跟直觉是相反的——你越是把光“关”在环里不让它跑出来环内的能量密度就越高。所以我才强调环形通量仿真本质上是在算一个间接量而不是单纯看透过率曲线。我的模拟场景里κ0.1光纤环长 5 m单模光纤衰减按 0.2 dB/km 算a≈0.9999谐振点上 P_circ 接近 38 倍 P_in。也就是说 1 mW 输入环内能扫出近 40 mW 的循环功率这个数值对后续选 EDFA、评估非线性效应阈值都很有参考价值。1.3 为什么要算通量而不是只看输出功率很多刚接触环形腔的人只盯着输出端功率直通端看到一个吸收谷就完事了。但环内通量决定了几件输出端功率看不到的事情非线性相移量、受激布里渊散射阈值、EDFA 的饱和状态以及环形传感腔的灵敏度。举个例子同样是 10 dB 的直通端消光比κ0.05 和 κ0.2 对应的环内通量放大倍数分别是 80 倍和 20 倍但光纤里的非线性效应却可能因此差出好几倍。做环内通量仿真目的就是在搭台子之前判断“环里到底会不会被非线性效应搞乱”这是只看光谱仪结果永远得不到的信息。2. 在 OptiSystem 里搭建环形腔顺序和参数2.1 器件选型与开环等效搭建我的做法是不追求一次成型先把环形腔拆成开环验证模型。组件清单如下CW Laser1550 nm1 mW线宽先设 1 MHz理想验证用2×2 定向耦合器耦合功率系数 κ0.1附加损耗 0.1 dB单模光纤段长度 0.005 km衰减 0.2 dB/km光功率计两个分别接直通端和环内输出点OSA用于粗看光谱但窄谐振峰不要依赖它后文会说原因。搭建顺序分四步CW Laser 输出口接到耦合器的 input port1耦合器的直通口接第一个光功率计耦合器的 tap 口接单模光纤段光纤末端接第二个光功率计先跑一次读出经过耦合器进环的功率、经过 5 米光纤衰减后的剩余功率得到单程往返损耗的实测等效值。为什么要用开环等效而不是直接把光纤末端接回耦合器的另一个输入口因为 OptiSystem 默认的数据流模式是按事件顺序把每个组件处理一遍直接闭环并不能自动迭代到稳态反而容易产生数值振荡这个坑我在后面专门讲。先测单程参数再用设计公式或者外部脚本叠加多次往返得到的稳态环内通量比硬连一个闭环可靠得多。2.2 全局参数设置长度单位与默认损耗OptiSystem 的全局参数默认是针对数字链路设置的做连续光仿真时特别容易踩单位换算的坑。光纤长度组件里的单位是 km5 m 环要填 0.005不是 5。我见过有人把 5 m 填成 5 km结果仿真里光纤损耗、色散和非线性效应全乱了出来的直通端曲线怎么看怎么怪。连续光仿真时Bit rate 保持默认 10 Gbps 即可Sample rate 我一般不动让软件自动算。真正要检查的是 OSA 的分辨率带宽和激光器线宽。谐振峰宽度很窄时OSA 的 Resolution bandwidth 如果设为 1 nm 级别峰谷会被展平看起来像没有谐振。建议先用 OSA 粗看再用参数扫描加光功率计逐点记录这才是能分辨亚 MHz 谐振峰的做法。我把常用的全局参数列在下面方便直接抄参数推荐设置说明Frequency193.1 THz对应 1550 nmBit rate10 Gbps影响采样时间窗连续光可保留Sequence length256时间窗长度谐振幅不够时可加大光纤长度0.005 km5 m 环单位是 km 不是 m光纤衰减0.2 dB/km1550 nm 单模常规值耦合器耦合比0.110% 进环耦合器附加损耗0.1 dB每个耦合器都会吃掉一点功率2.3 零散但关键的器件参数耦合器在 OptiSystem 组件库里对应 Optical Coupler 或 Directional Coupler重点检查两处Coupling ratio 和 Excess loss。部分自带模板的耦合器默认 Excess loss 是 0这个偏乐观你按 0.1 dB 设更接近实际。单模光纤组件里除了衰减还有 Dispersion 参数保留默认 16.75 ps/nm/km 即可连续光窄带场景下不影响通量。还有一个容易忽略的连接器损耗。某些模板会把 Connector loss 默认成 0.5 dB 一端一旦模型里出现多个连接器环内每往返一圈就额外掉 1 dB谐振增强直接从 40 倍掉到个位数。这是我在仿真里踩得最深的一个坑后面第 4 节专门展开。3. 参数扫描把谐振曲线一次跑透3.1 用 Parameter Sweep 扫入射波长环形通量仿真最有价值的部分是看“通量随失谐的变化曲线”也就是谐振曲线。在 OptiSystem 里我用 Parameter Sweep 把 CW Laser 的波长设成扫描变量从 1549.5 nm 扫到 1550.5 nm步长 0.01 nm。每跑完一步让光功率计记录当前波长下的读数最后把数据导出来画图。注意这里不用 OSA 作为主要读数工具。OSA 的分辨率带宽通常到不了亚 MHz 级别而 5 m 环的 3 dB 带宽可能只有几百 kHzOSA 会直接把谐振谷抹平。用功率计配合激光器逐点扫描等效分辨率取决于步长建议步长不超过 1 pm这样才能看清谐振峰的底。扫频时还要留意每步是否达到稳态。环形腔的充能时间大约等于 Q 值除以角频率精细度越高需要的点数就越多。这也是我为什么一直强调开环等效模型——闭环模型在扫频时每隔几步就出一次数值振荡非常让人头大。3.2 耦合系数对通量和带宽的影响把 κ 从 0.02 扫到 0.5你会看到一个经典 trade-off。κ 越小谐振点上环内通量放大倍数越大但谐振谷也越窄对波长和温度漂移越敏感κ 越大系统越容易对准但环内通量增强就没了。我用公式估算了一组理想值耦合系数 κ环内通量增强无源理想谐振峰宽度趋势工程意义0.02约 190 倍极窄高灵敏度传感难对准0.1约 38 倍中等论文最常用的折中0.5约 3.9 倍较宽容错好几乎无增强这组数据说明一件事环形通量仿真不是越大越好而是要根据系统容差选择耦合比。如果你后面还要接实际激光器激光器线宽本身有几十 MHz那拿 κ0.02 的环几乎等于宣判失败因为激光相位噪声就能把谐振峰盖掉。3.3 从仿真结果反推精细度仿真曲线出来之后我习惯顺手反推一下精细度 F用来跟理论值互相印证。精细度公式是F π√g / (1 - g)其中 g a√(1-κ)我的案例里 a≈0.9999κ0.1算出来 g≈0.9486F≈60。自由谱范围 FSR c/(n_eff L)1550 nm 下 n_eff 取 1.45L5 mFSR≈41.4 MHz。那么 3 dB 带宽就是 FSR/F≈0.7 MHz。如果你扫频步长能到 0.5 MHz 以下仿真曲线里就能看到一个清晰的尖锐峰如果步长比带宽还宽扫出来就是一个平滑大包啥也判断不了。这一步其实是在给仿真“验算”理论 F 和仿真峰宽若对不上基本可以断定模型里有损耗参数错了或者扫描步长不够细不用急着怀疑公式。4. 环形通量仿真最容易踩的四个坑4.1 直接把输出接回输入数据流模式下的振荡幻觉我第一次做环形通量仿真时图省事把光纤末端直接接回耦合器的另一个输入口形成一个真正的闭环。结果很有意思——环内功率随着 Sequence length 加大一路飙升看起来像激光器自激。但这不是真实的反馈而是 OptiSystem 数据流引擎按事件顺序把每个组件过了一遍信号在环里被反复累加本质是数值假象。解决方案就是回到第 2 节说的开环等效先测单程传输参数再用公式或者 MATLAB 脚本迭代多次往返。有人会问那 OptiSystem 是不是完全做不了闭环反馈我的经验是普通链路组件模式下最好别硬来非要闭环就要引入反馈控制和循环计数机制复杂度远高于收益。环形通量这种稳态量数学叠代比图形化闭环干净得多。4.2 双向通量混叠读数莫名偏高环形腔里如果两端都接了探测器有些人会把直通臂读数和环内功率相加当成总通量这是错的。全通环里直通臂和环内功率是同一光场的不同投影不是两个独立功率源。我第一次扫参时犯过这个错读出的峰值比公式算的高了一倍多还纳闷“增强怎么会超过理论极限”。正确做法是如果关心环内通量就只测环内环路中某个点的功率直通臂读数单独用于拟合透过率 T不要把两条支路的功率相加。环内通量是一个强度概念不是多端口总功率加法在这里没有任何物理意义。4.3 光纤长度单位与默认衰减单位错误是 OptiSystem 里最隐蔽的坑因为软件不报错只给你一堆看起来合理的数。光纤长度组件默认单位是 km5 m 环填 0.005 km但很多人填 5等效把环长放大 1000 倍。结果表面上看衰减不大但色散、非线性系数、群速度这些全按 5 km 算谐振峰位置和宽度全乱套。衰减参数也建议确认一下。OptiSystem 自带模板里单模光纤衰减常有 0.2 dB/km但如果你用的是通用光纤组件有的版本默认值是 0.3 dB/km5 m 影响不大可一旦环长几百米这个差异就不能忽略了。我建议每换一个模板就把所有无源器件的损耗参数展开看一遍别相信默认值等于物理值。4.4 耦合器附加损耗和连接器损耗决定 Finesse 的第一因素这是环形通量仿真里我踩过最深的一个坑。5 m 光纤本身的衰减只有 0.001 dB 量级在环内往返一圈几乎可以忽略。但耦合器 Excess loss 如果是 0.2 dB每过一次耦合器就掉 0.2 dB绕一圈损耗 0.4 dB这直接让场幅传输系数 a 从 0.9999 掉到 0.955精细度掉一个数量级。如果模板里还有连接器损耗默认 0.5 dB 每端绕一圈再扣 1 dB那 40 倍的通量增强就缩水到 4 倍左右。所以我做通量仿真时会把所有“看起来小”的损耗全部折算进去而不是只看光纤衰减。这也是仿真结果和实测总是对不上的头号原因。写一个检查习惯在布局里选中每一个无源器件右键看损耗把 Excess loss、Insertion loss、Connector loss 全部列出来加总成单圈损耗再代回公式验证仿真峰宽。这一步花十分钟能省下后面实测时三天排查时间。5. 把仿真结果带回工程两个可复用的场景5.1 场景一环形传感腔需要多少输入功率回到开头那个问题5 m 环、κ0.1、无源全通环想让 P_circ 稳定在 30 mW输入功率要多少按公式反推P_in≈30/38≈0.79 mW。但你不能直接按 0.79 mW 下单光源因为实际系统中失谐会偏移、耦合器附加损耗比仿真多、激光器线宽也会降低谐振深度。我自己的习惯是至少留 3 dB 余量也就是输入按 1.6 mW 以上选。如果环内还要插隔离器、滤波器每个附加元件的插损都要叠进这个计算里否则现场调半天就是不出通量。5.2 场景二MATLAB 联合验证与温度漂移OptiSystem 跑出来的通量曲线我会用 MATLAB 按第一节的公式复核一遍。方法很直接把 Parameter Sweep 导出的文本数据读进脚本与理论曲线对比。下面是一小段可供参考的验证脚本% 环形通量理论曲线 kappa 0.1; % 耦合功率系数 loss_ring_db 0.2 * 5e-3; % 5 m 光纤0.2 dB/km a 10^(-loss_ring_db/20); % 场幅往返传输系数 delta linspace(-2*pi, 2*pi, 2000); Pcirc kappa * a^2 ./ (1 (a^2)*(1-kappa) - 2*a*sqrt(1-kappa)*cos(delta)); [peak, idx] max(Pcirc); fprintf(峰值增强倍数: %.2f\n, peak);这段脚本不考虑耦合器损耗和连接器损耗只是理想光。如果你把实际器件的插损全部折算进 a就会发现峰值增强明显下降这时候你就知道该去调哪个器件的参数了。温度漂移也一样温度变化会改变光纤长度和折射率换算成等效失谐后继续用扫波长的方法在 OptiSystem 里重跑即可。1550 nm 光纤环的温度调谐系数大约在百 MHz 到 GHz 每摄氏度量级对 5 m 环、精细度 60 的结构来说1℃ 温度漂移足以让通量从峰值掉到半功率点以下。这个结论不算仿真误差而是环形腔的物理属性设计时一定要把温控预算考虑进去。最后说一个仿真之外的经验OptiSystem 里没有坏参数只有你没看的默认参数。环形通量这种间接量最容易出错的永远不是公式而是器件损耗、长度单位和测量端口选错。我第一次跑出来的 40 倍增强实际搭出来只有 9 倍回头一查就是连接器损耗默认值吃掉了大半。仿真给的是物理关系的放大镜真正决定结果的还是你对环内每一个 dB 的把握。祝大家的环都能稳定输出想要的通量。