CC BY
30
МЕТОД ФОКУСИРОВКИ ДЛЯ ИЗМЕРЕНИЯ ГЕОМЕТРООПТИЧЕСКИХ ПАРАМЕТРОВ СВЕТОВОДОВ
Андросик Андрей Борисович
канд. техн. наук, доцент МГОУ, РФ, г. Москва
Воробьев Сергей Андреевич канд. техн. наук, профессор МГОУ, РФ, г. Москва Мировицкая Светлана Дмитриевна
канд. техн. наук, доцент МГОУ, РФ, г. Москва E-mail: scotchwood@yandex. ru
THE FOCUSING METHOD FOR MEASUREMENT OF GEOMETRICAL AND OPTICAL PARAMETERS OF WAVEGUIDES
Audrey Androsik
candidate of Technical Sciences, Associate Professor of MSOU, Russia Moscow
Sergei Vorobev
candidate of Technical Sciences, Professor of MSOU, Russia Moscow
Svetlana Mirovitskaya
candidate of Technical Sciences, Associate Professor of MSOU, Russia Moscow
АННОТАЦИЯ
Рассмотрены особенности одного из наиболее удобных неразрушающих методов измерения основных геометрооптических характеристик волоконных световодов. Получено выражение для нахождения распределения показателя преломления. Приведены оценки точности измерений методом фокусировки.
ABSTRACT
The peculiarities of one of the most convenient non-destructive methods of measuring geometrical and optical parameters of optical fibers are reviewed. Expression for the location of the distribution of the refractive index is obtained. The accuracy of measurements by the focusing method is estimated.
Ключевые слова: волоконный световод; метод фокусировки; профиль показателя преломления; узкий пучок.
Keywords: the fiber waveguide; the focusing method; the profile of the refractive index; narrow beam.
Наиболее важными характеристиками волоконных световодов (ВС) являются геометрооптические — диаметр сердцевины, внешняя форма, распределение показателя преломления. При этом ВС должен быть однородным по всей длине. Отсюда вытекает необходимость разработки оперативных и высокоточных методов измерения формы и распределения показателя преломления оптических волокон (ОВ).
Метод фокусировки [1, с. 52, 2, с. 75] удобен для точного контроля профиля показателя преломления заготовок, а также оптических волокон и основан на фокусирующих свойствах сердцевины волокна, представляющей собой цилиндрическую линзу, при поперечном его освещении.
Рисунок 1. Схема измерения параметров ОВ методом фокусировки
На сердцевину под прямым углом к оси волокна, помещенного в иммерсионную жидкость, направляется параллельный некогерентный световой пучок постоянной интенсивности. Луч, входящий в сердцевину волокна (или заготовки) на расстоянии l от оптической оси (рис. 1), преломляется и выходит из сердцевины под углом, отличающимся от угла входа. В плоскости изображения волокна зондирующий луч попадает в точку на расстоянии у=у(1)
от оптической оси. Количественная связь профиля показателя преломления Ап
2 2
и координаты у(() описывается для параксиальных лучей уравнением d у^х = (1/п1) дп/ду, откуда наклон луча dy/dx на выходе из сердцевины определяется интегрированием:
(И) — ±Г Ё14х. (1)
Лх)х=ю п1 •'—ю ду
\йх/х=ж п1'}-жду
Для рассматриваемых осесимметричных объектов распределение показателя преломления является функцией радиальной координаты г,
связанной с переменной интегрирования х соотношением х — Vг2 — і2, т. е.
dx — гйг / Vг2 — I2, а следовательно, (dn/dy) = (дп/дг) (дг/ду) = (дп/дг) (у/г) ~ (дп/дг) (1/г). Внутри сердцевины волокна наблюдается слабое отличие у от I. Угол наклона луча с учетом (1) выражается интегральным соотношением
(Лу\ 21 гж дп дг
\йх)х=ж п1 ^1 дг ^г2—2 ( )
Вертикальная координата у луча в плоскости изображения (рис. 1)
у (I) — 1 + 1 (йу/йх )х=ю. (3)
Из (2) и (3) интегральное уравнение для искомой функции п(г), т. е. радиального распределения показателя преломления, имеет вид:
Сдд;-ші = шіу(і) — Ц (4)
Если функция [у(1)-1] известна, то распределение показателя преломления является решением интегрального уравнения (4). С введением обозначений
С(г) — п(г) — щ . р(і) — (5)
интегральное уравнение принимает вид:
Р(і) — С™ д° йг ґ(1) Ь дг(г2 — 12)1/2
Здесь верхний предел интеграла да, поскольку G(r) = 0 для г>а. Функция G (r) имеет вид:
G(r) — fj° g(u) cos (ur)du, (7)
тогда
F(l) — - j С ид (u)J0 (ul) du, (8)
где: J0 — функция Бесселя нулевого порядка.
Уравнение (8) является преобразованием Ханкеля известного обратного интегрального преобразования:
g(u) — -lf“lF(l)J0(ul)dl. (9)
Подставляя (9) в (7), можно получить развернутое интегральное представление G (r) через F(l):
G(r) — ~^С lF(l) [f0°° cos (ur)J0(ul) du]dl.
При интегрировании по u окончательно получается
G(r) = —^f;°(;ii^F(l)dl ; (10)
n(r) — n1 — — f™ l. dl. (11)
V j 1 nLJr v J
Входящая в основное интегральное уравнение (11) функция у(1) определяется путем прямого интегрирования плотности потока излучения Р(1).
Действительно, из закона сохранения энергии в потоке излучения следует, что мощность Р(у^у, переносимая через элемент длины dy, не зависит от координаты х. Поэтому Р(^1=Р(у^у, где Р(1) — плотность потока на входе в ОВ, а Р(у) — в плоскости изображения. При нормировании плотности потока излучения таким образом, чтобы Р(1) — 1, получается дифференциальное выражение dl=P(y)dy, после интегрирования которого
определяется обратная функция. Из нее находится и требуемая у(1), т. е. прямая функция, входящая в уравнение (11). При измерениях плоскость изображения определяется предметной плоскостью, на которую сфокусирован микроскоп.
Уравнение (11) справедливо, если объектив микроскопа, наведенный на плоскость изображения, согласован с показателем преломления путем погружения в иммерсионную жидкость. При отсутствии иммерсионной жидкости, т. е. когда объектив граничит с воздухом, необходимо учесть отражение и преломление света на границе жидкость — воздух. При этом (11) преобразуется к виду
Плоскость изображения необходимо удалить от оси сердцевины на расстояние до точки пересечения рефрагировавших лучей, покинувших ее. Наличие пересечения лучей нарушает однозначное соответствие между у и I, выражаемое соотношением (12). Оптимальным является расположение плоскости изображения непосредственно за сердцевиной, где Ь = а. При этом точка пересечения лучей, получаемая при продолжении входного и выходного направлений луча в сердцевину, расположена вблизи х = 0 (рис. 1).
1(У) — 1у Р(У'^у'
(12)
п(г) — п1 —
(13)
Интеграл (13), по которому рассчитывается распределение показателя преломления п (г), несколько напоминает аналогичный интеграл, получаемый с помощью метода интерферометрии при поперечном освещении. Однако между ними имеются два существенных различия. Во-первых, в случае фокусировки первичную информацию (плотность потока излучения Р(у)) в соответствии с (12) интегрируют, а не дифференцируют, что является значительным преимуществом, поскольку интегрирование приводит к сглаживанию случайных флуктуаций, вызванных шумами приборов, в то время как дифференцирование вызывает усиление сигнала с шумами. Во-вторых, хотя подынтегральное выражений в (11) имеет полюс при 1=г, его вклад в ошибку не стремится к бесконечности в точке г = 0, поскольку I - у (I) также обращается в нуль в этой точке. Погрешность метода фокусировки максимальна на оси волокна, но остается конечной.
Вклад в результирующую погрешность вносит также допущение, что у = I.
Отклонение луча внутри сердцевины определяется с учетом параболического распределения показателя преломления п (г) — п0[1 — (г/а)2А],где п0 — показатель преломления на оси; А — параметр профиля. Решение уравнения параксиального луча имеет вид у ~ 1[1 — (г/а)2А].
Координата г параллельна координате х (рис.1), но её начало сдвинуто таким образом, что в точке, где луч входит в сердцевину, г= 0. В точке выхода г — г (а2 — /2)1/2максимальное отклонение
Отсюда максимальное смещение луча в сердцевине 1 — 4^1 — ^) А.
^ Ут.
Среднее смещение луча по траектории внутри сердцевины --------------------------
У стандартных волокон А = 0,01 и смещение луча мало. Распределение плотности потока в стандартном волокне показано на рис. 2 сплошной кривой. На рис. 3.а сплошной линией изображен профиль показателя преломления ОВ, а штриховой — профиль того же волокна, определяемый методом интерферометрии среза, откуда видно, что центральный провал плохо разрешен методом фокусировки, если плоскость изображения выбрана на расстоянии
Рискнок 2. Распределение плотности потока в стандартном волокне
Ь=а. Снижение разрешающей способности вызвано тем, что центральный провал представляет собой сильную линзу диаметром несколько микрометров и лучи, проходящие через нее, пересекаются еще до того, как достигнут плоскости изображения. На рис. 3.б сплошной линий показан профиль, также определяемый методом фокусировки, плоскость изображения которого находилась в этом случае на расстоянии Ь= 16 мкм от оси волокна с радиусом сердцевины а =28 мкм. Наблюдаемый центральный провал углубляется при смещении плоскости изображения внутрь сердцевины. Как видно из рис. 3.б, размещение плоскости изображения в непосредственной близости от локализованных областей позволяет повысить, разрешение резко локализованных областей волокна со значительными градиентами показателя преломления (штриховой линией для сравнения показан профиль, полученный интерференционным методом). Однако повышение разрешающей способности в малых областях дается ценой потери информации о других областях.
Рисунок 3. Профили показателя преломления, измеренные методами фокусировки и интерферометрии среза
При вытягивании волокна из заготовки форма профиля показателя преломления сохраняется, трансформируется только радиальный размер, что обеспечивает достоверное представление о распределении показателя преломления волокна. Необходимо учитывать, что на измеренном с достаточно высокой точностью профиле показателя преломления заготовок видна микроструктура отдельных слоев сердцевины, в то время как в волокне эта структура размыта, поскольку ее характерные размеры соизмеримы с длиной волны зондирующего излучения.
Определение абсолютной точности метода измерения всегда является сложной задачей из-за наличия многих скрытых источников погрешностей. Однако уравнения (12) и (13) позволяют установить погрешность случайных и систематических флуктуаций уровня плотности потока излучения Р(у). Определив
п(г) —п1— Ап, (14)
можно записать, что Ап0 — решение уравнения (14) при отсутствии погрешности. Кроме того, у(1) — Уо(1) — Ау, где у0 — значение при отсутствии погрешности. Из (13) для усредненного значения (определяемого скобками ( )) можно получить
((Ап — Ап0)2) — (£) £ (I £ (V
(у(1)Ау(1'))
(12_г2)1/2(1'2_г2)1/2‘
С использованием соотношения dyldl= 1/Р, которое следует из (12), с помощью дифференцирования можно записать
Ау — А1/Р0, (15)
где: Р0 — идеальное значение. Из (15) можно получить
(АуАу(1')) = Шт
Ро(1)Ро(1'У
Функция корреляции для быстрых случайных флуктуаций плотности
([Р(у) — Ро(у)][Р(у') — Ро(у')]) — ([Р(у) — Ро(у)]2)°&(у — у'), (16)
где: о — длина корреляции (малая величина). Используя (16), можно переписать (12) в виде
(мм') — \°([р(у) — Ро(у)]2)у>у > у';
( ) 1с([Р(у) — Р0(у)]2)у',у<у':
Из рис. 2 видно, что Р(у) колеблется вокруг значений единицы, а у(1) — значения I. Поэтому приближенно можно считать Р0(у) =1 и у = I. Тогда порядок погрешности показателя преломления
[((Ап — Ап0)2)]1/2 — ^1{2о(а — г)([Р(у) — Р0(у)]2)}1/2.
Уравнение определяет погрешность профиля, вызванную быстрыми случайными пространственными флуктуациями плотности потока излучения.
Однако помимо случайных флуктуаций с очень малой длиной корреляций необходимо рассмотреть изменение Р(у) во времени. Поскольку облучение по полю наблюдения микроскопом не обязательно постоянно, необходимо замерить распределение интенсивности в отсутствие волокна и использовать эти уравнения для калибровки измеренной интенсивности излучения в присутствии сердцевины волокна. Любые флуктуации источника излучения в течение времени между калибровкой и измерениями будут искажать результаты определения профиля показателя преломления.
Длительные флуктуации излучения рассматриваются как систематическая погрешность 5Р плотности. Поскольку у(1) определяется по Р(у) из (12), он также вносит систематическую погрешность 5у. Погрешность показателя преломления может быть получена из выражения (13):
Используя (15), можно переписать (12) в виде
51 = !У5Р(у')с(у = у5Р,
и получить 5у = у[(5Р)/Р0]. Грубая оценка получается при использовании приближения Р0=1 и у=1. Интегрирование приводит к 5(Ап)=[(п1)/(л^)](а2—г2)12 5Р. Флуктуации повторных замеров показателя преломления при отсутствии усреднения 5(Ап)=0,003, а в случае усреднения по 100 замерам интенсивности излучения 5(Ап)=0,0005. Эту погрешность можно полностью отнести за счет временных флуктуаций излучения при 5Р=0,7 % без усреднения и 5Р=0,1 % с усреднением.
Описанный метод фокусировки справедлив только для объектов, помещенных в иммерсионную жидкость. Вместо иммерсионной жидкости так же можно ввести цилиндрические согласующие линзы. Эта модификация
метода фокусировки используется для контроля заготовок одномодовых волокон, но не применима для многомодовых, что обусловлено двумя факторами:
1. световой пучок внутри заготовки не может быть сколлимирован, что ограничивает область ввода излучения;
2. аберрации согласующих линз уменьшают внутри заготовки область, в которой коллимация пучка может быть достигнута оптическими методами.
На практике только половина центральной области поперечного сечения заготовки применима для метода фокусировки с использованием согласующих линз вместо иммерсионной жидкости. Преимуществом подобной модификации является устранение возможности разрушения поверхностей заготовки в процессе очистки при удалении остатков иммерсионной жидкости, что приводит к снижению прочности волокна. Кроме того, эта методика позволяет получить информацию о специальных заготовках с высокими показателями преломления или большими длинами волн облучения, для которых невозможно подобрать подходящую иммерсионную жидкость.
Список литературы:
1. Андросик А.Б., Воробьев С.А., Мировицкая С.Д. Рефракционный метод исследования волоконных световодов. Lambert Acamic Publishing 2012 — 183 с.
2. Лазарев Л.П., Мировицкая С.Д. Исследование оптических схем реализации измерителей геометрических характеристик оптических волокон и капилляров. ЦНИИТЭИ приборостр., Измерения, контроль, автоматизация — 1986, — № 3, — с. 52—86.
