Научная статья на тему 'ИЗМЕРЕНИЕ ВЕЛИЧИНЫ ДВУЛУЧЕПРЕЛОМЛЕНИЯ ВОЛОКНА ТИПА PANDA С ПОМОЩЬЮ ИНТЕРФЕРОМЕТРА САНЬЯКА'

ИЗМЕРЕНИЕ ВЕЛИЧИНЫ ДВУЛУЧЕПРЕЛОМЛЕНИЯ ВОЛОКНА ТИПА PANDA С ПОМОЩЬЮ ИНТЕРФЕРОМЕТРА САНЬЯКА Текст научной статьи по специальности «Нанотехнологии»

CC BY
55
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
двулучепреломление / волокна с сохранением поляризации / волокно Panda / интерферометр Саньяка / перестраиваемый лазер / birefringence / polarization maintaining fibers / PANDA fiber / Sagnac interferometer / tunable laser

Аннотация научной статьи по нанотехнологиям, автор научной работы — Константин Александрович Овчинников, Виктор Владимирович Криштоп, Даниил Георгиевич Гилев, Виталий Александрович Максименко, Анатолий Викторович Перминов

Представлена методика измерения величины двулучепреломления оптических волокон с применением интерферометра Саньяка и источника излучения c перестраиваемой длиной волны. Рассмотрены физические принципы, лежащие в основе представленной методики. Для устранения влияния нелинейности скорости перестройки длины волны источника излучения на точность измерений применен вспомогательный интерферометр с известными параметрами. Полученные значения двулучепреломления волокна совпадают с паспортными, что подтверждает эффективность предложенной методики измерения.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по нанотехнологиям , автор научной работы — Константин Александрович Овчинников, Виктор Владимирович Криштоп, Даниил Георгиевич Гилев, Виталий Александрович Максименко, Анатолий Викторович Перминов

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MEASURING THE BIREFRINGENCE VALUE OF A PANDA FIBER USING A SAGNAC INTERFEROMETER

A technique for measuring the birefringence of optical fibers using a Sagnac interferometer and a radiation source with a tunable wavelength is presented. Physical principles that form the basis the presented methodology are considered. To eliminate the influence of the nonlinearity of the radiation source wavelength tuning rate on the measurement accuracy, an auxiliary interferometer with known parameters is used. The obtained values of measured fiber birefringence coincide with the passport values, which confirms the proposed measurement technique effectiveness.

Текст научной работы на тему «ИЗМЕРЕНИЕ ВЕЛИЧИНЫ ДВУЛУЧЕПРЕЛОМЛЕНИЯ ВОЛОКНА ТИПА PANDA С ПОМОЩЬЮ ИНТЕРФЕРОМЕТРА САНЬЯКА»

ОПТИЧЕСКИЕ И ОПТИКО-ЭЛЕКТРОННЫЕ ПРИБОРЫ И КОМПЛЕКСЫ OPTICAL AND OPTO-ELECTRONIC INSTRUMENTS AND COMPLEXES

УДК 535.321.9

DOI: 10.17586/0021-3454-2023-66-10-845-851

ИЗМЕРЕНИЕ ВЕЛИЧИНЫ ДВУЛУЧЕПРЕЛОМЛЕНИЯ ВОЛОКНА ТИПА PANDA С ПОМОЩЬЮ ИНТЕРФЕРОМЕТРА САНЬЯКА

К. А. Овчинников u*, В. В. Криштоп1,2,3, Д. Г. Гилев1,2, В. А. Максименко2, А. В. Перминов2

'Пермская научно-производственная приборостроительная компания

Пермь, Россия [email protected]

2 Пермский национальный исследовательский политехнический университет

Пермь, Россия

3 Пермский государственный национальный исследовательский университет

Пермь, Россия

Аннотация. Представлена методика измерения величины двулучепреломления оптических волокон с применением интерферометра Саньяка и источника излучения c перестраиваемой длиной волны. Рассмотрены физические принципы, лежащие в основе представленной методики. Для устранения влияния нелинейности скорости перестройки длины волны источника излучения на точность измерений применен вспомогательный интерферометр с известными параметрами. Полученные значения двулучепреломления волокна совпадают с паспортными, что подтверждает эффективность предложенной методики измерения.

