Спектральная фильтрация интерферометром Фабри-Перо лидарных откликов вращательного комбинационного рассеяния света тема автореферата и диссертации по физике, 01.04.05 ВАК РФ

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

Сериков Илья Борисович

СПЕКТРАЛЬНАЯ ФИЛЬТРАЦИЯ ИНТЕРФЕРОМЕТРОМ ФАБРИ-ПЕРО ЛИДАРНЫХ ОТКЛИКОВ ВРАЩАТЕЛЬНОГО КОМБИНАЦИОННОГО РАССЕЯНИЯ СВЕТА

Специальность 01.04.05 - оптика

АВТОРЕФЕРАТ

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

Томск - 2005

Работа выполнена в Институте оптики атмосферы СО РАН и Томском государственном университете

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

старший научный сотрудник, Юрий Филимонович Аршинов

Официальные оппоненты: доктор технических наук,

профессор Артур Александрович Мицель

Ведущая организация: Институт мониторинга климатических и экологических систем СО РАН

Защита состоится "14" апреля 2005 г. в 14 час. 30 мин. на заседании диссертационного совета Д 212.267.04 при Томском государственном университете (634050, г. Томск, пр. Ленина, 36), ауд. 119.

С диссертацией можно ознакомиться в Научной библиотеке Томского государственного университета.

кандидат физико-математических наук, Алексей Анатольевич Земляное

Автореферат разослан « » марта 2005

г.

Ученый секретарь диссертационного совета

/Б.Н.Пойзнер/

щ-ч $.тт

ти

Актуальность темы:

Температура, характеризующая среднюю для ансамбля молекул скорость их хаотического теплового движения, является необходимым исходным параметром в задачах прогнозирования состояния атмосферы и тенденций развития протекающих в ней процессов. Это задачи метеопрогноза, прогнозирования смогообразования, задачи оценки интенсивности переноса в атмосфере антропогенных примесей, задачи исследования динамики развития атмосферных процессов, к примеру, динамики формирования облачных образований, динамики фазовых переходов. Во всех этих, и многих других задачах подобного плана вертикальный профиль температуры атмосферы, и его изменение во времени, является важным входным параметром. Необходимые для получения этой информации круглосуточные измерения температуры атмосферы до сих пор ведутся с применением контактных методов, которые не обеспечивают необходимого временного разрешения и повторяемости пространственной привязки профиля измеряемых параметров. Альтернативой контактным методам являются лидарные методы измерений, не имеющие обозначенных недостатков. Анализ существующих на сегодняшний день лвдарных методов измерения температуры позволяет выделить, как наиболее перспективный, метод измерения температуры атмосферы по вращательным спектрам комбинационного рассеяния (КР) света молекулами атмосферного азота и кислорода. Экспериментальные исследования, проводившиеся различными научными группами, показали высокую эффективность данного метода при проведении измерений в отсутствии помехи рассеянного фона дневного неба (ночных измерений) [1-4, 3*] А. В случае дневных измерений наличие помехи фона дневного неба существенно ограничивает высотный диапазон применимости метода. Один из подходов, позволяющих если не устранить совершенно, то, по крайней мере, значительно ослабить это ограничение, построен на фильтрации помехи рассеянного фона неба интерферометром Фабри-Перо [5]. Основная идея подхода такова: поскольку вращательный спектр КР света азотом (и кислородом) имеет вид гребенки линий, расположенных в шкале частот практически эквидистантно [6, 7], для выделения линяй и подавления фона между ними может быть использован интерферометр Фабри-Перо (ИФП), представляющий собой частотный гребенчатый фильтр [8-11]. Эта идея фильтрации была высказана в работе [5] лишь в общем виде, ее физические основы не были изучены вплоть до момента проведения представляемой работы, не был построен соответствующий математический аппарат, не были выполнены экспериментальные исследования. Последнее обуславливает актуальность данной диссертационной работы, целью которой является исследование физических основ метода фильтрации интерферометром Фабри-Перо вращательного спектра комбинационного рассеяния света молекулами атмосферного азота (и кислорода) в задаче измерения температуры атмосферы. Для достижения поставленной цели

были сформулированы следующие задачи: .......

•'«М. НАЦИОНАЛЫ!*» 1

--I БИБЛИОТЕКА I

А символом * отмечены ссылки на список публикаций Серикове И.Б С.Ятр<ш А/,- (

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

2. Исследовать сравнительную эффективность измерений температуры атмосферы в зависимости от положения спектрального интервала, в котором возбуждаются спектры рассеяния.

3. Построить математическую модель интерферометра Фабри-Перо, позволяющую учесть физические аспекты фильтрации рассеянного фона дневного неба в задаче измерения температуры атмосферы лидаром по вращательным спектрам комбинационного рассеяния света.

Методы исследования:

В работе применен комплексный подход: используется теория рассеяния Рэлея [в задаче № 2], рассеяния Ми [в задаче № 2], теория спонтанного комбинационного рассеяния [в задачах № 1-3], теория многолучевой интерференции [в задаче № 3], физическое и численное моделирование, методы математической статистики, натурный эксперимент.

На защиту выносятся следующие положения:

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

2. Интерферометр Фабри-Перо, обеспечивающий максимальное снижение времени измерений температуры при отношении сигнал-фон, равном 0,1 , обеспечивает, при отношении сигнал-фон в интервале от 1 до 0,01 , снижение времени измерений не менее 87 % от предельного уровня, достижимого при выборе параметров интерферометра для каждого конкретного значения отношения сигнал-фон.

3. В лидарной схеме измерений наличие упорядоченных по апертуре дефектов поверхности зеркал интерферометра Фабри-Перо приводит к зависимости пропускания спектральных каналов лидара от дальности зондирования.

Достоверность и обоснованность результатов диссертационной работы:

На основе теории спонтанного КР выполнены оценки характеристик температурного канала лидара, найденные значения подтверждаются оценками, полученными из эксперимента. В частности, значение температурной чувствительности отношения вращательных откликов КР, полученное из расчета, с точностью в 3 % согласуется с величиной, полученной по результатам натурных измерений [1*]. Вертикальные профили температуры, восстановленные по результатам лидарных измерений, с хорошей точностью (1-2 градуса) ложатся на профили температуры, полученные при помощи метеорологических зондов, запуакавШйхеЯГпарй&йеЛвйо с лидарными измерениями [3*, 6*, 7*, 9*].

| г*.-'

* * **

Результаты расчета параметров полос пропускания температурного канала лидара, выполненного в соответствии с предложенным в работе критерием, совпадают с результатами подобных расчетов, представленных в работе [5], выполненных для случая измерений в отсутствии помех в соответствии с критерием минимизации ошибки измерений.

Математическая модель интерферометра Фабри-Перо построена на основе теории многолучевой интерференции. Оценка уровня подавления помехи рассеянного фона дневного неба интерферометром, полученная с применением принятой модели, с хорошей точностью (2 %) совпадает со значением оценки этого параметра, найденным из натурных измерений. Пропускание интерферометра для лидарных откликов в каналах вращательного КР, рассчитанное по модели, хорошо согласуется с данными, полученными в натурном эксперименте (отличие 5 %).

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

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

Научная новизна работы заключается в следующем:

Впервые исследованы физические основы метода фильтрации интерферометром Фабри-Перо помехи рассеянного фона дневного неба в спектральных промежутках между линиями вращательного спектра комбинационного рассеяния. Разработана методика расчета параметров интерферометра Фабри-Перо, предназначенного для фильтрации помехи рассеянного фона дневного неба в спектральных промежутках между линиями вращательного спектра КР. На основании расчета по данной методике впервые определены границы применимости метода фильтрации рассеянного фона дневного неба с использованием интерферометра Фабри-Перо, получена оценка эффективности метода. [Пункты 1-2 списка защищаемых положений].

Впервые описано влияние дефектов поверхности зеркал интерферометра, упорядоченных по апертуре, на зависимость пропускания спектральных каналов лидара от дальности зондирования. [Пункт 3 списка защищаемых положений].

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

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

Научная ценность результатов диссертационной работы:

Использование выводов и рекомендаций, изложенных в работе, при построении температурных лидаров сделает возможными дневные лидарные измерения высотных профилей температуры атмосферы с адекватным временным и пространственным разрешением. [Пункты 1-3 списка защищаемых положений].

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

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

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

Практическая значимость результатов диссертационной работы:

Предложенная в работе простая методика определения базы интерферометра может быть использована при настройке интерферометров Фабри-Перо, предназначенных для выделения линий вращательного спектра КР азота и кислорода.

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

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

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

Внедрение результатов работы:

Описанная в работе методика фильтрации помехи рассеянного фона дневного неба в промежутках между линиями вращательного спектра комбинационного рассеяния азота с применением интерферометра Фабри-Перо опробована (ноябрь 1999 года) и внедрена (апрель 2001 года) на экспериментальном полигоне Института тропосферных исследований (Лейпциг, Германия). С апреля 2001 года интерферометр, установленный в температурном канале лидара, находится в непрерывной эксплуатации.

Апробация работы:

Основные результаты диссертационной работы докладывались и обсуждались на 19-ой Международной конференции по лазерному и акустическому зондированию (1998 г., Аннаполис, США), 5-ом Международном симпозиуме по оптике атмосферы и океана (1998 г., Томск), 20-ой Международной конференции по лазерному и акустическому зондированию (2000 г., Виши, Франция), 5-ом Международном симпозиуме по тропосферному профай-лингу (2000 г., Аделаида, Австралия), 11-ой Международной школе по квантовой электронике (2000 г, Варна, Болгария), 21-ой Международной конференции по лазерному и акустическому зондированию (2002 г., Квебек, Канада), 6-ом Международном симпозиуме по тропосферному профайлингу (2003 г., Лейпциг, Германия), 22-ой Международной конференции по лазерному и акустическому зондированию (2004 г., Матера, Италия).

Структура и объем работы:

Диссертация состоит из ведения, четырех глав, заключения, списка литературы, включающего 91 наименование, и приложения, содержит 236 страниц, 153 рисунка и 5 таблиц.

Во введении обоснована актуальность темы, сформулирована цель работы, изложены новые научные результаты и положения, выносимые на защиту.

В первой обзорной главе дан анализ достоинств и недостатков методов дистанционного измерения температуры, известных на момент проведения исследования. Показано, что из числа существующих методов наиболее оптимальным для построения на его основе системы, способной обеспечить возможность проведения круглосуточных дистанционных измерений профилей температуры атмосферы с адекватным временным и пространственным разрешением, является лидарный метод измеригия температуры по вращательным спектрам КР света молекулами атмосферного азота и кислорода [14]. По сравнению с методами, построенными на других принципах, данный метод отличается тем, что для него: - не существует принципиальных ограничений высотного диапазона применимости;

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

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

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

- метод является самодостаточным с точки зрения калибровки системы.

Эффективность метода при проведении измерений в отсутствии мешающего влияния рассеянного фона дневного неба подтверждается работами многих исследовательских групп [1-4, 11*, 13*]. К примеру, результат, полученный нашей группой совместно с лидарной группой Института тропосферных исследований (Лейпциг, Германия): при полуторачасовом усреднении ли-дарных откликов отклонение рассчитанного по ним профиля температуры атмосферы от профиля температуры, полученного при помощи метеорологического шара-зонда, не превысило одного градуса до высоты в 14 км [3*].

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

В той же первой главе кратко охарактеризованы существующие варианты снижения уровня помехи рассеянного фона дневного неба в каналах вращательного КР температурного лидара - помехи, ограничивающей возможности методики при проведении дневных измерений. Такой вариант, как переход в "солнечно слепую" область спектра, при необходимой дальности зондирования не менее 15 километров, рассматривается как неприемлемый, поскольку одни и те же атмосферные составляющие, участвующие, с одной стороны, в поглощении солнечного излучения и обеспечивающие снижение уровня рассеянного фона, с другой стороны, участвуют в ослаблении зондирующего излучения [15]. Отличие в ослаблении сигнала и фона определяется отличием оптических толщ, а потому пропадает при их равенстве, делая переход в "солнечно слепую" область спектра резонным лишь для области высот порядка нескольких километров. В качестве метода снижения уровня помехи рассеянного фона дневного неба выбран метод, для которого не существует принципиальных ограничений по дальности зондирования, а именно, метод, построенный на фильтрации рассеянного фона дневного неба в промежутках между выделяемыми линиями вращательного спектра КР с применением интерферометра Фабри-Перо [5,4*].

Во второй главе работы рассматриваются вопросы оптимального (с точки зрения минимизации времени накопления, необходимого для обеспече-

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

К(Т)= * , (1)

Р3 п ДТ) - средняя частота следования фотоотсчетов в измерительном канале

лидара для линии вращательного спектра КР (азота или кислорода) с номером ] с учетом пропускания канала PJ .

