1. Введение
Нефтяные загрязнения представляют серьёзную угрозу для морских экосистем, качества водной среды и здоровья человека. Также наличие нефтяной пленки может оказать местное влияние на газо-, тепло- и влагообмен системы океан-атмосфера, что необходимо учитывать при численном моделировании атмосферных процессов. При прогнозировании таких явлений погоды, как осадки, туман, низкой слоистой облачности, необходимо учитывать, что частицы нефти, попадая в атмосферу могут представлять собой дополнительные ядра конденсации. Актуальность исследования обусловлена высокой уязвимостью прибрежных и мелководных зон, где даже локальные утечки нефти могут привести к долговременным негативным последствиям. В таких районах необходимы непрерывное наблюдение состояния окружающей среды и разработка доступных, надёжных методов обнаружения нефтяной плёнки
[1]
Цель настоящей работы заключается в разработке и апробации спектрального метода количественной оценки нефтяных загрязнений морской поверхности, основанного на адаптированной математической модели радиолокационного анализа, ранее применявшейся авторами в иной постановке задачи. Проведённые численные эксперименты показали согласованность предложенного подхода с результатами наблюдений
[2][3]
Метод базируется на ключевом физическом признаке нефтяных плёнок — подавлении капиллярно-гравитационных волн, что приводит к изменению спектральных характеристик отражённого радиосигнала. Эффективность такого подхода подтверждается результатами современных исследований с применением радиолокационного (
[4][5][6]
Научная новизна работы состоит в применении спектрального анализа для количественного определения характеристик нефтяных плёнок, теоретическая значимость — в адаптации модели к задачам дистанционного мониторинга загрязнённой поверхности, практическая — в возможности реализации метода на стационарных прибрежных радиолокационных станциях как экономичного и локально привязанного инструмента.
2. Методы и принципы исследования
При попадании нефтепродуктов на морскую поверхность формируется тонкая плёнка, которая подавляет капиллярно-гравитационные волны — основной источник рассеяния в радиолокационном диапазоне
[7]
В рамках численного моделирования были заданы параметры, характерные для прибрежной радиолокационной системы, размещённой на стационарной вышке. Высота антенны составляла 20 метров, углы облучения — 30–45°, что обеспечивает охват участка поверхности длиной около 10 метров. Рабочая частота радара — 10 ГГц (X-band). Схема размещения и область облучения представлены на рис. 1.
Схема радиолокационного облучения с вышки
Для моделирования воздействия загрязнения использовался коэффициент сглаживания, отражающий степень подавления микроволнения плёнкой. Диапазон коэффициента сглаживания k от 1,00 до 0,70 соответствовал различным уровням загрязнения, параметры заданы по [8].
Для количественного анализа спектра отражённого сигнала была использована модель, разработанная авторами, в которой морская поверхность представляется как совокупность независимых флуктуирующих рассеивателей, равномерно распределённых вдоль направления излучения. Шаг между ними определяется из длины облучаемого участка L и числа рассеивателей. Расстояние определяется углом падения излучения. Исходный сигнал моделируется в виде гармонической волны, отражённой от всей поверхности. Суммарный сигнал формируется как суперпозиция откликов от всех рассеивателей с учётом фазовых сдвигов. Его аналитическое выражение имеет следующий вид:
где
[LATEX_FORMULA]r_{n}(t)[/LATEX_FORMULA] флуктуационное отклонение n-го рассеивателя.
Для перехода к спектру мощности используется автокорреляционная функция и теорема Винера–Хинчина. Влияние учитывается коэффициентом сглаживания Значения k и соответствующие им характеристики плёнки были заданы по данным
[8]
3. Основные результаты
Для первичной оценки влияния загрязнения на отражённый радиосигнал была выполнена серия численных экспериментов, сравнивающих спектр при наличии нефтяной плёнки и в условиях чистой воды. На рис. 2 представлены спектры в двух случаях: эталон (без загрязнения) и загрязнённая поверхность при коэффициенте сглаживания k = 0,95, соответствующем умеренному загрязнению.
Подавление высокочастотных компонент спектра при наличии нефтяной плёнки
Для более детального анализа проведено параметрическое исследование зависимости спектральной амплитуды от степени сглаживания поверхности. Это позволяет установить, как различные уровни загрязнения влияют на отражённый сигнал в спектральной области.
На рис. 3 представлена нормализованная амплитуда спектра в зависимости от коэффициента сглаживания
Зависимость амплитуды спектра от коэффициента сглаживания
<0,75[9]
Для получения количественного соответствия между уровнем сглаживания и физическими характеристиками загрязнения была построена табл. 1. Она объединяет результаты численного моделирования спектральной амплитуды с данными из работы
[8]
Соответствие толщины плёнки и коэффициента подавления ряби
| Коэффициент подавления ряби | Примерная толщина нефтяной плёнки, мкм | Характеристика плёнки | Нормализованная амплитуда спектра (по данным графика) |
| 1,00 | 0 | Чистая вода, без загрязнения | 1,00 |
| 0,97 | 0,5–1 | Очень тонкая плёнка, минимальное воздействие | ≈ |
| 0,95 | 1–2 | Начало подавления капиллярных волн | ≈ |
| 0,90 | 2–3 | Умеренное подавление | ≈ |
| 0,85 | 3–4 | Существенное снижение микроволнения | ≈ |
| 0,80 | 4–5 | Почти полное подавление капиллярных волн | ≈ |
| 0,78 | 5–6 | Плотная нефтяная плёнка | ≈ |
| 0,75 | 6–7 | Толстая плёнка, близка к максимальному затуханию | ≈ |
| 0,70 | ≥8 | Очень толстая/вязкая плёнка, сильное затухание | ≈ |
Как видно из таблицы, наиболее выраженное затухание отражённого сигнала наблюдается при
[9][10]
4. Заключение
В работе представлен спектральный метод количественной оценки нефтяных загрязнений на морской поверхности, основанный на подавлении капиллярно-гравитационных волн. Разработанная модель позволила численно воспроизвести спектр отражённого сигнала при различных уровнях загрязнения и выявить ключевой диапазон чувствительности метода — при коэффициенте сглаживания
k ≈ 0.80