Вычислительные технологии
Том 9, № 6, 2004
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НЕЛИНЕЙНЫХ РЕЖИМОВ ПЕРЕДАЧИ ИНФОРМАЦИИ В ВЫСОКОСКОРОСТНЫХ ЛИНИЯХ ОПТИЧЕСКОЙ СВЯЗИ*
Ю.И. Шокин, М.П. Федорук Институт вычислительных технологий СО РАН, Россия e-mail: [email protected], [email protected]
А. Д. ШАПИРО Новосибирский государственный университет, Россия e-mail: [email protected]
Е. Г. Шапиро, С. К. Турицын Институт автоматики и электрометрии СО РАН, Россия e-mail: [email protected]
We present a novel nonlinear transmission regime at 40 Gbit/s using non-periodic carrier pulses with 50% duty cycle and sliding (with propagation) chirp-free points. The quasi-stable transmission of such carrier pulses is supported by large normal span average dispersion and misbalanced optical amplification.
Введение
Концепция управления дисперсией (dispersion managament) является ключевой технологией при разработке магистральных, сверхбыстрых волоконно-оптических линий связи (со скоростью передачи 40 Гбит/с и более в одном частотном канале). В дисперсионно-управляемых системах используются волоконные световоды с противоположной по знаку хроматической дисперсией, что позволяет контролировать дисперсионное уширение импульса. Если средняя дисперсия линии на одном периоде < D >= 0, то тогда в линейном случае в отсутствие затухания и шума форма сигнала будет полностью восстанавливаться в конце линии [1]. Уменьшение суммарной хроматической дисперсии позволяет подавить флуктуации положения отдельных битов — так называемый эффект Гордона — Хаусса [2]. Поэтому часто при проектировании современных волоконно-оптических линий связи фиксируют нулевую среднюю дисперсию.
Однако использование линейных режимов передачи информации по магистральным линиям связи со встроенными оптическими усилителями не является оптимальным с точ-
*Работа выполнена при поддержке Президента РФ (грант № НШ-2314.2003.1). © Институт вычислительных технологий Сибирского отделения Российской академии наук, 2004.
ки зрения минимизации коэффициента ошибки из-за сильного воздействия на маломощный оптический сигнал шумов усиленной спонтанной эмиссии. Эти шумы приводят к искажению формы сигнала в единичных битах и появлению паразитных возмущений поля в нулевых битах (на исходно пустых местах). Увеличение мощности оптических сигналов вызывает их деградацию вследствие нелинейности. При скоростях передачи данных 40 Гбит/с в одном частотном канале основным нелинейным фактором, ограничивающим пропускную способность линии, является взаимодействие соседних битов.
Одним из возможных путей уменьшения паразитного влияния нелинейных факторов является использование квазилинейных режимов передачи информации на основе оптических импульсов с длительностью, много меньшей чем длительность битового интервала [3, 4]. Однако короткая длительность несущих импульсов ведет к спектральному ушире-нию, что не всегда оптимально для режимов с высокой спектральной эффективностью в системах со спектральным уплотнением каналов.
Здесь мы продемонстрируем возможность квазиустойчивых нелинейных режимов передачи данных с длительностью оптических сигналов, составляющей 50% длительности битового интервала. Найденные квазиустойчивые импульсы существуют при большой нормальной дисперсии и разбалансировке оптического усиления и представляют новый тип нелинейных носителей информации. Показано, что достаточно большая нормальная средняя дисперсия приводит к существенному увеличению пропускной способности волоконной линии связи в режиме передачи данных N х 40 Гбит/с и 100 ГГц разделением между соседними частотными каналами.
1. Постановка задачи
Распространение оптических импульсов по волоконной линии связи с распределенной дисперсией описывается обобщенным нелинейным уравнением Шредингера для комплексной огибающей A электромагнитного поля [1]:
.OA в2(z) d2A . . , .|2 . . „ , ч .
^ - 2 ~а2 + * (z) |A|2A = iG (z) A (1)
Здесь z — расстояние вдоль линии; t — время; |A|2 — мощность сигнала; в2 — параметр дисперсии групповой скорости; а — коэффициент керровской нелинейности. Величины y, а и в2 представлены как функции от z, чтобы учесть изменения этих параметров при переходе от одного типа световода к другому. Коэффициент нелинейности а определяется , „ 2ПП2 „ „ Л
формулой а = -—-—, где n2 — нелинейный показатель преломления; Ао — несущая длина
A0Aeff
волны; Aeff — эффективная площадь собственной моды световода. Функция G(z) описывает процессы затухания и усиления оптического сигнала. Конкретный вид G(z) задается типом оптического усиления. Для численного решения уравнения (1) использовался метод операторной экспоненты [1].
Оценкой "качества" коммуникационной системы является величина коэффициента ошибки (Bit-Error Rate), которая определяет количество ошибочных битов к общему числу переданных битов [5].
Предположим, что плотности вероятностей нулей и единиц pi (i = 0,1) распределены по нормальному закону
(x - уi)2
Pi(x) = z1 exP V 2nai
2а2
где — средние значения; аг — дисперсии (среднеквадратичные отклонения). Тогда можно ввести стандартную величину — ^-фактор,
Q
— + ^0 ;
которая связана с BER следующим образом:
BER
2er/c (
exp
Q! ■ 2/2
В линиях, где используются теоретико-информационные методы коррекции ошибок (Forward Error Correction (FEC)) [6], приемлемыми значениями Q-фактора считаются 4 и более.
В расчетах передаточная способность волоконно-оптических линий связи анализировалась в терминах Q-фактора (Q > 6 для коэффициента ошибки BER < 10-9), полученного посредством усреднения 25 псевдослучайных последовательностей по 27 — 1 бит в каждой [7].
2. Конфигурации линий и параметры волоконных световодов
Нами без ограничения общности были рассмотрены волоконно-оптические линии связи, состоящие из 17 участков вида А. SMF +EDFA +DCF+ EDFA,
где SMF — стандартное одномодовое волокно и DCF — волокно с компенсацией дисперсии.
Длина SMF составляла 85 км, и для полной компенсации дисперсии длина DCF выбиралась равной 14.45 км. Для компенсации затухания сигнала в линии использовались два волоконных эрбиевых усилителя (EDFA) с коэффициентом шума 4.5 дБ и одинаковыми коэффициентами усиления 13.4 дБ. Заметим, что потери после прохождения сигнала в SMF не полностью компенсировались и полная компенсация потерь производилась после DCF, приводя к разбалансировке усиления и поддержанию квазиустойчивых, нелинейных режимов распространения.
Параметры волоконных световодов
SMF Затухание при 1550 нм 0.2 дБ/км
Эффективная площадь u • 104 80 мкм2
Дисперсия 17 пс/нм/км
Дисперсионный наклон u • 1040.07 пс/нм2/км
Нелинейный показатель преломления u • 1042.740-20 м2/Вт
DCF Затухание при 1550 нм 0.65 дБ/км
Эффективная площадь u • 10419 мкм2
Дисперсия — 100 пс/нм/км
Дисперсионный наклон u • 104 —0.41 пс/нм2/км
Нелинейный показатель преломления u • 1042.740-20 м2/Вт
Б. SMF + BRA + DCF + EDFA.
В этой конфигурации вместо эрбиевого волоконного усилителя после участка SMF для компенсации потерь использовался волоконный рамановский усилитель с обратной накачкой (BRA). Для изучения влияния распределенного усиления на устойчивость нелинейного режима выбирались различные соотношения между коэффициентами усиления волоконного рамановского и эрбиевого усилителей.
Средняя дисперсия линии (SAD) изменялась за счет небольшого изменения длины волокна с компенсацией дисперсии DCF. Численное моделирование было выполнено для восьми частотных каналов в диапазоне длин волн от 1550.1 до 1555.8 нм с 100 ГГц разделением между соседними каналами и скоростью передачи данных 40 Гбит/с в одном частотном канале (см. таблицу). В качестве единичных битов в расчетах использовались гауссовы импульсы с длительностью 12.5 пс и различной пиковой мощностью.
3. Результаты численного моделирования
Рассмотрим результаты численного моделирования. Для первой конфигурации, даже при относительно низких входных пиковых мощностях около 1 мВт, режим распространения сигнала является нелинейным. На рис. 1 показана величина ^-фактора после 1700 км распространения для 8 частотных каналов (вертикальные линии) в зависимости от средней дисперсии для начальной пиковой мощности 2.6 мВт. Видно, что оптимальным с точки зрения величины ^-фактора является режим с достаточно большой нормальной величиной средней дисперсии < Б >= -0.7 пс/нм/км. Большая величина средней дисперсии свидетельствует, что оптимальный режим распространения нелинейный. Для понимания физического механизма этого нелинейного режима рассмотрим более детально эволюцию
Рис. 1. Конфигурация Л. ^-фактор в 8 каналах после распространения сигнала на расстояние 1700 км в зависимости от средней дисперсии при входной пиковой мощности 2.6 мВт.
Рис. 2. Ширина импульса Тр^нм в зависимости от дистанции распространения з для трех различных величин нормальной средней дисперсии линии.
одиночного импульса в системе для различных значений средней дисперсии. Эволюция длительности оптического импульса с пиковой мощностью 2 мВт вдоль линии показана на рис. 2. Уширение импульса при нулевой средней дисперсии свидетельствует о влиянии нелинейности, и для стабилизации ширины импульса требуется конечная средняя дисперсия. Уширение несущих сигналов при нулевой средней дисперсии ведет к их деградации и уменьшению величины Q-фактора. Это уширение можно предотвратить уже для режимов с малой величиной нормальной средней дисперсии (-0.05 пс/нм/км), использование которых ведет к стабилизации ширины импульса. Отметим, что в этом случае стабилизация происходит при ширинах больше начальной ширины импульса TFWнм = 12.5 пс, как это видно из рис. 2.
Использование большой нормальной дисперсии (^ -0.7 пс/нм/км) приводит к стабилизации начальной ширины импульса и даже к его незначительному сжатию. Физический механизм стабилизации или даже слабой компрессии импульса подобен механизму, описанному в [8]. Распространение в нелинейном режиме приводит к появлению чирпа (chirp) у импульса. Далее, благодаря разбалансировке усиления распространение импульса по световоду с противоположным знаком хроматической дисперсии происходит в линейном режиме (с низкой мощностью), что приводит к сжатию импульса. Хотя режим распространения нелинейный, важно отметить, что данный эффект не является солитонным, поскольку сигнал, распространяющийся в линии, не имеет характеристик дисперсионно-управляемого солитона [9]. Наиболее важное отличие от дисперсионно-управляемого со-литона состоит в том, что точка нулевого чирпа в стандартном одномодовом волокне сдвигается на расстояние Az = L(D)/Dsmf после каждой секции, как это показано на рис. 3. Здесь L есть длина периодического участка линии. Для данной входной мощности после конечного числа секций точка нулевого чирпа исчезает и оптический импульс начинает уширяться от секции к секции. Данное обстоятельство приводит к усилению эффекта
Рис. 3. Положение точки Zo нулевого чирпа оптического импульса в зависимости от номера секции N.
Рис. 4. Зависимость ^-фактора от входной мощности сигнала: сплошная линия — < О >=0, штрихпунктирная — < О >= -0.65 пс/нм/км.
межсимвольного взаимодействия.
Мощность оптического сигнала имеет некое оптимальное значение, поскольку при малой мощности шумы усиленной спонтанной эмиссии эрбиевых усилителей увеличивают коэффициент ошибки в передаче данных. Использование более мощных импульсов увеличивает отношение сигнал/шум, но при этом усиливается паразитная роль нелинейных эффектов, которые ведут к деградации сигнала. Поэтому имеется такое значение пиковой мощности входных импульсов с фиксированной шириной, которое обеспечивает наилучший баланс между эффектами шума и нелинейностью с точки зрения величины ^-фактора. Поскольку режимы с относительно большой величиной нормальной средней дисперсии допускают квазиустойчивый нелинейный режим распространения без значительного "расплывания" оптических импульсов, в таких режимах уменьшается взаимодействие соседних битов. Данное обстоятельство позволяет при выбранной ширине импульса использовать более мощные входные сигналы, что в свою очередь улучшает отношение сигнал/шум без существенного ухудшения качества сигнала из-за нелинейности.
На рис. 4 приведена зависимость ^-фактора в худшем из каналов от пиковой мощности входного импульса при средних дисперсиях < Б >= 0 и < Б >= -0.65 пс/нм/км. Оптимальное значение пиковой мощности входного сигнала при < Б >= -0.65 пс/нм/км больше, чем при нулевой дисперсии, однако уменьшение отрицательной роли нелинейных эффектов в этом случае приводит к уменьшению коэффициента ошибки (увеличению ^-фактора).
Необходимо подчеркнуть, что найденный нелинейный режим со скользящей точкой нулевого чирпа не зависит от конкретной конфигурации линии и может быть реализован для различных систем. Подобные результаты получены для второй конфигурации оптической линии связи. Эволюция ширины импульса для случаев, когда волоконный ра-мановский усилитель компенсирует 1 и 10 дБ затухания сигнала для различных средних
Рис. 5. Зависимость длительности импульса Т^нм от дистанции распространения з для коэффициента рамановского усиления 1 дБ в гибридной схеме усиления Б.
Рис. 6. Зависимость длительности импульса Т^НМ от дистанции распространения з для коэффициента рамановского усиления 10 дБ в гибридной схеме усиления Б.
<D> -0.2
-0.3
-0.4
-0.5
-0.6
9
10
11
12 т 13
1 T7VU1M
FVHM
Рис. 7. Конфигурация A. Величина средней дисперсии линии в зависимости от начальной длительности входных сигналов при средней мощности Pave = -1.5 дБм.
дисперсий линии показана на рис. 5 и 6 соответственно. Как можно видеть, тип усиления фактически не влияет на установление квазиустойчивого нелинейного режима, хотя, как показывают расчеты, оптимальная пиковая мощность сигнала зависит от типа усиления. Минимальная величина нормальной средней дисперсии, которая приводит к стабилизации ширины импульса, показана на рис. 7. Можно видеть, что более широкие импульсы требуют большей величины нормальной средней дисперсии для стабилизации уширения.
Заключение
Методами математического моделирования исследованы нелинейные режимы распространения оптических импульсов в режиме передачи данных 40 Гбит/с с длительностью оптического импульса, составляющей 50% битового интервала, и скользящей (с дистанцией распространения) точкой нулевого чирпа. Квазиустойчивое распространение таких импульсов возможно в режиме большой нормальной средней дисперсии и разбалансировки оптического усиления. Достигнутая в нелинейных режимах передачи информации величина Q >4 позволяет применять для таких линий теоретико-информационные методы коррекции ошибок (Forward Error Correction (FEC)).
Список литературы
[1] Agrawal G.P. Nonlinear Fiber Optics. N.Y.: Acad. Press, 2001.
[2] Gordon J.P., Haus H.A. Random walk of coherently amplified solitons in optical fiber transmission // Opt. Lett. 1986. Vol. 11. P. 665-667.
[3] Park S.-G., Gnauck A.H., Wiesenfeld J.M., Garrett L.D. 40-Gb/s transmission over multiple 120-km spans of conventional single-mode fiber using highly dispersed pulses // IEEE Photon. Technol. Lett. 2000. Vol. 12. P. 1085-1087.
[4] Gnauck A.H., Park S.-G., Wiesenfeld J.M., Garrett L.D. Highly dispersed pulses for 2 x 40 Gbit/s transmission over 800 km of conventional singlemode fibre // Electron. Lett. 1999. Vol. 35. P. 2218-2219.
[5] Agrawal G.P. Fiber-Optic Communication Systems. Second edition. N.Y.: John Wiley & Sons, Inc., 1997.
[6] Grover W.D. Forward error correction in dispersion-limited lightwave systems //J. Lightwave Technol. 1988. Vol. 6. P. 643.
[7] Shapiro E.G., Fedoruk M.P., Turitsyn S.K. Numerical estimate of BER in optical systems with strong patterning effects // Electron. Lett. 2001. Vol. 37, N. 19.
[8] Fisher R.A., Bishel W.K. Pulse compression for more efficient operation of solid-state laser amplifier chains // Appl. Phys. Lett. 1974. Vol. 24. P. 468.
[9] Turitsyn S.K., Shapiro E.G., Medvedev S.B., Fedoruk M.P., Mezentsev V.K. Physics and mathematics of dispersion-managed optical solitons // Comptes Rendus. Physique. 2003. Vol. 4. P. 145-161.
Поступила в редакцию 27 октября 2004 г.