Например, Бобцов

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

КОСМИЧЕСКАЯ ИНФОРМАТИКА
УДК 621.383.7
В. Ф. ФАТЕЕВ, В. Н. ЛАГУТКИН, Ю. В. СЛЫНЬКО
ИСПОЛЬЗОВАНИЕ МАЛЫХ КОСМИЧЕСКИХ АППАРАТОВ СО СПЕКТРОСКОПИЧЕСКИМИ ПРИБОРАМИ
ДЛЯ МОНИТОРИНГА ГАЗООБРАЗНЫХ ВЫБРОСОВ В АТМОСФЕРЕ
Рассмотрены вопросы применения инфракрасных спектроскопических приборов на малых космических аппаратах (МКА) для мониторинга газообразных выбросов в атмосфере. Приведены оценки требуемого спектрального разрешения спектроскопических приборов. Показано, что более предпочтительно использовать на МКА интерферометры с субнанометровым спектральным разрешением, настраиваемые на наиболее интенсивные спектральные линии газовых образований. Приведены примеры спектров излучения газовых образований, пропускания атмосферы, фона; определены особенности задачи адаптивного выбора спектральных линий.
Ключевые слова: малые космические аппараты, мониторинг атмосферы, спектроскопия газов, интерферометры, фон.
В рамках исследования направлений развития и областей применения малых космических аппаратов (МКА) значительное внимание уделяется проблеме создания многоспектральных датчиков для разнообразных систем мониторинга, управления их работой и автоматической обработки полученных с их помощью данных. В работах [1, 2] эта проблема рассматривалась применительно к задачам дистанционного зондирования Земли (ДЗЗ) с МКА и, в частности, исследовались методы выбора количества спектральных каналов и их расстановки на оси длины волны. Благодаря тому что спектральные образы наблюдаемых при ДЗЗ объектов являются непрерывными и гладкими (в той или иной степени) функциями длины волны, оказалось возможным при решении задачи определения положения спектральных каналов на оси длины волны использовать методы ортогональных разложений функций [2]. Задача определения ширины полос спектральных каналов детально не рассматривалась, поскольку применительно к ДЗЗ значение этого параметра не является критическим.
Спектральные образы газовых образований в атмосфере как естественного, так и искусственного происхождения кардинально отличаются от образов наземных объектов, так как спектры газов являются линейчатыми и как функции длины волны представляют собой совокупность иглообразных пиков разной высоты. В настоящей работе эта особенность спектров газов используется как базовая при выработке методического подхода к определению принципов организации аппаратуры для мониторинга газообразных выбросов в атмосфере.
Как известно, взаимодействие оптического излучения с газообразной средой описывается уравнением переноса излучения [3]. Решение этого уравнения для случая отсутствия рассеяния можно представить в виде
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 4

8 В. Ф. Фатеев, В. Н. Лагуткин, Ю. В. Слынько

S
∫Bν (r(S), l) = Bν0 (l)Hν (0, S; l) + EνBB (T (r(P)))dHν (P, S; l)dP , 0

(1)

где Bν (r, l) — спектральная плотность яркости излучения как функция координат r и направления наблюдения l, Р — длина пути вдоль направления l (dР>0), S — полный путь от наблюдателя до объекта, ν = 1/λ — частота волны излучения (в см–1), λ — длина волны, EνBB (T (r)) — спектральная плотность излучения абсолютно черного тела с температурой Т, определяемая формулой Планка, Hν (P, S; l) — спектральная функция пропускания газообразной среды на пути от точки r(P) до точки r(S) (положение наблюдателя) вдоль направ-

ления наблюдения, определяемая соотношением



(P,

S

;

l)

=

exp

⎡⎢−

S
∫βν

(r(S′))dS′⎥⎤

,

⎣P ⎦

(2)

βν (r) — спектральный показатель поглощения.
Заметим, что функция Hν (P, S; l) является неубывающей функцией Р, причем
0 ≤ Hν (P, S;l) ≤ 1. Функция Bν0 (l) определяет спектральную плотность потока излучения внешнего источ-
ника, находящегося на линии наблюдения за газообразным объектом, т.е. фона подстилающей поверхности, облаков или звезд.
Выражение для спектрального показателя поглощения излучения в многокомпонентном газе имеет вид

∑βν (r) = ρk (r)βν′ k (T (r), P(r)) ,

(3)

k

где ρk (r) — плотность k-й компоненты газа, βν′ k (T , P) — коэффициент поглощения k-й ком-

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

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

∑βν′ k (T , P) = Iki Fk (ν − νki , γki ) ,

(4)

i