Суммирование по ^ и ^ соответствует

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

Предложен критерий выбора положения и ширины выделяемых участков спектра, состоящий в минимизации вре-

I 1.0 £

щ о,а

* 0.6

I 0,4

Л

£ 0,2

Т = 220 К Т = 300 К

-200

20С

-100 0 100 О Частотный сдвиг, см1

Рис 1 Пример огибающих вращательного мени накопления, Дг(ДТ) выражение (2), спектра КР азота для двух температур газа

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

Д1(ДТ) =

г- - 2 Г 1

5>, "Г "Г ¿-и Л. Иы

1_ Ь J 1 Ли J

4,48 1

ДТ

сШ.'М

(2)

В соответствии с данным критерием рассчитаны положение и ширина участков вращательного спектра КР азота, выделение которых обеспечит минимальное время измерений (в рамках обсуждающейся методики измерений

температуры), в зависимости от уровня помехи несмещенного рассеяния, рис. 2, и помехи рассеянного фона дневного неба, рис. 3. 6

•"Г .....""! ' .....' (б)

.....|.....н -• - - -

-А-_-------

0 Ю'3 10"г 10"' ю° ю' юг Относительный уровень помехи

10° 10"г ю1 10° ю1 ю2 Относительный уровень помехи

Рис 2 Оптимальное положение границ выделяемых участков вращательного спектра КР азота с меньшим (а) и большим (б) частотным сдвигом относительно линии возбуждающего излучения в зависимости от относительного уровня помехи

Уровень помехи на рисунке 2 приведен в единицах отклика рассеяния для шестой линии стоксовой полосы вращательного спектра КР азота при температуре газа 220 К. При изменении уровня помехи в указанных на рисунке 2 пределах выбор положения и ширины выделяемых участков спектра с учетом уровня помехи позволяет снизить необходимое время накопления до двух раз по сравнению с вариантом, когда участки спектра выбираются для случая отсутствия помехи.

10~ 10" 10" 10'' 10'" 10° Спектральная плотность помехи

ю* ю-4 ю" ю" ю' ю" Спектральная плотность помехи

Рис 3 Оптимальное положение границ выделяемых участков вращательного спектра КР азота с меньшим (а) и большим (б) частотным сдвигом относительно линии возбуждающего излучения в зависимости от уровня спектральной плотности помехи

Спектральная плотность помехи на рисунке 3 приведена в единицах спектральной плотности отклика рассеяния для шестой линии стоксовой полосы вращательного спектра КР азота при температуре газа 220 К. При изменении уровня спектральной плотности помехи в указанных на рисунке 3 пределах, выбор положения и ширины выделяемых участков спектра с учетом уровня помехи позволяет снизить необходимое время накопления до трех раз по сравнению с тем вариантом, когда участки спектра выбираются для измерений в отсутствии помехи. Аналогичного рода зависимости (рис. 2, 3) получены в работе еще для двух вариантов формирования спектральных каналов:

выделения отдельно линий вращательного спектра КР кислорода и одновременного выделения линий спектра азота и спектра кислорода.

Для набора моделей, описывающих наиболее характерные атмосферные ситуации (чисто молекулярной атмосферы [16, 17], двух типов дымки, "прозрачной" и "плотной" (с дальностью видимости по горизонту в приземном слое 23 км и 5 км, для излучения с длиной волны 514,5 нм) [16, 18], атмосферы с присутствием кристаллических облаков [19, 20]), с учетом спектральной плотности освещенности прямой солнечной радиацией [13], выполнены оценки уровня рассеянного фона дневного неба, ожидаемого при проведении атмосферных измерений. С учетом полученных оценок исследована зависимость времени накопления (2), необходимого для обеспечения фиксированной точности измерений, от положения спектрального диапазона, в котором возбуждаются спектры рассеяния. Извлечение из этой части работы представлено на рисунке 4.

400

500 700 900

Длина волны, нм

1

0,7 300

500 700 900

Длина волны, нм

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

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

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

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

Для температурного канала комплексного Ми-КР-лидара Института тропосферных исследований (Лейпциг, Германия), на базе которого выполнены все представленные в работе экспериментальные исследования, сделаны оценки времени накопления, необходимого для обеспечения точности измерения температуры в 1 градус, в зависимости от дальности (высоты) зондирования, рис. 5. Параметры лидара приведены в таблице 1. На основании

I 5 10 15 20 25 30 Высота, км Рис 5 Время накопления как функция высоты измерений, при точности измерений температуры 1" Высотный ход температуры - модель среднеширотного лета [16] Пространственное разрешение 60 м. Зенитный угол Солнца 20° Кривые 1- молекулярная атмосфера (ночь), 2 - молекулярная атмосфера (день); 3 - кристаллическое облако, 4 -"прозрачная" дымка; 5 - "плотная" дымка

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

достаточной точности измерений Таблица! Параметры Ми-КР лидара Института (1-3 градуса) без дополнительных

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

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

Лазер

Длина волны, нм 532,075

Ширина линии по полувысоте, см"1 0,003

Энергия в импульсе, Дж 0,8

Частота повторения, Гц 30

Длительность импульса, не 10

Расходимость, мрад 0,5

Коллимация 1:10

Приемный телескоп

Диаметр, м 1

Эффективное фокусное расстояние, м 10

Угол поля зрения телескопа, мкрад 175

Фотодетектор

ФЭУ, Наташами Н7421-40

Квантовая эффективность, % 43

Длительность строба, не 400

Общее пропускание оптики, % 6

ветствии с моделью форма максимума функции пропускания интерферометра описывается соотношением

I 1-Я

агОД

