Информационно-аналитическая система прогнозирования процессов потребления природного газа в газотранспортной системе Украины Текст научной статьи по специальности «Математика»

УДК519.7:004+519.23

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

ТЕВЯШЕВ А.Д., ВЫХОДЦЕВ Е.И.,

ЩЕЛКАЛИНВ.Н., ИЕНАТОВА Ю.В.______________

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

Введение

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

Объектом автоматизации является служба режимов объединенного диспетчерского управления (ОДУ) дочерней компании (ДК) "У кртрансгаз. Одной из ее функций является обеспечение автоматизированного сбора, обработки, хранения и вычисления прогнозных значений суточных объёмов потребления природного газа по каждой газораспределительной станции (ГРС), находящейся в зоне ответственности ДК ”У кртрансгаз”.

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

1. Структура системы

Информационно-аналитическая система «Прогнозирование процессов потребления природного газа всеми внутренними потребителями ГТС Украины» (ПАС ППППГ) разработана как информационно-аналитический сервис, позволяющий службе режимов более эффективно и обоснованно принимать решения. В силу специфики функций и предложенной архитекту-рыИАС ППППГ обрабатываетинформацию из различных источников в автоматическом режиме и проводит расчеты без вмешательств пользователей системы.

На основании результатов системного анализа были определены основные функции системы:

- получение суточных фактических и плановых объемов потребления природного газа по каждому выходу ГРС;

- получение среднесуточной температуры воздуха и ее краткосрочных прогнозов в области нагрузки каждой ГРС;

- расчет прогнозных объемов суточного газопотребления по каждому выходу ГРС с разбивкой по выделенным категориям потребителей с учетом организационных, хронологических и метеорологических факторов в области нагрузки каждой ГРС;

- формирование информационных отчетов о расчётах суточных прогнозных объемов газопотребеления по различным категориям потребителей, областям и районам Украины, отдельно по каждому производственному отделению и в целом по ДК «Укртрансгаз».

Разработана модель бизнес-процесса «Прогнозирование процессов потребления природного газа внутренними потребителями ГТС У краины», связанного с решением задачи прогнозирования, которая представлена на рис. 1 в нотации BPMN [1].

Внешними системами к ПАС ППППГ являются электронный журнал диспетчера ОДУ ДК «Укртрансгаз», программный комплекс «Программа годовой отчетности ЛЧ МГ и ГРС» и сервис прогноза погоды Украинского гидрометеорологического центра (УкрГМЦ). На рис. 2 приведена схема потоков данных в разработанной системе.

ИАС ППППГ разработана с применением принципов сервис-ориентированной архитектуры (SOA), которая позволяет реализовать систему в виде модулей со слабой связанностью между компонентами системы. На рис. 3 приведена архитектура SOA-решения для ИАС ППППГ.

Процедура получения внешних данных заключается в программной реализации двух бизнес-правил: «Получение суточных фактических и плановых объемов потребления природного газа по выходам ГРС» и «Получение фактических значений температуры воздуха и её краткосрочного прогноза в зоне нагрузки ГРС». Дляреализации этих бизнес-правил разработаны два специализированных Window s Communication Foundation (WCF) сервиса GCF. WcfGdsDataCollector,

РИ, 2011, №3

92

о

Инициализация

параметров

системы,

восстановление

состояния

потока

выполнения

Запуск служб сбора данных

Получение суточных фактических и плановых объемов потребления природного газа по выходам ГРС

->

Получение среднесуточных фактических и краткосрочных прогнозных значений температуры воздуха

->

Суточные объемы потребления природного газа по каждому выходу ГРС

Данные прог * в целом по ДК \ * по производств* подразделения * по областям и р Украины Иилоь. кртрансгаз энным айонам

У

Среднесуточная темп-pa воздуха в зоне нагрузки ГРС

Анализ

данных

*

Ч/

*

____J

Вычисление прогнозов и доверительных границ объемов потребления природного газа для каждого выхода ГРС

