CC BY
53
УДК 681.785.574
А. Ю. Бойко, А. А. Григорьев, Г. В. Мацюк, А. Ю. Павлов, П. Е. Ш л ы г и н (Центральный научно-исследовательский испытательный институт), С. К. Дворук, М. В. Лельков, А. Н. Морозов, С. Е. Табалин, Г. В. Шишкин (МГТУ им. Н.Э. Баумана), В.Н.Корниенко(ИРЭРАН), И.В. Кочиков (НИВЦМГУим. М.В. Ломоносова), С. И. Светлич-н ы й (Институт энергетических проблем химической физики РАН)
ПРОБЛЕМА ИДЕНТИФИКАЦИИ И ОПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИЙ ЗАГРЯЗНЯЮЩИХ ВЕЩЕСТВ С ПОМОЩЬЮ ФУРЬЕ-СПЕКТРОРАДИОМЕТРА
На основе модели постоянства температур трассы и облака загрязняющих веществ предложен алгоритм решения задачи идентификации и определения концентраций загрязняющих веществ в атмосфере с помощью фурье-спектрорадиометра среднего класса разрешения. Рассмотрены процедуры обработки исходных интерферо-грамм и экспериментальных спектров, в том числе при малых яр-костных контрастах трасс. Приведены экспериментальные результаты лабораторных и натурных испытаний фурье-спектрорадио-метра.
Для лабораторных условий и стационарных фурье-спектрометров задача идентификации и определения концентрации веществ является методологически проработанной. Существует отлаженное программное обеспечение для решения подобных задач, включающее в себя набор математических процедур, позволяющих из экспериментального спектра выделить искомые спектральные компоненты и определить концентрации исследуемого вещества. Более того, имеются базы спектральных данных (постоянно пополняющиеся), с помощью которых решение задачи идентификации вещества, как правило, существенно упрощается. Естественно, что проведение подобных процедур требует достаточно высокой квалификации исследователя и занимает значительное время. Для мобильных фурье-спектрорадиометров (ФСР), предназначенных для работы в естественных условиях и в реальном масштабе времени, поставленная задача намного сложнее. Это обусловлено особенностями экспресс-анализа, жесткими временны ми рамками обработки и выдачи конечных результатов, а также непостоянством мишенной обстановки — вариациями температур подстилающей поверхности, изменением концентраций облаков загрязняющих веществ во времени и пространстве и т.п.
Цель настоящей работы — дeмoнcтpaция одного из возможных алгоритмов работы ФСР в реальном масштабе времени на основе накопленного опыта эксплуатации мобильного ФСР со спектральным разрешением около 4 см-1, рабочим диапазоном 700... 1400 см-1, временем скана менее 0,4 с и мгновенным полем зрения 1,5о х 1,5о. Фотоприемное устройство — фотосопротивление Cd-Hg-Te, охлаждаемое до 80 ^ с обнаружительной способностью около Б* = 1010 см- Гц1/2/Вт. ФСР снабжен электромеханической системой сканирования и высокочувствительным телевизионным визиром [1,2].
Основные соотношения, используемые при обработке спектров. При обработке спектров собственного излучения объектов в атмосфере можно исходить из наиболее общего уравнения переноса излучения, но для его решения, даже приближенного, требуются значительные вычислительные ресурсы, знание метеообстановки, привлечение специфических баз данных по излучательным свойствам подстилающих поверхностей и спектральным коэффициентам поглощения трасс. Поставленную задачу можно упростить, если предположить, что температуры облака загрязняющих веществ и трассы постоянны. Это предположение оправдано для приземных трасс наблюдений, где температурные перепады вдоль трассы не превышают нескольких градусов. Отсюда легко получить выражение для наблюдаемого спектра:
В = Р - ТрТЗВ(ВО - Р), (1)
где В = В(V) — экспериментальный спектр; Р = Р(^,Т) — функция Планка; ВО = ВО^) — спектр подстилающей поверхности (находящейся за облаком загрязняющих веществ по отношению к ФСР);
Ттр Ттр
(V) и тзв = тзв (V) — спектральные коэффициенты пропускания трассы и облака загрязняющих веществ соответственно; V — волновое число; Т — температура облака загрязняющих веществ и трассы. Дальнейшую обработку спектров можно проводить непосредственно согласно уравнению (1) составлением экспериментально наблюдаемого спектра из соответствующих спектров подстилающей поверхности и облака загрязняющих веществ. По-видимому, это не самый рациональный способ решения задачи для приземных трасс наблюдений, поскольку он требует обширной базы данных по спектральным свойствам различных подстилающих поверхностей и трасс наблюдения, но для ряда специфических задач данный подход оптимален.
Можно поступить иначе — предварительно снять спектр трассы (В0) при отсутствии облака загрязняющих веществ ( тзв = 1):
Во = Р - Ттр (ВО - Р). (2)
Комбинируя выражения (1) и (2), получим соотношения
B = P + тзв(Bo - P),
или
B - P
ТЗВ = B0-P ■ (3)
Можно получить соотношение, аналогичное выражению (3), но для условий, когда температура облака загрязняющих веществ (T1) и усредненная температура по трассе (T2) не совпадают. Проделав несложные алгебраические выкладки, получим в итоге
= B - P2 + Т-тр(р2 - Pi) (4)
ТЗВ Bo - P2 + Трр (P2 - Pi) ' ^
где P1 = P1(v,T1) и P2 = P2(v,T2) — функции Планка. Однако даже для такой сравнительно простой модели, описываемой выражением (4), подбор параметра т^ (v) требует дополнительных итерационных алгоритмов с привлечением баз данных типа HITRAN или GEISHA [3, 4]. Кроме того, требуется, чтобы были известны температуры трассы и облака загрязняющих веществ, что в натурных условиях не всегда возможно с необходимой точностью. По этим причинам в качестве основы алгоритма решения задачи было выбрано соотношение (3). Для стелющихся (приземных) облаков загрязняющих веществ выражение (3) описывает общую процедуру обработки спектров, полученных с помощью ФСР. При первом сканировании полученный спектр записывается в базу данных и рассматривается в дальнейшем как Bo, при последующих сканах регистрируются спектры B. Далее в соответствии с выражением (3) вычисляются спектральный коэффициент пропускания облака загрязнителей и, следовательно, спектральная оптическая плотность D:
n
D(v) = - ln(T3B(v)) = Y^ №в,
i=1
где n — число веществ в облаке загрязняющих веществ; о — сечение поглощения i-го вещества при заданном волновом числе (v); N/3B — интегральная концентрация загрязняющих веществ. Таким образом, чисто теоретически задача сводится к решению стандартной обратной задачи спектроскопии.
Предварительная обработка спектров. В отличие от классических спектральных приборов, фурье-спектрометр не дает спектральное распределение яркости излучения объекта наблюдения, а формирует автокоррелограмму, которая с учетом чувствительности фотоприемника S(v), спектральной функции пропускания оптического тракта F(v)
и функции фазовых искажений прибора ) может быть задана выражением
+5
I(ж) = Я0(V)В(V)ехр(г(2^ж - ф(и))) ¿V, (5)
-5
где I(ж) — регистрируемый сигнал двусторонней интерферограммы при разбалансе плеч интерферометра, равном ж; Я0(V) = 5(V(V); В (V) — спектр излучения объекта наблюдения; 6 — максимальное смещение подвижного зеркала интерферометра. При использовании в ФСР высокочувствительных охлаждаемых фотоприемников возникает необходимость учета собственного (фонового) излучения прибора, поскольку мощность входного сигнала от объекта исследования может быть сопоставима с собственным приборным излучением. B работе [5] показано, что регистрируемую интерферограмму объекта наблюдения с учетом фонового излучения можно представить в виде +5
I(ж)=2^УЯо^)(В(V) -)Р(V, Тиф))ехр(г(2^ж-^(у)))^, (6)
-5
где собственный фоновый (приборный) спектр является произведением эффективной излучательной способности £ИФ^) на функцию Планка Р(V, Тиф) при температуре Тиф, равной усредненной температуре деталей интерферометра. Естественно, что результат восстановления спектра по выражению (6) является комплексным:
В(V) = Яо^) (В(V) - ^иф(V)Р(V, Тиф)) ехр(-гф^)). (7)
Из выражений (6) и (7) следует, что предварительная обработка спектров должна включать в себя две обязательные стадии — фазовую коррекцию интерферограмм и алгоритм исключения из получаемого спектра фоновой составляющей излучения.
Коррекция фазы. Фазу ф^) можно определить непосредственно из соотношения (7) как аргумент комплексного спектра, но это определение, строго говоря, неоднозначно: из-за изменения фазовой функции на п рад или смены знака разностного сигнала могут получиться одинаковые спектры. Неоднозначность определения можно устранить с учетом того, что функция фазовых искажений — медленно изменяющаяся функция волнового числа. Таким образом, фаза спектра (7) должна представлять собой медленно изменяющуюся функцию с резкими изменениями на п рад в точках смены знака разности В(V) -гИФ(V)Р(V, Тиф). Функцию фазовой коррекции можно получить
выделением непрерывной составляющей фазы измеренного комплексного спектра. Это достижимо путем аппроксимации ) полиномом
2
невысокой степени (вданном случае — квадратичным: ) = ^ ¿),
¿=0
коэффициенты которого к определялись из условия минимума квадрата функции вида
^шах
Ф(к0,к1 ,к2) = У (ИеВ^)э1п(ф(^)) + 1шВ^)соэ(ф(^))) ¿V,
^шт
где ^ах — минимальное и максимальное волновые числа регистрируемого спектра.
Учет фоновой составляющей сигналов. Для исключения в исходных спектрах фоновой (приборной) составляющей необходимо, чтобы были известны функции Р0 (V), гИФ (V) и температура Тиф деталей интерферометра. Эти функции можно найти из измерений, проведенных для эталонных источников инфракрасного излучения при их различных температурах и различных собственных температурах интерферометра ФСР.
Пусть имеется несколько измерений спектров В^), выполненных для источников излучения с абсолютными температурами Т и излуча-тельными способностями = 1,...,п, причем при каждом измерении внутренним датчиком фиксировалась температура ТИФ корпуса ин-терферометра. Тогда с учетом переопределенности системы (при п > 2) получаем для каждой точки по волновым числам следующие уравнения для определения Р0 и в = £ИФ^0:
п \ / п \ п
^г2Р2(ТгН ^0 - X] ^Р№)Р(ТИФ) в = X] ^Р(ВД, ¿=1 ) V ¿=1 ! ¿=1
п \ / п \ п
^¿Р(Тг)Р(ТИФ) Р0 - ^ Р2(ТИФ) в = Р(ЗДВг.
¿=1 ) \г=1 ) ¿=1
(8)
Решая систему (8) для каждого значения V, нaйдем функции ) и £ИФ (V), с помощью которых можно легко исключить из спектров фоновую составляющую. На рис. 1 приведены искомые функции ) и £ИФ (V) экспериментального макета ФСР.
Финишная обработка спектров. Опишем алгоритм обработки спектров, необходимость которой вызвана непостоянством положения наблюдаемых спектров, малыми температурными контрастами и, как следствие, пересечением фоновых и экспериментальных спектров. Задача данного этапа обработки — ранжировать спектры, как это следует из соотношения (3), согласно которому спектры В облака загрязняю-
Рис. 1. Аппаратные функции
ФСР: 1 — функция спектральной чувствительности ФСР (Я0), 2 — эффективная излучательная способность(£ИФ)
0,6
/7,4
0,2
\7 1
2 \
__ / \ i ■ i
700
1100
1300i>. см'
,-1
щих веществ и В0 трассы без облака загрязняющих веществ, а также функция Планка Р (у, Т) должны удовлетворять соотношениям
В0 < В < Р либо Р < В < В0 (9)
для каждого волнового числа из рассматриваемого интервала (вид неравенства зависит от знака яркостного контраста).
В зависимости от режима работы прибора использовались следующие алгоритмы ранжирования. Если температура облака известна с достаточной точностью, вычислялся температурный контраст
ДТ = - 1), (10)
4 V (P,P)
где Т — абсолютная температура облака, а скалярные произведения рассчитывались по формуле
Zmax
(A,B) = J A(v)B(v)dv.
При достаточно большом положительном температурном контрасте ДТ > 2... 3 К спектр сигнала лежал выше спектра функции Планка , поэтому наблюдался спектр поглощения, а спектр фона лежал еще выше. При большом отрицательном контрасте спектр фона оказывался ниже спектра сигнала. Если спектры фона и сигнала пересекаются между собой, то спектр фона подстраивался так, чтобы он лежал выше или ниже спектра сигнала. Для этого спектр фона представлялся в виде
Во (а, Ь, V) = (а + Ьу )В о (у ),
где В0 (у) — первоначально измеренный спектр, а константы а, Ь выбираются из следующих соображений: спектр фона должен лежать выше
(или ниже) спектра сигнала; на краях спектрального участка спектры сигнала и фона не должны сильно различаться.
С помощью итерационной процедуры подстройки изменялся параметр а таким образом, чтобы (с известным допуском) спектр фона лежал выше или ниже спектра сигнала; далее при фиксированном значении а параметр Ь выбирался таким образом, чтобы минимизировать разность спектров сигнала и фона на краях спектрального участка.
Коррекция температуры облака загрязняющих веществ. В ситуации, когда температурный контраст мал или температура облака не могла быть измерена с достаточной точностью, процедура подстройки спектра фона сопровождалась процедурой подстройки температуры. В рамках этой процедуры вместо функции Планка использовалось выражение (1 + ^)Р(Т, V), где подстройке подлежала температура Т и параметр с. Критерии подстройки функции Планка аналогичны критериям подстройки спектра фона — функция Планка должна лежать выше (или ниже) спектра сигнала, а на краях рабочего диапазона спектр сигнала и функция Планка не должны сильно различаться. В данном случае значение температуры снова подбиралось таким образом, чтобы выполнялось первое из условий (9), а параметр с выбирался таким образом, чтобы на краях спектра значения сигнала и функции Планка были близки.
Вследствие того, что в подобных случаях нельзя заранее предсказать, будет ли функция Планка лежать выше или ниже спектра сигнала, анализ проводился для обоих вариантов, т.е. производилась подстройка и сверху, и снизу, причем спектр фона подстраивался таким образом, чтобы выполнялось соотношение (9), как это описано выше.
Определение температурного контраста и коррекция спектра. После подстройки спектра фона (и, возможно, температуры) проводилось вычисление коэффициента поглощения в облаке загрязняющих веществ на основе соотношений (3), а также заново вычислялся температурный контраст по формуле (10). Обычно рассматривалось несколько возможных вариантов расположения спектра фона и функции Планка относительно спектра сигнала, как это описано ранее. Выбор одного из вариантов производился на основе результатов решения задачи распознавания в каждом случае.
Алгоритм идентификации и определения концентрации загрязняющих веществ. Процедура отбора веществ из базы данных согласно коэффициентам корреляции. Практика расчетов показала, что результаты распознавания менее надежны, если применяется корреляционный алгоритм распознавания одновременно многих веществ. С целью повышения надежности процедуры было реализовано предвари-
тельное определение ограниченного набора веществ, для которых проводится количественный анализ.
На предварительном этапе рассчитывались корреляции спектров (коэффициентов пропускания) веществ из базы данных и полученного по формуле (3) спектра падающего излучения:
г. - , (т - ^ т - т> (11) 1 л/((т. - Т.)2>((т - т)2>' ^ '
где т — измеренный спектр пропускания, а т. — спектр пропускания г-го вещества из базы данных. Скалярные произведения вычислялись следующим образом:
{А, Б) = J A(v)Б(v)w(v) dv,
vmin
а средние значения — по формуле
Zmax
А = A(v)w(v) dv,
при этом весовая функция ) отражала погрешности определения коэффициента пропускания в каждой точке спектра и нормировалась таким образом, чтобы
^тах
J ) dv — 1.
^тш
Таким образом, определялось вещество, обладающее максимальным коэффициентом корреляции.
Нахождение второго и последующих компонентов смеси загрязняющих веществ. В связи с тем, что при наличии нескольких веществ в облаке загрязняющих веществ корреляции более слабых примесей могут быть сильно искажены влиянием основного вещества, проводилась процедура исключения первого найденного вещества из исходного спектра. Для этого решалась задача определения его концентрации: необходимо найти значение С такое, при котором достигается минимум функции невязки
Ф(С) = (exp(-CD(v)) - т(v))2 w(v) dv, (12)
V
V
где Б(у) — оптическая плотность для указанного вещества. Для найденного значения С * из исходной функции вычитался вклад, соответствующий данному веществу:
т>) = Т (У)
ехр(—С * б (у ))'
После этого коэффициенты корреляции остальных веществ пересчи-тывались заново, и процедура исключения повторялась для второго, третьего и т.д. веществ, обладающих наибольшими коэффициентами корреляции. Таким образом, для количественного анализа может быть использовано заранее заданное число веществ. Описанная процедура может привести к определению меньшего количества веществ, если наложить требование, чтобы коэффициенты корреляции для них были положительны и превосходили некоторый заданный минимальный уровень.
Определение концентраций компонентов смеси веществ. После отбора М веществ проводился количественный анализ с помощью формулы, аналогичной формуле (12), т.е. минимизировалась функция вида
^ / / м \ \ 2
Ф(С1,..., См )=у (ехр( — ^стбт(у и — т (у )\ ы(у ) Лу, (13)
где т — номер вещества из базы данных. Минимизация нелинейного функционала (13) проводилась с учетом ограничений на неотрицательность коэффициентов С > 0, причем ограничения накладывались и сверху, так как при очень больших концентрациях коэффициенты С с достаточной точностью не могут быть определены и, наконец, величины концентраций веществ в газовой фазе не могут превышать значений насыщенного пара при данной температуре.
Альтернативным методом является переход от коэффициентов пропускания к оптическим плотностям, в результате чего рассматривается линейная система уравнений относительно оптических плотностей, но ее решение также требует применения итерационных алгоритмов, если используются рассмотренные ограничения.
Фактически вместо интеграла в правой части выражения (13) стоит сумма по числу точек спектра (обычно по нескольким сотням точек). Переопределенность системы позволяет оценить погрешность вычисляемого решения по стандартным формулам, построив квадратичную аппроксимацию функционала (13) в точке минимума.
Определение списка распознанных веществ. В окончательный список распознанных веществ включались те из них, для которых значения
концентраций превосходят ошибку их определения. Кроме того, для каждого вещества в базе данных задавалась некоторая пороговая концентрация, которая должна быть превышена, в противном случае вещество исключалось из списка.
База спектральных данных. База спектральных данных загрязняющих веществ является одним из основных компонентов программного обеспечения комплекса ФСР. Только с ее помощью возможно проведение процедур идентификации и определение концентраций загрязняющих веществ, в том числе при многокомпонентном загрязнении. Очевидно, что качество базы данных играет ключевую роль в решении поставленной задачи. Прямое использование имеющихся баз данных, в том числе коммерческих, не всегда возможно по нескольким причинам. Конструктивная особенность подвески подвижного зеркала, примененного в данном ФСР, приводит к умножению исходной интерфе-рограммы на специфичную функцию аподизации и, как следствие, к изменению вида полос и линий поглощения, формы которых могут не совпадать с представленными в исходной базе данных. Не всегда концентрации веществ в базах данных близки по своим величинам к наблюдаемым, а условия их получения — к наблюдаемым в естественных условиях. Именно по этим причинам была создана спектральная база данных основных загрязняющих веществ и их имитаторов, для чего использовалась однопроходовая оптическая кювета длиной 8 ми сечением 0,6 м2 с принудительным механическим перемешиванием и вентилированием объема. В качестве источника подсветки применялся широкоапертурный термостабилизированный источник инфракрасного излучения, близкий по параметрам к абсолютно черному телу с яр-костной температурой от +4 до +50°С. Для контроля величин концентраций в кювете применялись как спектральная методика, так и пробо-отборные методы аналитической химии.
Регистрацию спектров проводили классическим методом. Сначала снимали спектр кюветы, заполненной атмосферным воздухом, потом спектр фонового абсолютно черного тела с температурой, равной комнатной, и затем в кювету вводили требуемое количество вещества. Спектральный коэффициент пропускания рассчитывали по соотношению (3) с последующей его коррекцией, заключающейся в том, что из конечного спектра Б = ) исключали постоянную составляющую в рабочем диапазоне 700... 1400 см-1, после чего спектр вносился в исходную базу данных. Данная процедура является обязательной, так как при обработке экспериментальных спектров в натурных условиях программа на заключительном этапе работы автоматически исключает постоянную составляющую, если таковая возникала при обработке.
Рис.2. Спектральные коэффициенты пропускания: а - аммиака (интегральная концентрация 0,34 г/м2); б - гептила (3 г/м2); в - метанола (1,2 г/м2); г - изопро-панола (4,9 г/м2)
На рис. 2 приведены спектральные коэффициенты пропускания для ряда веществ из созданной базы данных.
Результаты экспериментов. Первым этапом экспериментальных работ было лабораторное тестирование ФСР в режиме обнаружения с использованием оптической кюветы. После ввода заданного количества паров вещества (или смеси веществ) контролировался момент срабатывания ФСР при обнаружении вещества (смеси), правильность детектирования вещества и/или компонентов смеси и определялись концентрации обнаруженных компонентов. При испытаниях с помощью кюветы регистрировались предельные минимальные обнаружи-мые концентрации как отдельных веществ (табл. 1), так и смесей типа метанол-аммиак и метанол-аммиак-ацетон. Во всех проведенных экспериментах были зарегистрированы все компоненты смесей с определением их концентраций при яркостном контрасте, не превышающем 5 ...7 К.
Таблица 1
Минимальные обнаружимые концентрации веществ при яркостном контрасте не более 5 ... 7 К
№ Вещество Минимальная обнару-жимая интегральная концентрация, мг/м2 Минимальная обнаружи-мая концентрация при размере облака L « 100 м, мг/м3
1 Гептил 60 0,6
2 Метанол 20 0,2
3 Аммиак 10 0,1
4 Ацетон 150 1,5
5 Фенол 200 2
6 Бензол 300 3
7 Толуол 800 8
Экспериментальные работы в натуральных условиях были проведены утром при минимальных яркостных контрастах трасс, т.е. в наиболее жестких условиях работы для ФСР. Мишенная обстановка создавалась принудительным распылением в воздухе заданного вещества на высотах, не превышающих ^ 50 м. Минимальный объем при распылении вещества (или его раствора) составлял около 0,4 л. Испытания проводились на расстоянии 0,5 ... 6 км от ФСР до места распыления.
Тестирование работы ФСР проводилось в одном из двух возможных режимов срабатывания: в режиме срабатывания по первому факту обнаружения и в режиме перманентной (непрерывной) записи результатов работы в файл.
В режиме сканирования одного или нескольких секторов на расстояниях 1 и 2,5 км от ФСР надежно обнаруживались пары аммиака. На рис. 3 приведены результаты перманентной записи данных в файл при обнаружении облака паров аммиака на статической трассе длиной 0,5 км при одиночном заряде аммиака массой 80-90 г (0,45 л 20 %-ного раствора аммиачной воды). Наблюдение облака велось при неизменных углах возвышения и места. На этой же трассе успешно проведены серии экспериментов по обнаружению смесей паров этило-
200
о ю ю зо м 50 г, с
Рис.3. Временная зависимость изменения концентраций облака аммиака (суммарная масса аммиака в облаке — около 90 г)
вого и изопропилового спирта, ацетона и аммиака в различных комбинациях. На трассе длиной 2,5 км наблюдение и фиксация параметров облака паров аммиака при разовом распылении 0,7 кг аммиачной воды проводилось непрерывно в течение более 40 мин.
Наблюдение облаков загрязняющих веществ на трассах различной протяженности позволило определить минимально регистрируемые концентрации паров аммиака. Результаты экспериментов приведены в табл. 2, там же приведены усредненные яркостные контрасты трасс, которые измерялись с помощью калиброванного панорамного тепловизора при диапазоне 8 ... 12 мкм с точностью определения яркостных температур около 0,1 К.
Таблица 2
Минимально обнаружимые концентрации паров аммиака в зависимости от расстояния от ФСР до облака загрязняющих веществ
Расстояние, км
Оценка минимальной обнаружимой концентрации, мг/м2
Температурный контраст, К
0,5 1
1,5 2,5 4 6
12 15 25 30 20 30
5 ...6
4
2 ...4
2
2... 3
2
Заключение. В настоящей работе представлен один из апробированных вариантов алгоритма работы комплекса фурье-спектроради-ометра в естественных условиях. Полученные экспериментальные результаты подтверждают перспективность применения фурье-спектро-радиометра для решения задач оперативного контроля состава воздушной среды.
СПИСОК ЛИТЕРАТУРЫ
1. Д в о р у к С. К., К о ч и к о в И. В., Морозов А. Н., Н а з о л и н А. Л., Павлов Д. А., Поздняков В. А., Соловьев А. В., Светличный С. И., Табалин С. Е. Применение фурье-спектрорадиометра для определения концентраций загрязняющих веществ в атмосфере // Оптический журнал. - 2000. - Т. 67. - № 3. - С. 37-42.
2. Дворук С. К., Корниенко В. Н., Кочиков И. В., Морозов А. Н., Поздняков В. А., Светличный С. И., Соловьев А. В., Т а б ал и н С. Е., Ш и ш к и н Г. В. Применение фурье-спектрорадиометрии для анализа необратимых химических и физических процессов в природе // Вестник МГТУ им. Н.Э. Баумана. Сер. "Естественные науки". - 2001. - № 1. - С. 87-101.
3.Rothman J. et al. The HITRAN database 1986 // Appl. Opt. - 1987. - V. 26. -P. 4058-4097.
4. C h e d i n A. et al. The GEISHA Spectroscopic line parameter data bank in 1984 // Ann. Geophys. - 1984. - V. 4. - P. 2.
5. Дворук С. К., Корниенко В. Н., Кочиков И. В., Л е л ь к о в М. В., Морозов А. Н., Поздняков В. А., Светличный С. И., Т а б алин С. Е. Обработка двусторонних интерферограмм с учетом собственного фонового излучения фурье-спектрорадиометра // Оптика и спектроскопия. - 2002. - Т. 93. - № 6. - С. 1048-1053.
Статья поступила в редакцию 13.02.2003
Андрей Юрьевич Бойко родился в 1968 г., окончил в 1990 г. Саратовское высшее военное инженерное училище химической защиты. Канд. техн. наук, старший научный сотрудник, начальник отдела Центрального научно-исследовательского испытательного института. Автор более 20 научных работ в области спектроскопии, оптической локации и прецизионных измерений.
A.Yu. Boiko (b. 1968) graduated from the Saratov Military Higher School for Chemical Defense in 1990. Ph.D. (Eng.), senior researcher, head of department of the Central Research Institute for Testing. Author of over 20 publications in the field of spectroscopy, optical radar and precision measuring.
Александр Александрович Григорьев родился в 1951 г., окончил в 1975 г. Военную академию химической защиты. Канд. техн. наук, старший научный сотрудник, начальник отдела Центрального научно-исследовательского испытательного института. Автор более 100 научных работ в области математического моделирования, спектроскопии и оптической локации.
A.A. Grigoriev (b. 1951) graduated from the Military Academy for Chemical Defense in 1975. Ph.D. (Eng.), senior researcher, head of department of the Central Research Institute for Testing. Author of over 100 publications in the field of mathematical simulation, spectroscopy and optical radar.
Григорий Владимирович Мацюк родился в 1958 г., окончил в 1987 г. Военную академию химической защиты. Старший научный сотрудник Центрального научно-исследовательского испытательного института. Автор более 50 научных работ в области аналитической химии, оптической локации и прецизионных измерений.
G.V. Matsyuk (b. 1958) graduated from the Military Academy for Chemical Defense in 1987. Senior researcher, deputy director of the Central Research Institute for Testing. Author of over 50 publications in the field of analytical chemistry, optical radar and high precision measuring.
Александр Юрьевич Павлов родился в 1958 г, окончил в 1981 г. Саратовское высшее военное инженерное училище химической защиты. Д-р хим. наук, старший научный сотрудник, зам. начальника Центрального научно-исследовательского испытательного института. Автор более 100 научных работ в области аналитической химии, спектроскопии и прецизионных измерений.
A.Yu. Pavlov (b. 1958) graduated from the Saratov Military Higher School for Chemical Defense in 1981. D.Sc. (Chemistry), senior researcher, deputy director of the Central Research Institute for Testing. Author of over 100 publications in the field of analytical chemistry, spectroscopy and high precision measuring.
Петр Евгеньевич Шлыгин родился в 1969 г., окончил в 1991 г. Саратовское высшее военное инженерное училище химической защиты. Старший инженер Центрального научно-исследовательского испытательного института. Автор более 10 научных работ в области оптической локации и прецизионных измерений.
P.Ye. Shlygin (b. 1969) graduated from the Saratov Military Higher School for Chemical Defense in 1991. Senior engineer of the Central Research Institute for Testing. Author of over 10 publications in the field of optical radar and high precision measuring.
Станислав Казимирович Дворук родился в 1958 г., окончил в 1986 г. МВТУ им. Н.Э. Баумана. Ведущий инженер Центра прикладной физики МГТУ им. Н.Э. Баумана. Автор 5 научных работ в области прецизионных измерений.
S.K. Dvoruk (b. 1958) graduated from the Bauman Moscow High Technical School in 1986. Leading engineer of Applied Physics Center of the Bauman Moscow State Technical University. Author of 5 publications in the field of high precision measuring.
Максим Викторович Лельков родился в 1976 г., окончил в 2000 г. МГТУ им. Н.Э. Баумана. Аспирант кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор 4 научных работ в области расчета и создания оптических систем.
M.V. Lelkov (b. 1976) graduated from the Bauman Moscow State Technical University in 2000. Post-graduate of "Physics" department of the Bauman Moscow State Technical University. Author of 4 publications in the field of computation and creation of optical
Андрей Николаевич Морозов родился в 1959 г., окончил в 1981 г. МВТУ им. Н.Э. Баумана Д-р физ.-мат. наук, профессор, зав. кафедрой "Физика" МГТУ им. Н.Э. Баумана. Автор более 100 научных работ в области прецизионных измерений и физической кинетики.
A.N. Morozov (b. 1959) graduated from the Bauman Moscow Higher Technical School in 1981. D.Sc. (Phys.-Math.), professor, head of "Physics" department of the Bauman Moscow State Technical University. Author of over 100 publications in the field of high precision measuring and physical kinetics.
Сергей Егорович Табалин родился в 1961 г., окончил в 1984 г. МВТУ им. Н.Э. Баумана. Научный сотрудник НИИ ПММ МГТУ им. Н.Э. Баумана. Автор более 30 научных работ в области физической оптики и прецизионных измерений.
S.Ye. Tabalin (b. 1961) graduated from the Bauman Moscow Higher Technical School in 1984. Researcher of "PMM" research institute of the Bauman Moscow State Technical University. Author of over 30 publications in the field of physical optics and high precision measuring.
Григорий Владимирович Шишкин родился в 1979 г., окончил МГТУ им. Н.Э. Баумана. Аспирант МГТУ им. Н.Э. Баумана. Автор 4 научных работ в области систем управления.
G.V. Shishkin (b. 1979) graduated from the Bauman Moscow State Technical University. Post-graduate of the Bauman Moscow State Technical University. Author of 4 publications in the field of control systems.
systems.
о
II
Владимир Николаевич Корниенко родился в 1966 г., окончил в 1989 г. МГУ им. М.В. Ломоносова. Канд. физ.-мат. наук, старший научный сотрудник Института радиотехники и электроники РАН. Автор более 50 научных работ в области вакуумной микроволновой электроники.
V.N. Kornienko (b. 1966), graduated from the Moscow State University n.a. Lomonosov in 1989. Ph.D. (Phys.-Math.), senior researcher of the Institute for Radio Electronics of the Russian Academy of Sciences. Author of over 50 publications in the field of vacuum microwave electronics.
Игорь Викторович Кочиков родился в 1959 г., окончил в 1982 г. МГУ им. М.В. Ломоносова. Канд. физ.-мат. наук, старший научный сотрудник НИВЦ МГУ. Автор более 50 научных работ в области математической физики и вычислительной математики.
I.V. Kochikov (b. 1959) graduated from the Moscow State University n. a. Lomonosov in 1982. PhD (Phys.-Math.), senior researcher of "NIVTs" center of the Moscow State University n.a. Lomonosov. Author of over 50 publications in the field of mathematical physics and computation mathematics.
Сергей Иванович Светличный родился в 1950 г., окончил в 1973 г. Московский физико-технический институт. Канд. физ.-мат. наук, директор филиала Института энергетических проблем химической физики РАН. Автор более 60 научных работ в области химической физики и спектроскопии.
S.I. Svetlichny (b. 1950) graduated from Moscow Physical and Technical Institute in 1973. Ph.D. (Phys.-Math.), director of Affiliation of the Institute for Energy Problems of Chemical Physics of the Russian Academy of Sciences. Author of over 60 publications in the field of chemical physics, spectroscopy.
УДК 621.373+535
В. Г. Марченко
ПОЛЯ ОТКРЫТЫХ ШИРОКОАПЕРТУРНЫХ ОПТИЧЕСКИХ РЕЗОНАТОРОВ
Изложен подход с применением ключевых систем к расчету полей широкоапертурных лазеров с резонаторами, образованными плоскими, выпуклыми и вогнутыми зеркалами, основанный на точном решении полного волнового уравнения в скалярном представлении. Сформулирован физический смысл спецфункций разделяющих систем координат.
Открытые оптические резонаторы являются важным элементом квантовых оптических генераторов — лазеров. Их функции не ограничиваются созданием положительной обратной связи в оптических генераторах: с их помощью формируются поля с определенными свойствами в ближней и дальней зонах излучения.
Исторически лазеры прошли путь от генераторов с малыми апертурами резонаторов (с диаметром апертуры а порядка нескольких миллиметров, длиной резонатора L около метра, числами Френеля N = a2/4AL . 1, где А — длина волны) и — по мере возрастания
