ства. Его можно использовать при определении ранее неизвестных кристаллических структур химических соединений, полученных в поликристаллической форме.
Литература
1. David W.I.F., Shankland K. Structure determination from powder diffraction data // Acta Cryst. 2008. A64, pp. 52-64.
2. Yakimov Y.I., Semenkin E.S., Yakimov I.S. Two-level genetic algorithm for a fullprofile fitting of X-ray powder patterns //
Z. Kristallogr. Suppl. 2009. Vol. 30, pp. 21-26.
3. Якимов Я.И., Семёнкин Е.С. Автоматизированная система рентгеноструктурного анализа поликристаллов // Программные продукты и системы. 2009. № 4 (88). С. 82-84.
4. Sergienko R., Semenkin E. Competitive Cooperation for Strategy Adaptation in Coevolutionary Genetic Algorithm for Constrained Optimization // IEEE World Congress on Computational Intelligence (WCCI'2010). Barcelona, Spain, 2010, pp. 1626-1631.
5. Галушин П.В., Семёнкин Е.С. Асимптотический вероятностный генетический алгоритм // Вестн. СибГАУ им. акад. М.Ф. Решетнева. 2009. Вып. 4 (25). С. 37-42.
УДК 519.254
ПРОГРАММНЫЙ КОМПЛЕКС ОБНАРУЖЕНИЯ АНОМАЛИЙ ФОРМЫ РЕЛЬСОВЫХ ПУТЕЙ
Г.М. Соломаха, д.ф.-м.н. (Тверской государственный университет, gsolomakha@yandex.ru); В.И. Уманский, к.т.н. (ЗАО «ИнтехГеоТранс», г. Москва, umanvi@yandex.ru)
В статье представлен программный комплекс обнаружения аномалий формы рельсовых путей, реализующий статистические критерии обнаружения и оценки параметров пространственно-протяженных аномалий форм рельсовых нитей в условиях априорной неопределенности относительно их пространственного положения и времени возникновения.
Ключевые слова: программный комплекс, критерий, аномалия, рельсовая нить, моделирование.
Существует ряд зарубежных и отечественных систем контроля формы рельсовых нитей, в которых используются различные физические принципы. Как правило, качество измерений в этих системах обеспечивается специализированными и дорогостоящими передвижными лабораториями, а результатом измерений, проводимых во время специальных заездов, являются высокоточные данные о состоянии рельсовой колеи. Самостоятельное значение имеет задача оперативного мониторинга состояния железнодорожного пути, выполняемого измерительными средствами, размещенными на типовом подвижном составе. В этом случае следует говорить не о детальных высокоточных данных, а о фиксации некоторых аномальных состояний, угрожающих безопасности движения. Таким образом, контроль текущего состояния рельсовых путей (нитей) сводится, как это следует и из сущности задачи идентификации аномальных форм рельсовых нитей (АФРН), к решению задач обнаружения и оценки параметров аномальных участков.
Типовыми АФРН, обнаруживаемыми в первую очередь, являются аномалии следующих видов:
- односторонняя просадка (локальное вертикальное искривление) одной рельсовой нити относительно другой;
- одновременная и одинаковая по форме просадка (локальное вертикальное искривление) обеих рельсовых нитей;
- горизонтальное (в плане) одновременное и одинаковое искривление обеих рельсовых нитей.
Возможны комбинации этих простейших АФРН.
В результате предварительных исследований в созданной системе регистрации центральным элементом определен миниатюрный прибор ADIS 16364 фирмы Analog Device, включающий в себя по одному трехосевому датчику линейных ускорений (акселерометр) и угловых скоростей (гироскоп). Поэтому в данной работе учитываются только зависимости линейных ускорений движения вдоль каждой из трех осей декартовой системы координат и угловых скоростей поворота вокруг этих же осей.
Параметры прибора ADIS16364 были взяты из технического описания [1]. При имитации работы прибора осуществлена калибровка его шумов по результатам натурных записей с выхода прибора, полученных на движущемся локомотиве.
Для обоснования требований к системе регистрации моделирующий комплекс оснащен алгоритмами обнаружения аномалий формы рельсовых нитей и фильтрации помеховых и шумовых воздействий при оценке параметров аномалий. Моделирующий комплекс позволяет оценить требования по имитируемым и реальным потокам информации контрольно-измерительных приборов.
Кроме того, в комплексе предусмотрен режим полунатурного моделирования. В этом режиме
имитируемые данные используются совместно с реальной информацией, например, шумы в показаниях датчиков берутся реальные, а сигнал, формируемый в момент прохождения локомотивом участка пути с АФРН, - имитируемый. Необходимость полунатурного моделирования обусловлена тем, что в паспортных данных на прибор указываются лишь основные параметры датчиков, а параметры, кажущиеся разработчикам второстепенными, не приводятся (например, вид спектральной плотности помех). В результате сигналы датчиков, имитируемые в соответствии с паспортными параметрами, в определенных условиях могут существенно отличаться от реальных.
Вместе с тем моделирование по паспортным данным позволяет уже на начальном этапе разработки системы мониторинга отобрать из множества датчиков, реализующих разные физические принципы измерения, наиболее подходящие для решения поставленной задачи.
В основу моделирования геометрии рельсовых нитей взяты аналитические соотношения, приведенные в [2], используемые для описания железнодорожного полотна на местности. Наиболее важными участками пути считаются стрелочные переводы и повороты (далее моделируются только односекционные повороты). Это связано с тем, что сигналы, возникающие в датчиках линейных ускорений и угловых скоростей при прохождении локомотивом указанных участков пути, в определенной мере могут быть похожи на сигналы при наличии аномалий формы рельсовых нитей на относительно прямолинейных участках пути.
Сформулируем задачи, которые можно решить с помощью моделирующего комплекса в процессе разработки системы обнаружения АФРН.
1. Обоснование требований к системе регистрации АФРН. Проводится обоснование требований по составу, диапазону возможных значений, по точностным характеристикам, темпу выдачи-обновления информации, формируемой на выходе того или иного регистрирующего датчика о текущем состоянии рельсовых нитей на различных участках.
2. Сравнительная оценка информативности датчиков параметров движения, предлагаемых различными фирмами-изготовителями. Сравнение осуществляется путем моделирования показаний всевозможных датчиков для одних и тех же условий движения локомотива и типовых АФРН при воздействиях различного вида шумов и помех, связанных с кинематикой движения.
3. Сравнительный анализ алгоритмов обнаружения АФРН. Путем имитационного моделирования можно сформировать различные условия движения локомотива, в том числе специальные, которые могут не встретиться в ограниченных по времени и условиям натурных испытаниях. Это
позволяет всесторонне изучить сильные и слабые стороны различных подходов к выявлению АФРН.
4. Периодическая проверка работоспособности программно-алгоритмического обеспечения штатной системы обнаружения АФРН. Имитационное моделирование целесообразно использовать и для создания системы тестов, с помощью которых должна осуществляться периодическая проверка функционирования эксплуатируемой системы регистрации текущего состояния рельсовых нитей.
5. Конечная задача моделирования - оценка характеристик качества обнаружения АФРН.
Отметим, что показания датчиков линейных ускорений и угловых скоростей, установленных на локомотиве, в первую очередь зависят от геометрии пути (трехмерных координат Х(1), У®, Z(t) геометрического места точек обеих нитей железнодорожного полотна на участке моделирования), а также от функциональной зависимости У^) путевой скорости локомотива (величины модуля вектора скорости) от времени. Система координат прямоугольная, привязана к точке начала пути, ось ОУ направлена вдоль начального направления движения, ось OZ направлена вертикально вверх, а ось ОХ дополняет систему до правой. При моделировании рассчитываются показания датчиков для полностью исправного железнодорожного полотна, то есть без аномалий, а затем в показания датчиков на участках пути, где имеют место АФРН, добавляются соответствующие поправки, шумы и помехи.
В качестве базовой формы огибающей АФРН при моделировании используем гауссоиду. Выбор гауссоиды как огибающей АФРН объясняется тем обстоятельством, что, несмотря на сложную и изменяющуюся форму, которую принимает рельс в разные моменты времени в процессе движения локомотива и под воздействием нескольких точечных нагрузок (места касания колес), координаты точек касания колес практически всегда плавно смещаются вниз в аномальных зонах (например, под подвешенными шпалами) и возвращаются в исходное положение на безаномальных участках рельсовой колеи. Для АФРН также характерно наличие колебательной составляющей.
Приведем краткое описание основных элементов рассматриваемого программного комплекса, блок-схема которого представлена на рисунке.
В качестве исходных данных выступают параметры движения локомотива, характеристики приборов, входящих в систему регистрации, параметры обнаружителя, вид и параметры аномалии.
Моделирование движения локомотива.
Движение локомотива выбирается или по прямой линии, или в условиях односекционного поворота. Односекционный поворот состоит из трех характерных участков: начальная переходная кривая,
Ввод исходных данных для
моделирования
Моделирование АФРН
Моделирование обнаружения АФРН
Моделирование движения локомотива
Моделирование идеальных сигналов на выходе датчиков
1 г
Моделирование шумов и помех
1 г
Ввод данных —' натурного . / эксперимента /
Оценка точностных характеристик метода
Вывод результатов
круговой участок, конечная переходная кривая. На таком повороте поперечная компонента линейного ускорения локомотива плавно нарастает по модулю на участке начальной переходной кривой, постоянна на круговом участке, а затем плавно убывает до нуля. Отрицательность или положительность значения ускорения зависит от соотношения скорости локомотива, радиуса поворота и возвышения наружного рельса.
Моделирование аномалий формы рельсовых нитей. Локальные изменения формы рельсового пути моделируются гауссоидой по определенным направлениям в соответствии с задаваемыми исходными данными по виду АФРН и глубине просадки или бокового смещения.
Моделирование идеальных сигналов на выходе датчиков системы регистрации АФРН. Моделируется форма сигнала, возникающего в акселерометрах и гироскопах системы регистрации АФРН при отсутствии воздействия на приборы системы регистрации помех и шумов в условиях аномалий одного из трех указанных выше видов. На основе анализа результатов экспериментальных данных по формам огибающих сигналов выбрано следующее математическое описание эффектов от АФРН на выходе датчиков регистрации АФРН:
АЦ, 1) = z1je~I2i,^ , (1)
где А^, 1) - нелинейная функция компонент вектора Zj (] - номер вида аномалии) и времени (1 -номер отсчета); Zlj - амплитуда гауссоидальной аномалии вида j; Z2j - коэффициент затухания аномалии; Zзj - временная длительность аномалии или эквивалент пространственной протяженности; Z4j - несущая частота гауссоидальной аномалии типа j как гауссоидального радиоимпульса.
Моделирование шумов приборов системы регистрации. В модель заложены шумовые параметры реальных приборов, а также добавлен шум небольшой мощности, возникающий в процессе движения.
Моделирование помех. Моделируются помехи, вызванные различными причинами:
- ударными нагрузками на круги катания локомотива на стыках участков рельсовой колеи в стрелочных переводах и переходах;
- ударными нагрузками на круги катания вагонов на стыках рельсов;
- продольными толчками со стороны вагонов на локомотив при изменении скорости, рельефа местности и др.;
- резонансными квазигармоническими колебаниями в рельсовом волноводе с бегущими точечными нагрузками.
Моделирование реальных сигналов на выходе датчиков системы регистрации АФРН. На основе анализа вида шумов и помех в рассматриваемой задаче выбрана аддитивная модель влияния шумов и помех на полезный сигнал, соответствующий эффекту от АФРН на выходе соответствующих датчиков, вида y1=0jA(Zj, О+^+^+йь в которой 0j=O при отсутствии АФРН и 0j=1 при ее наличии; ^ - импульсное помеховое воздействие; й - узкополосное помеховое воздействие на прибор контроля рельсовых нитей; ^ - внутренний шум измерительного прибора; законы распределения вероятностей помех и шума известны - они установлены по статистическим данным показаний датчиков информации о возможных воздействиях на них при отсутствии аномалий рельсовых нитей; установлено также, что шум и помехи стохастически взаимно независимы.
При этом реальный сигнал является реализацией случайного процесса. Кроме того, предусмотрена возможность обработки данных натурных экспериментов, полученных рассматриваемой системой регистрации измерений, при прохождении локомотива по участкам, где заведомо имеется аномалия.
Моделирование процесса обнаружения АФРН. Решается задача определения наличия АФРН на исследуемом участке рельсового пути, а также оценки ее параметров в случае обнаружения аномалии.
Обнаружение и оценка параметров аномалий основываются на функциях правдоподобий порождения выборки измеренных датчиками значений
соответствующих параметров (линейных ускорений или угловых скоростей) при наличии или отсутствии аномалий в условиях помех и шумов, что предопределяет рисковые условия принятия решений по обнаружению, распознаванию и оцениванию параметров аномалий. Поэтому критерии обнаружения и оценки параметров АФРН следует формулировать в терминах теории статистических решений [3]. При этом для фиксированного вида аномалии (индекс ] вида аномалии опустим) выражение для отношения правдоподобия Л1 как критерия обнаружения аномалии в рекуррентной форме имеет вид:
1 Г (у, - ,1))2 )
■— ехр \--1--->
■ч/2яст I 2ст2
К =К
1
exp
y¡
2ст2
> я(а),
л/2лст
i=1, 2, n. (2)
Исходя из поступления измерений в последовательные моменты t¡, i=0, 1, 2, ..., n, когда параметры аномалии неслучайны и неизвестны, проведено осреднение по помехам и при нормальных функциях правдоподобий гипотезы 0=0 и альтернативы 0=1. С учетом всей выборки в целом эта структура сводится к следующему:
£ y¡A(z*, i) >1 (¿A2 (z*, i) + ln(rn)
2 V i = o
(3)
то есть к виду корреляционно-оценочного обнаружителя детерминированного сигнала с априори неизвестными параметрами, где ст - среднеквад-ратическое отклонение шума; z - вектор параметров аномалии, А^ , 1) имеет вид (1), а каждая выборка отсчетов уп={у1}, 1=0, 1, 2, ..., п, формируется по получаемым последовательностям измерений отдельно от каждого контрольно-измерительного прибора с помощью скользящего окна; длительность окна определяется по результатам анализа временной протяженности аномалии.
В (3) левая часть - результат раскрытия выражения для Л1, представляющий свертку выборки и сигнала от аномалии, или, иначе, согласованный фильтр, а правая - пороговый уровень, однозначно соответствующий пороговому уровню п обнаружения аномалии.
Здесь п - пороговый уровень принятия решения об обнаружении аномалии; устанавливается по допустимой вероятности ложного обнаружения аномалии (вероятности ложной тревоги) и условному закону распределения вероятностей статистики Лп при условии 0=0; вычисляется порог п по
да
выражению | у(Лп / 0 = 0)йЛ = а, где а - допус-
п
тимая вероятность ложной тревоги о наличии аномалии; у(Лп/0) - плотность распределения вероятностей статистики Лп при условии отсутствия
аномалии (0=0); эта плотность распределения определяется имитационным моделированием или аналитически.
При выполнении условия (3), то есть при обнаружении аномалии соответствующего вида, в качестве оценки ее параметров выбирается тот вектор z , для которого достигается максимальное значение в левой части неравенства (2).
Оценка точностных характеристик метода обнаружения аномалии с оценкой ее параметров. Осуществляется сравнением результатов моделирования с исходными данными об аномалии в соответствующей метрике. Данные оценки точностных характеристик являются выходными данными программно-моделирующего комплекса.
Представленный программный комплекс реализован в среде Delphi 7.
При математическом моделировании полагалось, что имеет место локальная аномалия рельсовой колеи и на выходе датчика по текущему состоянию колеи формируется последовательность измерений соответствующего признака. Эта последовательность подается на вход обнаружителя.
Например, для прямолинейных участков движения поезда «Сапсан» (250 км/час) результаты математического моделирования обнаружения аномалий представлены в таблице.
Характер сигнала в случае аномалий
Измеряемый датчиком признак Тип аномалии
Односторонняя просадка Двухсторонняя просадка Рыскание (отклонение пути в плане)
Амплитуда, мм
8 4 8 4 8 4
Линейное ускорение по оси ОХ СС СС
по оси ОУ - - - - - -
по оси О7 - СС СС - -
Угловая скорость поворота вокруг оси ОХ
вокруг оси ОУ ХС ПС - - - -
вокруг оси О7 - - - - СС ХС
Примечание: СС - сильный сигнал, ПС - пороговый сигнал, ХС - хороший сигнал.
Сильными считаются сигналы, обнаруживаемые практически достоверно, то есть вероятность пропуска и принятия решения о наличии аномалии из-за шумов по сути равна нулю.
Хорошими считаются сигналы, по которым аномалия обнаруживается путем согласованной фильтрации (3) и пороговой обработки с высокой вероятностью. Амплитуда данных сигналов превышает 3 -4 значения среднеквадратического отклонения шума.
И наконец, пороговыми сигналами от аномалий считаются те, амплитуда которых соизмерима
со среднеквадратическим отклонением шума. Для обнаружения аномалий в случае пороговых сигналов можно воспользоваться критерием проверки близких гипотез [4].
В результате можно сделать следующие выводы. На основе анализа причин и условий возникновения аномалий форм рельсовых нитей показано, что построение метода их обнаружения должно осуществляться по критерию совместного обнаружения и оценивания параметров альтернативных гипотез на различных участках движения локомотива при аддитивных воздействиях импульсных, узкополосных помех и гауссовых шумов с априори неизвестными мощностями. Критерий обнаружения и оценивания параметров аномалий представляется в виде отношения функ-
ций правдоподобия, которое максимизируется на множестве их альтернативных типов.
Представленный программный комплекс позволяет исследовать процессы обнаружения аномалий формы рельсового пути для данных, полученных имитационным моделированием или натурным экспериментом.
Литература
1. Six Degrees of Freedom Inertial Sensor ADIS 16364. Analog Devices, Inc. Norwood, MA 02062-9106, U.S.A., 2009.
2. Виноградов В.В. Расчеты и проектирование железнодорожного пути. М.: Маршрут, 2003.
3. Боровков А.А. Математическая статистика. М.: Лань, 2010.
4. Кудинов А.Н. [и др.]. Интегральный оператор дифференцирования случайных полей // Радиотехника (Журнал в журнале). 2008. № 14.
УДК 004.415.2
ПРОГРАММНЫЙ КОМПЛЕКС ОЦЕНКИ ХИМИЧЕСКОЙ ОБСТАНОВКИ ПРИ ВОЗНИКНОВЕНИИ ЧРЕЗВЫЧАЙНЫХ СИТУАЦИЙ
Ю.Н. Матвеев, к.т.н.
(Тверской государственный технический университет, fmas@tstu.tver.ru)
В статье представлены программный комплекс оценки химической обстановки при возникновении чрезвычайных ситуаций на объектах уничтожения химического оружия, интерфейсы системного администратора и конечного пользователя.
Ключевые слова: аварийная ситуация, чрезвычайная ситуация, принятие решений, химическое оружие, техническая безопасность, имитационное моделирование, интерфейс, семантика.
Программный комплекс оценки химической обстановки предназначен для обоснования технических решений по ликвидации последствий возможных чрезвычайных ситуаций (ЧС) при за-проектных авариях на объектах хранения и уничтожения химического оружия (ОХУ ХО) [1]. Комплекс обеспечивает выполнение следующих функций: ввод данных о запроектной аварии на объекте, оценка возможности возникновения ЧС и ее последствий, формирование отчетов, отображение оценки опасной зоны загрязнения на карте.
Процедура принятия решения о возможном развитии аварийной ситуации в ЧС заключается в следующем. Информация об аварии, содержащая данные о номере аварийного аппарата в соответствии с технологической схемой, о типе аварии, поступает в БД автоматизированной системы поддержки принятия решений (АСППР). Определяются координаты места аварии, тип аварийного химически опасного вещества (АХОВ) и его максимальная масса в каждом аппарате, модель процессов образования и распространения облака зараженного воздуха, необходимые константы и характеристики топографии местности. Эта ин-
формация, а также метеоданные и математические модели необходимы для расчета площади зоны возможного заражения местности. Величина оценки зоны возможного заражения и координаты границы санитарной защитной зоны поступают в анализатор ЧС, где определяется возможность развития аварийной ситуации в ЧС. В случае выхода зоны заражения за пределы санитарной защитной зоны производится оповещение сил, участвующих в локализации ЧС, персонала объекта, ответственных лиц и учреждений, осуществляется запуск автоматизированной системы оперативного управления в ЧС. В противном случае оповещаются службы безопасности ОХУ ХО о результатах моделирования аварийной ситуации. Отображение химической обстановки на топографической карте местности осуществляется постоянно с момента начала расчета зоны заражения с предельно допустимой концентрацией на границе зоны. Алгоритм оперативного управления в условиях ЧС [2], представленный на рисунке 1, работает следующим образом.
Информация о возможной ЧС поступает из системы анализа развития аварийной ситуации в ЧС. При необходимости может быть введена ин-