W

і

х

Ж

4—

Регистрация событий и операций системы

Журнал протоколирования операций системы

Рис. 1. Модель бизнес-процесса «Прогнозирование процессов потребления природного газа внутренними потребителями ГТС Украины»

Рис. 2. Модель потоков данных «Прогнозирование процессов потребления природного газа внутренними потребителями ГТС Украины»

GCF.WcfWeatherCollector, которые управляются интеллектуальным агентом GCF.Tasks. GCF.Tasks реализовано в виде WCF сервиса и использует планировщик заданий Quartz.NET (http:// auartznet.sourceforge.net).

Основными функциями сервиса получения среднесуточных значений температуры воздуха и её краткосрочных прогнозов в зоне нагрузки каждой ГРС ITC У краины GCF. WcfW eatherCollector являются:

- получение от УкрГМЦ официальных прогнозов погоды по всей территории Украины и привязка их к зоне нагрузки каждой ГРС ГТС Украины;

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

WCF сервис GCF.WcfWeatherCollector импортир\ ст данные из базы данных УкрГМЦ по установленному регламенту в таблицу WeatherCenter базы данных (БД) ИАС ППППГ. В случае успеха, выполняется процедура обновления данных среднесуточной температуры воздуха и краткосрочного прогноза погоды по всем ГРС, которые хранятся в таблице Daily Weather. В случае отсутствия или повреждения части импортируемых данных выполняется процедура восстановления недостающих данных.

РИ, 2011, №3

93

Рис. 3. Архитектура

Интеллектуальный агент GCF.Tasks для запуска процедуры импорта использует класс TasklmportDaily Weather, реализующий интерфейс IJob библиотеки Quartz.NET. В случае невозможности получения данных запускается задача, которая выполняет процедуру импорта с интервал ом в 15 минут, пока процедура импорта не закончится с успехом или пока не исчерпается временной лимит, отведенный на задачу импорта данных. Если временной регламент исчерпан, а новые данные не поступили в систему, то

НАС ППППГ

запускается процедура получения прогнозов погоды по их ретроспективной информации (рис. 4).

GCF.wcfGdsDataCollector запускается в автоматическом режиме с помощью службы Windows и управляется интеллектуальным агентом GCF.Tasks. GCF.wcfGdsDataCollector импортирует данные из текстовых файлов, поступающих в объединенное диспетчерское управление по каждому управлению магистральным газопроводом (УМГ) по установленному регламенту. Интеллектуальный агент GCF.Tasks

ОДУ

©

м

Бердичівське ЛВУМГ л Сивира Си вира Антонів Шамраівка Володаркв

ГРС Житомир АГНКЄ 1 Кор ОСТИ LLIIB Ступениця

ГРС Житом.(місто) Жил

Гуйва

С“ '-ТИ

' Озерянка Глубочиця

Висока Піч Бабічовка

Бердичів Котельня ВРТТ

Гришнівці

Райгородок

Маркуил

Рея (Галвчин)

Сапки . Чуднів

Б.Ка ровинці ївачопіль-побут

ФІГИНЦІ

Н.Чврторія

Липно

ф График температуры

Рис. 4. Оконная форма результатов прогнозирования с вероятностными пределами

П/

Воублівка

94

РИ, 2011, №3

для запуска процедуры импорта данных использует класс TasklmportDailyGdsData, который реализует интерфейс IJob библиотеки Quartz.NET и представляет задачу, запускаемую с суточной периодичностью. В случае невозможности получения данных запускается задача, которая реализует стратегию импорта данных, аналогичную стратегии, которая применяется для TasklmportDailyWeather.

Учёт фактического потребления природного газа и расчёт его суточного прогноза ведется по всем выходам каждой ГРС. Г ТС Украины имеет чёткую многоуровневую иерархическую структуру. Эта структура однозначно представлена в реляционной модели базы данных. На верхнем уровне находятся УМГ, которые связаны с линейно-производственными управлениями (ЛПУ) отношением «один ко многим». ЛПУ связаны с ГРС также отношением «один ко многим». В свою очередь ГРС связаны со своими выходами также отношением «один ко многим». Под объектом ГТС будем понимать множество технологических элементов, которые объединены по производственному, административному или территориальному признаку. Для того чтобы определить суточный прогноз потребления природного газа за определенные сутки, нужно выполнить операцию агрегирования вычисленных прогнозов по множеству объектов, подчинённых данному объекту7. Таким образом, прогноз потребления природного газа і-м объектом ГТС вычисляется по формуле:

Q°bi(i)=EQ,obl(j)(i)

j=i ’ '

где (1) - прогноз потребления природного газа і-м объектом ГТС Украины в момент времени t с упреждением 1 (1 = 1, 2 , 3, ...); ^ '(1) - прогноз

потребления природного газа j -м объектом і-й группы объектов ГТС Украины в момент времени t с упреждением 1.

Для реализации процедуры расчёта суточного прогноза потребления природного газа был разработан специализированный WCF сервис GCF.WcfForecast, который запускается в автоматическом режиме с помощью службы Windows и управляется интеллектуальным агентом GCF.Tasks. Интеллектуальный агент GCF.Tasks для запуска процедуры расчета использует класс TaskDailyForecasting, который реализует интерфейс IJob библиотеки Quartz.NET. Для прогнозирования процессов потребления природного газа по каждому выходу ГРС ГТС Украины используется методика построения моделей, адекватных исследуемым процессам потребления природного газа [2, 6]. Особенностью процедуры построения моделей процессов потребления природного газа является её строгая формализация и возможность полной автоматизации всех этапов построения и использования моделей.

На начало технологических суток (09:00 часов) для службы режимов ОДУ ДК «Укртрансгаз» в НАС П П1111Г подготавливаются информационные отчёты о фактических и прогнозных объемах газопотребления, о параметрах метеорологических условий:

- в целом по ДК «Укртрансгаз»;

- по производственным подразделениям ДК «Укртрансгаз»;

- по регионам и областям Украины;

- по различным категориям внутренних потребителей ГТС Украины.

При разработке НАС ППППГ применялись следующие технологии и средства разработки:

- СУБД - MS SQL Server 2008 R2;

- .NET Framework 4.0;

- компоненты сервера приложений на базе программного фреймворка Windows Communication Foundation;

- клиентские приложения - Silverlight 4.0;

- язык программирования - С#.

2. Математические модели прогнозирования процессов потребления природного газа

Математическая модель процессов потребления природного газа, зависящих от метеорологических и хронологических факторов в операторной форме, может быть представлена в виде модели сезонной авторегрессии - проинтегрированного скользящего среднего с экзогенными переменными (САРПССЭ) [2, 6]:

N со1. (в) . 0**(в)

О _ w сі і?і і Ч_____

Xt іГі5‘(в) t_bi Ф++(В)“1> (2)

1 p

где Qt, t = l,n - исходный или преобразованный (нормированный или прологарифмированный) центрированный временной ряд потребления природного газа; п - объём выборки; Ft‘, t = l,n, i = l,N -исходные или преобразованные (нормированные или прологарифмированные) центрированные временные ряды изменения значений внешних (экзогенных) факторов; N - количество экзогенных факторов; В -оператор сдвига по времени на одну единицу назад, такой что В‘Хк = Х|. ; Ь; - величина задержки і-го

экзогенного временного ряда по времени относительно прогнозируемого временного ряда Qt; 5j. (В) и и (В) -полиномы передаточных функций отВ степеней fj и Cj соответственно;

(b)=i-511b-512b2 -...-5;.вг> ;

со' (в)= COq -со}В-со2В2 —... — gL ВС‘ ;

Ф++ (в) - обобщенный оператор авторегрессии пор ns * ns

рядка р+ = р* + Z D;S;, Р = Z PiSi ;

i=l i=l

РИ, 2011, № З

95

Ф+ (в)=<i>\(b)v‘^ v‘)2... v^"s.