Ключевые слова: двулучепреломление, волокна с сохранением поляризации, волокно Panda, интерферометр Саньяка, перестраиваемый лазер

Благодарности: исследования выполнены при поддержке Министерства науки и высшего образования Российской Федерации (проект № FSNM-2023-0005).

Ссылка для цитирования: Овчинников К. А., Криштоп В. В., Гилев Д. Г., Максименко В. А., Перминов А. В. Измерение величины двулучепреломления волокна типа Panda с помощью интерферометра Саньяка // Изв. вузов. Приборостроение. 2023. Т. 66, № 10. С. 845—851. DOI: 10.17586/0021-3454-2023-66-10-845-851.

MEASURING THE BIREFRINGENCE VALUE OF A PANDA FIBER USING A SAGNAC INTERFEROMETER

К. А. Ovchinnikov1,2*, V. V. Krishtop1,2,3, D. G. Gilev1,2, V. А. Maksimenko2, А. V. Perminov2

1 Perm Research and Production Instrument-Making Company, Perm, Russia *[email protected] 2 Perm National Research Polytechnic University, Perm, Russia 3 Perm State University, Perm, Russia

Abstract. A technique for measuring the birefringence of optical fibers using a Sagnac interferometer and a radiation source with a tunable wavelength is presented. Physical principles that form the basis the presented methodology are considered. To eliminate the influence of the nonlinearity of the radiation source wavelength tuning rate on the measurement accuracy, an auxiliary interferometer with known parameters is used. The obtained values of measured fiber birefringence coincide with the passport values, which confirms the proposed measurement technique effectiveness.

© Овчинников К. А., Криштоп В. В., Гилев Д. Г., Максименко В. А., Перминов А. В., 2023

Keywords: birefringence, polarization maintaining fibers, PANDA fiber, Sagnac interferometer, tunable laser

Acknowledment: research was carried out with the support of the Ministry of Science and Higher Education of the Russian Federation (project No. FSNM-2023-0005).

For citation: Ovchinnikov K. A., Krishtop V. V., Gilev D. G., Maksimenko V. A., Perminov A. V. Measuring the birefringence value of a PANDA fiber using a Sagnac interferometer. Journal of Instrument Engineering. 2023. Vol. 66, N 10. P. 845—851 (in Russian). DOI: 10.17586/0021-3454-2023-66-10-845-851.

Введение. С развитием волоконно-оптических датчиков большую популярность получили PMF (polarization maintaining fiber) — оптические волокна, сохраняющие поляризацию, поскольку они способны передавать линейно-поляризованное излучение по всей длине волновода. Способность сохранять поляризацию в таких волокнах достигается за счет создания асимметричной структуры, которая приводит к возникновению двулучепреломления [1—3].

Двулучепреломление характеризует способность удерживать состояние поляризации при распространении излучения. Величина двулучепреломления определяет важные параметры PMF, такие как длина биений и дисперсия моды поляризованной волны. В ряде случаев при проектировании волоконно-оптических датчиков важно учитывать эти параметры и определять их значения [4—10].

В PMF различие показателей преломления из-за асимметричной структуры приводит к различию скоростей распространения света между двумя собственными поляризационными модами. Из-за различия скоростей распространения поляризационные моды имеют разные постоянные распространения, соответствующие „быстрой" (fast) и „медленной" (slow) поляризационным модам [11—13]. При этом величина двулучепреломления определяется как разница постоянных распространения между поляризационными модами:

АР = в, - в f, (1)

где Ps f = — ns f — постоянные распространения для „медленной" и „быстрой" поляризационных мод, X — длина волны излучения, ns,f — показатели преломления волокна вдоль „медленной" и „быстрой" осей. Это выражение можно также записать в виде:

Ав = f^n, (2)

где Дп = ns - nf— разность показателей преломления волокна вдоль „медленной" и „быстрой" осей.

Одним из основных параметров PMF является длина биений Л, определяющая пространственный период, через который две ортогональные поляризационные моды, распространяясь по собственным осям, приходят в фазе [14, 15]:

Л=^=А. (3)

Ав Дп

Таким образом, параметры Дв и Л характеризуют разность фазовых задержек двух поляризационных мод. Временная разность групповых задержек поляризационных мод для длины участка L:

==AnL (4)

г с

