全文HTML
--> --> -->目前, 空分复用已成为国内外研究热点. 在文献[3]中, 提出了一种新型空气沟槽辅助异质十二芯四模光纤, 该光纤结构可以显著降低串扰, 相较于普通沟槽异质光纤串扰可降低–66 dB/km. 在文献[4]中优化了一种七芯六模光纤, 成功获得低串扰, 低损耗特性, 可用于构建长途密集的空分复用传输链路. 在文献[5]中, 设计并制备出了一种低差分群时延折射率渐变型沟槽光纤, 可在C + L波段保持小于124 ps/km的低差分群时延下稳定传输六种LP模式.
多芯少模光纤的发展主要受到串扰、非线性系数和色散等问题的限制[6]. 其中芯间串扰和模间串扰对多芯少模光纤的限制最大, 过大的串扰会导致信号传输距离缩短和信号失真, 目前全世界的研究者都在致力于减小串扰以制造出实际可用的多芯少模光纤. 光纤的串扰可以通过设置结构差异来抑制, 例如控制相邻纤芯直径、纤芯间距、折射率分布及纤芯排布方式[6-9]. 然而受制于光纤的尺寸, 光纤设计的可能性就大大减小. 为抑制串扰, 目前多芯光纤的设计方案主要分为三种: 沟槽环绕型[10,11]、纤芯异质型[12]和空气孔辅助型[13]. 利用这三种结构可以将能量更好地限制在纤芯中从而减小纤芯之间的串扰. 三种结构相对比来说, 空气孔辅助型结构由于空气孔的存在, 在熔接时空气孔难以完全对准且拉制过程中空气孔容易塌陷造成光纤缺陷; 纤芯异质结构由于其纤芯各不相同在熔接问题上也较为困难; 沟槽环绕型结构在串扰抑制能力上具有一定的优势, 也利于优化光纤性能, 在制作工艺上也更为简单, 更符合实际应用. 对于采用模分复用的少模光纤, 单根纤芯中同时有多个独立且正交的模式, 为使光纤可以长距离稳定传输, 模式之间的串扰要求越小越好. 模式差分群时延(DMGD)是影响基于模分复用技术的少模光纤的重要因素之一[14,15]. DMGD是在少模光纤中由于不同模式具有不同的有效折射率, 各模式在纤芯中的传播速度不相同, 传输一段距离后形成模式间的时延[1]. 小DMGD的少模光纤模间串扰较大, 大DMGD的少模光纤模间串扰相对较小, 小的模间串扰能够保证模式间的相对独立, 而且较大的DMGD有利于接收端的解复用, 但DMGD不宜过大, 过大的DMGD不利于光纤长距离传输[1,13-19]. 在光纤的设计过程中, 设置合适的DMGD大小, 从而使光纤获得更优的性能. 上述方法对于设计多芯少模光纤, 实现长距离大容量稳定传输具有一定的指导意义.
本文设计了一种双沟槽环绕型十三芯五模光纤, 为了实现设计目标, 光纤设置的参数必须要以模式的传输数量为依据[6], 使得每个纤芯都可以稳定传输5个LP模式. 本文采用双沟槽辅助型结构, 相比于单沟槽结构可以将光能量更好地限制在纤芯中, 减小能量耦合; 优化双沟槽参数, 可以有效降低芯间串扰和模间串扰, 并利用最佳的DMGD减小模式间的相互影响. 仿真结果显示: 光纤在1.55 μm波长处传输60 km后, 对应于LP01, LP11, LP21, LP02, LP315个模式, 相邻纤芯间的串扰分别为–122.37, –114.76, –106.28, –100.68, –92.81 dB, 有效折射率差均大于1.026 × 10–3, 芯间和模间串扰可被有效抑制; 5个模式的非线性系数可保持在1.6 W–1·km–1以下. 综合光纤性能, 本文所设计的光纤可以应用于长距离大容量传输系统.
2.1.光纤结构
本文设计了一种双沟槽环绕型十三芯五模光纤, 其光纤截面图如图1所示. 该光纤基于空分-模分复用技术设计, 采用阶跃折射率分布结构, 光纤芯区使用掺锗石英材料提高芯区折射率, 沟槽包层采用掺氟石英材料降低沟槽折射率, 石英包层使用纯石英材料; 每一个芯区从内而外依次相间环绕2个石英包层和2个沟槽包层组成1个少模纤芯, 13个少模纤芯的参数相同.
Figure1. Schematic structure of a double-trench assisted 13-core 5-LP mode fiber.
相邻两个纤芯的折射率分布如图2所示. 其中: a为芯区半径; b1, b2为石英包层宽度; c1, c2为沟槽宽度; Λ为纤芯间距; R为光纤半径; 以上物理量单位为μm. 芯区与石英包层的相对折射率差为Δ1, 石英包层与沟槽包层的相对折射率差为Δ2.