p+V ’ р ' ' S>2 S>ns '

D;, і = 1, ns -порядок взятия разности Sj; Sj, і = 1, ns - период І-Й периодической компоненты, причем S ] = 1; ns - количество периодических компонент; Vs. и В 1 - упрощающие операторы такие, что

vSlxt=(i-Bs‘)-xt=xt-xt_Si.

Нередко такое интегрирование Vs. Xt =Xt -Xt_s. заменяется вычитанием из временного ряда авторегрессии первого или второго порядка, например,

?SiXt = Х t — (aSj + bSj Xt-Sj ) без предварительного вычитания из исходного временного ряда математического ожидания, а коэффициенты и Ь s.-

определяются по методу наименьших квадратов. Также достаточно часто для этих целей применяется модель экспоненциального сглаживания

VS.X, =aXt+(l-a)VSiXt_1,

а также модель двойного экспоненциального сглажи-

вания; Ф t (В) - обобщенный оператор авторегрес-

пч . „

* * s » С

сии порядка р вида ф (В)=ПФ (В 1 ) ; і Sj _____________ р* Г1 Р‘

Фр. (В 1 ), i = l,ns -ПОЛИНОМЫ ОТ BS| степеней Pj соответственно, определяющие составляющие авторегрессии периодических компонент с периодами Sj

