WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

загрузка...
   Добро пожаловать!

Pages:     || 2 |
Министерство образования и наук

и Федеральное агентство по образованию ГОУ ВПО РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ (РГГМУ) УДК 551.510.7

На правах рукописи

Саноцкая Надежда Александровна ВОССТАНОВЛЕНИЕ ХАРАКТЕРИСТИК АТМОСФЕРНОГО АЭРОЗОЛЯ ПО ДАННЫМ ОПТИЧЕСКИХ ИЗМЕРЕНИЙ Специальность 25.00.30 – метеорология, климатология и агрометеорология

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук Санкт-Петербург 2010

Работа выполнена в ГОУВПО «Российский государственный гидрометеорологический университет»

Научный консультант: доктор физико-математических наук Егоров Александр Дмитриевич

Официальные оппоненты: доктор физико-математических наук, профессор Ивлев Лев Семенович кандидат физико-математических наук, доцент Кашлева Лариса Владимировна ГУ «Главная геофизическая

Ведущая организация:

обсерватория им. А.И. Воейкова»

Защита состоится «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 системы уравнений Максвелла:

rot H = ikm2E, (1) rot E = - ikH, где k – волновое число, c – скорость света, m = n - i, – комплексный показатель преломления частицы.

Обозначим через r1 радиус частицы, r – расстояние от центра, = kr – волновое расстояние, 1 = kr1 – приведенный радиус частицы, m1 – показатель преломления атмосферы. Показатель преломления m() частицы задается функцией от (0 1); частица характеризуется показателем преломления в виде:

m (0) = m0 0 0, m () 0 < 1, (2) m (1) = m1 > 1.

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

Уравнения, которыми описываются потенциалы D, B, можно записать в виде:

-1 - (D cos ) + (D m2 - 2 m / m) cos = 0, (3) -1 - (B sin ) + Bm2 sin = 0.

Уравнения (3) решаются методом разделения переменных, в результате:

D = Dj() Pj1(cos ), (4) B = Bj() Pj1(cos ), где Pj1(cos ) – присоединенная функция Лежандра, функции Dj(), Bj() удовлетворяют уравнениям:

d -/ m2 ( D m2) + (m2 – j(j+1) ) Dj = 0, j d (5) - B + (m2 – j(j+1) ) Bj = 0.

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

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

Для решения уравнений (5) для потенциалов D, B их можно преобразовать введением функций:

j D () Yj ()= m2() D (), j (6) j B () Z ()= j, B () j в дифференциальные уравнения Риккати:

j( j +1) j ()=-m2 ()Yj2 ()+ 2m2 ()-1, (7) j(j +1) j ()=-2 ()+ -m2() j.

Yj (), Z () находятся путем численного решения j Функции дифференциальных уравнений (7).

Эти уравнения решаются для неоднородной оболочки (0 1), причем задаются следующие начальные условия:

j (m00) Y ( )=, j m0 j (m00 ) (8) j (m00 ) (0)= mj, j (m00 ) учитывающие, что величина m(0 - 0) – показатель преломления ядра, может не совпадать со значением m(0), где 0 – приведенный радиус ядра частицы.

Исходные для параметров рассеянной волны коэффициенты aj, bj находятся по формулам:

m1Yj(1)j (m1 )-j(m1 ) 1 aj =-, j m1Yj(1) (m1 )- (m1 ) j 1 j (9) ( ) (m1 )- m1 (m1 ) j 1 j 1 j b = -, j j ( ) (m1 )- m1 (m1 ) j 1 j 1 j которые получаются с учетом граничных условий.

Здесь aj =(-i)j E0, = j, j (11) mE0 – амплитуда вектора электрического поля облучающей волны, = + i. (11) j j j () удовлетворяют рекуррентным Функции (), j j соотношениям следующего вида:

2 j + (m) = (m) - -. (12) j+1 j j m Производные этих функций выражаются через сами функции. Для обеих функций имеет место формула:

j (m) = j -1 - j j (m).

(13) m При этом (m) = cos (m), -(m) = - sin (m), -(14) (m) = sin (m), (m) = cos (m).

0 Для расчета коэффициента направленного рассеяния, зависящего от вектора Пойнтинга, вводится сумма (i1 + i2), пропорциональная этому коэффициенту, где 2 j +(ajS + bjQj) j i1 =, (15) j(j +1) 2 j +(ajQj + bjS ) j i2 =, (16) j(j +1) Sj, Qj – угловые функции, связанные с присоединенными функциями Лежандра:

2 (j +1) Qj+2 -(j + 2) Qj Sj+1 =, 2 j +(2 j +3)Qj+1cos -(j + 2)Qj Qj+2 =, (17) j +Q1 = 1, Q2 = 3 cos, S1= cos.

Факторы эффективности ослабления и рассеяния находятся по формулам:

k = 2(Re j +1)aj + Re j +1)bj ) / 12. (18) (2 (1 2 ks = 2( (a + bj )(2 j +1))/ j. (19) В третьей главе анализируются данные оптических измерений характеристик атмосферного аэрозоля, осуществляется моделирование направленного рассеяния, приводятся результаты моделирования.

Результаты анализа направленного рассеяния излучения (угол рассеяния 1,57, показатель преломления в ядре из сажи m(0) = 1,82-0,64i, приведенный размер частицы 30) представлены на рис. 1, 2. Здесь L = k(r1 - r0) – приведенная толщина покрытия, IR – индикатриса рассеяния, отнесенная к её значению для L = 0, j =0,...,4.

IR 0,0,0,00,511,L Рис. 1. Результаты анализа направленного рассеяния излучения ( = 30) 1 – IR0, 2 – IR1, 3 – IRIR 0,0,0,00,511,L Рис.2. Результаты анализа направленного рассеяния излучения = 30, 1 – IR3, 2 – IR- индекс модельного рассеивающего центра (j = 0,1,означает присутствие воды в оболочке, j = 3,4 – ее отсутствие).

Перечень характеристик рассеивающего центра представлен в табл. 1, где d = ( - 0)/(1 - 0).

Таблица 1 – Характеристики рассеивающего центра j m() 0 1,1 m0 -(0,49-0,64i)d 2 (0,98-1,28i)/(1+d)+0,84+0,64i 3 m0 -(0,82-0,64i)d 4 (1,64-1,28i)/(1+d)+0,18+0,64i Как следует из рис. 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).

1,1,1,1,0 1 02 03 04 Приведенный размер Рис. 3. Значение G величины g, соответствующее максимальному уменьшению относительной индикатрисы рассеяния при наличии покрытия 1,0,0,0,0,0 Приведенный размер Рис. 4. Максимальное уменьшение индикатрисы рассеяния за счет неоднородности частицы Толщина слоя G IR(1)/IR(1,33) Таким образом, проблема интерпретации данных оптических измерений распределений частиц по размерам достаточно сложна.

С другой стороны, известен результат, когда в лабораторных условиях удалось радикально уменьшить влияние неоднородности частицы на результаты измерений, уменьшая угол рассеяния до 30о.

В табл. 2 приведены результаты модельных расчетов для характеристик рассеивающего центра, представленных в табл. 1, для частиц с приведенным радиусом ядра Ra = 30.

Pages:     || 2 |






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»