•5{.тах

•ехр

872

2 а

я,а,в у

(3)

где - смещение по частоте от положения максимума пропускания, Я и а, соответственно, коэффициенты отражения и поглощения (по интенсивности) зеркал интерферометра, рц, Р, и р, - отражательная, угловая и дефектная

резкости функции пропускания интерферометра, соответственно

Ря =

1-Я Ду

г11*

ук59

Р =1^—1— 11 п 2 2л/21п2ст,

(4)

(5)

(6)

ду _ - область свободной дисперсии, й - база интерферометра; - положе-2д

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

- среднеквадратическое отклонение в модели распределения отклонений от среднего расстояния между зеркалами интерферометра (по апертуре зеркал)

' N

Р(8) =

1

л/2го

-ехр

гот.

2ст

(7)

л /

Резкость функции пропускания интерферометра с идеально плоскими зеркалами р„ 0 связана с отражательной и угловой резкостью соотношением

1 1 1

= + (8)

р2 Р

Величина аЯ(1в определяется как

1

а

Av

к.<1,е

2л/21п2

Фактор gf™"1

LR,e

'<" У

и фактор 6F

vFd у

характеризующие, соответственно,

снижение амплитуды и уширение максимума пропускания вследствие неплоскостности зеркал интерферометра, описываются аппроксимациями

gfn*x

1,450

0,298

0,760

1 +

R.9

Г

1 +

R,0

5F,

0,974

0,527

1 +

0,390

1 R.e

F< у 0,917

ч

1 +

R.0

ч у

N 2

1+

R.fl

1 + -

d У

1 + -

(10)

(11)

d У

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

С учетом доплеровского уширения линии вращательного спектра КР и гауссовой формы контура линии возбуждающего излучения пропускание интерферометра для отдельной линии имеет вид

С = J

1

, л/2тса

-ехр

1 (У-Ур -АУ,)'

2\

•fd.e(v)dy

(12)

интегрирование ведется в пределах ширины выделяемой спектральной линии 5vRaman' vo " частота> определяющая положение линии возбуждающего излучения. Ду( - частотный сдвиг линии с номером j вращательного спектра комбинационного рассеяния (выбранного типа молекул). CTyL

среднеквадратическое отклонение, характеризующее ширину линии вращательного спектра КР. fpp| (v) - функция пропускания интерферометра.

Уровень помехи рассеянного фона дневного неба в каналах КР, Bto и

Bta в выражении (2), определяется следующим образом

®io = РГ^АУ|0, В,; = P^'nf^AVh,, (13)

Ду|о, Avlo определяют ширину полосы пропускания соответствующего канала. п'ку - спектральная плотность средней частоты следования фотоотсчетов для фотонов рассеянного фона в каналах вращательного комбинационного рассеяния лидара. р™ - пропускание интерферометра для помехи рассеянного фона

sky

дневного неба, которое, в приближении спектральной однородности фона и постоянства величины квантовой эффективности фотоприемника, имеет вид

"M'-PiJir (14>

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

На основании анализа зависимости фактора снижения необходимого времени накопления при фильтрации фона интерферометром от величины отношения сигнал-фон определены границы применимости рассматриваемого метода фильтрации фона. Показано, что применение интерферометра оправдано в том случае, если отношение сигнал-фон лежит в пределах примерно от 1 до 0,01 , рис. 6. Так же показано, что для системы, предназначенной для проведения измерений температуры в условиях измене-

10"' 10'3 10'2 10"' 10° 10* о Отношение сигнал-фон

Рис 6 Фактор снижения времени накопления при абсолютном подавлении помехи рассеянного фона (пунктирная кривая) и фильтрации фона интерферометром (сплошная кривая) в зависимости от величины отношения сигнал-фон (в канале с большим частотным сдвигом полосы пропускания) ния величины отношения сигнал-фон в пределах от 1 до 0,01 , параметры интерферометра следует выбирать для отношения сигнал-фон, равного 0,1.

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

Оценка влияния качества поверх-

Ю"4 10 3 102 101 10° ю' о Отношение сигнал-фон

Рис 7 Фактор снижения времени накопления в зависимости от величины отношения сигнал-фон при фильтрации помехи рассеянного фона интерферометром Фаб-ри-Перо (сплошная кривая) и "идеальным"

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

кости порядка 100. Предел снижения времени накопления ограничен вследствие конечной ширины выделяемых спектральных линий. При качестве поверхности зеркал интерферометра, соответствующем дефектной резкости функции пропускания порядка 20, фильтрация фона интерферометром малоэффективна, рис. 8. Показано, что при изменении величины дефектной резкости функции пропускания интерферометра от 20 до 100 оптимальная величина отражательной резкости функции пропускания меняется в пределах примерно от 14 до 28, и, соответственно, область оптимальных значений коэффициента отражения зеркал интерферометра лежит в пределах от 0,8 до 0,9.

к

X

а>

8 со

X

X О) 2 1>

7 6 5 4

3 2 1 0

ю'2

Р,= 100 0 ЮО 26.4) Р„= 70 ((у=24,7) Р„= 50 (Рв= 22,4) Р„»40 (^=20,6) Р„=30 ^„=18,0) — 25 <Р„=1в,3) 20 (Рв- 14,3)

- !

-

Отношениесигнал-фон

10"

Рис 8 Фактор снижения необходимого времени накопления в зависимости от величины отношения сигнал-фон с учетом фильтрации интерферометром линий спектра кислорода Отражательная резкость функции пропускания при каждом конкретном значении дефектной резкости и заданной угловой резкости соответствует оптимуму для отношения сигнал-фон, равного одной десятой.

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

В работе рассмотрены различные схемы фильтрации фона. Показано, что использование двух интерферометров (одного для выделения линий вращательного спектра КР азота, второго - для выделения линий вращательного спектра КР кислорода) позволит снизить время накопления до полутора раз по сравнению с вариантом использования одного интерферометра и выделения только линий спектра азота (рис. 8). Также показано, что использование двух интерферометров (по описанной выше схеме) с относительно невысоким качеством поверхности зеркал (с дефектной резкостью порядка 30) позволяет добиться того же снижения времени накопления, что с одним интерферометром идеального качества (настроенным для выделения линий спектра азота).

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

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

Описаны схема и параметры температурного канала комплексного Ми-КР-лидара Института тропосферных исследований (Лейпциг, Германия), на базе которого проводилась апробация методики фильтрации рассеянного фона дневного с применением интерферометра Фабри-Перо. Опи-

Модель TL-38 (Burleigh Instruments, Inc )

Световой диаметр зеркал, мм 32

База интерферометра, мкм 628+-4,25

Резкость функции пропускания

при диаметре пучка 2 мм 22

при диаметре пучка 30 мм 16,5

Коэффициент отражения зеркал 0,885

Коэффициент поглощения зеркал 0

Эксперимент —■—"lo" —•—"hi"

сана конструкция и параметры использовавшегося интерферометра, таблица 2.

Предложена методика настройки интерферометра (определения базы) по уникальному совпадению в интерференционной картине, получающейся при освещении интерферометра рассеянным излучением ртутной лампы и гелий неонового лазера, диаметров интерференционных колец для излучения с длинами волн 632,817 нм (Не-Ые лазер), 546,073 нм (зеленая линия излучения ртутной лампы) и 576,960 нм (желтая линия излучения ртутной лампы). Описана методика автоподстройки базы интерферометра, обеспечивающая стабильность средней (в пределах апертуры зеркал интерферометра) базы на уровне 2 ангстрем, и отклонение зеркал от параллельности не более 0,05 микрорадиана. Предложенная схема автоподстройки позволяет осуществлять контроль и коррекцию настройки интерферометра параллельно с проведением измерений. Дано качественное описание программного модуля автоподстройки интерферометра.

Приведены результаты экспериментальной проверки оптимальности настройки, рис. 9. Измерено относительное пропускание интерферометра в нескольких порядках интерференции в окрестности рабочего порядка интерференции, в котором, в соответствии с расчетом, пропускание интерферометра должно быть максимально. Результаты

0.48

2360,5 2362,5 2364,5 2366,5 Порядок интерференции Рис 9 Расчетное и экспериментально измеренное пропускание интерферометра для канала с меньшим (1о) и большим (Ы) частотным сдвигом полосы пропускания для порядков интерференции в окрестности рабочего порядка 2364,5 (относительно Х = 532,075 нм)- Доверительные интервалы оценок рассчитаны по уровню 0,95

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

Проведены последовательные измерения без интерферометра и с интерферометром в температурном канале лидара для оценки уровня подавления интерферометром помехи рассеянного фона дневного неба, рис. 10. Для контроля уровня рассеянного фона использовался дополнительный спектральный канал - канал упругого рассеяния на частоте возбуждения откликов вращательного КР, использующихся для измерения температуры. Результаты измерений показали хорошее согласие (2 %) теории и эксперимента - оценка уровня подавления интерферометром рассеянного фона дневного неба по результатам измерений оказалась равной 16, в то время как по расчету эта величина, для используемйго интерферометра, должна быть равна 16,3.

0 4 8 12 16 20 24 28 32 36 40 Высота, км.

4 8 12 16 20 24 28 32 36 40 Высота, км.

Рис. 10. Лидарные отклики (число фотоотсчетов) в каналах вращательного КР (сплошные кривые) и канале упругого рассеяния (пунктирная кривые), без интерферометра (а) и с интерферометром (б) в температурном канале лидара Измерения без интерферометра 29 11 1999, 14 4515-26, с интерферометром 30 11 1999, 15 01-15 42, [4*]

С целью демонстрации принципиальных возможностей методики измерения температуры атмосферы по вращательным спектрам КР атмосферных азота и кислорода, и того качества данных, к которому должно быть приближено качество дневных измерений, приведен пример вертикального профиля температуры, полученного при помощи лидара в серии ночных измерений (без интерферометра) [3*], в сравнении с профилем температуры, полученным при помощи метеорологического зонда, рис. 11 (а). На рисунке 11 (б) приведен пример дневных измерений профиля температуры атмосферы (из серии круглосуточных лидарных измерений, выполненных с применением для фильтрации рассеянного фона дневного неба интерферометра Фабри-Перо [7*]), в сравнении с профилем температуры, полученным при помощи метеорологического зонда. Профили лидарной и зондовой температуры, представленные на рисунке 11 (б), хорошо согласуются друг с другом для высот от 2 до 10 км. Выше 10 км полезная составляющая сигнала существенно меньше помехи рас-

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

