УДК 543.5/543.26
Белозерцев1 А.И., Эль-Салим2 С.З.
1АО «Научно-исследовательский институт физических измерений» 2ООО «Омега», Санкт-Петербург, Россия
ДЕТЕРМИНИРОВАННАЯ МОДЕЛЬ ПОВЫШЕНИЯ НАДЕЖНОСТИ АНАЛИТИЧЕСКИХ СИСТЕМ
Рассмотрена модель повышения достоверности, а следовательно, надежности аналитических измерений паров и газов. Исследованы факторы, связанные с результативным показателем, выраженным математической зависимостью.
Для задач, которые решаются в масштабе реального времени и не требуют высокой точности результатов моделирования, построена детерминированная модель, позволяющая уточнить длительность реальной экспозиции для повышения точности количественных расчетов
Ключевые слова:
модель, энтропия, детерминированные процессы, причина, следствие
Одной из важных задач теории надежности является активизация исследований по разработке не столько математических, сколько физических моделей отказов и причин, вызывающих их возникновение.
Следовательно, необходимо разработать такие физические или математические модели и требования, которые позволяют достаточно полно характеризовать надёжность системы независимо от ее природы. При этом необходимо знать причину, а не следствие отказа. Следовательно, необходим инженерный анализ и расчёт надёжности будущего изделия, а также алгоритма его работы.
Существующие модели, основанные на вероятностном подходе, не удовлетворяют условиям контроля готовой продукции, поскольку далеки от практики производства. Многочисленные исследования говорят о том, что в основе обеспечения надёжности лежат физико-химические процессы, протекающие в используемых материалах на стадии изготовления и последующей их эксплуатации [1, 2].
Поскольку физические явления описываются физическими теориями, на базе которых строятся математические модели, то и надёжность как процесс потери или сохранения работоспособности объекта за счет изменения внутреннего состояния описывается физической теорией, а затем уже возможно создание математической теории.
Вероятностно-статистический подход, основанный на теории вероятностей и математической статистике несовершенен, так как для установления закономерностей возникновения отказов необходимо проводить опыты, приводящие к изменению свойств и разрушению деталей. Кроме того, он не дает точного анализа поведения объекта при дальнейшей эксплуатации. В соответствии с этим подходом возможно с определенной степенью достоверности предсказать момент наступления отказа, а не предупредить его и обеспечить надёжность работы объекта.
Детерминированный (причинно-следственный)
подход позволяет оценить надёжность как отдельного элемента системы, так и всей системы в целом при одновременном воздействии на объект внешних и внутренних факторов. Поэтому такой подход к изучению надёжности является более приемлемым [3]. Детерминированная модель повышения надежности аналитических систем на основе причинно-следственных связей, прежде всего, выстраивается на анализе процессов, эту связь формирующих [4].
Существование «чистых» паров или газов в воздухе практически недостижимо в силу природного молекулярного многообразия. Молекулярный состав окружающего воздуха формируется молекулами (веществами), которые, в свою очередь, претерпевают физические и химические трансформации. Молекулы целевого (исследуемого) вещества могут сохранить структуру, или трансформироваться, сохраняя основные свойства. Также целевые вещества могут вступать во взаимодействие с другими компонентами, содержащимися в воздухе. При этом исходная концентрация целевого вещества может изменяться Для точного определения момента начала взаимодействия сенсора с газом необходимо установить параметр, который однозначно откликается на факт взаимодействия вещества газовой фазы с поверхностью. Так как для изменения электропроводности
при формировании аналитического сигнала необходима активация, в данном случае - подведение энергии для нагрева до определенной температуры, то целесообразно по измеренным и преобразованным параметрам рассчитать степень упорядоченности гетерогенной зоны.
в ту или иную сторону. То есть, состав окружающего воздуха можно считать динамически переменным. Также необходимо учитывать аэродинамическое перемешивание и способ (метод) отбора проб [7].
В основе адсорбционно-кинетического анализа [полупроводники-сенсоры] лежит измерение изменения электропроводности вследствие хемосорбции вещества из газовой фазы за определенное время. Очевидно, что даже при низких концентрациях число молекул велико. Так, количество молекул водорода в одном кубическом метре при нормальных условиях для концентрации 10-3 г/м3 равно Ы = -^^3-Ю20, где NAv - число Авогадро.
Практика изучения случайных процессов показывает, что хотя результаты отдельных наблюдений, проведенных в одинаковых условиях, могут сильно отличаться, в то же время средние результаты для достаточно большого числа наблюдений устойчивы и слабо зависят от результатов отдельных наблюдений.
Для построения феноменологической модели формирования аналитического сигнала необходимо построить цепь причинно-следственных связей и определить уравнения изменения состояния в соответствии со временем течения всего процесса измерения [5]. На рисунке 1 приведена общая схема событий в зависимости от их наступления.
Причинно-следственная модель состоит из двух практически независимых физически, но связанных идеологически блоков. К физическим блокам следует отнести:
- процессы, происходящие в объеме полупроводника при его нагреве;
- процессы, протекающие в гетерогенной зоне;
- конкурирующие процессы адсорбции и десорбции при контакте молекул газовой фазы с поверхностью.
По времени процессы располагаются по приоритету причинности:
1. Отбор и доставка вещества в гетерогенную зону газочувствительного слоя;
2. Взаимодействие молекул газовой фазы с поверхностью полупроводника;
3. Транспорт электронов адсорбата с учетом знака направления перехода;
4. Поверхностные процессы, формирующие значения изменения электропроводности.
Для повышения достоверности качественных и количественных измерений, проводимых с помощью полупроводниковых сенсоров, имеющих развитую структуру поверхности (рисунок 2), крайне важно определить момент начала взаимодействия вещества газовой фазы с поверхностным слоем газочувствительного элемента [7].
Так как данный метод анализа является косвенным, то перед формированием аналитического сигнала проводится первичная обработка, которая включает процедуру сглаживания, расчета текущего среднего значения и расчета нормированной дис-
персии
г(1) = [и(1)-(и(1))] .
Рисунок
сигнала
Рисунок 2 - Растровое изображение чувствительного слоя полупроводникового сенсора с увеличением
х50000
Так как в процессе изменения электропроводности (сопротивления) участвует большое количество частиц, сигнал П(Ь ,) распределен нормально. В этом случае интегральная функция и функция плотности вероятности имеют вид:
iy(t) =
1
V2na 1
1 (u-y)2
dt и
.К
2 ст2
где, ^ - мат. ожидание, а - стандартное отклонение величины и(1) на интервале [11, 12].
Очевидно, что приращения сигнала Д /(¿) = /(£2) — //(¿1), ¿2 > ^также нормально распределены с ^ — 0, и интегральное распределение определится как 1 г*2 .
P(AU < U) = Г'2 ехр (- ±^4) dût.
^ г V 2 Д
В этом случае полная дисперсия сигнала и(1) на интервале [11, 12] рассчитывается как О[/(^) — /(¿1)] = а2^2 — для любых значений внутри интервала. Видно, что рассчитанная дисперсия зависит только от разности Д и не зависит от самих значений сигнала и(1).
Для дисперсии и стандартного отклонения при-
lim D(»fe)-»fa)) = g
Л
At
ращений справедливо равенство
то есть, при Д£ = 1, О = а2. При изменении приращения стандартное отклонение флуктуирует относительно Б, но с бесконечно малым порядком малости: |О — а2|<£ — 0. В этом случае справедливо уравнение для мат. ожидания ¿П/^) —/(¿01] =
то есть приращение сигнала, соответствующего модели броуновского движения обладает свойством статистического самоподобия:
//(£ + ДО — /7(0 ^ -1 (/(£ + у ДО — У(0), у ем+.
Самоподобие сигнала и(1) и его эквивалентность модели броуновского движения позволяет рассчитать размерность результатов измерений на интервале [11, 12].
Стандартная процедура расчета фрактальной размерности определяется следующим образом: количество элементарных зон Ы, полностью «покрывающих» зависимость и(1), пропорционально М(Д 0 ос
1п(МД 0 ш
Л i 2 или dim U = — lim ■
В зависимости от
ДС-0 ^ДС 2у ш
частоты и у = фрактальная размерность сигнала
и(1) изменяется как при этом 8 определяется
характером взаимодействия вещества из газовой фазы с газочувствительным слоем. При этом большую роль играет состав материала, тип проводимости и мощность нагрева [2,5,6].
Таким образом, процесс изменения электропроводности вследствие хемосорбции, то есть обмен электронами между зоной проводимости и газовой фазой в гетерогенной зоне, является самоподобным по отношению к измеряемому падению напряжения в цепи газочувствительного слоя на интервале [11, 12].
Флуктуации падения напряжения в результате хемосорбции, соответствующие модели броуновского движения с точностью 8, предполагают некоторые дополнительные преобразования сигнала и(1) для повышения стабильности (однозначности) обработки результатов измерения для более точного определения аналитических параметров. Для аналитического применения полупроводниковых
2
а
сенсоров необходимо разработать модель, а на ее основе - алгоритм обработки сигнала и(1), позволяющий реализовать прямую аналитическую задачу - идентификацию вещества и определение его концентрации по измеренным значениям измеряемого напряжения.
Энтропия по разбросу значений состояния определяется как 5 = 21п2, полагая, что Z соответствует сумме предыдущих значений. В этом случае
изменение энтропии по времени равно
— = Ы! —,
& г
и, следовательно, производство энтропии:
Анализ производства энтропии в соответствии с теоремой Пригожина показывает: если 1 < 1 —1 — 21п^>0,и, следовательно, система находится в хаотичном состоянии, обмен электронами со средой (с газовой фазой) минимален. В этом случае потенциал, вызванный хемосорбцией, определяется только молекулами кислорода. Если же 2>1 —1 — 2Ы2 < 0, то обмен электронами со средой упорядочен, имеет градиент потока и обусловлен хемо-сорбцией молекул примесей, которые содержатся в газовой фазе [7].
Действительно, расчеты средних значений и дисперсий позволяют рассчитать энтропию динамического состояния и функцию динамического изме-¿г
нения энтропии —, являющуюся определяющей как
Измерения, приведенные в таблице 1, показывают, что при расхождении измеряемого сигнала (наличие мешающих примесей, наличие флуктуаций на поверхности чувствительного слоя) расчет текущих средних значений, а затем мгновенной дисперсии переводит процесс изменения электропроводности вследствие хемосорбции окислителя и восстановителя (А и НДМГ соответственно) из случайного в детерминированный. Расчет динамики и динамического изменения энтропии определяет момент начала реакции - истинной экспозиции. При подаче пробы в измерительный тракт инертность аналитической системы, определяется длиной газоотборной трубки, пульсациями побудителя расхода, температурой, влажностью, плотностью и давлением окружающего воздуха, поэтому аппаратно установленная экспозиция вносит погрешность в расчет концентраций. Расчет динамики энтропии 1п2 и динамического изменения — позволяет доста-
2
точно точно определить момент начала взаимодействия, соответственно время экспозиции уточняется и повышается точность расчета концентраций.
Таким образом, самоподобные преобразования позволяют повысить стабильность исходных измерений, уточнить время аппаратно установленной экспозиции, тем самым повысить точность количественного определения по величине дисперсии потенциала, обусловленного хемосорбцией вещества-аналита.
для скорости изменения порядка, так и для определения производства энтропии (таблица 1).
Первичные преобразования и динамика энтропии взаимодействия
Таблица 1
Измерение напряжения при хемосорбции НДМГ и АТ Концентрация НДМГ: 1 и 3,3 мг/м3, АТ: 3 мг/м3
Накопление средних значений за время экспозиции
Динамика дисперсии за время экспозиции
Динамика энтропии
В основе детерминированного моделирования лежит возможность построения тождественного преобразования для исходных данных по теоретически предполагаемым прямым связям с факторами. Детерминированное моделирование факторных систем -
это простое и эффективное средство для количественной оценки роли отдельных факторов в динамике изменения обобщающего показателя (рисунок 3).
Рисунок 3 - Переход от вероятностных моделей: равномерного распределения 1 и нормального распределения 2, 3 к детерминированной модели 4
При детерминированном анализе модель изучаемого явления не изменяется по изучаемым объектам и периодам. При необходимости сравнения результатов измерений в отдельные периоды возникает лишь вопрос о сопоставимости выявленных на основе модели количественных аналитических результатов.
Детерминированный анализ представляет собой методику исследования влияния факторов, связь которых с результативным показателем носит функциональный характер, то есть может быть выражен математической зависимостью.
Следовательно, предположение о детерминированном характере физической системы и описание
ее детерминированной математической моделью является идеализацией реальной системы. Детерминированная математическая модель отображает реальные физические процессы в усредненном смысле.
В тех задачах, где не требуется высокая точность результатов моделирования, предпочтение отдается детерминированной модели. Это объясняется тем, что реализация и анализ детерминированной математической модели много проще, чем стохастической.
Детерминированная модель недопустима в ситуациях, когда случайные процессы со^) соизмеримы с детерминированными х^). Результаты, полученные с помощью детерминированной математической модели, будут неадекватными реальным процессам.
ЛИТЕРАТУРА
1. Миллер Б.М., Панков А.Р. Случайные процессы в примерах и задачах. - М.: Изд-во МАИ. - 2001.
2. Вентцель А.Д. Курс теории случайных процессов. - М.: Наука. - 1975.
3. Вентцель Е.С., Овчаров Л.А. Теория случайных процессов и ее инженерные приложения. -М.: Наука. - 1991.
4. Обухова Л.В., Молдовская З.Я., Князева В.Ф. Теория вероятностей, математическая статистика и случайные процессы в примерах и задачах. - Киев: УМКВО. -1991.
5. Бендат Дж., Пирсол А. Прикладной анализ случайных данных. - М.: Мир. - 1989.
6. Тюрин Ю.Н., Макаров А.А. Стохастический анализ данных на компьютере. - М.: Инфра. -1997.
7. Громов В.Ф., Герасимов Г.Н., Белышева Т.В., Трахтенберг Л.И. Механизмы сенсорного эффекта в кондуктометрических датчиках на основе диоксида олова для детектирования газов-восстановителей // Ж. Рос.хим. об-ва им. Д.И. Менделеева, вып. 5. - 2008. — С. 80—87.
8. Ван-Желен В. Физическая теория надежности. - Симферополь: Изд-во «Крым». - 1998.
УДК 681.586.72
Николаев1 А.В., Эль-Салим2 С.З.
1АО «Научно-исследовательский институт физических измерений», Пенза, Россия 2ООО «Омега», Санкт-Петербург, Россия
КЛАССИФИКАЦИЯ СИГНАЛОВ В ПРОЦЕССЕ АНАЛИЗА ГАЗОВЫХ СМЕСЕЙ
В работе рассмотрены преобразования динамических измерений, их классификация и разработка алгоритма обработки кинетических данных. Процесс первичной обработки представлен в виде ряда процедур с операциями, принадлежащими пространству рациональных чисел. Показаны процедуры фильтрации данных, корректировка базовой линии, шумоподавление, предварительная статистическая оценка измеряемых параметров. Определена необходимость выбора весовых функций линейных операторов над полем обобщенных значений результатов измерений Ключевые слова:
алгоритм, суперпозиция, функционал Гаусса, детерминированные процессы, модель Лоренца и Гаусса