Оценка степени влияния природных и технологических факторов на различные типы аварий магистральных газопроводов

(блок 6) и алгоритм формирования мероприятий по ликвидации последствий ЧС (блок 7). Функ-

ционирование общего алгоритма обеспечивает анализ возможности возникновения ЧС при за-проектной аварии на ОХ и УХО, а также поддержку принятия решения должностными лицами ОХ и УХО в ЧС, а именно:

- по управлению силами и средствами при ликвидации последствий аварии на ОХ и УХО (локализация и ликвидация источника, ликвидация последствий);

- по управлению силами и средствами при ликвидации последствий ЧС.

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

В данном случае следует уточнить термин «Авария запроектная», который трактуется как событие, порожденное какими-либо внешними факторами (стихийные бедствия и катастрофы: ураган, землетрясение, падение летательного аппарата и т.д.), приведшими к нарушению нормального режима функционирования объекта и повлекшими за собой выделение токсичных веществ в окружающую среду [2].

В заключение отметим, что алгоритм реализован в виде программных модулей ПТК «КРАБ» для объекта «Горный» (Саратовская обл.).

Литература

1. Батырев В.В. [и др.]. Обоснование размеров зоны защитных мероприятий вокруг объектов по хранению и уничтожению химического оружия: отчет по НИР «Вагонетка-О». Этап 2. Новогорск: АГЗ, 1998. 214 с.

2. Нормы специального проектирования объектов 1281, 1282, 1596, 1597, 1721, 1728, 1729 по уничтожению химического оружия. НСПО1-99/МО РФ. Изд. официальное (утв. НВ РХБЗ МО РФ 17.06.1999 г.). М., 1999. 60 с.

УДК 656.56.338: 656.56.658

ОЦЕНКА СТЕПЕНИ ВЛИЯНИЯ ПРИРОДНЫХ И ТЕХНОЛОГИЧЕСКИХ ФАКТОРОВ НА РАЗЛИЧНЫЕ ТИПЫ АВАРИЙ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ

А.С. Жуков

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

Ключевые слова: трубопроводы, распознавание образов, коррозия, теория вероятностей, прогнозирование, факторы, аварии, методика.

Существующая в России Единая газотранспортная система (ЕСГ) - структурно сложная

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

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

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

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

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

1) выявление в системе элементов и подсистем, находящихся в наихудшем техническом состоянии в сочетании с наиболее тяжелыми условиями работы, повышенной важностью этих элементов;

2) планирование мероприятий, связанных с ремонтом и реконструкцией;

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

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

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

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

Для обеспечения крупномасштабного (регионального или по ЕСГ в целом) прогноза технического состояния и аварийности газопроводов в условиях изменения нагрузок, активных техногенных и природно-климатических воздействий на систему, физического старения системы и прочего прежде всего необходимо обеспечить возможность автоматизированного прогнозирования аварийности магистральных газопроводов по результатам анализа ретроспективных данных об аварийности.

Поставленную задачу целесообразно решать для газотранспортной системы в целом, однако она имеет смысл и для укрупненных фрагментов системы. Для таких подсистем можно вычислить фактические (имеющие место, например, с 1991 по 2009 гг.) удельные показатели аварийности и с учетом динамики этого показателя по годам экстраполировать их на ближайший период (например, до 2012-2015 гг.).

Затем, исходя из существующей (по наблюдениям) структуры и прогнозируемых темпов динамики аварийности, можно установить, например, путем экспертного анализа, на какую долю прогнозируемых аварий повлияют те или иные планируемые программой ремонтные мероприятия, предполагаемые меры по реконструкции ЕСГ, темпы физического старения системы, масштабы распространения стресс-коррозии и другие крупномасштабные воздействия на газотранспортную систему.

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

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

Ввиду большого разнообразия наблюдаемых и подлежащих распознаванию объектов их обычно объединяют абстрактным термином «образ» и говорят о «распознавании образов». Каждая задача распознавания сугубо индивидуальна и по-разному решается в различных областях науки и техники.

Постановку задачи распознавания применительно к предсказанию видов аварий трубопроводов реализуем на простейшем условном примере.