называется поляризационной модовой дисперсией (PMD).

В настоящей работе представлена методика определения величины двулучепреломления с применением интерферометра Саньяка [16] и перестраиваемого по длине волны источника излучения.

Описание методики. На рис. 1, a представлена схема эксперимента по определению двулучепреломления оптических волокон с применением интерферометра Саньяка. В данном случае линейно-поляризованное излучение от перестраиваемого лазерного источника (TLS) поступает на волоконный делитель 1x2 (C1), где разделяется по двум направлениям. Часть излучения поступает на вспомогательный интерферометр 2, представляющий собой несимметричный интерферометр Маха—Цандера, образованный двумя волоконными делителями C3 и С4, с известной временной задержкой Td между плечами интерферометра. Другая часть излучения поступает в интерферометр Саньяка, образованный волоконным делителем C2 и волоконным контуром 1 (FUT, Fiber Under Test — исследуемое волокно). Все оптические элементы схемы изготовлены из PMF типа PANDA. Для соединения волоконного контура с C2 использовано изотропное одномодовое волокно для нарушения соосности PMF волокон исследуемого контура и C2, с целью „разбиения" линейной поляризации по двум ортогональным осям, как показано на рис. 1, б и в.

интерферометр Рис. 1

Частота оптического излучения от перестраиваемого лазерного источника изменяется согласно линейному закону [17]:

& = yt +ю0, (5)

