WO2024055360A1 - Multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus - Google Patents

Multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus Download PDF

Info

Publication number
WO2024055360A1
WO2024055360A1 PCT/CN2022/122044 CN2022122044W WO2024055360A1 WO 2024055360 A1 WO2024055360 A1 WO 2024055360A1 CN 2022122044 W CN2022122044 W CN 2022122044W WO 2024055360 A1 WO2024055360 A1 WO 2024055360A1
Authority
WO
WIPO (PCT)
Prior art keywords
mode
core
electric field
inter
range
Prior art date
Application number
PCT/CN2022/122044
Other languages
French (fr)
Chinese (zh)
Inventor
向练
芮昊宇
Original Assignee
苏州大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 苏州大学 filed Critical 苏州大学
Publication of WO2024055360A1 publication Critical patent/WO2024055360A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation

Definitions

  • the invention relates to the field of optical fiber technology, and in particular, to a method, device and computer storage medium for inter-mode crosstalk detection of multi-core, few-mode optical fiber.
  • SMF Single-Mode Fiber
  • SDM is implemented through Multiple-input Multiple-output (MIMO) transmission, using the spatial mode of Multi-mode Fiber (MMF) or multiple single-mode fiber cores as channels.
  • MMF Multi-mode Fiber
  • FMF single-mode fiber
  • DSP Digital Signal Processing
  • FM-MCF has become a hot topic of research in recent years.
  • FM-MCF can bring ultra-high transmission capacity, it still has two important problems that will affect transmission performance.
  • One is that different modes have different transmission paths in the fiber core, which will lead to time differences in each signal received at the receiving end, eventually resulting in bit errors. This is called modal dispersion.
  • the second reason is that coupling between modes occurs between different paths, thereby reducing transmission performance.
  • modal dispersion can be alleviated through DSP technology or the introduction of dispersion-compensating optical fibers, but the evaluation of mode coupling is not yet complete. Therefore, accurate evaluation of mode coupling facilitates the preparation of multi-core few-mode fibers.
  • the author uses a mode electric field with polarization mode coupling effects, which further increases the computational complexity. And when dealing with the coupling coefficient, it is regarded as a constant that has nothing to do with the transmission distance, which is not consistent with the actual transmission situation.
  • the theory proposed by the ISMXT model it is proved that the crosstalk between supermodes caused by polarization mode coupling can be ignored. Therefore, when calculating the crosstalk between supermodes, the theory in the Monte Carlo model can be further simplified. .
  • the author only considered the impact of longitudinal perturbation, eliminated some redundant terms, and calculated inter-mode crosstalk directly from the mode power.
  • the technical problem to be solved by the present invention is to overcome the problem in the prior art that the processing of mode coupling coefficients is inconsistent with the actual situation.
  • the present invention provides a multi-core few-mode optical fiber inter-mode crosstalk detection method, which includes:
  • inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory is replaced with the inter-mode coupling mode coefficient, and the optical fiber inter-mode crosstalk value is calculated.
  • the fiber parameters include:
  • the calculation of the electric field intensity within the core range and the electric field intensity within the cladding range based on the optical fiber parameters includes:
  • E z1 and E z2 respectively represent the amplitude of the longitudinal electric field within the fiber core and the amplitude of the longitudinal electric field within the cladding.
  • A jUC/ ⁇ a
  • j is an imaginary unit
  • C is a system constant.
  • n 1 is the core refractive index
  • is the core propagation constant
  • a is the core radius
  • represents the light wavelength
  • J f (x) represents the first kind of Bessel function of order f
  • K f (x) represents the second kind of modified Bessel function of order f
  • n 2 is the refractive index of the cladding
  • r represents the radial direction in the cylindrical coordinate system
  • represents the circumferential direction in the cylindrical coordinate system
  • E r1 and E r2 respectively represent the amplitude of the transverse electric field within the core range and the amplitude of the transverse electric field within the cladding range
  • J' f (x) represents the first derivative of J f (x)
  • the electric field strength within the core is calculated:
  • Em and En represent the electric field intensity within the core range and the cladding range respectively, and the amplitude of the electric field in the circumferential direction e r represents the unit vector in the r direction, express The unit vector in the direction, e z represents the unit vector in the z direction.
  • the method before calculating the inter-mode coupling mode coefficient according to the electric field intensity within the core range, the electric field intensity within the cladding range and the injection power, the method includes:
  • the single-core transverse polar coordinate system (r, ⁇ ) is expanded to the dual-core transverse polar coordinate system (R, ⁇ ), D represents the core distance.
  • the inter-mode coupling mode coefficient is calculated based on the electric field intensity within the core range, the electric field intensity within the cladding range and the injection power:
  • is the angular frequency
  • ⁇ 0 is the vacuum dielectric constant
  • E m and H m represent the electric field intensity and magnetic field intensity in the core range respectively
  • E n represents the electric field intensity in the cladding range
  • e z represents The unit vector in the z direction
  • n 1 and n 2 represent the refractive index of the fiber core and the refractive index of the cladding respectively
  • the symbol * represents the conjugate operation of the matrix
  • P represents the injected power
  • the inter-core crosstalk calculation equation based on coupled mode theory is:
  • N is the number of segments after segmenting the optical fiber
  • d is the segment length after segmenting the optical fiber
  • ⁇ eq,mn,i is the equivalent phase mismatch of the i-th segment
  • k i (d) is the coupling mode coefficient between fiber cores.
  • the inter-core coupling mode coefficient k i (d) is replaced by the inter-mode coupling mode coefficient k mn,i , and the optical fiber inter-mode crosstalk value is calculated:
  • g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
  • the invention also provides a multi-core few-mode optical fiber inter-mode crosstalk detection device, which includes:
  • Fiber parameter and injection power acquisition module used to obtain fiber parameters and injection power
  • An electric field strength calculation module used to calculate the electric field strength within the core range and the electric field strength within the cladding range based on the optical fiber parameters
  • An inter-mode coupling mode coefficient calculation module used to calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
  • An optical fiber inter-mode crosstalk value calculation module is used to replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
  • the multi-core few-mode optical fiber inter-mode crosstalk detection device is applied to heterogeneous or homogeneous optical fibers.
  • the present invention also provides a computer-readable storage medium.
  • a computer program is stored on the computer-readable storage medium.
  • the computer program is executed by a processor, the above-mentioned method for detecting inter-mode crosstalk in multi-core, low-mode optical fiber is implemented. step.
  • the inter-mode crosstalk detection method of multi-core few-mode optical fiber first calculates the electric field distribution of each mode according to the fiber parameters, then calculates the inter-mode coupling coefficient according to the electric field distribution and injection power, and finally re-derives based on the coupled mode theory Compared with the Monte Carlo model and the ISMXT analytical model, the most prominent advantage of this method for calculating the amount of crosstalk between supermodes is that it takes into account the mode coupling coefficient that changes with distance, which further ensures the accuracy of the calculation and eliminates The redundancy brought by polarization mode coupling is more in line with the actual fiber laying situation and has a wider application range.
  • Figure 1 is a flow chart for the implementation of the inter-mode crosstalk detection method for multi-core few-mode optical fibers according to the present invention
  • Figure 2 is a schematic diagram of the longitudinal change of inter-mode crosstalk
  • Figure 3 is a schematic diagram of the dual-core polar coordinate system
  • Figure 4 is a schematic diagram of the crosstalk estimation model based on coupled mode theory
  • Figure 5 is a crosstalk comparison flow chart
  • Figure 6 shows the lateral electric field distribution of three modes: LP01, LP11a and LP11b;
  • Figure 7 is a schematic diagram of the transverse structure of a dual-core optical fiber
  • Figure 8 is a schematic diagram of nonlinear crosstalk changing with power
  • Figure 9 is a schematic diagram of the change of mode coupling coefficient with core spacing
  • Figure 10 is a schematic diagram of the inter-mode coupling coefficient changing with the core diameter ratio
  • Figure 11 is a schematic diagram of the inter-mode coupling coefficient changing with the refractive index difference
  • Figure 12 is a structural block diagram of a multi-core few-mode optical fiber inter-mode crosstalk detection device provided by an embodiment of the present invention.
  • the core of the present invention is to provide a method, device and computer storage medium for inter-mode crosstalk detection of multi-core, few-mode optical fibers, which take into account the mode coupling coefficient changing with distance and are more in line with the actual optical fiber laying situation.
  • Figure 1 is a flow chart of the implementation of the inter-mode crosstalk detection method for multi-core few-mode optical fibers provided by the present invention; the specific operation steps are as follows:
  • step-type multi-core few-mode optical fiber is to change the angle of light incident on the fiber core, so that each mode is transmitted along a different fold line trajectory in a single fiber core.
  • the specific schematic diagram is shown in Figure 2;
  • the optical fiber parameters include core refractive index, cladding refractive index, core radius, core propagation constant, core spacing and optical wavelength.
  • E z1 and E z2 respectively represent the amplitude of the longitudinal electric field within the fiber core and the amplitude of the longitudinal electric field within the cladding.
  • A jUC/ ⁇ a
  • j is an imaginary unit
  • C is a system constant.
  • n 1 is the core refractive index
  • is the core propagation constant
  • a is the core radius
  • represents the light wavelength
  • J f (x) represents the first kind of Bessel function of order f
  • K f (x) represents the second kind of modified Bessel function of order f
  • n 2 is the refractive index of the cladding
  • r represents the radial direction in the cylindrical coordinate system
  • represents the circumferential direction in the cylindrical coordinate system
  • E r1 and E r2 respectively represent the amplitude of the transverse electric field within the core range and the amplitude of the transverse electric field within the cladding range
  • J' f (x) represents the first derivative of J f (x)
  • the mode electric field can be generated into a cylindrical polar coordinate form, so that the electric field vector can be expressed as a superposition of three dimensions:
  • Em and En represent the electric field intensity within the core range and the cladding range respectively, and the amplitude of the electric field in the circumferential direction e r represents the unit vector in the r direction, express The unit vector in the direction, e z represents the unit vector in the z direction.
  • the mode coupling coefficient is closely related to the electric field distribution of each mode
  • the single-core transverse polar coordinate system (r, ⁇ ) in the electric field intensity within the cladding range is expanded to the dual-core transverse polar coordinate system (R , ⁇ ), D represents the core spacing;
  • is the angular frequency
  • ⁇ 0 is the vacuum dielectric constant
  • E m and H m represent the electric field intensity and magnetic field intensity in the core range respectively
  • E n represents the electric field intensity in the cladding range
  • e z represents The unit vector in the z direction
  • n 1 and n 2 represent the refractive index of the fiber core and the refractive index of the cladding respectively
  • the symbol * represents the conjugate operation of the matrix
  • P represents the injected power
  • the inter-mode crosstalk detection method of multi-core few-mode optical fiber first calculates the electric field distribution of each mode according to the fiber parameters, then calculates the inter-mode coupling coefficient according to the electric field distribution and injection power, and finally re-derives based on the coupled mode theory Compared with the Monte Carlo model and the ISMXT analytical model, the most prominent advantage of this method for calculating the amount of crosstalk between supermodes is that it takes into account the mode coupling coefficient that changes with distance, which further ensures the accuracy of the calculation and eliminates The redundancy brought by polarization mode coupling is more in line with the actual fiber laying situation and has a wider application range. On this basis, we can study the crosstalk characteristics in different communication systems. Based on the transmission characteristics of crosstalk in different working areas, the physical geometric characteristics of the fiber itself and the relationship between crosstalk and fiber parameters, we can then study ways to reduce crosstalk between different transmission modes. Theoretical methods.
  • step S104 further explains step S104:
  • This invention calculates inter-mode crosstalk based on coupled mode theory. First, it is necessary to introduce the principles of coupled mode theory and important solution methods:
  • the coupled mode equation is:
  • a n and A m are the slowly changing electric field complex amplitudes of core n and core m respectively
  • N is the number of fiber cores
  • K mn is the coupling coefficient from core n to core m
  • ⁇ m and ⁇ n are the propagation constants of m and n cores respectively
  • j is the imaginary unit
  • z represents the longitudinal propagation direction.
  • a general semi-analytic mathematical model is proposed, which accurately solves the linear coupled mode equation according to the distribution of phase matching points (Phase Machining Points, PMP).
  • K mn,i (z) is The coupling mode coefficient of the i-th section between core m and core n.
  • K mn,i (z) can be regarded as a constant K mn,i in the i-th section.
  • T is the coefficient matrix of the solution, which can be expressed as:
  • g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
  • the invention can provide a fast and accurate crosstalk estimation calculation method for multi-core and few-mode optical fiber communication systems, and its application range is wider.
  • This method is not only applicable to heterogeneous optical fibers, but also to real homogeneous optical fibers and ideal completely homogeneous optical fibers.
  • This method is based on the coupled mode theory that reflects the physical characteristics of optical fiber transmission, and re-derives the inter-mode coupling coefficient based on the random disturbances present in the fiber core. Therefore, it can more accurately reflect the changing process of inter-mode crosstalk than previous models. Based on this model, we can have a more forward-looking understanding of inter-mode crosstalk in few-mode multi-core optical fibers.
  • this embodiment verifies the accuracy of the above method through simulation, as follows:
  • transverse electric field distributions of the three modes LP01, LP11a and LP11b are shown in Figure 6, where LP01 is the transmission fundamental mode, and LP11a and LP11b are two mutually orthogonal polarization modes.
  • Figure 7 shows the lateral arrangement of the two fiber cores;
  • the influence of the model begins to be reflected, and the mode coupling coefficient calculated by the other two methods is an approximate empirical value and cannot accurately estimate the impact of random disturbances.
  • the inter-mode crosstalk results of the present invention have increased volatility. This is because the phase matching points between each mode are randomly distributed with distance, and the point with fluctuation indicates the occurrence of a phase matching point.
  • Figure 9 shows the process as the core spacing changes.
  • the core radius a 2.5 ⁇ m.
  • the inter-mode coupling coefficient gradually decreases, which shows that appropriately increasing the core spacing is one of the methods to effectively reduce inter-mode crosstalk.
  • Figure 9 considers three cases where the intrinsic effective refractive index difference is 0.21%, 0.31% and 0.41% respectively. It can be found that as the refractive index difference becomes larger, the inter-mode coupling coefficient will gradually decrease, and in the core When the spacing is small, the reduction is more obvious.
  • Figures 10 and 11 show the changes with core diameter ratio (ratio of core spacing to core radius) and intrinsic effective refractive index difference, respectively.
  • core diameter ratio ratio of core spacing to core radius
  • intrinsic effective refractive index difference respectively.
  • Figure 12 is a structural block diagram of a multi-core few-mode optical fiber inter-mode crosstalk detection device provided by an embodiment of the present invention; the specific device may include:
  • Optical fiber parameter and injection power acquisition module 100 used to acquire optical fiber parameters and injection power
  • the electric field strength calculation module 200 is used to calculate the electric field strength within the core range and the electric field strength within the cladding range according to the optical fiber parameters;
  • the inter-mode coupling mode coefficient calculation module 300 is used to calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
  • the optical fiber inter-mode crosstalk value calculation module 400 is used to replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
  • the inter-mode crosstalk detection device for multi-core and few-mode optical fiber in this embodiment is used to implement the aforementioned inter-mode crosstalk detection method for multi-core and few-mode optical fiber. Therefore, the specific implementation of the multi-core and few-mode optical fiber inter-mode crosstalk detection device can be found in the multi-core and few-mode optical fiber mentioned above.
  • the embodiment part of the few-mode optical fiber inter-mode crosstalk detection method for example, the optical fiber parameter and injection power acquisition module 100, the electric field strength calculation module 200, the inter-mode coupling mode coefficient calculation module 300, and the optical fiber inter-mode crosstalk value calculation module 400 are respectively used.
  • steps S101, S102, S103, and S104 in the inter-mode crosstalk detection method for multi-core few-mode optical fibers the specific implementation methods can be referred to the descriptions of the corresponding embodiments, and will not be described again here.
  • Specific embodiments of the present invention also provide a computer-readable storage medium.
  • a computer program is stored on the computer-readable storage medium.
  • the computer program is executed by a processor, the above-mentioned multi-core few-mode optical fiber inter-mode crosstalk is realized. Steps of the detection method.
  • embodiments of the present application may be provided as methods, systems, or computer program products. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment that combines software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk memory, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
  • computer-usable storage media including, but not limited to, disk memory, CD-ROM, optical storage, etc.
  • These computer program instructions may also be stored in a computer-readable memory that causes a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction means, the instructions
  • the device implements the functions specified in a process or processes of the flowchart and/or a block or blocks of the block diagram.
  • These computer program instructions may also be loaded onto a computer or other programmable data processing device, causing a series of operating steps to be performed on the computer or other programmable device to produce computer-implemented processing, thereby executing on the computer or other programmable device.
  • Instructions provide steps for implementing the functions specified in a process or processes of a flowchart diagram and/or a block or blocks of a block diagram.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Optical Communication System (AREA)