°. -20 | -зо

а -40 | -500) -60^ -70

20

10

оО 0

л -10

-20

го о. -30

1" -40

® -50

1- -60

6 8 10 12 Высота, км.

4 6В Высота, км.

10 12

Рис 11 Вертикальные профили температуры атмосферы по лидарным и зондовым измерениям (а) ночные измерения без интерферометра, выполнены 15 11 1999, 21 54-23 27 [3*] (б) дневные измерения с интерферометром, выполнены 23 04 2001, 14.55-15 55 [7*]. Вертикальными и горизонтальными отрезками на профилях "лидарной" температуры показаны, соответственно, доверительные интервалы по уровню 0,95 и пространственное разрешение

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

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

1 2 3 4 5 6 Расстояние, км Рис 12 Расчетов смещение восстанавливаемого по лидарным сигналам значения температуры от исходного, заданного значения (0°С), в зависимости от расстояния от лидара до рассеивающего объема.

щих падающий пучок, производит

вп

п г~

-

с

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

Е..............

Рис 13 Схема преобразования пучка с переменным углом расходимости в пучок с постоянным углом расходимости (а) и схема преобразования пучка с переменными угловыми и апертурными параметрами в пучок с постоянными угловыми и апертурными парамеграми (б) О диафрагма, [_ линза, Я] и Бц световоды Вертикальными пунктирными линиями показано положение фокальных плоскостей линзы

усреднение интенсивности по углу входного пучка, падающего на диафрагму D, в случае схемы рис. 13 (а), и падающего на входной торец световода 81, в случае схемы рис. 13 (б).

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

4 6 8 10 Высота, км. Рис. 14 Вертикальные профили температуры атмосферы, восстановленные по результатам измерений откликов рассеяния без скремблера (пунктирная кривая) и со скремблером (сплошная кривая) в температурном канапе лидара Измерения выполнены 30 марта 2003 г, без скремблера 30.03 2003, 21 01-21 31, со скремблером 30.03.2003, 22 47-23.52 [10*, 12*]

апертуре интерферометра дефектов поверхности зеркал.

В заключении приводятся основные результаты и выводы диссертационной работы:

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

2. Исследованы принципиальные возможности метода снижения необходимого времени накопления за счет фильтрации помехи рассеянного фона дневного неба в промежутках между линиями вращательных спектров комбинационного рассеяния молекул атмосферного воздуха с применением интерферометра Фабри-Перо. Построена математическая модель интерферометра, позволяющая сократить время вычислений в 100 раз при адекватном учете физических аспектов задачи фильтрации. Выработан критерий согласования параметров интерферометра Фабри-Перо, и в соответствии с ним определены границы применимости метода фильтрации рассеянного фона дневного неба в промежутках между линиями вращательного спектра комбинационного рассеяния атмосферных молекул с использованием интерферометра Фабри-Перо. Определены области значений параметров интерферометра и диапазон значений относительной спектральной плотности помехи, при которых эффективна фильтрация фона интерферометром. Показано, что:

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

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

- При измерении температуры атмосферы по вращательным спектрам комбинационного рассеяния атмосферного азота (или кислорода) использовать интерферометр Фабри-Перо для фильтрации помехи рассеянного фона дневного неба следует в том случае, если отношение сигнал-фон лежит в пределах диапазона от 1 до 0,01.

- Эффективность снижения времени накопления, отличающаяся не более чем на 9 % от той, что позволил бы получить интерферометр с идеально плоскими зеркалами, реализуется при дефектной резкости функции пропускания интерферометра больше 100.

- Фильтрация помехи рассеянного фона интерферометром при дефектной резкости функции пропускания интерферометра меньше 20 малоэффективна (при отношении сигнал-фон 0,01 эффективность снижения времени накопления в 2,4 раза меньше предельно возможной).

- При изменении величины дефектной резкости функции пропускания интерферометра от 20 до 100 оптимальная величина отражательной резкости

функции пропускания меняется в пределах примерно от 14 до 28, соответственно, область оптимальных значений коэффициента отражения зеркал интерферометра лежит в пределах от 0,8 до 0,9.

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

- При дефектной резкости функции пропускания интерферометра равной 100 и угловой резкости функции пропускания интерферометра равной 100 фильтрация интерферометром Фабри-Перо рассеянного фона дневного неба в промежутках меду линиями вращательного спектра комбинационного рассеяния азота позволит снизить время измерений температуры атмосферы в 6,6 раза при отношении сигнал-фон 0,01 , в 3,8 раза при отношении сигнал-фон 0,1 ив 1,34 раза при отношении сигнал-фон 1.

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

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

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

3. Проведены экспериментальные исследования, подтверждающие теоретические оценки эффективности фильтрации интерферометром Фабри-Перо помехи рассеянного фона дневного неба в каналах вращательного комбинационного рассеяния. Полученные экспериментальные данные показали, что применение интерферометра Фабри-Перо позволяет значительно снизить уровень помехи рассеянного фона дневного неба в каналах вращательного КР (с использовавшимся в эксперименте интерферометром, до 16 раз), а так же существенно улучшить отношение сигнал-фон (с использовавшимся интерферометром, до 6 раз).

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

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

Основные результаты опубликованы в следующих работах:

1*. AnsmannA., ArshinovYu, Bobrovnikov S., Mattis I, Serikovl, WandingerU. Double grating monochromator for a pure rotational Raman-lidar // Fifth International symposium on atmospheric and ocean optics / Zuev V.E., Matvienko G.G. (Editors), Proceedings of SPIE. -1998. - V.3583. - P. 491-497.

2* Wandinger U., Mattis I., Ansmann A., Arshinov Yu , Bobrovnikov S., Serikov I Tropospheric temperature profiling based on detection of Stokes and Anti-Stokes rotational Raman lines at 532 nm // 19th International laser radar conference- 6-10 July 1998 - Annapolis, 1998, P. 297-299.

3*. Althausen D., Arshinov Y., Bobrovnikov S , Serikovl., AnsmannA., Mattis I., Wandinger U. Temperature Profiling in the Troposphere Using a Pure Rotational Raman Lidar // Needs and Technology: Vth International Symposium on Tropospheric Profiling. - Adelaide, Australia, 2000. P. 31-33.

4* Arshinov Yu , Bobrovnikov S., Serikov I, Ansmann A., Althausen D., Mattis I., Wandinger U. Spectrally absolute instrumental approach to isolate pure rotational Raman lidar returns from nitrogen molecules of atmosphere // Advances in laser remote sensing / Editors- Dabas A , Loth C„ Pelon J. - Vichy, France, 2000, P. 121-124.

5*. ArshinovYu, Bobrovnikov S., Serikovl, Althausen D., Mattis I., WandingerU., Ansmann A, Temperature profiling in the atmosphere using lidars // 11th International school on quantum electronics- laser physics and applications / Atanasov P.A., Cartaleva S. (Editors), Proceedings of SPIE. - 2000 - V.4397. - P. 453-460.

6*. Mattis I., Ansmann A , Althausen D., Jaenish V., Wandinger U., Mueller D , Arshinov Yu , Bobrovnikov S., Serikovl Raman lidar observations of relative humidity // 21st International laser radar conference: 8-12 July 2002. - Quebec, 2002, P. 27-31.

7*. Bobrovnikov S., Arshinov Yu., Serikov I, Althausen D., Ansmann A., Mattis I., Wandinger U. Daytime temperature profiling in the troposphere with a pure rotational Raman lidar // 21st International laser radar conference- 8-12 July 2002. - Quebec, 2002, P. 717-720.

8*. Serikov I, Arshinov Yu , Bobrovnikov S , Althausen D., Mattis I. Distortion of the temperature profile of the atmosphere acquired with a pure rotational Raman lidar due to sphericity of the Fabry-Perot interferometer plates // 21st International laser radar conference: 8-12 July 2002. - Quebec, 2002, P. 721-724.

9*. Mattis I., Ansmann A., Althausen D., Jaenish V., Wandinger U., Mueller D., Arshinov Yu.F., Bobrovnikov S.M., Serikov I.B. Relative humidity profiling in the troposphere with a Raman lidar // Applied Optics. - 2002. - V.41. - №30. - P. 6451-6462

10*. ArshinovYu., Althausen D., AnsmannA., Bobrovnikov S., Mattis I, Mueller D., Serikov I, Wandinger U Round-the-clock temperature profiling in the troposphere with Raman lidar // Sixth International symposium on tropospheric profiling: Needs and technologies: 14-20 September 2003. - Leipzig, Saxony, 2003, P. 142-144.

II*. Balin I., Simeonov V., Serikovl, Bobrovnikov S., Calpini В., ArshinovYu., van den Berg H Simultaneous measurement of temperature, water vapor, aerosol extinction and backscatter by Raman lidar // Reviewed and revised papers presented at 22-nd International laser radar conference- 12-16 July 2004 - Matera, Italy, 2004, P 139-142

12* ArshinovYu, Bobrovnikov S., Serikovl., Althausen D , Ansmann Л., Mattis I, Mueller D., Wandinger U. Optic-fiber scramblers and Fourier transform lens as a mean to tackle the problem on the overlap factor of lidar // Reviewed and revised papers presented at 22nd International laser radar conference. 12-16 July 2004. - Matera, Italy, 2004, P. 227-230.