Iki — интенсивность излучения i-й линии, Fk (ν − νki , γki ) — контур линий, νki — несущая

частота, γki — полуширина контура. Наиболее часто используемой аппроксимацией контура

спектральных линий является лоренцевский контур (подробнее см. [4])

Fk

(ν − νki , γki )

=

1 π



γ ki − νki )2

+

γ 2ki

.

(5)

Интенсивность и полуширина линий зависят от температуры и давления. Для атмосфер-

ных газов типичная полуширина линий в нижней атмосфере составляет γ~10–1 см–1 (∆λ~10–1 нм

для средневолнового ИК-диапазона), в верхней — γ~10–2—10–3 см–1 (∆λ~10–2—10–3 нм).

Соотношения (1)—(5) представляют собой количественную модель сигналов и помех

для различных условий мониторинга газообразных объектов в атмосфере с МКА. Основную

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

рая, согласно выражениям (3)—(5), зависит от пространственного распределения состава,

температуры и давления газообразных объектов. В случае неразреженных газовых образова-

ний эта функция, в основном принимая значения около единицы, имеет очень узкие, порядка

ширины контуров спектральных линий, и глубокие, вплоть до нуля, „провалы“. В соответст-

вии с уравнением (1) это означает, что принимаемое излучение от газовых образований скон-

центрировано только в узких полосах (линиях), вне которых принимается излучение внешних

источников, являющееся в этом случае фоном наблюдения.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 4

Использование МКА для мониторинга газообразных выбросов в атмосфере

9

В случае разреженных газовых образований спектральная функция пропускания может не иметь глубоких „провалов“ даже в линиях с большой интенсивностью поглощения газов, при этом их собственное тепловое излучение может быть мало по сравнению с излучением фона. Тем не менее, если пространственное распределение фона относительно плавно изменяется или в случае, когда оно заранее известно, газовые образования можно наблюдать „на просвет“, выделяя области потемнения на изображениях фона, полученных в линиях поглощения газов.
Наблюдение и распознавание газообразных выбросов в атмосфере можно осуществлять на базе МКА, оснащенных спектроскопическими приборами. Каждый МКА в зависимости от высоты орбиты может осуществлять мониторинг газообразных выбросов в атмосфере в полосе шириной ~6000 км вдоль орбиты. Использование развитой орбитальной группировки МКА позволит обеспечить, во-первых, оперативный контроль атмосферы на всех широтах и, вовторых — возможность триангуляционных измерений, на основе которых можно определить пространственное положение выбросов и оценить их пространственную конфигурацию.
В качестве возможных спектроскопических приборов для МКА следует рассматривать фурье-спектрометры и перестраиваемые интерферометры [5]. Оба типа имеют высокое спектральное разрешение, до ~10–2 см–1 и меньше [6]. Сравнительный анализ этих приборов показал, что применение на МКА перестраиваемых интерферометров более предпочтительно. Во-первых, перестраиваемым интерферометром можно управлять, поочередно настраивая его на наиболее интенсивные спектральные (известные) линии газовых образований. Это обеспечивает большее быстродействие, что очень важно для реализации оперативного кругового обзора пространства с низкоорбитальных МКА. Во-вторых, интерферометры, настраиваемые на интенсивные спектральные линии газовых образований, имеют более высокую чувствительность, так как влияние фона в узких линиях сводится к минимуму. В-третьих, интерферометры могут работать в более широком поле зрения и объединяться с многоэлементными фотоприемниками.
При адаптивном определении набора спектральных линий, за которыми должен следить спектральный прибор, необходимо учитывать следующую информацию: возможный состав наблюдаемых газообразных объектов, спектральную функцию пропускания атмосферы, спектральную функцию яркости фона. На основе известных составов газовых образований можно заранее сформировать каталог интенсивных спектральных линий.
На рис. 1 проиллюстрированы расчетные фрагменты спектра излучения продуктов сгорания углеводородного топлива в верхней атмосфере в двух спектральных интервалах средневолновой части ИК-диапазона. Расчет выполнен на модели, описанной в работе [7].
На рис. 1, а хорошо видны отдельные линии и яркий дуплет линий паров воды на частоте ν≈2994,5 см–1, ширина линий ~0,1 см–1. Фрагмент спектра, показанный на рис. 1, б, интересен тем, что на нем в интервалах частот 3688—3689 и 3691—3692 см–1 наблюдаются группы плотно расположенных интенсивных линий с общей шириной ~0,5 см–1. Из рис. 1 видно, что для выделения излучения в спектральных линиях спектроскопические приборы должны обеспечивать разрешение ~(1—5)·10–1 см–1, т.е. на уровне десятых долей нанометра для средневолнового ИК-диапазона.
На рис. 2 приведен график спектральной функции пропускания атмосферы, рассчитанной на модели [8] в предположении, что газовое образование находится на высоте 10 км и зенитный угол наблюдения равен 45°.
„Провалы“ функции пропускания обусловлены линиями поглощения излучения молекулярными атмосферными составляющими, главным образом, водяным паром и углекислым газом, расположенными на пути от объекта до МКА.
На рис. 3 приведены графики расчетной спектральной плотности яркости фона облаков с высотой верхней границы 5 км в полосе частот 3680—3710 см–1 (в дБ к Вт/(м2·мкм·ср)). Расчет проведен на модели [9] с модернизированным модулем пропускания атмосферы [8].

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 4

