CC BY
45
УДК 621.396.
П.П. Гейко, С.С. Смирнов
Анализ содержания оксидов галогенов в атмосфере методом дифференциальной оптической спектроскопии в УФ-области спектра
Метод дифференциальной оптической абсорбционной спектроскопии является эффективным средством для трассовых измерений концентраций атмосферных и примесных газов. В газоанализаторах на его основе традиционно в качестве источника излучения используется ксеноновая лампа высокого давления, требующая высоковольтного питания, что делает такого рода газоанализаторы громоздкими и энергозависимыми. Благодаря последним достижениям в области создания ультрафиолетовых светодиодов и технологии оптических волокон удалось создать макет портативного энергонезависимого газоанализатора. Газоанализатор включает в себя два ультрафиолетовых светодиода, излучающих в ближней УФ-области спектра, систему световодов, приемопередающий телескоп, спектрометр и систему обработки. В статье приводятся описание методики и результаты измерений оксидов хлора и брома с помощью газоанализатора на основе ультрафиолетовых LED-диодов. Ключевые слова: дифференциальная оптическая спектроскопия, сечение поглощения, ультрафиолетовые све-тодиоды, оксид хлора, оксид брома. doi: 10.21293/1818-0442-2017-20-4-33-38
Интерес к дистанционному контролю содержания оксидов галогенов в воздухе обусловлен многочисленными приложениями, в частности, быстро развивающейся в последние годы химией энергоемких соединений. Энергоемкие вещества могут отдавать энергию практически мгновенно, например в камере сгорания реактивного двигателя, или с регулируемой скоростью, например в химических источниках тока. Эти соединения играют важную роль в вулканологии и экологии. Основными газами, выбрасываемыми из вулканического жерла, являются водяной пар, СО2, 802, И28 и галогениды водорода, оксиды хлора и брома. Состав выбрасываемых из вулканического жерла газов содержит информацию о магме, ее давлении и температуре.
Среди разработок аппаратуры дистанционного газоанализа лидирующее положение занимают оптические методы контроля. Одним из наиболее перспективных оптических методов является метод дифференциальной оптической абсорбционной спектроскопии (ДОАС) [1].
Перечислим основные достоинства ДОАС-метода. Во-первых, он позволяет проводить измерения концентраций большого числа газов в режиме реального масштаба времени без предварительных и пост-процессов. Здесь под пост-процессами понимаются предварительные операции с исследуемыми пробами. Во-вторых, ДОАС-газоанализаторы способны измерять средние концентрации на атмосферной трасе протяженностью до нескольких километров и в принципе позволяют создать трехмерную карту загрязнения. Наконец, метод ДОАС разделяет преимущества большинства других спектроскопических методов, включая высокую чувствительность вплоть до нескольких рр1 (триллионных долей). Кроме того, существует потенциальная возможность измерять распределенные значения концентрации и даже строить 2Б- и 3Б-карты загрязнений, однако это потребует нескольких контр-
отражателей, установленных в требуемых участках, автоматизированной системы сканирования приемопередающего устройства.
Цель настоящей работы - изучение возможно -стей практической реализация метода ДОАС для определения концентрации некоторых оксидов галогенов с помощью разработанного трассового газоанализатора на основе УФ-светодиодов.
Описание методики. Рассмотрим методические особенности измерений концентрации газов широкополосным источником излучения. С развитием вычислительной техники традиционный метод дифференциального поглощения существенно видоизменился, для идентификации и определения концентраций газовых загрязнителей вместо двух линий «в» и «вне» полосы поглощения стала возможной работа с большим количеством оцифрованных спектральных компонент. Обработка полученного спектра состоит в удалении спектра поглощения трассы и сравнении остаточного спектра с индивидуальными спектрами поглощения интересующих газов на основе процедур статистических испытаний. Подобная методика позволяет идентифицировать одновременно уже десятки и более веществ, повысить аппаратное быстродействие, обеспечить высокую точность результатов.
Математически интенсивность принятого излучения на заданной длине волны описывается уравнением [1]
1(1) = Л(1) /0(Х)х
(
м
(1)
хехр\-Ь £ст](Х)Су +стМе(Х)СМе
I
где Л(Х) - коэффициент пропускания оптической системы; 1о (X) - исходная интенсивность; ст j (X) -сечение селективного поглощения газами ст^{е(Х); стяау (X) - сечения рассеяния Ми и Рэлея соответст-
34
венно; Су - концентрация ¡-го газа; 3 - число измеряемых газов; См^, С^ау - концентрации аэрозольных и молекулярных рассеивателей. Введем коэффициенты ослабления за счет рассеяния Ми
8Ме(Х> = стМ1е(Х>Ске и Рэлея БЯау(^) = стЯау(^)СЯау
соответственно.
Эффективность рассеяния Рэлея и Ми изменяется плавно с длиной волны, их спектральные зависимости можно представить в виде полинома малого порядка. Поперечное сечение поглощения газами может быть разделено на высокочастотную (дифференциальную) и низкочастотную части
а} (X) = а} (X) + ст} о(Х), (2)
где ауо(Х),ау (X) - низкочастотная и высокочастотная части соответственно.
Подставим выражение (2) в формулу (1), получим
I (X) = А(Х) ■ 10 (X) • ехр{-1 X а' у (Х)Су }х у=1
хехр{-Ь( ¿а yо(X)Cy +eRay(X) + Sмie(X))}. (3) у=1
Введем интенсивность
ТО = A(X)■ ^ехрЬД ¿аjо(XX)Cj +eMe(X))},
у=1
имеющую смысл интенсивности в отсутствие селективного поглощения.
Перепишем выражение для 1о( X), представив низкочастотную (медленно меняющуюся) часть в виде полинома 1-го порядка.
I
Го(X),=Ы>^ЬрXр . (4)
р
Тогда, выражение (1) перепишется в виде
3 I
I (X)=1о (X) ■ ехр{-£ X а / (X)C^ + ¿Ьр X р }. (5)
у=1 р
Прологарифмируем уравнение (5) и получим выражение для оптической толщи -О(Х) (6), зависящее от дифференциальных сечений поглощения, концентраций измеряемых газов, длины атмосферной трассы и полинома:
ЭД = 1п (^О = -1 ■ £ а у (X)Cу■ + ¿Ьр Xр . (6)
^ М^») у=1 р
Уравнение (6) можно представить в дискретном виде:
3 I
^) = У^))^)Су +^р . (7)
у=1 р
Здесь к = 1, п; п - число дискретных отсчетов сигнала (определяется числом пикселей используемого детектора). По сути (7) представляет собой систему линейных уравнений, которую можно переписать в матричном виде и использовать тот или иной чис-
ленный метод, обеспечивающий лучшую устойчивость.
Экспериментальная установка
ДОАС-газоанализаторы обычно состоят из источника света, спектрального прибора, приемопередающего рефлектора и ретрорефлектора (моностатическая схема) [2]. Ранее в наших исследованиях использовался газоанализатор, в котором источником излучения выступала мощная ксеноновая лампа высокого давления [3, 4]. Недостатками дуговых ксеноновых ламп являются существенное снижение их характеристик после 2оо ч работы и флуктуации мощности разряда, для снижения которых требуется введение обратной отрицательной связи между источником питания лампы и ее излучением. Кроме того, для поджига дугового разряда требуется высокое напряжение (десятки кВ), которое вызывает сильные электромагнитные помехи.
Успехи в разработке источников излучения на базе гетероструктур послужили созданию мощных светодиодов в широком диапазоне спектра, которые являются перспективной альтернативой газоразрядным лампам при использовании в составе газоанализаторов [5, 6]. На рис. 1. приведена принципиальная схема разработанного газоанализатора.
Рис. 1. Принципиальная схема ДОАС-газоанализатора
Световой пучок излучения светодиодов формируется с помощью сферического зеркала и направляется на открытую атмосферную трассу длиной от десятков до сотен метров, отражается уголковым кварцевым отражателем и поступает в приемный канал телескопа. Размер и вес установки удалось значительно уменьшить, заменив кассегреновскую схему на коаксиальную конструкцию с оптоволоконной связью. Источник излучения связан с рефлектором многожильным оптико-волоконным световодом (ОВС), один конец которого помещен в фокальную плоскость сферического зеркала, выполняющего роль как передающего, так и принимающего устройства. Оптико-волоконный световод состоит из шести волокон диаметром 100 мкм для передачи и центрального диаметром 200 мкм для приема излучения. Это позволяет связать два светодиода в оптический канал, предварительно разделив волокна для передачи на две группы. Предусмотрена возможность подключения третьего источника для расширения спектра излучения или для увеличения
отношения сигнала к шуму. Использование такой схемы оптоволоконной связи исключает необходимость в каких-либо вторичных зеркалах и позволяет отделить источник излучения от рефлектора, что делает возможным облегчить конструкцию телескопа в целом. Такая конструкция в последние годы стала чрезвычайно популярна при создании автономных газоанализаторов [7-9]. В работе [7] подобный газоанализатор использовался для измерения реальных вулканических выбросов вулкана Ла-Фосса в Италии.
По световоду излучение попадает на входную щель монохроматора. В плоском поле изображения устанавливается линейка фотодиодов, сигнал с которой оцифровывается и поступает в компьютер для дальнейшей обработки. В качестве детектора использовался линейный датчик изображения Toshiba TCD 1205, имеющий 2048 пикселей с размером пикселя 0,014^0,2 мм2. Длина светочувствительной зоны детектора составляет 28,6 мм, при этом регистрируемый спектральный диапазон 190-500 нм определяется используемой дифракционной решеткой. Детектор обладает свойством, исключающим перетекание зарядов из засвеченных пикселей в соседние пиксели. Минимальное время накопления 4,2 мс, максимальное время накопления 4 с, максимальным временем накопления считается время, при котором темновой сигнал достигает 10% динамического диапазона.
Спектроскопические данные. Спектральные характеристики исследуемых газов являются одними из ключевых компонентов программного обеспечения метода ДОАС. Прямое использование имеющихся баз спектральных данных невозможно по следующим причинам. Конструктивные особенности конкретного дистанционного оптического анализатора накладывают отпечаток на регистрируемый спектр из-за индивидуальной аппаратной функции прибора, что приводит к изменению вида полос и линий поглощения, формы которых могут отличаться от представленных в исходной базе спектральных данных. Также не всегда условия получения (температура, давление и т.д.) спектральных данных, в частности сечений поглощения, совпадают с реальными полевыми условиями. Адаптация данных спектральных параметров атмосферных газов проводилась путем свертки спектров высокого разрешения с аппаратной функции нашего прибора.
В качестве источника излучения апробировались УФ-светодиоды, излучающие в области 365 и 395 нм, с мощностью излучения 20 и 30 мВт соответственно [http://www.hamamatsu.su/]. Относительно малый интервал излучения длин волн светодио-дов ограничивает измерения до трех газов: диоксида азота (NO2), монооксида брома (BrO) и диоксида хлора (ClO2) (рис. 2).
На рис. 2 представлены спектральные характеристики доступных нам светодиодов. Ввиду большой угловой расходимости излучения светодиодов было решено использовать двойную коллимацию
выходного излучения. Ультрафиолетовые светодио-ды устанавливались в мини-коллиматор прожекторного типа. Коллиматор использовался в качестве предварительного коллимационного устройства. Дальнейшая коллимация излучения осуществлялась основным зеркалом с диаметром 0,2 м и фокусным расстоянием 1 м.
УФ-светодиод (365 нм) 2 - УФ-светодиод (395 нм)
350 360 370 380 390 400 410 420 430 440 450 Длина волны, нм
Рис. 2. Относительное спектральное распределение интенсивностей УФ-светодиодов
Дифференциальные сечения исследуемых газов поглощения газов представлены на рис. 3. Расчет сечений поглощения газов, попадающих в область излучения используемых диодов, проводился с использованием аппаратной функции прибора на основе спектроскопической информации базы ИГТОАК [https://www.cfa.harvard.edu/hitran/].
Г ) Г. \ J ClO2
xN
-0,1
355 360 365 370 375 380 385 390 395 400 405 410 415 420 425
Длина волны, нм
Рис. 3. Дифференциальные сечения поглощения
диоксида хлора, оксида брома и диоксида азота
Видно, что интересующие нас соединения имеют выраженные спектральные особенности в области излучения используемых светодиодов и, следовательно, могут быть зарегистрированы при концентрациях, превышающих некий порог, который следует определить.
Минимально детектируемые концентрации. Для оценки возможностей газоанализатора проводи-
0
лись оценки минимально детектируемых концентраций (МДК) в режиме реального времени.
Ключевой характеристикой, определяющей чувствительность метода ДОАС, является минимально детектируемая толща -Оо(Х), определяемая статистикой фотоэлектронов (дробовой шум) 51, тепловым и фотонным шумами, включая прочие шумы 52:
где 51 = N "1/2 чения для
52 = В *Ш-4
А)^) = 4 + ^, (8)
а N - число фотонов, типичные зна-твердотельных фотоприемников Число фотонов легко оценить:
Р■ гт •-п^
N =-
ко
(9)
здесь Р - мощность излучения; гт - время накопления; ^ - квантовая эффективность; X - длина волны; к - постоянная Планка; о - скорость света. Полагая X = 37о нм; ^ = о,5; гт = 15о с и варьируя Р в пределах от 1о-6 до 1о2 Вт, можно увидеть, что дробовой шум не является определяющим.
С учетом этого получаем оценку минимальной детектируемой толщи:
-1 + В2 =
ко
-+В2«5 ■Ю
-4
(1о)
\р ■ т
Проведем оценки МДК с учетом уравнения (7). Рассчитанная МДК для трассы длиной 5оо м составляет о,о9 мкг/м3 для С102 (диоксид хлора); о,о7 для ВгО (оксид брома); 5,4 для М02 (диоксид азота) соответственно. Также проведен расчет МДК и для некоторых других соединений. Так, для СбИ60 (фенол) МДК составила о,7; для С9И12 (триметилбен-зол) - 2о; для С9И130 (триметилфенол) - 3о; для С8И80 (толуилальдегид) - 4; для С8И1о0 (диметил-фенол) - 7; для С7И80 (крезол) - 4; для С8И1о (ксилол) - 2о соответственно.
Экспериментальные результаты. Проведена серия измерений концентраций газов монооксида брома (Вг0), диоксида хлора (С102) и азота (N0^ в южной части г. Томска. Протяженность трассы составила 5оо м (25о м до ретрорефлекора).
Узкополосные спектральные структуры, полученные из измеренной оптической толщи, сравнивались с узкополосными структурами сечений поглощения искомых газов, взятыми из базы данных HITRAN и подвергнутыми процедуре сверки с аппаратной функцией спектрометра.
Процедура подгонки состояла из комбинации нелинейного метода Левенберга-Марквардта и стандартного линейного метода наименьших квадратов [1о, 11].
2
(11)
^ Ш1П,
п 3 1
X ^к ) - X ау аук -¿ЬрXр к=1 ]=1 р
где а у = СуЬ .
При расчетах использовались алгоритм и про грамма определения концентраций многокомпо нентных газовых смесей [12].
Для повышения точности подгонки проводился сдвиг зарегистрированного спектра в пределах ±о,3 нм с шагом спектрального разрешения спектрометра, что позволило компенсировать изменения, связанные с картированием длины волны - пиксели спектрометра, вызванные флуктуациями температуры окружающей среды.
Оценка зарегистрированных спектров поглощения по остаточным структурам и рассчитанная гипотетическая оптическая толща, обусловленная поглощением трех газов (монооксида брома, диоксида хлора и азота), представлена на рис. 4.
<и л и
1,о
о,5 о,о -о,5 -1,о 1,5 1,о о,5 о,о -о,5 -1,о
А/1 <д\ | '/{ Ч ;/ ^Чл \\ 1 \ л 1Лл , л' N. /Ч м С = 7,7 мкг/м3 С102 Г \ 1 1 ! ч>\ 1 V / М 1* ч.Ч V
хЮ"2
/Т\ /-Ц, . / ч 1 / ч г / \ /7\А 1' [ -'Г^Х. и \ л /У * / \Ч*~' \А / \ ! С = 14,8 мкг/м3 ( у ■
1,5 "1
1,о Вг0
о,5 Н --.Л
о,о / К 'г' \ ч/ VV /
-о,5 К / УЛ / /
-1,о С = 23,5 мкг/м3 уу
365 366 367 368 369 37о 371 372 373 374 375 376 377 378
Длина волны, нм Рис. 4. Оценка измеренных спектров поглощения для С102, N02 и Вг0
Суть его состоит в том, что из реальной дифференциальной оптической толщи, полученной с трассы (сплошная кривая на верхнем рис. 4, с выраженными пиками на 368,6 и 377,5 нм варьированием концентрации С102, методом наименьших квадратов подбирается дифференциальная оптическая толща, обусловленная поглощением С102 (пунктирная кривая на верхнем рис. 4).
Спектры поглощения других более 4о газов предварительно анализируются на предмет поглощения в этих областях. Затем из исходной оптической толщи убирается оптическая толща, обусловленная поглощением С102, получается сплошная линия на среднем рис. 4. Процедура последовательно повторяется для N02 и Вг0. В рассматриваемой спектральной области наибольший вклад в поглощение вносят три газа: С102, N0^ Вг0, но имеется
поглощение H2CO (формальдегид), и HNO2 (сечения поглощения порядка 10-22 см2) и димеров молекулы кислорода О4 (10-46 см2). Возможно, неучет этого поглощения и рассеянного солнечного света приводит к значительным остаткам. Мы рассчитывали погрешности не на всем спектральном интервале 365-378 нм, приведенном на рис. 4, а на несколько меньших интервалах. К примеру, для BrO интервал выбирался от 366,96 до 377,75 нм.
Рассчитанные значения концентраций при использовании светодиода с максимумом 371 нм составляют: диоксид хлора 7,7±1,5 мкг/м3 (3,85±0,75 ppb); оксид брома 23,5±2 мкг/м3 (36,4±3,2 ppb); диоксид азота 14,8±1,8 мкг/м3 (7,4±0,9 ppb).
Следует отметить, что местоположение измерительной трассы характеризуется как фоновое для указанных газов. Полученные значения концентраций для диоксида азота (NO2) хорошо согласуются с данными, полученными со станции «Фоновой» (60 км западнее г. Томска), где проводятся почасовые измерения диоксида азота. Мы сравнили полученные данные для монооксида брома, диоксида хлора и азота с результатами измерений газоанализатором на основе дуговой ксеноновой лампы высокого давления (Osram XBO 150 W). Оказалось, что значения концентраций неплохо согласуются между собой (различия не более 20%), но при этом ошибки измерений газоанализатором на основе ксеноновой лампы больше в среднем в 1,5 раза. Это обусловлено, прежде всего, различием спектральных характеристик источников излучения. Увеличить отношение сигнал/шум можно использованием двух или более светодиодов с одинаковыми спектральными характеристиками, объединенных в один оптический канал. Точность измерений можно повысить путем увеличения трассы за счет возрастания оптической толщи измеряемых газов.
Заключение
В статье приводятся описание методики и результаты расчетов минимально детектируемых концентраций оксидов хлора и брома с помощью газоанализатора на основе ультрафиолетовых светодио-дов. Создан и протестирован макет портативного энергонезависимого дистанционного трассового газоанализатора. Проведенные расчеты и эксперименты показали возможность дистанционного детектирования оксидов хлора и брома с высокой чувствительностью при использовании коммерчески доступных светодиодов УФ-диапазона спектра. Направление дальнейших работ должно быть связано с метрологической аттестацией газоанализатора.
Литература
1. Platt U. Differential optical absorbtion spectroscopy / U. Platt, J. Stutz. - Berlin; Heidelberg: Springer Verlag, 2008. - 593 p.
2. Lee J.S. Development of a Differential Optical Absorption Spectroscopy (DOAS) System for the Detection of Atmospheric Trace Gas Species: NO2, SO2, and O3 / J.S. Lee, B.J. Kuk, Y.J. Kim // J. of the Korean Phys. Soc. - 2002. -Vol. 41, No. 5. - P. 693-698.
3. Гейко П.П. Многокомпонентный дистанционный газоанализ атмосферы в УФ-области спектра / П.П. Гейко, С.С. Смирнов, И.В. Самохвалов // Изв. вузов. Физика. -
2015. - Т. 58, № 8/3. - С. 218-221.
4. Geiko P.P. Detection of concentration small gas components of atmosphere by DOAS method / P.P. Geiko, S.S. Smirnov, I.V. Samokhvalov // Optical Memory and Neural Networks (Information Optics). - 2015. - Vol. 24, No. 2. -P. 152-158.
5. Applicability of light-emitting diodes as light sources for active differential optical absorption spectroscopy measurements / C. Kern, S. Trick, B. Rippel, U. Platt // Applied Opt. - 2006. - Vol. 45, No. 9. - P. 2077-2068.
6. Катаев М.Ю. Анализ содержания метана атмосферы с помощью инфракрасных LED-диодов // Доклады ТУСУРа. - 2015. - Т. 38, № 4. - С. 155-160.
7. Vita F. Development of a portable active long-path differential optical absorption spectroscopy system for volcanic gas measurements / F. Vita, C. Kern, S. Inguaggiato // J. Sens. Syst. - 2014. - Vol. 3, No. 1. - P. 355-367.
8. An instrument for measurements of BrO with LED-based Cavity-Enhanced Differential Optical Absorption Spec-troscopy / D.J. Hoch, J. Buxmann, H. Sihler et al. // Atmos. Meas. Tech. - 2014. - No. 7. - P. 199-214.
9. A novel dual-LED based long-path DOAS instrument for the measurement of aromatic hydrocarbons / J. Stutz, S. Hurlock, S. Colosimo et al. // Atmospheric Environment. -
2016. - Vol. 147, No. 1. - P. 121-132.
10. Marquardt D.W. An algorithm for least-squares estimation of nonlinear parameters // J. Soc. Indust. Appl. Math. -1963. - Vol. 11, No. 2. - P. 431-441.
11. Stutz J. Numerical analysis and estimation of the statistical error of differential optical absorption spectroscopy measurements with least-squares methods / J. Stutz, U. Platt // Appl. Optics. - 1996. - Vol. 35. - P. 6041-6053.
12. Свидетельство о государственной регистрации программы для ЭВМ №2016611377. Программа определения концентраций многокомпонентных газовых смесей в приземной атмосфере по их селективному поглощению излучения в УФ-области спектра / С.С. Смирнов, И.Д. Брюханов, Е.В. Ни. - Заявка №2015662194. Дата поступления 14 декабря 2014 г. Зарегистрировано в Реестре программ для ЭВМ 01 февраля 2016 г.
Гейко Павел Пантелеевич
Д-р физ.-мат. наук, вед. науч. сотрудник Института мониторинга
климатических и экологических систем СО РАН, профессор каф. оптико-электронных систем и дистанционного зондирования (ОЭСиДЗ) Национального исследовательского Томского государственного университета Тел.: 8-905-991-9762, (382-2) 41-38-87 Эл. почта: ppg@imces.ru
Смирнов Сергей Сергеевич
Аспирант каф. ОЭСиДЗ
Тел.: 8-952-895-17-30, (382-2) 41-39-84
Эл. почта: sssmirnov@sibmail.com
Geiko P.P., Smimov S.S.
Atmosphere halogen oxides concentration retrieving by differential optical spectroscopy in the UV region of the spectrum
The method of differential optical absorption spectroscopy is an effective tool for open-path measurements of atmospheric gases and impurities. The gas analyzers based on this method traditionally include a high pressure xenon lamp as the radiation source, that require high voltage, which makes gas analyzers bulky and power-inefficient. Recent developments in
fiber-coupling telescope technology and the availability of ultraviolet light emitting diodes have now allowed us to construct a portable, long path gas analyzer. The gas analyzer comprises two UV-LEDs emitting in the near UV-region of the spectrum, system of optical fibers, receiving-transmitting telescope, spectrometer and treatment system. The article presents a description of the methodology and results of the measurement of chlorine and bromine oxides with gas analyzer based on UV-LEDs.
Keywords: differential optical spectroscopy, absorption cross-section ultraviolet LEDs, chlorine dioxide, bromine oxide.
