CC BY
98
УДК 551.510.532; 551.501.816
ПРИМЕНЕНИЕ ПАРАМЕТРИЧЕСКОГО ГЕНЕРАТОРА СВЕТА ДЛЯ ДИСТАНЦИОННОГО МОНИТОРИНГА АТМОСФЕРЫ В ДИАПАЗОНЕ СПЕКТРА 3-4 МКМ
Романовский Олег Анатольевич1,2,
roa@iao.ru
Садовников Сергей Александрович1,
sadsa@iao.ru
Суханов Александр Яковлевич1,3,
say@iao.ru
Харченко Ольга Викторовна1,
olya@iao.ru
Яковлев Семён Владимирович1,2,
ysv@iao.ru
1 Институт оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук, Россия, 634021, г. Томск, пл. Академика Зуева, 1.
2 Томский государственный университет, Россия, 634050, г. Томск, пр. Ленина, 36.
3 Томский государственный университет систем управления и радиоэлектроники, Россия, 634050, г. Томск, пр. Ленина, 40.
Актуальность работы обусловлена необходимостью развития методов и систем дистанционного зондирования атмосферы для обеспечения оперативного мониторинга состояния окружающей среды.
Цель работы: разработка методики планирования и проведения лазерного дистанционного зондирования малых газовых составляющих атмосферы посредством подхода, объединяющего преимущества метода дифференциального поглощения (МДП) и дифференциальной оптической абсорбционной спектроскопии (ДОАС), поиск информативных длин волн в среднем ИК-диапа-зоне, численное моделирование лидарного зондирования малых газовых составляющих атмосферы на выбранных длинах волн. Методы исследования: метод дифференциального поглощения, метод численного моделирования лидарных сигналов. Результаты. Исследованы возможности применения лазерной системы с параметрической генерацией света на основе нелинейного кристалла КТЮАзОц для лидарного зондирования атмосферы в диапазоне спектра 3-4 мкм. Разработана методика лидарных измерений газовых компонент атмосферы с помощью методов дифференциального поглощения и дифференциальной оптической абсорбционной спектроскопии. Проведена апробация разработанной методики для оценки возможностей лидарного зондирования малых газовых составляющих атмосферы. Представлены результаты моделирования лидарных измерений малых газовых составляющих атмосферы в диапазоне 3-4 мкм, показывающие возможность восстановления лидарного сигнала на вертикальных трассах до 5 км, на горизонтальных трассах до 10 км при зондировании метана, формальдегида, бромоводо-рода и хлороводорода, используя лазерное излучение с шириной аппаратной функции 2 см-'.
Вывод. Лазер с параметрической генерацией света на основе нелинейного кристалла КТЮАзОц с рабочим диапазоном 3-4 мкм является перспективным источником излучения для дистанционного зондирования рассматриваемых в работе малых газовых составляющих атмосферы на приземных тропосферных трассах с применением разработанной методики.
Ключевые слова:
Атмосфера, лидарное зондирование, МДП, ДОАС, газовые составляющие атмосферы, нелинейные кристаллы.
Развитие лазерных бесконтактных методов диагностики среды (атмосфера, океан) требует разработки лазерных источников излучения, способных генерировать излучение в определенном диапазоне длин волн с малым шагом перестройки по частоте. Современные лидары дифференциального поглощения, как правило, позволяют контролировать лишь несколько газовых компонент атмосферы. Используемые в их составе мощные источники излучения, такие как СО2, DF (ОТ) лазеры, имеют ограниченное число спектральных совпадений линий излучения и линий поглощения атмосферных
газов [1]. При этом одновременный спектрально-временной анализ многочастотного отклика атмосферы в принципе не возможен из-за необходимости перестройки частоты излучения. Этого недостатка лишены ИК фурье-газоанализаторы, однако, из-за отсутствия мощных широкополосных источников когерентного излучения в их составе используются маломощные тепловые источники излучения. Они позволяют проводить лишь интегральные измерения с длиной измерительных трасс до нескольких сот метров с использованием специальных отражателей.
Современные источники лазерного излучения для решения задач лидарного зондирования созданы на основе широкодиапазонных ИК молекулярных лазеров, а также параметрических преобразователей частоты на основе нелинейных кристаллов, позволяющих путем генерации обертонов, гармоник, суммарных и разностных частот излучения лазеров перекрыть спектральный диапазон от 2 до 18 мкм. Для перекрытия ближней и средней ИК-областей используется излучение параметрических генераторов света (ПГС) на нелинейных кристаллах [1, 2].
Наземные [3] и самолетные [4] измерения по методу дифференциального поглощения (МДП) использовались для получения вертикальных профилей озона в тропосфере и стратосфере [5, 6], тропосферного водяного пара [7], для картирования сильного загрязнения приземного слоя этаном и N0^ определения горизонтального распределения в выбросах дизельного двигателя [8]. Также МДП может быть использован в лидарном зондировании метеорологических параметров атмосферы [9-12].
Однако в стандартном МДП должно выполняться требование монохроматичности лазерного излучения. Широкий спектральный профиль лазерного излучения уменьшает эффективный коэффициент поглощения, что снижает чувствительность измерений. Наличие аэрозоля в больших количествах также ведёт к существенным ошибкам МДП. Наконец, проблема стандартных измерений МДП, использующих только две длины волны, означает, что наличие мешающего поглощения сторонними газами не учитывается, что приводит к ошибкам, вызванным неточностью априорного знания коэффициентов поглощения.
Этих недостатков позволяет избежать метод дифференциальной оптической абсорбционной спектроскопии (ДОАС) [13-16]. ДОАС позволяет проводить спектрально разрешенные измерения в широкой полосе, что обеспечивает идентификацию нескольких газов даже в случае перекрывания полос поглощения. Ещё одним преимуществом метода является независимость от аэрозольного и молекулярного рассеяния благодаря технике фильтрации высоких частот, используемой для получения спектров. Однако возможности метода ДОАС для получения вертикальных профилей ограничены, в настоящее время проведены только усреднённые по трассе измерения.
Перспективным путем решения проблемы является метод, сочетающий преимущества обоих методов - пространственное разрешение МДП и идентификацию газов ДОАС. Новый метод использует широкополосное излучение и спектрометр с ПЗС-детектором, что обеспечивает измерение сигналов обратно рассеянного излучения с разрешением как по высоте, так и по длине волны одновременно.
Целью работы является разработка методики лидарного зондирования малых газовых соста-
вляющих (МГС) атмосферы, совмещающей МДП и ДОАС, и ее апробация в численном эксперименте для оценки возможностей лидарного зондирования газового состава атмосферы с помощью созданной лазерной системы на основе ПГС в диапазоне спектра 3-4 мкм.
Методика лидарного зондирования МГС атмосферы
с помощью МДП и ДОАС
Основной принцип действия лидара ДП заключается в том, что часть обратно рассеянного лазерного излучения преобразуется молекулами и аэрозольными частицами при прохождении через атмосферу. Лазерный пучок ослабляется (в соответствии с законом Ламберта-Бера) молекулярным поглощением и молекулярным и аэрозольным рассеянием. Таким образом, сигнал, зарегистрированный детектором, является функцией этого ослабления, доли обратно рассеянных фотонов, длины пути и мощности лазерного излучения. Сигнал PR из рассеивающего слоя толщиной Az можно представить в виде
PR (z, 1) = P0(1) At nWO( z) Azp( z, 1)e(1)
где P0(1) - мощность лазерного излучения; AD -ширина аппаратной функции; 0(z) - область перекрытия между лазерным лучом и полем зрения приемника; /3(z,1) - объемный коэффициент обратно рассеянного излучения; r(1) - эффективность приемо-передающей системы; Az - пространственное разрешение по трассе зондирования; ф,А) - объемный коэффициент ослабления.
МДП заключается в определении концентрации атмосферного газа по разнице лидарных сигналов на длинах волн on- и off- с разным молекулярным поглощением. Концентрация исследуемого газа определяется следующим образом:
, , 1 , { Pr (z, К„ )PR (z + Az, Xoff )^
n( z ) =-ln I-—\, (2)
2AoabsAz \ Pr (z, 1ff )Pr (z + Az, 1 )\ ^
где Aaabs- сечение поглощения.
Это выражение справедливо только в случае, если мешающее поглощение другими газами мало и рассеивающие свойства аэрозольных частиц не изменяются в диапазонах A1 и Az. Следствием влияния мешающего поглощения и неоднородности аэрозольного слоя являются большие ошибки восстановления профилей газов.
Недостатки МДП обусловлены неточностью априорного знания коэффициентов поглощения на двух и более длинах волн. Метод ДОАС позволяет избежать их посредством проведения мониторинга пропускания атмосферы с использованием широкополосных источников излучения с полушириной на полувысоте от единиц до сотен нм. Молекулярное поглощение (шириной порядка нескольких нм) в ДОАС отделяется от рассеяния по закону Ламберта-Бера, которое слабо зависит от длины
волны и может быть удалено с помощью высокочастотного фильтра. Кроме того, возможно, путем подбора (например, методом наименьших квадратов) рассчитанного молекулярного поглощения к наблюдаемому, одновременное измерение нескольких поглощающих газов. Измерение поглощения той или иной газовой составляющей на большом количестве длин волн позволяет решить проблему применения априорной информации, вносящей искажения в результаты экспериментальных исследований содержания газов в атмосфере.
В случае проведения ДОАС-измерений в первую очередь необходимо найти отношение наблюдаемого спектра (Рш) к референтному (Рш), которое определяется из спектра лазерного источника (Р0), полученного одним и тем же детектором. В качестве РВЕТ обычно используется спектр атмосферы при известном содержании поглощающих газов, например пассивный метод спектроскопии зенитного неба. Рассчитав отношение (Р0ВЯ/РдЕ?) и проведя минимизацию выражения (3) путём подгонки переменных коэффициентов, можно вычислить значения концентраций газов, входящих в диапазон зондирования:
Г Ров8 (у., 7,.)^
OBS V j J IS I
l Pref (vj, z ) J
= exp I -2Az^ a(Vj, zk) + ]n(ß(Vj, zt)) J x
V k=0 /
(• v+Av/2 X J v-Av/2 f (V' Vj ) X
f . Я Na f ng (Zk)Kg (V, Zk) + ^
x exp
-Az Hl
k=0 g V
+ng (Zk+1)Kg (V. Zk +1>
JJ *
(3)
где а(ул) - коэффициент ослабления анализируемой среды; Ду],г1) - коэффициент обратного рассеяния, включающий молекулярное и аэрозольное рассеяние в направлении назад; Ду - разрешающая способность приемника; Ау,у) - аппаратная функция приемника на частоте V/, Ые - число учитываемых газов; Кг(у,гк) - коэффициент поглощения соответствующего газа; пе(гк) - концентрация газовой компоненты.
Отношение (Р0ВЯ/РДЕ?) не зависит от спектра излучения лазерного источника или спектральной зависимости принимающей оптики, спектрографа, приемника, что является важным преимуществом МДП-ДОАС методики. Измеряемая таким образом величина составляет разницу содержания поглощающих газов в референтном случае и в случае реальной атмосферы.
Незначительно изменяющееся рассеяние Ду) обычно удаляется с помощью фильтра высоких частот. Для этого используются различные методы, но самый простой - применение полинома второго или третьего порядка, который подгоняется к разнице в оптической толще, а затем удаляется, оста-
вляя так называемый дифференциальный спектр или дифференциальную оптическую толщу. Необходимо отметить, что именно этот процесс фильтрации позволяет МДП-ДОАС нивелировать проблемы ослабления неизвестным или некорректно учитываемым рассеянием (аэрозольным или Рэлеевским), присущим МДП.
Обычно дифференциальный спектр восстанавливается по нескольким сотням точек, при этом число параметров подгонки составляет не более шести. Таким образом, уравнение (3) становится переопределённым и успешно решается методом наименьших квадратов. В свою очередь, коэффициенты подгонки варьируются для получения наилучшего совпадения спектра подгонки и наблюдаемого спектра. В общем случае разрешение метода соответствует молекулярному поглощению, при котором выполняется закон Ламберта. В этом случае, если известны коэффициенты поглощения, можно определить концентрацию каждого газа. Важным преимуществом МДП-ДОАС по сравнению с МДП является возможность одновременного определения концентрации нескольких газов.
Разработанную методику планирования и проведения лидарных МДП-ДОАС измерений можно представить в виде трёх частей, приведённых на рис. 1: 1) поисковая часть (белые блоки); 2) расчётная (серые); 3) блок обработки экспериментальных результатов лидарного зондирования.
Известно, что выражение для минимально об-наружимых концентраций в экспериментах с использованием ДП при лидарном зондировании имеет следующий вид [16]:
N (z) =
In P(V0N,z) - In p(yOFF, z) 2AK(V)z :
(4)
где P(von,z) и P(voff,z) - сигналы, регистрируемые фотоприемником с трассы зондирования длиной z на частотах VON и VOFF,
расположенные в узком спектральном диапазоне линии поглощения исследуемого газа и в соседнем окне прозрачности, а AK(v) - дифференциальный коэффициент поглощения:
AK(V) = K(Von) - K(Voff ), (5)
где K(von) и K(voff) - коэффициенты поглощения в и вне линии поглощения.
Для большинства реальных устройств регистрации и цифровой обработки, применяемых в экспериментах, значение разности логарифмов в выражении (4) должно удовлетворять условию [17]:
ln P(Von , z) - ln P(Voff , z) * 0,02. (6)
С учётом выражения (6) значения минимальных интенсивностей линий поглощения можно определить путём моделирования и сравнения спектров пропускания атмосферы анализируемой газовой компоненты и мешающих составляющих с концентрациями соответствующими ПДК для трассы длиной 10 км [12].
СП
мдп _г
Входные данные
МДП-ДОАС
Л_
Поиск длин волн излучения по критериям: 1. Оптимальные интенсивности линий поглощения 2. Отсутствие мешающего поглощения
3. Наличие окна прозрачности атмосферы
4. Совпадение центров линий поглощения
и лазерного излучения
Отбор диапазонов для широкополосного зондирования по критериям: 1. Оптимальные интенсивности линий и полос поглощения . Мощность зондирующего лазерного излучения 3. Отношение сигнал/шум
Ш
Расчет коэффициентов поглощения и спектров пропускания атмосферы
Моделирование лазерного зондирования с применением
статистических распределений атмосферных характеристик
Моделирование систематических и случайных погрешностей восстановления газовых профилей
Требования к параметрам лидара
Приём лидарных сигналов
Обработка лидарных сигналов
Профили концентраций зондируемых газовых составляющих
Рис. 1. Схема методики планирования и проведения лидарных МДП-ДОАС измерений Fig. 1. Scheme of planning technique and DIAL-DOAS lidar measurements
В случае реализации МДП-ДОАС систем оптимальные интенсивности линий и полос поглощения выбираются исходя из характеристик приёмного устройства лидара. Минимальная интенсивность линий поглощения должна соответствовать пространственному разрешению фотоприёмного прибора, величине его суммарных шумов, а также концентрациям исследуемых газовых компонент. Данные критерии определяются спецификой метода, сущность которого заключается в извлечении информации о концентрациях газовых компонент атмосферы посредством нахождения отношения мощностей двух сигналов - опорного (или референтного) и сигнального.
В качестве входных данных выступают параметры трассы зондирования, тип трассы (горизонтальная, вертикальная или наклонная), характеристики источника излучения (ширина линии излучения, мощность) и приёмной аппаратуры (спектральное и пространственное разрешение детектируемых сигналов).
Поисковый блок основывается на использовании базы спектроскопических данных [18], метеомоделей [19], априорной информации о коэффициентах аэрозольного ослабления, обратного рассеяния [20] и молекулярного рассеяния [21].
Расчётная часть алгоритма необходима для решения следующих задач:
• расчет коэффициентов поглощения и спектров пропускания атмосферы методом «line-by-line»;
• моделирование лазерного зондирования с применением статистических моделей распределения атмосферных характеристик;
• моделирование систематических и случайных погрешностей восстановления газовых профилей;
• выработка требований к параметрам лидара. Результаты моделирования и поиска информативных пар длин волн для МДП и диапазонов длин волн для МДП-ДОАС методик формируются в табличной и графической формах.
Методика применима для реализации систем как с традиционным функционалом МДП и ДОАС подходов по отдельности, так и для перспективного объединённого МДП-ДОАС подхода с применением одной и той же аппаратуры.
Лазерная система с ПГС для дистанционного зондирования атмосферы
В работе рассматривается лазерная система, входящая в состав лидара дифференциального поглощения и обеспечивающая перестраиваемую генерацию наносекундных импульсов излучения в спектральном диапазоне 3-4 мкм.
Лазерная система включает:
• Nd: YAG импульсный лазер LQ529B;
• преобразователь излучения с перестройкой длины волны 3-4 мкм;
• управление длиной волны от шаговых двигателей (ШД);
• контроллер для управления ШД;
• спектрометр S100;
• общее основание для лазера и преобразователя с системой заводки излучения накачки в преобразователь.
В табл. 1 и 2 представлены основные характеристики лазера накачки и преобразователя излучения.
Таблица 1. Спецификация лазера накачки LQ529B Table 1. Specification of LQ529B pumping laser
Частота следования импульсов Pulse frequency 10 Гц (Hz)
Выходная энергия: на 1064 нм Output energy: at 1064 nm 350 мДж (mJ)
Длительность импульса на 1064 нм, FWHM Pulse length at 1064 nm, FWHM 10-13 нс (ns)
Диаметр луча на 1064 нм Beam diameter at 1064 nm <6 мм (mm)
Расходимость на 1064 нм Divergence angle at 1064 nm ~1,5 мрад (mrad)
Стабильность энергии импульсов на 1064 нм Stability of pulse energy at 1064 nm ±2,5 %
Таблица 2. Спецификация преобразователя Table 2. Specification of radiation converter
Диапазон перестройки длины волны Wavelength tuning range 3-4 мкм (,um)
Ширина линии излучения Radiation line width 1-5 см-1 (cm-1)
Энергия в импульсе, в максимуме перестроечной кривой Pulse energy, in the tuning curve peak >6 мДж (mJ)
Частота следования импульсов Pulse frequency 10 Гц (Hz)
Расходимость излучения Radiation divergence angle <2 мрад (mrad)
Управление перестройкой длины волны Wavelength tuning control от 3 ШД (step motor)
На рис. 2 представлена перестроечная кривая параметрического генератора света на кристалле кТЮАя04 (КТА). Видно, что в диапазоне длин волн 3-4 мкм ПГС на кристалле КТА имеет довольно высокую энергию в импульсе на выходе системы (в максимуме перестроечной кривой её значения достигают величины > 6 мДж).
Разработанная МДП-ДОАС методика измерений МГС была апробирована для оценки уровня лидар-ных сигналов с использованием характеристик описанной выше лазерной системы с ПГС на основе нелинейного кристалла КТА. На основе характеристик лазера проведена оценка возможностей зондирования некоторых газовых примесей атмосферы на приземных трассах в исследуемой области спектра.
7 -|—■——г—-1——р--I—-1—-i——г—-1—-1——р—-1——I——р——г
J
Ч 4
= з
s =
lT>
II
■si
ОС О* Ä — MC.^irj'Är^OCO^Ä'-f'in гчгч
Длина волны КТА-ОРО, нм
Рис. 2. Перестроечная кривая ПГС на кристалле KTA в спектральном диапазоне 3-4 мкм
Fig. 2. Tuning curve of the KTA crystal-based optical parametric oscillator (OPO) in spectral range 3-4 ¡im»
Результаты моделирования лидарных измерений МГС атмоферы в диапазоне 3-4 мкм
Проведено численное моделирование зондирования некоторых газовых составляющих в атмосфере на приземных тропосферных трассах (на высотах до 5 км и при длинах горизонтальных трасс до 10 км). При моделировании использовалась стандартная модель атмосферы лета средних широт [19]. Учитывалось мешающее поглощение всех основных атмосферных газовых компонент. Входные данные для проведения численного моделирования представлены в табл. 3.
Таблица 3. Входные данные для проведения численного моделирования лазерного зондирования Table 3. Input data for numerical simulation of laser sounding
Параметр лидарной системы Lidar system parameters Значение параметра Parameter value
Площадь приемника А„риеи„.Ш=0,3 м) Receiver area Arec.(D=0,3 m) 7-10-8км2 (km2)
Ширина аппаратной функции Instrumental function width 2 см-1 (cm-1)
Эффективность приемной системы Receiving system efficiency 0,3
Пространственное разрешение AR Spatial resolution AR 1 км (km)
Максимальная энергия в импульсе Pulse energy maximum 5 мДж (mJ)
Частота следования импульсов Pulse frequency 10 Гц (Hz)
Длительность импульса Pulse length 10 нс (ns)
Расходимость излучения Radiation divergence 2 мрад (mrad)
Диапазон перестройки лазера Tuning range of the laser 3-4 мкм (,um)
NEP фотодетектора Photodetector NEP 1 • 10-12 Вт (W)
1E-5 1E-6 1E-7 Е-
и
1E-8 ^
1E-9 1E-10
3,39 3,40 3,41 3,42 3,43 3,44 3,45 3,46 3,47
Длина волны, мкм
on-line
У! 1 N>
HBiv
7
Z^
— / \ / \
2
3
3,95 4,00 4,05
Длина волны, мкм
Длина волны, мкм
on-line
1Е-6
я
1E-7 РЦ
3,50 3,52 3,54 3,56
Длина волны, мкм
1E-5
С. 1E-6 1E-7 ^
1E-8 ^ 1E-9
1E-5 1E-6
1E-7 03
1E-8
1E-9
Рис. 3. Пространственно и спектрально разрешенныелидарные сигналы вертикального зондирования метана, формальдегида, бромоводорода и хлороводорода в области работы ПГС на основе кристалла KTA (ширина аппаратной функции 2 см-)
Fig. 3. Spatially and spectrally resolved lidar signals of methane, formaldehyde hydrogen bromide and hydrogen chloride vertical sounding in the region of KTA-based OPO operation (instrumental function width is 2 cm-)
3,39 3,40 3,41 3,42 3,43 3,44 3,45 3,46 3,47
Длина волны, мкм
3,45 3,46 3,47 3,48
Длина волны, мкм
on-line
HBr ,
111! У У IUI
✓
3,90 3,92 3,94 3,96 3,98 4,00 4,02 4,04 4,06 4,08
Длина волны, мкм
3,50 3,52 3,54 3,56
Длина волны, мкм
E-5 E-6
E-7 н
E-8 И E-9 ^ E-10 E-11
vV
Рис. 4. Пространственно и спектрально разрешенные лидарные сигналы горизонтального зондирования метана, формальдегида, бромоводорода и хлороводорода в области работы ПГС на основе кристалла KTA (ширина аппаратной функции 2 см
Fig. 4. Spatially and spectrally resolved lidar signals of methane, formaldehyde hydrogen bromide and hydrogen chloride horizontal sounding in the region of KTA-based OPO operation (instrumental function width is 2 cm--)
on-line
3,45
3,46
3,47
1E-5
1E-8
4
5
3,48
4,10
on-line
on-line
on-line
3,48
В табл. 4 приведены информативные длины волн, пригодные для зондирования исследуемых в работе малых газовых составляющих атмосферы методом МДП-ДОАС.
Таблица 4. Длины волн, пригодные для зондирования метана, формальдегида, бромоводорода и хлорово-дорода
Table 4. Informative wavelengths chosen for sounding methane, formaldehyde, hydrogen bromide and hydrogen chloride
Газ Gas Яга-л., мкм (в воздухе) ^absorption, Mm (in the air) v„om, см-' (в воздухе) ^absorption., см-' (in the air) Тгаза Tgas Тмеш. погл. Tint. absorption
CH4 3,46490 - on-line 3,46395 - off-line 2886,078 2886,877 0,86 0,90 0,91 0,93
3,42746 - on-line 3,42962 - off-line 2917,607 2915,774 0,60 0,98 0,72 0,98
3,41537 - on-line 3,41716 - off-line 2927,932 2926,406 0,46 0,93 0,61 0,94
3,39101 - on-line 3,38844 - off-line 2948,970 2951,210 0,32 0,97 0,46 0,97
H2CO 3,72324 - on-line 3,72100 - off-line 2687,449 2685,832 0,85 0,92 0,99 0,98
3,46827 - on-line 3,46696 - off-line 2883,273 2884,37' 0,61 0,78 0,98 0,96
3,45081 - on-line 3,44928 - off-line 2897,864 2884,37' 0,42 0,83 0,91 0,98
HBr 3,95955 - on-line 3,96269 - off-line 2525,539 2523,538 0,82 0,99 0,99 0,99
3,98730 - on-line 3,99005 - off-line 2507,962 2506,234 0,79 0,99 0,99 0,99
4,01634 - on-line 4,01279 - off-line 2489,829 2492,03' 0,79 0,99 0,99 0,99
4,04642 - on-line 4,05281 - off-line 2471,320 2467,423 0,81 0,98 0,97 0,98
4,07744 - on-line 4,07180 - off-line 2452,5'9 2455,9'6 0,85 0,99 0,97 0,99
HCl 3,48950 - on-line 3,49548 - off-line 2865,740 2860,837 0,64 0,98 0,89 0,98
3,51576 - on-line 3,51985 - off-line 2844,335 284',030 0,57 0,95 0,91 0,96
3,54318 - on-line 3,54770 - off-line 282',829 28'8,727 0,52 0,96 0,85 0,96
3,57191 - on-line 3,57690 - off-line 2799,622 2795,7'6 0,60 0,97 0,93 0,97
На рис. 3 и 4 приведены рассчитанные для вертикальной (рис. 3) и горизонтальной (рис. 4) тропосферной трассы пространственно и спектрально разрешенные лидарные сигналы для СН4 (а), Н2СО (б),
СПИСОК ЛИТЕРАТУРЫ
1. Васильев Б.И., Маннун У.М. ИК-лидары дифференциального поглощения для экологического мониторинга окружающей среды // Квантовая электроника. - 2006. Т. 36. - № 9. -С. 801-820.
2. Mid-IR DIAL for high-resolution mapping of explosive precursors / V. Mitev, S. Babichenko, J. Bennes, R. Borelli, A. Dolfi-Bouteyre, L. Fiorani, L. Hespel, T. Huet, A. Palucci, M. Pistilli, A. Puiu, O. Rebane, I. Sobolev // SPIE Remote Sensing. International Society for Optics and Photonics. - 2013. - V. 8894. -P. 88940S-88940S-13.
HBr (в) и HCl (г) в области информативных длин волн МГС атмосферы, представленных в табл. 4.
На рисунках видно, что во всем рассматриваемом диапазоне высот 0-5 км и длин горизонтальных трасс 0-10 км уровень лидарных сигналов превышает уровень эквивалентной мощности шума фотоприемника №Р=10-12Вт.
Результаты численного моделирования показывают возможность восстановления лидарного сигнала на вертикальных трассах до 5 км, на горизонтальных трассах до 10 км при зондировании излучением ПГС на кристалле KTA в диапазоне 3-4 мкм метана, формальдегида, бромоводорода и хлороводорода.
Заключение
Применение разработанной методики лидарно-го зондирования МГС атмосферы, совмещающей МДП и ДОАС, и ее апробация в численном эксперименте подтверждают перспективность использования выбранных информативных длин волн для лидарного зондирования газового состава атмосферы с помощью лазерной системы на основе ПГС в диапазоне спектра 3-4 мкм. Проведенное численное моделирование показывает, что лазер с ПГС на основе нелинейного кристалла KTA является перспективным источником излучения для дистанционного зондирования методом МДП-ДОАС рассматриваемых в работе МГС атмосферы на приземных тропосферных трассах. Причем конструкция лазерной системы предусматривает возможность сужения линии лазерного излучения в пределах 0,01-0,05 см-1. Возможность такой модернизации, наряду с малым шагом перестройки линии лазерного излучения и наличием в рассматриваемом спектральном диапазоне линий поглощения других газовых составляющих и, в частности, атмосферных загрязнителей, делает рассматриваемый лазерный источник уникальным прибором для создания лидара дифференциального поглощения наземного базирования.
Работа выполнена при поддержке Российского научного фонда (Соглашение № 14-27-00022 в части разработки методики планирования и проведения лидарных МДП-ДОАС измерений) и Российского фонда фундаментальных исследований (Грант № 16-45-700722в части проведения численного моделирования лидарных измерений МГС атмосферы).
3. Sunesson J.A., Apituley A., Swart D.P.J. Differential absorption lidar system for routine monitoring of tropospheric ozone // Applied Optics. - 1994. - V. 33. - № 30. - P. 7045-7058.
4. Browell E.V. Differential absorption lidar sensing of ozone // Proceeding of the IEEE. - 1989. - V. 77. - № 3. - P. 419-432.
5. Differential absorption lidar sensing of ozone / T.J. McGee, M. Gross, U.N. Singh, J.J. Butler, P.E. Kimvilakani // Optical Engineering. - 1995. - V. 34. - № 5. - P. 1421-1430.
6. Лидарное зондирование озона в верхней тропосфере - нижней стратосфере: методика и результаты измерений / В.Д. Бурлаков, С.И. Долгий, А.А. Невзоров, А.В. Невзоров, О.А. Рома-
новский, О.В. Харченко // Известия Томского политехнического университета. - 2015. - Т. 326. - № 9. - С. 124-132.
7. Airborne differential absorption lidar system for measurements of atmospheric water vapor and aerosols / N.S. Higdon, E.V. Browell, P. Ponsardin, B.E. Grossmann, C.F. Butler, T.H. Chyba, M. Neale Mayo, R.J. Allen, A.W. Heuser, W.B. Grant, S. Ismail, S.D. Mayor, A.F. Carter // Applied Optics. - 1994. - V. 33. - № 27. - P. 6422-6438.
8. Toriumi R., Tai H., Takechi N. Tunable solid-state blue laser differential absorption lidar system for NO2 monitoring // Optical Engineering. - 1996. - V. 35. - № 8. - P. 2371-2375.
9. Харченко О.В. Методика планирования и проведения лидар-ных измерений профилей метеорологических параметров атмосферы // Оптика атмосферы и океана. 2012. - Т. 25. -№06.- С. 523-528.
10. Результаты моделирования лидарных измерений профилей метеопараметров с помощью обертонного СО-лазера / Г.Г. Матвиенко, О.А. Романовский, О.В. Харченко, С.В. Яковлев // Оптика атмосферы и океана. - 2014. - Т. 27. - № 02. -С. 123-125.
11. Романовский О.А., Харченко О.В., Яковлев С.В. Применение многоволновых ИК-лазеров для лидарных и трассовых измерений метеорологических параметров атмосферы // Известия вузов. Физика. - 2014. - Т. 57. - № 10. - С. 74-80.
12. Лидарный спектроскопический газоанализ атмосферы / С.М. Бобровников, Г.Г. Матвиенко, О.А. Романовский, И.Б. Сериков, А.Я. Суханов. - Томск: Изд-во ИОА СО РАН, 2014. - 510 с.
13. Platt U., Perner D., Patz H.W. Simultaneous measurement of atmospheric CH2O, O3, and NO2 by differential optical absorption // J. Geophys. Res. - 1979. - V. 84. - Ibis. C10. - P. 6329-6335.
14. Platt U. Air monitoring by spectroscopic techniques // Chemical Analysis Series. - 1994. - V. 127. - P. 27-84.
15. Platt U., Stutz J. Differential optical absorption spectroscopy. -New-York, Berlin, Heidelberg: Springer-Verlag, 2008. - 593 p.
16. Fourier-transform lidar / M. Douard, R. Bacis, Р. Rambaldi, A. Ross, J.-P. Wolf, G. Fabre, R. Stringat // Optics Letters. -1995. - V. 20. - № 20. - P. 2140-2143.
17. Коллис Р.Т.Х., Рассел П.Б. Лидарные измерения аэрозольных частиц и газов посредством упругого рассеяния назад и дифференциальное поглощение // Лазерный контроль атмосферы / под ред. Е.Д. Хинкли. - М.: Мир, 1979. - С. 91-180.
18. The HITRAN2012 molecular spectroscopic database / L.S. Rothman, I.E. Gordon, Y. Babikov, A. Barbe, Benner Ch.D., P.F. Bernath, M. Birk, L. Bizzocchi, V. Boudon, L.R. Brown,
A. Campargue, K. Chance, E.A. Cohen, L.H. Coudert, V.M. Devi,
B.J. Drouin, A. Fayt, J.-M. Flaud, R.R. Gamache, J.J. Harrison, J.-M. Hartmann, C. Hill, J.T. Hodges, D. Jacquemart, A. Jolly, J. Lamouroux, R.J. Le Roy, G. Li, D.A. Long, O.M. Lyulin,
C.J. Mackie, S.T. Massie, S. Mikhailenko, H.S.P. Muller, O.V. Naumenko, A.V. Nikitin, J. Orphal, V. Perevalov, A. Per-rin, E.R. Polovtseva, C. Richard, M.A.H. Smith, E. Starikova, K. Sung, S. Tashkun, J. Tennyson, G.C. Toon, Vl.G. Tyuterev, G. Wagner // Journal of Quantitative Spectroscopy and Radiative Transfer. - 2013. - V. 130. - P. 4-50.
19. Зуев В.Е., Комаров В.С. Статистические модели температуры и газовых компонент атмосферы. - Л.: Гидрометеоиздат, 1986.- 264 с.
20. Креков Г.М., Рахимов Р.Ф. Оптико-локационная модель континентального аэрозоля. - Новосибирск: Наука, 1982. - 199 с.
21. Optical properties of atmosphere. Report AFCRL-71-0297 / R.A. McClatchey, R.W. Fenn, J.E.A. Selby, F.E. Volz, J.S. Ga-ring. - Bedford: Mass, 1971. - 86 p.
Поступила 10.06.2016 г.
Информация об авторах
Романовский О.А., доктор физико-математических наук, зав. научно-образовательным центром Института оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук; профессор каф. управления инновациями факультета инновационных технологий Национального исследовательского Томского государственного университета.
Садовников С.А., младший научный сотрудник научно-образовательного центра Института оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук.
Суханов А.Я., кандидат физико-математических наук, старший научный сотрудник лаборатории лидарных методов Института оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук; доцент каф. автоматизированных систем управления Томского государственного университета систем управления и радиоэлектроники.
Харченко О.В., кандидат физико-математических наук, старший научный сотрудник центра лазерного зондирования атмосферы Института оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии
Яковлев С.В., кандидат физико-математических наук, научный сотрудник лаборатории лидарных методов Института оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук; доцент кафедры управления инновациями факультета инновационных технологий Национального исследовательского Томского государственного университета.
UDC 551.510.532; 551.501.816
APPLICATION OF OPTICAL PARAMETRIC OSCILLATOR FOR REMOTE MONITORING OF ATMOSPHERE IN THE 3-4 ¡m SPECTRAL RANGE
Oleg A. Romanovskii1'2,
roa@iao.ru
Sergey A. Sadovnikov1,
sadsa@iao.ru
Alexander Ya. Sukhanov1,3,
say@iao.ru
Olga V. Kharchenko1,
olya@iao.ru
Semyon V. Yakovlev1,2,
ysv@iao.ru
1 V.E. Zuev Institute of Atmospheric Optics SB RAS,
1, Academician Zuev square, Tomsk, 634021, Russia.
2 National Research Tomsk State University, 36, Lenin Avenue, Tomsk, 634050, Russia.
3 Tomsk State University of Control Systems and Radioelectronics, 40, Lenin Avenue, Tomsk, 634050, Russia.
Relevance of the research is caused by the need of developing the atmosphere remote sensing methods and systems for real-time monitoring of the environment.
The main aim of the study is to develop the methodology of planning and carrying out laser remote sensing of atmospheric trace gases through the approach combining the advantages of differential absorption lidar (DIAL) and the differential optical absorption spectroscopy (DOAS), search for informative wavelengths in the middle infrared range, numerical simulation of lidar sensing of atmospheric trace gases at the selected wavelengths.
The methods used in the study: differential absorption lidar, differential optical absorption spectroscopy.
The results. The authors have studied the opportunities of applying laser system with parametric light generation based on a nonlinear crystal KTiOAsO4 for atmosphere lidar sensing in the spectral range of 3-4 ¡m. The technique of lidar measurements of atmospheric gaseous components by the differential absorption lidar and differential optical absorption spectroscopy was developed. The authors tested the developed method to assess the feasibility of lidar sensing of atmospheric trace gases. The paper introduces the results of atmospheric lidar measurements modeling in the range of 3-4 ¡m, which demonstrate the possibility of lidar signal restoring on vertical paths up to 5 km and on horizontal paths up to 10 km when probing methane, formaldehyde, hydrogen bromide and hydrogen chloride using laser beam with a width of the instrumental function of 2 cm-1.
Conclusion. Laser with parametric light generation based on nonlinear KTiOAsO4 crystal with the operating range of 3-4 ¡m is the advanced light source for remote sensing the atmospheric trace gases, considered in the work, by the DIAL-DOAS technique at the ground level tropospheric paths.
Key words:
Atmosphere, lidar sounding, DIAL, DOAS, atmospheric trace gas, nonlinear crystals.
The research was supported by the Russian Science Foundation (Agreement no. 14-27-00022 as related to the development of the planning technique and DIAL-DOAS lidar measurement) and the Russian Foundation of Fundamental Research (Grant no. 16-45-700722 as related to numerical simulation of lidar measurements of atmosphere trace gas).
REFERENCES
1. Vasilev B.I., Mannun O.M. IK lidary differentsialnogo po-glashcheniya dlya ekologicheskogo monitoring okruzhayushchey sredy [IR differential-absorption lidars for ecological monitoring of the environment]. Quantum Electronics, 2006, vol. 36, no. 9, pp.801-820.
2. Mitev V., Babichenko S., Bennes J., Borelli R., Dolfi-Bouteyre A., Fiorani L., Hespel L., Huet T., Palucci A., Pistilli M., Puiu A., Rebane O., Sobolev I. Mid-IR DIAL for high-resolution mapping of explosive precursors. SPIE Remote Sensing. International So-
ciety for Optics and Photonics, 2013, vol. 8894, pp. 88940S-88940S-13.
3. Sunesson J.A., Apituley A., Swart D.P.J. Differential absorption lidar system for routine monitoring of tropospheric ozone. Applied Optics, 1994, vol. 33, no. 30, pp. 7045-7058.
4. Browell E.V. Differential absorption lidar sensing of ozone. Proceeding of the IEEE, 1989, vol. 77, no. 3, pp. 419-432.
5. McGee T.J., Gross M., Singh U.N., Butler J.J., Kimvilakani P.E. Differential absorption lidar sensing of ozone. Optical Engineering, 1995, vol. 34, no. 5, pp. 1421-1430.
6. Burlakov V.D., Dolgii S.I., Nevzorov A.A., Nevzorov A.V., Romanovskii O.A., Kharchenko O.V. Lidar sensing of ozone in the upper troposphere - lower stratosphere: technique and results of measurements. Bulletin of the Tomsk Polytechnic University, 2014, vol. 325, no. 2, pp. 127-136. In Rus.
7. Higdon N.S., Browell E.V., Ponsardin P., Grossmann B.E., Butler C.F., Chyba T.H., Neale Mayo M., Allen R.J., Heuser A.W., Grant W.B., Ismail S., Mayor S.D., Carter A.F. Airborne differential absorption lidar system for measurements of atmospheric water vapor and aerosols. Applied Optics, 1994, vol. 33, no. 27, pp. 6422-6438.
8. Toriumi R., Tai H., Takechi N. Tunable solid-state blue laser differential absorption lidar system for NO2 monitoring. Optical Engineering, 1996, vol. 35, no. 8, pp. 2371-2375.
9. Kharchenko O.V. The method for planning and carrying out the lidar measurements of profiles of atmospheric meteorological parameters. Atmospheric and Oceanic Optics, 2012, vol. 25, no. 06, pp. 523-528. In Rus.
10. Matvienko G.G., Romanovskii O.A., Kharchenko O.V., Yakovlev S.V. Simulation of lidar measurements of profiles of atmospheric meteorological parameters using overtone CO laser. Atmospheric and Oceanic Optics, 2014, vol. 27, no. 02, pp. 123-125. In Rus.
11. Romanovskii O.A., Kharchenko O.V., Yakovlev S.V. Application of multiwave IR laser for lidar and path measurements of the meteorological parameters of the atmosphere. Russian Physics Journal, 2015, vol. 57, no. 10, pp. 1380-1387. In Rus.
12. Bobrovnikov S.M., Matvienko G.G., Romanovskii O.A., Seri-kov I.B., Sukhanov A.Ya. Lidarny spektroskopichesky gazoanaliz atmosfery [Lidar spectroscopic gas analysis of the atmosphere]. Tomsk, IOA SO RAN Publ., 2014. 510 p.
13. Platt U., Perner D., Patz H. W. Simultaneous measurement of atmospheric CH2O, O3, and NO2 by differential optical absorption. J. Geophys. Res., 1979, vol. 84, Iss. C10, pp. 6329-6335.
14. Platt U. Air monitoring by spectroscopic techniques. Chemical Analysis Series, 1994, vol. 127, pp. 27-84.
15. Platt U., Stutz J. Differential optical absorption spectroscopy. New-York, Berlin, Heidelberg, Springer-Verlag, 2008. 593 p.
16. Douard M., Bacis R., Rambaldi P., Ross A., Wolf J.-P., Fabre G., Stringat R. Fourier-transform lidar. Optics Letters, 1995, vol. 20, no. 20, pp. 2140-2143.
17. Collis R.T.H., Russell P.B. Lidar measurement of particles and gases by elastic backscattering and differential absorption. Laser monitoring of the atmosphere. Berlin, Heidelberg, Springer-Verlag, 1976. pp. 71-151.
18. Rothman L.S., Gordon I.E., Babikov Y., Barbe A., Chris Benner D., Bernath P.F., Birk M., Bizzocchi L., Boudon V., Brown L.R., Campargue A., Chance K., Cohen E.A., Cou-dert L.H., Devi V.M., Drouin B.J., Fayt A., Flaud J.-M., Ga-mache R.R., Harrison J.J., Hartmann J.-M., Hill C., Hodges J.T., Jacquemart D., Jolly A., Lamouroux J., Le Roy R.J., Li G., Long D.A., Lyulin O.M., Mackie C.J., Massie S.T., Mikha-ilenko S., Müller H.S.P., Naumenko O.V., Nikitin A.V., Orphal J., Perevalov V., Perrin A., Polovtseva E.R., Richard C., Smith M.A.H., Starikova E., Sung K., Tashkun S., Tennyson J., Toon G.C., Tyuterev Vl.G., Wagner G. The HITRAN2012 molecular spectroscopic database. Journal of Quantitative Spectroscopy and Radiative Transfer, 2013, vol. 130, pp. 4-50.
19. Zuev V.E., Komarov V.S. Statisticheskie modeli temperatury i ga-zovykh komponent atmosfery [Statistic models of temperature and gas components of the atmosphere]. Leningrad, Gidrometeo-izdat Publ., 1986. 264 p.
20. Krekov G.M., Rakhimov R.F. Optiko-lokatsionnaya model kon-tinentalnogo aerozolya [Optical location model of continental aerosol]. Novosibirsk, Nauka Publ., 1982. 199 p.
21. McClatchey R.A., Fenn R.W., Selby J.E.A., Volz F.E., Garing J.S. Optical properties of atmosphere. Report AFCRL-71-0297. Bedford, Mass, 1971. 86 p.
Received: 10 May 2016.
Information about the authors
Oleg A. Romanovskii, Dr. Sc., head of the center, V.E. Zuev Institute of Atmospheric Optics SB RAS; professor, National Research Tomsk State University.
Sergey A. Sadovnikov, junior researcher, V.E. Zuev Institute of Atmospheric Optics SB RAS.
Alexander Ya. Sukhanov, Cand. Sc., senior researcher, V.E. Zuev Institute of Atmospheric Optics SB RAS; associate professor, Tomsk State University of Control Systems and Radioelectronics.
Olga V. Kharchenko, Cand. Sc., senior researcher, V.E. Zuev Institute of Atmospheric Optics SB RAS.
Semyon V. Yakovlev, Cand. Sc., researcher, V.E. Zuev Institute of Atmospheric Optics SB RAS; associate professor, National Research Tomsk State University.