Abstract

The present invention relates to the technical field of optical fibers, in particular to a multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus, and a computer storage medium. The multi-core few-mode fiber inter-mode crosstalk measurement method comprises: first calculating the electric field distribution in each mode according to optical fiber parameters; then calculating an inter-mode coupling coefficient according to the electric field distribution and injection power; and finally re-deriving a calculation way of an inter-supermode crosstalk quantity on the basis of the coupled mode theory. Compared with the Monte Carlo model and the ISMXT analysis model, the method has the most prominent advantage of taking into account the mode coupling coefficient varying with the distance, which further ensures the calculation accuracy and eliminates the redundancy caused by polarization mode coupling, thereby getting closer to actual optical fiber laying conditions and achieving a wider application range.

Description

一种多芯少模光纤模间串扰检测方法及装置A method and device for detecting inter-mode crosstalk in multi-core, few-mode optical fiber 技术领域Technical field
本发明涉及光纤技术领域,尤其是指一种多芯少模光纤模间串扰检测方法、装置及计算机存储介质。The invention relates to the field of optical fiber technology, and in particular, to a method, device and computer storage medium for inter-mode crosstalk detection of multi-core, few-mode optical fiber.
背景技术Background technique
随着互联网数据流量的持续指数型增长,基于单模光纤(Single-Mode Fiber,SMF)传输的信息主干网正在迅速接近其容量极限。过去,SMF传输系统的容量增加是通过利用各种维度,包括偏振和波分复用,以及先进的调制格式和相干传输技术来实现的。然而,即将到来的容量紧缩意味着我们需要新的传输技术以支持社会容量需求。为了缓解使用额外SMF带来的与线性容量扩展相关的相应成本和增加的能源需求,单光纤内的空分复用(Space Division Multiplexing,SDM)可以提供一种解决方案。通过引入额外的正交复用维度,可以降低容量、频谱和能量效率,从而降低每比特的成本,这对于维持主要网络利益相关者的商业模式至关重要。为了实现可持续发展管理的承诺,设想了一种新的范式,这使得SMF的容量增加了两个数量级。SDM通过多输入多输出(Multiple-input Multiple-output,MIMO)传输实现,使用多模光纤(Multi-mode Fiber,MMF)的空间模式或多个单模纤芯作为信道。近些年,有研究开发了一种独特的MMF,即少模光纤(Few-mode Fiber,FMF),用于共传输三个或六个线性极化(LP)模式。在高速电子技术快速增强的驱动下,数字信号处理(Digital Signal Processing,DSP)MIMO技术可以快速地恢复混合传输信道,从而在空间信道占用相同波长时提高频谱效率。最先进的单载波FMF传输实验表明,通过利用六种空间模式,单光纤中的容量增加达到32bit/(s·Hz)的频谱效率。如果使用多芯传输,有研究使用12个单模纤芯实现了频谱效率为109bit/(s·Hz)的传输。在另一项工作中,有研究结合了多芯和少模的特性,使用7芯3模少模多芯光纤 (Few-mode Multicore Fiber,FM-MCF),并在1km的传输中实现了255Tbit/s的最高速率。As Internet data traffic continues to grow exponentially, the information backbone network based on single-mode fiber (Single-Mode Fiber, SMF) transmission is rapidly approaching its capacity limit. In the past, capacity increases in SMF transmission systems have been achieved by leveraging various dimensions, including polarization and wavelength division multiplexing, as well as advanced modulation formats and coherent transmission techniques. However, the coming capacity crunch means we need new transmission technologies to support society's capacity needs. To alleviate the consequential costs and increased energy requirements associated with linear capacity expansion using additional SMFs, Space Division Multiplexing (SDM) within a single fiber can provide a solution. By introducing additional orthogonal multiplexing dimensions, capacity, spectrum and energy efficiency can be reduced, thereby reducing cost per bit, which is critical to sustaining the business models of key network stakeholders. To realize the promise of sustainable development management, a new paradigm is envisioned, which increases the capacity of SMFs by two orders of magnitude. SDM is implemented through Multiple-input Multiple-output (MIMO) transmission, using the spatial mode of Multi-mode Fiber (MMF) or multiple single-mode fiber cores as channels. In recent years, research has developed a unique MMF, namely few-mode fiber (Few-mode Fiber, FMF), which is used to transmit a total of three or six linear polarization (LP) modes. Driven by the rapid enhancement of high-speed electronic technology, Digital Signal Processing (DSP) MIMO technology can quickly restore mixed transmission channels, thereby improving spectral efficiency when spatial channels occupy the same wavelength. State-of-the-art single-carrier FMF transmission experiments show that by utilizing six spatial modes, the capacity in a single fiber is increased to a spectral efficiency of 32bit/(s·Hz). If multi-core transmission is used, some research has used 12 single-mode fiber cores to achieve a transmission with a spectral efficiency of 109 bit/(s·Hz). In another work, there is research that combines the characteristics of multi-core and few-mode, using 7-core 3-mode few-mode multi-core fiber (Few-mode Multicore Fiber, FM-MCF), and achieving 255Tbit in 1km transmission /s maximum rate.
因此,FM-MCF成为了近年来研究的热点话题。同时,尽管FM-MCF可以带来超高的传输容量,但其依旧存在两个会影响传输性能的重要问题。其一是由于不同模式在纤芯中传输的路径不同,这会导致在接收端接收到的每路信号存在时间差,最终产生误码,这被称为模式色散。其二是由于不同路径之间会发生模式之间的耦合,从而降低传输性能。目前有大量研究表明模式色散可以通过DSP技术或是引入色散补偿光纤进行缓解,而对于模式耦合的评估还并不完善。因此,准确地评估模式耦合有助于多芯少模光纤的制备。Therefore, FM-MCF has become a hot topic of research in recent years. At the same time, although FM-MCF can bring ultra-high transmission capacity, it still has two important problems that will affect transmission performance. One is that different modes have different transmission paths in the fiber core, which will lead to time differences in each signal received at the receiving end, eventually resulting in bit errors. This is called modal dispersion. The second reason is that coupling between modes occurs between different paths, thereby reducing transmission performance. There are currently a large number of studies showing that modal dispersion can be alleviated through DSP technology or the introduction of dispersion-compensating optical fibers, but the evaluation of mode coupling is not yet complete. Therefore, accurate evaluation of mode coupling facilitates the preparation of multi-core few-mode fibers.
在Monte Carlo模型中,作者采用了带有偏振模耦合影响的模式电场,这进一步使得计算复杂度提升。并且在对于耦合系数的处理上,看作是一个与传输距离无关的常数,这与实际传输情况并不符合。此外,在ISMXT模型所提出的理论中,证明了偏振模耦合带来的对于超模间的串扰可以忽略不计,因此在计算超模间的串扰时,Monte Carlo模型中的理论可以得到进一步的简化。在ISMXT解析模型中,作者只考虑了纵向扰动带来的影响,消除了一部分冗余项,直接从模式功率出发计算模间串扰,从最终的结果来看,Monte Carlo模型中的偏振模耦合带来的影响确实可以忽略不计。但是ISMXT解析模型中把最重要的模式耦合系数做了近似,由一般经验值计算得到的模间串扰显然是不够严谨的。In the Monte Carlo model, the author uses a mode electric field with polarization mode coupling effects, which further increases the computational complexity. And when dealing with the coupling coefficient, it is regarded as a constant that has nothing to do with the transmission distance, which is not consistent with the actual transmission situation. In addition, in the theory proposed by the ISMXT model, it is proved that the crosstalk between supermodes caused by polarization mode coupling can be ignored. Therefore, when calculating the crosstalk between supermodes, the theory in the Monte Carlo model can be further simplified. . In the ISMXT analytical model, the author only considered the impact of longitudinal perturbation, eliminated some redundant terms, and calculated inter-mode crosstalk directly from the mode power. From the final result, the polarization mode coupling band in the Monte Carlo model The impact is indeed negligible. However, the ISMXT analytical model approximates the most important mode coupling coefficient, and the inter-mode crosstalk calculated from general empirical values is obviously not rigorous enough.
发明内容Contents of the invention
为此,本发明所要解决的技术问题在于克服现有技术中模式耦合系数的处理与实际情况不符的问题。To this end, the technical problem to be solved by the present invention is to overcome the problem in the prior art that the processing of mode coupling coefficients is inconsistent with the actual situation.
为解决上述技术问题,本发明提供了一种多芯少模光纤模间串扰检测方法,包括:In order to solve the above technical problems, the present invention provides a multi-core few-mode optical fiber inter-mode crosstalk detection method, which includes:
获取光纤参数和注入功率;Obtain fiber parameters and injection power;
根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;Calculate the electric field intensity within the core range and the electric field intensity within the cladding range according to the optical fiber parameters;
根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;Calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range and the injected power;
将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。The inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory is replaced with the inter-mode coupling mode coefficient, and the optical fiber inter-mode crosstalk value is calculated.
优选地,所述光纤参数包括:Preferably, the fiber parameters include:
纤芯折射率、包层折射率、纤芯半径、纤芯传播常数、芯间距和光波长。Core refractive index, cladding refractive index, core radius, core propagation constant, core spacing and optical wavelength.
优选地,所述根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度包括:Preferably, the calculation of the electric field intensity within the core range and the electric field intensity within the cladding range based on the optical fiber parameters includes:
利用亥姆霍兹方程计算纤芯范围内和包层范围内的纵向电场的幅值:Use the Helmholtz equation to calculate the amplitude of the longitudinal electric field within the fiber core and within the cladding:
Figure PCTCN2022122044-appb-000001
Figure PCTCN2022122044-appb-000001
其中,E z1和E z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数,
Figure PCTCN2022122044-appb-000002
n 1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J f(x)表示f阶的第一类贝塞尔函数,K f(x)表示f阶的第二类修正贝塞尔函数,
Figure PCTCN2022122044-appb-000003
n 2是包层折射率,r表示柱坐标系中的半径方向,θ表示柱坐标系中的圆周方向;
Among them, E z1 and E z2 respectively represent the amplitude of the longitudinal electric field within the fiber core and the amplitude of the longitudinal electric field within the cladding. A=jUC/βa, j is an imaginary unit, and C is a system constant.
Figure PCTCN2022122044-appb-000002
n 1 is the core refractive index, β is the core propagation constant, a is the core radius, k = 2π/λ is the wave number, λ represents the light wavelength, J f (x) represents the first kind of Bessel function of order f , K f (x) represents the second kind of modified Bessel function of order f,
Figure PCTCN2022122044-appb-000003
n 2 is the refractive index of the cladding, r represents the radial direction in the cylindrical coordinate system, and θ represents the circumferential direction in the cylindrical coordinate system;
将所述纤芯范围内和包层范围内的纵向电场的幅值代入麦克斯韦方程,得到纤芯范围内和包层范围内的横向电场的幅值:Substituting the amplitude of the longitudinal electric field within the core range and the cladding range into Maxwell's equation, the amplitude of the transverse electric field within the core range and the cladding range is obtained:
Figure PCTCN2022122044-appb-000004
Figure PCTCN2022122044-appb-000004
其中,E r1和E r2分别表示纤芯范围内的横向电场的幅值和包层范围内的横向电场的幅值,J' f(x)表示J f(x)的一阶导数; Among them, E r1 and E r2 respectively represent the amplitude of the transverse electric field within the core range and the amplitude of the transverse electric field within the cladding range, and J' f (x) represents the first derivative of J f (x);
根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度:
Figure PCTCN2022122044-appb-000005
According to the amplitude of the longitudinal electric field within the core and the amplitude of the transverse electric field within the core, the electric field strength within the core is calculated:
Figure PCTCN2022122044-appb-000005
根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度:
Figure PCTCN2022122044-appb-000006
Calculate the electric field strength within the cladding range based on the amplitude of the longitudinal electric field within the cladding range and the amplitude of the transverse electric field within the cladding range:
Figure PCTCN2022122044-appb-000006
其中,E m和E n分别表示在纤芯范围内和包层范围内的电场强度,电场圆 周方向幅值
Figure PCTCN2022122044-appb-000007
e r表示r方向上的单位向量,
Figure PCTCN2022122044-appb-000008
表示
Figure PCTCN2022122044-appb-000009
方向上的单位向量,e z表示z方向上的单位向量。
Among them, Em and En represent the electric field intensity within the core range and the cladding range respectively, and the amplitude of the electric field in the circumferential direction
Figure PCTCN2022122044-appb-000007
e r represents the unit vector in the r direction,
Figure PCTCN2022122044-appb-000008
express
Figure PCTCN2022122044-appb-000009
The unit vector in the direction, e z represents the unit vector in the z direction.
优选地,所述根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数前包括:Preferably, before calculating the inter-mode coupling mode coefficient according to the electric field intensity within the core range, the electric field intensity within the cladding range and the injection power, the method includes:
将所述包层范围内的电场强度中,单芯横向极坐标系(r,θ)扩展为双芯的横向极坐标系(R,Θ),
Figure PCTCN2022122044-appb-000010
D表示芯间距。
In the electric field intensity within the cladding range, the single-core transverse polar coordinate system (r, θ) is expanded to the dual-core transverse polar coordinate system (R, Θ),
Figure PCTCN2022122044-appb-000010
D represents the core distance.
优选地,所述根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及注入功率计算模式间耦合模系数:Preferably, the inter-mode coupling mode coefficient is calculated based on the electric field intensity within the core range, the electric field intensity within the cladding range and the injection power:
计算纤芯范围内的电场强度、包层范围内的电场强度和纤芯折射率和包层折射率之差点积后的二重积分与角频率以及真空介电常数的乘积,再除以4倍的注入功率,得到模式间耦合模系数:Calculate the double integral of the electric field strength within the fiber core, the electric field strength within the cladding and the difference between the core refractive index and the cladding refractive index, the angular frequency and the vacuum dielectric constant, and then divide by 4 times Injection power, the inter-mode coupling mode coefficient is obtained:
Figure PCTCN2022122044-appb-000011
Figure PCTCN2022122044-appb-000011
其中,ω是角频率,ε 0是真空介电常数,E m和H m分别代表在纤芯范围内的电场强度和磁场强度,而E n代表在包层范围内的电场强度,e z表示z方向上的单位向量,n 1和n 2分别表示纤芯的折射率和包层的折射率,符号*表示对矩阵进行共轭操作,P表示注入的功率, Among them, ω is the angular frequency, ε 0 is the vacuum dielectric constant, E m and H m represent the electric field intensity and magnetic field intensity in the core range respectively, and E n represents the electric field intensity in the cladding range, and e z represents The unit vector in the z direction, n 1 and n 2 represent the refractive index of the fiber core and the refractive index of the cladding respectively, the symbol * represents the conjugate operation of the matrix, P represents the injected power,
Figure PCTCN2022122044-appb-000012
Figure PCTCN2022122044-appb-000012
优选地,所述基于耦合模理论的芯间串扰计算方程为:Preferably, the inter-core crosstalk calculation equation based on coupled mode theory is:
Figure PCTCN2022122044-appb-000013
Figure PCTCN2022122044-appb-000013
其中,N为将光纤分段后的段数,d是将光纤分段后的段长,Δβ eq,mn,i是第i段的等效相位失配,
Figure PCTCN2022122044-appb-000014
为第i段修正后的纤芯间耦合模系数,k i(d)为纤芯间的耦合模系数。
Among them, N is the number of segments after segmenting the optical fiber, d is the segment length after segmenting the optical fiber, Δβ eq,mn,i is the equivalent phase mismatch of the i-th segment,
Figure PCTCN2022122044-appb-000014
is the modified coupling mode coefficient between fiber cores in the i-th section, k i (d) is the coupling mode coefficient between fiber cores.
优选地,所述将纤芯间耦合模系数k i(d)替换为所述模式间耦合模系数k mn,i,计算得到光纤模间串扰值: Preferably, the inter-core coupling mode coefficient k i (d) is replaced by the inter-mode coupling mode coefficient k mn,i , and the optical fiber inter-mode crosstalk value is calculated:
Figure PCTCN2022122044-appb-000015
Figure PCTCN2022122044-appb-000015
其中,g mn,i为第i段修正后的模式间耦合模系数。 Among them, g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
本发明还提供了一种多芯少模光纤模间串扰检测装置,包括:The invention also provides a multi-core few-mode optical fiber inter-mode crosstalk detection device, which includes:
光纤参数和注入功率获取模块,用于获取光纤参数和注入功率;Fiber parameter and injection power acquisition module, used to obtain fiber parameters and injection power;
电场强度计算模块,用于根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;An electric field strength calculation module, used to calculate the electric field strength within the core range and the electric field strength within the cladding range based on the optical fiber parameters;
模式间耦合模系数计算模块,用于根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;An inter-mode coupling mode coefficient calculation module, used to calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
光纤模间串扰值计算模块,用于将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。An optical fiber inter-mode crosstalk value calculation module is used to replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
优选地,所述多芯少模光纤模间串扰检测装置应用于异质或同质光纤。Preferably, the multi-core few-mode optical fiber inter-mode crosstalk detection device is applied to heterogeneous or homogeneous optical fibers.
本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质上 存储有计算机程序,所述计算机程序被处理器执行时实现上述一种多芯少模光纤模间串扰检测方法的步骤。The present invention also provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, the above-mentioned method for detecting inter-mode crosstalk in multi-core, low-mode optical fiber is implemented. step.
本发明的上述技术方案相比现有技术具有以下优点:The above technical solution of the present invention has the following advantages compared with the existing technology:
本发明所述的多芯少模光纤模间串扰检测方法,首先根据光纤参数计算每个模式的电场分布,然后根据电场分布和注入功率计算模式间耦合系数,最后基于耦合模理论,重新推导了超模之间串扰量的计算方式,相比于Monte Carlo模型和ISMXT解析模型,该方法最突出的优势就是考虑了随距离变化的模式耦合系数,这进一步保证了计算的准确性,且剔除了偏振模耦合带来的冗余,更加贴合实际光纤铺设情况,并且应用范围更广。The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to the present invention first calculates the electric field distribution of each mode according to the fiber parameters, then calculates the inter-mode coupling coefficient according to the electric field distribution and injection power, and finally re-derives based on the coupled mode theory Compared with the Monte Carlo model and the ISMXT analytical model, the most prominent advantage of this method for calculating the amount of crosstalk between supermodes is that it takes into account the mode coupling coefficient that changes with distance, which further ensures the accuracy of the calculation and eliminates The redundancy brought by polarization mode coupling is more in line with the actual fiber laying situation and has a wider application range.
附图说明Description of drawings
为了使本发明的内容更容易被清楚的理解,下面根据本发明的具体实施例并结合附图,对本发明作进一步详细的说明,其中:In order to make the content of the present invention easier to understand clearly, the present invention will be further described in detail below based on specific embodiments of the present invention and in conjunction with the accompanying drawings, wherein:
图1是本发明多芯少模光纤模间串扰检测方法的实现流程图;Figure 1 is a flow chart for the implementation of the inter-mode crosstalk detection method for multi-core few-mode optical fibers according to the present invention;
图2是模间串扰的纵向变化示意图;Figure 2 is a schematic diagram of the longitudinal change of inter-mode crosstalk;
图3是双芯极坐标系示意图;Figure 3 is a schematic diagram of the dual-core polar coordinate system;
图4是基于耦合模理论的串扰估计模型原理图;Figure 4 is a schematic diagram of the crosstalk estimation model based on coupled mode theory;
图5是串扰对比流程图;Figure 5 is a crosstalk comparison flow chart;
图6是LP01、LP11a和LP11b三种模式的横向电场分布;Figure 6 shows the lateral electric field distribution of three modes: LP01, LP11a and LP11b;
图7是双芯光纤的横向结构示意图;Figure 7 is a schematic diagram of the transverse structure of a dual-core optical fiber;
图8是非线性串扰随功率变化示意图;Figure 8 is a schematic diagram of nonlinear crosstalk changing with power;
图9是模式耦合系数随芯间距变化示意图;Figure 9 is a schematic diagram of the change of mode coupling coefficient with core spacing;
图10是模间耦合系数随芯径比变化示意图;Figure 10 is a schematic diagram of the inter-mode coupling coefficient changing with the core diameter ratio;
图11是模间耦合系数随折射率差变化示意图;Figure 11 is a schematic diagram of the inter-mode coupling coefficient changing with the refractive index difference;
图12为本发明实施例提供的一种多芯少模光纤模间串扰检测装置的结构框图。Figure 12 is a structural block diagram of a multi-core few-mode optical fiber inter-mode crosstalk detection device provided by an embodiment of the present invention.
具体实施方式Detailed ways
本发明的核心是提供一种多芯少模光纤模间串扰检测方法、装置及计算机存储介质,考虑了随距离变化的模式耦合系数,更加贴合实际光纤铺设情况。The core of the present invention is to provide a method, device and computer storage medium for inter-mode crosstalk detection of multi-core, few-mode optical fibers, which take into account the mode coupling coefficient changing with distance and are more in line with the actual optical fiber laying situation.
为了使本技术领域的人员更好地理解本发明方案,下面结合附图和具体实施方式对本发明作进一步的详细说明。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to enable those skilled in the art to better understand the solution of the present invention, the present invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. Obviously, the described embodiments are only some of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present invention.
请参考图1,图1为本发明所提供的多芯少模光纤模间串扰检测方法的实现流程图;具体操作步骤如下:Please refer to Figure 1, which is a flow chart of the implementation of the inter-mode crosstalk detection method for multi-core few-mode optical fibers provided by the present invention; the specific operation steps are as follows:
阶跃型多芯少模光纤的原理是改变光入射纤芯的角度,从而使每个模式在单个纤芯中按照不同的折线轨迹进行传输,其具体的原理图如图2所示;The principle of step-type multi-core few-mode optical fiber is to change the angle of light incident on the fiber core, so that each mode is transmitted along a different fold line trajectory in a single fiber core. The specific schematic diagram is shown in Figure 2;
由于模式之间的串扰发生会非常频繁,我们可以采用分段思想,假设在长度较小的一段里,串扰只会发生一次,并且在每一小段中模式耦合系数可以看作是一个常数。Since crosstalk between modes will occur very frequently, we can adopt the segmentation idea, assuming that crosstalk will only occur once in a small segment, and the mode coupling coefficient can be regarded as a constant in each small segment.
S101:获取光纤参数和注入功率;S101: Obtain fiber parameters and injection power;
所述光纤参数包括纤芯折射率、包层折射率、纤芯半径、纤芯传播常数、芯间距和光波长。The optical fiber parameters include core refractive index, cladding refractive index, core radius, core propagation constant, core spacing and optical wavelength.
S102:根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;S102: Calculate the electric field intensity within the core range and the electric field intensity within the cladding range according to the optical fiber parameters;
利用亥姆霍兹方程计算纤芯范围内和包层范围内的纵向电场的幅值:Use the Helmholtz equation to calculate the amplitude of the longitudinal electric field within the fiber core and within the cladding:
Figure PCTCN2022122044-appb-000016
Figure PCTCN2022122044-appb-000016
其中,E z1和E z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数,
Figure PCTCN2022122044-appb-000017
n 1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J f(x)表示f阶的第一类贝塞尔函数,K f(x)表示f阶的第二类修正贝塞尔函数,
Figure PCTCN2022122044-appb-000018
n 2是包层折射率,r表示柱坐标系中的半径方向,θ表示柱坐标系中的圆周方向;
Among them, E z1 and E z2 respectively represent the amplitude of the longitudinal electric field within the fiber core and the amplitude of the longitudinal electric field within the cladding. A=jUC/βa, j is an imaginary unit, and C is a system constant.
Figure PCTCN2022122044-appb-000017
n 1 is the core refractive index, β is the core propagation constant, a is the core radius, k = 2π/λ is the wave number, λ represents the light wavelength, J f (x) represents the first kind of Bessel function of order f , K f (x) represents the second kind of modified Bessel function of order f,
Figure PCTCN2022122044-appb-000018
n 2 is the refractive index of the cladding, r represents the radial direction in the cylindrical coordinate system, and θ represents the circumferential direction in the cylindrical coordinate system;
将所述纤芯范围内和包层范围内的纵向电场的幅值代入麦克斯韦方程,得到纤芯范围内和包层范围内的横向电场的幅值:Substituting the amplitude of the longitudinal electric field within the core range and the cladding range into Maxwell's equation, the amplitude of the transverse electric field within the core range and the cladding range is obtained:
Figure PCTCN2022122044-appb-000019
Figure PCTCN2022122044-appb-000019
其中,E r1和E r2分别表示纤芯范围内的横向电场的幅值和包层范围内的横向电场的幅值,J' f(x)表示J f(x)的一阶导数; Among them, E r1 and E r2 respectively represent the amplitude of the transverse electric field within the core range and the amplitude of the transverse electric field within the cladding range, and J' f (x) represents the first derivative of J f (x);
由于光纤是一个圆柱体结构,可以将模式电场产开成圆柱极坐标形式,这样电场向量就可以表示成三个维度的叠加:Since the optical fiber is a cylindrical structure, the mode electric field can be generated into a cylindrical polar coordinate form, so that the electric field vector can be expressed as a superposition of three dimensions:
根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度:
Figure PCTCN2022122044-appb-000020
Calculate the electric field strength within the fiber core according to the amplitude of the longitudinal electric field within the fiber core and the amplitude of the transverse electric field within the fiber core:
Figure PCTCN2022122044-appb-000020
根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度:
Figure PCTCN2022122044-appb-000021
Calculate the electric field strength within the cladding range based on the amplitude of the longitudinal electric field within the cladding range and the amplitude of the transverse electric field within the cladding range:
Figure PCTCN2022122044-appb-000021
其中,E m和E n分别表示在纤芯范围内和包层范围内的电场强度,电场圆周方向幅值
Figure PCTCN2022122044-appb-000022
e r表示r方向上的单位向量,
Figure PCTCN2022122044-appb-000023
表示
Figure PCTCN2022122044-appb-000024
方向上的单位向量,e z表示z方向上的单位向量。
Among them, Em and En represent the electric field intensity within the core range and the cladding range respectively, and the amplitude of the electric field in the circumferential direction
Figure PCTCN2022122044-appb-000022
e r represents the unit vector in the r direction,
Figure PCTCN2022122044-appb-000023
express
Figure PCTCN2022122044-appb-000024
The unit vector in the direction, e z represents the unit vector in the z direction.
S103:根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;S103: Calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
模式耦合系数与每个模式的电场分布密切相关;The mode coupling coefficient is closely related to the electric field distribution of each mode;
如图3,为考虑两根纤芯之间的模式耦合,将所述包层范围内的电场强度中,单芯横向极坐标系(r,θ)扩展为双芯的横向极坐标系(R,Θ),
Figure PCTCN2022122044-appb-000025
D表示芯间距;
As shown in Figure 3, in order to consider the mode coupling between the two cores, the single-core transverse polar coordinate system (r, θ) in the electric field intensity within the cladding range is expanded to the dual-core transverse polar coordinate system (R ,Θ),
Figure PCTCN2022122044-appb-000025
D represents the core spacing;
计算纤芯范围内的电场强度、包层范围内的电场强度和纤芯折射率和包 层折射率之差点积后的二重积分与角频率以及真空介电常数的乘积,再除以4倍的注入功率,得到模式间耦合模系数:Calculate the double integral of the electric field strength within the fiber core, the electric field strength within the cladding and the difference between the core refractive index and the cladding refractive index, the angular frequency and the vacuum dielectric constant, and then divide by 4 times Injection power, the inter-mode coupling mode coefficient is obtained:
Figure PCTCN2022122044-appb-000026
Figure PCTCN2022122044-appb-000026
其中,ω是角频率,ε 0是真空介电常数,E m和H m分别代表在纤芯范围内的电场强度和磁场强度,而E n代表在包层范围内的电场强度,e z表示z方向上的单位向量,n 1和n 2分别表示纤芯的折射率和包层的折射率,符号*表示对矩阵进行共轭操作,P表示注入的功率, Among them, ω is the angular frequency, ε 0 is the vacuum dielectric constant, E m and H m represent the electric field intensity and magnetic field intensity in the core range respectively, and E n represents the electric field intensity in the cladding range, and e z represents The unit vector in the z direction, n 1 and n 2 represent the refractive index of the fiber core and the refractive index of the cladding respectively, the symbol * represents the conjugate operation of the matrix, P represents the injected power,
Figure PCTCN2022122044-appb-000027
用数值求解的方式得到最终的结果。
Figure PCTCN2022122044-appb-000027
The final result is obtained by numerical solution.
S104:将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。S104: Replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
本发明所述的多芯少模光纤模间串扰检测方法,首先根据光纤参数计算每个模式的电场分布,然后根据电场分布和注入功率计算模式间耦合系数,最后基于耦合模理论,重新推导了超模之间串扰量的计算方式,相比于Monte Carlo模型和ISMXT解析模型,该方法最突出的优势就是考虑了随距离变化的模式耦合系数,这进一步保证了计算的准确性,且剔除了偏振模耦合带来的冗余,更加贴合实际光纤铺设情况,并且应用范围更广。在此基础上,我们可以研究不同通信系统中的串扰特性,根据串扰在不同工作区间的传输特性、光纤本身的物理几何特性和串扰与光纤参数的关系,进而研究降低不同传输模式之间串扰的理论方法。The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to the present invention first calculates the electric field distribution of each mode according to the fiber parameters, then calculates the inter-mode coupling coefficient according to the electric field distribution and injection power, and finally re-derives based on the coupled mode theory Compared with the Monte Carlo model and the ISMXT analytical model, the most prominent advantage of this method for calculating the amount of crosstalk between supermodes is that it takes into account the mode coupling coefficient that changes with distance, which further ensures the accuracy of the calculation and eliminates The redundancy brought by polarization mode coupling is more in line with the actual fiber laying situation and has a wider application range. On this basis, we can study the crosstalk characteristics in different communication systems. Based on the transmission characteristics of crosstalk in different working areas, the physical geometric characteristics of the fiber itself and the relationship between crosstalk and fiber parameters, we can then study ways to reduce crosstalk between different transmission modes. Theoretical methods.
基于以上实施例,本实施例对步骤S104进行进一步说明:Based on the above embodiment, this embodiment further explains step S104:
本发明是基于耦合模理论来计算模间串扰,首先需要介绍耦合模理论的原理以及重要的求解方法:This invention calculates inter-mode crosstalk based on coupled mode theory. First, it is necessary to introduce the principles of coupled mode theory and important solution methods:
耦合模方程为:
Figure PCTCN2022122044-appb-000028
The coupled mode equation is:
Figure PCTCN2022122044-appb-000028
其中,A n,A m分别为芯n和芯m缓慢变化的电场复振幅,N为纤芯数,K mn是芯n到芯m耦合系数。β m和β n分别为m和n芯的传播常数,j是虚数单位,z表示纵向传播方向。为了基于耦合模理论求解线性串扰,提出了一种通用半解析数学模型,其根据相位匹配点(Phase Machining Point,PMP)的分布来精确求解线性耦合模方程,其原理如图4所示,其中d是将光纤分段后的段长,P 0表示注入功率,P n和P m分别表示n芯和m芯的输出功率。基于以上的分段模型,我们可以将式(1)看作两个分量形式: Among them, A n and A m are the slowly changing electric field complex amplitudes of core n and core m respectively, N is the number of fiber cores, and K mn is the coupling coefficient from core n to core m. β m and β n are the propagation constants of m and n cores respectively, j is the imaginary unit, and z represents the longitudinal propagation direction. In order to solve linear crosstalk based on coupled mode theory, a general semi-analytic mathematical model is proposed, which accurately solves the linear coupled mode equation according to the distribution of phase matching points (Phase Machining Points, PMP). The principle is shown in Figure 4, where d is the segment length after segmenting the optical fiber, P 0 represents the injection power, P n and P m represent the output power of n core and m core respectively. Based on the above segmented model, we can regard equation (1) as two component forms:
Figure PCTCN2022122044-appb-000029
Figure PCTCN2022122044-appb-000029
其中Δβ eq,mn,i是第i段的等效相位失配,可表示为Δβ eq,mn,i=β eq,m,ieq,n,i;K mn,i(z)为纤芯m和纤芯n之间的第i段耦合模系数,当分段段长足够小的时候,K mn,i(z)在第i段可以看作是常数K mn,i。当光功率在纤芯n中进行传输时,通过求解式(2)可以得到纤芯m和n在每一段末尾处的电场解析解为: where Δβ eq,mn,i is the equivalent phase mismatch of the i-th segment, which can be expressed as Δβ eq,mn,ieq,m,ieq,n,i ; K mn,i (z) is The coupling mode coefficient of the i-th section between core m and core n. When the segment length is small enough, K mn,i (z) can be regarded as a constant K mn,i in the i-th section. When the optical power is transmitted in the fiber core n, the analytical solution of the electric field of the fiber core m and n at the end of each section can be obtained by solving equation (2) as:
Figure PCTCN2022122044-appb-000030
Figure PCTCN2022122044-appb-000030
其中T是解的系数矩阵,可以表示为:where T is the coefficient matrix of the solution, which can be expressed as:
Figure PCTCN2022122044-appb-000031
Figure PCTCN2022122044-appb-000031
Figure PCTCN2022122044-appb-000032
Figure PCTCN2022122044-appb-000032
Figure PCTCN2022122044-appb-000033
Figure PCTCN2022122044-appb-000033
Figure PCTCN2022122044-appb-000034
Figure PCTCN2022122044-appb-000034
式中
Figure PCTCN2022122044-appb-000035
为第i段修正后的耦合模系数,并且假设纤芯间的耦合模系数k i(d)=K mn,i(d)=K nm,i(d)。如此,每一段电场幅值都可以通过前一段的值得到,同时也可以得到每段串扰的变化,将其叠加就可以计算出总串扰值:
in the formula
Figure PCTCN2022122044-appb-000035
is the modified coupling mode coefficient of the i-th section, and it is assumed that the coupling mode coefficient between fiber cores k i (d) = K mn,i (d) = K nm,i (d). In this way, the electric field amplitude of each segment can be obtained from the value of the previous segment, and the change of crosstalk in each segment can also be obtained. By superimposing them, the total crosstalk value can be calculated:
Figure PCTCN2022122044-appb-000036
Figure PCTCN2022122044-appb-000036
我们只需要将其中的k i(d)替换为模式之间的耦合系数k mn,i,就可以将式(5)更改为模式之间的串扰值公式: We only need to replace k i (d) with the coupling coefficient k mn,i between modes, and then equation (5) can be changed into the crosstalk value formula between modes:
Figure PCTCN2022122044-appb-000037
Figure PCTCN2022122044-appb-000037
其中,g mn,i为第i段修正后的模式间耦合模系数。 Among them, g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
本发明可以为多芯少模光纤通信系统提供一个快速、准确的串扰估计计算方法,其应用范围更加广泛。该方法不仅适用于异质光纤情况,同样适用于真实的同质光纤和理想的完全同质光纤。该方法基于能体现光纤传输物理特性的耦合模理论,并且针对纤芯中存在的随机扰动重新推导了模间耦合系数,因此相较于以往模型能更准确地反应模间串扰的变化过程。基于该模型,我们可以对少模多芯光纤中的模间串扰有更前瞻性的认识。The invention can provide a fast and accurate crosstalk estimation calculation method for multi-core and few-mode optical fiber communication systems, and its application range is wider. This method is not only applicable to heterogeneous optical fibers, but also to real homogeneous optical fibers and ideal completely homogeneous optical fibers. This method is based on the coupled mode theory that reflects the physical characteristics of optical fiber transmission, and re-derives the inter-mode coupling coefficient based on the random disturbances present in the fiber core. Therefore, it can more accurately reflect the changing process of inter-mode crosstalk than previous models. Based on this model, we can have a more forward-looking understanding of inter-mode crosstalk in few-mode multi-core optical fibers.
如图5,基于以上实施例,本实施例通过仿真验证上述方法的准确性,具体如下:As shown in Figure 5, based on the above embodiment, this embodiment verifies the accuracy of the above method through simulation, as follows:
LP01、LP11a和LP11b三种模式的横向电场分布如图6所示,其中LP01是传输基模,LP11a和LP11b是两个相互正交的偏振模。图7展示两根纤芯的横向排布;The transverse electric field distributions of the three modes LP01, LP11a and LP11b are shown in Figure 6, where LP01 is the transmission fundamental mode, and LP11a and LP11b are two mutually orthogonal polarization modes. Figure 7 shows the lateral arrangement of the two fiber cores;
考虑一个纤芯半径a=2.5μm,包层折射率n 2=1.45,纤芯折射率约为n 1=1.4545,芯间距为D nm=30μm,光脉冲波长为1550nm的双芯三模光纤,首先对比传输距离为z=10m情况下三种串扰计算方法的区别,比较结果如图8所示,图8给出了ISMXT随着传输距离增加的变化情况。从整体的变化趋势看,本发明与MonteCarlo模型和ISMXT解析模型基本保持一致,随着传输距离逐渐变长,本发明逐渐与另外两者保持3dB左右的差距,这主要是因为随机扰动因素带来的影响开始体现,而另外两种方法参与计算的模式耦合系数是一个近似的经验值,并不能准确估计随机扰动带来的影响。本发明的模间串扰结果存在波动性上升的情况,这是因为各个模式之间的相位匹配点会随着距离进行随机分布,存在波动的那个点就表示出现了相位匹配点。 Consider a dual-core three-mode fiber with core radius a = 2.5 μm, cladding refractive index n 2 = 1.45, core refractive index approximately n 1 = 1.4545, core spacing D nm = 30 μm, and optical pulse wavelength 1550 nm. First, compare the differences between the three crosstalk calculation methods when the transmission distance is z=10m. The comparison results are shown in Figure 8. Figure 8 shows the changes in ISMXT as the transmission distance increases. From the overall change trend, the present invention is basically consistent with the MonteCarlo model and the ISMXT analytical model. As the transmission distance gradually becomes longer, the present invention gradually maintains a gap of about 3dB with the other two. This is mainly due to random disturbance factors. The influence of the model begins to be reflected, and the mode coupling coefficient calculated by the other two methods is an approximate empirical value and cannot accurately estimate the impact of random disturbances. The inter-mode crosstalk results of the present invention have increased volatility. This is because the phase matching points between each mode are randomly distributed with distance, and the point with fluctuation indicates the occurrence of a phase matching point.
图9展示了随芯间距变化的过程,纤芯的半径a=2.5μm。当芯间距逐渐变大时,模间耦合系数逐渐减小,这说明适当增加芯间距是有效降低模间串扰的方法之一。此外,图9分别考虑了本征有效折射率差为0.21%、0.31%和0.41%三种情况,可以发现随着折射率差的变大,模间耦合系数也会逐渐减小,并且在芯间距较小的情况下,减小地更明显。Figure 9 shows the process as the core spacing changes. The core radius a = 2.5 μm. As the core spacing gradually increases, the inter-mode coupling coefficient gradually decreases, which shows that appropriately increasing the core spacing is one of the methods to effectively reduce inter-mode crosstalk. In addition, Figure 9 considers three cases where the intrinsic effective refractive index difference is 0.21%, 0.31% and 0.41% respectively. It can be found that as the refractive index difference becomes larger, the inter-mode coupling coefficient will gradually decrease, and in the core When the spacing is small, the reduction is more obvious.
我们还根据推导出的模型做了其他物理特性的仿真。图10和图11分别显示了与芯径比(芯间距与纤芯半径的比值)和本征有效折射率差的变化。如图10所示,从整体上看芯径比的增大会使得模间耦合系数变小,但是在折射率差较大的情况下,这种影响并不明显,相反在折射率差较小的情况下影响较为明显,这说明起更主要作用的还是折射率差。图11我们将横坐标设置为了折射率差,最终发现,无论芯径比的取值是大是小,折射率差的增大会使得模式耦合系数显著减小。基于以上结论,我们发现芯间距、芯径比和有效折射率差等物理参数会对模式耦合系数产生影响,从而影响模间串扰的变化。 基于以上结论,我们可以通过制作具有合适的物理特性的多芯少模光纤来进一步减小模间串扰。We also performed simulations of other physical properties based on the derived model. Figures 10 and 11 show the changes with core diameter ratio (ratio of core spacing to core radius) and intrinsic effective refractive index difference, respectively. As shown in Figure 10, overall, an increase in the core diameter ratio will make the inter-mode coupling coefficient smaller. However, when the refractive index difference is large, this effect is not obvious. On the contrary, when the refractive index difference is small, The effect is more obvious in this case, which shows that the refractive index difference plays a more important role. In Figure 11, we set the abscissa to the refractive index difference, and finally found that no matter whether the value of the core diameter ratio is large or small, the increase in the refractive index difference will significantly reduce the mode coupling coefficient. Based on the above conclusion, we found that physical parameters such as core spacing, core diameter ratio and effective refractive index difference will affect the mode coupling coefficient, thereby affecting the changes in inter-mode crosstalk. Based on the above conclusions, we can further reduce inter-mode crosstalk by making multi-core few-mode optical fibers with suitable physical properties.
请参考图12,图12为本发明实施例提供的一种多芯少模光纤模间串扰检测装置的结构框图;具体装置可以包括:Please refer to Figure 12. Figure 12 is a structural block diagram of a multi-core few-mode optical fiber inter-mode crosstalk detection device provided by an embodiment of the present invention; the specific device may include:
光纤参数和注入功率获取模块100,用于获取光纤参数和注入功率;Optical fiber parameter and injection power acquisition module 100, used to acquire optical fiber parameters and injection power;
电场强度计算模块200,用于根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;The electric field strength calculation module 200 is used to calculate the electric field strength within the core range and the electric field strength within the cladding range according to the optical fiber parameters;
模式间耦合模系数计算模块300,用于根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;The inter-mode coupling mode coefficient calculation module 300 is used to calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
光纤模间串扰值计算模块400,用于将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。The optical fiber inter-mode crosstalk value calculation module 400 is used to replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
本实施例的多芯少模光纤模间串扰检测装置用于实现前述的多芯少模光纤模间串扰检测方法,因此多芯少模光纤模间串扰检测装置中的具体实施方式可见前文多芯少模光纤模间串扰检测方法的实施例部分,例如,光纤参数和注入功率获取模块100,电场强度计算模块200,模式间耦合模系数计算模块300,光纤模间串扰值计算模块400,分别用于实现上述多芯少模光纤模间串扰检测方法中步骤S101,S102,S103,S104,所以,其具体实施方式可以参照相应的各个部分实施例的描述,在此不再赘述。The inter-mode crosstalk detection device for multi-core and few-mode optical fiber in this embodiment is used to implement the aforementioned inter-mode crosstalk detection method for multi-core and few-mode optical fiber. Therefore, the specific implementation of the multi-core and few-mode optical fiber inter-mode crosstalk detection device can be found in the multi-core and few-mode optical fiber mentioned above. The embodiment part of the few-mode optical fiber inter-mode crosstalk detection method, for example, the optical fiber parameter and injection power acquisition module 100, the electric field strength calculation module 200, the inter-mode coupling mode coefficient calculation module 300, and the optical fiber inter-mode crosstalk value calculation module 400 are respectively used. In order to implement steps S101, S102, S103, and S104 in the inter-mode crosstalk detection method for multi-core few-mode optical fibers, the specific implementation methods can be referred to the descriptions of the corresponding embodiments, and will not be described again here.
本发明具体实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述一种多芯少模光纤模间串扰检测方法的步骤。Specific embodiments of the present invention also provide a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, the above-mentioned multi-core few-mode optical fiber inter-mode crosstalk is realized. Steps of the detection method.
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘 存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will understand that embodiments of the present application may be provided as methods, systems, or computer program products. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment that combines software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk memory, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each process and/or block in the flowchart illustrations and/or block diagrams, and combinations of processes and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine, such that the instructions executed by the processor of the computer or other programmable data processing device produce a use A device for realizing the functions specified in one process or multiple processes of the flowchart and/or one block or multiple blocks of the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that causes a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction means, the instructions The device implements the functions specified in a process or processes of the flowchart and/or a block or blocks of the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device, causing a series of operating steps to be performed on the computer or other programmable device to produce computer-implemented processing, thereby executing on the computer or other programmable device. Instructions provide steps for implementing the functions specified in a process or processes of a flowchart diagram and/or a block or blocks of a block diagram.
显然,上述实施例仅仅是为清楚地说明所作的举例,并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。Obviously, the above-mentioned embodiments are only examples for clear explanation and are not intended to limit the implementation. For those of ordinary skill in the art, other changes or modifications may be made based on the above description. An exhaustive list of all implementations is neither necessary nor possible. The obvious changes or modifications derived therefrom are still within the protection scope of the present invention.