13* Balin I., Serikovl, Bobrovnikov S , SimeonovV, Calpini В., Arshinov Y., van den Bergh H. Simultaneous measurement of atmospheric temperature, humidity, and aerosol extinction and backscatter coefficients by a combined vibrational-pure-rotational Raman lidar // Applied Physics B. - 2004. - № 79. P. 775-782.

Цитируемая литература:

1 Arshinov Yu.F., Bobrovnikov S.M., Zuev V.E., Mitev V.M. Atmospheric temperature measurements using a pure rotational Raman lidar4 Aplied Optics. - 1983 - V 22 -№ 19.-P. 2984-2990.

2 Nedeljkovich D . Hauchecorne A., Chanin M.L. Rotational Raman lidar to measure temperature from the groung to 30 km // ГЕЕЕ Trans. Geosci. Remote Sens -1993 -V.31.-№l.-P 90-101.

3. Vaughan G., Wareing D.P., Pepler S.J., Thomas L., Mitev V. Atmospheric temperature measurements made by rotational Raman scattering// Applied Optics. - 1993. - V 32, №¡15, P. 2758-2764.

4. Behrendt A., Reichardt J. Atmospheric temperature profiling in the presence of clouds with a pure rotational Raman lidar by use of an interference-filter-based polychromator // Applied Optics. - 2000. - V 39. - №9. - P. 1372-1378.

5 Аршинов Ю.Ф Разработка и исследование лидарного метода измерения температуры атмосферы по вращательным спектрам КР молекул воздуха- Диссертация на соискание ученой степени кандидата физико-математических наук. - Томск, 1982. - 139 с

6. Плачек Г Рэлеевское рассеяние и Раман-эффект - Харьков: ГНТИ Укр.ССР, 1935. - 173 с.

7. Butcher R.J., Willetts D.V., Jones W J. On the use of Fabry-Perot etalon for the determination of rotational constants of simple molecules - the pure rotational Raman spectra of oxygen and nitrogen //Proc. Roy. Soc. Lon A - 1971 - V. 324, P 231-245.

8. Малышев В И Введение в экспериментальную спектроскопию. - М Наука, 1979 - 478 с.

9. Тарасов К И. Спектральные приборы - Ленинград: Машиностроение, 1977, - 387 с

10. Нагибина И.М. Интерференция и дифракция света. - Ленинград' Машиностроение, 1974 - 359 с.

11. Зайдель А.Н. Основы спектрального анализа. - М.: Наука, 1965. - 322 с.

12. Monzon J.J., Sanchez-Soto L.L., Bernabeu Е. Influence of coating thickness on the performance of a Fabry-Perot interferometer // Applied Optics. - 1991. - V.30. - №28. - P 4126-4132.

13. Valley S.L. Handbook of Geophysics and Space Environments. - Belford, Massachusetts- Air Force Cambridge Research Laboratories, 1965. - 692 p

14. CooneyJ. Measurement of atmospheric temperature profiles by Raman backscatter // J. Appl. Meteorol. - 1972. - V.ll. -№1.-P. 108-112.

15. Аршинов Ю.Ф., Бобровников C.M., Попов А.Г., Шелефонтюк Д.И., Шумский В.К. Дальность обнаружения малых загрязняющих компонентов атмосферы в солнечнослепой области спектра // Оптика атмосферы и океана. - 1994. - т. 7. - № 8. - С. 1132-1138

16. McClatchey R.A , Fenn R.W , Selby J E.A., Volz F E , Garing J S Optical properties of the atmosphere // Environmental research papers - 1972. - №411 - p. 108.

17. Мак-Картни Э. Оптика Атмосферы. - M.: Мир, 1979. - 421 с.

18. Дейрменджан Д. Рассеяние электромагнитного излучения полидисперсньши частицами: Пер. с англ. -М.: Мир, 1971. - 165 с.

19. Measurements with lidar systems during the international cirrus experiment 1989 - Hamburg. Max-Planck-lnstitut fuer Meteorologie, №60, 1990 - 158 p

20 Волковицкий O.A., Павлова Л.Н , Петрушин А Г. Оптические свойства кристаллических облаков. - Ленинград: Гидрометеоиздат, 1984. - 198 с,

21. Kaufmann К. Photomultipliers detect rare events // Laser focus world. - 2000. -V.36. -№10. - P. 77-81.

Отпечатано на участке оперативной полиграфии Редакционно-издательского отдела I I У Лицензия ПД №00208 от 20 декабря 1999 г.

Заказ № 59 от" ± "

2005 г. Тираж 100 экз.

43 84

РНБ Русский фонд

2006-4 14711

Введение

Глава 1. Методы дистанционного измерения температуры

Глава 2. Спектральная фильтрация в лидарных измерениях температуры атмосферы по вращательным спектрам комбинационного рассеяния молекул атмосферного воздуха

2.1. Вращательный спектр комбинационного рассеяния азота и кислорода

2.2. Методика лидарных измерений температуры атмосферы по вращательным спектрам комбинационного рассеяния

2.3. Выбор положения и ширины спектральных участков

2.3.1. Критерий выбора участков спектра

2.3.2. Измерения в отсутствии помех

2.3.2.1. Выбор отдельных линий спектра

2.3.2.2. Выбор участков спектра

2.3.3. Измерения при конечном отношении сигнал-шум

2.3.3.1. Помеха несмещенного рассеяния

2.3.3.2. Помеха фонового излучения

Глава 3. Фильтрация помехи фона в каналах вращательного комбинационного рассеяния с применением интерферометра Фабри-Перо

3.1. Функция пропускания интерферометра Фабри-Перо

3.1.1. Пропускание интерферометра для пучка конечной расходимости

3.1.2. Зависимость пропускания интерферометра от угла падения пучка

3.1.3. Влияние дефектов поверхности зеркал

3.2. Согласование функции пропускания интерферометра Фабри-Перо и вращательного спектра комбинационного рассеяния

3.3. Подавление помехи рассеянного фона дневного неба

3.4. Фактор повышения отношения сигнал-фон

3.5. Подавление помехи несмещенного рассеяния

3.6. Выбор отражательной резкости функции пропускания

3.6.1. Критерий выбора отражательной резкости функции пропускания

3.6.2. Идеальный интерферометр в плоскопараллельном падающем пучке

3.6.3. Выбор отражательной резкости с учетом дефектов поверхности зеркал

3.6.4. Выбор отражательной резкости с учетом дефектов поверхности зеркал и конечной расходимости падающего пучка

3.7. Учет фильтрации интерферометром отклика рассеяния кислорода

3.8. Фильтрация фона в промежутках между линиями вращательных спектров комбинационного рассеяния азота и кислорода

3.9. Выбор ширины линии возбуждающего излучения

Глава 4. Практические вопросы реализации методики фильтрации помехи рассеянного фона дневного неба с применением интерферометра Фабри-Перо

4.1. Температурный канал комплексного Ми-КР-лидара

4.2. Настройка базы интерферометра

4.2.1. Методика определения базы интерферометра

4.2.2. Автоподстройка базы интерферометра

4.2.3. Экспериментальная проверка настройки

4.3. Фильтрация помехи рассеянного фона дневного неба с применением интерферометра Фабри-Перо

4.4. Влияние крупномасштабных дефектов поверхности зеркал интерферометра на точность измерений температуры

Температура, характеризующая среднюю для ансамбля молекул скорость их хаотического теплового движения, является необходимым исходным параметром в задачах прогнозирования состояния атмосферы и тенденций развития протекающих в ней процессов. Это задачи метеопрогноза, прогнозирования смогообразования, задачи оценки интенсивности переноса в атмосфере антропогенных примесей, задачи исследования динамики развития атмосферных процессов, к примеру, динамики формирования облачных образований, динамики фазовых переходов. Во всех этих, и многих других задачах подобного плана вертикальный профиль температуры атмосферы, и его изменение во времени, является важным входным параметром. Необходимые для получения этой информации круглосуточные измерения температуры атмосферы до сих пор ведутся с применением контактных методов, которые не обеспечивают необходимого временного разрешения и повторяемости пространственной привязки профиля измеряемых параметров. Альтернативой контактным методам являются лидарные методы измерений, не имеющие обозначенных недостатков. Анализ существующих на сегодняшний день лидарных методов измерения температуры позволяет выделить, как наиболее перспективный, метод измерения температуры атмосферы по вращательным спектрам комбинационного рассеяния света молекулами атмосферного азота и кислорода. Экспериментальные исследования, проводившиеся различными научными группами, показали высокую эффективность данного метода при проведении измерений в отсутствии помехи рассеянного фона дневного неба (ночных измерений) [1-5]. В случае дневных измерений наличие помехи фона дневного неба существенно ограничивает высотный диапазон применимости метода. Один из подходов, позволяющих если не устранить совершенно, то, по крайней мере, значительно ослабить это ограничение, построен на фильтрации помехи рассеянного фона неба интерферометром Фабри-Перо [6]. Основная идея подхода такова: поскольку вращательный спектр комбинационного рассеяния азота (и кислорода) имеет вид гребенки линий, расположенных в шкале частот практически эквидистантно [7, 8], для выделения линий и подавления фона между ними может быть использован интерферометр Фабри-Перо, представляющий собой гребенчатый фильтр [9-12]. Эта идея фильтрации была высказана в работе [6], но лишь в общем виде, ее физические основы не были изучены вплоть до момента проведения представляемой работы, не был построен соответствующий математический аппарат, не были выполнены экспериментальные исследования. Последнее обуславливает актуальность данной диссертационной работы, целью которой является исследование физических основ метода фильтрации интерферометром Фабри-Перо вращательного спектра комбинационного рассеяния света молекулами атмосферного азота (и кислорода) в задаче измерения температуры атмосферы. Для достижения поставленной цели были сформулированы следующие задачи:

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

