CC BY
114
Литература
1. Путилин Э.С. Оптические покрытия: Учебное пособие. - СПб: СПбГУ ИТМО, 2010. - 227 с.
2. Фурман Ш.А. Тонкослойные оптические покрытия. - Л.: Машиностроение, 1977. - 264 с.
3. Телен А. Конструирование многослойных интерференционных светофильтров // Физика тонких пленок. - М.: Мир, 1972. - Т. 5. - С. 46-83.
4. Виногоров В.М., Исаев А.Г., Симонов В.П., Учайкин А.Г. Оптический интерференционный спектро-делитель. Патент в02Б5/28, 2006 [Электронный ресурс]. - Режим доступа: http://bd.patent.su/2365000-2365999/pat/servl/servlet3200.html, св. Яз. рус. (дата обращения 01.12.2011).
Губанова Людмила Александровна - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор технических наук, профессор, [email protected]
Путилин Эдуард Степанович - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор технических наук, профессор, зав. кафедрой, eputilin @ yandex.ru
УДК 621.3.038.8.004.14
ВЫБОР ПАРАМЕТРОВ МНОГОВОЛНОВОГО АЭРОЗОЛЬНОГО ЛИДАРА ДЛЯ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ АТМОСФЕРЫ
Н.Н. Волков
Описывается возможность дистанционного зондирования атмосферы при помощи лидара, позволяющего получать коэффициенты обратного рассеяния и экстинкции аэрозоля в диапазоне длин волн от 355 до 1064 нм. Оцениваются параметры многоволнового лидара для мониторинга аэрозоля в высотном диапазоне 1-10 км. Ключевые слова: аэрозоль, лидар, рамановское рассеяние.
Введение
Увеличение количества антропогенных аэрозольных выбросов в атмосферу, повышение концентрации С02 и других парниковых газов оказывают воздействие на радиационный баланс планеты, и, как следствие, на ее климат. Следует отметить, что среди перечисленных факторов воздействие аэрозоля на радиационный баланс изучено в наименьшей степени. Воздействие атмосферного аэрозоля на радиационный баланс Земли двояко. Во-первых, аэрозоль влияет на радиационный баланс путем рассеяния и поглощения солнечной радиации (прямое воздействие). Во-вторых, это влияние происходит за счет модификации свойств облаков (косвенное воздействие). Для уменьшения соответствующих погрешностей в климатических моделях необходима глобальная информация об основных микрофизических параметрах аэрозоля, таких как радиус частиц, концентрация и комплексный показатель преломления.
Одним из наиболее перспективных инструментов, способных дистанционно проводить измерения параметров частиц, являются лидары. Малая длина волны излучения, сопоставимая с размером аэрозолей, и высокое пространственное разрешение делают лидары уникальным инструментом для исследования аэрозолей и облаков. Возможности метода значительно расширяются при проведении зондирования одновременно на нескольких длинах волн. Многоволновые системы позволяют получать информацию не только о спектральной зависимости коэффициентов рассеяния аэрозоля, но и о микрофизических параметрах частиц, таких как размер, концентрация и комплексный показатель преломления. Целью работы является оценка параметров многоволнового лидара, предназначенного для исследования аэрозолей различного типа: индустриального, пустынного и вулканического происхождения.
Существующие методы и средства для измерения параметров аэрозольных частиц атмосферы
Теоретические и экспериментальные исследования последних лет показывают, что измеренные коэффициенты обратного рассеяния и экстинкции аэрозоля могут быть инвертированы в микрофизические параметры частиц путем решения обратной задачи зондирования [1]. Наиболее удобными с практической точки зрения являются лидары на основе №:УАв-лазера с генератором третьей гармоники. Такая система позволяет измерять три коэффициента обратного рассеяния (Р) на длинах волн 355, 532, 1064 нм и два коэффициента экстинкции (а) на длинах волн 355, 532 нм. Вычисление коэффициентов экстинкции может быть проведено только с использованием рамановского рассеяния на молекулах азота или кислорода [1]. В измерениях, как правило, используется рассеяние на азоте, и длины волн соответствующих рамановских компонент для второй и третьей гармоники №:УАв составляют 608 нм и 387 нм. Излучаемая длина волны составляет 355 нм, в то время как после рассеяния длина волны составляет 387 нм, поэтому соответствующая экстинкция пересчитывается на 355 нм с использованием коэффициента Ангст-
Н.Н. Волков
рема. Соответствующая погрешность вычисления коэффициента экстинкции, связанная с неопределенностью выбора параметра Ангстрема, не велика и не превосходит 5%.
Сечение рамановского рассеяния почти на три порядка меньше сечения рэлеевского рассеяния, поэтому для надежной регистрации рамановского сигнала в нижней тропосфере необходимо использовать мощные лазерные источники излучения и приемные телескопы значительной апертуры. В работе проводится оценка параметров лазера и системы регистрации, требуемых для создания многоволнового аэрозольного лидара.
Процесс детектирования рассеянного излучения описывается статистикой Пуассона, и погрешность измерения может быть оценена как е = ^ , где N - количество зарегистрированных фотонов [2].
л/ЛТ
При измерении лидарного сигнала системой эта ошибка должна составлять не более 0,5%. Таким образом, из одного высотного слоя система должна регистрировать не менее 4 104 фотонов.
Количество детектируемых фотонов может быть оценено из лидарного уравнения [2]:
1Г Л, ■ 5■А■Т■ "П-В -2Г а^л
N-2-— ■ е 1о , (*)
г
где N - количество фотонов, принимаемых системой; Л, - количество фотонов, излучаемых в атмосферу; 5 - площадь приемного зеркала телескопа; г - дистанция зондирования; Т - коэффициент пропускания оптической системы; п - квантовая эффективность фотоприемника. Коэффициент обратного рассеяния молекул определяется как р = ст ■ нто1, где с - дифференциальное сечение рассеяния (рэлеевского либо рамановского), а пто1 - концентрация молекул. Коэффициент молекулярной экстинкции а связан с коэффициентом обратного рассеяния как а = ■8П ■ р .
Задача, описываемая уравнением (1), является неопределенной, поскольку это уравнение содержит произведение неизвестных Л, ■ 5 ■ Аг ■ Т ■ ^, а требуемое количество фотонов N может быть обеспечено как за счет увеличения энергии лазера, так и за счет увеличения апертуры телескопа, либо же за счет увеличения толщины слоя Аг.
Выбор и расчет параметров многоволнового аэрозольного лидара
При выборе параметров системы должны приниматься во внимание такие факторы, как стоимость, габариты, энергопотребление и т.д. Расчетные лидарные сигналы, обусловленные молекулярным рассеянием длин волн 355, 532, 1064 нм, показаны на рисунке. Результаты расчетов количества зарегистрированных фотонов за один лазерный импульс и соответствующие частоты счета фотонов на высотах 1 км и
N
10 км приведены в таблице. Частота счета V связана с количеством фотонов N как V = ^^ • В расчетах
с
апертура телескопа принималась равной 400 мм. Телескоп данной апертуры компактен, обеспечивает надежную регистрацию рамановского сигнала и не приводит к необоснованному удорожанию системы. Высотное разрешение Аг =100 м позволяет профилировать аэрозольный состав атмосферы и является достаточным для большинства климатологических задач [3]. Максимальная высота зондирования предполагалась равной 10 км. Энергия излучения составляла 100 мДж на каждой длине волны. Результаты расчетов могут быть легко пересчитаны для любых других параметров системы.
Длина волны приемного канала X, нм Тип рассеяния Скорость счета фотонов, принимаемых системой V, МГц Количество фотонов, принимаемых системой N
1 км 10 км 1 км 10 км
355 Рэлеевское 1,62106 3,05 • 103 10,81 • 105 2,04103
532 0,50106 1,71 • 103 3,35-105 1,41 • 103
1064 0,062 106 0,24103 0,41105 0,16103
387 Рамановское на N 1,61103 3,38 10,74102 2,26
608 0,43 • 103 1,51 2,87^102 1,01
Таблица. Параметры многоволновой аэрозольной лидарной системы
При выборе параметров системы следует рассматривать сигнал рамановского рассеяния как наиболее слабый. Из таблицы можно заключить, что на длине волны 608 нм регистрируется один фотон за импульс с высоты 10 км, т.е. для обеспечения требуемой точности измерения должны проводиться в течение 30 мин (при частоте повторения импульсов лазера 20 Гц). Такое временное разрешение достаточно для отслеживания временной эволюции вертикального распределения аэрозоля.
Помимо количества фотонов, определяемых по формуле (*), при выборе параметров системы необходимо учитывать динамический диапазон изменения лидарного сигнала (который для современных систем регистрации составляет четыре-пять порядков), а также максимальную и минимальную частоту счета фотонов фотодетектором. Лучшие фотоэлектронные умножители (ФЭУ), позволяющие проводить одновременное детектирование в аналоговом режиме и в режиме счета фотонов, надежно работают в диапазоне частот счета 103-10-2 МГц. На высоте 10 км скорость счета составляет 1,5 МГц для 608 нм, что достаточно для проведения надежных измерений (собственный шум ФЭУ составляет порядка 1 кГц). С уменьшением высоты скорость счета растет и на высоте 1 км достигает 1,6 103 МГц для 387 нм. Следует отметить, что в расчетах не рассматривалось влияние атмосферного аэрозоля. При типичных величинах экстинкции в пограничном слое атмосферы это может приводить к двукратному уменьшению сигнала. Таким образом, скорости счета фотонов в рамановских каналах соответствуют характеристикам детектора. В то же время рэлеевский сигнал на три порядка превосходит рамановский, поэтому для проведения измерений в каналы упругого рассеяния должны быть введены ослабляющие фильтры с коэффициентами ослабления от 100 до 1000.
В расчетах не рассматривался геометрический фактор перекрытия, характеризующий высоту, на которой лазерный пучок полностью входит в угловое поле телескопа. Для телескопа апертурой 400 мм с угловым полем 0,25 мрад эта высота составляет около 400 м при расходимости лазерного пучка 0,15 мрад. Таким образом, система позволяет производить измерения на высотах от 400 м, однако для этого часть апертуры телескопа должна быть заблокирована, чтобы обеспечить линейный режим работы ФЭУ. На рисунке показаны зависимости скорости счета фотодетектора от высоты зондирования для
2000 4000 6000 8000 10 000
Высота, м
Рисунок. Зависимость частоты счета детектора от высоты зондирования для молекулярного рассеяния
на длинах волн 355 нм (1), 532 нм (2) и 1064 нм (3)
Заключение
Проблема определения параметров лидарной системы является многопараметрической, и выбор конкретной конфигурации определяется спецификой поставленной задачи. В работе была проведена оценка параметров многоволнового рамановского лидара для мониторинга аэрозоля в высотном диапазоне 1-10 км. Проведенные вычисления показывают, что для уверенного приема сигнала рамановского рассеяния азота достаточной является система с апертурой телескопа 400 мм при энергии лазера не менее 100 мДж на длинах волн 355 и 532 нм. Для обеспечения требуемой точности на высоте 10 км измерение должно проводиться в течение 30 мин (при частоте повторения импульсов лазера 20 Гц). На основании проведенных расчетов спроектирована лидарная система, которая должна начать работу в 2012 г. в Центральной аэрологической обсерватории.
А.В. Демин, М.И. Моисеева
Литература
1. Ansmann A., Muller D. Lidar and atmospheric aerosol particles. - Springer Series in Optical Sciences, 2005 [Электронный ресурс]. - Режим доступа: http://www.springerlink.com/content/j75258113x74r713/iulltext.pdf, св. Яз. англ. (дата обращения 26.09.2011).
2. Межерис Р. Лазерное дистанционное зондирование: Пер. с англ. - М: Мир, 1987. - 510 с.
3. Althausen D., Muller D., Ansmann A., Wandinger U., Hube H., Clauder E., Zorner S. Scanning 6-wavelength 11-channel aerosol lidar // J. Atmos. and Oceanic Technol. - 2000. - V. 17. - P. 1469-1482.
Волков Николай Николаевич - Московский государственный университет геодезии и картографии
(МИИГАиК), аспирант, [email protected]
УДК 001.891.57
ИНВАРИАНТНАЯ МОДЕЛЬ ДЛЯ ОЦЕНКИ КОЭФФИЦИЕНТА ПРОПУСКАНИЯ АТМОСФЕРЫ ПРИ МОНИТОРИНГЕ ОБЪЕКТОВ В ОПТИЧЕСКОМ ДИАПАЗОНЕ СПЕКТРА А.В. Демин, М.И. Моисеева
Рассматривается задача учета влияния атмосферы на пропускание излучения в среднем и дальнем диапазоне инфракрасного спектра с целью повышения достоверности результатов при мониторинге летательных аппаратов. Предложена аналитическая апостериорная модель пропускания излучения атмосферой. Показана общая структура такой модели. Приводятся результаты идентификации значения коэффициента пропускания атмосферы на основе экспериментальных данных.
Ключевые слова: аналитическая модель, идентификация, коэффициент пропускания атмосферы, мониторинг объектов, экспериментальные данные.
Введение
Проблема достоверности результатов при дистанционном зондировании среды в среднем и дальнем диапазоне инфракрасного (ИК) спектра излучения (ИКСид) является актуальной при решении задачи пеленгации высокоскоростных летательных аппаратов (ЛА). Зондирование среды представляет собой процесс определения состояния среды по измеренным параметрам с целью выявления в ней отклонений от некоего «эталона», характеризуемого совокупностью физико-химических и оптико-физических параметров. Отклонения среды от «эталона» возникают как при наличии в среде объектов, отличных от нее (например, летательный аппарат), так и вследствие техногенного вмешательства.
На достоверность определения пеленгационных характеристик ЛА при его мониторинге влияет атмосфера. Последнее особенно сказывается при пеленгации высокоскоростных ЛА, которая возможна только в ИКСид. Это объясняется их температурой, лежащей в пределах 270-1000 К.
Для учета состояния атмосферы, характеризуемого коэффициентом пропускания (входящим в уравнение пеленгации), необходимо наличие базы данных со сведениями об атмосфере при любых возможных условиях. Последнее практически невозможно на момент мониторинга без соответствующей аппаратуры. На рис. 1 изображена структурная схема процесса мониторинга объекта с помощью оптико-электронных систем. В этом случае общим уравнением для измеряемой яркости будет ^изм = /атм ◦ /объекта, где /атм и /объекта - функции, характеризующие яркость атмосферы и объекта мониторинга; (◦) - композиция функций. Учет атмосферной компоненты /атм особенно важен тогда, когда необходима высокая вероятность принятия верного решения.
Известны работы [1, 2], посвященные измерениям коэффициента пропускания атмосферы в различных спектральных диапазонах. Следует заметить, что имеющиеся практические результаты получены при конкретных условиях, которые в общем случае могут отличаться от текущих. В настоящей работе предлагается получение апостериорной аналитической модели для учета коэффициента пропускания атмосферы ха(Х), инвариантной относительно условий наблюдения.
Аналитическая модель для оценки коэффициента пропускания атмосферы
При построении аналитической модели для оценки коэффициента пропускания атмосферы в ИКСид необходимо получить расчетное соотношение для прохождения ИК излучения в атмосфере.
Известно, что аналитическая модель экспериментальных данных, полученных в результате исследования определенного физического явления или технического объекта, может быть представлена композицией функций вида
Z = ЗД, r2,..., rk) о F2(S,, S2,..., Sm) ,
