随着光纤通信的快速发展,有关光纤通信信息容量的需求也在逐渐扩大。少模光纤的模分复用[1]技术利用光纤中多个模式传递信息,可有效提高光纤网络的容量;同时,其大模场面积通过降低光功率密度,也能减小光纤激光器中的光学非线性效应[2]。为实现光纤激光器的单模工作,以改善光束质量[3],需要对光纤模式的功率进行分解。因此,对少模光纤模式进行精确测量是一项关键任务。此外,光纤模式的分析在色散补偿器[4]、光纤模式转换器[5]、偏振控制仪[6]等方面具有重要作用。
多种用于光纤模式测量的技术被陆续提出和采用。最早采用的是侧发射技术[7]。该技术通过蚀刻光纤来控制包层厚度和等效势垒厚度。这种控制方法既可激发光纤中的高阶模,也能测量其相对功率。2007年出现的低相干反射技术[8]不仅能测量模式的功率,还能测量少模光纤的群延迟以及色散值。随后出现的三镜环形谐振器[9]通过挤压光纤产生高阶模,使各光纤模式可以表示为谐振器的几个本征模式的叠加,通过测量谐振腔模式就可以计算出基模和高阶模的功率。有研究提出了基于空间和频谱分辨成像法(fiber mode measurement method based on spatial and spectral resolution),简称S2法[10-13]。它可以用来确定激发模式的横向分布、相对相位。该方法可以在光纤参数未知的情况下,采集被测光纤输出端横截面上的光强值,然后对该区域每个位置的光强值进行频率扫描,再根据扫描所得的光谱曲线做傅里叶变换[14-15]。利用傅里叶变换后的振幅时延图和相位时延图中的数据信息对各个模式进行求解,即可得基模和高阶模的光强分布,以及高阶模与基模的相位差[16-17]。但现有工作无法判断光纤模式是否存在简并,也无法将简并模式从光纤模场中有效分离出来[18-19]。
本文提出利用S2法来求解简并模式的测量方法。利用该方法可以判断某一光纤模式是否存在简并,并能够对存在简并的模式进行分析测量。这是对现有S2法的改进和扩展。本研究提出的简并模式分解方法在保持模式正交性的同时实现简并模式的有效分解,显著提升了分解精度。该方法为光纤模式测量[20]提供了新思路,有望推动该领域的发展。该方法首先对入射光场进行调控,完成不同光纤入射模场下S2法的测量;接着,以这些测量数据作为样本进行奇异值分解;最后,对光纤模式是否存在简并进行判断,并完成简并模式的分解。搭建了基于空间和频谱分辨测量的光纤模式测量实验装置。选用少模光纤进行实验测试,结果表明:采用基于空间和频谱分辨测量的光纤模式测量方法,可以确定少模光纤中激发模式的横向分布;基于机器学习理论,利用奇异值分解可判断模式是否存在简并模式、确定简并模式的数量,并求解其光场分布。
1 实验内容 1.1 实验理论根据电磁波理论,光纤中的导模电场为
| $ {E_{mn}}(x,y,\omega ,t) = {E_{mn}}(x,y) \cdot \exp ({\mathrm{i}}\omega \cdot {n_{mn}} \cdot L - {\mathrm{i}}\omega t) $ | (1) |
式中:
| $\begin{split} E(x,y,\omega ,t) =& \displaystyle\sum\limits_{mn} {{c_{mn}} \cdot {E_{mn}}(x,y) \cdot \exp ({\mathrm{i}}\omega \cdot {n_{mn}} \cdot L - {\mathrm{i}}\omega t)} + \\&{c_{01}} \cdot {E_{01}}(x,y) \cdot \exp ({\mathrm{i}}\omega \cdot {n_{01}} \cdot L - {\mathrm{i}}\omega t)\\[-1pt] \end{split} $ | (2) |
式中:
将基模的光强定义为参考量,高阶模的光强采用参考量来表示。求解出参考量及参考量和其他高阶模之间的关系,就能够求得各个高阶模的光强分布。从基模和高阶模之间的光谱干涉公式出发,推导出用于光纤模式测量的S2法的理论公式。其前提条件是基模能量占主导。在此条件下,高阶模与高阶模的干涉项可近似忽略。测量得到的光纤模式相位是高阶模与基模之间的相对相位。假设基模光场分布在光纤截面上是等相位的,可对光纤中不同模式的相位分布进行测量。
设基模振幅
| $ {a_{mn}}(x,y) = \frac{{{c_{mn}} \cdot {E_{mn}}(x,y)}}{{{c_{01}} \cdot {E_{01}}(x,y)}} $ | (3) |
利用式(3),可以求出基模与高阶模之间光强关系,具体关系式为
| $ c_{mn}^2 \cdot {I_{mn}}(x,y) = {\left| {{a_{mn}}(x,y)} \right|^2} \cdot c_{01}^2 \cdot {I_{01}}(x,y) $ | (4) |
将式(3)代入到式(2),可得
| $ \begin{split} E(x,y,\omega ,t) =& {c_{01}} \cdot {E_{01}}(x,y) \cdot \exp ({\mathrm{i}}\omega \cdot {n_{01}} \cdot L - {\mathrm{i}}\omega t)\\&\left[1 + \displaystyle\sum\limits_{mn} {{a_{mn}}(x,y) \cdot \exp ({\mathrm{i}}\omega \cdot {\tau _{mn}})} \right] \end{split}$ | (5) |
式中,
| $ \begin{split} I(x,y,\omega) =& E(x,y,\omega) \cdot E^*(x,y,\omega)= c_{01}^2 \cdot I_{01}(x,y) \cdot \left[ 1 + \displaystyle\sum_{(m,n) \neq (0,1)} a_{mn}(x,y) \cdot a_{mn}^*(x,y) \right] + \\ & c_{01}^2 \cdot I_{01}(x,y) \cdot \left[ \displaystyle\sum_{\substack{(m,n) \neq (0,1), \\ (p,q) \neq (0,1)}} a_{pq}^*(x,y) \cdot \exp[{\mathrm{i}} \omega \cdot \tau_{pq}] \right] + c_{01}^2 \cdot I_{01}(x,y) \cdot \left[ \displaystyle\sum_{\substack{(p,q) \neq (0,1), \\ (m,n) \neq (0,1)}} a_{mn}(x,y) \cdot \exp[-{\mathrm{i}} \omega \cdot \tau_{mn}] \right] + \\ & c_{01}^2 \cdot I_{01}(x,y) \cdot \displaystyle\sum_{\substack{(m,n) \neq (0,1), \\ (p,q) \neq (0,1), \\ (m,n) \neq (p,q)}} a_{mn}(x,y) \cdot a_{pq}^*(x,y) \cdot \exp[{\mathrm{i}} \omega (\tau_{pq} - \tau_{mn})] \end{split} $ | (6) |
式中:
| $ \begin{split} F(x,y,\tau) &= c_{01}^2 \cdot F_{01}(x,y,\tau) \cdot \left[ 1 + \displaystyle\sum_{(m,n) \neq (0,1)} |a_{mn}(x,y)|^2 \right] + c_{01}^2 \cdot \displaystyle\sum_{\substack{(m,n) \neq (0,1), \\ (p,q) \neq (0,1)}} a_{pq}^*(x,y) \cdot F_{01}(x,y,\tau - \tau_{pq}) + \\ & c_{01}^2 \cdot \displaystyle\sum_{\substack{(p,q) \neq (0,1), \\ (m,n) \neq (0,1)}} a_{mn}(x,y) \cdot F_{01}(x,y,\tau + \tau_{mn}) + c_{01}^2 \cdot \displaystyle\sum_{\substack{(m,n) \neq (0,1), \\ (p,q) \neq (0,1), \\ (m,n) \neq (p,q)}} a_{mn}(x,y) \cdot a_{pq}^*(x,y) \cdot F_{01}(x,y,\tau - \tau_{pq} + \tau_{mn}) \end{split} $ | (7) |
式中:
|
图 1 少模光纤仿真得到的干涉结果 Figure 1 The interference results obtained from the simulation of few-mode fibers |
图1(b)中群时延差
| $ {f_{mn}}(x,y) = \frac{{\left| {F(x,y,{\tau _{mn}})} \right|}}{{\left| {F(x,y,0)} \right|}} = \frac{{\left| {{a_{mn}}(x,y)} \right|}}{{1 + \displaystyle\sum\limits_{(m,n) \ne (0,1)} {{{\left| {{a_{mn}}(x,y)} \right|}^2}} }} $ | (8) |
将上式化简得
| $ \begin{split} \left| {{a_{mn}}(x,y)} \right| =& a{{{f_{mn}}(x,y)\left[1 - \sqrt {1 - 4\displaystyle\sum\limits_{(p,q) \ne (0,1)} {f_{pq}^2(x,y)} } \right]\cdot}}\\&{{ \left[2 \cdot\displaystyle\sum\limits_{(p,q) \ne (0,1)} {f_{pq}^2(x,y)}\right]^{-1} }}\\[-1pt] \end{split} $ | (9) |
利用式(9)可以求出基模与高阶模之间空间位置系数的模值。通过式(4)和式(9)可以求出少模光纤中所有模式在空间位置上每个点的光强值。具体关系式为
| $ \begin{split} \displaystyle\sum\limits_\omega {c_{01}^2 \cdot {I_{01}}(x,y,\omega )} =& \frac{1}{{1 + \displaystyle\sum\limits_{mn} {{{\left| {{a_{mn}}(x,y)} \right|}^2}} }} \cdot\\& \displaystyle\sum\limits_\omega {I(x,y,\omega )} \end{split} $ | (10) |
| $ \begin{split} \displaystyle\sum\limits_\omega {c_{mn}^2 \cdot {I_{mn}}(x,y,\omega )} =& \frac{{{{\left| {{a_{mn}}(x,y)} \right|}^2}}}{{1 + \displaystyle\sum\limits_{mn} {{{\left| {{a_{mn}}(x,y)} \right|}^2}} }} \cdot\\& \displaystyle\sum\limits_\omega {I(x,y,\omega )} \end{split} $ | (11) |
式中,
所采用的光纤为圆柱形。由于圆形具有对称性,且测试光纤是少模光纤,为了便于分析讨论,考虑设定一种特定的少模光纤。该光纤一共具有
| $ \begin{split} E(x,y,L,\omega ,t) = & [c_{11}^{\mathrm{o}} \cdot E_{11}^{\mathrm{o}}(x,y) + c_{11}^{\mathrm{e}} \cdot E_{11}^{\mathrm{e}}(x,y)] \cdot\\& \exp (i\omega \cdot {n_{11}} \cdot L - i\omega t) \end{split} $ | (12) |
式中:
| $ \left| {{a_{11}}(x,y)} \right| = \left| {a_{11}^{\mathrm{o}}(x,y) + a_{11}^{\mathrm{e}}(x,y)} \right| $ | (13) |
式中,
| $ \displaystyle\sum\limits_\omega {c_{11}^2 \cdot {I_{11}}(x,y,\omega )} = \frac{{{{\left| {a_{11}^{\mathrm{o}}(x,y) + a_{11}^{{{\mathrm{e}}}}(x,y)} \right|}^2} \cdot \displaystyle\sum\limits_\omega {I(x,y,\omega )} }}{{1 + {{\left| {{a_{02}}(x,y)} \right|}^2} + {{\left| {a_{11}^{\mathrm{o}}(x,y) + a_{11}^{\mathrm{e}}(x,y)} \right|}^2} + {{\left| {a_{21}^{{\mathrm{o}}}(x,y) + a_{21}^{\mathrm{e}}(x,y)} \right|}^2}}} $ | (14) |
式(14)与式(11)的区别为,用
本文介绍了一种数据处理方法−奇异值分解(singular value decomposition,SVD)法。它一般适用于数据降维、图像压缩、信号处理等领域,其计算复杂度与矩阵的规模有关,取决于矩阵的行数和列数。SVD法的优势在于它能够严格保证简并模式之间的正交性。为实现入射光场的调控,假设在窄频率带宽范围内光纤模式保持不变,在此范围内取多个子频带进行S2模式测量。将每次测得的结果进行归类,对不同入射模场条件下求得的同一模式的不同模场数据进行整理。将这些光场数据重新排列成长度为
| $ {\boldsymbol{A }}= {\left( {\begin{array}{*{20}{c}} {x_1^1}&{x_1^2}& \cdots &{x_1^M} \\ {x_2^1}&{x_2^2}& \cdots &{x_2^M} \\ \vdots & \vdots & \ddots & \vdots \\ {x_{{l^2}}^1}&{x_{{l^2}}^2}& \cdots &{x_{{l^2}}^M} \end{array}} \right)_{{l^2} \times M}} $ | (15) |
矩阵
基于空间和频谱分辨的光纤模式测量方法的实验装置如图2所示。其主要由可调谐激光器(TSL–550,Santec,1 480~1 630 nm)、裸纤适配器、少模光纤(纤芯直径d=18.5 μm,数值孔径0.13,长度L=50 m)、物镜(40×,数值孔径0.65)、偏振片(通光孔径ϕ=21.5 mm)、CCD相机(Bobcat–640–GE,Xenics,分辨率640×512)等组成。选择物镜数值孔径(0.65)大于测试光纤数值孔径(0.13)的目的是确保从测试光纤出射的光能够被物镜全部收集。激光器发出的光经法兰耦合进入待测少模光纤后,从光纤另一端出射的光依次通过物镜进行准直和放大,再经偏振片将各模式光的偏振态调整一致,最后由CCD相机采集各点的光强值。通过频率扫描,测得光强随频率变化曲线。对该曲线作傅里叶变换分析,求解出光纤中的各个模式。实验中,频率的扫描范围为193.050 1~193.422 5 THz,频率间隔为0.000 7 THz,频率采样点数为532个。选择该频率范围及间隔的目的是为了在频域信号中能够清晰分辨出各个模式的群时延差,从而提高群时延差测量的准确性。
|
图 2 基于空间和频谱分辨的光纤模式测量方法的实验装置示意图 Figure 2 Schematic diagram of the experimental setup based on a spatially and spectrally resolved optical fiber pattern measurement method |
实验采用四模阶跃型光纤,长50 m,纤芯直径为18.6 μm,数值孔径为0.13。其在1 500~1 600 nm共有4个模式:
|
图 3 基于空间和频谱分辨光纤模式测得的干涉结果 Figure 3 Interference results measured in the fiber mode measurement experiment based on spatial and spectral resolution |
测得光纤出射端面上某一点的光强随光源频率变化的曲线,结果如图3(a)所示。对该曲线作傅里叶变换,得到如图3(b)所示的振幅群时延差曲线。图3(b)所示曲线中,横坐标为0处的峰是由干涉光强中不随频率变化的常量经傅里叶变换后得到的,而另外3个比较高的峰则是由基模与高阶模的干涉效应所产生。由于基模能量占主导,所以高阶模与基模干涉形成的峰比较高,而高阶模与高阶模之间的干涉形成的峰比较低。选取振幅群时延差曲线中4个最高峰值,再结合
|
(a)~(d)为 |
图4展示了各模式的光强分布及相较于基模的相位差,与理论结果基本一致。
对于简并模式的处理,本文假设在窄频率带宽范围内光纤模式保持不变,实验中频率为188.09~193.42 THz(对应波长1 551~1 590 nm)。在该范围内取5个子频带(188.09~188.68 THz、189.27~189.87 THz、191.69~192.31 THz、192.31~192.92 THz和193.05~193.42 THz),分别进行5次S2模式测量。在整个测量过程中,保持测量装置和光纤静止,以提高测量系统的稳定性和测量结果的可靠性。将每次测量的光纤模式光强和相位结果中同模式的光场信息放入一个矩阵中,得到4个矩阵,然后对每个矩阵进行奇异值分解,可求出每个矩阵的奇异值和对应的向量,结果如图5所示。
|
图 5 简并模式分解实验结果 Figure 5 Results of experimental data in degenerate modes |
从图5可知,
为了验证该方法对噪声的鲁棒性,对
|
图 6 加入噪声(13 dB)后SVD分解结果 Figure 6 SVD decomposition results after adding noise (13 dB) |
在扫描过程中引入大小为13 dB的噪声后,利用重构算法对焦点场进行重构,并分析噪声对重构精度的影响。结果显示,在一定噪声范围内,重构算法的抗干扰能力比较理想。
4 讨 论采用一种基于空间和频谱分辨的光纤模式测量方法对少模光纤模式进行了测量和分析。利用光纤中传输的激发模式具有不同群时延的特点,建立了基模与高阶模之间干涉的理论模型,给出了模式重建分析的方法。本文的创新点在于:利用奇异值分解法对简并模式进行特征提取,根据提取的奇异值和对应向量来判断简并模式的个数,并实现简并模式的分解,对现有S2光纤模式测量方法进行了扩展。同时,利用相应的光学仪器搭建了一套完备的实验装置,对少模光纤的模式分布进行了实验研究。结果表明,采用本文提出的方法可以实现光纤模式的测量。与现有方法相比,该方法可以对模式是否存在简并进行判断,且能在保持正交性的同时对简并模式进行有效分解。另外,对光纤中可能存在的有效折射率相近的近简并模式而言,在传统S2测量方法中,需要采用更为密集的频率分格扫描来完成模式的分析。本文提出的方法有助于在较粗略的频率分格测量中判断光纤中是否存在近似简并模式,保证光纤模式实验分析的完备性。实验所面临的问题是,难以使基模能量占据主导地位。为解决这一问题,将单模光纤与少模光纤熔接,成功实现了少模光纤中基模能量占主导的状态,有效解决了上述实验难题。后续将针对不同偏振态(如相互正交偏振态等)开展模式测量工作,为相关研究和应用奠定基础。对不同偏振态模式的测量,能为通信网络和传感网络领域的研究提供关键数据支持,为该领域研究提供重要依据。
| [1] | BERDAGUÉ S, FACQ P. Mode division multiplexing in optical fibers[J]. Applied Optics, 1982, 21(11): 1950–1955. DOI:10.1364/AO.21.001950 |
| [2] | CAI P J, LIAO T Q. Multi-wavelength fiber laser based on highly erbium-doped fiber and nonlinear effects in SMF[J]. Optik, 2020, 206: 164257. DOI:10.1016/j.ijleo.2020.164257 |
| [3] | WU F Y, ZHONG S, WANG X, et al. Improving the beam quality of single-grating spectrally combined fiber lasers with adaptive optics[J]. Applied Optics, 2023, 62(31): 8248–8260. DOI:10.1364/AO.499838 |
| [4] | POOLE C D, WIESENFELD J M, DIGIOVANNI D J, et al. Optical fiber-based dispersion compensation using higher order modes near cutoff[J]. Journal of Lightwave Technology, 1994, 12(10): 1746–1758. DOI:10.1109/50.337486 |
| [5] | RAMACHANDRAN S, WANG Z Y, YAN M. Bandwidth control of long-period grating-based mode converters in few-mode fibers[J]. Optics Letters, 2002, 27(9): 698–700. DOI:10.1364/OL.27.000698 |
| [6] | RAMACHANDRAN S, DAS M, WANG Z, et al. High extinction, broadband polarisers using long-period fibre gratings in few-mode fibres[J]. Electronics Letters, 2002, 38(22): 1327–1328. DOI:10.1049/el:20020926 |
| [7] | SZCZEPANEK P S, BERTHOLD III J W. Side launch excitation of selected modes in graded-index optical fibers[J]. Applied Optics, 1978, 17(20): 3245–3247. DOI:10.1364/AO.17.003245 |
| [8] | HAMEL P, JAOUËN Y, GABET R, et al. Optical low-coherence reflectometry for complete chromatic dispersion characterization of few-mode fibers[J]. Optics Letters, 2007, 32(9): 1029–1031. DOI:10.1364/OL.32.001029 |
| [9] | ANDERMAHR N, THEEG T, FALLNICH C. Novel approach for polarization-sensitive measurements of transverse modes in few-mode optical fibers[J]. Applied Physics B, 2008, 91(2): 353–357. DOI:10.1007/s00340-008-3011-9 |
| [10] | NICHOLSON J W, YABLON A D, RAMACHANDRAN S, et al. Spatially and spectrally resolved imaging of modal content in large-mode-area fibers[J]. Optics Express, 2008, 16(10): 7233–7243. DOI:10.1364/OE.16.007233 |
| [11] | NICHOLSON J W, YABLON A D, FINI J M, et al. Measuring the modal content of large-mode-area fibers[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2009, 15(1): 61–70. DOI:10.1109/JSTQE.2008.2010239 |
| [12] | NGUYEN D M, BLIN S, NGUYEN T N, et al. Modal decomposition technique for multimode fibers[J]. Applied Optics, 2012, 51(4): 450–456. DOI:10.1364/AO.51.000450 |
| [13] | OTTO H J, JANSEN F, STUTZKI F, et al. Improved modal reconstruction for spatially and spectrally resolved imaging(S2)[J]. Journal of Lightwave Technology, 2013, 31(8): 1295–1299. DOI:10.1109/JLT.2013.2242430 |
| [14] | 胡丽荔, 冯国英, 董哲良. 基于空间和频谱分辨的光纤模式测量方法[J]. 红外与激光工程, 2015, 44(8): 2517–2522. DOI:10.3969/j.issn.1007-2276.2015.08.047 |
| [15] | HU H F, LI Z H, ZHU Y L, et al. Non-degenerate mode power decomposition in optical fiber without prior knowledge[J]. Optik, 2020, 206: 164354. DOI:10.1016/j.ijleo.2020.164354 |
| [16] | Gao H, Hu H, Zhan Q. Generation and applications of spectral-spatially correlated principal mode in multimode fibers[J]. Optica, 2024, 11(2): 205–213. DOI:10.1364/OPTICA.510202 |
| [17] | Sharma V, Brandmüller D, Bütow J, et al. Universal photonic processor for spatial mode decomposition[J]. Nature Communications, 2025, 16(1): 7982. DOI:10.1038/s41467-025-63359-5 |
| [18] | Sánchez-Montes A R, Singh S, Márquez A, et al. Nonlinear modal decomposition of structured light[J]. Optics Express, 2025, 33(19): 41261–41270. DOI:10.1364/OE.567825 |
| [19] | Mangini F, Ferraro M, Tonello A, et al. High-temperature wave thermalization spoils beam self-cleaning in nonlinear multimode GRIN fibers[J]. Optics Letters, 2023, 48(18): 4741–4744. DOI:10.1364/OL.497917 |
| [20] | Hu X, Wu J, Gao Y. A review of structural health monitoring for flexible composite materials[J]. Applied Composite Materials, 2025, 32(2): 431–471. DOI:10.1007/s10443-024-10285-2 |
2026, Vol. 48
Issue (2): 51-59


