CC BY
40ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2018, том 28, № 4, c. 103-109
- РАБОТЫ С КОНФЕРЕНЦИИ -
УДК 551.508.953
© П. П. Гейко, Д. В. Петров, С. С. Смирнов
РЕАЛИЗАЦИЯ МЕТОДА ДИФФЕРЕНЦИАЛЬНОЙ ОПТИЧЕСКОЙ АБСОРБЦИОННОЙ СПЕКТРОСКОПИИ ДЛЯ ИЗМЕРЕНИЯ ВУЛКАНИЧЕСКИХ ГАЗОВЫХ ВЫБРОСОВ
Метод дифференциальной оптической абсорбционной спектроскопии является эффективным средством измерений интегральной концентрации атмосферных и примесных газов на длинных трассах. Предложен макет портативного энергонезависимого газоанализатора, который может использоваться для дистанционных трассовых измерений дегазации от вулканических выбросов (диоксида серы, сероуглерода и оксидов хлора и брома). Газоанализатор включает в себя ультрафиолетовые светодиоды, излучающие в ближней УФ-области спектра, световоды, приемо-передающий телескоп, спектрометр и систему обработки. В статье приводятся описание методики и некоторые результаты полевых измерений содержания диоксида серы и оксидов хлора и брома.
Кл. сл.: дифференциальная оптическая абсорбционная спектроскопия, ультрафиолетовые светодиоды, диоксид серы, диоксид хлора, оксид брома
ВВЕДЕНИЕ
Метод дифференциальной оптической абсорбционной спектроскопии (ДОАС) является эффективным средством для трассовых измерений атмосферных и примесных газов [1, 2]. При этом в измерительных системах, реализующих его, традиционно в качестве источника излучения используется ксеноновая лампа высокого давления, требующая высоковольтного питания, что делает такого рода газоанализаторы громоздкими и энергозависимыми [3].
Основными достоинствами метода является то, что он позволяет проводить измерения концентраций нескольких газов в режиме реального масштаба времени с высокой чувствительностью, также имеется возможность измерять средние концентрации на трассе протяженностью до нескольких километров. Кроме того, существует возможность измерять распределенные значения концентрации и даже строить 2D- и 3D-карты загрязнений, однако это требует нескольких контротражателей, установленных в требуемых участках, и автоматизированной системы сканирования приемо-пере-дающего устройства. Основными газами, выбрасываемыми из вулканического жерла, являются водяной пар, СО2, SO2, Н^ и галогениды водорода, оксиды хлора и брома. Состав выбрасываемых из вулканического жерла газов содержит информацию о магме, ее давлении и температуре. Поэтому разработка аппаратуры дистанционного газоанализа вулканических выбросов является актуальной.
Целью настоящей работы является изучение возможностей практической реализация метода ДОАС для определения концентрации некоторых оксидов галогенов с помощью разработанного макета трассового газоанализатора на основе УФ-светодиодов.
МЕТОДИКА И МАТЕМАТИЧЕСКИЙ ФОРМАЛИЗМ ДОАС
Традиционный метод дифференциального поглощения в последние годы существенно видоизменился: для идентификации и определения концентраций газовых загрязнителей вместо двух линий в и вне полосы поглощения стала возможной работа с большим количеством оцифрованных спектральных компонент. Обработка полученного спектра состоит в удалении спектра поглощения трассы и сравнении остаточного спектра с индивидуальными спектрами поглощения интересующих газов на основе процедур статистических испытаний. Подобная методика позволяет идентифицировать одновременно уже десятки газов, повысить аппаратное быстродействие, обеспечить высокую точность результатов.
Рассмотрим методические особенности измерений концентрации газов широкополосным источником излучения, какими являются светодиоды и лампы высокого давления. Поскольку эффективности рассеяния Рэлея и Ми изменяются плавно с длиной волны, их спектральные зависимости можно представить в виде полинома малого по-
рядка. Поперечное сечение поглощения газами разделяется на высокочастотную (дифференциальную) и низкочастотную части:
а, (Я) = а'. (Я) + а°(Я),
(1)
xexp
-Ь •Ха0(Я)С7 +%(Я) + £М(Я)
7=1
, (3)
имеющую смысл интенсивности в отсутствие селективного поглощения. Перепишем выражение
для /0 (Я), представив низкочастотную (медленно меняющуюся) часть в виде полинома 1-го порядка:
10(Я)=10(Я)£ Ьряр. р=0
Тогда выражение (2) перепишется в виде
(4)
I (Я) = 10(Я) X
X
exp Ь^а 7 (Я)С. +Х ЬрЯ
7=1 Р=0
где (Г0 (Я), о( (Я) — низкочастотная и высокочастотная части соответственно.
Основным выражением для любого абсорбционного спектроскопического метода является закон Бугера—Ламбета—Бера, записанный в той или иной форме. Запишем выражение для зарегистрированной детектором интенсивности
I (X):
I (Я) = 10(Я) • exp |-Ь (Я)С. х exp|-Ь • 2 о"(Я)С,. + £м(Я) + % (Я)|, (2)
где 10(Я) — исходная интенсивность; ¿м (Я), £К(Я) — коэффициенты ослабления за счет
рассеяния Ми и Рэлея соответственно, С. —
концентрация .-го газа. Коэффициенты ослабления за счет рассеяния Ми можно определить как
^м(Я) =ом(Я)СМ и Рэлея — £к(я) = °К(Я)СК соответственно.
Поскольку интенсивность рассеяния Рэлея
(~ Я 4) и Ми (~ Я ", П = 1 ^ 4) изменяются плавно с длиной волны, их спектральные зависимости можно представить в виде полинома относительно невысокого порядка. Введем интенсивность
10 (Я) = А(Я) • 10(Я) X
J
Прологарифмируем уравнение (5) и включим эффективность рассеивания в полином, получим выражение (6), зависящее от оптической толщи, полинома 1-го порядка, концентрации и длины атмосферной трассы:
D(Я) = 1п
' I (Я) л
V10(Я),
J
= -Ь^о. (Я)С. + 2 ЬрЯР, (6)
7=1 Р=0
где Ь* — коэффициенты полинома, с включенной
эффективностью рассеивания.
В общем виде оптическая толща записывается как Б(Я) = Ь • С •((Я), а с учетом выражения (1) перепишем ее в виде:
БДЯ) = Ь • С • о. (Я) = Ь • С х
х(о' (Я) + о0(Я)) = Б/ (Я) + Б 0(Я), (7)
где Б 7(Я), Б°(Я) — высокочастотная и низкочастотная части оптической толщи соответственно.
Для определения дифференциальной опти-
ческой толщи
Б' (Я)
мы использовали раз-
личные процедуры фильтрации. Эти процедуры включали приближение полиномами,
аппроксимацию сплайнами, цифровое сглаживание, аппроксимацию отрезками ряда Фурье.
Процедура фильтрации позволяет разделить спектры на дифференциальную и плавную составляющие. Дифференциальная составляющая атмосферного спектра Бэ (Як) сравнивается с дифференциальными составляющими искомых J I
газов Б(Як) = 20}ка] + 2ЬрЯР (процедура
7=1 р=0
подгонки дифференциальных частей спектров), что позволяет идентифицировать искомый газ и определить его концентрацию:
2
к=1
Бэ (Як)-2(7ка7-2 ЬрЯ
7=1 р=0
^ тт,
(8)
2
где а^ = С]Х , k — число отсчетов (определяется
числом пикселей используемого детектора, J — число измеряемых газов.
В итоге решено было остановиться на приближении полиномами, поскольку данная процедура фильтрации обеспечивает наименьшее количество коэффициентов и простоту реализации, что способствует уменьшению систематической ошибки и времени обработки экспериментальных данных. В программном обеспечении, в качестве процедуры фильтрации, реализовано приближение полиномами с выбором степени полинома (I = 3-6).
ОПИСАНИЕ УСТАНОВКИ
Ранее в наших исследованиях использовался газоанализатор, в котором источником излучения выступала мощная ксеноновая лампа высокого давления [4]. Недостатками дуговых ксеноновых ламп являются существенное снижение их характеристик после 200 ч работы и флуктуации мощности разряда, для снижения которых требуется введение обратной отрицательной связи между источником питания лампы и ее излучением. Кроме того, для поджига дугового разряда требуется высокое напряжение (десятки кВ), которое вызывает сильные электромагнитные помехи. Успехи последних лет в разработке источников излучения на базе гетероструктур послужили созданию мощных светодиодов в широком диапазоне спектра, которые являются перспективной альтернативой газоразрядным лампам при использовании в составе газоанализаторов [5]. ДОАС-газо-анализаторы обычно состоят из источника света,
спектрального прибора, приемо-передающего рефлектора и ретрорефлектора (моностатическая схема). На рис. 1. приведена принципиальная схема разработанного газоанализатора.
Пучок излучения светодиодов формируется с помощью сферического зеркала и направляется на открытую атмосферную трассу длиной от десятков до сотен метров, отражается уголковым кварцевым отражателем и поступает в приемный канал телескопа. Ввиду значительной угловой расходимости излучения светодиодов использовалась двойная коллимация выходного излучения. Ультрафиолетовые светодиоды устанавливались в мини-коллиматор прожекторного типа. Коллиматор использовался в качестве предварительного коллимационного устройства. Последующая коллимация осуществлялась сферическим зеркалом диаметром 20 см. Размер и вес установки удалось значительно уменьшить, заменив кассегреновскую схему на коаксиальную конструкцию с оптоволоконной связью. Источник излучения связан с рефлектором многожильным оптико-волоконным световодом, один конец которого помещен в фокальную плоскость сферического зеркала, выполняющего роль как передающего, так и принимающего устройства. Оптико-волоконный световод состоит из шести волокон диаметром 100 мкм для передачи и центрального — диаметром 200 мкм для приема излучения. Это позволяет связать два све-тодиода в оптический канал, предварительно разделив волокна для передачи на две группы. Предусмотрена возможность подключения третьего источника для расширения спектра излучения или для увеличения отношения сигнала к шуму.
(нсточншс света)
Спектрограф
Рис. 1. Оптическая схема трассового газоанализатора
Использование такой схемы оптоволоконной связи исключает необходимость в каких-либо вторичных зеркалах и позволяет отделить источник излучения от рефлектора, что делает возможным облегчить конструкцию телескопа в целом. Такая конструкция в последние годы стала чрезвычайно популярна при создании автономных газоанализаторов. В работе [6] подобный газоанализатор использовался для измерения реальных вулканических выбросов вулкана Ла-Фосса в Италии.
По световоду излучение попадает на входную щель монохроматора. В плоском поле изображения устанавливается линейка фотодиодов, сигнал с которой оцифровывается и поступает в компьютер для дальнейшей обработки. В качестве детектора использовался линейный датчик изображения Toshiba TCD 1205, имеющий 2048 пикселей с размером пикселя 0.014 х 0.2 мм. Длина светочувствительной зоны детектора составляет 28.6 мм, при этом регистрируемый спектральный диапазон 190-500 нм определяется используемой дифракционной решеткой. Детектор обладает свойством, исключающим перетекание зарядов из засвеченных пикселей в соседние пиксели. Минимальное время накопления 4.2 мс, максимальное время накопления 4 с, максимальным временем накопления считается время, при котором темновой сигнал достигает 10 % динамического диапазона. В качестве источника излучения апробировались УФ-све-тодиоды, излучающие в области 370, 380 и 396 нм, с мощностью излучения от 20 до 50 мВт соответственно [http://www.hamamatsu.su/]. Кроме того, для определения концентрации SO2 использовался светодиод, излучающий в области 280 нм.
РЕЗУЛЬТАТЫ ЭКСПЕРИМЕНТОВ
Проведена серия измерений концентраций газов монооксида брома (ВгО), диоксида хлора (С1О2) и серы ^О2) в южной части г. Томска. Протяженность трассы изменялась от 300 до 500 м.
Спектральные характеристики исследуемых газов являются одним из ключевых компонентов программного обеспечения метода ДОАС. Прямое использование имеющихся баз спектральных данных невозможно по следующим причинам. Конструктивные особенности конкретного дистанционного оптического анализатора накладывают отпечаток на регистрируемый спектр из-за индивидуальной аппаратной функции прибора, что приводит к изменению вида полос и линий поглощения, формы которых могут отличаться от представленных в исходной базе спектральных данных. Также не всегда условия получения (температура, давление и т. д.) спектральных данных, в частности сечений поглощения, совпадают с реальными полевыми условиями. Адаптация данных спектральных параметров атмосферных газов проводилась путем свертки спектров высокого разрешения с аппаратной функцией нашего прибора.
Расчет сечений поглощения газов, попадающих в область излучения используемых диодов, проводился с использованием аппаратной функции прибора на основе спектроскопической информации базы ШTRAN [https://www.cfa.harvard.edu/hitran/]. Спектральные структуры, полученные из измеренной оптической толщи, сравнивались с узкополосными структурами сечений поглощения искомых газов, взятыми из базы данных HITRAN
Ш 368 370 372 374 376 37S
Длина волны, нм
с О
рсз
т
и
Д
2,2 S 2.24 2,23 2.22
2 21 2,20
2.4 1.6
0,8 0.0 -0.S -1,6
N -1.1, 1-ТП . I ■ 1 _
xlOz
г—.—с—.—1—.— 1 SJ 1 ' j ■ ■ ■ - 1 ■ ■ -i
390 392 394 3% 398 400 402 404
Длина волны, нм
Рис. 2. Разделение оптической толщи.
a — оптическая толща для излучения светодиода с максимумом 370 нм; б — оптическая толща для излучения светодиода с максимумом 396 нм; 1 — оптическая толща, 2 — сглаживающий полином
а
б
и подвергнутыми процедуре свертки с аппаратной функцией спектрометра. На рис. 2 приведен пример разделения оптической толщи.
Процедура подгонки состояла из комбинации нелинейного метода Левенберга—Марквардта [7] и стандартного линейного метода наименьших квадратов.
Для повышения точности подгонки проводился сдвиг зарегистрированного спектра в пределах ± 0.3 нм с шагом спектрального разрешения спектрометра, что позволило компенсировать изменения, связанные с картированием (длина волны— пиксель спектрометра), вызванные флуктуациями температуры окружающей среды.
Рассчитанные значения концентраций составляют: диоксид хлора 7.7 ± 1.5 мкг/м3 (3.85 ± ± 0.75 ррЬ); оксид брома 23.5 ± 2 мкг/м3 (36.4 ± 3.2 ррЬ); диоксид азота 14.8 ± 1.8 мкг/м3 (7.4 ± 0.9 ррЬ). На рис. 3 приведен пример измерений концентраций диоксида серы в г. Томске.
Местоположение измерительной трассы характеризуется как фоновое для указанных газов. Мы сравнили полученные данные для монооксида брома, диоксида хлора и азота с результатами измерений газоанализатором на основе дуговой ксе-ноновой лампы высокого давления (Osram ХВО 150 W). Оказалось, что значения концентраций неплохо согласуются между собой, но при этом ошибки измерений газоанализатором на основе ксеноновой лампы больше в среднем в 1.5 раза. Это обусловлено прежде всего различием спектральных характеристик источников излучения. Увеличить отношение сигнал/шум можно использованием двух или более светодиодов с одинаковыми спектральными характеристиками, объединенных в один оптический канал. Точность измерений можно повысить путем увеличения трассы за счет возрастания оптической толщи измеряе-
мых газов. На рис. 3 приведен пример часового хода измерений концентраций диоксида серы. Минимально детектируемая концентрация диоксида серы составляет 0.1 ррЬ на трассе длиной 500 м. Одновременно проводились измерения диоксида серы на параллельной трассе (удаление 20 м) с помощью сертифицированного газоанализатора ДОАС 1-М. Результаты измерений оказались в пределах допустимых погрешностей 10 %.
ВЫВОДЫ
Методика ДОАС может быть успешно использована для измерения дегазации от вулканических выбросов, прежде всего диоксида серы, азота и оксидов хлора и брома. Благодаря достижениям последних лет в области создания ультрафиолетовых светодиодов и технологии оптических волокон нам удалось создать макет портативного энергонезависимого газоанализатора для дистанционных трассовых измерений, например для измерения дегазации от вулканических выбросов. Апробировались УФ-светодиоды, излучающие в области 280, 365, 380 и 395 нм, с мощностью излучения от 1 до 50 мВт соответственно. Проведенные расчеты и эксперименты показали возможность дистанционного детектирования двуокиси серы, оксидов хлора и брома с высокой чувствительностью (единицы и доли ррЬ) при использовании коммерчески доступных светодиодов УФ-диапазона спектра. Параллельные измерения диоксида серы с помощью сертифицированного газоанализатора ДОАС 1-М и нашего газоанализатора показали хорошую корреляцию, различие показаний не превышало 10 %. Направление дальнейших работ должно быть связано с метрологической аттестацией газоанализатора.
Рис. 3. Временной ход содержания диоксида серы на трассе длиной 500 м
СПИСОК ЛИТЕРАТУРЫ
1. Platt U. Differential optical absorbtion spectroscopy / U. Platt, J. Stutz. Berlin, Heidelberg: Springer Verlag, 2008. 593 p.
2. Stutz J., Hurlock S., Colosimo S., Tsai C., Cheung R., Fes-ta J., Pikelnaya O., Alvarez S., Flynn J., Erickson M., Olaguer E. A novel dual-LED based long-path DOAS instrument for the measurement of aromatic hydrocarbons // Atmospheric Environent. 2016. Vol. 147, no. 1. P. 121-132. Doi: 10.1016/j.atmosenv.2016.09.054.
3. Geiko P.P, Smirnov S.S., Samokhvalov I.V. Detection of concentration small gas components of atmosphere by DOAS method // Optical Memory and Neural Networks (Information Optics). 2015. Vol. 24, no. 2. P. 152-158.
4. Смирнов С.С., Гейко П.П. Многокомпонентный дистанционный газоанализ атмосферы в УФ-области спектра // Известия ВУЗов. Физика. 2015. Т. 58, № 8/3. С. 218-221.
5. Kern C., Trick S., Rippel B., Platt U. Applicability of light-emitting diodes as light sources for active differential optical absorption spectroscopy measurements // Applied Opt. 2006. Vol. 45, no. 9. Р. 2077-2068. Doi: 10.1364/AO.45.002077.
6. Vita F., Kern C., Inguaggiato S. Development of a portable active long-path differential optical absorption spec-troscopy system for volcanic gas measurements // J. Sens. Syst. 2014. Vol. 3, no. 1. Р. 355-367. Doi: 10.5194/jsss-3-355-2014.
7. 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.
Институт мониторинга климатических и экологических систем СО РАН, г. Томск
(Гейко П.П., Петров Д.В., Смирнов С.С.)
Национальный исследовательский Томский государственный университет, г. Томск
(Гейко П.П., Петров Д.В., Смирнов С.С.)
Контакты: Гейко Павел Пантелеевич, [email protected]
Материал поступил в редакцию 26.06.2018
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2018, Vol. 28, No. 4, pp. 103-109
IMPLEMENTATION OF THE METHOD OF DIFFERENTIAL OPTICAL ABSORPTION SPECTROSCOPY FOR MEASUREMENTS OF VOLCANIC GAS EMISSIONS
P. P. Geiko1,2, D. V. Petrov1,2, S. S. Smirnov1,2
1 Institute of Monitoring of Climatic and Ecological Systems SB RAS, Tomsk, Russia 2Tomsk State University, Tomsk, Russia
An active method of differential optical absorption spectroscopy (DOAS) is an effective tool for long-path measurements of atmospheric gases and impurity. A model of a portable non-volatile gas analyzer was developed, which can be used for remote trace measurements of degassing from volcanic emissions (sulfur dioxide, carbon disulfide and chlorine and bromine oxides). The gas analyzer includes ultraviolet LEDs emitting in the near UV region of the spectrum, fiber bundle, a receiving-transmitting telescope, a spectrometer and a processing system. The article describes the method and some results of field measurements of sulfur dioxide and chlorine and bromine oxides.
Keywords: differential optical absorption spectroscopy, ultraviolet LEDs, sulfur dioxide, chlorine dioxide, bromine oxide
REFERENСES
1. Platt U. Differential optical absorption spectroscopy. Eds.: U. Platt, J. Stutz. Berlin, Heidelberg, Springer Verlag, 2008. 593 p.
2. Stutz J., Hurlock S., Colosimo S., Tsai C., Cheung R., Festa J., Pikelnaya O., Alvarez S., Flynn J., Erickson M., Olaguer E. A novel dual-LED based long-path DOAS instrument for the measurement of aromatic hydrocarbons. Atmospheric Environent, 2016, vol. 147, no. 1, pp. 121132. Doi: 10.1016/j.atmosenv.2016.09.054.
3. Geiko P.P., Smirnov S.S., Samokhvalov I.V. Detection of concentration small gas components of atmosphere by DOAS method. Optical Memory and Neural Networks (Information Optics), 2015, vol. 24, no. 2, pp. 152-158. Doi: 10.3103/S1060992X15020034.
4. Smirnov S.S., Geiko P.P. [Multicomponent remote gas analysis of the atmosphere in the UF-area of a range]. Iz-vestiya VUZov. Fizika [News of higher education institu-
Contacts: Geiko Pavel Panteleevitch, [email protected]
tions. Physics], 2015, vol. 58, no. 8/3, pp. 218-221. (In Russ.).
5. Kern C., Trick S., Rippel B., Platt U. Applicability of light-emitting diodes as light sources for active differential optical absorption spectroscopy measurements. Applied Opt, 2006, vol. 45, no. 9, pp. 2077-2068.
Doi: 10.1364/AO.45.002077.
6. Vita F., Kern C., Inguaggiato S. Development of a portable active long-path differential optical absorption spec-troscopy system for volcanic gas measurements. J. Sens. Syst, 2014, vol. 3, no. 1, pp. 355-367. Doi: 10.5194/jsss-3-355-2014.
7. Marquardt D.W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Indust. Appl. Math., 1963, vol. 11, no. 2, pp. 431-441. Doi: 10.1137/0111030.
Article received in edition 26.06.2018
