Примерами термодинамически неоднородных газовых сред являются: атмосфера, как в своей вертикальной структуре, так и вдоль горизонтальных трасс, выбросы промышленных предприятий, факелы реактивных двигателей, области вблизи нагретых объектов и пожаров, и т. п. Интерес к исследованию их газового состава неслучаен. Изменения газового состава атмосферы Земли, обусловленные антропогенными факторами, приводят к различным экологическим последствиям. Рост содержания парниковых газов (С02, метан и др.) изменяет радиационные свойства атмосферы и, как следствие, приводит к изменению климата планеты [1]. Увеличение содержания фреонов и ряда других малых газовых составляющих атмосферы вызывает снижение концентрации стратосферного озона, что приводит к увеличению уровня жесткой ультрафиолетовой радиации, достигающей поверхности Земли. Наряду с этим замечен рост концентрации тропосферного озона, являющегося сильным окислителем. Исследование характеристик газового состава атмосферы является одной из центральных задач современной физики и химии атмосферы. Анализ термодинамических характеристик и химического состава газообразных продуктов горения углеводородов также важен для мониторинга и разработки средств управления процессом. Из всех способов определения концентраций газов наиболее привлекательны оптические дистанционные методы.
Анализ литературных данных по оптическим методам дистанционных исследований состава и термодинамических характеристик газовых сред позволяет сделать следующие выводы:
1. В настоящее время разработаны оптические методы исследования газового состава сред, включающие в себя разнообразные по используемому спектральному диапазону, источнику поля излучения, задаваемой модели среды, методам разделения сигнала и обращения данных подходы к анализу концентрационных полей неоднородных газовых сред. Это позволяет говорить о возможности комплексного использования различных методов для взаимного контроля полученных данных и извлечения дополнительной информации об исследуемом объекте.
2. Задачи восстановления газового состава обладают некорректностью в классическом смысле и требуют использования априорной информации и регуляризирующего алгоритма, критерием возможности использования которых на сегодняшний день выступает вычислительный эксперимент, проводимый по замкнутой схеме.
3. Все существующие методы восстановления концентрационных полей опираются на термодинамическую модель газовой среды. В связи с этим большое внимание уделяется созданию методов дистанционного измерения термодинамических характеристик среды.
4. Все существующие алгоритмы обращения данных рассматриваемых методов используют численные методы восстановления, что налагает требования оценки устойчивости и дополнительно возникающих погрешностей дискретизации.
5. Современные тенденции развития рассматриваемых методов связаны с созданием аппаратуры высокого спектрального разрешения, прецизионными измерениями оптических характеристик взаимодействия молекул с полем излучения, совершенствованием физических основ дистанционных измерений, и, в частности, разработкой методов расчета характеристик взаимодействия излучения с молекулярной газовой средой, адекватно описывающих экспериментальные данные, и разработкой новых алгоритмов обращения данных, направленных на все большую детализацию восстанавливаемой пространственной структуры концентрационного поля и минимизацию влияния закладываемой модели среды.
Все вышесказанное свидетельствует об актуальности разработки новых алгоритмов обращения данных оптических дистанционных измерений характеристик газовых сред.
Основная цель исследования разработка новых подходов к восстановлению термодинамических параметров и пространственного распределения концентрации неоднородной v газовой среды по измеренным спектральным зависимостям характеристик ее взаимодействия с оптическим излучением.
Для достижения указанной цели необходимо решить следующие задачи:
1. Создать новый математический аппарат обращения данных измерений спектральной зависимости прозрачности неоднородной газовой среды, базирующийся на особенностях пространственного распределения концентраций отдельных газовых компонент.
2. Разработать логическую, функциональную схемы и реализовать их в виде алгоритмов и программ для ЭВМ для следующих задач: модели пирометрических измеренийметодики расчетов коэффициентов поглощения и излучательной способности газовых сред, — предложенного метода обработки и интерпретации экспериментальных данных.
3. Проанализировать влияние на результаты моделирования следующих параметров задачи: аппаратной функции измерительного устройстваотклонений в задании исходной информациивариаций задаваемых моделей средразличных вариантов методов численного решения.
4. Предложить методы определения температуры нагретой газовой среды по данным пирометрических измерений для различных диапазонов оптической толщи, и, следовательно, концентраций молекул.
Методы исследования.
В качестве основных методов исследования выбраны метод математического моделирования и вычислительный эксперимент, проводимый по замкнутой схеме, опирающиеся на: математический аппарат статистических моментовпрямой метод роасчета прозрачности неоднородной газовой среды по стратифицированной моделиаппроксимационные методыметод наименьших квадратов в двух вариантах (нормальные уравнения и ортогонализация Грама — Шмидта) — метод Гаусса для решения систем линейных уравнений и метод дихотомии.
Защищаемые положения:
1. Восстановление начальных моментов пространственного распределения концентрации компоненты термодинамически неоднородной газовой среды при обработке данных прозрачности в спектральном интервале, содержащем центры линий поглощения только этого газа, можно осуществить путем аппроксимации полиномом зависимости коэффициента поглощения этой компоненты от пространственной переменной. Подстановка этой аппроксимации в интегральное уравнение для оптической толщи преобразует это уравнение в систему линейных алгебраических уравнений относительно искомых моментов.
2. В температурно-неоднородных газовых средах для молекулярной малоатомной компоненты в пределах изолированной полосы поглощения возможно ослабление температурной зависимости эффективного коэффициента поглощения при регистрации прибором с аппаратной функцией с центром в точке локального минимума температурной зависимости поглощения за счет частичной компенсации различных по знаку изменений интенсивностей линий поглощения вблизи центра и на краю полосы для оптически тонкого слоя между линиями поглощения и оптически толстого слоя в их центрах. В частности, в неоднородных средах с изменением температуры в пределах 400 -1400 К, эффективный коэффициент поглощения угарного газа, полученный прибором с дисперсионной аппаратной функцией с центром 2189 см" 1 и полушириной 12 см" 1 не зависит от температуры (с точностью до 5%) и является монотонной функцией от общего содержания угарного газа в пределах изменения 0.05 — 20 атм-см (STP), что позволяет при наличии данных об ослаблении излучения другими факторами (т0) по результату измерения прозрачности среды прибором с вышеуказанной аппаратной функцией оценить общее содержание СО с относительной ошибкой Асо, связанной с погрешностью измерений прозрачности, А соотношением, А =Т° + А.
3. Определение кинетической температуры в газовой среды пирометрическими методами может осуществляться по монотонно зависящей от 9 линейной комбинации радиационных температур, измеренных вблизи центра фундаментальной полосы поглощения СО и на ее крыле, коэффициенты которой зависят от выбора аппаратной функции измерительного устройства и остаются постоянными при изменении концентрации излучающих молекул в широких пределах. В частности, для угарного газа при температурах 300 — 1400 К коэффициенты остаются постоянными при изменении содержания СО в пределах трех порядков. Достоверность:
• расчетов коэффициентов поглощения и прозрачности атмосферы подтверждается сравнением параметров спектральных линий и исследуемых величин с экспериментальными данными и применением прямого метода расчета;
• моделирования пирометрических измерений подтверждается совпадением результатов расчетов в случаях оптически тонкого слоя и больших оптических толщ с теоретическими асимптотическими зависимостями;
• функционирования используемых программных комплексов подтверждается многократным тестированием их в других задачах.
Научная новизна:
• Предложен метод поиска особых точек пространственного распределения концентрации компоненты газовой среды, базирующийся на математическом аппарате статистических моментов, что дает возможность дистанционного определения и интерпретации характеристик пространственных профилей концентраций газов термодинамически неоднородной газовой среды, являющийся альтернативой существующих методов.
• Показана возможность определения общего содержания компонентов температурно-неоднородной газовой среды по измерениям ее прозрачности, основанная на компенсации температурной зависимости коэффициента поглощения в максимуме и на крыле полосы поглощения выбором аппаратной функции измерительного устройства.
• Разработан способ дистанционного измерения средней температуры нагретой газовой среды методами пирометрии. Предлагаемая схема опирается на выявление спектральных интервалов с минимальным и максимальным значениями температур, и показано, что линейная комбинация соответствующих радиационных температур однозначно связана с кинетической температурой и сохраняет свое поведение в широком диапазоне общего содержания излучающих молекул.
• Создан не имеющий аналогов программный комплекс по моделированию пирометрических измерений излучения газовой среды и проведены расчеты излучательной способности нагретой (300−3000 К) газовой среды.
Научная ценность:
• Применение аппроксимации полиномом ядра интегрального уравнения Фредгольма 1 рода, являющегося математической основой формулировки многих обратных задач оптики газовых сред, позволяет записать его в виде линейной комбинации начальных моментов целевой функции, что, при наличии набора значений правой части дает возможность сфор-мировать систему линейных алгебраических уравнений относительно указанных величин.
• Выбор аппаратной функции регистрирующего устройства позволяет минимизировать температурную зависимость коэффициента поглощения посредством компенсации изменений коэффициентов поглощения вблизи центра и на крыле полосы.
• Предложен дистанционный метод определения кинетической температуры газовой среды, базирующийся на интерпретации стандартных пирометрических измерений;
• Разработанный программный комплекс RADTEMP позволяет предсказать спектральную излучательную способность газовой смеси СО, С02, Н20 в диапазоне температур 300−3000 К с точностью в пределах погрешностей экспериментов.
Практическая значимость:
— Применение алгоритма (положение 1) к спектру прозрачности атмосферы по вертикальной трассе, регистрируемого с относительной ошибкой не более 3% прибором разрешающим не менее 20 спектральных точек в пределах характерной ширины линии поглощения и возможностью сканирования по спектральному отрезку, содержащему центры линий одного из основных атмосферных газов, позволяет восстановить с отклонениями не выше 10% следующие характеристики высотного распределения концентрации этого газа:
• общее содержание в толще атмосферы;
• высоту, выше и ниже которой содержания исследуемой компоненты равны;
• ширину высотного профиля.
— Измерение прозрачности температурно-неоднородных газовых сред прибором с конкретной аппаратной функцией (положение 2) может быть использовано для оперативного контроля содержания угарного газа.
— Использование разности радиационных температур (положение 3), измеренных вблизи центра полосы поглощения и на ее крыле позволяет в ряде практически важных случаев оценить кинетическую температуру газовой среды без измерения концентрации излучающих молекул.
Содержание диссертации отражает личный вклад автора в опубликованных работах. Они выполнены в соавторстве с научным руководителем Войцеховской О. К. В работе использованы программные комплексы ОРТАТМ и HOTGAS, разработанные совместно с Шатровым Ф. Г. и Шелудяковым Т. Ю. соответственно.
Апробация работы. Материалы, вошедшие в диссертацию, обсуждались на следующих конференциях: 1) Международная конференция «Импульсные лазеры на переходах атомов и молекул», Томск, 22 — 26 сентября 1997; 2) международная конференция по сопряженным задачам механики и экологии, Томск, 6−10 июля 1998 г.- 3Международная конференция «Импульсные лазеры на переходах атомов и молекул», Томск, сентября 1999 г.- 4) 1-я школа-семинар молодых ученых «Современные проблемы физики и технологии», Томск, 1−3 февраля 2000 г.- 5) научная конференция «Фундаментальные проблемы охраны окружающей среды и экологии природно-территориальных комплексов Западной Сибири». — Горно-Алтайск, 27−30 июня 2000 г.- 6) международная конференция «Сопряженные задачи механики и экологии», Томск, 4−9 июля 2000 г.- 7) VII международный симпозиум «Оптика атмосферы и океана», Томск, 16 — 19 июля 2000 г.- 8) 2-я школа-семинар молодых ученых «Современные проблемы физики и технологии», Томск, 5−7 февраля 2001 г.
Содержание диссертационной работы. Диссертация изложена на 121 странице, содержит 36 рисунков и 6 таблиц, и состоит из введения, четырех глав, заключения, и списка используемой литературы и работ соискателя из 87 наименований.
ВЫВОДЫ ГЛАВЫ 3.
1. Теоретически обоснована возможность минимизации температурной зависимости поглощения малоатомными молекулярными компонентами специальным выбором аппаратной функции с центром в окрестности нелокальной точки пересечения спектральных зависимостей коэффициентов поглощения при различных температурах в случае оптически тонкого слоя, реализующегося вне линий поглощения и оптически толстого слоя в центрах линий поглощения.
2. Произведен выбор аппаратной функции, минимизирующей температурную зависимость эффективного коэффициента поглощения СО в динамическом диапазоне 400 — 1400 К. Выбрана дисперсионная аппаратная функция с центром 2189 см" 1 и полушириной 12 см" 1.
3. Предложена номограмма определения общего содержания угарного газа в температурно-неоднородных средах по измерениям прозрачности с указанной аппаратной функцией.
4. Выполнен вычислительный эксперимент по восстановлению общего содержания СО в температурно-неоднородных средах, показавший, что предложенная методика может использоваться для контроля концентрации угарного газа вдоль трассы с общим содержанием СО 0.05 — 20 атм-см в указанном температурном диапазоне, с отклонениями не более 10%.
5. Предложен способ оценки средней температуры, при которой находится поглощающий газ вдоль оптической трассы.
ГЛАВА 4. ОПРЕДЕЛЕНИЕ ТЕМПЕРАТУРЫ ГАЗОВОЙ СРЕДЫ ПО ДАННЫМ ПИРОМЕТРИЧЕСКИХ ИЗМЕРЕНИЙ.
В главе предложена оригинальная методика измерения кинетической температуры нагретого газового объема по данным пирометрических измерений его собственного теплового излучения без измерения концентрации излучающего газа, основанная на сравнении радиационных температур, измеренных на крыле и вблизи центра колебательно-вращательной полосы поглощения. Следует отметить, что специфика газовой среды, в общем случае, не позволяет использовать результаты пирометрических измерений в качестве оценки кинетической температуры, поскольку коэффициенты поглощения невелики и могут меняться при незначительном изменении частоты на несколько порядков. Для газовых сред невозможно ввести коэффициент серости. В инженерных методиках решение обратных задач, связанных с изучением характеристик нагретых газов базируется на феноменологических подходах, когда экспериментальные данные по характеристикам излучения различных объектов в атмосфере аппроксимируются простыми полиномиальными (или иными, например, функцией Планка) соотношениями, коэффициенты которых определяются методами минимизации [81]. Такой подход не позволяет исследовать физическую картину процессов, происходящих в излучающей газовой среде и не дает оснований для решения обратной задачи по контролю за составом и термодинамикой среды [82,83]. В данном подходе вычисление излучательной способности газов осуществляется по точной методологии прямых расчетов (line by line).
4.1. РАСЧЕТ РАВНОВЕСНОГО ТЕПЛОВОГО ИЗЛУЧЕНИЯ ГАЗОВОЙ СРЕДЫ.
Расчет излучательной способности газовой среды может быть осуществлен по формуле, являющейся прямым следствием законов Кирхгофа и БугераЛамберта: e (v)= $B (v, e{l))a (v, l) expa (v, l) dl [dl.
4.1.1) где v — спектральная координата, / - текущая точка оптической трассы луча, В (у, 6) — функция Планка, 0 — кинетическая температура, a (v, l) — коэффициент поглощения, зависимость которого от I обусловлена изменением вдоль трассы как концентраций излучающих газов, так и термодинамических параметров.
Заметим, что d dl exp< - ja (v, T) dT dl.
4.1.2).
Это дает возможность взять интеграл в (2.2.1) по частям:
LrdB (v, 0(l)) Г ' e (v) = B (v, 0(O))-B (v, 9(L))T+ldl expa (vj)dT (4.1.3) о1*" L о где Г — пропускание среды. В случае однородной газовой субстанции dB (v, 0(l)) dl 0, и формула (4.1.1) переходит в закон Кирхгофа.
Для численного моделирования процессов излучения целесообразно использовать стратифицированную модель среды в виде N однородных слоев. Для этого случая dB (v, 6(l)) n-1 dl I,{B (v, eM)-B (vAM-h)>
4.1.4) к=1 где 6к — температура ктого слоя, 1к — его граница, — функция Дирака, и формула (4.1.1) приобретает вид: n ety) = 2 B (v, е{)(1 — Tt (у) (4.1.5) i=l i.
Имеется также возможность сделать более тонкую аппроксимацию температурно неоднородных сред путем представления их в виде слоев с линейно меняющейся температурой:.
0(1) = вк 9к~1 l = P + Ql, le [lkx, lk ].
U-i.
4.1.6) к-1.
В этом случае можно аппроксимировать Q IQ.
А/ (4.1.7).
Тогда выражение для спектральной излучательной способности принимает вид:.
Использование приемника излучения предполагает осреднение по спектральной координате с некоторой аппаратной функцией A (v, v'). При пирометрических исследованиях излучающих объектов измеряемыми величинами являются эквивалентные температуры, получаемые сравнением одной из характеристик излучения с аналогичной характеристикой излучения АЧТ, например энергия излучения. Соответствующая эквивалентная температура вк называется радиационной и может быть найдена из решения следующего нелинейного уравнения: где пределы интегрирования определяются областью ненулевых значений аппаратной функции..
В качестве параметров среды в слоях задаются температура, концентрации газов и длина оптической трассы. Ограничением на интервал температур является требование, чтобы используемая база данных по параметрам спектральных линий была сконструирована в том же интервале температур. Предельной концентрацией газа является одна атмосфера. Длина пути любая. e (v) = B (v, 0q) — B (v, 6>o + cL) expj — ja (v, l) dl + 0.
4.1.8).
4.1.1.Моделирование пирометрических измерений.
4.1.9).
1600−1 я 1400Н.
CL I.
О.1200Н Ф с 1000 о: го i 800 о s =г g 600 сС го Q..
400 гггр 1.
Концентрация СО 0.05 атм.
I ill.
10 100 1000 Длина излучающего слоя, см.
Рис. 4.2.1. Зависимость регистрируемой радиационной температуры от общего содержания СО. Штриховая линия — кинетическая температура слоя..
4.2. МОДЕЛИРОВАНИЕ СПЕКТРАЛЬНОГО ИЗЛУЧЕНИЯ НЕОДНОРОДНОЙ НАГРЕТОЙ.
ГАЗОВОЙ СРЕДЫ.
Моделирование спектра теплового излучения нагретой газовой среды проводилось с целью изучения влияния полуширины и формы аппаратной функции приемного устройства, газового состава смеси, наличия температурных неоднородностей на результаты пирометрических измерений и выявления возможностей определения физических параметров среды по спектральным измерениям ее излучательной способности..
Излучательная способность сред моделировалась программным комплексом RADTEMP. Рассматривалась среды, представляющие собой смеси газов СО, С02 и Н20 в различных соотношениях концентраций и при различных температурах [84−86]. В расчетах использовались общепринятые допущения (диффузность излучения и наличие локального термодинамического равновесия). Проведенные расчеты позволяют выявить наиболее важные факторы, влияющие на достоверность вычислений излучательной способности газового объема (ширина и контур излучения, геометрии трассы, спектрального разрешения, температурного распределения и т. д.) и оценить адекватность экспериментальным значениям результатов расчета излучательных характеристик газовой фазы факела [65−67]..
Возможность регистрации излучения газообразных продуктов сгорания углеводородов при его распространении в атмосфере Земли оценивалась в [87]..
Как и ожидалось, столь многопараметрическая задача трудноразрешима относительно получения однозначных пространственных распределений концентраций газов. Определяющей в нашей схеме оказалась оптическая толща, являющаяся произведением концентрации газа на длину пути. На рис. 4.1. показана на примере СО зависимость радиационной температуры от значения оптической толщи — при больших значениях последней слой излучает как абсолютно черное тело (АЧТ) (на рис. 4.1. это величина 50 атм. см), что и следовало ожидать..
Анализ роли аппаратной функции, в данном случае определяемой спектральными характеристиками приемного устройства, при определении радиационной температуры проиллюстрирован на рис. 4.2.2.-4.2.3. Как следует из рис. 4.2.4, профиль радиационной температуры обнаруживает полную корреляцию со спектральной картиной коэффициента поглощения СО приведенной на рис. 3.1.2, 3.3.1, 3.3.2. Очевидно, для анализа с целью получения информации о физических параметрах среды лучше использовать данные, измеренные прибором с меньшей шириной аппаратной функции, поскольку это приводит к большей детализации спектра измеряемых радиационных температур..
Из рис. 4.2.3. следует, что наиболее близкие к термодинамической температуре слоя значения радиационной температуры получаются при прочих равных условиях при треугольной форме спектральной зависимости чувствительности регистрирующего устройства, что соответствует наиболее используемой в реальных экспериментах форме. В данном случае это связано с.
98 тем, что центр аппаратной функции был выбран в максимуме поглощения СО, и при усреднении треугольной функцией этот участок входил с максимальными весами. Опираясь на эту же зависимость, можно сказать, что для пирометрических измерений газовых сред не имеет смысла использовать приборы с шириной аппаратной функции, превышающей ширину полосы поглощения..
ЗАКЛЮЧЕНИЕ.
1. Выполнен анализ существующих оптических методов определения характеристик неоднородных газовых сред, позволяющий говорить о необходимости разработки новых алгоритмов обращения данных дистанционных измерений с целью повышения устойчивости решения к заданию модели среды и погрешностям в экспериментальных данных..
2. Разработан метод определения характеристик одномерного концентрационного поля неоднородной газовой среды по данным измерений спектра ее прозрачности:.
— предложен алгоритм обработки данных, основанный на математическом аппарате статистических моментов-.
— модифицирован алгоритм метода наименьших квадратов для применения к анализу спектра оптической толщи, нелинейно связанной с измеряемой функцией пропускания-.
— проведен вычислительный эксперимент на примере определения характеристик высотных профилей концентраций атмосферных газов, подтвердивший применимость предложенного метода к анализу высотной структуры содержания газов в атмосфере-.
— анализ влияния ошибок в задании экспериментальных данных свидетельствует о повышенной устойчивости предложенного алгоритма обращения данных..
3. Предложен способ дистанционного определения концентраций молекулярных газов в температурно-неоднородных газовых средах по измерению прозрачности на участках с минимальной температурной зависимостью коэффициента поглощения:.
— теоретически обоснована возможность минимизации температурной зависимости поглощения малоатомными молекулярными компонентами специальным выбором аппаратной функции с центром в окрестности нелокальной точки пересечения спектральных зависимостей коэффициентов поглощения при различных температурах в случае оптически тонкого слоя, реализующегося вне линий поглощения и оптически толстого слоя в центрах линий поглощения-.
— показана возможность определения in situ концентрации угарного газа в температурно-неоднородной среде, основанная на выборе аппаратной функции, компенсирующей различную температурную зависимость коэффициента поглощения вблизи центра фундаментальной полосы и на ее краю-.
— предложена номограмма, связывающая логарифм прозрачности среды с искомой концентрацией угарного газа-.
— проведен вычислительный эксперимент подтверждающий возможность восстановления концентрации СО с точностью до 10% в средах с изменением температур в пределах 400 — 1400 К и содержанием угарного газа от 0.05 до 20 атм-см (STP) —.
— предложен способ оценки средней температуры, при которой находится поглощающий газ вдоль оптической трассы..
4. Рассмотрена концентрационная зависимость регистрируемой пирометром радиационной температуры нагретого газового объема. Предложен пирометрический способ определения температуры нагретой газовой среды по ее спектральному излучению в случае неопределенной концентрации излучающего газа, основанный на сравнении радиационных температур в максимуме полосы поглощения и на ее крыле..
5. Созданы концептуальная, логическая и функциональная схемы моделирования пирометрических измерений собственного теплового излучения нагретой газовой среды, реализованные в программном комплексе RADTEMP. Проанализировано влияние возможных факторов (центра, ширины и формы аппаратной функции, наличия температурных неоднородностей, диапазона изменения концентрации) на регистрируемую.
Ill излучательную способность газового объема, для оптимизации эксперимента по определению термодинамических характеристик неоднородных газовых сред. В главу 2 вошли исследования, выполненные по соглашению с Национальным космическим агентством Японии по исследованию космического пространства (NASDA), проект № 1110 «Создание информационно-поисковой системы для анализа ослабления излучения С02-лазера атмосферными газами» (1994;2000 г.), а также поддержанные грантом Регионального конкурса РФФИ «Р 98 Сибирь (Томская область, Республика Алтай)» (рег.№ 98−01−3 003), и опубликованы в [13, 14, 62, 68−70].
В главу 3 вошли материалы, выполненные по программе «Университеты России — фундаментальные исследования» (грант per. № 990 695), и опубликованные в [63, 75]..
В главу 4 вошли исследования, поддержанные программой «Научные исследования высшей школы по экологии и рациональному природопользованию» (грант рег.№ 08.01.09), и опубликованные в [63−67].