10 В. Ф. Фатеев, В. Н. Лагуткин, Ю. В. Слынько
Спектр фона облаков имеет „обратную“ линейчатую структуру, в которой узкие и глубокие провалы (на 15—20 дБ) обусловлены линиями поглощения атмосферных газов.
а) I, (Вт⋅см)/ср 150
120
90
60
30

б) I, (Вт⋅см)/ср
200

2991 2992 2993 2994 2995 2996 2997 2998 ν, см–1

150

100

50

3686 3687 3688 3689 3690 3691 3692 3693 ν, см–1 Рис. 1
Нν
0,8
0,6
0,4
0,2

В, дБ –20

2975

2980 2985 Рис. 2

2990

2995 ν, см–1

–30

–40
–50 3685 3690 3695 3700 3705 ν, см–1 Рис. 3
Спектры пропускания атмосферы и яркости фона зависят от состояния атмосферы, параметров ее влажности и облачности. Совместный анализ спектральных функций излучения газовых образований, фона и пропускания атмосферы, подобных представленным на рис. 1—3,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 4

Использование МКА для мониторинга газообразных выбросов в атмосфере

11

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

СПИСОК ЛИТЕРАТУРЫ

1. Фатеев В. Ф., Миньков С. А. Новое направление развития МКА дистанционного зондирования Земли // Изв. вузов. Приборостроение. 2004. Т. 47, № 3. С. 18—22.

2. Мальцев Г. Н., Козинов И. А., Фатеев В. Ф. Методы выбора наиболее информативных спектральных каналов при дистанционном зондировании Земли с малых космических аппаратов // Изв. вузов. Приборостроение. 2007. Т. 50, № 6. С. 23—31.

3. Чандрасекар С. Перенос лучистой энергии. М.: ИЛ, 1953. 431 с.

4. Зуев В. Е. Распространение видимых и инфракрасных волн в атмосфере. М.: Сов. радио, 1970. 496 с.

5. Курбатов Л. Н. Оптоэлектроника видимого и инфракрасного диапазонов спектра. М.: Изд-во МФТИ, 1999. 320 с.

6. Справочник по инфракрасной технике / Под ред. У. Вольфа и Г. Цисиса. Т. 4. М.: Мир, 1999. 472 с.

7. Лагуткин В. Н., Слынько Ю. В. Метод расчета спектральных яркостей потока смеси газа и частиц // Тр. XLX науч. конф. МФТИ „Современные проблемы фундаментальных и прикладных наук“. Ч. 1. М.: МФТИ, 2007. C. 119—121.

8. Слынько Ю. В. Модель расчета переноса излучения на основе открытой базы данных HITRAN // Вопросы радиоэлектроники. Сер. РТЛ. 2007. Вып. 4. С. 5—11.

9. Лагуткин В. Н., Лукьянов А. П., Подкорытов Е. Н., Репин В. Г. Трехмерная динамическая модель полей яркости фона, использующая расчетные индикатрисы многократного рассеяния излучения в облаках и изображения, полученные метеорологическими спутниками // Вопросы радиоэлектроники. Сер. РТЛ. 2004. Вып. 1. С. 101—117.

Вячеслав Филиппович Фатеев Владимир Николаевич Лагуткин Юрий Вячеславович Слынько

Сведения об авторах — д-р техн. наук, профессор; президент ОАО МАК „Вымпел“, кафедра
информационных систем, Москва; E-mail: vimpel@vimpel.ru — канд. техн. наук, доцент; ОАО МАК „Вымпел“, кафедра информа-
ционных систем, Москва; E-mail: vlagutkin@mail.ru — инженер; ОАО МАК „Вымпел“, кафедра информационных систем,
Москва; E-mail: urassl@gmail.com

Рекомендована Ученым советом ВКА им. А. Ф. Можайского

Поступила в редакцию 20.10.08 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 4