соответственно; 0** (в)-обобщенныйоператорсколь-

* ns

зящего среднего порядка q = Zqisi вида

і—1

е**(В)= ГІ0І. (BSi); Oq.(BSi), і = І^Г -по-q і=0 q. q‘

g.

линомы от В 1 степеней q, соответственно, определяющие составляющие скользящего среднего периодических компонент с периодами Sj соответственно; а t - остаточные ошибки модели.

Решение задач стру ктурной идентификации параметров с j. г,. b j. і = 1. N осуществляется автоматически в результате анализа функции отклика на единичный импульс. Автоматизированная идентификация параметров р j, і = 1, ns - путём анализа частной автокорреляционной функции (ЧАКФ), а параметров q;, j = 1. ns - путём анализа автокорреляционной функции (АКФ). Параметры ns, Dj ,Sj также идентифицируются путём анализа АКФ и ЧАКФ процесса.

В15] представлены другие аддитивные, мультипликативные и смешанные модели прогнозирования целевых продуктов в инженерных сетях. Там же предложены различные модели для моделирования и прогнозирования трендовой, сезонной, недельной, оста-

точной составляющих процессов. Среди предлагаемых моделей прогнозирования отдельных составляющих процесса используются экспоненциальное сглаживание, модель Хольта-Винтерса, двойного экспоненциального сглаживания, фильтр Калмана, ряды Фурье, модели спектрального анализа, экспоненциально взвешенного скользящего среднего и пр. Идеи методов построения приведенных моделей развиваются также в статье [7].

Класс линейных дискретных передаточных функций с рациональной структурой вида:

со1 (В) ___

Q1 =-^—64'. i = LN (3)

s; (В)

используется для учета связи процессов потребления природного газа с метеорологическими факторами, где Qj, і = 1, N - составляющая процесса Q t, связанная с і-м метеорологическим или организационным фактором Ft‘ - i = l,N-Перепишем (3) в виде:

C0q сО]В-аьВ 1-о!В-5‘2В2

-...-со' вс'

---------^-----BblFt‘

...-5' Bri n

і = 1, N •

В данной работе используется только один экзогенный фактор - изменение температуры окружающего воздуха в зоне нагрузки каждой ГРС, поэтому (3) можно записать следующим образом:

2

со0 — C0jB — со2В 1-5]В-52В2 -

— ИдВ1-

...-5ГВГ

ВЬТ,

(4)

