Выпуск #7-8/2020
А.Кривошеев, Ф.Барков, Ю.Константинов
КОНТРОЛЬ ВОЛС МЕТОДОМ БРИЛЛЮЭНОВСКОЙ РЕФЛЕКТОМЕТРИИ: ПРОБЛЕМЫ И ВОЗМОЖНЫЕ РЕШЕНИЯ
КОНТРОЛЬ ВОЛС МЕТОДОМ БРИЛЛЮЭНОВСКОЙ РЕФЛЕКТОМЕТРИИ: ПРОБЛЕМЫ И ВОЗМОЖНЫЕ РЕШЕНИЯ
Просмотры: 2279
DOI: 10.22184/2070-8963.2020.92.7-8.38.44
Рассмотрена актуальность проблемы контроля состояния волоконно-оптических линий в условиях температурно-деформационных воздействий. Смоделировано влияние погрешности калибровки сенсора ВОЛС и приборной погрешности на точность определения температуры и деформации. Продемонстрирован новый метод постобработки оптического спектра отраженного бриллюэновского излучения, позволяющий получать более точные значения температуры или деформации волокна. Результаты моделирования и эксперимента показали достоверность предлагаемых решений для одновременного контроля температуры и деформации оптического кабеля.
Рассмотрена актуальность проблемы контроля состояния волоконно-оптических линий в условиях температурно-деформационных воздействий. Смоделировано влияние погрешности калибровки сенсора ВОЛС и приборной погрешности на точность определения температуры и деформации. Продемонстрирован новый метод постобработки оптического спектра отраженного бриллюэновского излучения, позволяющий получать более точные значения температуры или деформации волокна. Результаты моделирования и эксперимента показали достоверность предлагаемых решений для одновременного контроля температуры и деформации оптического кабеля.
Теги: deformation of an fiber optic cables fiber optic lines reflected brillouin optical spectrum temperature-deformation effects бриллюэновское излучение волоконно-оптические линии деформация оптического кабеля оптический спектр температурно-деформационное воздействие
КОНТРОЛЬ ВОЛС методом бриллюэновской рефлектометрии: проблемы и возможные решения
А.Кривошеев, м.н.с. Пермского федерального исследовательского центра
Уральского отделения Российской академии наук (ПФИЦ УрО РАН) / antokri@ya.ru,
Ф.Барков, к.ф.-м.н., ст. науч. сотр. ПФИЦ УрО РАН / theodor1@mail.ru,
Ю.Константинов, к.т.н., заведующий лабораторией ПФИЦ УрО РАН / yuri.al.konstantinov@ro.ru
УДК 53.084, DOI: 10.22184/2070-8963.2020.92.7-8.38.44
Введение
В настоящее время достаточно актуален вопрос температурного и деформационного мониторинга волоконно-оптических линий связи (ВОЛС). Телекоммуникационные рефлектометры временной области (OTDR) позволяют контролировать проблемные участки ВОЛС, но не обеспечивают возможность предупреждения возможных обрывов оптического волокна.
Важная задача контроля состояния ВОЛС – именно обнаружение и предупреждение появления аварийных участков, например, с измененными температурой и натяжением оптического кабеля. Причин может быть множество, в частности, обледенение кабеля в холодное время года, несанкционированный доступ к ВОЛС и др.
Для контроля целостности различных технических систем хорошо зарекомендовали себя распределенные волоконно-оптические датчики. Особенно удобно, когда и контролируемый, и контролирующий объекты – один и тот же оптический кабель.
Разработчики активно ведут исследования по улучшению точности и быстродействия современных систем распределенного мониторинга. Не остается без внимания и метод бриллюэновской рефлектометрии (BOTDA/BOTDR). Точность определения температуры и деформации напрямую зависит от точности получения бриллюэновского спектра (BGS – Brillouin Gain Spectrum), а также сдвига частоты бриллюэновского спектра (BFS – Brillouin Frequency Shift).
Проблематика описанной задачи охватывает две области работ: первая – раздельное измерение температуры и деформации в одних и тех же точках кабеля, вторая – повышение точности и достоверности сенсорного сигнала в целом.
Решение задачи одновременного и независимого измерения температуры и деформации в оптическом волокне для многих систем мониторинга до сих пор остается достаточно нетривиальным. Современные методы получения информации о температурных и деформационных воздействиях могут включать использование двух различных волокон с записанными в них брэгговскими решетками [1, 2], внесение в конструкцию волоконного световода двух бессердцевинных волоконных тейперов разного диаметра [3], поляризационно-бриллюэновскую рефлектометрию [4, 5], фазочувствительную рефлектометрию [6], применение волокна с двумя сердцевинами [7], получение рефлектограмм для двух поляризационных мод в волокне, сохраняющем состояние поляризации вводимого излучения, с последующим корреляционным анализом рэлеевского спектра [8]. Характерными являются точностные характеристики в пределах нескольких градусов температуры (°С) и нескольких микрострейн (µε) деформации.
Помимо упомянутых методов, успешно применяется поляризационная бриллюэновская рефлектометрия. Например, в работах [4, 5] описывается метод анализа динамической акустической решетки в волокне, сохраняющем состояние поляризации, которая генерируется во время вынужденного бриллюэновского рассеяния. Метод сложен с точки зрения аппаратной части сенсора, но обеспечивает высокую точность определения температуры и деформации порядка 0,08 °С и 3 µε соответственно.
Резюмируя представленные методы, необходимо отметить, что для раздельного измерения двух величин нужны два независимых канала с разными откликами на воздействия – это могут быть как два различных телекоммуникационных волокна, так и один световод, поддерживающий две поляризационные моды (например, типа PANDA). Описанные ниже подходы проецируются на оба случая.
Что касается второй задачи (повышение достоверности сигнала сенсора), то точность системы в случае бриллюэновской рефлектометрии напрямую зависит от качества определения BFS. Исследованиям источников шумов, методов шумоподавления и получения более точных данных о BGS и BFS посвящены многочисленные работы [9–20]. Условно можно выделить два основных подхода к получению реального значения BFS.
Первый подход опирается на использование аппаратных мер, как с применением цифровых фильтров, так и обособленно, позволяющих определить максимум BGS при обработанном сигнале. Работы авторов [13, 14, 15 из MDPI] посвящены применению низкочастотной фильтрации сигнала, модуляции интенсивности излучения накачки и зондирующей волны, а также модуляции длины волны зондирующего сигнала. В работах [16, 17] описано успешное применение вейвлет-фильтрации. Эти достаточно эффективные для поставленных задач методы, в то же время слишком сложные для рядового использования распределенных датчиков в ряде приложений.
Второй подход – реконструкция спектра бриллюэновского рассеяния. В работе [18] представлен вариант на основе аппроксимации спектра при помощи известных аналитических функций. В работе [19] авторы описали использование лоренцевой функции для восстановления BGS, что позволило удвоить точность определения BFS по сравнению с традиционными методами. Тем не менее этот подход ограничен формой спектра рассеяния, поскольку возможные шумы, сбои оцифровки сигнала, фотон-фононные взаимодействия излучения в конкретных волокнах могут приводить к существенным искажениям формы BGS [20, 21].
Настоящая работа посвящена исследованию погрешностей, возникающих при анализе данных, получаемых при помощи рефлектометрических систем, на примере бриллюэновского анализатора, а также новому методу, направленному на борьбу с такими погрешностями.
Теоретические основы
Как уже отмечалось, направлению одновременного измерения температуры и деформации при помощи оптического волокна посвящено достаточное количество работ. Среди рассматриваемых подходов выделяется общая идея использования двух систем с различными чувствительностями к температуре и деформации.
С этой точки зрения интерес представляет волокно, сохраняющее поляризацию. Такое волокно способно выступать каналом одновременно для двух оптических волн, состояния поляризации которых повернуты относительно друг друга на 90°. При этом волны не будут взаимодействовать друг с другом. Конструкция таких волокон позволяет предположить, что две волны будут отличаться по чувствительности к внешним воздействиям.
Основы бриллюэновской рефлектометрии без использования поляризационной чувствительности изложены в [22].
В общем виде зависимость сдвига бриллюэновского спектра от температуры и деформации можно записать следующим образом:
Δf = a1ΔT + a2Δε, (1)
где Δf – величина сдвига бриллюэновского спектра по частоте; Δε – величина изменения деформации волокна; ΔT – величина изменения температуры волокна; a1, a2 – параметры, характеризующие чувствительность оптического волокна к воздействию температуры и деформации соответственно.
В случае поляризационной бриллюэновской рефлектометрии в распределенном анализе участвуют две спектральные компоненты излучения, соответствующие двум поляризационным осям волокна. При этом у каждой поляризационной оси оптического волокна своя чувствительность к воздействию температуры и деформации. В таком случае выражение (1) можно преобразовать следующим образом:
, (2)
где Δfx и Δfy – величины сдвига бриллюэновского спектра для каждой поляризационной оси.
Для нахождения коэффициентов чувствительности aij необходимо провести калибровочный эксперимент для конкретного оптического волокна, поскольку для различных волокон коэффициенты могут отличаться, что приведет к некорректным данным о состоянии оптического кабеля.
В оптическом волокне, сохраняющем состояние поляризации, как уже было сказано, возможно независимое распространение двух оптических волн, состояния поляризации которых повернуты относительно друг друга на 90°. Сохранение стабильности поляризации вводимого излучения поддерживается благодаря наведенному в волокне двулучепреломлению, которое фактически представляет собой разницу показателей преломления двух ортогональных осей в оптическом волокне. Эта разница обеспечивается конструкцией волоконного световода. Типичный порядок величины двулучепреломления для волокон, сохраняющих состояние поляризации, типа PANDA составляет 10–4.
Исходя из вышеуказанного, можно сделать вывод, что весомый вклад в измерение температуры и деформации вносит точность измерительной системы. Точность измерения бриллюэновских сдвигов fx и fy определяется рядом факторов [10, 12, 23]. При частотном сканировании спектр подвергается дискретизации, что заведомо вносит ошибку в определение бриллюэновского сдвига частоты. Помимо этого, даже при бесконечно малом шаге сканирования неточность определения максимума бриллюэновского спектра может быть вызвана зашумленностью сигнала. Для увеличения отношения сигнал/шум (SNR) используются различные приемы, например фильтрация низких частот, модуляция интенсивности волны накачки и зондирующей волны, модуляция длины волны зондирующего сигнала [14, 15, 24].
Наряду с инструментальными погрешностями известен еще один источник ошибок – различие между действительным и расчетным значениями чувствительности оптического волокна к внешним воздействиям. В таком случае даже идеальное определение величины бриллюэновских сдвигов fx и fy не даст достоверного результата об истинных температуре и деформации в оптическом волокне.
Оценка погрешности калибровочного измерения для разделения температуры и деформации
Для оценки влияния ошибок, обусловленных погрешностями, в расчетах был поставлен калибровочный эксперимент с целью нахождения коэффициентов чувствительности aij. Процедура калибровочного измерения предусматривала два этапа. На первом этапе волокно в свободной намотке прошло термоциклирование c одночасовой выдержкой на температурах 25, 30, 40, 50 и 60 °С для двух состояний поляризации. При этом для всей длины волокна снимались показания бриллюэновского сдвига при пространственном разрешении 0,1 м и шаге частотного сканирования в 5 МГц.
На втором этапе часть свободно намотанного волокна была помещена в стенд, обеспечивающий заданные продольные натяжения. На участке волокна 1 м при постоянной температуре окружающей среды 25 °С было осуществлено продольное растяжение с усилиями до 5 Н. На указанном участке волокна также фиксировались бриллюэновские сдвиги fx и fy .
На рис.1 представлена зависимость значений бриллюэновских сдвигов fx и fy от температуры. Зависимость для опыта с растяжением волокон имеет аналогичный характер, но иной наклон по отношению к оси абсцисс.
По тангенсам угла наклона кривых были получены значения коэффициентов матрицы:
. (2)
Эксперимент дал понимание того, что матрица (2) невырожденная, а значит, раздельное измерение температуры и деформации ВОЛС таким способом возможно. Предполагается, что при наличии двух разных волокон различие между соответствующими элементами матриц будет еще больше.
Затем было промоделировано влияние возможных погрешностей в определении коэффициентов aij на ошибку, возникающую при пересчете температуры и деформации. На рис.2 представлена зависимость ошибки в определении температуры от ошибки расчета коэффициентов чувствительности aij.
Зависимость погрешности определения деформации носит аналогичный характер.
В работе [25] авторы показали также вклад в погрешность измерения шага сканирования измерительного прибора и соотношение сигнал/шум получаемых спектров. По результатам моделирования, даже относительно небольшой шумовой вклад (SNR ~ 20dB) может приводить к значительным отклонениям показаний прибора. Эта часть погрешности закладывается уже на этапе присвоения калибровочных коэффициентов, а значит, ее вклада в погрешность в описываемом методе разделения температурных и деформационных воздействий не избежать. Второе весомое, но все же контролируемое, слагаемое общей погрешности – непосредственно погрешность измерения для разных методов определения бриллюэновских сдвигов – рассмотрим ниже.
Точность измерений и способ ее повышения
Предположим, что спектр, полученный в результате измерения, в виде дискретного набора выборок, состоящего из 2N+1 пар [f0 + i*Δf, Pi], где i – номер выборки в спектре, изменяющийся от 0 до 2N; f0 – минимальная частота спектра; Δf – шаг сканирования частоты, определяемый аппаратной дискретизацией; Pi – зарегистрированная плотность мощности обратного рассеяния на частоте fi = f0 + i*Δf (рис.3). Pi состоит из двух частей – собственно полезного сигнала и шума: Pi = Psi + Pni. При отсутствии шумовой составляющей, например при Pni = 0, нахождение максимальной мощности спектра дало бы центральную частоту fb с точностью до ошибки выборки, которая в среднем составляет Δf / 4. Однако в работе [25] показано, что наличие даже небольшого шума (SNR <20 дБ) приводит к тому, что погрешность в основном определяется шумом в спектре и слабо уменьшается по мере снижения Δf.
Для оценки эффективности метода была проведена серия численных экспериментов в диапазоне сканирования от 10 400 до 10 800 МГц. Для данного шага сканирования и соотношения SNR был смоделирован спектр в соответствии с выражением (3). При этом величина R задавалась случайным образом, имитируя шумовую составляющую при получении спектров бриллюэновских сдвигов.
, (3)
где w – ширина бриллюэновского спектра (при моделировании принимается за постоянную величину 40 МГц); Pn – амплитуда шумовой составляющей. Pin был выбран случайным образом равновероятно в диапазоне R:
. (4)
Для моделируемого спектра максимум был найден простым методом (по максимуму – FindMax), а также описанным способом (по методу "обратной корреляции"). Разница между фактическим значением BFS и найденными значениями была единичной погрешностью измерения для каждого метода. Затем была выбрана новая центральная частота для каждого спектра, и процесс повторился. Было сгенерировано более 10 тыс. спектров. Результаты для шага частотного сканирования 4 МГц приведены на рис.4.
Из полученной зависимости видно, что сделанное предположение верно – в целом при низких SNR корреляционный метод дает меньшую погрешность в определении значения сдвига бриллюэновской частоты. Точку пересечения C двух кривых на рис.4 (18 дБ) уместно называть "критическим SNR".
Зависимость "критического SNR" от шага частотного сканирования показана на рис.5. Когда SNR находится ниже критической линии, метод "обратной корреляции" обеспечивает более высокую точность. Из рис.5 следует, что чем меньше шаг сканирования, тем больше диапазон SNR, где корреляционный метод дает лучшие результаты.
В работе [26] представленный способ также сравнивался с методом восстановления спектра по лоренцевой функции и предварительно показал схожую точность регистрации физических величин.
Выводы
В работе рассмотрены и оценены различного рода погрешности, как измерительные, так и инструментальные, которые вносят вклад в точность получения информации о температуре и деформации. Раздельное измерение температуры и деформации с использованием двух измерительных каналов теоретически и практически реализуемо, но при этом большое внимание необходимо уделять калибровочному эксперименту, а именно нахождению калибровочных коэффициентов. Если нет возможности получения этих величин более точным способом, авторы предлагают прибегнуть к более длительному временному или пространственному накоплению. Первое возможно в условиях соблюдения стабильности всех внешних факторов, второе – в случае оптической однородности сенсора, контролирующего ВОЛС. Второе условие выглядит вполне реализуемым ввиду применения в ВОЛС довольно хорошо известных равномерных волокон, производимых способами OVD/VAD.
Разработанный авторами метод "обратной корреляции", предназначенный для получения точных значений сдвига бриллюэновского спектра с использованием рефлектометров и анализаторов, продемонстрировал свои преимущества перед простейшим способом расчета максимума спектра – как по моделируемым спектрам, так и по данным, полученным в эксперименте. Более того, метод также сравнивался с алгоритмом восстановления спектра по лоренцевой функции и предварительно показал схожую точность регистрации физических величин. Тот факт, что метод эффективен на реальных данных с достаточно высоким уровнем шума, использует типовой алгоритм свертки, скорее всего, указывает на его основную роль в потенциальном удешевлении рефлектометров и бриллюэновских анализаторов, поскольку становится возможным применение менее дорогостоящих фотоприемников и менее производительной вычислительной мощности (ввиду простоты реализации).
Проведенное исследование профинансировано в рамках выполнения государственного задания АААА-А19-119042590085-2.
ЛИТЕРАТУРА
Shishkin V., Terentyev V., Kharenko D., Dostovalov A., Wolf A., Simonov V., Fedotov M., Shienok A., Shelemba I., Babin S. Experimental Method of Temperature and Strain Discrimination in Polymer Composite Material by Embedded Fiber-Optic Sensors Based on Femtosecond-Inscribed FBGs // Journal of Sensors. 2016.
Sivanesan P., Sirkis J., Murata Y., and Buckley S. Optimal wavelength pair selection and accuracy analysis of dual fiber grating sensors for simultaneously measuring strain and temperature // Optical Engineering. 1 October 2002. 41(10).
Andre R.M., Biazoli C.R., Silva S.O., Marques M.B., Cordeiro C.M.B., and Frazao O. Strain-Temperature Discrimination Using Multimode Interference in Tapered Fiber // IEEE Photonics Technology Letters. 15 January 2013. V. 25. No. 2. P. 155–158.
Zou W., He Z., Hotate K. Complete discrimination of strain and temperature using Brillouin frequency shift and birefringence in a polarization-maintaining fiber // Optics Express. 2009. V. 17. Р. 1248–1255.
Hotate K., Zou W., Yamashita R. et al. Distributed discrimination of strain and temperature based on Brillouin dynamic grating in an optical fiber // Photonic Sens. 2013. V. 3. Р. 332–344.
Froggatt M., Gifford D., Kreger S., Wolfe M., and Soller B. Distributed Strain and Temperature Discrimination in Unaltered Polarization Maintaining Fiber // Optical Fiber Sensors. 23–27 October 2006.
Zaghloul M.A.S., Wang M., Milione G., Li M., Li S., Huang Y., Wang T., Chen K.P. Discrimination of Temperature and Strain in Brillouin Optical Time Domain Analysis Using a Multicore Optical Fiber // Sensors. 2018. V. 18. P. 1176.
Lu X., Soto M., Thévenaz L. Temperature-strain discrimination in distributed optical fiber sensing using phase-sensitive optical time-domain reflectometry // Optics Express. 2017. V. 25. P. 16059.
Urricelqui J., Soto M., Thévenaz L. Sources of noise in Brillouin optical time-domain analyzers // Proceedings of the 24th International Conference on Optical Fibre, Sensors, Curitiba, Brazil. 28 September 2015.
Zheng H., Fang Z., Wang Z., Lu B., Cao Y., Ye Q., Qu R., Cai H. Brillouin frequency shift of fiber distributed sensors extracted from noisy signals by quadratic fitting // Sensors. 2018. V. 18. P. 409.
Zaslawski S., Yang Z., Soto M., Thévenaz L. Impact of Fitting and Digital Filtering on Signal-to-Noise Ratio and Brillouin Frequency Shift Uncertainty of BOTDA Measurements. Optical Society of America: Washington, DC, USA, 2018.
Soto M., Thévenaz L. Modeling and evaluating the performance of Brillouin distributed optical fiber sensors // Optics Express. 2013. V. 21. Р. 31347–31366.
Feng C., Preussler S., Kadum J., Schneider T. Measurement accuracy enhancement via radio frequency filtering in distributed brillouin sensing // Sensors. 2019. V. 19. P. 2878.
Gyger F., Yang Z., Soto M., Yang F., Tow K., Thévenaz L. High Signal-to-Noise Ratio Stimulated Brillouin Scattering Gain Spectrum Measurement. – Optical Society of America: Washington, DC, USA. ThE69, 2018.
Urricelqui J., Sagues M., Loayssa A. Synthesis of Brillouin frequency shift profiles to compensate non-local effects and Brillouin induced noise in BOTDA sensors // Opt. Express. 2014. V. 22.
Feng Z., Gan J., Lv H., Cui L. Application of wavelet analysis in distributed optical fiber Brillouin temperature strain monitoring system // IOP conference series. Earth Environ. Sci. 2018. V. 189.
Zhang Z., Hu W., Yan J., Zhang P. The research of optical fiber Brillouin spectrum denoising based on wavelet transform and neural network // Proc. SPIE Int. Soc. Opt. Eng. 2013.
Minardo A., Caccavale M., Coscetta A., Esposito G., Matano F., Sacchi M., Somma R., Zeni G., Zeni L. Monitoring test of crack opening in volcanic tuff (Coroglio Cliff, Italy) using distributed optical fiber sensor. Geophysics: Principles, Applications and Emerging Technologies. Hauppauge, NY: Nova Science Publishers, 2016.
Feng C., Lu X., Preussler S., Schneider T. Gain spectrum engineering in distributed brillouin fiber sensors // Journal of Lightwave Technology. 2019. V. 37. Р. 5231–5237.
Jia X., Rao Y., Wang Z., Zhang W., Ran Z., Deng K., Yang Z. Non-local effect in Brillouin optical time-domain analyzer based on Raman amplification // Qiangjiguang Yu Lizishu. High power laser and particle Beams. 2012. V. 24. Р. 1667–1671.
Dragic P., Kucera C., Furtick J., Guerrier J., Hawkins T., Ballato J. Brillouin spectroscopy of a novel baria-doped silica glass optical fiber // Opt. Express. 2013. 21. Р. 10924–10941.
Листвин А.В., Листвин В.Н. Рефлектометрия оптических волокон. М.: ЛЕСАРарт, 2005.
Zaslawski S., Yang Z., Soto M.A., and Thévenaz L. Impact of Fitting and Digital Filtering on Signal-to-Noise Ratio and Brillouin Frequency Shift Uncertainty of BOTDA Measurements // Proc. 26th International Conference on Optical Fiber Sensors. 24–28 September 2018.
Feng C., Preussler S., Emad Kadum J., Schneider T. Measurement Accuracy Enhancement via Radio Frequency Filtering in Distributed Brillouin Sensing // Sensors. 2019. V. 19. P. 2878.
Barkov F.L., Yu A., Konstantinov V., Burdin V., Krivosheev A.I. Theoretical and experimental estimation of the accuracy in simultaneous distributed measurements of temperatures and strains in anisotropic optical fibers using polarization-brillouin reflectometry // Instrum. Exp. Tech. 2020. V. 20.
Barkov F.L., Konstantinov Y.A., Krivosheev A.I. A Novel Method of Spectra Processing for Brillouin Optical Time Domain Reflectometry // Fibers. 2020. V. 8. Р. 60.
А.Кривошеев, м.н.с. Пермского федерального исследовательского центра
Уральского отделения Российской академии наук (ПФИЦ УрО РАН) / antokri@ya.ru,
Ф.Барков, к.ф.-м.н., ст. науч. сотр. ПФИЦ УрО РАН / theodor1@mail.ru,
Ю.Константинов, к.т.н., заведующий лабораторией ПФИЦ УрО РАН / yuri.al.konstantinov@ro.ru
УДК 53.084, DOI: 10.22184/2070-8963.2020.92.7-8.38.44
Введение
В настоящее время достаточно актуален вопрос температурного и деформационного мониторинга волоконно-оптических линий связи (ВОЛС). Телекоммуникационные рефлектометры временной области (OTDR) позволяют контролировать проблемные участки ВОЛС, но не обеспечивают возможность предупреждения возможных обрывов оптического волокна.
Важная задача контроля состояния ВОЛС – именно обнаружение и предупреждение появления аварийных участков, например, с измененными температурой и натяжением оптического кабеля. Причин может быть множество, в частности, обледенение кабеля в холодное время года, несанкционированный доступ к ВОЛС и др.
Для контроля целостности различных технических систем хорошо зарекомендовали себя распределенные волоконно-оптические датчики. Особенно удобно, когда и контролируемый, и контролирующий объекты – один и тот же оптический кабель.
Разработчики активно ведут исследования по улучшению точности и быстродействия современных систем распределенного мониторинга. Не остается без внимания и метод бриллюэновской рефлектометрии (BOTDA/BOTDR). Точность определения температуры и деформации напрямую зависит от точности получения бриллюэновского спектра (BGS – Brillouin Gain Spectrum), а также сдвига частоты бриллюэновского спектра (BFS – Brillouin Frequency Shift).
Проблематика описанной задачи охватывает две области работ: первая – раздельное измерение температуры и деформации в одних и тех же точках кабеля, вторая – повышение точности и достоверности сенсорного сигнала в целом.
Решение задачи одновременного и независимого измерения температуры и деформации в оптическом волокне для многих систем мониторинга до сих пор остается достаточно нетривиальным. Современные методы получения информации о температурных и деформационных воздействиях могут включать использование двух различных волокон с записанными в них брэгговскими решетками [1, 2], внесение в конструкцию волоконного световода двух бессердцевинных волоконных тейперов разного диаметра [3], поляризационно-бриллюэновскую рефлектометрию [4, 5], фазочувствительную рефлектометрию [6], применение волокна с двумя сердцевинами [7], получение рефлектограмм для двух поляризационных мод в волокне, сохраняющем состояние поляризации вводимого излучения, с последующим корреляционным анализом рэлеевского спектра [8]. Характерными являются точностные характеристики в пределах нескольких градусов температуры (°С) и нескольких микрострейн (µε) деформации.
Помимо упомянутых методов, успешно применяется поляризационная бриллюэновская рефлектометрия. Например, в работах [4, 5] описывается метод анализа динамической акустической решетки в волокне, сохраняющем состояние поляризации, которая генерируется во время вынужденного бриллюэновского рассеяния. Метод сложен с точки зрения аппаратной части сенсора, но обеспечивает высокую точность определения температуры и деформации порядка 0,08 °С и 3 µε соответственно.
Резюмируя представленные методы, необходимо отметить, что для раздельного измерения двух величин нужны два независимых канала с разными откликами на воздействия – это могут быть как два различных телекоммуникационных волокна, так и один световод, поддерживающий две поляризационные моды (например, типа PANDA). Описанные ниже подходы проецируются на оба случая.
Что касается второй задачи (повышение достоверности сигнала сенсора), то точность системы в случае бриллюэновской рефлектометрии напрямую зависит от качества определения BFS. Исследованиям источников шумов, методов шумоподавления и получения более точных данных о BGS и BFS посвящены многочисленные работы [9–20]. Условно можно выделить два основных подхода к получению реального значения BFS.
Первый подход опирается на использование аппаратных мер, как с применением цифровых фильтров, так и обособленно, позволяющих определить максимум BGS при обработанном сигнале. Работы авторов [13, 14, 15 из MDPI] посвящены применению низкочастотной фильтрации сигнала, модуляции интенсивности излучения накачки и зондирующей волны, а также модуляции длины волны зондирующего сигнала. В работах [16, 17] описано успешное применение вейвлет-фильтрации. Эти достаточно эффективные для поставленных задач методы, в то же время слишком сложные для рядового использования распределенных датчиков в ряде приложений.
Второй подход – реконструкция спектра бриллюэновского рассеяния. В работе [18] представлен вариант на основе аппроксимации спектра при помощи известных аналитических функций. В работе [19] авторы описали использование лоренцевой функции для восстановления BGS, что позволило удвоить точность определения BFS по сравнению с традиционными методами. Тем не менее этот подход ограничен формой спектра рассеяния, поскольку возможные шумы, сбои оцифровки сигнала, фотон-фононные взаимодействия излучения в конкретных волокнах могут приводить к существенным искажениям формы BGS [20, 21].
Настоящая работа посвящена исследованию погрешностей, возникающих при анализе данных, получаемых при помощи рефлектометрических систем, на примере бриллюэновского анализатора, а также новому методу, направленному на борьбу с такими погрешностями.
Теоретические основы
Как уже отмечалось, направлению одновременного измерения температуры и деформации при помощи оптического волокна посвящено достаточное количество работ. Среди рассматриваемых подходов выделяется общая идея использования двух систем с различными чувствительностями к температуре и деформации.
С этой точки зрения интерес представляет волокно, сохраняющее поляризацию. Такое волокно способно выступать каналом одновременно для двух оптических волн, состояния поляризации которых повернуты относительно друг друга на 90°. При этом волны не будут взаимодействовать друг с другом. Конструкция таких волокон позволяет предположить, что две волны будут отличаться по чувствительности к внешним воздействиям.
Основы бриллюэновской рефлектометрии без использования поляризационной чувствительности изложены в [22].
В общем виде зависимость сдвига бриллюэновского спектра от температуры и деформации можно записать следующим образом:
Δf = a1ΔT + a2Δε, (1)
где Δf – величина сдвига бриллюэновского спектра по частоте; Δε – величина изменения деформации волокна; ΔT – величина изменения температуры волокна; a1, a2 – параметры, характеризующие чувствительность оптического волокна к воздействию температуры и деформации соответственно.
В случае поляризационной бриллюэновской рефлектометрии в распределенном анализе участвуют две спектральные компоненты излучения, соответствующие двум поляризационным осям волокна. При этом у каждой поляризационной оси оптического волокна своя чувствительность к воздействию температуры и деформации. В таком случае выражение (1) можно преобразовать следующим образом:
, (2)
где Δfx и Δfy – величины сдвига бриллюэновского спектра для каждой поляризационной оси.
Для нахождения коэффициентов чувствительности aij необходимо провести калибровочный эксперимент для конкретного оптического волокна, поскольку для различных волокон коэффициенты могут отличаться, что приведет к некорректным данным о состоянии оптического кабеля.
В оптическом волокне, сохраняющем состояние поляризации, как уже было сказано, возможно независимое распространение двух оптических волн, состояния поляризации которых повернуты относительно друг друга на 90°. Сохранение стабильности поляризации вводимого излучения поддерживается благодаря наведенному в волокне двулучепреломлению, которое фактически представляет собой разницу показателей преломления двух ортогональных осей в оптическом волокне. Эта разница обеспечивается конструкцией волоконного световода. Типичный порядок величины двулучепреломления для волокон, сохраняющих состояние поляризации, типа PANDA составляет 10–4.
Исходя из вышеуказанного, можно сделать вывод, что весомый вклад в измерение температуры и деформации вносит точность измерительной системы. Точность измерения бриллюэновских сдвигов fx и fy определяется рядом факторов [10, 12, 23]. При частотном сканировании спектр подвергается дискретизации, что заведомо вносит ошибку в определение бриллюэновского сдвига частоты. Помимо этого, даже при бесконечно малом шаге сканирования неточность определения максимума бриллюэновского спектра может быть вызвана зашумленностью сигнала. Для увеличения отношения сигнал/шум (SNR) используются различные приемы, например фильтрация низких частот, модуляция интенсивности волны накачки и зондирующей волны, модуляция длины волны зондирующего сигнала [14, 15, 24].
Наряду с инструментальными погрешностями известен еще один источник ошибок – различие между действительным и расчетным значениями чувствительности оптического волокна к внешним воздействиям. В таком случае даже идеальное определение величины бриллюэновских сдвигов fx и fy не даст достоверного результата об истинных температуре и деформации в оптическом волокне.
Оценка погрешности калибровочного измерения для разделения температуры и деформации
Для оценки влияния ошибок, обусловленных погрешностями, в расчетах был поставлен калибровочный эксперимент с целью нахождения коэффициентов чувствительности aij. Процедура калибровочного измерения предусматривала два этапа. На первом этапе волокно в свободной намотке прошло термоциклирование c одночасовой выдержкой на температурах 25, 30, 40, 50 и 60 °С для двух состояний поляризации. При этом для всей длины волокна снимались показания бриллюэновского сдвига при пространственном разрешении 0,1 м и шаге частотного сканирования в 5 МГц.
На втором этапе часть свободно намотанного волокна была помещена в стенд, обеспечивающий заданные продольные натяжения. На участке волокна 1 м при постоянной температуре окружающей среды 25 °С было осуществлено продольное растяжение с усилиями до 5 Н. На указанном участке волокна также фиксировались бриллюэновские сдвиги fx и fy .
На рис.1 представлена зависимость значений бриллюэновских сдвигов fx и fy от температуры. Зависимость для опыта с растяжением волокон имеет аналогичный характер, но иной наклон по отношению к оси абсцисс.
По тангенсам угла наклона кривых были получены значения коэффициентов матрицы:
. (2)
Эксперимент дал понимание того, что матрица (2) невырожденная, а значит, раздельное измерение температуры и деформации ВОЛС таким способом возможно. Предполагается, что при наличии двух разных волокон различие между соответствующими элементами матриц будет еще больше.
Затем было промоделировано влияние возможных погрешностей в определении коэффициентов aij на ошибку, возникающую при пересчете температуры и деформации. На рис.2 представлена зависимость ошибки в определении температуры от ошибки расчета коэффициентов чувствительности aij.
Зависимость погрешности определения деформации носит аналогичный характер.
В работе [25] авторы показали также вклад в погрешность измерения шага сканирования измерительного прибора и соотношение сигнал/шум получаемых спектров. По результатам моделирования, даже относительно небольшой шумовой вклад (SNR ~ 20dB) может приводить к значительным отклонениям показаний прибора. Эта часть погрешности закладывается уже на этапе присвоения калибровочных коэффициентов, а значит, ее вклада в погрешность в описываемом методе разделения температурных и деформационных воздействий не избежать. Второе весомое, но все же контролируемое, слагаемое общей погрешности – непосредственно погрешность измерения для разных методов определения бриллюэновских сдвигов – рассмотрим ниже.
Точность измерений и способ ее повышения
Предположим, что спектр, полученный в результате измерения, в виде дискретного набора выборок, состоящего из 2N+1 пар [f0 + i*Δf, Pi], где i – номер выборки в спектре, изменяющийся от 0 до 2N; f0 – минимальная частота спектра; Δf – шаг сканирования частоты, определяемый аппаратной дискретизацией; Pi – зарегистрированная плотность мощности обратного рассеяния на частоте fi = f0 + i*Δf (рис.3). Pi состоит из двух частей – собственно полезного сигнала и шума: Pi = Psi + Pni. При отсутствии шумовой составляющей, например при Pni = 0, нахождение максимальной мощности спектра дало бы центральную частоту fb с точностью до ошибки выборки, которая в среднем составляет Δf / 4. Однако в работе [25] показано, что наличие даже небольшого шума (SNR <20 дБ) приводит к тому, что погрешность в основном определяется шумом в спектре и слабо уменьшается по мере снижения Δf.
Для оценки эффективности метода была проведена серия численных экспериментов в диапазоне сканирования от 10 400 до 10 800 МГц. Для данного шага сканирования и соотношения SNR был смоделирован спектр в соответствии с выражением (3). При этом величина R задавалась случайным образом, имитируя шумовую составляющую при получении спектров бриллюэновских сдвигов.
, (3)
где w – ширина бриллюэновского спектра (при моделировании принимается за постоянную величину 40 МГц); Pn – амплитуда шумовой составляющей. Pin был выбран случайным образом равновероятно в диапазоне R:
. (4)
Для моделируемого спектра максимум был найден простым методом (по максимуму – FindMax), а также описанным способом (по методу "обратной корреляции"). Разница между фактическим значением BFS и найденными значениями была единичной погрешностью измерения для каждого метода. Затем была выбрана новая центральная частота для каждого спектра, и процесс повторился. Было сгенерировано более 10 тыс. спектров. Результаты для шага частотного сканирования 4 МГц приведены на рис.4.
Из полученной зависимости видно, что сделанное предположение верно – в целом при низких SNR корреляционный метод дает меньшую погрешность в определении значения сдвига бриллюэновской частоты. Точку пересечения C двух кривых на рис.4 (18 дБ) уместно называть "критическим SNR".
Зависимость "критического SNR" от шага частотного сканирования показана на рис.5. Когда SNR находится ниже критической линии, метод "обратной корреляции" обеспечивает более высокую точность. Из рис.5 следует, что чем меньше шаг сканирования, тем больше диапазон SNR, где корреляционный метод дает лучшие результаты.
В работе [26] представленный способ также сравнивался с методом восстановления спектра по лоренцевой функции и предварительно показал схожую точность регистрации физических величин.
Выводы
В работе рассмотрены и оценены различного рода погрешности, как измерительные, так и инструментальные, которые вносят вклад в точность получения информации о температуре и деформации. Раздельное измерение температуры и деформации с использованием двух измерительных каналов теоретически и практически реализуемо, но при этом большое внимание необходимо уделять калибровочному эксперименту, а именно нахождению калибровочных коэффициентов. Если нет возможности получения этих величин более точным способом, авторы предлагают прибегнуть к более длительному временному или пространственному накоплению. Первое возможно в условиях соблюдения стабильности всех внешних факторов, второе – в случае оптической однородности сенсора, контролирующего ВОЛС. Второе условие выглядит вполне реализуемым ввиду применения в ВОЛС довольно хорошо известных равномерных волокон, производимых способами OVD/VAD.
Разработанный авторами метод "обратной корреляции", предназначенный для получения точных значений сдвига бриллюэновского спектра с использованием рефлектометров и анализаторов, продемонстрировал свои преимущества перед простейшим способом расчета максимума спектра – как по моделируемым спектрам, так и по данным, полученным в эксперименте. Более того, метод также сравнивался с алгоритмом восстановления спектра по лоренцевой функции и предварительно показал схожую точность регистрации физических величин. Тот факт, что метод эффективен на реальных данных с достаточно высоким уровнем шума, использует типовой алгоритм свертки, скорее всего, указывает на его основную роль в потенциальном удешевлении рефлектометров и бриллюэновских анализаторов, поскольку становится возможным применение менее дорогостоящих фотоприемников и менее производительной вычислительной мощности (ввиду простоты реализации).
Проведенное исследование профинансировано в рамках выполнения государственного задания АААА-А19-119042590085-2.
ЛИТЕРАТУРА
Shishkin V., Terentyev V., Kharenko D., Dostovalov A., Wolf A., Simonov V., Fedotov M., Shienok A., Shelemba I., Babin S. Experimental Method of Temperature and Strain Discrimination in Polymer Composite Material by Embedded Fiber-Optic Sensors Based on Femtosecond-Inscribed FBGs // Journal of Sensors. 2016.
Sivanesan P., Sirkis J., Murata Y., and Buckley S. Optimal wavelength pair selection and accuracy analysis of dual fiber grating sensors for simultaneously measuring strain and temperature // Optical Engineering. 1 October 2002. 41(10).
Andre R.M., Biazoli C.R., Silva S.O., Marques M.B., Cordeiro C.M.B., and Frazao O. Strain-Temperature Discrimination Using Multimode Interference in Tapered Fiber // IEEE Photonics Technology Letters. 15 January 2013. V. 25. No. 2. P. 155–158.
Zou W., He Z., Hotate K. Complete discrimination of strain and temperature using Brillouin frequency shift and birefringence in a polarization-maintaining fiber // Optics Express. 2009. V. 17. Р. 1248–1255.
Hotate K., Zou W., Yamashita R. et al. Distributed discrimination of strain and temperature based on Brillouin dynamic grating in an optical fiber // Photonic Sens. 2013. V. 3. Р. 332–344.
Froggatt M., Gifford D., Kreger S., Wolfe M., and Soller B. Distributed Strain and Temperature Discrimination in Unaltered Polarization Maintaining Fiber // Optical Fiber Sensors. 23–27 October 2006.
Zaghloul M.A.S., Wang M., Milione G., Li M., Li S., Huang Y., Wang T., Chen K.P. Discrimination of Temperature and Strain in Brillouin Optical Time Domain Analysis Using a Multicore Optical Fiber // Sensors. 2018. V. 18. P. 1176.
Lu X., Soto M., Thévenaz L. Temperature-strain discrimination in distributed optical fiber sensing using phase-sensitive optical time-domain reflectometry // Optics Express. 2017. V. 25. P. 16059.
Urricelqui J., Soto M., Thévenaz L. Sources of noise in Brillouin optical time-domain analyzers // Proceedings of the 24th International Conference on Optical Fibre, Sensors, Curitiba, Brazil. 28 September 2015.
Zheng H., Fang Z., Wang Z., Lu B., Cao Y., Ye Q., Qu R., Cai H. Brillouin frequency shift of fiber distributed sensors extracted from noisy signals by quadratic fitting // Sensors. 2018. V. 18. P. 409.
Zaslawski S., Yang Z., Soto M., Thévenaz L. Impact of Fitting and Digital Filtering on Signal-to-Noise Ratio and Brillouin Frequency Shift Uncertainty of BOTDA Measurements. Optical Society of America: Washington, DC, USA, 2018.
Soto M., Thévenaz L. Modeling and evaluating the performance of Brillouin distributed optical fiber sensors // Optics Express. 2013. V. 21. Р. 31347–31366.
Feng C., Preussler S., Kadum J., Schneider T. Measurement accuracy enhancement via radio frequency filtering in distributed brillouin sensing // Sensors. 2019. V. 19. P. 2878.
Gyger F., Yang Z., Soto M., Yang F., Tow K., Thévenaz L. High Signal-to-Noise Ratio Stimulated Brillouin Scattering Gain Spectrum Measurement. – Optical Society of America: Washington, DC, USA. ThE69, 2018.
Urricelqui J., Sagues M., Loayssa A. Synthesis of Brillouin frequency shift profiles to compensate non-local effects and Brillouin induced noise in BOTDA sensors // Opt. Express. 2014. V. 22.
Feng Z., Gan J., Lv H., Cui L. Application of wavelet analysis in distributed optical fiber Brillouin temperature strain monitoring system // IOP conference series. Earth Environ. Sci. 2018. V. 189.
Zhang Z., Hu W., Yan J., Zhang P. The research of optical fiber Brillouin spectrum denoising based on wavelet transform and neural network // Proc. SPIE Int. Soc. Opt. Eng. 2013.
Minardo A., Caccavale M., Coscetta A., Esposito G., Matano F., Sacchi M., Somma R., Zeni G., Zeni L. Monitoring test of crack opening in volcanic tuff (Coroglio Cliff, Italy) using distributed optical fiber sensor. Geophysics: Principles, Applications and Emerging Technologies. Hauppauge, NY: Nova Science Publishers, 2016.
Feng C., Lu X., Preussler S., Schneider T. Gain spectrum engineering in distributed brillouin fiber sensors // Journal of Lightwave Technology. 2019. V. 37. Р. 5231–5237.
Jia X., Rao Y., Wang Z., Zhang W., Ran Z., Deng K., Yang Z. Non-local effect in Brillouin optical time-domain analyzer based on Raman amplification // Qiangjiguang Yu Lizishu. High power laser and particle Beams. 2012. V. 24. Р. 1667–1671.
Dragic P., Kucera C., Furtick J., Guerrier J., Hawkins T., Ballato J. Brillouin spectroscopy of a novel baria-doped silica glass optical fiber // Opt. Express. 2013. 21. Р. 10924–10941.
Листвин А.В., Листвин В.Н. Рефлектометрия оптических волокон. М.: ЛЕСАРарт, 2005.
Zaslawski S., Yang Z., Soto M.A., and Thévenaz L. Impact of Fitting and Digital Filtering on Signal-to-Noise Ratio and Brillouin Frequency Shift Uncertainty of BOTDA Measurements // Proc. 26th International Conference on Optical Fiber Sensors. 24–28 September 2018.
Feng C., Preussler S., Emad Kadum J., Schneider T. Measurement Accuracy Enhancement via Radio Frequency Filtering in Distributed Brillouin Sensing // Sensors. 2019. V. 19. P. 2878.
Barkov F.L., Yu A., Konstantinov V., Burdin V., Krivosheev A.I. Theoretical and experimental estimation of the accuracy in simultaneous distributed measurements of temperatures and strains in anisotropic optical fibers using polarization-brillouin reflectometry // Instrum. Exp. Tech. 2020. V. 20.
Barkov F.L., Konstantinov Y.A., Krivosheev A.I. A Novel Method of Spectra Processing for Brillouin Optical Time Domain Reflectometry // Fibers. 2020. V. 8. Р. 60.
Отзывы читателей