Пусть при расследовании аварии трубопровода 25-летнего возраста, работающего на давление 7,5 МПа, по виду разрушения невозможно установить, что явилось ее причиной: обычная (классическая) коррозия или стресс-коррозия. Для упрощения задачи рассматривается только один сопутствующий признак - запись о величине рабочего давления в момент аварии.

Индексом «1» обозначим событие, заключающееся в реализации аварии по причине обычной коррозии, индексом «2» - событие, заключающееся в реализации аварии по причине стресс-коррозии.

Результаты появления образов разных классов представляют собой несовместимые и образующие полную группу события. Напомним, что конечная совокупность событий называется полной группой событий, если хотя бы одно из них обязательно появляется в результате опыта. В данном случае это реализация аварии по какой-либо причине из группы причин, образующих полную группу (коррозия, стресс-коррозия, СМР, брак металла, механическое повреждение, стихийное бедствие, диверсия).

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

Кроме записи о величине рабочего давления при разрыве, есть еще наблюдения, накопившиеся за время эксплуатации трубопроводов данного региона, работающих на давление 7,5 МПа:

• о количестве аварий по причине коррозии (3 аварии) и стресс-коррозии (55 аварий) и, соответственно, о соотношении (частости) аварий того и другого видов:

Pi - частость (априорная безусловная вероятность) коррозионных аварий (0,052),

Р2 - частость (априорная безусловная вероятность) стресс-коррозионных аварий (0,948);

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

Р(р/1) - частость (апостериорная условная вероятность) любого значения давления при наблюдавшихся коррозионных авариях;

Р(р/2) - частость (апостериорная условная вероятность) любого значения давления при наблюдавшихся стресс-коррозионных авариях.

На рисунке показаны плотности соответствующих распределений 1"(р/1) и 1"(р/2), полученные в результате реальной обработки сведений о значениях давления при авариях по причине коррозии (слева на рисунке) и стресс-коррозии (справа на рисунке) газопроводов. Обработка произведена распределением Крицкого-Менкеля, проверка -по критерию Колмогорова-Смирнова. Точка пересечения распределений - 6,75 МПа. Были рассмотрены все случаи аварий по указанным причинам газопроводов с рабочим давлением 7,5 МПа.

С помощью этой дополнительной информации необходимо принять решение о коррозионном или стресс-коррозионном происхождении аварии. Иначе по конкретной реализации признака (значения рабочего давления в момент конкретной аварии) с привлечением статических наблюдений за аварийностью и распределением давлений при авариях обоих типов требовалось найти решающее правило для отнесения аварии к 1 -му или 2-му типу, то есть найти модель распознавания.

Модель процесса распознавания представляет собой зависимость оценки 10, 20, ... номера класса образов 1, 2, ... от конкретного значения х наблюдаемой случайной величины X. Обычно наблюдаемая величина X представляет собой конечномерный случайный вектор, компонентами которого являются различные дискретные и непрерывные распределения случайных факторов (класс материала трубы, параметры режима эксплуатации, свойства грунтов, состояние изоляционного покрытия, тип изоляции и др.). Случайный характер наблюдаемых величин X приводит к тому, что одно и то же значение х может наблюдаться при разных классах образов (разных значениях номеров классов (1, 2, ..., п)). Построение модели распознавания сводится к нахождению этой зависимости, которую следует определить так, чтобы число ошибок при применении модели распознавания было как можно меньше.

Такое решающее правило в теории распознавания образов определяется соотношением Р1 АрЛ)=Р2 1"(р/2).

Если области значений признака графически изобразить в виде кривых распределения Р(А1)= =P1 f(p/l) и Р(А2)=Р2 f(p/2), получим, что оптимальная граница разбиения этих областей (значение давления р0) находится на пересечении кривых распределения, а вероятность ошибки отнесения аварии к несоответствующему классу измеряется площадью заштрихованного криволинейного треугольника (см. рис.). При значении р0 площадь заштрихованного треугольника наименьшая. Если граница областей р0 смещается влево или вправо, вероятность ошибки увеличивается в зависимости от площади заштрихованных криволинейных треугольников ABC или CDE.

Оптимальная модель распознавания вырабатывает конкретное (опытное) значение 1оп или 2оп по следующему правилу:

. f(P|i) р2 „ f(p|i) , р2

1 пп - п р и —— > — , 2 пп - при —— < — .

оп V fp2 Pi' оп F f(p|2) Pi

Отношение

f

называется отношением

правдоподобия. Отсюда видно, что решающую роль в построении модели распознавания играет отношение условных плотностей распределения признаков, в данном случае - плотностей распределения признака р, то есть ^р/1) и ^р/2).

Определяя по имеющемуся значению давления в момент конкретной аварии ординату плотности распределения ^рк/1) и ординату плотности ^рк/2), получаем конкретные результаты распознавания, то есть решение 1оп или 2оп.

Например, если по записи о величине рабочего давления в момент конкретной аварии получено значение Рав=7,4 МПа, то описанная процедура распознавания относит эту аварию к стресс-коррозионным, и т.д.

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

Г(р|2) Р!

Это дает возможность перейти к линейной операции суммирования плотностей распределения случайных величин, используя следующее правило:

1оп -1п f(p/l) + (- 1п f(p/2)) > (1п Р2 - 1п Р1),

2оп -1п f(p/l) + (- 1п А(р/2))< (1п Р2 - 1п Р1).

Охарактеризованная выше оптимальная модель, гарантирующая максимум правильного распознавания, работает по критерию максимума апостериорной вероятности события (появления того или иного номера класса образов при реализации случайного вектора признаков), поскольку произведения Р1 f(x|i) представляют собой переменную часть (делимое) в формуле вычисления апостериорных вероятностей появления образов

Р1

разных классов: Р (Y = i | х) =

2ÍLiP¡f(x I i)'

Поэтому при произвольном числе классов образов области оптимального решения определяются условием A k = (х : max ¡Pjf (х | i) = Pk f (х | k) }, k = 1 ,. . ., N , а вероятность правильного решения -условием

Соответственно скорректируются решающее правило и отношение правдоподобия. Графически это будет означать, что оптимальная граница (точка х0) определяется на линии пересечения распределений max Pi f(x|i) и Pk f(x|k).

Принципиально не меняет дело и то, что наблюдаемая величина X представляет собой конечномерный случайный вектор. Предположив, что все рассматриваемые признаки (составляющие случайного вектора Х19 Х2, ..., Хп) являются независимыми (а признаки для целей распознавания следует стараться выбирать с учетом фактора их взаимной независимости), определим совместную плотность распределения случайного вектора X по формуле f(X)=f(X), f(X2), ..., (Xn).

Последующее решение задачи соответствует вышеприведенной методике.

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

Поэтому организованная указанным способом машина распознавания помогает в выработке оптимальных стратегий (оптимальной очередности) капитальных ремонтов и реконструкции и в определении вероятности аварии.

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

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

Показанная на простейшем примере методика прогнозирования аварийности трубопроводов обосновывает возможность и перспективность

применения теории распознавания образов для решения описанной задачи.

Литература

1. Дедешко В.Н. Техническое состояние магистральных трубопроводов РАО «Газпром» и организация работ по внут-

ритрубной диагностике // Диагностика-98: сб. тр. Между-нар. деловой встречи. М.: ИРЦ Газпром, 1998. С. 3-31.

2. Стратегия развития газовой промышленности России; под общ. редак. Р.И. Вяхирева и А.А. Макарова. М.: Энерго-атомиздат, 1997.

3. Горелик А.Л., Скрипкин В.А. Об одном методе решения задачи классификации объектов и явлений. Изв. АН СССР // Техническая кибернетика. М., 1965. № 1.

УДК 004.946

ПРОБЛЕМЫ МОДЕЛИРОВАНИЯ ДЕЙСТВИЙ ОПЕРАТОРА В СОВРЕМЕННОМ ТРЕНАЖЕРЕ

Д.С. Ситалов (Донской филиал Центра тренажеростроения, г. Новочеркасск,,

сИтад;Ьз@гатЫег. ги)

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

Ключевые слова: тренажер, моделирование движений человека, геометрическое моделирование человека, скелетная анимация.

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

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

Модель скелета аватара представляет собой иерархически связанные в древовидную структуру элементы скелета. Здесь под элементом скелета

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