CC BY
173
'Въотпьк, ЯАТ£\Оу£\А
Уфа: УГАТУ. 2012___________________________________________^¿____________________Т. 16, №8(53). С. 29-35
ЭЛЕКТРОНИКА, ИЗМЕРИТЕЛЬНАЯ ТЕХНИКА, РАДИОТЕХНИКА И СВЯЗЬ
УДК 681.7:535.5
А. Х. Султанов, В. Х. Багманов, Р. В. Кутлуяров, С.В. Харитонов
МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ 1Ш-СИГНАЛОВ В ВОЛОКОННО-ОПТИЧЕСКОЙ ЛИНИИ ПЕРЕДАЧИ В УСЛОВИЯХ НЕЛИНЕЙНОСТЕЙ И ПОЛЯРИЗАЦИОННО-МОДОВОЙ ДИСПЕРСИИ
Рассматривается математическая модель волоконно-оптической линии связи на основе системы связанных нелинейных уравнений Шредингера. С использованием модели проведены расчёты параметров распространения сигнала в типовой волоконнооптической линии передачи. Волоконно-оптические линии передачи; Керр-нелинейности; ПМД; связанные уравнения Шредингера
ВВЕДЕНИЕ
Уплотнение каналов по длине волны (WDM) на сегодняшний день является лидирующей технологией в магистральных волоконно-оптических линиях связи благодаря высокой емкости и спектральной эффективности. Однако увеличение скорости передачи до 40 Гбит/с и выше в отдельном канале связано с целым комплексом сложностей. Ключевыми факторами, существенно ограничивающими возможности волоконно-оптических систем связи, являются хроматическая дисперсия, поляризационные эффекты и нелинейности в волокне [1].
Поляризационные эффекты представляют достаточную сложность для изучения, так как существенно зависят от случайного пространственного и временного изменения параметра двулучепреломления оптического волокна. Случайное изменение параметров среды передачи происходит во временных масштабах от миллисекунд до дней [2].
Особый интерес представляют поляризационные эффекты, возникающие в волокне в условиях нелинейной передачи. Проявление таких эффектов имеет также случайный характер, и при скоростях передачи выше 10 Гбит/с и на длинных (сотни километров) линиях связи должно обязательно приниматься во внимание. В ряде работ показано, что при взаимодействии поляризационно-модовой дисперсии (ПМД) и нелинейностей существенно снижается эффективность оптических и электрических схем компенсации ПМД первого и более высоких порядков [3, 4].
Таким образом, при проектировании высокоскоростных волоконно-оптических линий пе-
редачи и модернизации существующих линии необходим оценочный анализ параметров линии с учетом хроматической дисперсии, нелинейных и поляризационных эффектов. Эффективным способом проведения таких оценок является численное моделирование линий передачи.
ТЕОРЕТИЧЕСКИЙ АППАРАТ МОДЕЛИРОВАНИЯ ЛИНИИ ПЕРЕДАЧИ
Феномен нелинейности возникает, когда амплитуда приложенного электрического поля достигает значений, при которых отклик среды не может рассматриваться как линейный. Нелинейности описываются уравнениями Максвелла с включением нелинейных членов. Уравнения Максвелла могут быть использованы для получения волнового уравнения, которое описывает распространение света в оптическом волокне [1]:
„ „ 1 Э2Е Э2Р
VxVx Е = —г—г- ц0-
(1)
с2 dt2 0 dt2 ’
где E - вектор электрического поля, P -электрическая поляризация, с - скорость света в вакууме, и использовано соотношение Цо£о = = 1/с2.
В случае оптических волокон в диапазоне длин волн 1-2 мкм наведенная поляризация состоит из двух слагаемых:
Р(Г t) = PL (Г t) + PNL (Г tX (2)
где линейный и нелинейный член определяются, соответственно, как
Pl (r,0 = £о j-¥X(1)(i -1 ')• E(r, t ')dt ' (3)
Контактная информация: 8(347)273-06-89
(4)
(/' = 1, 2, ...) - восприимчивость 7-ГО по-
-?2, ? - ?з); Е(г, ?! )Е(г, ?2 )Е(г, ?з), где У
рядка, х® - тензор ранга / + 1, вводимый для учета поляризационных эффектов.
С использованием (2), уравнение (1) может быть записано как
• = Ио
+ ^0
Э2РЛ,
(5)
с- Э?2 г и Э?2 г и Э?2
Математическое описание нелинейных поляризационных эффектов основано на том факте, что одномодовое волокно поддерживает две ортогонально-поляризованные моды с одинаковым пространственным распределением. Реальные волокна имеют некоторую асимметрию из-за неточностей технологического процесса и/или механического воздействия на волокно после производства. Асимметрия нарушает вы-рожденность ортогонально поляризованных НЕП мод, приводя к двулучепреломлению -разности в фазовых и групповых скоростях двух мод. При этом величина этой разности и ориентация осей двулучепреломления случайным образом изменяются на расстояниях порядка десятков метров.
Для получения уравнений, описывающих эволюцию двух ортогонально-поляризованных волоконных мод, можно описатьэлектрическое поле следующим образом:
Е(г,о=^[ш, Ме>ф(-/чО
+
+уЕу (г,?)ехр(-7Ю/) + с.с.\
4 (м)=^ (х, >0 4 (^ 0 ехР (Фо^) > ^ (г, о=^ (^, ^) 0 ехр (г'Ро^) ’
(6)
(7)
(8)
где х,у - два ортогональных единичных вектора, поперечных направлению распространения, 1'(х. у) - пространственное распределение моды, поддерживаемое волокном, А(г. /) - медленно меняющаяся амплитуда, (30 - постоянная распространения, с. с. - комплексно сопряженная величина.
Медленно меняющиеся амплитуды удовлетворяют системе связанных уравнений [1]:
аа.
2 д(2
+
■у 4 =*'у(|4
+
21,12 +з! 'I
4 +у^42ехр(_2/АРг)’
(9)
а л ■ (\ л I2
т4=ог 4
|4Г 14
/у
(10)
3 44 с\р(2/Др2
где А(3 = (30х - Рп, - Ро, Рь Рг - коэффициенты разложения постоянной распространения моды Р в ряд Тейлора вблизи несущей частоты со0-При этом
&=-, (П)
V
g
где групповая скорость.
Неравенство величин рЬ; Ф Р |Г приводит к уширению импульса по мере распространения в волокне. Величина уширения импульса может быть оценена с использованием временной задержки АТ, возникающей между двумя поляризационными компонентами при распространении импульса. Для волокна длиной Ь задержка определяется как
Ь Ь
А Т =
(I2)
Таким образом, второй член в левой части уравнений (9)—(10) обуславливает эффекты, связанные с двулучепреломлением.
Уравнение (12) определяет величину временной задержки между двумя ортогонально-поляризованными составляющими импульса в случае, если оси двулучепреломления не изменяют своего положения в волокне, т. е., оно справедливо для относительно коротких участков волокна (здесь мы имеем дело с ПМД первого порядка). При рассмотрении длинных волокон необходимо учитывать, что локальные оси двулучепреломления меняют свою ориентацию случайным образом от одного участка волокна к другому, и групповая задержка будет увеличиваться не линейно с расстоянием, а пропорционально его квадратному корню. Среднеквадратичное значение групповой задержки
Ат
(13)
где Ьс - длина корреляции, Ахъ - групповая задержка, получаемая на расстоянии, равном длине биений Ьъ, Ь - длина волокна.
Параметр р2 определяет хроматическую дисперсию волокна, которая принимается поляризационно-независимой.
Последний член в уравнениях (9) и (10) появляется вследствие когерентного связывания между двумя поляризационными компонентами
и ведет к вырожденному четырехволновому смешению. Если длина волокна существенно больше длины биений, последний член в уравнениях (9) и (10) часто меняет знак и его усреднённый вклад равен нулю.
Длина биений ЬБ = Х/Ди, где Ди - разность эффективных показателей преломления двух ортогональных мод в волокне представляет собой расстояние, на котором накапливается разность фаз в 2п между двумя модами или, что эквивалентно, поляризация вращается по полному циклу. Экспериментальные измерения на установленных волокнах дали усредненные значения ЬБ порядка 40 м для стандартных ступенчатых волокон и порядка 20 м для волокон со смещённой дисперсией и ненулевой смещённой дисперсией [5].
При распространении на большие расстояния экспоненциальным членом можно пренебречь, и уравнения (9) и (10) принимают вид
дА
— + PiX —- +
dz dt
дА- , ¡p2 д2А-
2 dt2
а .
+—А =
=¡ч (| 41' +1| А,|1 ] А.
¡Р2 д2 А
(14)
^+Р, —+
dz dt 2
■ (\А Р 2,. |2
=¡ч II ау| + з\А-\
dt2
АУ •
у a .
У +—А. =
(15)
Эти уравнения известны как связанные нелинейные уравнения Шредингера для линейно-двулучепреломляющего волокна.
Указанные уравнения управляют эволюцией поля одного оптического канала и могут быть обобщены для WDM-сигнала для описания межканальных эффектов [6].
Для 1-го канала ^-канальной WDM-системы электрическое поле
Е (г, ^ =
F, (XУ)А- (Zt)exp[/' (Poi-z -®Mt)] + cc
Fl (XУ) А,У (Zt)eXP[(P 0lyZ W0,t)] + cc-0
v J
(F ( X, У ) А- (z,t) exp [ (Poi-Z - w0,t)] + c-c^
F( X, У) А,у (Zt) eXP [ (P0lyZ - W0lt)] + c-c.
0
(16)
где Л1р(2, ?) и рМр - медленно меняющаяся огибающая и константа распространения р-й (р =
= х, у) компоненты канала I, соответственно, ^(х, у) - поперечное распределение моды канала I, которое приблизительно эквивалентно ^(х, у) в случае одномодового волокна.
Нелинейная компонента поляризации может быть записана как
( N
РNL (Г t) = ^
XР- eXP [¡' (b0ixz - W0lt)] + c-c-
Л
XР,У exP [ (Р
0,yZ -®0,?)] + c.c.
. (17)
Наконец, система уравнений, описывающих ПМД и различные нелинейные эффекты, может быть получена в следующем виде:
дА,р I р дА,р I ¡Р2р д А,р + а а =
+ _ А1„ =
dz
+ Р
1,Р
= ¡4,
АР +
(18)
dt 2 dt2 2
N 2 N
W +2 х ki + 3 Xk
j=1,j*, 3 j=1 j
+ig, X А1Атр eXP\ip2jp (W -wj )
j,m, j^m,
wl=2wj -“m 2
+igl 3 X ^jAq eXP [¡' (p2jp К
j,m,j^m,
-Wj )2 + ^ (“m -Wj ) I Z
где p, q = x, y и p Ф q , у - нелинейный коэффициент для канала l.
Первое слагаемое в правой части уравнения (18) описывает фазовые самомодуляцию и кросс-модуляцию, тогда как следующие части описывают четырехволновое смешение.
Фазовая кросс-модуляция (ФКМ) - это нелинейный феномен, возникающий при одновременной передаче двух и более оптических каналов по одному волокну. Проявление эффекта ФКМ зависит от суммарной мощности сигналов в волокне и разности поляризаций двух сигналов и приводит к флуктуации интенсивности и деполяризации сигнала на временных промежутках порядка длительности одного бита. Такое быстрое изменение поляризации сигнала на выходе линии обуславливает существенное снижение эффективности компенсаторов ПМД.
Нелинейная модуляция фазы, наводимая посредством ФКМ, зависит от последовательности бит в индуцирующем канале. Известно, что взаимное влияние WDM-каналов при ФКМ зависит от длины сдвига (walk-off) между битовыми по-
0
2
2
I
2
l-^j шга
следовательностями соседних каналов, возникающей из-за разности групповых скоростей двух WDM-каналов. Однако в протяженных WDM-системах с управлением дисперсией оптическая линия состоит из чередующихся отрезков с нормальной и аномальной дисперсией, так что дисперсия и потери в волокне компенсируются после каждого периода карты дисперсии. Поэтому, хотя импульсы в соседних каналах сдвигаются относительно друг друга на протяжении пролета, они снова «возвращаются» при компенсации дисперсии и снова взаимодействуют с битами соседнего канала в следующем пролёте линии. В результате ФКМ увеличивается от одного пролета к другому и может достигать достаточно больших значений для длинных линий передачи.
РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ
Уравнения (14)-(15), (18) в общем случае не могут быть решены аналитически.
Общепринятым подходом к решению связанных уравнений Шредингера является их численное решение.
Учитывая случайный характер моделируемых поляризационных эффектов, для получения наглядной картины исследуемого процесса передачи необходимо проведение ряда вычислительных экспериментов с получением ансамбля реализаций случайного процесса при каждом сочетании исходных условий.
Надлежащее моделирование ПМД требует воспроизводства ее случайной природы. Простейшая модель ПМД отображает эффект первого порядка и представляет собой отрезок дву-лучепреломляющего волокна, дающего временную задержку между двух ортогонально-поляризованных компонент поля. Более реалистичные модели должны воспроизводить ПМД-эффекты более высокого порядка. Это может быть достигнуто использованием хорошо-известной модели, состоящей из последовательности двулучепреломляющих волокон, повернутых на случайные углы относительно оси распространения сигнала. Для учета нелинейных эффектов каждое волокно должно быть смоделировано как нелинейная среда распространения с использованием системы уравнений Шредингера. В разработанной вычислительной модели связанные уравнения Шредингера решаются с использованием фурье-метода расщепления по физическим факторам [1].
Описанная модель схематически изображена на рис. 1. ПМД моделируемой линии определяется распределением Др^).
В серии проведенных симуляций была смоделирована волоконно-оптическая линия передачи, состоящая из четырех пролетов, каждый из которых включает в себя участок волокна с ненулевой смещенной дисперсией длиной 80 км, волокно, компенсирующее дисперсию, длиной 3 км и идеальный усилитель на основе волокна, легированного эрбием (на рис. 2 обозначены как NZDSF, DCF, EDFA соответственно; там же WDM-мультиплексор обозначен как MUX WDM, WDM-демультиплексор - как DEMUX WDM).
При симулировании рассматривался WDM-сигнал, состоящий из двух каналов на длине волны порядка 1550 нм с межканальным интервалом 0,8 нм. В каждом канале передавалась последовательность бит, образованная супер-гауссовскими импульсами. Параметры волокна, использованные при расчете, соответствуют справочным параметрам для NZDSF-волокон.
NZDSF DCF
о
Рис. 2. Структурная схема моделируемой ВОЛП
При симулировании было рассмотрено несколько случаев.
На рис. 3, а изображен ансамбль реализаций осциллограммы одноканального сигнала с пиковой мощностью 1 мВт после прохождения линии передачи, в которой отсутствует ПМД, то есть волокно не имеет свойств двулучепрелом-ления. На рис. 3, б изображен тот же сигнал после прохождения линии с ПМД (моделировалось волокно с низким значением параметра ПМД - Бр = 0,05 пс/^км). По результатам моделирования видно, что ПМД в случае низкой мощности сигнала в волокне ухудшает параметры сигнала на приеме из-за случайного фазового дрожания импульсов.
На рис. 4, а изображен ансамбль реализаций осциллограммы одного из двух каналов WDM-сигнала, пиковая мощность которого составляет 1 мВт после прохождения линии передачи, в которой отсутствует ПМД. При этом мощ-
ность второго канала принята 100 мВт, то есть он играет роль сигнала накачки, создающего нелинейные эффекты в волокне и вызывающего существенные искажения сигнала малой мощности. Поляризации двух канальных сигналов совпадают, что является худшим случаем с точки зрения их взаимного влияния. На рис. 4, б изображен случай, когда при тех же условиях два сигнала имеют ортогональные поляризации. При этом видно, что влияние сигнала накачки ниже, чем в предыдущем случае. На рис. 4, в, изображена осциллограмма того же сигнала при прохождении волокна со слабой ПМД при одинаковой исходной поляризации канальных сигналов. Слабая ПМД способствует снижению нелинейных помех, наводимых сигналом накачки на сигнал малой мощности, что связано с частичной деполяризацией сигналов.
Б
о
х
-1 0 1 Нормализованное время
б
Рис. 3. Передача сигнала малой мощности: а - в линии без ПМД; б - со слабой ПМД. Штриховой линией показана форма исходного сигнала
Нормализованное время
б
Нормализованное время
в
Рис. 4. Передача двухканального сигнала с накачкой большой мощности: а - в линии без ПМД при одинаковой поляризации каналов; б - в линии без ПМД при ортогональной поляризации каналов; в - в линии со слабой ПМД. Штриховой линией показана форма исходного сигнала
ВЫВОДЫ
На основе системы связанных нелинейных уравнений Шредингера создана компьютерная модельволоконно-оптической линии связи, позволяющая рассчитывать параметры передачи WDM-сигналов. Результаты проведенных расчетов показывают, что при низкой мощности сигналов ПМД оказывает большее влияние на систему по сравнению с нелинейностями. Когда суммарная мощность WDM-сигнала велика, нелинейные эффекты становятся доминирующими. При этом наличие в линии слабой ПМД незначительно снижает искажение сигналов, вызванное нелинейными эффектами. Полученные
при моделировании результаты применимы к реальным WDM-сетям с большими скоростями передачи. Разработанная модель может применяться при анализе существующих и проектировании новых волоконно-оптических линий передачи.
СПИСОК ЛИТЕРАТУРЫ
1. Agrawal G. P. Nonlinear Fiber Optics. Elsevier Inc., USA, 2007.
2. Menyuk C. R., Marks B. S. Interaction of polarization mode dispersion and nonlinearity in optical fiber transmission systems // Journal of Lightwave Technology. 2006. Vol. 24. P. 2806-2826.
3. Bit-pattem-dependent polarization rotation in first-order PMD-compensated WDM systems / R. Kho-sravani [et al.] // Optics Communications. 2006. № 257. P. 191-196.
4. Effects of nonlinear crosstalk in optical PMD compensation / J. H. Lee [et al.] // IEEE Photonics Technology Letters. Vol.14, № 8.
5. Measurement of local beat length and differential group delay in installed single-mode fibers / A. Galtarossa [et al.] // Journal of Lightwave Technology. 2000. Vol. 18(10). P. 1389-1394.
6. Impact of polarization mode dispersion and nonlinear effect on 40 Gbit/s dense wavelength division multiplexing system / X. Sun [et al.] // Front. Electr. Electron. Eng. China. 2006. Vol. 3. P. 361-366.
ОБ АВТОРАХ
Султанов Альберт Ханович, проф., зав. каф. теле-коммуникац. систем. Дипл. инженер по многока-нальн. электросвязи (Новосибирск. электротехн. ин-т связи, 1973). Д-р техн. наук по управлению в техн. сист. (УГАТУ, 1996). Иссл. в обл. телекоммуникац. систем, оптоэлектр. аэрокосм. систем.
Багманов Валерий Хусаинович, проф. той же каф. Дипл. физик (МГУ, 1975). Д-р техн. наук по системн. анализу, управлению и обработке информации (УГАТУ, 2007). Иссл. в обл. математическ. моделирования и обработки сигналов.
Кутлуяров Руслан Владимирович, асп. той же каф. Дипл. инженер по многоканальн. телекоммуникац. системам (УГАТУ, 2008). Иссл. в обл. волоконнооптической связи.
Харитонов Святослав Валериевич, магистрант Ин-та прикл. физики (г. Йена, ФРГ). Дипл. инженер по многоканальн. телекоммуникац. системам (УГАТУ, 2010). Иссл. в обл. волоконно-оптической связи.
