CC BY
50
Аналіз теоретичних основ лазерних вимірювань складу повітряного середовища та можливості його застосування в надзвичайних ситуаціях (НС). Запропоновано метод розв’язання задачі оперативного моніторингу складу газів домішок і аерозолів в зоні НС застосуванням єдиного комплексу лазерної вимірювальної апаратури
Ключові слова: лазер, атмосфера, оперативний моніторинг
□-----------------------------------□
Анализ теоретических основ лазерных измерений состава воздушной среды и возможности его применения в чрезвычайных ситуациях (ЧС). Предложен метод решения задачи оперативного мониторинга состава газов примесей и аэрозолей в зоне ЧС применением единого комплекса лазерной измерительной аппаратуры
Ключевые слова: лазер, атмосфера, оперативный мониторинг
□-----------------------------------□
Analysis of the theoretical foundations of laser measurements of air pollution and its application in emergency situations (ES). We propose a method for solving the real-time monitoring of trace gases and aerosols in the emergency area using a single set of laser measuring devices
Key words: laser, atmosphere, real-time monitoring
УДК 621.373:772.99
ВОЗМОЖНОСТИ ПРИМЕНЕНИЯ ЛАЗЕРНЫХ ИССЛЕДОВАНИЙ АТМОСФЕРЫ ЗОНЫ ЧРЕЗВЫЧАЙНОЙ СИТУАЦИИ
Л.О. Черногор
Доктор физико-математических наук, профессор
Факультет радиофизики Харьковский национальный университет им. В.Н. Каразина пл. Свободы, 4, г. Харьков, 61077 Контактный тел.: (057) 707-55-61 E-mail: Leonid.F.Chernogor@univer.kharkov.ua
А.С. Рашкевич
Адъюнкт
Кафедра охраны труда и техногенно-экологической
безопасности
Национальный университет гражданской защиты
Украины
ул. Чернышевского, 94, г. Харьков, 61023 Контактный тел.: (057) 707-34-57 E-mail: acmilan 1899@mail.ru
Постановка проблемы
Стремительное развитие всех отраслей промышленности, энергетики, транспорта, увеличение численности населения, урбанизация и химизация всех сред деятельности человека приводят к нарушению и загрязнению биосферы, ее отдельных компонентов. Экологическая ситуация, сложившаяся в ряде промышленных центров, в районах добычи и переработки минерального сырья, строительства и эксплуатации промышленных объектов, часто близка к критической [1].
Осознание глобальной экологической катастрофы заставляет мировое сообщество искать пути выхода из кризисной ситуации. Вывод о необходимости перехода цивилизации к экологически сбалансированному развитию имеет непосредственное отношение к опасности, угрожающей человеческой цивилизации. Идея глобального мониторинга окружающей человека природной среды была выдвинута в 1972 г. на Стокгольмской конференции и нашла отклик в документах Конференции ООН в Рио-де-Жанейро (1992) [2,3].
Основной задачей экологического мониторинга является изучение изменений природной среды, возникающих в результате воздействия на нее человека,
получение как качественных, так и количественных характеристик изменений в природной среде.
Анализ последних исследований и публикаций
Наиболее перспективным методом диагностики загрязнений воздуха является дистанционные методы зондирования.
Среди дистанционных методов особое место занимают лазерные методы [4]. Одним из самых эффективных средств дистанционного мониторинга степени загрязнения окружающей среды являются лазерные мониторы, называемые также лидарами (по аналогии с радарами).
Основы лазерного мониторинга окружающей среды заложены в работах [4-6]. Их обобщение проведено в работах [7,8].
Суть методики кратко состоит в следующем.
Загрязненная атмосфера содержит не присущие ей газы (двуокись серы - SO2, окислы азота - N0 и N02, углеводороды - НХСХ и другие), продукты их реакции типа кислот и окислителей, а также твердые пылинки с размерами от 10-8 до 10-3 м (аэрозоли).
Молекулы газов - загрязнителей и аэрозолей поглощают и рассеивают лазерное излучение. На этом
базируются методы лазерного мониторинга загрязнений в атмосфере.
Цель экологического мониторинга - информационное обеспечение управления природоохранной деятельностью и экологической безопасностью.
Для достижения поставленной цели необходимо дать ответы на следующие вопросы:
• каково состояние природной среды в рассматриваемый отрезок времени в сравнении с предшествующим техногенезу состоянием (в относительной или абсолютной форме) и какие изменения (положительные, отрицательные) ожидаются в природной среде в прогнозируемый отрезок времени;
• в чем причины происшедших изменений и возможных изменений в будущем (в том числе нежелательных, губительных, критических) и что явилось, является или будет являться источником этих изменений (как правило, вредных техногенных воздействий);
• какие воздействия на данную локальную природную среду, определяемые исходя из выработанной для данного случая критериальной основы оценок функции полезности - вредности, являются вредными (нежелательными или недопустимыми);
• какой уровень техногенных воздействий, в том числе в совокупности с естественными или стихийными процессами и воздействиями, происходящими в рассматриваемой природной среде, является допустимым для природной среды и отдельных ее компонентов или комплексов (ценозов) и какие резервы имеются у природной среды для саморегенерации состояния, адекватного исходному, принятому за состояние экологического баланса;
• какой уровень техногенных воздействий на природную среду, отдельные ее компоненты и комплексы является недопустимым или критическим, после которого восстановление природной среды до уровня экологического баланса является неосуществимым.
Прикладная задачи и ее решение
Состав выбросов вредных веществ зависит от материалов, находящихся в зоне взрыва и горения. Этот состав априори неизвестен.
При горении нефти, бензина, природного газа основной состав продуктов горения хорошо известен. В него входит СО, СО2, Н20, S0, S02 и др.
Пожары сопровождаются горением сопутствующих материалов. Состав выбросов при этом может быть самым различным.
Одновременно с выбросами газов в зонах горения образуется аэрозоли, т.е. твердые несгоревшие частицы вещества с характерным размером от 10-8 до 10-3 м. Их состав и концентрация также априори неизвестны.
Несвойственные атмосфере вредные выбросы газов и аэрозолей угрожают жизни и здоровью людей: работникам производств, где возникла зона ЧС, жителям близлежащих районов, а также ликвидаторам ЧС.
Для уменьшения риска и угрозы людям требуется оперативный и достаточно точный мониторинг загрязнений в зоне ЧС.
В данной работе будут рассмотрены теоретические основы лазерного мониторинга атмосферы зоны чрезвычайной ситуации.
Основное уравнение лазерной локации. Исходные соотношения
Основным уравнением лазерной локации будем называть зависимость мощности отраженного (рассеянного) объектом зондирования излучения на заданной длине волны от расстояния до объекта и его параметров [9].
Уравнение может быть приведено к такому виду:
Р = ^ Є-2
5 Я2
(1)
где Р5 - мощность отраженного (рассеянного) излучения,
Р - мощность лазера,
S - эффективная площадь (апертура) приемного элемента,
П - его эффективность,
в - объемный коэффициент обратного рассеяния падающего излучения,
к
Г = |а(Я) dR - интегральный коэффициент погло-
0
щения лазерного излучения,
а - объемный коэффициент поглощения лазерного излучения,
К -расстояние до объекта зондирования.
В случае импульсного зондирования:
К = ст/2 ,
где с - скорость света,
т - задержка импульса при его движении от передающего объекта и от него до приемного устройства.
При использовании в качестве отражателя зондирующего пучка уголкового отражателя уравнение (1) модифицируется к виду:
р = РЗдС 2г(Ко)
п2 С ’
Я;
(2)
где С - аппаратная константа,
К0 - длина трассы,
Г(К0) - интегральный коэффициент поглощения при прохождении излучением всей трассы.
В уравнениях (1) и (2) зависит от а , (3 в свою очередь:
где ат и аа - коэффициенты, описывающие молекулярное и аэрозольное ослабления.
Ослабление лазерного излучения вследствие молекулярного рассеяния
Для ат имеем следующее соотношение:
а = ат +«а
€
= 8лвт’
(4)
где Рт - коэффициент молекулярного рассеяния. Он дается формулой Рэлея [10]:
В = 8 П(п2 -1)2 6+38 - 3 Ж4 6 - 78’
(5)
где п - показатель преломления воздуха,
8 - степень деполяризации рассеянного излучения. Для его описания преобразуем множитель в (5):
q = 0.585Я-3,
где Ят измеряется в км.
Для Ят , изменяющихся в пределах 6-50 км q = 1.3. При Ят > 50 км q = 1.6 [8].
Значение Ят для различных условий видимости приведен в табл. 2. Результаты расчета q, а и Га (К- ) для различных Ят приведена в табл. 3. Из табл. 3 видно, значение аа существенно уменьшается при улучшениивидимости,втожевремязначения Га (Ят) при увеличивающихся Ят остаются неизменными и равными 3.910-3. Иначе говоря, при Я = Ят поглощается около 0.78% мощности видимого излучения.
л = 6+38 = 1+38/ 6 8 = 6 - 78 = 1 - 78/6 .
Поскольку значение 8 при однократном молекулярном рассеянии принимается в пределах 0 - 9,7% (табл. 1). Тогда с учетом того, что 8 << 1
38 78 5
Л8 = 1 + — + — = 1 + -8 . 8 6 6 3
Значение Л8 отличается от 1 на 58/3 , т.е. не более, чем на 16%. Поэтому часто полагают Л8 = 1 Из (4) и (5) следует, что:
2 (п2 -1)2 6 + 38
ат = п--------г---------.
т Ж4 6 - 78
(5а)
При п2 -1 << 1 или то же самое при п = 1 представим в скобках в следующем виде
п2 -1 = (п- 1)(п +1) = 2(п-1).
Тогда из (5а) имеем:
2 (п2 -1)2
ат “ 4я-----.
т т4
Ослабление лазерного излучения вследствие аэрозольного рассеяния
Выражение для аа следует из теории аэрозольного рассеяния (теории Ми) [10, 11].
Однако более удобное для описания аэрозольного ослабления является эмпирическое соотношение [8]:
Таблица 1
Степень деполяризации основных атмосферных газов и примесей [8]
Газ 8 , %
Воздух 3.5
н2о 2
N2 3.6
02 6.5
со2 9.7
S02 3.1
H2S 0.3
СН4 0
С12 4.1
НС12 0.7
со 1.3
СпН2п+2 0
Таблица 2
Значения метеорологической видимости при длине волны
Х = 5.5 10-7 м [8]
Условия видимости К-
Плотный туман 3-5 м
Густой туман 5-50 м
Средний туман 50-500 м
Легкий туман 0,5-1 км
Слабый туман 1-2 км
Дымка 2-4 км
Легкая дымка 4-10 км
Ясно 10-20 км
3.912 ( X
Г-М'\
1,0.55 )
а Я- 10
где аа - измеряется в км-1,
X - измеряется в мкм,
Ят - метеорологическая дальность видимости в км при Х = 0.55 мкм,
q - показатель степени, который, вообще говоря, зависит от Ят.
При этом
Я
Г а =|а^Я.
0
Так, при Ят < 6 км [8]:
Таблица 3
Рассчитанные значения q , аа и а (Ят) при различных Ят
и Х = 5.5 10-7 м
Условия в зоне ЧС Я— , м q а , м-1 а Га (Я-)
Видимость очень плохая 1 0.059 3.9 10-3 3.9 10-3
Видимость плохая 5 0.100 0.810-3 3.910-3
Видимость средняя 10 0.126 3.9 10-4 3.9 10-3
Видимость удовлетворительная 30 0.182 1.310-4 3.910-3
Видимость хорошая 50 0.216 0.8 10-4 3.9 10-3
Видимость отличная 1000 0.585 3.910-6 3.910-3
Проведем расчеты аа и Га (Ят) для С02- и СО- лазеров, которые планируется использовать для экспериментальных исследований. Как известно, для них диапазон длин волн составляет 11.7-13 мкм и 5-6 мкм соответственно [8].
Результаты расчетов приведены в табл. 4. Считалось, что для этих лазеров в среднем Х1 = 12.3 мкм и Х2 = 5.5 мкм.
Таблица 4
Значения аа и Га (Ят) для ОД2- и ОТ- лазеров
Таблица 5
Результаты расчета максимальных размеров зоны ЧС для ОД2- и ОД- лазеров
Условия в зоне ЧС (Х1), км (Х2), км
Видимость очень плохая 0.354 0.338
Видимость плохая 1.75 1.61
Видимость средняя 3.79 3.42
Видимость удовлетворительная 13.5 11.65
Видимость хорошая 25 21.0
Видимость отличная 1560 1000
Условия в зоне ЧС Ят, м q «а (Х1) , м-1 «а (Х2 ) , м-1 Га (Х1) Га (Х2 )
Види- мость очень плохая 1 0.059 3.2510-3 3.4010-3 3.25-10-3 3.40-10-3
Види- мость плохая 5 0.100 5.73-10-4 6.2110-4 2.8610-3 3.11-10-3
Види- мость средняя 10 0.126 2.6410-4 2.9210-4 2.64-10-3 2.92-10-3
Види- мость удов- летво- ритель- ная 30 0.182 7.41-10-5 8.5810-5 2.22-10-3 2.57-10-3
Види- мость хоро- шая 50 0.216 4.010-5 4.7610-5 2.00-10-3 2.38-10-3
Види- мость отлич- ная 1000 0.585 6.4 10-7 1.01-10-6 0.64-10-3 1.01-10-3
Из табл. 4 видно, что ослабление лазерного излучения аэрозолями - относительно небольшое на расстоянии Ят . Оценим максимальный размер зоны ЧС, считая, что аэрозоли ослабляют мощность лазерного излучения в 10 раз.
При этом
е2г““ = 10.
Отсюда
1
Г =—1п10 = 1.15. 2
Тогда
Я. =
1.15 «а (X) .
(5б)
Результаты расчета приведены в табл. 5.
Из табл. 5 видно, максимальные размеры зоны ЧС, подлежащий диагностике, в условиях очень плохой и плохой видимости ограничиваются сотнями метров - одним (двумя) километрами.
При лучшей видимости главную роль уже играет не ослабление аэрозолями, а мощность лазера и удаление до зоны ЧС.
Заметим, что Ят описывает видимость в видимом диапазоне волн. Значения Ят изменяются от 1 до 1000 м. В то же время значения Ятах (Х1) изменяется в пределах: 354-1750 км, а Ятах (Х2) - в пределах: 0.338-
1.619 км. Отсюда следует исключительно важный для служб МЧС вывод.
В условиях визуальной невидимости (“слепоты”) лазеры на С02 и СО позволяют производить дистанционную диагностику загрязнений зон ЧС с размерами от 340 до 1750 м. При очень крупных авариях зоны ЧС могут достигать таких размеров.
Примером подобной аварии, точнее катастрофы, были взрывы военных складов вблизи г. Мелитополя в мае 2004 г.
Метод дифференциального поглощения. Однородное распределение загрязнений в атмосферном газе
Обратимся к уравнению (2). В этом уравнении только Г зависит от длины волны Х (частоты V ) электромагнитного излучения. Для исключения неизвестных (или известных, но с большей погрешностью) величин измерения следует вести на двух длинах волн Х12 (двух частотах V12 ). Одна из длин волн должна соответствовать максимуму в спектре поглощения, а вторая - должна быть вне этого спектра. Представим уравнение (2) в виде:
Р, (Х) = Ве
-2Г(Х)
(6)
РЗпС где В = ——.
К0
Тогда из (6) имеем: ЬР5 (Х) = 1пВ - 2Г(Х).
(7)
При измерении на двух длинах волн для каждой составляющей, загрязняющей атмосферу, получим:
ЬР5 (Х1 ) = 1пВ - 2Г(Х1),
ЬР5 (Х2 ) = 1пВ - 2Г(Х2).
После вычитания (5) из (8) имеем:
Г(Х )-Г(Х2 )=
2 Р, (Х2)'
(8)
(9)
(10)
Учтем, что:
Е
Г(Ч2 )=JN (R)o(R, і12 ) dR.
(11)
Важно, что интегрирование в (11) ведется вдоль трассы.
При однородном распределении атмосферных примесей:
г(х12 ) = no(X1j2 )R0.
Если размер облака выбросов l << R0, то: Г(Х,2 )= No(Xi,2 )l.
(12)
(13)
В соотношениях (12), (13) N - искомая концентрация примесей.
Из (12) следует, что:
ДГ = Г(Х2 )-Г(Х1 )= N (о(Х2 )-о(Х1 ))Я0 или
ДГ = ^оЯ0,
где До = о(Х2)-о(Х1).
С другой стороны, из (10) следует, что:
ДГ=11„Р<Ц.
2 Р, М
Тогда из (14), (15) получаем:
N = (2Д„Ко П-Ц).
Аналогично из (13) и (14) следует, что:
Р, (Х1)
N = (2Aol)ln
Ps (І2)'
(14)
(15)
(16)
(17)
0(v)=0(vL )
vL + (v-v0 )2’
(18)
v - текущее значение частоты, Vl - лоренцева частота.
Тогда:
. do . do
До = —Да = —Av , dA dv
где ДА = A2 - А1, Av = v1 - v2.
Неоднородное распределение загрязнений в атмосферном газе
В этом случае для нахождения N вдоль лазерного луча необходимо решать уравнение (11).
Уравнение (11) - интегральное уравнение Вольтер-ры. Для решения этого уравнения необходимо и правую части в (11) продифференцировать по верхнему пределу Я0. Тогда решение (11) примет вид:
N (R0 ) = о(І1,12)
-1 d^Rp, Іі, І2)
dR
Из (15) и (20) следует, что:
N (R0 Н2*»)-'^^
(20)
(21)
Длянахождения учтем,чтозависимость о(А) или o(v) описывается лоренцевой кривой, т.е.:
,,2
где v0 - центральная частота,
Соотношение (21) дает значение концентрации примесей на расстоянии Я0.
На практике при нахождении решения (21) возникает трудность, связанная с численным дифференцированием, которое приводит к неустойчивости (случайным колебаниям) искомой величины N(Я0).
Бороться с этой трудностью можно применением регуляризации по Тихонову. Для этой цели разработаны специальные регуляризующие алгоритмы [12].
Выводы
Дистанционный лазерный мониторинг позволяет решить при помощи одного комплекта аппаратуры двуединую задачу - оперативного мониторинга как газовых примесей, так и аэрозолей в зоне ЧС.
L
Литература
1. Беккер А.А. Охрана и контроль загрязнения природной среды. - Гидрометеоиздат, 1989.
2. Стокгольмская конференция по окружающей среде 1972 года.
3. Конференция ООН по окружающей среде и развитию в Рио-де-Жанейро 1992 года.
4. Межерис Р. Лазерное дистанционное зондирование. - М.: Мир, 1979.
5. Лазерный контроль атмосферы. под редакцией Э.Д. Хинкли. М.: Мир, 1979. 386 с.
6. В.Е. Зуев, В.В. Зуев. Дистанционное оптическое зондирование атмосферы. С.-П.: Гидрометеоиздат, 1992. 212 с.
7. P.L. Meyer, M.W. Sigrist. Atmospheris pollution using CO2-laser photoacoustic spectroscopy and other technigues. Rov. Sci. Instrum.,
1990, v. 61, № 7, pp. 1779-1807.
8. Б.И. Васильев, У.М. Маннун. ИК лидары дифференциального поглощения для экологического мониторинга окружающей
среды. Квантовая электроника, 2006, т. 36, № 9, с. 801-820.
9. Л.Ф. Черногор. Дистанционное радиозондирование атмосферы и космоса. Х.: ХНУ имени В.Н. Каразина, 2009. 348 с.
10. А. Исимару. Распространение и рассеяние волн в случайно -неоднородных средах. Т. 1. Однократное рассеяние и теория переноса. М.: Мир, 1981. 281 с.
11. Л.С. Ивлев, С.О. Андреев. Оптические свойства аэрозолей. Л.: ЛГУ, 1986. 278 с.
12. А.Н. Тихонов, В.Я. Арсенин. Методы решения некорректных задач. М.: Наука, 1979. 288 с.
3
