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 PDFInfo
- 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
Links
- 239000000835 fiber Substances 0.000 title claims abstract description 68
- 238000000691 measurement method Methods 0.000 title abstract 3
- 230000005684 electric field Effects 0.000 claims abstract description 115
- 230000008878 coupling Effects 0.000 claims abstract description 76
- 238000010168 coupling process Methods 0.000 claims abstract description 76
- 238000005859 coupling reaction Methods 0.000 claims abstract description 76
- 239000013307 optical fiber Substances 0.000 claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000002347 injection Methods 0.000 claims abstract description 25
- 239000007924 injection Substances 0.000 claims abstract description 25
- 238000003860 storage Methods 0.000 claims abstract description 11
- 238000005253 cladding Methods 0.000 claims description 67
- 238000001514 detection method Methods 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000009826 distribution Methods 0.000 abstract description 10
- 230000010287 polarization Effects 0.000 abstract description 9
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000005540 biological transmission Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting 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
本发明涉及光纤技术领域,尤其是指一种多芯少模光纤模间串扰检测方法、装置及计算机存储介质。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.
随着互联网数据流量的持续指数型增长,基于单模光纤(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:
其中,E
z1和E
z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数,
n
1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J
f(x)表示f阶的第一类贝塞尔函数,K
f(x)表示f阶的第二类修正贝塞尔函数,
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. 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, 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:
其中,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);
根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度:
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:
根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度:
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:
其中,E
m和E
n分别表示在纤芯范围内和包层范围内的电场强度,电场圆 周方向幅值
e
r表示r方向上的单位向量,
表示
方向上的单位向量,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 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.
优选地,所述根据所述纤芯范围内的电场强度、所述包层范围内的电场强度以及所述注入功率计算模式间耦合模系数前包括: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,Θ),
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, Θ), 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:
其中,ω是角频率,ε
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,
优选地,所述基于耦合模理论的芯间串扰计算方程为:Preferably, the inter-core crosstalk calculation equation based on coupled mode theory is:
其中,N为将光纤分段后的段数,d是将光纤分段后的段长,Δβ
eq,mn,i是第i段的等效相位失配,
为第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, 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:
其中,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.
为了使本发明的内容更容易被清楚的理解,下面根据本发明的具体实施例并结合附图,对本发明作进一步详细的说明,其中: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.
本发明的核心是提供一种多芯少模光纤模间串扰检测方法、装置及计算机存储介质,考虑了随距离变化的模式耦合系数,更加贴合实际光纤铺设情况。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:
其中,E
z1和E
z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数,
n
1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J
f(x)表示f阶的第一类贝塞尔函数,K
f(x)表示f阶的第二类修正贝塞尔函数,
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. 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, 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:
其中,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:
根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度:
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:
根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度:
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:
其中,E
m和E
n分别表示在纤芯范围内和包层范围内的电场强度,电场圆周方向幅值
e
r表示r方向上的单位向量,
表示
方向上的单位向量,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 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.
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,Θ),
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 ,Θ), 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:
其中,ω是角频率,ε
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,
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:
其中,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:
其中Δβ
eq,mn,i是第i段的等效相位失配,可表示为Δβ
eq,mn,i=β
eq,m,i-β
eq,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,i =β eq,m,i -β eq,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:
其中T是解的系数矩阵,可以表示为:where T is the coefficient matrix of the solution, which can be expressed as:
式中
为第i段修正后的耦合模系数,并且假设纤芯间的耦合模系数k
i(d)=K
mn,i(d)=K
nm,i(d)。如此,每一段电场幅值都可以通过前一段的值得到,同时也可以得到每段串扰的变化,将其叠加就可以计算出总串扰值:
in the formula 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:
我们只需要将其中的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:
其中,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)
- 一种多芯少模光纤模间串扰检测方法,其特征在于,包括: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.
- 根据权利要求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.
- 根据权利要求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:利用亥姆霍兹方程计算纤芯范围内和包层范围内的纵向电场的幅值: Use the Helmholtz equation to calculate the amplitude of the longitudinal electric field within the fiber core and within the cladding:其中,E z1和E z2分别表示纤芯范围内的纵向电场的幅值和包层范围内的纵向电场的幅值,A=jUC/βa,j是虚数单位,C是一个系统常数, n 1是纤芯折射率,β是纤芯传播常数,a是纤芯半径,k=2π/λ是波数,λ表示光波长,J f(x)表示f阶的第一类贝塞尔函数,K f(x)表示f阶的第二类修正贝塞尔函数, 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. 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, 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:其中,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);根据纤芯范围内的纵向电场的幅值和纤芯范围内的横向电场的幅值,计算纤芯范围内的电场强度: 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:根据包层范围内的纵向电场的幅值和包层范围内的横向电场的幅值,计算包层范围内的电场强度: 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:其中,E m和E n分别表示在纤芯范围内和包层范围内的电场强度,电场圆周方向幅值 e r表示r方向上的单位向量, 表示 方向上的单位向量,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 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.
- 根据权利要求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:
- 根据权利要求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:其中,ω是角频率,ε 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,
- 根据权利要求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:其中,N为将光纤分段后的段数,d是将光纤分段后的段长,Δβ eq,mn,i是第i段的等效相位失配, 为第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, 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.
- 根据权利要求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:其中,g mn,i为第i段修正后的模式间耦合模系数。 Among them, g mn,i is the modified inter-mode coupling mode coefficient of the i-th section.
- 一种多芯少模光纤模间串扰检测装置,其特征在于,包括: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.
- 根据权利要求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.
- 一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求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.
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)
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)
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 |
-
2022
- 2022-09-16 CN CN202211131009.7A patent/CN115455355B/en active Active
- 2022-09-28 WO PCT/CN2022/122044 patent/WO2024055360A1/en unknown
Patent Citations (6)
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 |