Министерство образования и наук
и
Федеральное агентство по образованию
ГОУ ВПО
РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ
УНИВЕРСИТЕТ
(РГГМУ)
УДК 551.510.7
На правах рукописи
Саноцкая Надежда Александровна
ВОССТАНОВЛЕНИЕ ХАРАКТЕРИСТИК АТМОСФЕРНОГО АЭРОЗОЛЯ
ПО ДАННЫМ ОПТИЧЕСКИХ ИЗМЕРЕНИЙ
Специальность 25.00.30 – метеорология, климатология и агрометеорологияАВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук Санкт-Петербург
Работа выполнена в ГОУВПО «Российский государственный гидрометеорологический университет»
Научный руководитель: доктор физико-математических наук Егоров Александр Дмитриевич
Официальные оппоненты: доктор физико-математических наук, профессор Ивлев Лев Семенович кандидат физико-математических наук, доцент Кашлева Лариса Владимировна ГУ «Главная геофизическая
Ведущая организация:
обсерватория им. А.И. Воейкова»
Защита состоится «25» февраля 2010 года в 1530 на заседании диссертационного совета Д 212.197.01 при ГОУВПО «Российский государственный гидрометеорологический университет» по адресу:
195196, Санкт-Петербург, Малоохтинский пр., 98.
С диссертацией можно ознакомиться в библиотеке ГОУВПО «Российский государственный гидрометеорологический университет».
Автореферат разослан « » января 2010 года.
Ученый секретарь диссертационного совета д-р физ.- мат. наук, профессор А.Д.Кузнецов
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Решение проблемы совершенствования методов расчета в области атмосферной оптики в настоящее время приобретает особую актуальность в связи с постоянно увеличивающейся загрязненностью атмосферы и со сложностью параметризации оптико-микроструктурных свойств аэрозольных частиц.
В настоящее время в практике определения концентраций загрязняющих веществ в воздухе используются, в основном, контактные методы измерений. В частности, оптическими счетчиками частиц и фильтровыми аспирационными устройствами измеряются параметры аэрозоля в непосредственной близости от этих приборов. Применение оптических методов позволяет автоматизировать процесс измерений.
Недостаток оптических методов состоит в том, что они не дают возможности непосредственного определения характеристик загрязнения атмосферы.
Исследование проблемы интерпретации результатов, получаемых в процессе определения характеристик атмосферного аэрозоля оптическими методами, позволит усовершенствовать эти методы.
Цель работы – повышение достоверности результатов интерпретации оптической информации об атмосферном аэрозоле на основе применения алгоритмов обработки экспериментальных результатов, базирующихся на теории рассеяния света.
Основные задачи исследования:
- анализ особенностей оптических методов, применяемых для определения параметров атмосферного аэрозоля;
- оценка систематических погрешностей, существенно влияющих на результаты оптических измерений, с учетом особенностей измерительной аппаратуры;
- разработка модели сигнала направленного рассеяния с учетом экспериментальных данных;
- оценка эффективности методов оптических измерений параметров атмосферного аэрозоля, предназначенных для повышения достоверности результатов измерений.
Методы исследования базируются на решении уравнений Максвелла.
Используются методы численного анализа, компьютерное моделирование с применением данных натурных экспериментов, статистический анализ.
Научная новизна работы:
- определено, что направленное рассеяние излучения частицей с приведенным размером, значительно превышающим единицу, может существенно зависеть от структуры частицы;
- установлено практическое отсутствие зависимости направленного рассеяния излучения частицей с приведенным размером порядка единицы от структуры частицы;
направленного рассеяния излучения частицей от структуры частицы, удается объяснить имеющиеся экспериментальные данные;
- установлено уменьшение зависимости направленного рассеяния излучения частицей от структуры частицы при уменьшении угла рассеяния.
Основные положения, выносимые на защиту:
- зависимость направленного рассеяния излучения частицей с приведенным размером, значительно превышающим единицу, от структуры частицы;
- отсутствие зависимости направленного рассеяния излучения частицей с приведенным размером порядка единицы от структуры частицы;
- основы метода оптических измерений аэрозольных характеристик, базирующегося на уменьшении направленного рассеяния излучения частицей при уменьшении угла рассеяния зависимости от структуры частицы.
Достоверность полученных результатов подтверждается хорошими результатами сопоставления теоретических и экспериментальных результатов определения характеристик атмосферного аэрозоля, а также использованием современных физико-математических методов анализа и обобщения данных.
Практическая значимость работы определяется, в соответствии с поставленной целью, повышением достоверности результатов интерпретации оптической информации об атмосферном аэрозоле. Полученные результаты использованы в учебном процессе в РГГМУ и могут быть использованы для целей совершенствования оптических технических средств, предназначенных для мониторинга загрязнения атмосферы.
Личный вклад автора Все основные результаты получены автором лично. Автор выполнил анализ современного состояния проблемы. Нашел новые результаты в области исследования рассеяния излучения частицей.
Апробация работы. Результаты исследований докладывались на международном симпозиуме «33 International symposium on remote sensing of environment», 2009, на Международном симпозиуме стран СНГ «Атмосферная радиация и динамика» МСАР, 2009, на итоговых сессиях ученого совета РГГМУ, 2008, 2009.
Публикации. Основные результаты исследований опубликованы в 5 статьях.
Структура и объем работы. Работа состоит из введения, 3 глав, заключения, списка литературы - 120 наименований, списка обозначений и сокращений. В ней содержится 121 страница текста, 8 таблиц, 25 рисунков.
СОДЕРЖАНИЕ РАБОТЫ
Во введении анализируется современное состояние рассматриваемой проблемы, обосновывается актуальность темы диссертации, указывается цель работы, формулируются основные задачи исследования и положения, выносимые на защиту, отмечается научная новизна работы, ее практическая значимость, аргументируется обоснованность и достоверность полученных результатов, дается краткая характеристика диссертации.
В первой главе осуществляется физическая и математическая постановка задачи интерпретации данных, получаемых в процессе оптических измерений характеристик атмосферного аэрозоля.
Задача установления оптико-микроструктурных связей аэрозольных частиц достаточно сложна из-за существования различных аэрозольных фракций. При ее решении следует учитывать существование различных источников аэрозоля, поставляющих в атмосферу частицы, различающиеся по размеру, химическому составу, оптическим характеристикам. Среди других факторов, усложняющих задачу, можно выделить значительную изменчивость микрооптических характеристик частиц, не укладывающуюся в рамки теории рассеяния света однородной частицей. В связи с этим в работе осуществляется моделирование микрооптических свойств частиц атмосферного аэрозоля с радиально изменяющимся показателем преломления в оболочке, покрывающей однородное ядро.
Во второй главе разрабатываются методы интерпретации данных оптических измерений.
Математическое моделирование рассеяния частицей электромагнитных волн с круговой частотой = kc и напряженностями электрического Eexp(it) и магнитного Hexp(it) полей предполагает решение относительно векторов E, H системы уравнений Максвелла:
где k – волновое число, c – скорость света, m = n - i, – комплексный показатель преломления частицы.
Обозначим через r1 радиус частицы, r – расстояние от центра, = kr – волновое расстояние, 1 = kr1 – приведенный радиус частицы, m1 – показатель преломления атмосферы. Показатель преломления m() частицы задается функцией от (0 1); частица характеризуется показателем преломления в виде:
С целью решения этой системы уравнений для моделирования рассеяния плоских волн облучаемыми частицами в работе вводятся потенциалы Дебая D, B, с которыми компоненты напряженностей, электрического и магнитного полей в сферической системе координат связаны известными соотношениями.
Уравнения, которыми описываются потенциалы D, B, можно записать в виде:
результате:
где P j 1 (cos ) – присоединенная функция Лежандра, функции D j (), B j () удовлетворяют уравнениям:
В работе находятся численные решения уравнений (5), что позволяет аэрозольными частицами. Эти уравнения решаются для неоднородной оболочки, и используется представление потенциалов плоской волны.
На поверхностях раздела: ядро – оболочка, оболочка – среда задаются составляющих векторов напряженностей полей.
Для решения уравнений (5) для потенциалов D, B их можно преобразовать введением функций:
в дифференциальные уравнения Риккати:
дифференциальных уравнений (7).
Эти уравнения решаются для неоднородной оболочки (0 1), причем задаются следующие начальные условия:
учитывающие, что величина m(0 - 0) – показатель преломления ядра, может не совпадать со значением m(0), где 0 – приведенный радиус ядра частицы.
Исходные для параметров рассеянной волны коэффициенты aj, bj находятся по формулам:
которые получаются с учетом граничных условий.
E0 – амплитуда вектора электрического поля облучающей волны, соотношениям следующего вида:
Производные этих функций выражаются через сами функции. Для обеих функций имеет место формула:
Для расчета коэффициента направленного рассеяния, зависящего от вектора Пойнтинга, вводится сумма (i1 + i2), пропорциональная этому коэффициенту, где Sj, Qj – угловые функции, связанные с присоединенными функциями Лежандра:
Факторы эффективности ослабления и рассеяния находятся по формулам:
В третьей главе анализируются данные оптических измерений характеристик атмосферного аэрозоля, осуществляется моделирование направленного рассеяния, приводятся результаты моделирования.
рассеяния 1,57, показатель преломления в ядре из сажи m(0) = 1,82-0,64i, приведенный размер частицы 30) представлены на рис. 1, 2. Здесь L = k(r1 - r0) – приведенная толщина покрытия, IR – индикатриса рассеяния, отнесенная к её значению для L = 0, j =0,...,4.
Рис. 1. Результаты анализа направленного рассеяния излучения Рис.2. Результаты анализа направленного рассеяния излучения означает присутствие воды в оболочке, j = 3,4 – ее отсутствие).
Перечень характеристик рассеивающего центра представлен в табл. 1, где d = ( - 0)/(1 - 0).
Таблица 1 – Характеристики рассеивающего центра Как следует из рис. 1, 2 индикатриса направленного рассеяния может существенно изменяться при изменении структуры аэрозольной частицы (в частности, при ее обводнении). Характер модельных изменений структуры рассеивающего центра (j = 4 – табл. 1, L = 2,0 - рис. 2) соответствует экспериментальным результатам.
Результаты расчетов относительной величины IR(1)/IR(1,33) (отношения индикатрисы рассеяния для покрытия, не содержащего воду, когда m1 = 1, к индикатрисе рассеяния для покрытия, содержащего воду, когда m1 = 1,33) для разных величин при заданном отношении g = (r1 - r0)/r показывают существенную изменчивость величины IR.
На рис. 3 представлено значение G величины g, соответствующее максимальному уменьшению относительной индикатрисы рассеяния при наличии покрытия (рис. 4).
Рис. 3. Значение G величины g, соответствующее IR(1)/IR(1,33) Рис. 4. Максимальное уменьшение индикатрисы Таким образом, проблема интерпретации данных оптических измерений распределений частиц по размерам достаточно сложна.
С другой стороны, известен результат, когда в лабораторных условиях удалось радикально уменьшить влияние неоднородности частицы на результаты измерений, уменьшая угол рассеяния до 30о.
характеристик рассеивающего центра, представленных в табл. 1, для частиц с приведенным радиусом ядра Ra = 30.
Здесь I, I(j) – индикатрисы рассеяния однородной и неоднородной частицей сажи (R1 = 32). Их отношение I(j)/I может быть очень малой величиной (несколько процентов для j = 4) при рассеянии под углом = 90о.
Оно существенно увеличивается при меньших углах рассеяния (составляет 30 % при = 30о для j = 4).
Найденное изменение, однако, существенно превышает вариации индикатрисы для различных неоднородностей частиц, найденные в работе в лабораторных условиях.
Таблица 2 – Относительная индикатриса для различных Настоящее исследование относится к условиям натурных измерений в атмосфере, когда оптические свойства частиц изменяются в существенно более широких пределах. Тем не менее, и при измерениях характеристик атмосферного аэрозоля можно достичь существенного уменьшения влияния оптических свойств частиц на результаты измерений за счет выбора угла рассеяния излучения.
Важно подчеркнуть, что в рассматриваемом диапазоне углов имеет место факт меньшего влияния оптических свойств частиц на направленное рассеяние по сравнению с их влиянием на ослабление излучения. Кроме того, влияние меньше для частиц меньших размеров. Особенно это касается частиц аккумулятивной фракции. Данное обстоятельство, а также относительная стабильность распределения числа частиц атмосферного аэрозоля по размерам объясняют существование корреляции между коэффициентом ослабления и концентрацией частиц. Это, в свою очередь, расширяет возможности методов оптических измерений, в частности, лидарного зондирования атмосферного аэрозоля.
В заключении сформулированы основные выводы работы:
- найдено, что направленное рассеяния излучения частицей с приведенным размером, значительно превышающим единицу, может существенно зависеть от структуры частицы;
- установлено практическое отсутствие зависимости от структуры частицы направленного рассеяния излучения частицей с приведенным размером порядка единицы;
- найденными результатами, описывающими зависимость от структуры частицы направленного рассеяния излучения частицей, удается объяснить имеющиеся экспериментальные данные;
- установлено уменьшение зависимости от структуры частицы направленного рассеяния излучения частицей при уменьшении угла рассеяния.
ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ
Егоров А.Д., Ржонсницкая Ю.Б., Потапова И.А., Саноцкая Н.А.Определение характеристик атмосферного аэрозоля по данным оптических измерений // Естественные и технические науки, 2009, №1, с.44 – 54.
Егоров А.Д., Потапова И.А., Ржонсницкая Ю.Б., Саноцкая Н.А. Методы лидарного определения характеристик атмосферного аэрозоля, Сборник тезисов Международного симпозиума стран СНГ «Атмосферная радиация и динамика» (МСАРД – 2009) 22 – 45.
Егоров А.Д., Перельман А.Я., Куликов В.Н., Саноцкая Н.А., микроструктуры по данным оптических измерений, Сборник тезисов Международного симпозиума стран СНГ «Атмосферная радиация и динамика» (МСАРД – 2009) 22 – 26 июня 2009 г., Санкт – Петербург, СПбГУ, МСАРД, 2009, с. 84 – 85.
Sanotskaya N.A. Remote Sensing of Atmospheric Aerosols, International Symposium on Remote Sensing of Environment, Stresa, Italy, 2009 [electronic resource] http://www.isprs.org/publications/related/ISRSE/html/papers/191.pdf.
Yegorov A.D., Potapova I.A., Sanotskaya N.A. Lidar probing of Atmospheric Aerosols, 33 International Symposium on Remote Sensing of Environment, Stresa, Italy, 2009 [electronic resource] http://www.isprs.org/publications/related/ISRSE/html/papers/898.pdf.