Detection of oil spills on the sea surface based on analysis of the spectrum of reflected radio signals from coastal radar stations
Detection of oil spills on the sea surface based on analysis of the spectrum of reflected radio signals from coastal radar stations
Abstract
The article suggests a spectral method for the detection and quantification of oil pollution on the sea surface, based on the analysis of changes in the spectrum of the reflected radar signal from coastal radar stations when capillary-gravitational waves are suppressed by an oil film. The method is implemented using an adapted mathematical model of radar analysis that takes into account the physical properties of the film and the conditions of Bragg scattering. Numerical simulations were carried out for parameters typical of X-range radar systems, followed by an evaluation of the dependence of the spectral amplitude on the degree of microwave smoothing of the surface. It is shown that the spectral amplitude varies non-monotonically as the film thickens: at a certain thickness, the method’s sensitivity reaches a maximum, after which its effectiveness decreases due to saturation and interference distortions. A generalised table has been compiled linking the physical characteristics of the film to changes in spectral parameters, which makes it possible to apply the method for remote monitoring of pollution levels. The results are consistent with recent research in the field of radar, hyperspectral and infrared remote sensing, and the proposed approach can be implemented at fixed coastal radar stations as a cost-effective tool for regular monitoring of environmental conditions and analysis of the impact of oil slicks on atmospheric processes and the formation of various weather phenomena.
1. Введение
Нефтяные загрязнения представляют серьёзную угрозу для морских экосистем, качества водной среды и здоровья человека. Также наличие нефтяной пленки может оказать местное влияние на газо-, тепло- и влагообмен системы океан-атмосфера, что необходимо учитывать при численном моделировании атмосферных процессов. При прогнозировании таких явлений погоды, как осадки, туман, низкой слоистой облачности, необходимо учитывать, что частицы нефти, попадая в атмосферу могут представлять собой дополнительные ядра конденсации. Актуальность исследования обусловлена высокой уязвимостью прибрежных и мелководных зон, где даже локальные утечки нефти могут привести к долговременным негативным последствиям. В таких районах необходимы непрерывное наблюдение состояния окружающей среды и разработка доступных, надёжных методов обнаружения нефтяной плёнки
.Цель настоящей работы заключается в разработке и апробации спектрального метода количественной оценки нефтяных загрязнений морской поверхности, основанного на адаптированной математической модели радиолокационного анализа, ранее применявшейся авторами в иной постановке задачи. Проведённые численные эксперименты показали согласованность предложенного подхода с результатами наблюдений
, .Метод базируется на ключевом физическом признаке нефтяных плёнок — подавлении капиллярно-гравитационных волн, что приводит к изменению спектральных характеристик отражённого радиосигнала. Эффективность такого подхода подтверждается результатами современных исследований с применением радиолокационного (SAR) зондирования
, гиперспектрального мониторинга с беспилотных летательных аппаратов и инфракрасного спектрального анализа .Научная новизна работы состоит в применении спектрального анализа для количественного определения характеристик нефтяных плёнок, теоретическая значимость — в адаптации модели к задачам дистанционного мониторинга загрязнённой поверхности, практическая — в возможности реализации метода на стационарных прибрежных радиолокационных станциях как экономичного и локально привязанного инструмента.
2. Методы и принципы исследования
2.1. Физические основы метода
При попадании нефтепродуктов на морскую поверхность формируется тонкая плёнка, которая подавляет капиллярно-гравитационные волны — основной источник рассеяния в радиолокационном диапазоне X-band (около 3 см). Это приводит к снижению отражённой мощности радиосигнала, особенно при углах облучения 30–45°, что обусловлено ростом вязкости и снижением поверхностного натяжения; максимальная чувствительность достигается при условии брегговского рассеяния, когда длина радиоволны соизмерима с длиной поверхностных волн
.
2.2. Исходные условия численного эксперимента
В рамках численного моделирования были заданы параметры, характерные для прибрежной радиолокационной системы, размещённой на стационарной вышке. Высота антенны составляла 20 метров, углы облучения — 30–45°, что обеспечивает охват участка поверхности длиной около 10 метров. Рабочая частота радара — 10 ГГц (X-band). Схема размещения и область облучения представлены на рис. 1.
Рисунок 1 - Схема радиолокационного облучения с вышки
2.3. Математическая модель
Для количественного анализа спектра отражённого сигнала была использована модель, разработанная авторами, в которой морская поверхность представляется как совокупность независимых флуктуирующих рассеивателей, равномерно распределённых вдоль направления излучения. Шаг между ними определяется из длины облучаемого участка L и числа рассеивателей. Расстояние определяется углом падения излучения. Исходный сигнал моделируется в виде гармонической волны, отражённой от всей поверхности. Суммарный сигнал формируется как суперпозиция откликов от всех рассеивателей с учётом фазовых сдвигов. Его аналитическое выражение имеет следующий вид:
где
c — скорость распространения электромагнитной волны;
Для перехода к спектру мощности используется автокорреляционная функция и теорема Винера–Хинчина. Влияние учитывается коэффициентом сглаживания Значения k и соответствующие им характеристики плёнки были заданы по данным
.3. Основные результаты
Для первичной оценки влияния загрязнения на отражённый радиосигнал была выполнена серия численных экспериментов, сравнивающих спектр при наличии нефтяной плёнки и в условиях чистой воды. На рис. 2 представлены спектры в двух случаях: эталон (без загрязнения) и загрязнённая поверхность при коэффициенте сглаживания k = 0,95, соответствующем умеренному загрязнению.
Рисунок 2 - Подавление высокочастотных компонент спектра при наличии нефтяной плёнки
Для более детального анализа проведено параметрическое исследование зависимости спектральной амплитуды от степени сглаживания поверхности. Это позволяет установить, как различные уровни загрязнения влияют на отражённый сигнал в спектральной области.
На рис. 3 представлена нормализованная амплитуда спектра в зависимости от коэффициента сглаживания k, моделирующего степень микроволнового подавления.
Рисунок 3 - Зависимость амплитуды спектра от коэффициента сглаживания
Для получения количественного соответствия между уровнем сглаживания и физическими характеристиками загрязнения была построена табл. 1. Она объединяет результаты численного моделирования спектральной амплитуды с данными из работы
, где установлена зависимость между коэффициентом подавления ряби и толщиной плёнки. На основе этой связи удалось отразить, как изменение k влияет одновременно и на спектр отражённого сигнала, и на физическую структуру загрязнения.Таблица 1 - Соответствие толщины плёнки и коэффициента подавления ряби
Коэффициент подавления ряби | Примерная толщина нефтяной плёнки, мкм | Характеристика плёнки | Нормализованная амплитуда спектра (по данным графика) |
1,00 | 0 | Чистая вода, без загрязнения | 1,00 |
0,97 | 0,5–1 | Очень тонкая плёнка, минимальное воздействие | ≈ 0,87 |
0,95 | 1–2 | Начало подавления капиллярных волн | ≈ 0,75 |
0,90 | 2–3 | Умеренное подавление | ≈ 0,55 |
0,85 | 3–4 | Существенное снижение микроволнения | ≈ 0,40 |
0,80 | 4–5 | Почти полное подавление капиллярных волн | ≈ 0,05 |
0,78 | 5–6 | Плотная нефтяная плёнка | ≈ 0,03 |
0,75 | 6–7 | Толстая плёнка, близка к максимальному затуханию | ≈ 0,15 |
0,70 | ≥8 | Очень толстая/вязкая плёнка, сильное затухание | ≈ 0,55 (возможно, локальный минимум) |
Как видно из таблицы, наиболее выраженное затухание отражённого сигнала наблюдается при k ≈ 0,80, что соответствует толщине плёнки порядка 4–5 мкм. Именно в этом диапазоне достигается максимум чувствительности метода, поскольку капиллярные волны максимально подавлены, а интенсивность отражения радиоволн от поверхности минимальна. При дальнейшем снижении k, несмотря на увеличение толщины плёнки, амплитуда спектра начинает возрастать, что связано с потерей чувствительности спектрального метода. Такой эффект подтверждает, что метод надёжно работает лишь в ограниченном диапазоне: при слишком тонкой плёнке влияние на спектр слабо выражено, а при чрезмерной толщине — спектральная картина искажается из-за насыщения и интерференции, как показано в
, .4. Заключение
В работе представлен спектральный метод количественной оценки нефтяных загрязнений на морской поверхности, основанный на подавлении капиллярно-гравитационных волн. Разработанная модель позволила численно воспроизвести спектр отражённого сигнала при различных уровнях загрязнения и выявить ключевой диапазон чувствительности метода — при коэффициенте сглаживания k ≈ 0.80, соответствующем толщине плёнки 4–5 мкм. Полученные зависимости спектральной амплитуды согласуются с результатами современных исследований и отражают нелинейную природу взаимодействия радиоволн с загрязнённой поверхностью. Итоговая таблица связывающая коэффициент подавления с физическими параметрами плёнки, подтверждает, что метод позволяет оценивать степень загрязнения на основе радиолокационных данных. При этом зафиксированы ограничения спектрального подхода: при избыточной толщине наблюдается эффект насыщения, приводящий к снижению чувствительности. Метод может применяться на прибрежных радиолокационных станциях для организации непрерывных наблюдений за состоянием окружающей среды в целях обеспечения экологической безопасности, а также исследования влияния нефтяных пленок на развитие атмосферных процессов и формирования различных явлений погоды.