Figure2. Refractive index profile of adjacent fiber-core.
2
2.2.光纤参数设计
对于传统阶跃折射率光纤, 纤芯中传输的模式数量由归一化频率V决定[1]:本文在参数设计过程中, 着重考虑了光纤抑制串扰的能力、有效模场面积大小和非线性系数. 对于多芯少模光纤, 抑制串扰的能力决定了光纤是否可以长距离稳定传输光能量, 非线性系数决定了光纤容量大小. 非线性系数与有效模场面积成反比, 通过增大有效模场面积可以有效降低非线性系数, 而增大有效模场面积最直接的方法就是增大纤芯直径, 因此本文会合理设置纤芯直径从而抑制非线性系数, 优化光纤参数使串扰达到目标水平. 表1给出了光纤初始参数.
a/μm | b1/μm | b2/μm | c1/μm | c2/μm | Λ/μm | R/μm | Δ1 | Δ2 |
6 | 2 | 2 | 5 | 4 | 42 | 100 | 0.017 | –0.01 |
表1光纤初始参数
Table1.The initial fiber parameters
本文使用基于全矢量有限元分析方法的Comsol Multiphysics仿真软件, 结合完美匹配层(perfectly matched layer, PML)边界条件对设计的光纤结构进行仿真分析. 控制光纤直径为200 μm, 优化调整光纤其他参数, 获得低串扰低非线性系数的双沟槽环绕型十三芯五模光纤.
3.1.芯间串扰计算方法
芯间串扰是多芯少模光纤中必须要解决的问题, 特别是用于长距离大容量传输的多芯少模光纤, 串扰必须要足够低, 光纤才能稳定传输. 目前, 计算串扰的方法主要有模式耦合理论和功率耦合理论, 其中功率耦合理论更适用于多芯少模光纤的串扰计算[2], 所以本文将采用功率耦合理论计算串扰, 并在仿真模拟时沿横向引入随机波动的传播常数, 以此更准确的计算纤芯之间的串扰. 图3给出了具有横向波动的两相邻纤芯之间的串扰计算模型.
Figure3. Crosstalk between two adjacent cores with random fluctuation along longitudinal direction.
相邻纤芯之间的串扰(XT, 单位为dB)计算表达式为[2,6]:

由于光纤在制作过程中会不可避免地产生制造误差, 本文假设相邻纤芯制造的相对误差为δa/a = 1%[10]; 由于次邻近纤芯间距较大, 次邻近芯间串扰足够小可以被忽略[8], 所以本文只计算中间纤芯和相邻纤芯的串扰值.
2
3.2.模间串扰的表征
模间串扰是影响光纤传输的重要参数之一. 对于微结构光纤, 纤芯与包层的折射率差要大于传统光纤, 这会使得微结构光纤的模间串扰小于传统光纤[13]. 目前有效抑制模间串扰的方法主要有两种[19]: 1)通过设计较大的有效折射率差, 减小模式间的能量耦合, 抑制模间串扰; 2)保持模式间的差分群时延足够小, 使用多输入多输出技术在光纤的接受端对信号进行补偿处理. 两种方法各有优劣, 本文将采用第一种方法控制模间串扰.根据文献[1, 13]所述, 同一纤芯中各模式之间的有效折射率相差较大且大于10–3时, 纤芯内的模间串扰可以被忽略, 而模式间的有效折射率差与MDGD关系密切, 具有较大的有效折射率差同时会具有较大的MDGD, 因此MDGD较大的少模光纤具有更低的模间串扰, 但是过大的MDGD会影响光纤的传输距离. 目前, 长飞公司制造出的通信用阶跃型少模光纤在1.55 μm处LP31与基模的最大MDGD为15 ps/m. 本文将综合考虑光纤性能, 合理设计适宜的MDGD.
MDGD的计算公式为[3]
2
3.3.沟槽结构对串扰的影响
首先本文对比了芯区外无沟槽、单沟槽和双沟槽三种结构对光纤芯间串扰的影响, 除沟槽数量外三种光纤的其他结构参数完全相同. 针对LP01模式三种结构的芯间串扰如图4所示.
Figure4. No/single/double trench structure crosstalk contrast of LP01 mode.
由图4可以看出, 单沟槽结构的串扰小于无沟槽结构, 双沟槽结构对串扰的抑制效果又明显优于单沟槽, 采用双沟槽结构可以增强光纤抑制串扰的能力. 确定双沟槽结构优势后, 继续研究了沟槽宽度和沟槽相对折射率差Δ2对串扰的影响. 计算了光纤在1.55 μm处传输60 km, 串扰随沟槽宽度和Δ2的变化趋势. 图5给出了Δ2 = –0.01时芯间串扰与沟槽宽度的关系, 图6给出了沟槽宽度c1和c2都等于4 μm时, 串扰和Δ2的关系.