где у — скорость перестройки циклической частоты, ю0 — начальная частота оптического излучения. Данное излучение после прохождения С2 разделяется по двум направлениям обхода волоконного контура: по часовой стрелке (С^ и против (СС^. В местах соединения С2 с волоконным контуром линейно-поляризованное излучение „разбивается" на две собственные ортогональные поляризационные моды. Вследствие наличия двулучепреломления эти моды как будто распространяются по двум различным оптическим контурам для каждого из направлений. При повторном прохождении областей соединения волоконного контура и С2 поляризационные моды для каждого направления „смешиваются" между собой, при этом происходит интерференция полей для обоих направлений CW и СС^ Результат интерференции регистрируется фотоприемником РБ1, сигнал с которого записывается при помощи осциллографа.

Интенсивность излучения, приходящего на фотоприемник, можно описать как:

I « 1 + cos (Аф(*)), (6)

где Дф(0 = ф^(0 - фХО — разность фаз между поляризационными модами. Поскольку частота излучения изменяется во времени согласно (5), то разность фаз не является постоянной величиной и зависит от времени. Рассмотрим данное обстоятельство подробнее. Фаза сигнала определяется как интеграл от частоты, поэтому, учитывая (5):

ф(?) =Jra(t)dt = J(yt + ш0)dt = -f yt2 + Q0t + ф0 . (7)

Поляризационная мода, распространяющаяся по „быстрой" оси, опережает моду, распространяющуюся по „медленной" оси, на величину поляризационной модовой дисперсии т,, так что ф/t) = 9s(t - т,). Учитывая это, а также (7), уравнение (6) можно записать как:

1«1 + cos(yx,t -yx,2 -Ш0Т,). (8)

Применяя к (8) преобразование Фурье, можно вычленить частоту периодического сигнала f которая в данном случае определяется как f = ут,/2л. Вычислив данную частоту и учитывая (4), можно найти разность показателей преломления:

An =

2/ /X2

(9)

yL vL

где v = dk/dt — скорость перестройки длины волны лазера. Определив разность показателей преломления по (9) и учитывая (2)—(4), можно найти основные параметры PMF.

Сигнал вспомогательного интерферометра формируется аналогично рассмотренному выше сигналу, только время задержки Td определяется разницей оптических путей между плечами интерферометра. Если известна величина Td, можно, проанализировав сигнал вспомогательного интерферометра, определить реальную скорость перестройки; далее по алгоритму, представленному в [18—20], можно скомпенсировать нелинейность перестройки лазера для повышения точности определения измеряемых величин.

Результаты эксперимента. Для проведения эксперимента согласно рис. 1 была собрана схема. В работе использовался TLS 81606A производства Keysight с параметрами перестройки k = 1550 нм — центральная длина волны, Дк = 4 нм — диапазон перестройки лазера, v = 100 нм/с — скорость перестройки, фотоприемники FPD 610-FC NIR производства MenloSystems, осциллограф Tektronix DPO 7254. В эксперименте использовался волоконный контур длиной L = 500 м из PMF производства ПАО „ПНППК".

На рис. 2 представлен пример нормированного выходного сигнала интерферометра Саньяка.

А,о.е. 0,75 0,50 0,25 0

-0,25 -0,50 -0,75 -1,00

0,0022 0,0024 0,0026 0,0028 t, с

Рис. 2

В ходе эксперимента снят ряд интерференционных сигналов. Для их обработки написана программа на языке Python. Результат обработки в виде частотного спектра представлен на рис. 3.

А, дБ 0

-10 -2-30 -40 -50

0 10 20 30 40 50 60 f кГц

Рис. 3

В данном эксперименте спектры на рис. 3 имеют пик на частоте f = 15,342 кГц. Полученное значение частоты позволяет рассчитать основные параметры PMF согласно (2)—(4), (9). Значения полученных параметров приведены ниже:

— разность показателей преломления An = 7,4 10-4;

3 -1

— двулучепреломление Ар = 3 10 м ;

— длина биений Л = 2 мм;

— поляризационная модовая дисперсия тр = 1,23 нс;

— удельная поляризационная модовая дисперсия тр /L = 2,46 нс/км.

Измеренная разница показателей преломления An = 7,4 10-4 попадает в диапазон паспортных значений An = (7,1—7,7)10-4, что подтверждает работоспособность представленной методики.

Заключение. Представлена методика определения величины двулучепреломления с применением перестраиваемого лазерного источника и интерферометра Саньяка. По результатам эксперимента данная методика показала высокую точность определения параметров волокна типа Panda, в том числе за счет использования вспомогательного интерферометра. Измеренная с применением предложенной методики разность показателей преломления вдоль „быстрой" и „медленной" осей составила 7,4 10-4, что попадает в диапазон паспортных значений.

СПИСОК ЛИТЕРАТУРЫ

1. Fujii Y. and Sanos K. Polarization transmission characteristics of optical fibers with elliptical cross section // Electron. Comm. Jpn. 1980. Vol. 63. Pt. I. P. 87—93. https://doi.org/10.1002/ecja.4400630812.

2. Kumar A., Varshney R. K., Thyagarajan K. Birefringence calculations in elliptical-core optical fibres // Electronics Letters. 1984. Vol. 20. Р. 112.

3. Noda J., Okamoto K. and Sasaki Y. Polarization-maintaining fibers and their applications // Journal of Lightwave Technology. 1986. Vol. 4, N 8. P. 1071—1089. DOI: 10.1109/JLT.1986.1074847.

4. Medeiros A., Barcelos S., Rigon E., Rando R., Seminario J., Santos M., and Oliveira R. PMD Characterization of Installed Fiber Networks - Compromise Between Result Accuracy and Measurement Time // Optical Fiber Communication Conference and Exposition and The National Fiber Optic Engineers Conference, Technical Digest (CD). Optica Publishing Group, 2006. Art. no. NTuC4.

5. Ten S., Edwards M. An Introduction to the Fundamentals of PMD in Fibers // White Paper. 2006. Vol. WP5051.

6. Gilev D. G., Valiushina P. M., Maksimenko V. A., and Krishtop V. V. Determination of induced birefringence in fiber-optic resonator from frequency difference between main and additional resonance peaks // Opt. Continuum. 2022. Vol. 1, N 3. Р. 487—493.

7. Drozdov I. R., Ovchinnikov K. A., Boychuk E. S., & Krishtop V. V. Two-stage method for comparing the lengths of optical fibers using OFDR // Intern. Conf. Laser Optics (ICLO). June 2022.

8. Гилев Д. Г. Оценка динамического диапазона датчика угловой скорости на основе волоконно-оптического резонатора // Прикладная фотоника. 2023. Т. 10, № 1. С. 116—130.

9. Денисов И. В., Лисовский Н. В. Систематизация чувствительных к изгибу волоконных световодов // Прикладная фотоника. 2023. Т. 10, № 1. С. 131—148.

10. Краузин П. В., Трефилов И. А., Сагирова А. Р., Оглезнев А. А., Кондрашов А. Н., Юрина А. Д., Мельников Р. М. Алгоритм калибровки распределенного датчика температуры // Прикладная фотоника. 2023. Т. 10, № 2. С. 68—87.

11. Kumar A., Ghatak A. Polarization of light with applications in optical fibers. SPIE, 2011. 246 р. ISBN: 9780819482150. https://doi.org/10.1117/3.861761.

12. Rogers A. Polarization in Optical Fibers. Artech house, 2008. 250 р. ISBN: 9781580535342.

13. Lefevre H. Fiber-optic gyroscopes. Artech house, 2022. 416 р.

14. Бурдин В. А. Измерения длины биений оптических волокон при зондировании импульсами увеличенной длительности // T-Comm-Телекоммуникации и Транспорт. 2009. № S1. С. 175—177.

15. Моршнев С. К. и др. Измерение длины биений в двулучепреломляющих волоконных световодах // Фотоника. 2018. Т. 12, №. 6. С. 616—633.

16. Малыкин Г. Б. Эффект Саньяка. Корректные и некорректные объяснения // Успехи физических наук. 2000. Т. 170, №. 12. С. 1325—1349.

17. Song J. Optical frequency domain reflectometry: Sensing range extension and enhanced temperature sensitivity: M.Sc. Thesis. University of Ottawa, 2014. 113 р.

18. Song Jia, Li Wenhai, Lu Ping, Xu Yan,ing, Chen Liang & Bao Xiaoyi. Long-Range High Spatial Resolution Distributed Temperature and Strain Sensing Based on Optical Frequency-Domain Reflectometry // Photonics Journal. 2014. Vol. 6. Р. 1—8. doi 10.1109/JPHOT.2014.2320742.

19. Gilev D. G., Ovchinnikov K. A., Krishto, V. V. et al. Fiber Optic Resonators for Angular Rate Sensors // Bull. Russ. Acad. Sci. Phys. 2022. Vol. 86, Suppl. 1. Р. S75—S80. https://doi.org/10.3103/S1062873822700423.0vchinnikov K. A., Gilev D. G., Krishto, V.V. et al. Application of Optical Frequency Domain Reflectometry for the Study of Polarization Maintaining Fibers // Bull. Russ. Acad. Sci. Phys. 2022. Vol. 86, Suppl. 1. Р. S156—S162. https://doi.org/10.3103/S1062873822700599.

20. Ovchinnikov K. A., Gilev D. G., Krishto, V. V. et al. Application of Optical Frequency Domain Reflectometry for the Study of Polarization Maintaining Fibers // Bull. Russ. Acad. Sci. Phys. 2022. Vol. 86, Suppl. 1. Р. S156—S162. https://doi.org/10.3103/S1062873822700599.

Сведения об авторах

Константин Александрович Овчинников — аспирант; Пермский национальный исследовательский

политехнический университет; Пермская научно-производственная приборостроительная компания, начальник лаборатории; E-mail: [email protected]

Виктор Владимирович Криштоп — д-р физ.-мат. наук, профессор; Пермский государственный

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

национальный исследовательский университет; Пермский национальный исследовательский политехнический университет; Пермская научно-производственная приборостроительная компания, главный научный сотрудник; E-mail: [email protected]

Даниил Георгиевич Гилев — канд. техн. наук; Пермская научно-производственная при-

боростроительная компания, начальник бюро; Пермский национальный исследовательский политехнический университет, научный сотрудник; E-mail: [email protected]

Виталий Александрович Максименко — канд. физ.-мат. наук, доцент; Пермский национальный

исследовательский политехнический университет, кафедра общей физики; доцент; E-mail: [email protected]

Анатолий Викторович Перминов — д-р физ.-мат. наук, доцент; Пермский национальный ис-

следовательский политехнический университет, кафедра общей физики; заведующий кафедрой; E-mail: [email protected]

Поступила в редакцию 07.06.23; одобрена после рецензирования 19.06.23; принята к публикации 28.08.23.

REFERENCES

1. Fujii Y. and Sanos K. Electron. Comm. Jpn., pt. I, 1980, vol. 63, pp. 87-93, https://doi.org/10.1002/ ecja.4400630812.

2. Kumar A., Varshney R.K., Thyagarajan K. Electronics Letters, 1984, vol. 20, pp. 112.

3. Noda J., Okamoto K. and Sasaki Y. Journal of Lightwave Technology, 1986, no. 8(4), pp. 1071-1089, DOI: 10.1109/JLT.1986.1074847.

4. Medeiros A., Barcelos S., Rigon E., Rando R., Seminario J., Santos M., and Oliveira R. Optical Fiber Communication Conference and Exposition and The National Fiber Optic Engineers Conference, Optica Publishing Group, 2006, paper NTuC4.

5. Ten S., Edwards M. White Paper, 2006, vol. WP5051.

6. Gilev D.G., Valiushina P.M., Maksimenko V.A., and Krishtop V.V. Opt. Continuum, 2022, no. 1(3), pp. 487-493.

7. Drozdov I.R., Ovchinnikov K.A., Boychuk E.S., & Krishtop V.V. International Conference Laser Optics (ICLO), 2022, pp. 1-1.

8. Gilev D.G. Applied photonics, 2023, no. 1(10), pp. 116-130. (in Russ.)

9. Denisov I.V., Lisovsky N.V. Applied photonics, 2023, no. 1(10), pp. 131-148. (in Russ.)

10. Krauzin P.V., Trefilov I.A., Sagirova A.R., Ogleznev A.A., Kondrashov A.N., Yurina A.D., Melnikov R.M. Applied photonics, 2023, no. 2(10), pp. 68-87. (in Russ.)

11. Kumar A., Ghatak A. Polarization of light with applications in optical fibers, SPIE, 2011, 246 p., ISBN: 9780819482150, https://doi.org/10.1117/3.861761.

12. Rogers A. Polarization in Optical Fibers, Artech house, 2008, 250 p., ISBN: 9781580535342.

13. Lefevre H. Fiber-optic gyroscopes, Artech house, 2022, 416 p.

14. Burdin V.A. T-Comm: Telecommunications and transport, 2009, no. S1, pp. 175-177. (in Russ.)

15. Morshnev S.K., Gubin V.P., Starostin N.I., Prshiyalkovsky Ya.V., Sazonov A.I. Photonics Russia, 2018, no. 6(12), pp. 616-633. (in Russ.)

16. Malykin G.B. Physics-Uspekhi (Advances in Physical Sciences), 2000, no. 12(43), pp. 1229-1252.

17. Song J. Optical frequency domain reflectometry: Sensing range extension and enhanced temperature sensitivity: M.Sc. Thesis, Université d'Ottawa, 2014, 113 p.

18. Song Jia, Li Wenhai, Lu Ping, Xu Yanping, Chen Liang, & Bao Xiaoyi. Photonics Journal, 2014, vol. 6, pp. 1-8, DOI: 10.1109/JPHOT.2014.2320742.

19. Gilev D.G., Ovchinnikov K.A., Krishtop V.V. et al. Bull. Russ. Acad. Sci. Phys., 2022, vol. 86, suppl. 1, pp. S75-S80, https://doi.org/10.3103/S1062873822700423.

20. Ovchinnikov K.A., Gilev D.G., Krishtop V.V. et al. Bull. Russ. Acad. Sci. Phys., 2022, vol. 86, suppl. 1, pp. S156-S162, https://doi.org/10.3103/S1062873822700599.

Konstantin A. Ovchinnikov

Victor V. Krishtop

Daniil G. Gilev

Vitaliy A. Maksimenko

Anatoliy V. Perminov

Data on authors

Post-Graduate Student; Perm National Research Polytechnic University; Perm Research and Production Instrument-Making Company, Head of Laboratory; E-mail: [email protected]

Dr. Sci., Professor; Perm State University; Perm National Research Polytechnic University; Perm Research and Production Instrument-Making Company, Chief Researcher; E-mail: [email protected]

PhD; Perm Research and Production Instrument-Making Company; Perm National Research Polytechnic University, Researcher; E-mail: [email protected]

PhD, Associate Professor; Perm National Research Polytechnic University, Department of General Physics, Associate Professor; E-mail: [email protected]

Dr. Sci., Associate Professor; Perm National Research Polytechnic University, Department of General Physics, Head of the Department; E-mail: [email protected]

Received 07.06.23; approved after reviewing 19.06.23; accepted for publication 28.08.23.

i Надоели баннеры? Вы всегда можете отключить рекламу.