Уфа : УГАТУ, 2010_ОУ ' <и_Т. 14, №5(40). С. 110-116
УПРАВЛЕНИЕ, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И ИНФОРМАТИКА
УДК 622.3:004.023
А. В. Жевак, В. Ю. Арьков
МОДЕЛИРОВАНИЕ И ОПТИМИЗАЦИЯ ПРОЦЕССА СБОРА ДАННЫХ В СИСТЕМЕ МОНИТОРИНГА ГАЗОВЫХ СКВАЖИН И ГАЗОСБОРНЫХ КОЛЛЕКТОРОВ
От эффективности сбора данных в системе мониторинга, основанной на автономных сенсорных узлах, собирающих измеряемые на контрольных точках газосборных коллекторов геофизические параметры в контрольный пункт, зависит среднее энергопотребление узлов, а следовательно, и затраты на замену источников питания в сложных климатических условиях Крайнего Севера. В статье предложена новая модель процесса сбора данных в системе мониторинга, в отличие от известных, учитывающая возможности современных радиомодулей изменять мощность передатчика (что существенно влияет на энергопотребление). Рассматривается цифровой алгоритм оптимизации, учитывающий специфику предложенной модели. Мониторинг; оптимизация; моделирование; транспортная задача; расписание
Отсутствие оперативной информации об устьевых параметрах скважин и состоянии газосборных коллекторов делает процесс поддержания объемов добычи газа весьма проблематичным. Это связано, прежде всего, с несвоевременным обнаружением процесса развития жидкостных пробок, снижающих пропускную способность шлейфов и провоцирующих остановку скважин, что может длительное время оставаться незамеченным [6].
Наиболее очевидным решением являлась бы электрификация и оборудование фонда скважин кабельными системами телеметрии. Однако, в связи с удаленностью объектов добычи газа, затраты на такую автоматизацию будут очень большими и вряд ли скоро окупятся, особенно на завершающем этапе эксплуатации месторождений.
В связи с этим, для решения задач телеметрии устьевых параметров скважин и контроля за состоянием газосборных коллекторов предлагается использовать автономную информационно-измерительную систему (по перепаду давления в разных точках коллектора возможно выявить жидкостную пробку, как показано на рис. 1) [1].
На рис. 2 схематически представлена система мониторинга скважин и газосборных коллекторов, развернутая на УКПГ-5 ООО «Газпром добыча Уренгой».
Ограниченная емкость источников питания сенсорных узлов системы ведет к необходимости оптимизации процесса сбора данных по критерию энергопотребления. Сложность натурного эксперимента делает оправданным разработку комплекса предметно-ориентированных программ для проведения вычислительного эксперимента.
Для решения задачи оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов целесообразным является применение математического моделирования, численных методов и комплексов программ. Математическая модель системы мониторинга (технического объекта) построена с применением методов теории графов. Алгоритм оптимизации сбора данных реализован в виде комплекса проблемно-ориентированных программ для проведения вычислительного эксперимента. Таким образом, проведено комплексное исследование научно-технической проблемы мониторинга газовых скважин и газосборных коллекторов с применением современной технологии математического моделирования и вычислительного эксперимента.
1. СОСТОЯНИЕ ВОПРОСА
Известные модели процесса сбора данных в автономных системах мониторинга можно разбить на 2 основных класса: модели сбора данных по расписанию и адаптивные модели. Важно отметить, что нельзя рассматривать алгоритм оптимизации процесса сбора данных отдельно от модели процесса сбора данных и наоборот. Основной идеей алгоритмов сбора данных по расписанию является создание некоторого расписания передачи данных для узлов системы. Расписание жестко прописывает для каждого узла, что он должен делать в каждый момент времени. При этом время может быть разбито на интервалы. Как правило, расписание строится таким образом, чтобы исключить дублирование данных [2]. Адаптивные модели реализуют идею так называемых самоорганизующихся узлов.
Контактная информация: (347) 272-89-81
Место скопления конденсата
\\
ЛР = |Р1-Р2|
Рис. 1. Определения жидкостной пробки по перепаду давления на контрольных точках
Рис. 2. Система мониторинга скважин и газосборных коллекторов, развернутая на УКПГ-5 ООО
«Газпром добыча Уренгой»
Основной смысл заключается в том, что в системе не существует никакого заданного извне расписания. Узлы системы самостоятельно определяют кому передавать данные в тот или иной момент времени. В таблице приведена сводка основных характеристик существующих моделей информационных потоков.
Характеристики моделей процесса сбора данных в автономных системах мониторинга
Класс Время Энер- Под- Отказо-
моделей сбора гопот- держка устойчи-
данных ребле- мобиль- вость
ние ности
Основан- малое, низкое, нет Низкая
ные фиксиро- посто-
на распи- ванное янное
сании
Адаптив- большое, высо- да Высокая
ные без верх- кое,
ней гра- непред-
ницы пред-
ска-
зуемое
Ни одна из известных моделей процесса сбора данных в автономных информационно-измерительных системах не отвечает в полной
мере требованиям, предъявляемым к системе мониторинга газовых скважин и газосборных коллекторов, а именно: отсутствует формула расчета энергопотребления для каждого узла системы, не учитывается возможность переключения передатчиков узлов между несколькими уровнями мощности, что могло бы снизить энергопотребление.
Таким образом, актуальной и практически востребованной является задача моделирования и оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов.
2. ПОСТАНОВКА ЗАДАЧИ
Целью работы является повышение эффективности сбора данных в системе мониторинга газовых скважин и газосборных коллекторов путем создания эффективных моделей и алгоритмов оптимизации, а также оценка эффективности предложенных решений. Для достижения указанной цели поставлены и решены следующие задачи:
1) разработать математическую модель процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов;
2) разработать алгоритм оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов;
3) разработать алгоритм генерации тестовых графов, описывающих процесс сбора данных, для проведения вычислительного эксперимента;
4) пеализовать разработанные алгоритмы в виде комплекса программ для проведения вычислительного эксперимента.
3. МЕТОДИКА ИССЛЕДОВАНИЯ
При решении поставленных задач использовались методы теории графов, теории вероятностей и математической статистики, объектно-ориентированного программирования, проведения вычислительного эксперимента.
Система мониторинга, состоящая из п узлов (п-1 сенсорного и 1 центрального), есть связный неориентированный граф G(V, Е), где
V = (уь у2, ..., уп} - множество вершин графа G, причем: у1, у2, ..., vn_l - сенсорные узлы; уп - центральный узел; Е = {е1, е2, ..., ет} -множество ребер графа G, причем: {и, V} е Е, если между и и V возможна прямая передача данных; *(и, V) - минимальный уровень мощности передатчика, при котором возможна прямая передача данных между и и V.
Представление системы мониторинга в виде графа встречалось в работах других авторов, однако в модель был также введен новый элемент - функция *(и, V), указывающая минимальный уровень мощности передатчика, необходимый для установления прямой связи между узлами и и V. Такое нововведение обусловлено наличием у большинства современных радиомодулей возможности переключения между несколькими уровнями мощности при работе в режиме передачи. При этом чем выше мощность, тем выше потребляемый ток и радиус действия. Так как не всегда для связи между 2 узлами требуется использовать самый высокий уровень мощности, то за счет снижения мощности может быть достигнута экономия энергии, а как следствие - увеличения срока службы узлов без замены источника питания. Так, используемые в системе мониторинга газовых скважин и газосборных коллекторов радиомодули ХЕ1203 фирмы Хешю8 имеют 4 дискретных уровня мощности передатчика [3].
Каждый сенсорный узел генерирует пакет данных, включающий в себя значения измеряемых им параметров. Все такие пакеты с каждого из сенсорных узлов должны быть переданы в центральный узел. Сбор данных в системе мониторинга осуществляется по расписанию. Время разбивается на единичные временные
интервалы. Период, за который данные со всех сенсорных узлов поступают в центральный узел, называется циклом сбора данных. Цикл сбора данных состоит из Ь единичных временных интервалов. Каждый узел в течение 1 единичного временного интервала выполняет 1 из 3 возможных действий: 1) принимает данные от другого узла, 2) передает данные другому узлу, 3) не осуществляет прием/передачу данных.
В последнем случае узел расходует гораздо меньше энергии за счет временного отключения приемопередатчика, пребывания микропроцессора в спящем режиме и тому подобных мер.
Расписание сбора данных описывается математически как упорядоченное множество {51, 52, ..., 5Ь }, где 57 описывает передачу данных в 7-й единичный временной интервал и представляет собой множество упорядоченных пар вершин: 57 = {(аь Ъ{), (а2, 62), ..., (ат, Ьт)}, где а7 - узел, передающий информацию в течение 7-го единичного временного интервала узлу Ь7, а Ь7, соответственно, узел, принимающий информацию от узла а7, в течение 7-го временного интервала.
Расход энергии вычисляется исходя из расписания сбора данных и технических характеристик используемого оборудования (токов, потребляемых аппаратными компонентами сенсорного узла в различных режимах). Формула расчета расхода энергии сенсорным узлом системы, приведенная ниже, используется для прогнозирования срока службы каждого сенсорного узла системы без замены источника питания. Предложенная в работе модель является первой в своем роде моделью процесса сбора данных, позволяющей строить такой прогноз. Формула расчета расхода энергии сенсорным узлом 7 за 1 цикл передачи данных приведена ниже:
е7 = Ье + *т (ет + е* - ел )+ *с (е* - е) +
+ (г + 1)**(ек -е,)+Г*р(е* + ег -es)+
/л N
Г (г V + Vc)+ кг¥т + Vc + V, + Vp)
+
Х(ег + е* - е6, ) +
2 Г - фг +
В
+
х
+
г V + V, + Vp)
В
+1 2 г -1
+
Щ + Гг - е, )) +
л
ГV +(кг+ЛУт + Vc ) В
х
+ 1 2 г-1
х
Х(*г + е* - е, )■
Г
Здесь
tnl - время, необходимое АЦП для измерения всех необходимых параметров;
tc - время, необходимое микропроцессору для обработки значений, полученных от АЦП;
tw - время перехода микропроцессора в рабочий режим (из спящего режима);
tp - максимальное время ожидания сигнала приемником (интервал времени между включением узла-приемника на прием и началом передачи данных узлом-передатчиком);
tr - время перехода радио-модуля в режим приемника;
tt - время перехода радио-модуля в режим передатчика;
es - ток, потребляемый микропроцессором в спящем режиме;
em - ток, потребляемый АЦП во время измерения;
ew - ток, потребляемый микропроцессором в рабочем режиме;
er - ток, потребляемый приемопередатчиком в режиме приемника;
et - ток, потребляемый приемопередатчиком в режиме передатчика;
B - скорость обмена данными между двумя узлами;
Vp - объем служебной информации, передаваемой при каждой передаче данных (кроме передачи отклика);
Vc - размер пакета отклика; Vt - размер 1 пакета синхронизации времени;
Vm - размер 1 пакета измерений 1; Lr - реальная длительность цикла сбора данных;
г - ожидаемое среднее количество попыток связи;
г/ - количество узлов, у которых узел / принимает данные в течение 1 цикла передачи данных;
^ - количество пакетов измерений, принимаемых узлом / в течение 1 цикла передачи данных;
Ь - сумма значений энергопотребления радиопередатчика /-го узла при передаче откликов во все узлы, передающие данные напрямую узлу / (с учетом необходимых уровней мощности);
^ - энергопотребление радиопередатчика /-го узла при передаче пакетов измерений (с учетом требуемого уровня мощности).
На основании формулы (1) определена формула расчета общего расхода энергии (т. е. суммарного расхода энергии всеми сенсорными узлами системы за 1 цикл сбора данных).
E = (n - 1)Lres + (n -А» (em + - es) + ,
+ (n - l)(c (e - es) + (2n - q - 1>w (ew - es) +
+ (n - q)tp {er + ew - es) + + (n - l)f 2 r-1\ (er + ew - es) +
+ (n -1)2rtr(er + eK -e^) +
A n-1
+ 2rtt£(z, + ew -es) +
+ 12 r -1 \t
£ z + (n - qtew- es)
V itV/Q
A V + V + V
+ (n -1)r--p-- (er + eu
B
+
r
+ — B
- es) +
n-1 Л
(n - q)(Vp + Vm + Vc)+ Vm £k
t=1 J
(e
+r
A n- Vp + Vc +(k, + 1)Vm r + ew - es )+ r £ B ч
, =1 X(Z, + ew - es)
Vt + V + V ( ^
' " ■ £zt +(n-qXew-es)
[teV/Q
(2)
B
Общий расход энергии E является критерием оптимизации в задаче оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов, постановка которой приведена ниже.
Дано:
G(V, E) - граф, описывающий систему мониторинга;
B V V V V, t t t t t t, e e e e C - характеристики системы мониторинга.
Найти:
S - расписание сбора данных для заданной системы мониторинга, такое, что E ® min (общий расход энергии минимален).
Здесь C = (ei, c2, ..., cm} - множество, описывающее уровни мощности передатчика; c, -ток, потребляемый РМ в режиме передачи при ,-м уровне мощности; m - количество уровней мощности.
Задача оптимизации процесса сбора данных, посредством определенных манипуляций, сводится к задаче поиска максимального потока минимальной стоимости в транспортной сети. Таким образом, разработанный алгоритм состоит из 3 основных шагов:
1) построение транспортной сети по заданному графу системы мониторинга;
2) нахождение максимального потока минимальной стоимости в транспортной сети, построенной на предыдущем шаге;
3) построение расписания передачи данных для исходной системы мониторинга по потоку в транспортной сети, найденному на предыдущем шаге.
,=1
х
X
+
Рис. 3. Граф системы мониторинга (а) и соответствующая ему транспортная сеть (б)
Второй шаг решается с помощью модифицированного алгоритма Форда - Фалкерсона [5]. Шаги 1 и 3, в свою очередь, представляют собой новые задачи, для которых не существует известных решений.
На первом шаге алгоритма по исходному графу системы мониторинга G(V, E) строится транспортная сеть G'(V, E') таким образом, что V = V и 5 и t , где 5 - исток, t - сток; E' = E и и (5, 1) и (5, 2) и ... и (5, n - 1) и (n, t) - множество дуг; p(5, v) = 1, для всех v, таких что (5, v) е E'; p(u, v) = да, для всех u, v, таких что (u, v) е E' и и Ф 5; a(u, v) = 0, для всех и, v, таких что (u, v) е E' / E; a(u, v) = с^, v) + er + 2ew - 2es, для всех u, v, таких что (u, v) е E, u Ф n, v Ф n; a(u, n) = cw(u, v) + ew - e5, для всех u, таких что (u, n) е E. Здесь p(u, v) - пропускная способность дуги (u, v); a(u, v) - цена прохождения потока по дуге (u, v).
На рис. 3 показан пример построения транспортной сети по графу системы мониторинга согласно приведенному выше алгоритму.
Третий шаг алгоритма включает в себя операцию «объединения потоков», корректность которой доказывается в теореме, формулировка которой приведена ниже.
Теорема об объединении потоков: если в транспортной сети, построенной по графу системы мониторига, для некоторого узла u существует несколько узлов v1, v2, ..., vk, таких что из u в каждый из них проистекают потоки ненулевой величины, то объединение этих потоков в единый поток из u в v, (для любого i е [1, k]) не изменит ни суммарной величины потока в графе, ни его стоимости.
Сложность проведения натурного эксперимента приводит к необходимости создания комплекса компьютерных программ для имитационного моделирования [4]. Для проведения последнего спроектирован (с помощью унифицированного языка моделирования UML (Unified
Modeling Language)) и реализован (на языке программирования C++) комплекс программ для ЭВМ. Данный комплекс позволяет тестировать алгоритм оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов всевозможных конфигураций не прибегая к дорогостоящему и длительному натурному эксперименту.
Для проведения вычислительного эксперимента было необходимо предложить алгоритм генерации тестовых графов систем мониторинга (входных данных алгоритма оптимизации). Постановка задачи генерации тестовых графов систем мониторинга приведена ниже.
Дано: n - количество узлов системы; D = = {db d2, ..., dm}, где dt - количество связей между узлами с минимальным уровнем мощности передатчика равным i. Требуется сгенерировать граф G( V, E), такой что:
G(V, E) - связный неориентированный граф;
|V = n; w(u, v) е {1, 2, ..., m} V {u, v} е E; количество {u, v} е E и w(u, v) = i, равно di для всех i е {1, 2, ..., m}. Алгоритм генерации тестовых графов систем мониторинга, решающий поставленную задачу:
Задать V = 0, E = 0.
Задать V' = {v1, v2, ., vn} - множество всех допустимых вершин.
Задать E' = \(u,v) - множество всех до-
uеE',vеE'
пустимых ребер.
Построить случайное остовое дерево:
Задать u = случайная вершина из V.
Задать V = V и u.
Задать u = случайная вершина из V.
Задать v = случайная вершина из V.
Задать V = V и u, V' = V' / u.
Задать E = E и (u, v), E' = E' / (u, v).
Задать w(u, v) = случайный элемент из {1, 2,., m}.
dw(u,v) = dw(u,v) - 1 .
б
а
Если IV < п, перейти на шаг 4с.
Дополнить граф случайным образом до достижения заданного числа ребер:
Пока Х^г- > 0, выполнять:
Задать (и, V) = случайное ребро из Е'.
Задать Е = Е и (и, V), Е' = Е' / (и, V).
Задать w(u,v) = случайный элемент из {1, 2,..., т}.
dw(u,v) = dw(u,v) - 1.
Примечание к шагам 4g и 5с: вероятность генерации элемента / равна
т ]=1
Новизна данного алгоритма заключается в его способности случайным образом генерировать графы систем мониторинга Е) с заданным числом вершин учетом и количеством ребер по каждому уровню мощности радиопередатчиков (т. е. с учетом новой функции w(u, V)).
4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
На рис. 4 и 5 показаны результаты экспериментов для графов различной плотности, состоящих из 100 узлов.
общий расход энергии, мАс 35
30
25
20
Е
15
10 20 30 40 50 60 70 80 90 100 Ь-плотность сети, %
Рис. 4. Общий расход энергии (Е) по результатам вычислительного эксперимента для графов систем мониторинга, состоящих из 100 узлов
На рис. 5 по оси ординат отложен средний срок службы сенсорного узла системы до замены источника питания (К), вычисляемый по
следующей формуле: К =-
1
-нЕ -
24 • 365,25 4 средний срок службы сенсорного узла системы в годах. На рис. 4 и 5 сплошной линией показан график результатов для алгоритма оптимизации с использованием возможности варьирования уровня мощности радиопередатчика, а пунктирной - без использования такой возможности. Из результатов вычислительного эксперимента видно, что варьирование уровней мощности
радиопередатчиков позволяет снизить общий расход энергии (значение целевой функции Е) в среднем на 14-16%.
ю с
3,0 -|
2,5 -
£Р
О
01
2,0 -
1,5
10 20 30 40 50 60 70 80 90 100 d - плотность сети, % Рис. 5. Средний срок службы сенсорного узла системы мониторинга (К) по результатам вычислительного эксперимента для графов систем мониторинга, состоящих из 100 узлов
Таким образом, результаты эксперимента доказывают целесообразность добавления нового параметра в модель процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов.
5. ПРИЛОЖЕНИЕ РЕЗУЛЬТАТОВ
Разработанный и реализованный комплекс программ для проведения вычислительного эксперимента позволяет протестировать алгоритм оптимизации сбора данных менее чем за 1 минуту, не прибегая к натурному эксперименту длительностью в 2-3 года.
Разработанный алгоритм оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов позволяет составлять расписание сбора данных для системы, обеспечивающее функционирование сенсорных узлов системы без замены источника питания в среднем в течение 2-3 лет, что превышает период изменения конфигурации системы мониторинга. Использование в процессе оптимизации такого нового фактора, как дискретные уровни мощности передатчиков, дает эффект в виде повышения среднего срока службы сенсорного узла системы без замены источника питания на 14-16%. Данная модель также позволяет осуществить прогноз срока службы каждого сенсорного узла без замены источника питания с точностью 4-6%. Такая погрешность является незначительной на фоне вариации емкости источника питания в зависимости от температурного режима (более 50%).
ВЫВОДЫ
1) Разработана математическая модель процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов, основанная на идее сбора данных по расписанию и, в отличие от существующих, учитывающая такой фактор, как наличие у передатчиков узлов системы нескольких дискретных уровней мощности. Данная модель может быть использована для расчета расхода энергии для каждого узла системы. Данная модель позволяет избежать конфликтов и дублирования при передаче данных, обеспечивает сбор данных со всех сенсорных узлов за фиксированное время (менее 2 минут для системы из 100 узлов) и не требует изменения расписания сбора данных до тех пор, пока не изменится конфигурация системы мониторинга. Прогноз срока службы каждого сенсорного узла без замены источника питания по математической модели отличается от экспериментального показателя в среднем на 4-6%, что является незначительным на фоне вариации емкости источника питания в зависимости от температурного режима (более 50%).
2) Предложен алгоритм оптимизации процесса сбора данных в системе мониторинга газовых скважин и газосборных коллекторов, основанный на разработанной модели, использующий модифицированный алгоритм Форда-Фалкерсона для поиска максимального потока минимальной стоимости в графе. Предложенный алгоритм содержит в себе такие новые шаги, как преобразование графа системы мониторинга в транспортную сеть и построение расписания сбора данных для системы мониторинга по потоку в упомянутой транспортной сети. Алгоритм позволяет составлять расписание сбора данных для системы, обеспечивающее функционирование сенсорных узлов системы без замены источника питания в среднем в течение 2-3 лет. Использование в процессе оптимизации такого нового фактора, как дискретные уровни мощности передатчиков, дает эффект в виде повышения среднего срока службы сенсорного узла системы без замены источника питания в среднем на 14-16%.
3) Предложен алгоритм генерации тестовых графов, описывающих процесс сбора данных, для проведения вычислительного эксперимента. Алгоритм включает в себя такие шаги, как построение случайного остового дерева и дополнение графа до достижения требуемого количества ребер по каждому уровню мощности передатчика. В отличие от ранее известных, данный алгоритм позволяет генерировать графы с заданным числом узлов (до 100) и количеством
связей по каждому уровню мощности передатчика.
4) Разработан комплекс программ для проведения вычислительного эксперимента, реализованный на языке программирования C++, реализующий новую модель процесса сбора данных и новые алгоритмы оптимизации и генерации тестовых графов систем мониторинга. Данный комплекс программ позволяет протестировать алгоритм оптимизации на 2000 различных графах систем мониторинга менее чем за 1 мин., не прибегая к натурному эксперименту длительностью в 2-3 года.
СПИСОК ЛИТЕРАТУРЫ
1. Жевак А. В., Арьков В. Ю., Фридлянд А. М.
Оптимизация сбора данных в беспроводных сенсорных сетях с использованием нейронной сети с градиентным алгоритмом обучения // Нейрокомпьютеры: разработка, применение. М., 2007. Вып. 10. С. 47-49.
2. Коловертнов Г. Ю., Ланчаков Г. А. Беспроводная информационно-измерительная система контроля технологических параметров скважин и газосборных коллекторов // Газовая промышленность. М., 2008. Вып. 8. С. 26-28.
3. Жевак А. В. Распределенная автономная информационно-измерительная система мониторинга состояния трубопровода // Сб. тезисов второй науч.-техн. конф. молодых спец. Уфа, 2006. С. 70-71.
4. Трусов П. В. Введение в математическое моделирование: Учеб. пособие. М.: Логос, 2004. 307 с.
5. Акимов О. Е. Дискретная математика: логика, группы, графы. М.: Лаборатория Базовых Знаний, 2001. 352 с.
6. Hauge E., Aamo O. M. Model-Based Monitoring and Leak Detection in Oil and Gas Pipelines // SPE Projects, Facilities & Construction. 2009. Vol. 4. P. 53-60.
ОБ АВТОРАХ
Жевак Алексей Владимирович, инженер. Дипл. инж. по прикладной информатике в экономике (УГАТУ, 2005). Иссл. в обл. моделирования и оптимизации сбора данных в системах мониторинга.
Арьков Валентин Юльевич,
проф. каф. автоматизир. систем упр-я. Дипл. инж. по пром. электронике (УГАТУ, 1986). Д-р техн. наук по системному анализу, управлению и обработке информации (УГАТУ, 2002). Иссл. в обл. идентификации и моделирования сложных систем.