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

Вовна А.В., Рак А.И., Хламов М.Г.

Донецкий национальный технический университет, кафедра ЭТ

     Abstract.
    Vovna A.V.,, Khlamov M.G., Rak A.I. The mathematical model of an optical channel of concentration measurer of carbon oxide of infrared gas analisator. The mathematical model of an absorbing spectrum of infrared radiation by carbon oxide is offered. The model can be used in the defining algorithm of the carbon oxide concentration in automobile exhaust gases. The introduction of an additional channel and a developed algorithm provide termination of soot influence in the working optical channel, and also the influence calculation of temperature, pressure and concentration of dioxide on the measurement result of carbon oxide concentration.

     Метод, положенный в основу построения быстродействующего устройства измерения концентрации оксида углерода, основан на способности атмосферных газов избирательно поглощать инфракрасное излучение в спектральной области, и относится к бездисперсионным методам абсорбционной спектроскопии. В соответствии с экспериментальными данными о молекулярном поглощении [1] спектр поглощения оксида углерода имеет сложную структуру. Спектр поглощения оксида углерода, полученный при высокой разрешающей способности анализатора, включает в себя спектральные полосы, каждая из которых содержит электронно-колебательно-вращательные компоненты и их комбинационные составляющие, отличающиеся местом положения в спектральной области и интенсивностью поглощения. Наиболее интенсивная колебательно-вращательная полоса поглощения лежит в области с центром l = 4,66 мкм. Спектр структуры носит характер длинно- и коротковолновых не полностью симметричных крыльев при центральной вращательной составляющей колебательно-вращательные полосы поглощения, представляющие собой последовательности линий поглощения, контур каждой из которых определяется межмолекулярными соударениями и описывается коэффициент поглощения K(l) полосы, который выражается суммой коэффициентов поглощения отдельных линий:

Коэффициент поглощения (1)
     где Интегральная интенсивность i-й линии - интегральная интенсивность i-й линии; di - полуширина линии; l0i - центр линии.

     Для вычисления ряда (2-6) вводились различные допущения о характере изменения ai, di и n0i в пределах колебательно-вращательной полосы. У. Эльзассер рассматривал полосу поглощения как бесконечную последовательность одинаковых по интенсивности (ai=const) и полуширине (di=d) равноотстоящих линий. Для упрощения расчетов была рассмотрена модель полосы поглощения, в которой равноотстоящие линии равной полуширины d имеют распределение интенсивности в каждой из двух ветвей полосы, соответствующее распределению молекул по вращательным уровням энергии. На рис. 1 приведен график функции K(n), вычисленные для полосы поглощения СО с использованием принятой модели полосы.

Спектр поглощения оксида углерода

     Рис. 1 - Спектр поглощения оксида углерода

     В общем виде поглощение излучения в газе описывается законом Бугера-Ламберта:

Закон Бугера-Ламберта (2)

     где Ф0 - поток излучения на границе тестируемого вещества; a - показатель поглощения; l - длина пути излучения в тестируемом веществе.

     Оптическая плотность вещества D устанавливается законом Беера в виде:

Закон Беера (3)

     где С - концентрация активной компоненты вещества, поглощающего излучение; k - удельный показатель поглощения, l - длина кюветы.

     Вычисление концентрации поглощающих частиц возможно путем измерения величины k(n), характеризующей интенсивность линии поглощения, и параметров контура линии поглощения. В реальных условиях прямое исследование контура линии поглощения представляет весьма сложную задачу. Поэтому на практике при определении концентраций атомов и молекул измеряют интегральную интенсивность линий (полос) поглощения. Аналитический сигнал в этом случае определяется разностью интенсивностей зондирующего излучения до и после кюветы с поглощающей газовой средой. Аналитическую связь между изменением интенсивности зондирующего излучения и концентрацией поглощающих частиц находят экспериментально и используют в виде градуировочных графиков. Для окрестности ln0 т.е. l0+Dl/2 процесс поглощения частотно зависим и определяется функцией интегрального коэффициента поглощения. Эта зависимость учитывается следующим образом:

Интегральный коэффициент поглощения Оптическая плотность (4,5)

     Коэффициент передачи оптического канала КОК(l,С) определяется как:

Коэффициент передачи оптического канала (6)

     Как известно, в инфракрасной спектрометрии влияние давления и температуры на оптическая плотность исследуемого вещества учитывается следующим образом: от давления сильно зависит полуширина спектральной линии поглощения и в меньшей от изменения температуры, огибающая спектра полосы тоже прямо пропорциональна изменению давления и обратно пропорциональна квадратному корню изменения температуры. Наибольшее значение интегрального коэффициента поглощения оксида углерода от давления и температуры, рассчитано с использованием аппроксимирующего полинома, который получен на основании статистических данных [2]:

Аппроксимирующий полином зависимости коэффициента поглощения оксида углерода от давления и температуры (7)

      Зависимость коэффициента поглощения оксида углерода определяется с учетом возмущающих факторов по следующей зависимости (см. рис. 2):

Нормированный коэффициент поглощения оксида углерода с учетом возмущающих факторов (8)

Коэффициент поглощения оксида углерода с учетом возмущающих факторов (9)

Нормированный коэффициент поглощения оксида углерода с учетом возмущающих факторов
Рис.2 - Зависимость коэффициента поглощения оксида углерода от давления и температуры
1 - КСО (l, Р=101325, Т=273)
2 - КСО (l, Р=101325, Т=393)
3 - КСО (l, Р=125000, Т=273)
4 - КСО (l, Р=125000, Т=393)

     Поток излучения, создаваемый светодиодом и направляемый на облучение ОК, представляет собой входной сигнал ОК. Поток излучения, прошедший слой тестируемого вещества, образует выходной сигнал ОК. Этот сигнал содержит информацию о концентрации ингредиента в тестируемом веществе.

     Описывается он в следующем виде:

Выходной сигнал оптического канала. (10)

     Обе функции представляются графически, причем график функции SвхОК(l) имеет вид SвыхОК(l, C=0) , а график функции SвыхОК(l, C) для различных значений концентрации ингредиента приведен на рис. 3. Полученные выражения спектральных плотностей потоков излучения используются для построения спектров входного и выходного сигналов ОК. Спектр выходного сигнала строится при максимальном значении концентрации ингредиента.

 Спектр выходного сигнала ОК при различных концентрациях ингредиента
Рис.3 - Спектр выходного сигнала ОК при различных концентрациях ингредиента
1 - SвыхОК(l, C=0,00)
2 - SвыхОК(l, C=0,05)
3 - SвыхОК(l, C=0,50)
4 - SвыхОК(l, C=0,75)
5 - SвыхОК(l, C=1,00)

     Информация о концентрации ингредиента в выходном сигнале представлена в виде изменения формы спектра. Входной величиной всех фотоприемников является поток излучения, а не отдельные спектральные составляющие. Выходной сигнал фотоприемника пропорционален интегральной характеристике сигнала:

Выходной сигнал фотоприемника Выходной сигнал оптического канала (11,12)

Выходной сигнал оптического канала (13)

     Результатом измерений должна быть концентрация, а не величина выходного потока. Для этого следует установить обратную зависимость:

Концентрация - обратная зависимость от выходного сигнала оптического канала (14)

     Одной из задач представленной работы является установление теоретической функции преобразования "концентрация в поток" и ей обратная "поток в концентрацию". Последняя предназначена для градуировки оптического канала и решения измерительных задач. Зависимость концентрации от потока излучения оказывается более наглядной и воспринимаемой, если она устанавливается не для потока излучения, а его относительного изменения:

Зависимость концентрации от относительного значения потока излучения (15)

     Эти зависимости представляются в графической форме. Обратная функция W-1 представляет собой решение уравнения относительно С. В этом выходной поток оптического канала входит в виде интегрального соотношения со сложной зависимостью от длины волны. Это значительно осложняет поиск решения в аналитическом виде, но за это решение может быть принято табличное представление зависимости, где аргумент и функцию легко поменять местами.

     В разрабатываемой математической модели учтено подавляющее большинство факторов, влияющих на прохождение инфракрасного излучения по кюветам газоанализатора.