2. Исследовать сравнительную эффективность измерений температуры атмосферы в зависимости от положения спектрального диапазона, в котором возбуждаются спектры рассеяния.

3. Построить математическую модель интерферометра Фабри-Перо, позволяющую учесть физические аспекты фильтрации рассеянного фона дневного неба в задаче измерения температуры атмосферы лидаром по вращательным спектрам комбинационного рассеяния света.

Методы исследования:

В работе применен комплексный подход: используется теория рассеяния Рэлея, рассеяния Ми, теория спонтанного комбинационного рассеяния, теория многолучевой интерференции, физическое и численное моделирование, методы математической статистики, натурный эксперимент.

Научная новизна работы заключается в следующем:

Впервые исследованы физические основы метода фильтрации интерферометром Фаб-ри-Перо помехи рассеянного фона дневного неба в спектральных промежутках между линиями вращательного спектра комбинационного рассеяния. Разработана методика расчета параметров интерферометра Фабри-Перо, предназначенного для фильтрации помехи рассеянного фона дневного неба в спектральных промежутках между линиями вращательного спектра комбинационного рассеяния. На основании расчета по данной методике впервые определены границы применимости метода фильтрации рассеянного фона дневного неба с использованием интерферометра Фабри-Перо, получена оценка эффективности метода.

Впервые описано влияние дефектов поверхности зеркал интерферометра, упорядоченных по апертуре, на зависимость пропускания спектральных каналов лидара от дальности зондирования.

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

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

Научная ценность результатов диссертационной работы:

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

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

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

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

Практическая значимость результатов диссертационной работы:

Предложенная в работе простая методика определения базы интерферометра может бьггь использована при настройке интерферометров Фабри-Перо, предназначенных для выделения линий вращательного спектра комбинационного рассеяния азота и кислорода.

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

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

Внедрение результатов работы:

Описанная в работе методика фильтрации помехи рассеянного фона дневного неба в промежутках между линиями вращательного спектра комбинационного рассеяния азота с применением интерферометра Фабри-Перо опробована (ноябрь 1999 года) и внедрена (апрель 2001 года) на экспериментальном полигоне Института тропосферных исследований (Лейпциг, Германия). С апреля 2001 года интерферометр, установленный в температурном канале лидара, находится в непрерывной эксплуатации.

Апробация работы:

Основные результаты диссертационной работы докладывались и обсуждались на 19-ой Международной конференции по лазерному и акустическому зондированию (1998 г., Аннаполис, США), 5-ом Международном симпозиуме по оптике атмосферы и океана (1998 г., Томск), 20-ой Международной конференции по лазерному и акустическому зондированию (2000 г., Виши, Франция), 5-ом Международном симпозиуме по тропосферному профайлингу (2000 г., Аделаида, Австралия), 11-ой Международной школе по квантовой электронике (2000 г., Варна, Болгария), 21-ой Международной конференции по лазерному и акустическому зондированию (2002 г., Квебек, Канада), 6-ом Международном симпозиуме по тропосферному профайлингу (2003 г., Лейпциг, Германия), 22-ой Международной конференции по лазерному и акустическому зондированию (2004 г., Матера, Италия).

Структура работы:

Диссертация состоит из ведения, четырех глав, заключения, списка литературы, включающего 91 наименование, и приложения, содержит 235 страниц, 153 рисунка и 5 таблиц.

Во введении обоснована актуальность темы, сформулирована цель работы, изложены новые научные результаты и положения, выносимые на защиту.

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

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

В третьей главе проводится исследование физических основ метода фильтрации интерферометром Фабри-Перо вращательного спектра комбинационного рассеяния света молекулами азота (и кислорода) в задаче измерения температуры атмосферы. Построена математическая модель интерферометра. Рассматриваются вопросы согласования функции пропускания интерферометра и вращательного спектра комбинационного рассеяния. Дан анализ влияния неидеальной плоскостности зеркал интерферометра и величины угла расходимости падающего на интерферометр пучка на качество фильтрации помехи рассеянного фона. Выполнены оценки предела снижения необходимого времени накопления. Описан расчет параметров интерферометра по критерию минимизации времени измерений (за счет снижения уровня помехи рассеянного фона), предложенному во второй главе работы. Определяются границы применимости метода фильтрации рассеянного фона дневного неба в промежутках между линиями вращательного спектра комбинационного рассеяния атмосферных молекул с использованием интерферометра Фабри-Перо. Рассматриваются различные варианты формирования спектральных каналов и различные варианты схемы фильтрации фона.

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

В заключении приводятся основные выводы по результатам диссертационной работы.

На защиту выносятся следующие положения:

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

2. Интерферометр Фабри-Перо, обеспечивающий максимальное снижение времени измерений температуры при отношении сигнал-фон, равном 0.1, обеспечивает, при отношении сигнал-фон в интервале от 1 до 0.01, снижение времени измерений не менее 87 % от предельного уровня, достижимого при выборе параметров интерферометра для каждого конкретного значения отношения сигнал-фон.

3. В лидарной схеме измерений наличие упорядоченных по апертуре дефектов поверхности зеркал интерферометра Фабри-Перо приводит к зависимости пропускания спектральных каналов лидара от дальности зондирования.

Достоверность и обоснованность результатов диссертационной работы:

На основе теории спонтанного комбинационного рассеяния выполнены оценки характеристик температурного канала лидара, найденные значения подтверждаются оценками, полученными из натурного эксперимента. В частности, значение температурной чувствительности отношения вращательных откликов комбинационного рассеяния, полученное из расчета, хорошо (3 %) согласуется с величиной, полученной по результатам натурных измерений [14]. Вертикальные профили температуры, восстановленные по результатам ли-дарных измерений, с хорошей точностью (1-2 градуса) ложатся на профили температуры, полученные при помощи метеорологических шаров-зондов, запускавшихся параллельно с проведением лидарных измерений [5, 15-17].

Результаты расчета параметров полос пропускания температурного канала лидара, выполненного в соответствии с предложенным в работе критерием, совпадают с результатами подобных расчетов, представленных в работе [6], выполненных для случая измерений в отсутствии помех в соответствии с критерием минимизации ошибки измерений.

Математическая модель интерферометра Фабри-Перо построена на основе теории многолучевой интерференции. Оценка уровня подавления помехи рассеянного фона дневного неба интерферометром, полученная с применением принятой модели, с хорошей точностью (2 %) совпадает со значением оценки этого параметра, найденным из натурных измерений. Пропускание интерферометра для лидарных откликов в каналах вращательного комбинационного рассеяния, рассчитанное по модели, хорошо согласуется с данными, полученными в натурном эксперименте (отличие 5 %).

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

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

Основные результаты, представленные в диссертации, опубликованы в работах [5], [14], [15], [16], [51], [52], [56], [57], [79], [86], [87], [88].

Автор выражает глубокую благодарность кандидатам физико-математических наук Ю.Ф. Аршинову и С.М. Бобровникову за руководство работой, кандидатам физико-математических наук Б.В. Каулю и А.И. Надееву за участие в обсуждении полученных результатов, директору Института тропосферных исследований (Лейпциг, Германия) профессору Й. Хайценбергу и сотрудникам лидарной группы института докторам А. Ансманну, У. Вандингер, Д. Альтхаузену, И. Маттис, Д. Мюллеру и инженеру X. Хубэ, сотрудникам лидарной группы Федерального политехнического института (Лозанна, Швейцария) докторам В. Симеонову, И. Балину, Б. Калпини, X. ван ден Бергу за плодотворное сотрудничество при проведении совместных экспериментальных исследований.

226

223 Заключение

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

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

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

- При относительной спектральной плотности помехи в пределах от 0 до 1 (в единицах спектральной плотности отклика рассеяния для шестой линии стоксовой полосы вращательного спектра комбинационного рассеяния азота при температуре газа 220 К) минимальное время накопления, необходимое для достижения заданной точности измерений температуры по соотношению интенсивности линий вращательного спектра комбинационного рассеяния азота, реализуется при выборе параметров полос пропускания каналов для спектральной плотности помехи 3.5 -Ю"3. При тех же условиях, но измерении температуры по соотношению интенсивности линий вращательного спектра комбинационного рассеяния кислорода, параметры полос пропускания каналов следует выбирать для относительной спектральной плотности помехи, равной 2.6-10"3. При выделении в спектральных каналах линий спектра азота и линий спектра кислорода, параметры полос пропускания каналов следует выбирать для относительной спектральной плотности помехи 5.6 • 10-3.

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

2. Исследованы принципиальные возможности метода снижения необходимого времени накопления за счет фильтрации помехи рассеянного фона дневного неба в промежутках между линиями вращательных спектров комбинационного рассеяния молекул атмосферного воздуха с применением интерферометра Фабри-Перо. Построена математическая модель интерферометра, учитывающая физические аспекты задачи фильтрации фона. Выработан критерий согласования параметров интерферометра Фабри-Перо, и в соответствии с ним определены границы применимости метода фильтрации рассеянного фона дневного неба в промежутках между линиями вращательного спектра комбинационного рассеяния атмосферных молекул с использованием интерферометра Фабри-Перо. Определены области значений параметров интерферометра и диапазон значений относительной спектральной плотности помехи, при которых эффективна фильтрация фона интерферометром. Показано, что:

- При измерении температуры атмосферы по вращательным спектрам комбинационного рассеяния атмосферного азота (или кислорода) использовать интерферометр Фабри-Перо для фильтрации помехи рассеянного фона дневного неба следует в том случае, если отношение сигнал-фон лежит в пределах диапазона от 1 до 0.01. При отношении сигнал-фон больше единицы эффект от фильтрации фона интерферометром дает выигрыш во времени накопления не более 34 %, при отношении сигнал-фон меньше ста фактор эффективности снижения времени накопления выходит в насыщение, что связано с ограничением степени подавления помехи рассеянного фона за счет конечной ширины выделяемых спектральных линий.

- При выделении интерферометром Фабри-Перо линий вращательного спектра комбинационного рассеяния (азота или кислорода), максимальная эффективность снижения времени измерений температуры газа, по интенсивности выделяемых линий, составляет 68 %, 86 % и 96 % при отношении сигнал-фон 0.01, 0.1 и 1, соответственно, от предельного уровня, достижимого при идеальном совмещении каждой из выделяемых спектральных линий с максимумом функции пропускания. Соответственно: а) Отсутствие возможности непрерывной перестройки длины волны возбуждающего излучения (что характерно для широкого ряда лазерных источников) некритично с точки зрения методики фильтрации рассеянного фона с применением интерферометра Фабри-Перо - необходимое согласование положения спектральных линий и максимумов функции пропускания интерферометра можно осуществить без существенной потери пропускания перестройкой базы интерферометра. б) Принципиально неустранимое рассогласование гребенки линий вращательного спектра комбинационного рассеяния и гребенки максимумов функции пропускания интерферометра Фабри-Перо, являющееся следствием эффекта центробежного растяжения рассеивающих молекул, так же может быть скомпенсировано без существенной потери пропускания перестройкой базы интерферометра.

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

- При отношении сигнал-фон в пределах от 0 до 0.01, интерферометр, обеспечивающий максимальное снижение времени измерений при отношении сигнал-фон 0.1, обеспечивает снижение времени измерений не менее 87 % от предельного, достижимого при выборе параметров интерферометра для каждого значения отношения сигнал-фон.

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

- Эффективность снижения времени накопления, отличающаяся не более чем на 9 % от той, что позволил бы получить интерферометр с идеально плоскими зеркалами, реализуется при дефектной резкости функции пропускания интерферометра больше 100.

- Фильтрация помехи рассеянного фона интерферометром при дефектной резкости функции пропускания интерферометра меньше 20 малоэффективна (при отношении сигнал-фон 0.01 эффективность снижения времени накопления в 2.4 раза меньше предельно возможной).

- При изменении величины дефектной резкости функции пропускания интерферометра от 20 до 100 оптимальная величина отражательной резкости функции пропускания меняется в пределах примерно от 14 до 28, соответственно, область оптимальных значений коэффициента отражения зеркал интерферометра лежит в пределах от 0.8 до 0.9.

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

- При дефектной резкости функции пропускания интерферометра равной 100 и угловой резкости функции пропускания интерферометра равной 100 фильтрация интерферометром Фабри-Перо рассеянного фона дневного неба в промежутках меду линиями вращательного спектра комбинационного рассеяния азота позволит снизить время измерений температуры атмосферы в 6.6 раза при отношении сигнал-фон 0.01, 3.8 раза при отношении сигнал-фон 0.1 ив 1.34 раза при отношении сигнал-фон 0.1.

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

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

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

3. Проведены экспериментальные исследования, подтверждающие теоретические оценки эффективности фильтрации интерферометром Фабри-Перо помехи рассеянного фона дневного неба в каналах вращательного комбинационного рассеяния. Полученные экспериментальные данные показали, что применение интерферометра Фабри-Перо позволяет значительно снизить уровень помехи рассеянного фона дневного неба в каналах вращательного КР (с использовавшимся в эксперименте интерферометром, до 16 раз), а так же существенно улучшить отношение сигнал-фон (с использовавшимся интерферометром, до 6 раз).

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

1. Arshinov Yu.F., Bobrovnikov S.M., ZuevV.E., Mitev V.M. Atmospheric temperature measurements using a pure rotational Raman lidar // Aplied Optics. 1983. - V. 22. - # 19. -P. 2984-2990.

2. Nedeljkovich D., Hauchecorne A., Chanin M.L. Rotational Raman lidar to measure temperaturefrom the groung to 30 km // IEEE Trans. Geosci. Remote Sens. 1993. - V.31. - #1. - P. 90101.

3. Vaughan G., Wareing D.P., Pepler S.J., Thomas L., Mitev V. Atmospheric temperature measurements made by rotational Raman scattering // Applied Optics. 1993. - V.32, #15, P. 2758-2764.

4. Behrendt A., Reichardt J. Atmospheric temperature profiling in the presence of clouds with apure rotational Raman lidar by use of an interference-filter-based polychromator // Applied Optics. 2000. - V.39. - #9. - P. 1372-1378.

5. Аршинов Ю.Ф. Разработка и исследование лидарного метода измерения температуры атмосферы по вращательным спектрам КР молекул воздуха: Диссертация на соискание ученой степени кандидата физико-математических наукю Томск, 1982. - 139 с.

6. ПлачекГ. Рэлеевское рассеяние и Раман-эффект. Харьков: ГНТИ Укр.ССР, 1935.173 с.

7. Butcher R.J., Willetts D.V., Jones W.J. On the use of Fabry-Perot etalon for the determinationof rotational constants of simple molecules the pure rotational Raman spectra of oxygen and nitrogen // Proc. Roy. Soc. Lon. A. - 1971. - V. 324, P. 231-245.

8. Малышев В.И. Введение в экспериментальную спектроскопию. М.: Наука, 1979. - 478 с.

9. Тарасов К.И. Спектральные приборы. Ленинград: Машиностроение, 1977, - 387 с.

10. Нагибина И.М. Интерференция и дифракция света. Ленинград: Машиностроение, 1974. -359 с.

11. Зайдель А.Н. Основы спектрального анализа. М.: Наука, 1965. - 322 с.

12. MonzonJ.J., Sanchez-Soto L.L., BernabeuE. Influence of coating thickness on the performance of a Fabry-Perot interferometer // Applied Optics. 1991. - V.30. - #28. - P 4126-4132.

13. Mattis I., Ansmann A., Althausen D., JaenishV., Wandinger U., Mueller D., Arshinov Yu., Bobrovnikov S., Serikov I. Raman lidar observations of relative humidity //21st International laser radar conference: 8-12 July 2002. Quebec, 2002, P. 27-31.

14. Mattis I., Ansmann A., Althausen D., JaenishV., Wandinger U., Mueller D., Arshinov Yu.F., Bobrovnikov S.M., Serikov I.B. Relative humidity profiling in the troposphere with a Raman lidar //Applied Optics. 2002. - V.41. - #30. - P. 6451-6462.

15. Valley S.L. Handbook of Geophysics and Space Environments. Belford, Massachusetts: Air Force Cambridge Research Laboratories, 1965. - 692 p.

16. Хинкли Э.Д. Лазерный контроль атмосферы. М.: Мир, 1979. - 416 с.

17. Hauchecorne A., Chanin M.-L. Density and temperature profiles obtained by lidar between 35 and 75 km // Geophys. Res. Lett. 1980. - #7, P. 565-568.

18. ShibataT., Kobuchi M., Maeda M. Measurement of density and temperature profiles in the middle atmosphere with a XeF lidar // Applied Optics. 1986. - V.25. - #5. - P. 685-688.

19. Evans K.D, Melfi S.H., Ferrare R.A., Whiteman D.N. Upper tropospheric temperature measurements with the use of a Raman lidar // Applied Optics. 1997. - V.36. - #12. - P. 2594-2602

20. Gross M.R., McGee T.J., Ferrare R.A., Singh U.N., Kimvilakani P. Temperature measurements made with a combined Rayleigh Mie and Raman lidar // Applied Optics. -1997. V.36. - #24. P. 5987-5995.

21. Конингстайн Г. Введение в теорию комбинационного рассеяния. М.: Мир, 1976. -192 с.

22. ВеберА. Спектроскопия комбинационного рассеяния света в газах и жидкостях. М.: Мир, 1982.-373 с.

23. Gu Y.Y., Gardner Ch.S., Castleberg P A., Papen G.C., Kelley M.C. Validation of the Lidar InSpace Technology Experiment: stratospheric temperature and aerosol measurements // Applied Optics. 1997. - V.36. - #21. - P. 5148-5157

24. Bills R.E., Gardner C.S., She C.-Y. Narrowband lidar technique for sodium temperature and Dopier wind observation of the upper atmosphere // Optical Engineering. 1991. - V.30. -#1.-P. 13-21.

25. Fricke K.H., von Zahn U. Mesopause temperatures derived from probing the hyperfine structure of the D2 resonance line of sodium by lidar // J. Atmos. Terr. Phys. 1985. - V.47. -P. 499-512.

26. Neuber R., von der Gathen P., von Zahn U. Altitude and temperature of the mesopause at 69° N latitude in winter // J. Geophys. Res. 1988. - V.93. - P. 11093-11098.

27. Grund C.J., Eloranta E.W. University of Wisconsin high spectral resolution lidar // Optical Engineering. 1991. - V.30. - #1. - P. 6-12.

28. She C.-Y., YuJ.R., LatifiH., Bills R.E. High-spectral-resolution fluorescence light detection and ranging for mesospheric sodium temperature measurements // Applied Optics. 1992. -V.31.-#12.-P. 2095-2106

