COMPUTING OF PARAMETERS FOR SCHEME OF MDM CHANNEL PRECISION POSITIONING AT THE FEW-MODE OPTICAL FIBER END

Full Text

На сегодняшний день сильное проявление нелинейных эффектов в волоконно-оптических линиях передачи (ВОЛП) является ключевым фактором, ограничивающим реализацию транспортных сетей связи нового поколения, ориентированных на передачe данных со скоростью сотни Тбит/с и выше [1]. Поэтому переход на маломодовые оптические волокна (ОВ) (FMF - Few Mode Fibers) с экстремально увеличенным, по сравнению с традиционными телекоммуникационными коммерческими одномодовыми [2] и маломодовыми [3] (ОВ), диаметром сердцевины, обеспечивающим практически полное устранение нелинейности самого волоконного световода, в сочетании с пространственным мультиплексированием (SDM - Space Division Multiplexing) на базе технологий модового уплотнения каналов (MDM - Mode Division Multiplexing) и MIMO (Multiple Input Multiple Output - множественные входы/множественные выходы), является одним из перспективных подходов для решения данной задачи [4 - 7]. Очевидно, что увеличение диаметра сердцевины FMF неизбежно приводит к появлению еще большего числа новых модовых составляющих высших порядков, удовлетворяющих условию отсечки рассматриваемого ОВ. В результате требуется проведение дополнительных мероприятий по уменьшению дифференциальной модовой задержкой (ДМЗ) - основного фактора линейных искажений оптического сигнала при распространении по ОВ в маломодовом режиме [8]. Неслучайно при разработке конструкций FMF ключевой проблемой является поиск компромисса между эффективной площадью сечения ОВ, непосредственно связанной с диаметром сердцевины, который, в том числе, определяет модовый состав волоконного световода, и минимизацией ДМЗ на заданной оптической несущей или в определенном выделенном диапазоне длин волн [6, 7, 9, 10]. В то ж время, из теории связи мод хорошо известно [11], что перераспределение мощности оптического сигнала, поступаемого с выхода когерентного источника в торец возбуждаемого волоконного световода, осуществляется между ограниченным числом мод этого ОВ, количество и порядок которых определяются геометрией и параметрами возбуждаемого световода, исходным модовым составом (поперечными модами) излучения лазера и непосредственно самими условиями ввода. Так, например, при центрированном вводе, а также вводе с прецизионным радиальным рассогласованием, мощность поступающих с выхода источника поперечных мод оптического сигнала передается модам возбуждаемого ОВ только аналогичного азимутального порядка [11, 12]. Безусловно, при дальнейшем распространении сигнала по ОВ ввиду его нерегулярности, а также наличия микро- и макроизгибов, неизбежно возникающих на практике при инсталляции реальных кабельных линий, эта мощность будет также частично передаваться новым модам, но вновь ближайших азимутальных и радиальных порядков [13 - 15]. Однако если изначально подобрать условия ввода сигнала с выхода лазера в торец волокна таким образом, чтобы практически вся мощность поступаемого излучения передавалась непосредственно одной определенной моде заданного порядка, то и амплитуды новых, возбуждаемых в процессе распространения по ОВ, компонентов будут пренебрежимо малы, по сравнению с ней [16, 17]. В данной работе предложена методика расчета параметров схемы прецизионного позиционирования каналов пространственного мультиплексора, реализующего технологию MDM, на торце возбуждаемого FMF ВОЛП. На основании представленной методики получены некоторые предварительные результаты теоретического исследования потенциальных возможностей выбора оптимальных параметров передачи вводимых мод и их позиционирования на торце кварцевого FMF с увеличенным диаметром сердцевины 42 мкм, с точки зрения уменьшения ДМЗ за счет передачи мощности вводимого излучения каналов MDM ограниченной группе направляемых мод возбуждаемого ОВ. Методика расчета коэффициентов связи мод на вводе в ОВ, выполненном с прецизионным радиальным рассогласованием На сегодняшний день расчет коэффициентов связи по-прежнему остается одним из наиболее известных и простых подходов к анализу перехода оптического излучения от одной моды к другой. Данный метод широко применяется при оценке вносимых потерь и коэффициента отражения на соединении ОВ одинаковой конфигурации с разбросом отдельных технологических параметров - например, [12, 18 - 20] и др., моделировании и исследовании влияния условий ввода оптического излучения на возбуждение оптических волноводов ([11, 21 - 23] и др.), в том числе и при моделировании распространения оптических сигналов по ОВ в маломодовом режиме ([8, 13, 14, 24, 25] и др.). В общем случае коэффициент связи моды m, возбуждающей моду с индексом n, определяется через интеграл перекрытия полей взаимодействующих мод p и q (LPlm(p) и LPlm(q)), который в скалярной постановке задачи для цилиндрической системы координат имеет следующий вид [8, 11 - 14, 18 - 25]: , (1) где и - радиальные распределения полей взаимодействующих - мод поступающей с выхода источника моды LPlm(p) и возбуждаемой в ОВ моды LPlm(q). Очевидно, что наиболее простая структура поля относится к основной моде LP01. Это позволяет перейти непосредственно к аналитической форме записи интеграла перекрытия полей (1), используя, например, для представления поля соответствующие аппроксимирующие выражения или разложение в ряд в базисе заданных функций. Подобный подход активно применяется, например, для анализа стыков одномодовых ОВ, что и отражено в перечисленных выше публикациях - в частности, [11, 18, 19, 21, 23] и др. Напротив, структура поля высших мод носит более сложный характер. По этой причине в работах, ориентированных на анализ взаимодействия мод высших порядков - например, [8, 24, 25] и др. - интеграл числителя выражения (1) предлагается брать численно или использовать простые аппроксимирующие выражения, соответствующие точным решениям волнового уравнения. В частности, в рамках приближения Гаусса [11] для волокон с идеальным неограниченным параболическим профилем показателя преломления, в состав которого входит удобная для интегрирования функция Лагерра-Гаусса: , (2) где l и m - азимутальный и радиальные порядки моды LPlm, соответственно; - нормированная радиальная координата; r - радиальная координата; a - радиус сердцевины; - нормированный радиус пятна моды; r0 - радиус пятна моды; - полином Лагерра. Однако и в этом случае возникают определенные трудности перехода к аналитическим выражением для коэффициента связи мод неединичного радиального порядка: при m>1 в подынтегральном выражении числителя формулы (1) фигурирует произведение полиномов Лагерра разного порядка, что, в общем случае, не относится к «табличным» интегралам. По этой причине, например, в монографии [11] представлены решения только для основной моды и высших мод 1го радиального порядка m=1 LPl1. Кроме того, представление профиля показателя преломления ОВ с помощью набора простых гладких степенных функций или, тем более, только одной степенной функции - в частности, неограниченной параболы - для перехода к описанию распределения поля моды не является корректным как с точки зрения расчета параметров модового состава реальных промышленных образцов ОВ, градиентный профиль показателя преломления которых существенно отличатся от гладких модельных a-профилей [26, 27], так и тем более специализированных профилей FMF, обеспечивающих уменьшенную ДМЗ. Для решения описанной проблемы предлагается воспользоваться ранее разработанной модификацией приближения Гаусса (ОМПГ) [15], обобщенной на случай расчета спектральных характеристик дисперсионных параметров направляемых мод произвольного порядка, распространяющихся в слабонаправляющих ОВ с произвольным осесимметричным профилем показателя преломления. Предложенный метод базируется на сочетании приближения Гаусса [11], модифицированного для расчета многомодовых ОВ [15], и методе стратификации [28]. Здесь, в отличие от известных решений , базирующихся на классическом приближении Гаусса [11] и его различных модификациях [29 - 35] предлагается рассматривать исследуемый слабонаправляющий волоконный световод с произвольным осесимметричным профилем показателя преломления, ограниченным одной внешней сплошной оболочкой, как слабонаправляющий волоконный световод с многослойным профилем показателя преломления, представляемый в области сердцевины ОВ в виде конечного числа N слоев, в пределах которых значение показателя преломления остается постоянным. В результате при задании исходных данных обеспечивается детализированное воспроизведение профиля показателя преломления исследуемого ОВ, что существенно снижает погрешность расчетов [36]. Таким образом, с точки зрения предложенного приближенного метода ОМПГ, подстановку выражения (2), описывающее распределение поля моды, в интеграл перекрытия полей (1) можно полагать корректным. Действительно, ОМПГ использует для представления поля моды заданного порядка простое аппроксимирующее выражение (2), однако при этом учитывается реальная структура профиля исследуемого градиентного ОВ, что и обеспечивает малую погрешность расчетов [36]. Как было отмечено выше, при подстановке обобщенного на произвольный порядок моды аппроксимирующего выражения (2) в интеграл перекрытия полей (1) для высших мод возникают определенные трудности интегрирования числителя (1). В первую очередь это связано с интегралами от произведения полиномов Лагерра разной степени и порядка, которые, в общем случае, не относятся к табличным. В рамках данной работе эта проблема преодолевается путем использования формулы явного представления полинома Лагерра в виде конечного степенного ряда соответствующего вида [37]. Все это позволило получить достаточно громоздкие, содержащие несколько конечных вложенных сумм, но при этом аналитические формулы для оценки коэффициентов связи мод произвольного порядка на стыке двух слабонаправляющих ОВ с неодинаковыми параметрами, выполненном с осевым смещением. В результате обобщенное аналитическое выражение для коэффициентов связи мод LPlm и LPln произвольного порядка на стыке ОВ с осевым смещением d будет иметь вид: , (3) где ; ; ; ; ; ; . . (4) - соответствующие коэффициенты степенного ряда формулы явного выражения полинома Лагерра . Для перехода от эквивалентных, полученных в результате расчета на основании ОМПГ, значений радиусов пятна моды к фактическим, предлагается воспользоваться хорошо известной интегральной формой записи формулы поправки Петерманна-2 для радиуса пятна моды в ближнем поле () [39], которое в результате соответствующих преобразований после подстановки аппроксимирующего выражения для радиального распределения поля моды (2) приводится к аналитической формуле (4). Таким образом, методика расчета коэффициентов связи поперечных мод оптического излучения, поступающих с выхода портов MDM-мультиплексора, на торец FMF ВОЛП включает в себя следующие этапы. Предварительно задается профиль показателя преломления и базовые геометрические параметры исследуемого ОВ. Далее с помощью ОМПГ осуществляется расчет эквивалентных, в рамках приближения Гаусса, значений радиусов пятна мод модового состава FMF. Затем выполняется переход от эквивалентных к фактическим значениям радиусов пятна мод с использованием выражения для поправки Петерманна-2 (4). После чего для оценки коэффициента связи в формулу (3) осуществляется подстановка полученных значений радиусов пятна моды соответствующего порядка FMF, этого же параметра соответствующей поперечной моды излучения канала MDM и непосредственно самого прецизионного радиального смещения ввода этой моды относительно центра сердцевины исследуемого ОВ. Результаты расчетов: выбор оптимальных параметров схемы прецизионного позиционирования каналов MDM на торце FMF 42/125 В качестве примера в данном разделе представлены результаты анализа распределения коэффициентов связи мод FMF соответствующих азимутальных порядков в заивисимости от комбинации радиуса пятна моды 1го радиального порядка LPl1 излучения канала MDM и прецизионного радиального смещения ввода d этой моды относительно центра сердцевины возбуждаемого ОВ. Рис. 1. Градиентный профиль показателя преломления FMF 42/125 с умееьшенной ДМЗ Рассматривается кварцевое слабонаправляющее FMF с диаметром сердцевины 42 мкм и оболочкой 125 мкм. Данное ОВ поддерживает распространение до 20 направляемых мод LPlm - до 6го азимутального l=0…6 и 4го радиального m=1…4 порядков включительно в «С»-диапазоне длин волн. Полученный в результате оптимизиции специализированный профиль показателя преломления, внешний вид которого приведен на рис. 1, обеспечивает значение ДМЗ относительно основной моды менее 30 пс/км на оптической несущей l=1550 нм. Предполагается, что при определенной комбинации радиуса пятна моды LPl1 соответствующего канала MDM и ее положения на торце, смещенного относительно центра сердцевины возбуждаемого ОВ, можно добиться передачи практически всей ее мощности одной определенной моде FMF этого же азимутального порядка. Для проверки данной гипотезы на основании представленной методики анализа ввода оптического излучения каналов MDM в торец ОВ был проведен расчет коэффициентов связи мод FMF и, соответственно, вводимых мод в зависимости от радиуса пятна последних и их положения относительно центра сердцевины возбуждаемого ОВ. На первом этапе предлагается рассмотреть 5-канальную систему MDM, оснащенную мультиплексором, обеспечивающим ввод 5 мод 1го радиального порядка LP0m…LP4m. Рассматривался диапазон значений радиуса пятна мод каналов MDM =0,5…8,0 мкм и осевого смещения d=0,0…18,0 мкм. Полученные результаты расчетов коэффициентов связи позволили перейти к построению диаграмм зависимости минимального разброса нормированных амплитуд возбуждаемых мод FMF в зависимости от перечисленных выше параметров. В качестве примера, на рис. 2 представлены указанные диаграммы для мод 0го (а), 3го (б) и 4го (в) азимутальных порядков. Анализ построенных диаграмм позволил выделить две локации торца FMF, соответствующих непосредственно центру сердцевины ОВ d=0 мкм и ее периферийной части, отстоящей от последнего на расстоянии d=4 мкм, при которых для всех 5ти перечисленных модовых групп LP0m…LP4m минимальное отклонение нормированных амплитуд внутри одной группы идентичного азимутального порядка l=0…4 достигает значения от 0,85 и более. (а) (б) (в) Рис. 2. Диаграмма зависимости минимального разброса нормированных амплитуд возбуждаемых мод FMF в зависимости от радиуса пятна вводимой моды LPl1 и ее радиального смещения d относительно центра ОВ: (а) моды LP0m; (б) моды LP3m ; (в) моды LP4m Далее для указанных локаций были построены кривые зависимости минимального разброса нормированных амплитуд возбуждаемых мод FMF от радиуса пятна вводимой моды LPl1 при фиксированном радиальном смещении d относительно центра ОВ, представленные на рис. 3. (а) (б) Рис. 3. Кривые зависимости минимального разброса нормированных амплитуд возбуждаемых мод FMF в зависимости от радиуса пятна вводимой моды LPl1 при фиксированном радиальном смещении d относительно центра ОВ: (а) центрированный ввод d=0,0 мкм; (б) d=4,0 мкм Это позволило выделить искомые значения радиуса пятна вводимых мод LPl1 каналов MDM, при которых достигается искомое максимальное отклонение нормированных амплитуд возбуждаемых мод одного азимутального порядка, составляющих более 0,85. А последующее сопоставление распределения расчетных значений коэффициентов связи мод внутри одной группы - определить, какой именно моде передается основная часть мощности сигнала канала MDM LPl1. Далее на основании анализа диаграммы ДМЗ, были отобраны оптимальные, с точки зрения, минимизации ДМЗ по основной моде LP01, комбинации возбуждаемой моды FMF заданного порядка, радиуса пятна вводимой моды канала MDM и ее позиционирования на торце исследуемого ОВ d. Анализ полученных результатов показал следующее: центрированный ввод целесообразно осуществлять для моды LP21, при этом ее радиус пятна моды должен составлять =2,45 мкм. В этом случае практически вся мощность данного канала MDM перейдет моде FMF этого же порядка, при этом минимальный разброс нормированных амплитуд внутри модовой группы ОВ 2го азимутального порядка l=2 составит 0,99. Для остальных четырех мод каналов MDM можно рекомендовать позиционирование в периферийной зоне сердцевины возбуждаемого ОВ на расстоянии d=4 мкм от центра последней. В этом случае при вводе основной моды LP01 c радиусом пятна =2,08 мкм в FMF будет также возбуждаться основная мода LP01 нормированная амплитуда которой будет отличаться от аналогичного параметра мод LP0m этой же группы 0го азимутального порядка l=0 более, чем на 0,99. В свою очередь, ввод моды LP11 канала MDM при =0,88 мкм и аналогичном смещении d=4 мкм от центра сердцевины обеспечит возбуждение моды LP12 с отклонением нормированной амплитуды от мод этой же группы l=1 более 0,99. В свою очередь, при вводе LP31 с радиусом пятна =0,80 мкм и d=4 мкм практически вся мощность это моды переходит компоненте LP32 с отклонением нормированной амплитуды в этой группе l=3 не менее, чем на 0,89. Наконец, для группы мод LP4m ввод LP41 с радиусом пятна =0,58 мкм при таком же смещении d=4 мкм обеспечит возбуждении моды FMF LP42, нормированная амплитуда которой будет превышать остальные компоненты из той группы l=4 боле, чем на 0,99. Для оценивания потенциальных возможностей управления ДМЗ модового состава FMF 42/125 за счет выбора соответствующих параметров схемы позиционирования вводимых мод каналов MDM на торце исследуемого ОВ, был выполнен расчет спектральных характеристик ДМЗ в «С»-диапазоне длин волн для всех мод рассматриваемого FMF, удовлетворяющих условию отсечки, и 5 возбуждаемых мод (LP01, LP12, LP21, LP32, LP42) при контроле указанных условий ввода. Данные кривые представлены на рис. 4. Анализ полученных результатов показывает, что выбор соответствующей комбинации параметров схемы позиционирования вводимых мод каналов MDM на торце исследуемого FMF 42/125 позволяет уменьшить локальное значение ДМЗ на искомой оптической несущей 1550 нм с 27,70 пс/км до 19,8 пс/км, увеличить на 18,5 нм относительного исходного значения спектральный диапазон центральной области «С»-диапазона, соответствующий пониженным значения ДМЗ по основной моде до 60 пс/км и, в целом, более в чем в 2 раза уменьшить ее значение за пределами указанной области длин волн. Рис. 4. Спектральные зависимости ДМЗ по основной моде FMF 42/125 в пределах волнового «С»-диапазона без и с подбором параметров схемы позиционирования мод каналов MDM на торце исследуемого ОВ Заключение Представлена методика расчета параметров схемы прецизионного позиционирования каналов MDM на торце возбуждаемого ОВ ВОЛП. Проведен расчет параметров ввода 5 мод низших азимутальных порядков системы MDM в торец FMF с диаметром сердцевины 42 мкм и уменьшенной дифференциальной модовой задержкой (ДМЗ). Показано, что при соответствующем сочетании значений радиуса пятна моды канала MDM и радиального рассогласования на вводе удается добиться передачи практически всей мощности вводимой моды MDM отдельной направляемой моде ОВ этого же азимутального порядка. На основании результатов проведенных теоретических расчетов получены практические рекомендации по выбору параметров схемы позиционирования каналов MDM на торце ОВ 42/125. Продемонстрирована потенциальная возможность снижения ДМЗ в центральной области «С»-диапазона при соответствующем подборе комбинации параметров схемы позиционирования: радиуса пятна вводимой моды канала MDM, ее положения на торце FMF относительно центра сердцевины и порядок возбуждаемой моды ОВ, нормированная амплитуда которой превышает амплитуды остальных мод из этой же группы аналогичного азимутального порядка не менее, чем на 0,85.

About the authors

Anton Vladimirovich Bourdine

Povolzhskiy State University of Telecommunications and Informatics

Email: bourdine@yandex.ru

References

Supplementary files