Claims (10)

  1. 一种多芯少模光纤模间串扰检测方法,其特征在于,包括:A method for detecting inter-mode crosstalk in multi-core, few-mode optical fiber, which is characterized by including:
    获取光纤参数和注入功率;Obtain fiber parameters and injection power;
    根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;Calculate the electric field intensity within the core range and the electric field intensity within the cladding range according to the optical fiber parameters;
    根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;Calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range and the injected power;
    将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。The inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory is replaced with the inter-mode coupling mode coefficient, and the optical fiber inter-mode crosstalk value is calculated.
  2. 根据权利要求1所述的多芯少模光纤模间串扰检测方法,其特征在于,所述光纤参数包括:The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 1, characterized in that the optical fiber parameters include:
    纤芯折射率、包层折射率、纤芯半径、纤芯传播常数、芯间距和光波长。Core refractive index, cladding refractive index, core radius, core propagation constant, core spacing, and wavelength of light.
  3. 根据权利要求1所述的多芯少模光纤模间串扰检测方法,其特征在于,所述根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度包括:The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 1, wherein the calculation of the electric field intensity within the fiber core range and the electric field intensity within the cladding range based on the fiber parameters includes:
    利用亥姆霍兹方程计算纤芯范围内和包层范围内的纵向电场的幅值:
    Figure PCTCN2022122044-appb-100001
    Use the Helmholtz equation to calculate the amplitude of the longitudinal electric field within the fiber core and within the cladding:
    Figure PCTCN2022122044-appb-100001
    其中,E z1和E z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数,
    Figure PCTCN2022122044-appb-100002
    n 1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J f(x)表示f阶的第一类贝塞尔函数,K f(x)表示f阶的第二类修正贝塞尔函数,
    Figure PCTCN2022122044-appb-100003
    n 2是包层折射 率,r表示柱坐标系中的半径方向,θ表示柱坐标系中的圆周方向;
    Among them, E z1 and E z2 respectively represent the amplitude of the longitudinal electric field within the fiber core and the amplitude of the longitudinal electric field within the cladding. A=jUC/βa, j is an imaginary unit, and C is a system constant.
    Figure PCTCN2022122044-appb-100002
    n 1 is the core refractive index, β is the core propagation constant, a is the core radius, k = 2π/λ is the wave number, λ represents the light wavelength, J f (x) represents the first kind of Bessel function of order f , K f (x) represents the second kind of modified Bessel function of order f,
    Figure PCTCN2022122044-appb-100003
    n 2 is the refractive index of the cladding, r represents the radial direction in the cylindrical coordinate system, and θ represents the circumferential direction in the cylindrical coordinate system;
    将所述纤芯范围内和包层范围内的纵向电场的幅值代入麦克斯韦方程,得到纤芯范围内和包层范围内的横向电场的幅值:
    Figure PCTCN2022122044-appb-100004
    Substituting the amplitude of the longitudinal electric field within the core range and the cladding range into Maxwell's equation, the amplitude of the transverse electric field within the core range and the cladding range is obtained:
    Figure PCTCN2022122044-appb-100004
    其中,E r1和E r2分别表示纤芯范围内的横向电场的幅值和包层范围内的横向电场的幅值,J′ f(x)表示J f(x)的一阶导数; Among them, E r1 and E r2 respectively represent the amplitude of the transverse electric field within the core range and the amplitude of the transverse electric field within the cladding range, and J′ f (x) represents the first derivative of J f (x);
    根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度:
    Figure PCTCN2022122044-appb-100005
    Calculate the electric field strength within the fiber core according to the amplitude of the longitudinal electric field within the fiber core and the amplitude of the transverse electric field within the fiber core:
    Figure PCTCN2022122044-appb-100005
    根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度:
    Figure PCTCN2022122044-appb-100006
    Calculate the electric field strength within the cladding range based on the amplitude of the longitudinal electric field within the cladding range and the amplitude of the transverse electric field within the cladding range:
    Figure PCTCN2022122044-appb-100006
    其中,E m和E n分别表示在纤芯范围内和包层范围内的电场强度,电场圆周方向幅值
    Figure PCTCN2022122044-appb-100007
    e r表示r方向上的单位向量,
    Figure PCTCN2022122044-appb-100008
    表示
    Figure PCTCN2022122044-appb-100009
    方向上的单位向量,e z表示z方向上的单位向量。
    Among them, Em and En represent the electric field intensity within the core range and the cladding range respectively, and the amplitude of the electric field in the circumferential direction
    Figure PCTCN2022122044-appb-100007
    e r represents the unit vector in the r direction,
    Figure PCTCN2022122044-appb-100008
    express
    Figure PCTCN2022122044-appb-100009
    The unit vector in the direction, e z represents the unit vector in the z direction.
  4. 根据权利要求3所述的多芯少模光纤模间串扰检测方法,其特征在于,所述根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数前包括:The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 3, characterized in that the calculation is based on the electric field intensity within the fiber core range, the electric field intensity within the cladding range and the injection power. The inter-mode coupling mode coefficients include:
    将所述包层范围内的电场强度中,单芯横向极坐标系(r,θ)扩展为双芯的横向极坐标系(R,Θ),
    Figure PCTCN2022122044-appb-100010
    D表示芯间距。
    In the electric field intensity within the cladding range, the single-core transverse polar coordinate system (r, θ) is expanded to the dual-core transverse polar coordinate system (R, Θ),
    Figure PCTCN2022122044-appb-100010
    D represents the core distance.
  5. 根据权利要求4所述的多芯少模光纤模间串扰检测方法,其特征在于,所述根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及注入功率计算模式间耦合模系数包括:The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 4, wherein the inter-mode crosstalk is calculated based on the electric field intensity within the fiber core range, the electric field intensity within the cladding range and the injection power. Coupling mode coefficients include:
    计算纤芯范围内的电场强度、包层范围内的电场强度和纤芯折射率和包层折射率之差点积后的二重积分与角频率以及真空介电常数的乘积,再除以4倍的注入功率,得到模式间耦合模系数:Calculate the double integral of the electric field strength within the fiber core, the electric field strength within the cladding and the difference between the core refractive index and the cladding refractive index, the angular frequency and the vacuum dielectric constant, and then divide by 4 times Injection power, the inter-mode coupling mode coefficient is obtained:
    Figure PCTCN2022122044-appb-100011
    Figure PCTCN2022122044-appb-100011
    其中,ω是角频率,ε 0是真空介电常数,E m和H m分别代表在纤芯范围内的电场强度和磁场强度,而E n代表在包层范围内的电场强度,e z表示z方向上的单位向量,n 1和n 2分别表示纤芯的折射率和包层的折射率,符号*表示对矩阵进行共轭操作,P表示注入的功率, Among them, ω is the angular frequency, ε 0 is the vacuum dielectric constant, E m and H m represent the electric field intensity and magnetic field intensity in the core range respectively, and E n represents the electric field intensity in the cladding range, and e z represents The unit vector in the z direction, n 1 and n 2 represent the refractive index of the fiber core and the refractive index of the cladding respectively, the symbol * represents the conjugate operation of the matrix, P represents the injected power,
    Figure PCTCN2022122044-appb-100012
    Figure PCTCN2022122044-appb-100012
  6. 根据权利要求1所述的多芯少模光纤模间串扰检测方法,其特征在于,所述基于耦合模理论的芯间串扰计算方程为:The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 1, characterized in that the inter-core crosstalk calculation equation based on coupled mode theory is:
    Figure PCTCN2022122044-appb-100013
    Figure PCTCN2022122044-appb-100013
    其中,N为将光纤分段后的段数,d是将光纤分段后的段长,Δβ eq,mn,i是第i段的等效相位失配,
    Figure PCTCN2022122044-appb-100014
    为第i段修正后的纤芯间耦合模系数,k i(d)为纤芯间的耦合模系数。
    Among them, N is the number of segments after segmenting the optical fiber, d is the segment length after segmenting the optical fiber, Δβ eq,mn,i is the equivalent phase mismatch of the i-th segment,
    Figure PCTCN2022122044-appb-100014
    is the modified coupling mode coefficient between fiber cores in the i-th section, k i (d) is the coupling mode coefficient between fiber cores.
  7. 根据权利要求6所述的多芯少模光纤模间串扰检测方法,其特征在于,所述将纤芯间耦合模系数k i(d)替换为所述模式间耦合模系数k mn,i,计算得到光纤模间串扰值: The inter-mode crosstalk detection method of multi-core few-mode optical fiber according to claim 6, characterized in that the inter-core coupling mode coefficient k i (d) is replaced by the inter-mode coupling mode coefficient k mn,i , Calculate the crosstalk value between optical fiber modes:
    Figure PCTCN2022122044-appb-100015
    Figure PCTCN2022122044-appb-100015
    其中,g mn,i为第i段修正后的模式间耦合模系数。 Among them, g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
  8. 一种多芯少模光纤模间串扰检测装置,其特征在于,包括:A multi-core few-mode optical fiber inter-mode crosstalk detection device, which is characterized by including:
    光纤参数和注入功率获取模块,用于获取光纤参数和注入功率;Fiber parameter and injection power acquisition module, used to obtain fiber parameters and injection power;
    电场强度计算模块,用于根据所述光纤参数计算纤芯范围内的电场强度和包层范围内的电场强度;An electric field strength calculation module, used to calculate the electric field strength within the core range and the electric field strength within the cladding range based on the optical fiber parameters;
    模式间耦合模系数计算模块,用于根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数;An inter-mode coupling mode coefficient calculation module, used to calculate the inter-mode coupling mode coefficient according to the electric field intensity within the fiber core range, the electric field intensity within the cladding range, and the injection power;
    光纤模间串扰值计算模块,用于将基于耦合模理论的芯间串扰计算方程中的纤芯间耦合模系数替换为所述模式间耦合模系数,计算得到光纤模间串扰值。An optical fiber inter-mode crosstalk value calculation module is used to replace the inter-core coupling mode coefficient in the inter-core crosstalk calculation equation based on coupled mode theory with the inter-mode coupling mode coefficient, and calculate the optical fiber inter-mode crosstalk value.
  9. 根据权利要求8所述的一种多芯少模光纤模间串扰检测装置,其特征在于,应用于异质或同质光纤。A multi-core few-mode optical fiber inter-mode crosstalk detection device according to claim 8, characterized in that it is applied to heterogeneous or homogeneous optical fibers.
  10. 一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述一种多芯少模光纤模间串扰检测方法的步骤。A computer-readable storage medium, characterized in that a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the multi-core system according to any one of claims 1 to 7 is implemented. Steps of the inter-mode crosstalk detection method for few-mode optical fiber.
