НАУКА и ОБРАЗОВАНИЕ

УДК 551.501

Россия, НИИ Радиоэлектроники и лазерной техники МГТУ им. Н.Э. Баумана

belov@bmstu.ru

ekomonit@bmstu.ru

gorod@bmstu.ru

В настоящее время актуальной является проблема загрязнения нефтепродуктами прибрежных зон морей (см., например [1-4]). Именно в прибрежных зонах морей ярко проявляются отрицательные последствия активной хозяйственной деятельности человека.

Наиболее перспективными методами контроля акваторий с целью обнаружения нефтяных загрязнений на поверхности воды являются методы, основанные на активной лазерной локации водной поверхности. Лазерные методы  не требуют естественного освещения, могут использоваться в широком диапазоне атмосферных условий независимо от времени суток, обладают высокой пространственно-угловой разрешающей способностью, возможностью фильтрации и временного стробирования полезного сигнала на фоне помех (солнечных бликов, импульсов обратного рассеяния от атмосферных слоев и т.п.)

            Физической основой дистанционного обнаружения нефтяных загрязнений на водной поверхности методом активного лазерного зондирования является наличие контраста яркости лазерного излучения отраженного от чистой водной поверхности и поверхности, покрытой пленкой нефтепродуктов. Контраст обусловлен двумя причинами [5, 6]: нефтяная пленка имеет другой коэффициент отражения, чем граница раздела воздух-вода, и  нефтяная пленка сглаживает морское волнение.

            Обнаружение нефтяных загрязнений обычно проводят, регистрируя лазерный локационный сигнал, отраженный от морской поверхности, и определяя контраст между чистой морской поверхностью и морской поверхностью, покрытой пленкой нефтепродуктов (см., например, [7-10]).

Эти методы позволяют надежно обнаруживать нефтяные загрязнения, однако они имеют недостаток – они обнаруживают нефтяные загрязнения в схеме зондирования морской поверхности вертикально вниз или в направлениях очень близких к вертикальным (т.е. полоса контроля нефтяных загрязнений на морской поверхности является очень узкой – единицы и или десятки).

Ниже описан угловой лазерный метод обнаружения нефтяных загрязнений на взволнованной морской поверхности, использующий угловое сканирование и позволяющий при большой вероятности правильного обнаружения проводить контроль нефтяных загрязнений в широкой полосе на морской поверхности.

Считаем, что лазерный локатор для зондирования нефтяных загрязнений морской поверхности установлен на авиационном носителе. Схема зондирования морской поверхности имеет узкий лазерный пучок, сканирующий поперек линии полета (см. рис. 1).

Рис. 1. Cхема зондирования морской поверхности

На рисунке Л – лазерный локатор, расположенный на авиационном носителе, S - морская поверхность,  - угол расходимости лазерного источника,  - текущий угол между нормалью к плоскости z=0 и оптической осью лазерного локатора. Режим подсвета может быть как непрерывным, так и импульсным.

В качестве измеряемой величины используется массив значений средней (строго – по ансамблю реализаций морской поверхности, в практических задачах – по участкам морской поверхности много большим размера ее неровностей) мощности эхо-сигнала , регистрируемого при пространственном сканировании (поперек направления полета носителя) морской поверхности. При импульсном подсвете  - мощность в момент максимума регистрируемого импульсного локационного сигнала.

Мощность эхо-сигнала , регистрируемого приемником от взволнованной морской поверхности, определяется (импульсный подсвет при ; небольшие углы зондирования , ) следующей формулой (см., например, [8]):

                                 (1)

где:

;

;

L - расстояние от лидара до морской поверхности;  - коэффициент отражения Френеля для морской поверхности; - длительность зондирующего импульса; ,  - дисперсии наклонов морской поверхности поперек и вдоль линии полета  (в плоскостях  XOZ и YOZ.);

            Для оценки дисперсий наклонов морской поверхности без нефтяных загрязнений используются следующие соотношения [11, 12]:

;  ;

 ,

 ,  - дисперсии наклонов морской поверхности по направлению ветра и поперек направления ветра;  - угол между осью x и направлением ветра;  - скорость приводного ветра (м/с).

            Для прозрачной аэрозольной атмосферы величины , ,  равны:

;  ;  ;

 - угол расходимости источника излучения и угловое поле зрения приемной системы;  - мощность излучаемая источником;  - эффективный размер приемной апертуры;  - оптическая толщина атмосферы между лазерным локатором и морской поверхностью.