где Tt, t = 1. п -временной ряд изменения температуры окужающего воздуха. Такая передаточная функция описывает чувствительность процессов потребления природного газа к колебаниям температуры. Степень влияния температуры воздуха на прогнозы газопотребления характеризуется так называемым коэффициентом усиления К, который определяется путем замены оператора в и переменной Tt в модели передаточной функции (4) на 1, т.е.

£ _ СОр - СО! -С02 — ... — С0с

1-5] -52 -...-5Г

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

96

РИ. 2011. №3

мечено в зимний период с момента начала и до окончания отопительного сезона. В остальное время года изменение температуры воздуха не оказывает существенного влияния на потребление природного газа. Расчёт степени влияния температуры воздуха на потребление природного газа осуществляется анализом коэффициента детерминации модели (4) передаточной функции.

Таким образом, вне отопительного сезона для прогнозирования процессов потребления природного газа используется модель авторегрессии - проинтегрированного скользящего среднего:

е**(в)

Р, = Фtwa'- (5)

р

Прогноз Qt(l) процесса потребления природного газа Qt+i с упреждением 1 в момент времени t с минимальной ере днеквадратической ошибкой определяется как условное математическое ожидание Qt+i в момент времени!. Формулы для вычисления прогнозов Qt (l) выводятся в [2, 6].

Для повышения надёжности и живучести системы данная модель используется как для восстановления фактических объёмов потребления природного газа, таки для восстановления прогнозов изменения температуры воздуха в случае неполучения их непосредственно от У крГМЦ с заменой Qt на Tt.

В ИАС ППППГ реализована усовершенствованная процедура вычисления доверительных границ прогнозов. В [2, 6] представлен вывод формул для вычисления доверительных границ прогнозов

Q,+1(t)=Q,(i)±u£Vv®,

2

где 1 - период упреждения прогноза; Qt (l) - прогноз процесса потребления природного газа с упреждением 1; V(l) - дисперсия прогноза процесса потребления природного газа с упреждением 1; Qt+1 (t) - значения

вероятностных пределов в момент времени t+І; и Е -

1 Є 2 квантиль уровня І- — стандартного нормального рас-

иЁ

пределения. Численные значения — выбираютсярав-

ными0.68,1.65,1.96 или 2.58 в зависимости от того, находятся ли будущие значения Qt+i в интервале

Qt(0-u е л/у(0 с вероятностью 0.5,0.9,0.95 или 0,99 2

соответственно.

При вычислении прогнозов в соответствии с уравнением. выведенным из выражения (5), дисперсия прогнозов с упреждением 1 равна

РИ, 2011, №3

v(i)=^zVf,

где оа - среднее квадратичное отклонение процесса остаточных ошибок at модели (5), а формула для вычисления весов i|/j. і = 0.1 — 1 приведена в [2,3, 6].

При вычислении прогнозов по уравнению, выведенному из выражения (2), дисперсия прогнозов с упреждением 1 равна

N 1-і

v(i) = ZZ^<

i=lj=0 1

„ 1-1 „ 2 2 ■°a ХУ;

1—0

здесь oa - среднее квадратичное отклонение процесса остаточных ошибок at модели (2); o(tj . і = 1. N -средние квадратичные отклонения временного ряда остаточных ошибок моделей процессов изменения экзогенных факторов, а формулы для вычисления

весов V|/j,і =0,1-1 И Xjj, і = 1,N , j = 0,1-1 приведены в [2, 6].

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

.. « «22 учета изменения каждой из дисперсии , оц. ,

і = 1, N от времени и имеет вид:

о2 = \у + 0(в)а2 + ^>(в)о2 , (6)

где о2 - временной ряд изменения дисперсии остаточных ошибок at моделей (2) или (5); е(в) = 01В - 02В2 -... - 0рвт, ^в)= ^в - <?>2В2 -... - -

полиномы ОТ В модели (6), скользящего среднего и авторегрессии соответственно; at - остаточные члены модели. Модель GARCH(m.k) может быть записана через модель APCC(s, ш) следующим образом:

,2_ w + (l-p(B)) _

‘ = М^ЫВ))'1’