PCT/CN2022/122044 2022-09-16 2022-09-28 Multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus WO2024055360A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202211131009.7A CN115455355B (en) 2022-09-16 2022-09-16 Method and device for detecting crosstalk between modes of multi-core few-mode optical fibers
CN202211131009.7 2022-09-16

Publications (1)

Publication Number Publication Date
WO2024055360A1 true WO2024055360A1 (en) 2024-03-21

Family

ID=84304644

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/122044 WO2024055360A1 (en) 2022-09-16 2022-09-28 Multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus

Country Status (2)

Country Link
CN (1) CN115455355B (en)
WO (1) WO2024055360A1 (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706535A (en) * 2011-02-25 2012-10-03 住友电气工业株式会社 Measuring method for crosstalk between cores in multi-core optical fiber
JP2014206517A (en) * 2013-04-16 2014-10-30 日本電信電話株式会社 Evaluation method for crosstalk characteristic of multicore optical fiber, and system thereof
CN110445534A (en) * 2019-08-13 2019-11-12 中天宽带技术有限公司 A kind of method, system and equipment that multi-core optical fiber crossfire value determines
CN112803996A (en) * 2020-12-30 2021-05-14 中天通信技术有限公司 Method for detecting high nonlinear optical fiber coupling crosstalk
CN112859329A (en) * 2021-01-25 2021-05-28 苏州大学 Weak-coupling multi-core optical fiber crosstalk calculation method based on segmentation idea
CN113711094A (en) * 2019-04-25 2021-11-26 日本电信电话株式会社 Multi-core optical fiber and design method

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9195000B2 (en) * 2009-12-02 2015-11-24 Ofs Fitel, Llc. Techniques for reducing crosstalk in multicore fibers
WO2011071750A1 (en) * 2009-12-02 2011-06-16 Ofs Fitel Llc. A Delaware Limited Liability Company Techniques for manipulating crosstalk in multicore fibers
US10302463B2 (en) * 2015-11-19 2019-05-28 Corning Incorporated Distributed fiber sensors and systems employing multicore optical fibers
JP6334582B2 (en) * 2016-01-29 2018-05-30 株式会社フジクラ Multi-core fiber
JP6654104B2 (en) * 2016-02-26 2020-02-26 株式会社フジクラ Method and apparatus for measuring crosstalk of multi-core fiber
CN110927863B (en) * 2019-12-10 2020-08-25 东北大学 Multi-core few-mode micro-structure optical fiber used in space division-mode division multiplexing field
CN112198587B (en) * 2020-09-28 2021-12-31 北京邮电大学 Multi-core few-mode optical fiber and method for determining parameters thereof
CN112733073A (en) * 2020-12-30 2021-04-30 中天通信技术有限公司 Multi-core optical fiber crosstalk detection method based on coupling power theory

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706535A (en) * 2011-02-25 2012-10-03 住友电气工业株式会社 Measuring method for crosstalk between cores in multi-core optical fiber
JP2014206517A (en) * 2013-04-16 2014-10-30 日本電信電話株式会社 Evaluation method for crosstalk characteristic of multicore optical fiber, and system thereof
CN113711094A (en) * 2019-04-25 2021-11-26 日本电信电话株式会社 Multi-core optical fiber and design method
CN110445534A (en) * 2019-08-13 2019-11-12 中天宽带技术有限公司 A kind of method, system and equipment that multi-core optical fiber crossfire value determines
CN112803996A (en) * 2020-12-30 2021-05-14 中天通信技术有限公司 Method for detecting high nonlinear optical fiber coupling crosstalk
CN112859329A (en) * 2021-01-25 2021-05-28 苏州大学 Weak-coupling multi-core optical fiber crosstalk calculation method based on segmentation idea

Also Published As

Publication number Publication date
CN115455355A (en) 2022-12-09
CN115455355B (en) 2023-07-25

Similar Documents

Publication Publication Date Title
WO2022156762A1 (en) Weakly-coupled multi-core fiber crosstalk calculation method based on segmentation idea
Gan et al. Investigation of channel model for weakly coupled multicore fiber
Savović et al. Explicit finite difference solution of the power flow equation in W-type optical fibers
CN112803996B (en) Method for detecting high nonlinear optical fiber coupling crosstalk
CN110926758B (en) Device for measuring hybrid optical fiber mode
García et al. Bending and twisting effects on multicore fiber differential group delay
CN112733073A (en) Multi-core optical fiber crosstalk detection method based on coupling power theory
WO2024055360A1 (en) Multi-core few-mode fiber inter-mode crosstalk measurement method and apparatus
Sabitu et al. High dispersion four-mode fiber for mode-division multiplexing systems
Ferreira et al. Semi-analytical model for linear modal coupling in few-mode fiber transmission
Choutagunta et al. Designing high-performance multimode fibers using refractive index optimization
Anisimov et al. Fast multi step-index mode solver for analysis and optimization of optical fiber performance
Bhandari Orbital angular momentum (OAM) mode mixing in a bent step index fiber in perturbation theory
Mohapatra et al. Intermodal dispersion free few-mode (quadruple mode) fiber: A theoretical modelling
WO2023123651A1 (en) Multi-core fiber optic crosstalk detection method and apparatus, and storage medium
Ma et al. Finite-difference beam-propagation method for anisotropic waveguides with torsional birefringence
da Silva Ge‐doped microstrutured fiber design by genetic algorithm for directional coupling
Nesic et al. Transformation-optics modeling of 3D-printed freeform waveguides
Kazakov et al. Parallelization of the Generalized Multimode Nonlinear Schrödinger Equation Solver: A Performance Analysis
Kogan et al. Design and characterization of bend-insensitive multimode fiber
Bickham et al. Theoretical and experimental studies of macrobend losses in multimode fibers
Bourdine et al. Simulation and research of few-mode optical fiber DMD degradation due to geometry deviation from optimized form
Uvarov et al. GPU-assisted simulations of SDM systems
Velamuri et al. Study of bending losses in optical fibers using COMSOL
Kumar et al. Analysis of Crosstalk Performance and its Reduction Amount in Single-mode Homogeneous Trench-assisted Multicore Fiber

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22958550

Country of ref document: EP

Kind code of ref document: A1