29. Fioccio G., Beneditti-Michelangeli G., Maischenberger K., Madonna E. Measurement of temperature and aerosol to molecule ratio in the troposphere by optical radar // Nature (London) Phys. Sci. 1971. - #229. - P. 78-79.

30. Lading L., Jensen A.S. Estimating the spectral width of a narrowband optical signal // Applied Optics. 1980. - V. 19. - #16. - P. 2750-2756.

31. Schwiesow R.L., Lading L. Temperature profiling by Rayleigh-scattering lidar // Applied Optics. 1981.- V.20. - #11. — P. 1972-1978.

32. Shimizu H., Lee S.A., She C.-Y. High-spectral-resolution lidar system with atomic blocking filters for measuring atmospheric parameters // Applied Optics/ 1983. - V.22. - #9. P. 13731381.

33. Shimizu H., Noguchi K., She C.-Y. Atmospheric temperature measurement by a high-spectral-resolution lidar//Applied Optics. 1986. - V.25. - #9. - P. 1460-1466.

34. Lehmann F.J., Lee S.A., She C.-Y. Laboratory measurements of atmospheric temperature and backscatter ratio using high-spectral-resolution lidar technique // Optics Letters. 1986. -V.l 1. - #9. P. 563-565.

35. She C.-Y., Alvarez R.J. II, Cardwell L.M., Krueger D.A. High-spectral-resolution Rayleigh— Mie lidar measurements of aerosol and atmospheric profiles // Optics Letters. 1992. - V.l 7. -#7.-P. 541-543.

36. Voss E., Weitkamp C., Michaelis W. Lead-vapor filters for high-spectral-resolution temperature lidar// Applied Optics. 1994. - V.33. - #15. -P. 3250- 3260.

37. Mason J.B. Lidar measurement of temperature: a new approach // Applied Optics. 1975. -V.14. - #1. — P. 76-78.

38. Rosenberg A., Hogan D.B. Lidar technique of simultaneous temperature and humidity measurements: analysis of Mason's method // Applied Optics. 1981. - V.20. - #19. - P. 32863288.

39. Shwemmer G.K., Wilkerson T.D. Lidar temperature profiling: performance simulation of Mason's method // Applied Optics. 1979.- V. 18. - #21. - P. 3539-3541.

40. Wulfmeyer V. Ground-Based Differential Absorption Lidar for Water-Vapor and Temperature Profiling: Development and Specifications of a High-Performance Laser Transmitter // Applied Optics. 1998. - V.37. - #18. - P. 3804-3824.

41. Cooney J. Measurement of atmospheric temperature profiles by Raman backscatter // J. Appl. Meteorol. 1972. - V. 11. - # 1. - P.108-112.

42. Cohen A., Cooney J.A., GellerK.N. Atmospheric temperature profiles from lidar measurements of rotational Raman and elastic scattering // Applied Optics. 1976. - V.15. -#11.-P. 2896-2901.

43. Penney M., Peters RL.St., Lapp M. Absolute rotational Raman cross sections for N2, O2, and C02 // J.Opt.Soc.Am. -1974. V.64. - #5. - P. 712-716.

44. Булдаков M.A., Матросов И.И., Попова Т.Н. Определение анизотропии тензора поляризуемости молекул О2 и N2 // Оптика и спектроскопия. 1979. - т.46. - вып.5. - С. 867-869.

45. Кондиленко И.И., Коротков П.А., Клименко В.А., Голубева Н.Г. Абсолютные сечения КР вращательных линий азота и кислорода // Оптика и спектроскопия. 1980. - т.48. -вып.4.-С. 745-751.

46. Cooney J.A. Atmospheric temperature measurement using a pure rotational Raman lidar: comment // Applied Optics. 1984. -V.23. - #5. - P. 653-654.

47. Petri K., SalikA., Cooney J.A. Variable-wavelength solar-blind Raman lidar for remote measurement of atmospheric water-vapor concentration and temperature // Applied Optics.-1982.-V.21.-#7.-P. 1212-1218.

48. Zeyn J., Lahmann W., Weitkamp C. Remote daytime measurements of tropospheric temperature profiles with a rotational Raman lidar // Optics Letters. 1996. - V.21. - #16. — P. 1301-1303.

49. Аршинов Ю.Ф., Бобровников C.M., Попов А.Г., Шелефонтюк Д.И., Шумский В.К. Дальность обнаружения малых загрязняющих компонентов атмосферы в солнечнослепой области спектра // Оптика атмосферы и океана. 1994. - т. 7. - № 8. - С. 1132-1138.

50. McKay J.A. Single and tandem Fabry-Perot etalons as solar background filters for lidar // Applied Optics/ -1999. V.38. - #27. - P. 5851-5858.

51. Trauger J.T., Roesler F.L. Fabry-Perot spectrometer adjustment for the compensation of Doppler shift from rapidly rotating and rapidly flowing sources // Applied Optics. 1972. -V.l 1. - #9. - P. 1964-1969.

52. Daehler M., Roesler F.L. High contrast in a polyetalon Fabry-Perot spectrometer // Applied Optics. 1968. - V.7. - #6. - P. 1240-1241.

53. Skinner W.R., Hays P.B., Abreu V.J. Optimization of a triple etalon interferometer I I Applied Optics. -1987. V.26. - #14. - P. 2817-2827.

54. GehrsitzS., SiggH., Siegwart H., KriegerM., Heine C., MorfR., Reinhart F.K., Martin W., Rudigier H. Tandem triple-pass Fabry Perot interferometer for applications in the near infrared //Applied Optics. 1997. - V.36. - #22. - P. 5355-5361.

55. Serabyn E., Benford D.J., Wu S., Pardo J.R. Dual Fabry Perot Filter for Measurement of CO Rotational Spectra: Design and Application to the CO Spectrum of Venus // Applied Optics. -2000. V.39. - #34. - P. 6448-6452.

56. HaysP.B., RobieR.G. A technique for recovering Doppler line profiles from Fabry-Perot interferometer fringes of very low intensity // Applied Optics. 1971. - V.10. - #1. - P. 193200.

57. Skinner W.R., Hays P.B. Observation of molecular rotational bands with a single Fabry-Perot interferometer // Applied Optics. 1986. - V.25. - #23. - P. 4272-4279.

58. Arshinov Yu., Bobrovnikov S. Use of a Fabry-Perot interferometer to isolate pure rotational Raman spectra of diatomic molecules // Applied Optics. 1999. - V.38. - #21. - P. 4635-4638.

59. Armstrong R.L. Rotational Raman Interferometric technique to measure gas temperatures // Applied Optics. -1975. V.14. - #2. - P. 383-387.

60. Плачек Г. Рэлеевское рассеяние и Раман-эффект. Харьков: ГНТИ Укр.ССР, 1935. -173 с.

61. Ярив А. Квантовая электроника и нелинейная оптика. М.: Сов.Радио, 1973.-455 с.

62. Аршинов Ю.Ф., Бобровников С.М., Сапожников С.В. О методе лидарного измерения температуры атмосферы по отноошению сигналов чисто вращательного спектра КР N2 и 02 // ЖПС. 1980. - т. XXXII. - вып.4. - С. 725-731.

63. Камминс Г.,Пайк Э. Спектроскопия оптического смешения и корреляция фотонов. М.: Мир, 1978. -584 с.

64. Горяинов В.Т., Журавлев А.Г., Тихонов В.И. Статистическая радиотехника: Примеры и задачи. М.: Сов. Радио, 1980. - 544 с.

65. McCIatchey R.A., FennR.W., SelbyJ.E.A., VolzF.E., GaringJ.S. Optical properties of the atmosphere // Environmental research papers. 1972. - #411. - p. 108.

66. Measurements with lidar systems during the international cirrus experiment 1989. Hamburg: Max-Planck-Institut fuer Meteorologie, №60, 1990. - 158 p.

67. Дейрменджан Д. Рассеяние электромагнитного излучения полидисперсными частицами: Пер. с англ. -М.: Мир, 1971. 165 с.

68. Arshinov Yu.F., Bobrovnikov S.M., ZuevV.E., MitevV.M. Atmospheric temperature measurements using a pure rotational Raman lidar // Aplied Optics. 1983. - V. 22. - # 19. -P. 2984-2990.

69. Мак-Картни Э. Оптика Атмосферы. M.: Мир, 1979. - 421 с.

70. Волковицкий О.А., Павлова JI.H., Петрушин А.Г. Оптические свойства кристаллических облаков. Ленинград: Гидрометеоиздат, 1984. - 198 с.

71. Борн M., Вольф Э. Основы оптики.- M.: Наука, 1970. 856 с.

72. Миберн Дж. Обнаружение и спектрометрия слабых источников света. М.: Мир, 1979. -304 с.

73. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Гос.Изд.Физ.Мат.Лит., 1963. - 1100 с.

74. Таблицы спектральных линий / Зайдель А.Н., Прокофьев В.К., Райский С.М., Славный В.А., Шрейдер Е.Я. М.: Наука, 1969.- 782 с.

75. Каталог цветного стекла: Каталог / Машиностроение, М., 1967. 63 с.

76. Радиационые процессы в средней атмосфере / Атмосфера. Ленинград: Гидрометеоиздат, 1991. - 509 с.

77. Kaufmann К. Photomultipliers detect rare events // Laser focus world. 2000. - V.36. - #10. -P. 77-81.91. http ://optoe 1 ectronics.perkinelmer.com.231