Из формулы (1) получим формулу для отношения :

              (2)

Если  мало и измеряемые величины ,  для соседних участков морской поверхности относятся к одной однородной области (воды или нефти), то , ,  и из (2) получим:

                    (3)

            После преобразований из (3) имеем:

                                          (4)

            Таким образом, проводя измерения принимаемой мощности для небольших соседних участков морской поверхности, относящихся к одной однородной области (воды или нефти) можно измерить дисперсию  наклонов морской поверхности в этой области (а значит и в пределах всего сектора сканирования на морской поверхности).

Отметим, что если для соседних участков измеряемые величины ,  не относятся к одной однородной области (воды или нефти), то они будут сильно отличаться по величине (а значит, могут быть отбракованы при обработке данных измерений).

Для принятия решения о наличии нефтяных загрязнений на морской поверхности вычисляется  - отношение величины  исследуемого участка и чистой морской поверхности (величина  позволяет оценить эффект сглаживания морского волнения из-за нефтяных загрязнений):

,                                               (5)

где:

 - дисперсии наклонов исследуемого (текущего – регистрируемого при угле зондирования ) участка и чистой морской поверхности; величина  дисперсии наклонов чистой морской поверхности может либо вычисляться в бортовом вычислителе по данным измерений скорости приводного ветра, либо измеряться при полете над заведомо чистым (без нефтяных загрязнений) участком морской поверхности.

Контроль величины  позволяет относить исследуемый участок морской поверхности к случаю чистой воды или нефтяного загрязнения. Однако для надежного обнаружения нефтяных загрязнений одной величины  недостаточно (см., например, [8, 9]) – кроме сглаживания морского волнения из-за разлива нефтепродуктов необходимо контролировать и коэффициент отражения поверхности.

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

                         (6)

            Используя (1) из (6) имеем:

,                                               (7)

где:

;

 – величина, не зависящая от характеристик морской поверхности; .

            В левой части выражения (7) стоит величина, зависящая от данных измерений , результатов определения величины  (по данным обработки результатов измерений), дальности зондирования, характеристик источника и приемника. В правой части (7) стоит коэффициент отражения  морской поверхности (распределение которого в пределах всего участка сканирования на морской поверхности мы хотим определить). Величина  в правой части (7) умножена на некоторый  коэффициент , однако это не мешает контролировать  и относить полученную оценку  к случаю чистой воды или нефтяного загрязнения. Величина  увеличивается при переходе от чистой воды к нефтяному загрязнению в 2 раза в видимом и ближнем ИК диапазонах и в 4 раза в среднем ИК диапазоне, в то время как коэффициент  довольно слабо меняется при изменении характеристик волнения (см. рис. 2).

Рис. 2. Зависимость коэффициента  от  скорости приводного ветра

На рисунке 2 показана зависимость коэффициента  от скорости приводного ветра  в широком диапазоне условий волнения на морской поверхности  - от штиля (0,5 м/с) до шторма (18 м/с). Хорошо видно, что, по крайней мере, в среднем ИК диапазоне вариации коэффициента  не помешают отнести измеренную величину  к чистой воды или нефтяному загрязнению.

Если для принятия решения о наличии нефтяных загрязнений на морской поверхности использовать  (- отношение величины  исследуемого участка и чистой морской поверхности), то величина  вообще не будет зависить от коэффициента  (считаем, что для нефтяных пленок дисперсии наклонов  одинаково уменьшаются по сравнению с чистой водной поверхностью):

                        (8)

где:

 и - значение величины  и коэффициента отражения для исследуемого участка и чистой морской поверхности.

          Таким образом, величина  содержит информацию только о дисперсии наклонов исследуемого участка морской поверхности, а величина  - только о коэффициенте отражения исследуемого участка морской поверхности.

Для принятия решения о наличии нефтяных загрязнений на исследуемом (текущем – при угле зондирования ) участке морской поверхности необходимо определить величины  и  и сравнить их с некоторыми пороговыми значениями.

          На рисунках 3, 4 представлены зависимости величин ,  от скорости приводного ветра  в идеальном случае отсутствия шума измерения.

Рис. 3. Зависимость величины  от скорости приводного ветра.

Рис. 4. Зависимость величины  от скорости приводного ветра.

На рисунках 5, 6 представлены результаты математического моделирования – пространственные реализации величин , .

Рис. 5. Пространственная реализация величины .

Рис. 6. Пространственная реализация величины .

