COMBINED NEURAL NETWORK METHOD FOR DETERMINING THE BRILLOUIN SPECTRUM MAXIMUM IN DISTRIBUTED FIBER-OPTIC SENSORS
- Authors: Krivosheev A.I1, Konstantinov Y.A1, Pervadchuk V.P2, Barkov F.L1
- Affiliations:
- Perm Federal Research Center of the Ural Branch of the RAS
- Perm National Research Polytechnic University
- Issue: No 3 (2021)
- Pages: 95-106
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/amcs/article/view/2070
- DOI: https://doi.org/10.15593/2499-9873/2021.3.05
- Cite item
Abstract
Distributed fiber-optic sensors are becoming an increasingly common solution for monitoring and diagnosing extended information transmission lines, industrial devices, buildings and devices. One of the directions - Brillouin reflectometry, which allows diagnosing a fiber line for changes in ambient temperature or mechanical deformation, is becoming increasingly popular for engineers and researchers. However, modern standards impose increasingly strict requirements on diagnostic systems in terms of the accuracy of determined parameters. For Brillouin reflectometry, where the value of the environmental parameters is determined by the position of the maximum of the Brillouin gain spectrum, the task of more accurately determining this maximum becomes the main one. The paper considers modern computer-computational methods for detecting the maximum of the Brillouin gain spectrum in an optical fiber. The authors note that imperfections in the shape of the optical spectrum, such as the signal-to-noise ratio, as well as possible digital defects that occur during digitization, can significantly impair the accuracy of the system. The authors consider three approaches to detecting the maximum of the spectrum: the classical method of Lorentz curve fitting, the method of cross-correlation with the ideal Lorentzian function, as well as the method of inverse correlation developed by the authors earlier. To combine the results of the work of the three methods, a neural network model was developed that accepts the input data of each method, together with the parameters of noise and distortion of the spectrum. The model is presented in the form of a four-layer perceptron with two hidden layers. As a result, the authors achieved an increase in the accuracy of determining the position of the maximum of the spectrum by 10% on the model data.
Full Text
Исследовательский интерес к бриллюэновской рефлектометрии временной области с точки зрения распределенного волоконного зондирования значительно вырос со времени ее открытия [1]. Распределенное волоконное зондирование позволяет получать непрерывную информацию о состоянии волоконной линии длиной до нескольких километров c использованием одного блока обработки. Это является большим преимуществом перед локальными [2] и квазираспределенными [3, 4] системами мониторинга. В основе бриллюэновской рефлектометрии лежит вынужденное рассеяние Мандельштама - Бриллюэна (ВРМБ) в оптическом волокне, которое позволяет получить информацию о температуре и деформации с высокой точностью и разрешением [5-7]. Суть метода заключается во взаимодействии импульсной волны накачки со встречной непрерывной зондирующей волной. Когда разность частот между волной накачки и зондирующей совпадает со значением сдвига бриллюэновской частоты оптического волокна, который составляет порядка 10 ГГц для кварцевых волокон, зондирующая волна испытывает усиление за счет ВРМБ. Во время внешних возмущений, таких как температура или деформация, бриллюэновский сдвиг в этой точке волокна будет смещен. Величина этого смещения несет в себе информацию о величине изменения температуры или деформации, воздействующей на оптическое волокно [8, 9]. Спектр бриллюэновского усиления рассеянного в волокне излучения может быть восстановлен путем сканирования частоты между волной накачки и зондирующей волной. Затем распределение бриллюэновского сдвига вдоль волоконной линии можно вычислить путем определения центральной частоты в спектре усиления. Традиционно в системах мониторинга на основе бриллюэновского рассеяния центральную частотную компоненту в спектре усиления определяют путем аппроксимации функцией Лоренца, модель которой задается уравнением [10]: (1) где - значение бриллюэновского сдвига, - ширина спектра усиления, - амплитуда спектра усиления. Такой метод сопряжен с некоторыми трудностями. К ним относятся высокие вычислительные затраты, а также тот факт, что точность метода определяется заданным изначально параметром точности системы, по достижении которой цель аппроксимации будет достигнута [11, 12]. Предложены также методы, не требующие задания первоначальных параметров. К ним относятся методы кросс-корреляции [12, 13], использующие идеальную функцию Лоренца, заданную по формуле (1), для вычисления функции корреляции между измеренным спектром усиления и идеальной кривой Лоренца. Тем не менее для этого все равно требуется генерация дополнительного спектра. Кроме того, метод остается чувствительным к шумам системы и величине шага частотного сканирования. Недавно авторами был предложен новый метод [14], основанный на методе кросс-корреляции, но не требующий никаких дополнительных вводимых величин. Метод, названный методом обратной корреляции, использует тот же спектр, полученный в результате измерения, однако отраженный по оси абсцисс. Авторами было проведено исследование точности данного метода [15] в сравнении с традиционными аналогами. Показано, что метод имеет преимущества при наличии весомых несовершенств в измеренном спектре, таких как шумовая составляющая, а также различные искажения в спектре, которые могут быть вызваны дефектами оцифровки. На сегодняшний день известно множество исследований по внедрению машинного обучения в бриллюэновские распределенные датчики в качестве альтернативы традиционным методам получения информации о температуре и деформации. Машинное обучение было успешно применено в таких сферах, как обработка изображений, распознавание образов, классификации и оптики. Машинное обучение зарекомендовало себя благодаря способности точно предсказывать решение, имея лишь несколько входных данных, которые было бы слишком сложно организовать в простую логическую цепочку. Модель строится и обучается, используя известный набор данных, как входных, так и выходных. В области бриллюэновской рефлектометрии исследуются такие направления машинного обучения, как искусственные нейронные сети [16], глубинные нейронные сети [17], а также метод опорных векторов [18]. Однако, несмотря на достоинства таких подходов, классические методы обработки информации, полученной с помощью бриллюэновских датчиков, более распространены в типовых аттестованных системах, использующихся в производственных и лабораторных целях. Такие методы имеют как преимущества, так и недостатки в зависимости от условий проведения измерения. При этом каждый несет в себе часть полезной информации об истинном положении центральной компоненты бриллюэновского спектра. Данная работа посвящена использованию нейронной сети для объединения информации, полученной с бриллюэновского анализатора классическими способами, с целью увеличения точности определения максимума спектра бриллюэновского сдвига. В качестве входных параметров сети были использованы результаты трех методов поиска максимума спектра, а также значение отношения сигнал/шум спектра. Другой тип цифрового дефекта, помимо шума, хоть и нечасто встречающийся в современных датчиках, - локальное пропадание полезного сигнала спектра. Этот тип дефекта может иметь различную протяженность, поэтому ширина такого дефекта, в случае его наличия в спектре, также включена в набор входных параметров нейронной сети. Для обучения сети было сгенерировано 2000 спектров с различными соотношениями сигнал/шум, центральной частотой, а также положением и протяженностью цифрового дефекта. Соотношение объема данных при обучении и тестировании сети составило 1500:500 соответственно. Проектирование и верификация работоспособности сети проводились на задаче распознавания рукописного набора цифр базы данных MNIST. Точность распознавания при этом составила более 97 %, что соответствует результатам, представленным в литературе [19]. Путем последовательного перебора различных параметров нейронной сети было обнаружено, что наилучший результат определения максимума бриллюэновского спектра получается при использовании четырехслойного персептрона с двумя скрытыми слоями по 50 и 25 нейронов соответственно. В качестве функции активации нейронов использовалась сигмоидальная функция [20]. Конечная структура сети представлена на рис. 1. Сеть состоит из пяти входных нейронов, двух скрытых слоев и одного выходного нейрона. Рис. 1. Структура нейронной сети: BWC - результат метода обратной корреляции; CCM - классического метода кросс-корреляции; LCF - аппроксимация лоренцевой функцией; SNR - отношение сигнал/шум; Δart - ширина области дефекта Обучение сети производилось на сгенерированных спектрах по принципу обратного распространения ошибки. При этом способе обучения выходное значение сети сравнивается с правильным значением, а ошибка, представляющая собой разницу этих значений, возвращается через сеть. При этом корректируются значения весов связей между нейронами для последующего уменьшения ошибки. В связи с тем, что входные параметры сети имеют различные порядки и единицы измерения, каждый параметр был индивидуально нормализован в пределах от 0 до 1. На рис. 2 представлен график динамики ошибки сети в процессе обучения. Рис. 2. Динамика изменения ошибки определения максимума бриллюэновского спектра при помощи нейронной сети Оценка точности методов производилась путем расчета стандартного отклонения ошибки определения максимума спектра для каждого метода. Сравнительные характеристики методов расчета максимума для смоделированных спектров приведены в таблице. В результате применения такого подхода удалось добиться снижения ошибки определения максимума бриллюэновского спектра на величину порядка 10 %. Этот результат можно считать существенным, принимая во внимание то, что к спектру уже были применены современные методы поиска максимума, для которых дальнейшее увеличение точности повлекло значительное увеличение необходимой вычислительной мощности системы. Сравнение результатов расчета № п/п Метод Стандартное отклонение ошибки 1 BWC 0,0072 2 CCM 0,0071 3 LCF 0,0086 4 Комбинированный (нейросеть) 0,0064 Представленная работа демонстрирует потенциальную возможность увеличения точности классических распределенных систем бриллюэновского мониторинга без использования сложных алгоритмов обработки, а также без усложнения аппаратной части сенсора. Это открывает возможности для использования типового аттестованного оборудования для получения более точных сведений о состоянии волоконной линии. Помимо этого, такой принцип обработки позволит удешевить создание лабораторных систем, поскольку становится возможным использование менее дорогих фотоприемных устройств и вычислительной техники.About the authors
A. I Krivosheev
Perm Federal Research Center of the Ural Branch of the RAS
Yu. A Konstantinov
Perm Federal Research Center of the Ural Branch of the RAS
V. P Pervadchuk
Perm National Research Polytechnic University
F. L Barkov
Perm Federal Research Center of the Ural Branch of the RAS
References
- Horiguchi T., Tateda M. BOTDA-nondestructive measurement of single-mode optical fiber attenuation characteristics using Brillouin interaction: theory // Journal of Lightwave Technology. - 1989. - Vol. 7, iss. 8. - P. 1170-1176. doi: 10.1109/50.32378
- Behavior of femtosecond laser-induced eccentric fiber bragg gratings at very high temperatures / H. Chikh-Bled, K. Chah, Álvaro González-Vila, B. Lasri, C. Caucheteur // Optics Letters. - 2016. - Vol. 41, iss. 17. - P. 4048-4051. doi: 10.1364/OL.41.004048
- Fast and stable gratings inscription in pofs made of different materials with pulsed 248 nm krf laser / C.A.F. Marques, R. Min, A.L. Junior, P. Antunes, A. Fasano, G. Woyessa, K. Nielsen, H.K. Rasmussen, B. Ortega, O. Bang // Optics Express. - 2018. - Vol. 26, iss. 2 - P. 2013-2022. doi: 10.1364/OE.26.002013
- Compensation method for temperature cross-sensitivity in transverse force applications with FBG sensors in POFs / A.G. Leal-Junior, A. Theodosiou, C. Marques, M.J. Pontes, K. Kalli, A. Frizera // Journal of Lightwave Technology. - 2018. - Vol. 36, iss. 17 - P. 3660-3665. doi: 10.1109/JLT.2018.2848704
- Galindez-Jamioy C.A., López-Higuera J.M. Brillouin distributed fiber sensors: an overview and applications // Journal of Sensors. - 2012. - Vol. 2012. - Art. 204121. doi: 10.1155/2012/204121
- Zan M.S.D.B., Horiguchi T. A dual golay complementary pair of sequences for improving the performance of phase-shift pulse BOTDA fiber sensor // Journal of Lightwave Technology. - 2012. - Vol. 30, iss. 21. - P. 3338-3356. doi: 10.1109/JLT.2012.2215578
- Bao X., Chen L., Recent progress in distributed fiber optic sensors // Sensors. - 2012. - Vol. 12, iss. 12. - P. 8601-8639. doi: 10.3390/s120708601
- Kurashima T., Horiguchi T., Tateda M. Distributed-temperature sensing using stimulated Brillouin scattering in optical silica fibers // Optics Letters. - 1990. - Vol. 15, iss. 18. - P. 1038. doi: 10.1364/OL.15.001038
- Nikles M., Thevenaz L., Robert P.A. Measurement of the distributed-Brillouin gain spectrum in optical fibers by using a single laser source // OSA Technical Digest Series (Optical Society of America) [Conference on Optical Fiber Communication, San Jose, California United States, 20 February 1994]. - 1994. - Vol. 4. - Art. WF1. doi: 10.1364/OFC.1994.WF1
- Soto M.A., Thévenaz L. Modeling and evaluating the performance of Brillouin distributed optical fiber sensors // Optics Express. - 2013. - Vol. 21, iss. 25. - P. 31347. doi: 10.1364/OE.21.031347
- Li C., Li Y. Fitting of Brillouin spectrum based on LabVIEW // 2009 5th International Conference on Wireless Communications, Networking and Mobile Computing, 24-26 September 2009, Beijing, China. - IEEE, 2009. - P. 1-4. doi: 10.1109/WICOM.2009.5303692
- Farahani M.A., Castillo-Guerra E., Colpitts B.G. Accurate estimation of Brillouin frequency shift in Brillouin optical time domain analysis sensors using cross correlation // Optics Letters. - 2011. - Vol. 36, iss. 21. - P. 4275. doi: 10.1364/OL.36.004275
- Farahani M.A., Castillo-Guerra E., Colpitts B.G. A detailed evaluation of the correlation-based method used for estimation of the Brillouin frequency shift in BOTDA sensors // IEEE Sensors Journal. - 2013. - Vol. 13, iss. 12. - P. 4589-4598. doi: 10.1109/JSEN.2013.2271254
- Barkov F.L., Konstantinov Y.A., Krivosheev A.I. A novel method of spectra processing for Brillouin optical time domain reflectometry // Fibers. - 2020. - Vol. 8, iss. 9. - Р. 60. doi: 10.3390/fib8090060
- Сравнительный анализ точности определения сдвига бриллюэновской частоты в экстремально зашумленных спектрах различными корреляционными методами / А.И. Кривошеев, Ю.А. Константинов, Ф.Л. Барков, В.П. Первадчук // Приборы и техника эксперимента. - 2021. - № 2. - C. 84-89. doi: 10.31857/S0032816221050062
- Temperature extraction in Brillouin optical time-domain analysis sensors using principal component analysis based pattern recognition / A.K. Azad, F.N. Khan, W.H. Alarashi, N. Guo, A.P.T. Lau, C. Lu // Optics Express. - 2017. - Vol. 25, iss. 14. - P. 16534-16549. doi: 10.1364/OE.25.016531
- Deep neural networks assisted BOTDA for simultaneous temperature and strain measurement with enhanced accuracy / B. Wang, L. Wang, N. Guo, Z. Zhao, C. Yu, C. Lu // Optics Express. - 2019. - Vol. 27, iss. 3. - P. 2530. doi: 10.1364/OE.27.002530
- Support vector machine based differential pulse-width pair Brillouin optical time domain analyzer / H. Wu, L. Wang, Z. Zhao, C. Shu, C. Lu // IEEE Photonics Journal. - 2018. - Vol. 10, iss. 4. - P. 1-11. doi: 10.1109/JPHOT.2018.2858235
- Тарик Р. Создаем нейронную сеть: пер. с англ. - СПб.: Альфа-книга, 2017. - 272 с.
- Cybenko G. Approximiation by superpositions of sigmoidal functions // Mathematics of Control Signals and Systems. - 1989. - Vol. 2, no. 4. - P. 303-314.
Statistics
Views
Abstract - 69
PDF (Russian) - 35
Refbacks
- There are currently no refbacks.