Figure5. Relation between crosstalk and trench width at Δ2 = –0.01.

Figure6. Relation between crosstalk and Δ2 at c1, c2 = 4 μm.
由图5和图6可以得出芯间串扰随着沟槽宽度和Δ2的增加不断减小, 增加沟槽的宽度和Δ2可以有效的抑制串扰. 本文将选取最佳的沟槽宽度和Δ2, 使双沟槽十三芯五模光纤具有最优的串扰抑制能力.
2
3.4.纤芯分析
芯区半径对光纤的串扰和性能有着显著的影响, 分析芯区大小与光纤性能的关系是必不可少的步骤, 本文拟定对芯区半径a = 5—8 μm进行分析, 分析结果如图7所示.
Figure7. The relationship between crosstalk, Aeff, Δneff of five modes and core size at 1.55 μm: (a) The relationship between core size and crosstalk; (b) the relationship between core size and Aeff; (c) the relationship between core size and Δneff.
由图7(a)可以看出5个模式的芯间串扰呈现非线性变化, LP01模式的串扰随芯区大小变化具有一定的波动性, 芯区较小时LP01模式呈现低串扰特性, 但芯区半径超过7 μm后, LP01模式的串扰只有微小的变化; LP11模式的串扰呈现出随着芯区半径增大, 串扰缓慢增长的趋势; LP21和LP02两模式的串扰变化趋势基本相同, 但LP21模式串扰低于LP02, 随着芯区增大串扰在a = 6.5 μm处存在一最小值, 芯区半径超过7 μm后串扰变化趋于平缓; LP31模式的串扰随芯区半径变化最为明显, 随着芯区增大串扰迅速降低, 在芯区半径超过7 μm后, 串扰变化趋于平缓. 图7(a)中的串扰变化规律可以从模式的电磁能量分布角度解释, 同一大小的芯区对不同模式的能量束缚能力不同, 高阶模式的能量更容易泄露到纤芯外, 且高阶模的能量主要分布在芯区外围更容易引起能量耦合, 通过增大芯区可以增强对高阶模能量的束缚能力以达到减小高阶模式的芯间串扰的目的, 当芯区增大到一定程度后, 芯区对模式的束缚能力增长程度不在明显. 综上分析, 5个模式的串扰均有一个相同的特点: 芯区半径超过7 μm后, 串扰均能保持在较低水平, 并且串扰随芯区半径变化很小; 若想保持5个模式的串扰均保持在较低水平, 芯区半径最佳的选择范围是7—8 μm.
由图7(b)可知, 光纤中各模式的Aeff均与芯区半径成正比, 5个模式的Aeff增长趋势基本形同, 但LP02模式的Aeff大小却低于其他4个模式. 对于高折射率芯区, 其具有较强的集光能力, 光会被集中在芯区内传播, 当芯区半径较小时, 传输模式会被集中在小芯区内, 模式所占据的横向面积就会很小, 对应的有效模式面积也会较小; 当芯区半径增大时, 模式被集中在较大的芯区面积内, 模式的横向分布面积也会增大, 有效模式面积也会随之增大. 根据图7(b)分析, 若想获得较大的Aeff, 就要尽量扩大芯区半径.
图7(c)为相邻模式之间的有效折射率差, 随着芯区增大, 模式之间的有效折射率差减小. 从图7(c)中可以看出LP21和LP02模式之间的有效折射率很接近, 两者之间的差值相比于其他4个模式略低, 所以在考虑折射率差时应以LP21–LP02为标准. 根据上文所述, 各模式之间的有效折射率差大于10–3便可忽略模间串扰, 但过大的折射率差将会引起较大的MDGD, 所以应选择适宜的芯区半径以获得最佳的模间折射率差值, 从而获得足够小的模间串扰.
除芯区半径外, 芯区的折射率大小对光纤也有着重要的影响. 现对芯区的相对折射率差Δ1进行分析, 结果如图8所示.