Оптический канал спектрометрического абсорбционного газоанализатора
Рис.4 - Оптический канал спектрометрического абсорбционногогазоанализатора

      Оптическая неоднородность канала вызвана применением в канале материалов с различными оптическими свойствами. Оптическая неоднородность учитывается коэффициентами преломления. Потери излучения на границах раздела и в средах с различными коэффициентами преломления учитываются коэффициентами отражения и пропускания границами раздела, рассеяния и поглощения излучения в этих средах. В настоящей работе в качестве таких сред выступают оптические линзы, изготавливаемые из специальных материалов: специальных сортов оптического стекла, природных и синтетических кристаллов. Рассматриваемые линзы считаются тонкими (их толщина и диаметр много меньше фокусных расстояний) в связи, с чем потерями на рассеяние и поглощение можно пренебречь. По этой же причине будем считать, что потоки излучения падают на границы раздела сред по нормали к поверхности раздела.

     В диапазоне малых углов падения луча на границу раздела, коэффициент отражения практически остается постоянной величиной. Поэтому малое отклонение угла падения на границу раздела обусловленное кривизной линзы практически не влияет на коэффициент отражения. Сделанное допущение следует считать справедливым. Источник излучения имеет конечные размеры и поток выходящий из линзы оказывается расходящимся, поэтому при определении коэффициента ввода излучения в объектив учитывается угол расхождения луча.

     Эффективность ввода излучения в окно фотоприемного устройства оценивается коэффициентом ввода, при этом полагается, что интенсивность потока по сечению однородна. Выходной сигнал фотоприемника представляет собой сигнал напряжения, поскольку фотоприемное устройство - фоторезистор с интегральным усилителем напряжения, и определяется соотношением:

Выходной сигнал фотоприемника (16)

     где Scинт - интегральная чувствительность фотоприемника к сигналу излучателя; эта характеристика связана со спектральными характеристиками как фотоприемника, так и излучателя и отличается от интегральной чувствительности фотоприемника. Обе характеристики определяются следующим образом:

Sинт - интегральная чувствительность фотоприемника Sсинт - интегральная чувствительность фотоприемника к сигналу (17,18)

     Подынтегральная функция первой формулы содержит спектральную плотность чувствительности фотоприемника в то время, как вторая формула в качестве подинтегральных функций содержит нормированные функции спектральной плотности чувствительности фотоприемника и плотности потока излучения источника.

     На основе приведенной математической модели построена имитационная модель измерителя концентрации оксида углерода, учитывающая влияние комплекса возмущающих факторов, и предложены технические решения, обеспечивающие повышение точности и оперативности получения измерительной информации.

ВЫВОДЫ

     1. Предложенная математическая модель оптического канала измерителя концентрации оксида углерода учитывает комплекс возмущающих факторов, а именно: температуру и давление выхлопных газов, влияние концентрации диоксида углерода и сажи. Учет этих факторов позволяет повысить точность измерения концентрации оксида углерода.

     2. Методы математического и имитационного моделирования, а также современные оптические и микроэлектронные средства позволяют разработать инфракрасные газоанализаторы высокого быстродействия и точности.

ЛИТЕРАТУРА

  1. Зуев В.Е. Распространение видимых и инфракрасных волн в атмосфере. М.: Советское радио, 1970. - 496 с.
  2. Бреслер П.И. Оптические абсорбционные газоанализаторы и их применение. Л.: Энергия, 1980. - 164 с.


     Вовна А.В., Рак А.И., Хламов М.Г. Математическая модель оптического канала абсорбционного измерителя концентрации оксида углерода инфракрасного газоанализатора. - Научные работы Донецкого национального технического университета. Серия: Вычислительная техника и автоматизация: - Донецьк: ДонНТУ,-2003.