где s = max(r, m), V, = а,2 - с,2 •

Для обнаружения разладок процессов в системе используются метод отношения правдоподобия и метод кумулятивных сумм [8].

В ИАС ППППГ предусмотрены процедуры коррекции фактических данных соответствующим нерегулярным дням (государственным праздникам, выходным дням и т.п.) и соответствующие корректирующие коэффициенты, которые записываются в специальную таблицу и периодически адаптируются совместно с коэффициентами модели. Прогнозные значения, вычисленные по скорректированным данным, умно-

97

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

Выводы

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

Литература: 1. Business Process Modeling Notation VI. 1 / / OMG Available Specification. OMG Document Number: formal/2008-01-17. URL: http://www.omg.org/spec/BPMN/ 1.1/PDF. 2. Евдокимов А. Г., Тевяшев А..Д. Оперативное управление потокораспределением в инженерных сетях. Х.:Вищашк., 1980.144с. З.БоксДж.,ДженкинсГ. Анализ временных рядов. Прогноз и управление: Пер. с англ. М.: Мир, 1974. Вып. I. 406 с. 4. Бокс Дж. Дженкинс Г. Анализ временных рядов. Прогноз и управление: Пер. с англ. М.: Мир, 1974. Вып. II. 197 с. 5. Бэнн Д. В., Фармер Е. Д. Сравнительные модели прогнозирования электрической нагрузки: Пер. с англ. М.: Энергоатомиздат, 1987. 200 с.

6. Тевяшев А. Д, ЕЦелкалин В. Е1. Сравнительный анализ методов прогнозирования процессов потребления природного газа // Восточно-европейский журнал передо-

вых технологий. 2009. №42. С. 8 - 16. 7. ЕЦелкалин В. Н., Тевяшев А. Д. Трендовый и декомпозиционный подходы прогнозирования процессов потребления электроэнергии // Восточно-европейский журнал передовых технологий. 2011. № 5/4 (53). С. 30 - 37. 8. Жиглявский А. А., Красковский А. Е. Обнаружение разладки случайных процессов в задачах радиотехники. Л.: Изд-во Ленингр. ун-та. 1988. 224 с.

Поступила в редколлегию 28.07.2011

Тевяшев Андрей Дмитриевич, д-ртехн. наук, проф., зав. каф. прикладной математики ХНУРЭ. Научные интересы: математическое моделирование энергетических систем, теория стохастических моделей. Адрес: Украина, 61166,Харьков,пр. Ленина, 14,тел.: (057) 702-14-36,e-mail: tevj a shev@kture. kharkov. ua.

Выходцев Евгений Иванович, научный сотрудник кафедры прикладной математики ХНУРЭ. Научные интересы: информационно-аналитические системы, системы поддержки принятия решений, компьютерное моделирование. Адрес: Украина, 61166, Харьков, пр. Ленина, 14,тел.: (057)702-14-36,

e-mail: eugene.vykhoddtsev@hotmail. com.

Щелкалин Виталий Николаевич, инженер 1 -й категории кафедры прикладной математики ХНУРЭ. Научные интересы: математическое моделирование, прогнозирование, теория автоматического управления, искусственные нейронные сети. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел.: (057) 719-91-09, e-mail:

vitalii. shchelkalin@gmail.com.

Игнатова Юлия Васильевна, заместитель начальника отдела разработки и сопровождения информационных систем диспетчерского управления филиала объединенного диспетчерского управления ДК «Укртрансгаз». Научные интересы: автоматизированные системы управления, численные методы. Адрес: Украина, 01021, Киев, пер. Кловский-9/1, тел.: (044) 461-22-94, e-mail: j ulia @utg. com. ua.

98

РИ, 2011, №3