Figure8. . The relationship between crosstalk, Aeff, ?neff of five modes and Δ1 at 1.55 μm: (a) The relationship between Δ1 and crosstalk; (b) The relationship between Δ1 and Aeff; (c) The relationship between Δ1 and ?neff.
由图8(a)所示, LP01模式的芯间串扰在Δ1 = 0.01–0.016区间内, 串扰在–110 dB上下波动, Δ1超过0.016后, 串扰随着Δ1的增大开始降低; LP11模式在Δ1为0.01—0.017范围内, 串扰呈现降低趋势, 但Δ1大于0.017后串扰变化不在下降, 呈现平稳趋势; LP21, LP02和LP31模式的串扰随着Δ1的增加不断减小. 由图8(a)分析可得, 纤芯折射率越大对串扰的抑制能力越强. 但芯区折射率的增大要受实验上掺锗浓度的限制.
如图8(b)所示, 芯区的相对折射率Δ1增大, 会减小模式的有效模式面积, 但减小程度并不强烈. 若想获得较大的模式面积, 需减小芯区的相对折射率Δ1.
由图8(c)可知, 芯区的相对折射率Δ1对模式之间的有效折射率差具有一定的影响, 模式之间的有效折射率差随着Δ1增大而增大, 但相比芯区半径, 影响效果较小.
综上分析, 较大的Δ1具有更优的串扰抑制能力, 但同时会降低模式的有效模式面积, 增大光纤的非线性系数.
a/μm | b1/μm | b2/μm | c1/μm | c2/μm | Λ/μm | R/μm | Δ1 | Δ2 |
8 | 2 | 2 | 4 | 4 | 42 | 100 | 0.015 | –0.008 |
表2光纤优化参数
Table2.The optimal fiber performance.
Modes | Crosstalk/ (dB/60 km) | Aeff /μm2 | MDGD/ (ps·m–1) |
LP01 | –122.37 | 147 | 0 |
LP11 | –114.76 | 134 | 5.855 |
LP21 | –106.28 | 125 | 13.452 |
LP02 | –100.68 | 87 | 15.799 |
LP31 | –92.81 | 118 | 22.314 |
表35个LP模式的串扰、有效模面积和MDGD(LPmn–LP01)
Table3.Estimated values of crosstalk, effective area and MDGD(LPmn–LP01) for 5-LP modes at 1.55 μm.

Figure9. Transverse mode profile for 5-LP modes.
芯间串扰受波长影响较大, 如图10所示, 在1.3—1.7 μm波段, 芯间串扰随着波长逐渐增大. 在1.3 μm处传输60 km, LP31模式的芯间串扰可以保持在–140.61 dB左右, 但在1.7 μm处, LP31模式的芯间串扰却达到了–66.75 dB, 若要在传输过程中保持低串扰, 应尽量选择在较低波段进行传输.

Figure10. Relation between wavelength and core-to-core crosstalk for the double-trench assisted 13-core 5-LP mode fiber.
对于模间串扰, 根据图11分析, 优化后的光纤结构相邻模式之间可以保持较大的有效折射率差, 模式区分度大, 模间串扰可以被有效抑制, 当波长大于1.54 μm, 5个模式间的有效折射率差均大于10–3, 可以忽略模间串扰. 相邻模式间的模式差分群时延也可保持在适宜大小, 可以在接收端解复用.

Figure11. The relationship between ?neff, MDGD of five modes and wavelength: (a) The relationship between ?neff and wavelength; (b) the relationship between ?neff and wavelength.
通过仿真优化, 双沟槽十三芯五模光纤的芯间串扰和模间串扰均可以被有效抑制, 在长距离传输后, 串扰依旧可以保持在低水平.
5.1.各模式非线性系数分析
非线性系数是衡量光纤容量的重要参数, 光纤的非线性系数的计算公式为[21,22]:光纤的非线性系数与有效模场面积成反比关系, 通过增大有效模面积可以有效降低非线性系数. 目前通信用的G.652和G.655光纤的有效模式面积在80 μm2左右, 长飞公司生产的同质七芯光纤的有效模场面积在100 μm2左右, 对于多芯少模光纤, 模式的有效面积最好保持大于100 μm2[3,6,7,23]. 本文设计的光纤结构在1.55 μm波长处各模式的有效模面积如表3所示, 除LP02模式的有效模面积为86.5 μm2, 其他4个模式的有效模面积均保持在110 μm2以上, 可以使光纤保持低非线性系数.
图12(a)和图12(b)分别给出了光纤中各模式的有效模式面积Aeff和非线性系数与波长的关系. 从图12(b)可以看出, 非线性系数γ随着波长的增大逐渐减小, LP02模式的非线性系数在5个模式中最大, 但也未超过1.6 W–1·km–1, 在1.55 μm波长处各模式的非线性系数为0.74, 0.82, 0.88, 1.26, 0.93 W–1·km–1, 5个模式均可以保持低非线性传输.

Figure12. The relationship between Aeff, γ of five modes and wavelength: (a) The relationship between Aeff and wavelength; (b) the relationship between γ and wavelength.
2
5.2.色散分析
通过增大模场面积减小非线性系数会对光纤的色散产生影响, 在设计时必须要考虑到光纤传输的色散问题[10], 折射率高掺杂也会对色散有明显的影响. 色散(D)与有效折射率相对于波长的二阶导数成正比, 通常由(11)式表示[3,6]
Figure13. Relation between dispersion and wavelength for 5-LP modes.
2