Разработка метода лазерной спектроскопии и аппаратно-программного комплекса для экологического мониторинга подводных акваторий c использованием телеуправляемых необитаемых подводных аппаратов

Полный текст

Введение Миссия экологического мониторинга морских акваторий, в случае загрязнения нефтепродуктами, включает в себя не только исследование непосредственно самого процесса загрязнения нефтепродуктами, но и оценку результатов воздействия процессов загрязнения на фитопланктонное сообщество. Мониторинг непосредственно самого процесса загрязнения нефтепродуктами состоит из последовательного решения несколько задач: детектирование загрязнений, исследование динамики загрязнения, определение объемов загрязнения, определение сортов нефтепродуктов на загрязненной акватории и выявление источников загрязнения. В работе рассматривается загрязнение подводных морских акваторий, когда нефтепродукты находятся в морской воде в растворенном состоянии. В этом случае задача определения объемов загрязнения сводится к измерению концентрации растворенных нефтепродуктов. Такие ситуации возникают в случае сброса льяльных вод судами, в случае аварий на подводных нефтепроводах, при наличии на морской поверхности пленок нефтепродуктов, также происходит их растворение и поступление в более глубокие слои. Процесс определения ущерба, который наносится морской экосистеме, является сложным, включает в себя много направлений и, в каждом конкретном случае, эти задачи решаются по уже отработанным методикам [1-2][7]. Летальное воздействие загрязнений на подводную морскую флору и фауну можно оценить путем подсчета погибшей рыбы, птиц, моллюсков, подводных растений и т.д. Воздействие нефтяных загрязнения на фитопланктонные сообщества является предметом исследований. В работах [3; 4] представлены обзоры, наиболее полно описывающие результаты в этой области. В случае загрязнения больших масштабов используются спутниковые данные, которые позволяют при усреднении по большим площадям морской поверхности оценить изменения, которые происходят с фитопланктонными сообществами в самом приповерхностном слое океана[8]. Однако в настоящее время отсутствуют оперативные методы, которые позволяют делать оценки воздействия загрязнений на фитопланктонные сообщества insitu под водой в локальных масштабах. В работе описывается новый подход к разработке метода и технических средств, позволяющих решать задачи, связанные как с мониторингом непосредственно самого процесса загрязнения подводных акваторий нефтепродуктами, так и с оценкой изменения состояния фитопланктонных сообществ под воздействием этих загрязнений. Подход основан на использовании метода лазерной индуцированной флуоресценции (ЛИФ). Он развивается и используется для измерения концентрации растворенных в морской воде углеводородов нефти, концентрации хлорофилла А в морской воде и удельного воспроизводства растворенного органического вещества клетками фитопланктона [5-8]. Две последние характеристики могут использоваться для описания состояния фитопланктонных сообществ. Спектр ЛИФ содержит информацию о концентрации растворенного органического вещества (РОВ) (включая растворенные нефтепродукты) и концентрации хлорофилла. Проведенные ранее исследования показали, что измерение концентрации растворенных в морской воде нефтепродуктов возможно при использовании лазерного излучения с длиной волны в УФ-области спектра (менее 300 нм), так как наиболее интенсивная область ЛИФ от растворов нефтепродуктов находится в спектральном диапазоне от 300 до 450 нм [9; 10]. Вместе с тем в [7; 8; 11] было показано, что наиболее подходящий диапазон лазерного излучения для измерения концентрации хлорофилла А и РОВ, воспроизводимого клетками фитопланктона, лежит в «зеленой» области спектра. Кроме того, при облучении УФ-излучением часть клеток фитопланктона может погибать. Таким образом, для того чтобы осуществлять одновременный мониторинг состояния фитопланктонных сообществ и загрязнения морских вод в виде растворенных нефтепродуктов, необходимо использовать возбуждение спектров ЛИФ двумя длинами волн, одна из которых лежит в УФ-области, а другая - в зеленой области длин волн. Ранее были разработаны спектрометры, которые позволяли решать по отдельности задачи мониторинга состояния фитопланктонных сообществ и процессы загрязнения акваторий нефтепродуктами. В [10; 11] проведена разработка метода и спектрометра для измерения концентрации растворенных нефтепродуктов в морской воде, определения сортов нефтепродуктов с использованием УФ-излучения[9]. В [12-13] приведены результаты разработки погружаемого лазерного спектрометра, предназначенного для использования на телеуправляемом необитаемом подводном аппарате (ТНПА) для исследования распределения хлорофилла А и процессов воспроизводства РОВ клетками фитопланктона. В этом спектрометре использовалась вторая гармоника Nd:YAG лазера с длиной волны излучения 532 нм. В настоящей работе описываются результаты разработки малогабаритного лазерного спектрометра для ТНПА, который предназначен одновременно как для мониторинга загрязнения подводных акваторий растворами нефтепродуктов, так и для исследования состояния фитопланктонных сообществ на акваториях, подверженных загрязнению. Описание методики мониторинга Методика мониторинга состояния фитопланктонных сообществ основана на исследовании спектров ЛИФ от морской воды. На рис. 1 приведены спектры ЛИФ от морской воды, в которой содержатся клетки фитопланктона и не содержатся растворы нефтепродуктов. Рис.1. Спектр ЛИФ морской воды при возбуждении лазерным излучением с длиной волны 532 нм / Figure 1. LIF spectrum of sea water upon excitation by laser radiation with a wavelength of 532 nm По вертикальной оси - интенсивность спектра ЛИФ, по горизонтальной оси - длины волн в нанометрах. Максимум с центром на длине волны 648 нм соответствует комбинационному рассеянию (КР) морской воды, линия флуоресценции хлорофилла А наблюдается на длине волны 680 нм. Спектральное распределение, представленное на рис. 1 сплошной черной кривой, соответствует спектру ЛИФ от РОВ, которое воспроизведено клетками фитопланктона. Спектр ЛИФ был получен в лабораторных условиях на пробах морской воды из Амурского залива. Для измерения концентрации хлорофилла А определяется интенсивность линии флуоресценции на длине волны 680 нм (за вычетом величины интенсивности РОВ на этой длине волны). Затем проводится нормировка интенсивности линии хлорофилла А на интенсивность линии КР. Концентрация хлорофилла А (CchA) определяется по нормированной величине интенсивности линии флуоресценции хлорофилла А: СchlA = k (IchlA/Iкр), (1) где k - калибровочная константа, которая определяется в процессе калибровки метода на стандартный метод измерения [14; 15]. На рис. 2 приведен один из примеров калибровки метода ЛИФ на стандартный фотометрический метод измерения концентрации хлорофилла А, которая проводилась на пробах морской воды с различной концентрацией хлорофилла А в лабораторных условиях. Измерения выполнены на пробах морской воды, которые были отобраны в период цветения водорослей в различных точках Амурского залива. По вертикальной оси - значение концентрации хлорофилла А, которое получено в пробах воды стандартным фотометрическим методом. По горизонтальной оси - значение нормированной интенсивности ЛИФ. Подобные калибровки проточных и погружаемых спектрометров ЛИФ использовались для измерения концентрации хлорофилла А на спектрометрах ЛИФ в натурных условиях [16; 17]. В верхней части приведено уравнение линейной регрессии и коэффициент детерминации R2. Состояние фотосинтезирующего аппарата клеток фитопланктона можно анализировать не только по общему количеству хлорофилла А в морской воде, но и по темпам воспроизводства РОВ клетками фитопланктона [7; 8]. Чтобы оценить скорость воспроизводства РОВ клетками фитопланктона, необходимо проводить измерение не общего объема РОВ, содержащегося на момент измерения в морской воде, а только «молодого» РОВ (или «лабильной части» РОВ), т. е. той части РОВ, которая воспроизводится живыми клетками фитопланктона за период активной жизнедеятельности клеток (например, за период цветения водорослей) и не успела пройти процессы деградации. В периоды цветения водорослей такой подход к оценке состояния клеток наиболее нагляден. Объем РОВ, воспроизводимого клетками фитопланктона, описывается параметром Q, который определяется как площадь под кривой спектрального распределения ЛИФ от РОВ, в диапазоне от 560 до 740 нм (заштрихованная часть спектра на рис. 1). Q (2) Параметр Q пропорционален количеству РОВ и при соответствующей калибровке на стандартные методы измерения можно проводить измерения концентрации РОВ в абсолютных единицах. Рис. 2. Пример калибровочной зависимости для измерения концентрации хлорофилла А по интенсивности спектра ЛИФ / Figure 2. An example of a calibration dependence for measuring the concentration of chlorophyll A by the intensity of the LIF spectrum Измеряя в различных местах выбранной водной акватории (или на различных глубинах) значения концентрации хлорофилла А (СchlA) и значения параметра Q, можно восстанавливать диаграммы рассеяния Q - C, которые показывают функциональную связь между этими величинами, т.е. можно восстановить соотношения между конкретной концентрацией хлорофилла А и тем значением величины РОВ, которое клетки воспроизводят в результате своей жизнедеятельности. На рис. 3 приведены виды функциональных Q - C зависимостей в различных биопродуктивных районах Мирового океана в период цветения водорослей на рассматриваемых акваториях. Измерения проведены в периоды цветения водорослей, т.е. в периоды, когда клетки фитопланктона имеют максимальные концентрации хлорофилла А и воспроизводят максимальное количество РОВ. Диаграммы рассеяния Q-C хорошо аппроксимируются линейными функциями, что дает возможность определить темпы воспроизводства РОВ (ν = dQ/dC) в процессе цветения водорослей [7; 8; 18]. Для анализа удельного воспроизводства РОВ клетками фитопланктона (величины ν) достаточно определять величину РОВ в относительных единицах (Q). Сравнение величины ν в различные периоды жизнедеятельности клеток фитопланктона на данной акватории позволяет определять периоды цветения водорослей, стадии активного воспроизводства РОВ или состояния, в котором процессы жизнедеятельности клеток подавляются наличием растворенных нефтепродуктов. Рис. 3. Примеры Q - C диаграмм для различных биопродуктивных районов Мирового океана в периоды цветения водорослей или близкие к ним: 1 - Охотское море; 2 - Северное море; 3 - Юго-западная часть Атлантического океана / Figure 3. Examples of Q - C diagrams for various bio productive regions of the world's oceans during or close to algal blooms: 1 - Sea of Okhotsk; 2 - North Sea; 3 - Southwestern part of the Atlantic Ocean Таким образом, для оценки состояния фитопланктонного сообщества на данной акватории, в рамках данного метода, используются две характеристики: концентрация хлорофилла А - СchlA и величина ν, которая характеризует удельное воспроизводства РОВ на анализируемой акватории. При наличии растворенных углеводородов нефтепродуктов величины СchlA и ν будут меняться в результате воздействия загрязнения на аппарат фотосинтеза клеток фитопланктона. Задача оценки изменения состояния фитопланктонных сообществ сводится в этом случае к сравнению средних концентраций и темпов воспроизводства РОВ, которые имели место непосредственно перед загрязнением на данных акваториях, и тех значений, которые определяются в период загрязнения. В процессе мониторинга необходимо проводить измерение концентраций растворенных нефтепродуктов на исследуемой акватории. В [8; 9] описаны методика определения концентраций растворенных нефтепродуктов, определение сорта нефтепродуктов с использованием методики ЛИФ. На рис. 4 продемонстрированы спектры ЛИФ растворов судового топлива DMA в морской воде с различными концентрациями. Для возбуждения ЛИФ использовалось излучение лазера с длиной волны 278 нм. По вертикальной оси - интенсивность спектра ЛИФ в относительных единицах, по горизонтальной оси - длины волн в нанометрах. Справа приведены значения концентраций растворенного DMA и их соответствие цветам на графиках. На рис. 5 приведена калибровка метода измерения концентрации раствора DMA по спектрам ЛИФ на стандартный метод измерения концентрации растворенных нефтепродуктов в морской воде с использованием ИК Фурье-спектроскопии. По вертикальной оси отложено значение интегрального параметра Q (см. формулу 2), где интегрирование производится в диапазоне длин волн от 300 до 525 нм, по горизонтальной оси отложена концентрация DMA, определенная стандартным методом ИК Фурье-спектроскопии. Аналогичные калибровочные линейные регрессии были построены и для других видов судового топлива, включительно и сырой нефти. Для каждого вида топлива спектры ЛИФ имеют характерный вид, что позволяет идентифицировать типы нефтепродуктов. Создание библиотеки спектров ЛИФ и разработка процедуры компьютерного обучения позволили создать элементы искусственного интеллекта, которые в автономном режиме способны определять концентрацию растворенных нефтепродуктов и их типы [10]. Рис. 4. Спектры ЛИФ растворов судового топлива DMA в морской воде / Figure 4. LIF spectra of DMA marine fuel solutions in sea water Рис. 5. Калибровочная зависимость при определении концентрации растворенного DMA в морской воде / Figure 5. Calibration dependence for determining the concentration of dissolved DMA in sea water Разработка двухчастотного спектрометра ЛИФ для ТНПА На рис. 6а, б приведена функциональная схема малогабаритного спектрометра ЛИФ, разработанного для ТНПА. а б Рис. 6. Схема малогабаритного, двухчастотного спектрометра ЛИФ (а) и фото спектрометра в герметичном корпусе (б): 1 - спектрометр Ocean Optics Maya 2000Pro; 2 - процессорная плата ITX-N29; 3 - светодиод LED - LEUVA66H70HF00 (LG); 4 - светодиод LED KSK-304; 5 - объектив 74-ud (Ocean Optics); 6, 7 - фильтры Fgl550 Thor labs; 8 - кварцевые иллюминаторы; 9 - герметичный корпус / Figure 6. Schematic diagram of a small-sized, dual-frequency spectrometer LIF (a) and photo spectrometer in a sealed case (b): 1 - Ocean Optics Maya 2000Pro spectrometer; 2 - ITXN29 processor board; 3 - LED - LEUVA66H70HF00 (LG); 4 - LED KSK304; 5 - lens 74ud (Ocean Optics); 6, 7 - Thor labs Fgl550 filters; 8 - quartz portholes; 9 - sealed housing При использовании спектрометра ЛИФ на ТНПА существенным являются массогабаритные параметры аппаратуры и энергопотребление. Для их уменьшения, в качестве излучателей для возбуждения ЛИФ, использовались светодиоды (LED). Для мониторинга растворов нефтепродуктов - LED типа LEUVA66H70HF00 c длиной волны излучения 278 нм (3 на рис. 6). Для возбуждения спектров ЛИФ от клеток фитопланктона и «молодого» РОВ использовался LED типа KSK- 340 (4 на рис. 6) с длиной волны излучения 532 нм. Излучение светодиодов (3, 4) через кварцевые иллюминаторы (8), направлялись в воду, флюоресценция регистрировалась объективом (5). Фильтры 6, 7 использовались для подавления упругого рассеяния в каждом случае. При измерении спектра ЛИФ, возбуждаемого длиной волны 278 нм, излучал только светодиод (3), перед объективом 5, располагался фильтр для подавления длины волны упругого рассеяния на 278 нм. В случае измерения спектра на 532 нм в канал регистрации вводился фильтр, подавляющий эту длину волны. Фильтры размещались перед объективом автоматически. Излучение, собранное объективом 5, направлялось на входную цель спектрометра 1. Управление работой всего спектрометра и обработка спектров осуществлялось процессором 2. Спектрометр размещался в герметичном корпусе 9. Подключение спектрометра осуществляется к бортовому питанию и линии связи ТНПА. Передача данных и управление осуществляется при помощи Ethernet-канала ТНПА. Подводное соединение спектрометра выполнено в виде герморазъема (10), позволяющего в любой момент отсоединить спектрометр от ТНПА и демонтировать устройство. Лабораторные эксперименты Разработка двухчастотного спектрометра ЛИФ позволила проводить эксперименты по детектированию растворенных нефтепродуктов в морской воде, определять их тип и концентрацию и одновременно исследовать их воздействие на состояние фотосинтезирующего аппарата клеток фитопланктона. На рис. 7 и 8 продемонстрированы результаты измерения спектров ЛИФ в пробах чистой морской воды и в пробах, где присутствовал раствор судового топлива DMS. Данные спектры были получены в лабораторных условиях. Кювета, содержащая пробы морской воды, располагалась непосредственно перед выходными иллюминаторами спектрометра (иллюминаторы 8 на рис. 6а). После измерения спектров в кювету добавлялся раствор судового топлива MGO в таком количестве, чтобы итоговая концентрация топлива в кювете составляла порядка 50 ppm. В зеленом цвете на обоих графиках представлены спектры ЛИФ при возбуждении излучением в зеленой области спектра (длина волны 532 нм), в фиолетовом цвете - спектры, возбуждаемые в ультрафиолете (278 нм). Для наглядности сравнения величины сигналов, вертикальные шкалы слева показывают величины сигналов ЛИФ для УФ-излучения, шкала справа демонстрирует величину для возбуждения 532 нм. Сигнал флуоресценции хлорофилла А, при возбуждении длиной волны 523 нм, значительно превышал флуоресценцию при возбуждении длиной волны 278 нм, при одинаковой концентрации хлорофилла А в пробе. Поэтому показания относительной интенсивности по правой шкале приведены к значениям интенсивности на левой шкале, но масштаб шкалы необходимо умножить на 10. Для случая чистой морской воды (рис. 7) сравнение интенсивностей спектра ЛИФ на центральной длине волны линии флуоресценции хлорофилла А (680 нм) показывает, что флуоресценция наиболее эффективно возбуждается излучением в зеленой области спектра, т. е. длиной волны 532 нм. Это относится также и к интенсивности флуоресценции лабильной части РОВ. Наблюдается значительное превышение интенсивности широкополосной флуоресценции, возбуждаемой излучением 532 нм, в интервале от 620 до 660 нм над интенсивностью этого интервала при возбуждении излучением 278 нм. В то время как в коротковолновой области спектра, где флуоресцирует РОВ, прошедшее стадию разложения, в диапазоне длин волн от 350 до 500 нм наблюдается только интенсивность флуоресценции, которая возбуждается УФ-излучением. На рис. 8 приведены спектры ЛИФ в случае, когда в пробу морской воды добавлялся раствор судового топлива MGO. Концентрация раствора в кювете равнялась 50 ppm. В этом случае наблюдается значительное снижение линии флуоресценции хлорофилла А. Уменьшение интенсивности линии флуоресценции можно оценить по графику на рис. 9. рис 6 Рис. 7. Спектры ЛИФ, возбуждаемые УФ и зеленым излучением в чистой морской воде с фитопланктоном / Figure 7. LIF spectra excited by UV and green radiation in pure sea water with phytoplankton Рис. 8. Спектры ЛИФ возбуждаемые, УФ и зеленым излучением в морской воде, содержащей раствор судового топлива / Figure 8. LIF spectra excited by UV and green radiation in sea water containing a marine fuel solution В этом эксперименте раствор добавлялся в кювету сразу после измерения спектра ЛИФ от чистой воды. Затем, после 10 минут, было проведено измерение спектра ЛИФ. Интенсивность линии флуоресценции уменьшилась более чем на порядок. Это может быть вызвано двумя причинами. Первая - это гибель клеток фитопланктона и разрушения молекул хлорофилльных комплексов. Вторая - это уменьшение квантового выхода флуоресценции за счет воздействия на хлорофилл загрязнения нефтепродуктами. В любом случае, наблюдается прямое воздействие загрязнение на клетки фитопланктона. Нам представляется более вероятным первый механизм воздействия, который приводит к прямой гибели клеток. Об этом свидетельствует и тот факт, что в области флуоресценции «молодого» РОВ наблюдается увеличение интенсивности спектральных линий. То есть происходит увеличение РОВ за счет гибели клеток фитопланктона. Подобное поведение интенсивностей линий хлорофилла А и «молодого» РОВ наблюдалось в работе [19], где гибель клеток хлорофилла А вызывалась нагреванием проб морской воды. Рис. 9. Сравнение интенсивности флуоресценции хлорофилла А в случае чистой морской воды и при загрязнении раствором судового топлива MGO / Figure 9. Comparison of the intensity of chlorophyll A fluorescence in the case of pure sea water and when contaminated with MGO marine fuel solution Заключение Результаты разработки метода ЛИФ и малогабаритного аппаратно-программного комплекса лазерного спектрометра для ТНПА создают методическую и аппаратурную базу для оперативного мониторинга загрязнения подводных акваторий растворенными нефтепродуктами в локальных масштабах. Использование двухчастотного спектрометра ЛИФ на ТНПА позволяет проводить insitu измерение параметров, характеризующих объемы загрязнений их тип, исследовать реакцию фитопланктонного сообщества на эти загрязнения. Использование спектроскопии ЛИФ на ТНПА позволяет исследовать реакцию различных видов фитопланктона не только на наличие растворенных нефтепродуктов, при различных параметрах морской воды, а также и на диспергенты, которые используются для ликвидации нефтяных пленок.

Об авторах

Дмитрий Юрьевич Прощенко

Морской государственный университет имени адмирала Г.И. Невельского

Email: dima.prsk@mail.ru
ORCID iD: 0000-0001-8378-8208

кандидат физико-математических наук, заведующий лабораторией морской робототехники научно-инновационного комплекса «Морской технопарк»

Российская Федерация, 690003, г. Владивосток, ул. Верхнепортовая, д. 50а

Илья Олегович Букин

Морской государственный университет имени адмирала Г.И. Невельского

Автор, ответственный за переписку.
Email: il_bukin@mail.ru
ORCID iD: 0000-0002-2929-890X

младший научный сотрудник лаборатории морской робототехники научно-инновационного комплекса «Морской технопарк»

Российская Федерация, 690003, г. Владивосток, ул. Верхнепортовая, д. 50а

Список литературы