Для рисунков 5, 6 по оси абсцисс отложен i - номер пространственного дискрета на линии сканирования; среднеквадратическое значение шума измерения равно 0,5 %; =0,12; =0,003; скорость приводного ветра – 3 м/с; высота полета носителя – 2000 м; длина волны лазера 1,06 мкм; ; чистая морская поверхность: ; толстая нефтяная пленка: . Пространственные реализации величин ,  получены по данным «измерений» по описанному выше методу. Принималось, что для толстых нефтяных пленок дисперсии наклонов  уменьшаются по сравнению с чистой водной поверхностью в 3 раза [6, 11], а для слика (участка морской поверхности со сглаженным ветровым волнением) - в 10 раз. В области слика по разным причинам происходит очень сильное сглаживание волнения на морской поверхности и его можно спутать с нефтяным загрязнением, однако коэффициент отражения в области слика такой же как и у чистой поверхности воды.

На рисунках 3 и 5 линии 1, 2, 3 соответствуют величинам  “нефтяное загрязнение - чистая морская поверхность”, “слик - чистая морская поверхность”, “чистая морская поверхность - чистая морская поверхность”, линии 4 и 5 - пороговые значения =1,5, =5. На рисунках 4 и 6 обозначения те же (как и на рисунках 3, 5) с одним отличием:  кривая 5  соответствует пороговому значению =3.

            Из рисунков 3 - 6 хорошо видно, что определяя величины ,  можно уверенно обнаруживать нефтяные загрязнения и отличать их от сликов на морской поверхности.             Математическое моделирование показывает, что одновременное использование условий  и  позволяет принимать решения о наличии нефтяных загрязнений на морской поверхности с вероятностью правильного обнаружения 0,9935 при вероятности ложных тревог 0,0025. Причем, при используемых для расчета параметрах схемы зондирования ширина полосы контроля нефтяных загрязнений на морской поверхности была не менее 500 м.

Таким образом, угловой лазерный метод обнаружения нефтяных загрязнений на взволнованной морской поверхности, позволяет при большой вероятности правильного обнаружения и приемлемой вероятности ложных тревог проводить обнаружение нефтяных загрязнений в широкой (сотни метров) полосе на морской поверхности.

Литература

1. Патин С.А. Нефть и экология континентального шельфа. М.: ВНИРО, 2001. 247 с.

2. Техногенное загрязнение природных вод углеводородами и его экологические последствия / В.М. Гольдберг [и др.] М.: Наука, 2001. 125 с.

3. Алешин И.В. Экологический мониторинг Мирового океана. СПб.: Изд-во СПГМТУ, 1997. 76 с.

4. Современное состояние прибрежных экосистем морей Российской Федерации / Ю.А.Израэль [и др.] // Метеорология и гидрология. 1995. N9. С. 6-21.

5. Гуревич И.Я., Шифрин К.С. Отражение видимого и ИК излучения нефтяными пленками на море // Оптические методы изучения океанов и внутренних водоемов. Новосибирск: Наука, 1979. С. 166-176.

6. Гардашов Р.Г., Гуревич И.Я., Шифрин К.С. Отражение оптического излучения от взволнованной морской поверхности, покрытой нефтяной пленкой // Оптика атмосферы и океана. Баку: ЭЛМ, 1983. С. 33-44.

7. Оптико-физические средства исследования океана / Под ред. Е.Г. Пащенко. Л.: Судостроение, 1984. 264 с.

8. Белов М.Л., Городничев В.А., Козинцев В.И. Трехлучевой лазерный метод обнаружения нефтяной пленки на морской поверхности. // Вестник МГТУ. сер. Приборостроение. 2000. № 3. С. 7-12.

9. Дистанционное обнаружение нефтяных загрязнений на взволнованной морской поверхности с помощью двухспектрального метода / Белов М.Л. [и др.]  // Оптика атмосферы и океана. 2006.Т.19. N10. С. 872-874.

10. Белов М.Л. Спектрально-угловой лазерный метод дистанционного обнаружения нефтяных загрязнений на морской поверхности // Информационно-измерительные и управляющие системы. 2005. Т.3. N4-5. С.60-64.

11. Cox C.,Munk W. Measurement of the roughness of the sea surface from photographs of the sun's glitter // J.O.S.A. 1954. V.44. N.11. P.838-850.

12. Басс Ф.Г., Фукс И.М. Рассеяние волн на статистически неровной поверхности. М.:Наука, 1972. 424 с.