УДК 504.064.37
DOI: 10.18698/1812-3368-2019-3-34-45
КОНТРОЛЬ ДОСТОВЕРНОСТИ ПОКАЗАНИЙ СПЕКТРОРАДИОМЕТРА В УСЛОВИЯХ ЗАПЫЛЕННОЙ АТМОСФЕРЫ
Р.Н. Садовников И.В. Кудымова А.С. Самородов
krz525@mail.ru irina.kudymova.83@mail.ru
ФГБУ «33 ЦНИИИ» Минобороны России, Вольск-18, Саратовская обл., Российская Федерация
Aннотация
Рассмотрена возможность контроля достоверности результата работы инфракрасного спектрорадио-метра при регистрации токсичного химиката в запыленной атмосфере. Предложен метод оценки интегральной концентрации аэрозольных частиц на основе определения изменения контраста тестового объекта относительно фона в видимом диапазоне длин волн при прохождении газоаэрозольного облака. Критическое значение контраста, соответствующее концентрации аэрозоля, когда уже невозможно срабатывание спектрорадиометра на наличие в воздухе токсичного химиката, предложено определять на основе численного решения задачи переноса излучения в среде с заданными параметрами. Изложены результаты решения задачи для мешающего аэрозоля, состоящего из капель воды диаметром 5...20 мкм. Показано, что при значении отношения сигнал/шум 100 возможна регистрация токсичного химиката при концентрации аэрозоля до 5 г/м2. Реализация предложенного метода предполагает использование видеокамеры для наблюдения за тестовым объектом и компьютера для обработки получаемых изображений
Ключевые слова
Аварийный выброс, достоверность результатов, газоаэрозольное облако, контраст, спектрорадио-метр, тестовый объект
Поступила 17.04.2018 © Автор(ы), 2019
Введение. В настоящее время появление заражения воздуха и отслеживание перемещения облака токсичного химиката возможно практически в режиме реального времени. Одними из наиболее эффективных технических средств, которые позволяют реализовать такую возможность, являются спектрорадиометры, работающие в инфракрасном диапазоне длин волн [1]. Применение таких приборов для контроля обстановки на
предприятиях химической промышленности позволит в случае аварии обеспечить принятие неотложных мер по защите персонала и населения.
Аварийные выбросы могут сопровождаются взрывами, в результате которых происходит бризантное дробление окружающих объектов [2-4]. Вследствие этого, даже если токсичный химикат в условиях свободной атмосферы является газообразным веществом, в результате аварии с высокой вероятностью сформируется газоаэрозольное облако.
Содержащийся в облаке аэрозоль, не представляя угрозы для людей, является одним из мешающих факторов распознавания наличия в воздухе токсичного химиката. Результаты анализа маскирующей роли аэрозольного компонента облака при распознавании газообразной токсичной примеси с использованием инфракрасного спектрорадиометра приведены в работе [5]. Показано, что отсутствие срабатывания спектрорадиометра, нацеленного на облако выброса, еще не означает, что в воздухе отсутствуют токсичные химикаты, опасные для человека. Для окончательного вывода необходимо решить задачу оценки концентрации в облаке мешающего аэрозоля, а затем определить возможность срабатывания спектрорадиометра.
В настоящее время определение концентрации аэрозоля в воздухе возможно несколькими способами. Наиболее простой и распространенный способ заключается в применении ловушек, использующих методы фильтрации, импакции и аспирации [6-8]. Однако такой подход нецелесообразен вследствие его низкой оперативности. Действительно, значения концентрации могут быть получены только после прохождения облака, съема ловушек и обработки полученных данных.
Оценка параметров аэрозольного компонента облака может быть осуществлена практически в режиме реального времени с помощью ли-дарной системы. Использование многочастотного зондирования позволяет оценить микрофизические характеристики аэрозоля и вычислить его концентрацию [9, 10]. Однако несмотря на привлекательность такого способа, он обладает определенными недостатками. В первую очередь это относится к тому, что в настоящее время серийно не выпускают подобных устройств, а разработка единичного образца потребует значительных затрат. Кроме того, затруднительно обеспечить непрерывную работу ли-дара в течение длительного периода времени, так как это потребует достаточно частого обслуживания, предполагающего наличие специального высококвалифицированного персонала.
Указанные трудности обусловливают необходимость разработки метода для определения параметров аэрозольного компонента облака, об-
ладающего высокой оперативностью и не требующего использования оригинальной дорогостоящей аппаратуры.
Постановка задачи. Пылевой компонент газоаэрозольного облака оказывает маскирующее воздействие при регистрации спектров поглощения газового компонента за счет рассеяния и поглощения излучения в инфракрасной области спектра. Следует отметить, что рассеяние и поглощение излучения будет происходить и в видимой области спектра.
Следовательно, концентрацию аэрозоля можно оценивать по снижению контраста в оптическом диапазоне изображения некоторого объекта до аварии и изображения, полученного в момент, когда облако выброса расположено между фиксируемым объектом и регистрирующей аппаратурой. Далее, зная оценку концентрации аэрозоля, можно сделать вывод о возможности срабатывания спектрорадиометра, работающего в инфракрасном диапазоне.
В качестве регистрирующей аппаратуры можно использовать серийно выпускаемый образец цифровой фото- или видеокамеры. Подключение регистрирующей аппаратуры к компьютеру позволит периодически получать и обрабатывать изображение тестового объекта в целях определения его контраста, например, на фоне открытого неба. Схема расположения аппаратуры и тестового объекта приведена на рисунке.
Расположение регистрирующей аппаратуры и тестового объекта относительно
газоаэрозольного облака выброса: 1 — аварийный объект; 2 — тестовый объект; 3 — газоаэрозольное облако, содержащее токсичный химикат; 4 — видеокамера; 5 — спектрорадиометр
Если затем решить обратную задачу распространения света в рассеивающей среде, то будет получено значение концентрации аэрозоля в об-
лаке, что позволит сделать вывод о возможности срабатывания спектро-радиометра при наличии токсичного вещества в виде газового компонента облака выброса.
Анализ процесса снижения контраста изображения облаком выброса. Если яркость объекта обозначить через Bо, а яркость — фона через Bф, то контраст между ними можно описать формулой [6]:
Во - Вф
L = ■
Во
В более общем случае, когда соотношение яркостей отдельных частей изображения может быть произвольным, выражение для определения контраста имеет вид
^ _ 1В° _ Вф|
" Во + Вф .
Наличие аэрозоля между объектом и наблюдателем приводит к ослаблению световых лучей, идущих от предмета и фона, и к рассеянию аэрозолем света, поступающего со всех других направлений. Яркость объекта Во, перед которым находится облако аэрозоля, можно оценить по соотношению
Во (И) = Во ехр(~оСмИ), (1)
где а — сечение рассеяния излучения одной частицей; Сы — счетная концентрация аэрозольных частиц; И — толщина слоя аэрозоля. В этом случае яркость фона
Вф (И) = Вф ехр (-аСмИ). (2)
Внешний свет, рассеянный аэрозолем, создает свечение, яркость которого О накладывается на яркость объекта и фона. В этих условиях контраст [6]
ь = Во (И)-Вф (И) = Во - Вф | ехр (~оСмИ) (3)
Во (И) + Вф (И) (Во + Вф )ехр (-аСмИ) + О '
В рамках рассматриваемой задачи можно полагать, что светорассеива-ющие частицы равномерно заполняют все пространство и их число в единице объема достаточно велико. В этом случае интенсивности являются (приближенно) непрерывными функциями толщины светорассеивающего слоя, и можно применять методы дифференциального исчисления.
Отражение, рассеяние вперед и поглощение света слоем единичной толщины можно описать с использованием коэффициентов
г = Сыст^ (т,й); I = Сы) (т, й); ц = СN&(з) (т,й),
где о^"* (т, й)— интегральное сечение рассеяния частицей электромагнитного излучения в диапазоне углов тс / 2...тс; а^ (т, й) — интегральное сечение рассеяния частицей электромагнитного излучения в диапазоне углов 0...тс/2; а(я) — интегральное сечение поглощения частицей электромагнитного излучения; й — диаметр частицы; т — комплексный показатель преломления вещества аэрозоля.
Счетную концентрацию легко определить, зная массовую концентрацию монодисперсного аэрозоля:
Сы =-3 >
жрй
где С — массовая концентрация аэрозоля в воздухе; р — плотность вещества аэрозоля.
Очевидно, что суммарное ослабление прямого излучения будет определяться суммой всех коэффициентов а = t + г + д. При этом коэффициент а представляет собой Сыо из соотношений (1)-(3).
Электромагнитное излучение в рассеивающем слое можно разложить на два потока: первый распространяется по направлению падающего излучения, второй — навстречу ему. Первый поток состоит из нерассеянного света и света, испытавшего одно, два или более рассеяний; второй создается только рассеянным светом. Представим рассеянное излучение в виде суммы потоков однократно и более раз рассеянного света. Обозначим световые потоки к-й кратности, рассеянные по направлению и
навстречу падающему излучению, через 1(1) и I(2) соответственно, где к = 1, 2, ... [11].
Рассмотрим прохождение света через бесконечно тонкий слой йг, находящийся на глубине г. Если интенсивность света на глубине г равна 4», то при прохождении через слой йг часть его рассеется, часть поглотится, а некоторая часть пройдет без рассеяния. Предположим, что рассеяние вперед и назад (отражение), а также поглощение слоем йг можно полагать пропорциональным толщине йг. Тогда доля света, отраженного слоем йг, равна ^гйг, рассеянного вперед — поглощенного —
а прошедшего без рассеяния — - (г + t + |д)].
Изменение интенсивности прямого света в слое йг, обусловленное рассеянием и поглощением, имеет вид
^ .-«IО). йг
Интенсивность однократно рассеянного вперед света уменьшается в результате рассеяния и поглощения и увеличивается за счет рассеяния прямого света в слое йг:
йт (1)
Ё^ = _а7(1) +1 (1).
йг
Изменение интенсивности рассеянного к-й раз света происходит по тем же причинам, что и однократного, но световой поток дополнительно пополняется в результате отражения рассеянного света (к—1)-й кратности, распространяющегося назад:
Ё1к г« , г(2) , .Л1)
-г- = -о1к +Пк-1 +Ик-1.
Получим систему уравнений для света, рассеянного навстречу падающему излучению:
, (2)
_ ^ =_а71(3) + Г1 (1);
йг
^ -al(2) + rl(1) + tl(2)
йг
Граничные условия:
- падающее излучение на передней границе рассеивающего слоя имеет начальную интенсивность 1о: I(1) (0) = То;
- рассеянное вперед произвольное число раз излучение на передней границе рассеивающего слоя полностью отсутствует: 1(1) (0) = 0;
- рассеянное назад произвольное число раз излучение на задней границе рассеивающего слоя полностью отсутствует: I(2) (Н) = 0, где Н — толщина рассеивающего слоя.
Последовательное решение приведенных уравнений необходимо вести вплоть до учета ^-кратного рассеяния света, пока не будет выполнено условие
I(1) (Н(Н(0)« 10. к=1 к=1
Будем полагать, что в качестве фонового объекта используется квадрат черного цвета. Поскольку яркость такого объекта равна нулю, выражение (3) принимает вид
L = Дф ехр (-оШ) . (4)
Bф ехр (-апИ) + О
Учитывая, что цветопередача видеоаппаратуры близка к идеальной, выражение (4) можно переписать следующим образом:
I(1) (И)
L = ■
N
I(1) (И)+2 к (И) к=1
Вариант решения задачи контроля показаний спектрорадиометра.
Результаты расчетов значений контраста тестового объекта при различных концентрациях аэрозоля приведены в таблице. Ширина факела выброса выбрана равной 50 м, плотность вещества аэрозоля — 1 г/см3. При проведении расчетов предполагалось, что коэффициент преломления составлял 1,33. Факторы эффективности рассеяния и поглощения электромагнитной энергии выбраны равными [12]:
Ку = 4а(/У(ттй2 ) = 2,1; Кп = 4а(п5)/(ттй2 ) = 0,1; Ка = 4а(я)/(ттй2 ) = 0.
Результаты расчетов значений контраста тестового объекта при различных концентрациях аэрозоля
Концентрация аэрозоля Са, г/м2 Диаметр аэрозольных частиц, мкм
5 10 20
0,25 0,854 0,924 0,961
0,5 0,730 0,854 0,924
1 0,532 0,730 0,854
2,5 0,206 0,455 0,674
5 0,042 0,206 0,455
10 0,002 0,042 0,206
Данные, приведенные в таблице, можно аппроксимировать аналитической зависимостью
L = exp
-3'15 d J •
(5)
Преобразование выражения (5) относительно Са дает результат
Са =-0,317 ln(L) d. (6)
Оценочная формула для определения минимальной интегральной концентрации газообразного токсичного химиката, позволяющей распознать наличие опасности при заданном значении концентрации аэрозоля в облаке, предложена в работе [5]:
_ __ 2-107 1 + 1 - 4ехр [ 3 Са К/№р)]ц{
— /„ ч |п , (7)
аг (Яр )Ыл 2
где Шц — молярная массы газообразного соединения; аг (Л,р ) — сечение
поглощения на резонансной частоте Яр для одной молекулы; Ыл — число Авогадро; й — медианное значение диаметра аэрозольных частиц; ^ — отношение сигнал/шум для спектрорадиометра; К — фактор эффективности рассеяния излучения аэрозольной частицей.
Из формулы (7) следует, что возрастание концентрации аэрозоля Са требует наличия более высокой концентрации газообразной примеси Сшш. При этом максимальное значение концентрации аэрозоля, при котором может быть обнаружено наличие газовой примеси, определяется уравнением, получаемым из (7):
1 - 4ехр [ 3Са К/ (2йр)]^_1 = 0.
Решение имеет вид
с(кр)=^ - е). (8)
Формула (8) получена при условии, что отношение интенсивности излучения на резонансной частоте и интенсивность фонового излучения отличаются в 2 раза. Если интенсивности будут отличаться в большее число раз, то критическая концентрация аэрозоля будет больше. Следовательно, с помощью формулы (8) можно оценить снизу максимальную концентрацию, при которой существует вероятность, что спектрорадио-метр не обнаружит в воздухе токсичный химикат.
Подставляя выражение (8) в (6) и выполняя тождественные преобразования, получаем формулу для определения критического значения контраста изображения черного квадрата:
Ькр = exp
2Лр. in'
K
,3
(9)
Поскольку плотность вещества аэрозоля выбрана равной 1 г/см3, а К = Ку + Кп = 2,2, формула (9) для выбранных параметров принимает окончательный вид
, л-0,995
Lкр=[ 4J . (10)
Согласно (10), при ц = 10 получаем L^ = 0,402, а при ц = 100 Lкр = = 0,041. Сопоставляя вычисленные критические значения c данными, приведенными в таблице, устанавливаем, что при ^ = 10 срабатывание спектрорадиометра возможно при интегральной концентрации мешающего аэрозоля до 1 г/м2, а при ^ = 100 возможно срабатывание и при концентрации до 5 г/м2.
Заключение. Предложенный метод позволяет осуществлять оценку интегральной плотности аэрозоля в аварийном облаке выброса в режиме, близком к реальному времени, что делает возможным повышение достоверности показаний спектрорадиометра о наличии газообразных токсичных примесей при работе в наиболее сложных условиях. Определение невозможности срабатывания спектрорадиометра вследствие запыленности атмосферы позволит внести коррективы в принимаемые меры по защите персонала и населения, а также снизить вероятность поражения людей.
Для реализации метода требуется наличие цифровой видеокамеры, подключенной к компьютеру, а также специализированное программное обеспечение для обработки получаемых изображений тестового объекта. Это обусловливает его практическую значимость, так как делает возможным не только создание более совершенных систем контроля химической обстановки, но и достаточно простую модернизацию уже существующих систем в целях повышения их эффективности.
ЛИТЕРАТУРА
[1] Морозов А.Н., Светличный С.И. Основы фурье-спектрорадиометрии. М., Наука, 2014.
[2] Лебская Т.А. Прогнозирование чрезвычайных ситуаций на потенциально опасных производственных объектах. Инженерный вестник, 2016, № 2.
URL: http://engbul.bmstu.ru/doc/833854.html
[3] Соколов Э.М., Панарин В.М., Горюнкова А.А. и др. Совершенствование теории развития чрезвычайных ситуаций техногенного характера. Известия ТулГУ. Технические науки, 2010, № 4-2, с. 241-251.
[4] Матвеев Ю.Н., Карельская К.А., Стукалова Н.А. и др. Структуры математических моделей аварийных и чрезвычайных ситуаций на химически опасных объектах. Науковедение, 2016, т. 8, № 5.
URL: http://naukovedenie.ru/PDF/89TVN516.pdf
[5] Садовников Р.Н., Бойко А.Ю., Шлыгин П.Е. Оценка влияния особенностей состава газоаэрозольного облака на спектр фонового излучения открытой атмосферы. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2013, № 4, с. 69-78.
[6] Грин Х., Лейн В. Аэрозоли — пыли, дымы и туманы. Л., Химия, 1972.
[7] Толчинский А.Д., Сигаев В.И., Варфоломеев А.Н. Современные тенденции развития технических средств пробоотбора и анализа биоаэрозолей. Экологические системы и приборы, 2008, № 11, с. 30-48.
[8] Назаров Б.И., Абдуллаев С.Ф., Маслов В.А. и др. Исследование функции распределения аэрозольных частиц и вариация сажевого аэрозоля в атмосфере. Доклады академии наук Республики Таджикистан, 2010, т. 53, № 5, с. 358-363.
[9] Зуев В.Е., Зуев В.В. Современные проблемы атмосферной оптики. Т. 8. Дистанционное оптическое зондирование атмосферы. СПб., Гидрометеоиздат, 1992.
[10] Боровлев А.Э., Кунгурцев С.А. Лидарный аппаратно-программный комплекс как элемент геоинформационной системы г. Белгорода. Экологические приборы и системы, 2008, № 11, с. 56-59.
[11] Чекалинская Ю.И. О кратности рассеяния света. Известия АН СССР, 1957, т. 21, № 11, с. 1494-1499.
[12] Борн М., Вольф Э. Основы оптики. М., Наука, 1973.
Садовников Роман Николаевич — д-р техн. наук, профессор, главный научный сотрудник ФГБУ «33 ЦНИИИ» Минобороны России (Российская Федерация, 412918, Саратовская обл., Вольск-18, ул. Краснознаменная, д. 1).
Кудымова Ирина Владимировна — научный сотрудник ФГБУ «33 ЦНИИИ» Минобороны России (Российская Федерация, 412918, Саратовская обл., Вольск-18, ул. Краснознаменная, д. 1).
Самородов Александр Сергеевич — канд. техн. наук, старший научный сотрудник ФГБУ «33 ЦНИИИ» Минобороны России (Российская Федерация, 412918, Саратовская обл., Вольск-18, ул. Краснознаменная, д. 1).
Просьба ссылаться на эту статью следующим образом:
Садовников Р.Н., Кудымова И.В., Самородов А.С. Контроль достоверности показаний спектрорадиометра в условиях запыленной атмосферы. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2019, № 3, с. 34-45. DOI: 10.18698/1812-3368-2019-3-34-45
MONITORING RELIABILITY OF SPECTRORADIOMETER READINGS IN DUSTY ATMOSPHERIC CONDITIONS
R.N. Sadovnikov krz525@mail.ru
I.V. Kudymova irina.kudymova.83@mail.ru
A.S. Samorodov
33rd Central Scientific Research Institute of the Ministry of Defense of the Russian Federation, Volsk-18, Saratov Region, Russian Federation
Abstract
The paper considers the possibility of monitoring reliability of infrared spectroradiometer operation when detecting the presence of a toxic chemical in a dusty atmosphere. We present a method for estimating cumulative concentration of aerosol particles based on evaluating variations in contrast of a test object against the background in the visible wavelength range during gas and aerosol cloud propagation. In order to determine the critical contrast value that corresponds to the aerosol concentration designating further inability of the spectroradiometer to detect the toxic chemical in the air, we propose to numerically solve the problem of radiation transfer in a medium with known parameters. We report the results of solving this problem for an interfering aerosol consisting of water droplets with a diameter of 5-20 |^m. We show that for a signal-to-noise ratio of 100 it is possible to detect a toxic chemical at aerosol concentrations of up to 5 g/m2. Implementing the method proposed involves using a video camera to monitor the test object and a computer to process the images recorded
Keywords
Accidental pollution, result reliability, gas and aerosol cloud, contrast, spectroradiometer, test object
Received 17.04.2018 © Author(s), 2019
REFERENCES
[1] Morozov A.N., Svetlichnyy S.I. Osnovy fur'ye-spektroradiometrii [Fundamentals of Fourier-spectroscopy]. Moscow, Nauka Publ., 2014.
[2] Lebskaya T.A. Forecasting emergency at potentially dangerous manufacturing sites. Inzhenernyy vestnik [Engineering Bulletin], 2016, no. 2 (in Russ.).
Available at: http://engbul.bmstu.ru/doc/833854.html
[3] Sokolov E.M., Panarin V.M., Goryunkova A.A., et al. Improvement of the theory of development emergency. Izvestiya TulGU. Tekhnicheskie nauki [News of the Tula State University. Technical Science], 2010, no. 4-2, pp. 241-251 (in Russ.).
[4] Matveev Yu.N., Karel'skaya K.A., Stukalova N.A., et al. Structures of mathematical models of contingency and emergency situations on chemically dangerous objects. Nau-kovedenie, 2016, vol. 8, no. 5 (in Russ.).
Available at: http://naukovedenie.ru/PDF/89TVN516.pdf
[5] Sadovnikov R.N., Boyko A.Yu., Shlygin P.E. Estimation of influence of features of the gas-aerosol cloud composition on the background radiation spectrum of the exposed atmosphere. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2013, no. 4, pp. 69-78 (in Russ.).
[6] Green H.L., Lane W.R. Particulate clouds: dusts, smoke, and mists. Their physics and physical chemistry and industrial and environmental aspects. Princeton, N.J., Van Nostrand, 1964.
[7] Tolchinskiy A.D., Sigaev V.I., Varfolomeev A.N. Modern trends in development of technical means for bioaerosol sampling and analysis. Ekologicheskie sistemy i pribory [Ecological Systems and Devices], 2008, no. 11, pp. 30-48 (in Russ.).
[8] Nazarov B.I., Abdullaev S.F., Maslov V.A., et al. Studies of the distribution of an aerosol particles and variation soot aerosols in the atmosphere. Doklady akademii nauk Respubliki Tadzhikistan, 2010, vol. 53, no. 5, pp. 358-363 (in Russ.).
[9] Zuev V.E., Zuev V.V. Sovremennye problemy atmosfernoy optiki. T. 8. Distantsion-noe opticheskoe zondirovanie atmosfery [Today problems of atmosphere optics. Vol. 8. Optical remote sensing of the atmosphere]. St. Petersburg, Gidrometeoizdat Publ., 1992.
[10] Borovlev A.E., Kungurtsev S.A. The hardware-software lidar complex as the element of geoinformation analytical system of city of Belgorod. Ekologicheskie sistemy i pribory [Ecological Systems and Devices], 2008, no. 11, pp. 56-59 (in Russ.).
[11] Chekalinskaya Yu.I. On degree of light scatter. Izvestiya AN SSSR, 1957, vol. 21, no. 11, pp. 1494-1499 (in Russ.).
[12] Born M., Wolf E. Principles of optics: electromagnetic theory of propagation, interference and diffraction of light. Oxford, Pergamon Press, 1965.
Sadovnikov R.N. — Dr. Sc. (Eng.), Professor, Chief Research Scientist, 33rd Central Scientific Research Institute of the Ministry of Defense of the Russian Federation (Krasnoznamennaya ul., Volsk-18, Saratov Region, 412918 Russian Federation).
Kudymova I.V. — Researcher Scientist, 33rd Central Scientific Research Institute of the Ministry of Defense of the Russian Federation (Krasnoznamennaya ul., Volsk-18, Saratov Region, 412918 Russian Federation).
Samorodov A.S. — Cand. Sc. (Eng.), Senior Researcher, 33rd Central Scientific Research Institute of the Ministry of Defense of the Russian Federation (Krasnoznamennaya ul., Volsk-18, Saratov Region, 412918 Russian Federation).
Please cite this article in English as:
Sadovnikov R.N., Kudymova I.V., Samorodov A.S. Monitoring reliability of spectro-radiometer readings in dusty atmospheric conditions. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2019, no. 3, pp. 34-45 (in Russ.). DOI: 10.18698/1812-3368-2019-3-34-45