Научная статья на тему 'Математическое моделирование процесса теплообмена в регенераторе газоперекачивающего агрегата с использованием аппарата обратных задач'

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

CC BY
141
33
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГАЗОПЕРЕКАЧИВАЮЩИЙ АГРЕГАТ / ТЕХНИЧЕСКОЕ СОСТОЯНИЕ / РЕГЕНЕРАТОР / ДИАГНОСТИРОВАНИЕ / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ОБРАТНАЯ ЗАДАЧА / ТЕМПЕРАТУРНОЕ ПОЛЕ / GAS PUMPING UNIT / TECHNICAL STATE / REGENERATOR / DIAGNOSING / MATHEMATICAL MODEL / INVERSE PROBLEM / TEMPERATURE FIELD

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Замиховский Л. М., Петрив С. Я.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The necessity of technical state control of regenerators during operation of gas pumping unit was proved in the article. Іt was proposed to develop a new method based on the use of methods of mathematical modeling of heat distribution on the surface of the regenerator and hardware methods to determine its temperature. It is considered the regularization algorithm of incorrect inverse problem of heat conduction in the material of regenerator design using values of temperature fields, which were defined experimentally

Текст научной работы на тему «Математическое моделирование процесса теплообмена в регенераторе газоперекачивающего агрегата с использованием аппарата обратных задач»

УДК 622.691.4.052 :(517.956+519.876.5) DOI: 10.15587/2313-8416.2015.40764

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ТЕПЛООБМЕНА В РЕГЕНЕРАТОРЕ ГАЗОПЕРЕКАЧИВАЮЩЕГО АГРЕГАТА С ИСПОЛЬЗОВАНИЕМ АППАРАТА ОБРАТНЫХ ЗАДАЧ

© Л. М. Замиховский, С. Я. Петрив

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

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

The necessity of technical state control of regenerators during operation of gas pumping unit was proved in the article. It was proposed to develop a new method based on the use of methods of mathematical modeling of heat distribution on the surface of the regenerator and hardware methods to determine its temperature. It is considered the regularization algorithm of incorrect inverse problem of heat conduction in the material of regenerator design using values of temperature fields, which were defined experimentally

Keywords: gas pumping unit, technical state, regenerator, diagnosing, mathematical model, inverse problem, temperature field

1. Введение

Вопросам оценки технического состояния газотранспортной системы Украины, учитывая ее длительную эксплуатацию (более 30 лет), всегда уделялось и уделяется значительное внимание. Сегодня этот вопрос является особенно актуальным, поскольку в период сложившейся экономической ситуации в стране от эффективной и надежной эксплуатации оборудования компрессорных станций (КС) - газоперекачивающих агрегатов (ГПА) зависит выполнение контрактов по транзиту газа в Европу и обеспечение газом собственных нужд. Одним из важных элементов ГПА являются подогреватели воздуха, который поступает в камеру сгорания - регенераторы. От целостности регенератора в значительной степени зависит мощность ГПА. Наличие утечек воздуха из регенератора из-за появления неплотностей приводит к снижению мощности ГПА до 30 %. В случае появления сквозных трещин в корпусе регенератора продукты сгорания снова попадают в камеру сгорания, в результате чего возрастает температура рабочего тела, а это, в свою очередь, приводит к увеличению концентрация оксида углерода и азота в воздухе вследствие чего ухудшается состояния окружающей среды в зоне КС [1].

2. Постановка проблемы

Сегодня из известных типов теплообменных аппаратов - рекуперативных, регенеративных (пластинчатых, трубчатых) и смешанного типа широкое применение на ГПА нашли теплообменные аппараты рекуперативного типа [2]. Среди них пластинчатый регенератор является одним из наиболее эффективных, поскольку степень регенерации достигает до 0,92. Однако они имеют ряд недостатков, которые обусловливают сегодня широкое внедрение на ГПА

трубчатых регенераторов. Основными из них является достаточно сложная технология изготовления, которая не позволяет при сборке пакетов обеспечить высокую герметичность значительного числа сварных швов, а также высокая стоимость исходного материала для изготовления регенераторов (нержавеющая сталь 1Х18Н9Т). Проведенный в [2] анализ причин и факторов, обусловливающих изменение технического состояния регенераторов показал, что основным дефектом являются трещины сварных швов, которые приводят к потере герметичности его воздушного тракта. Так, простая и надежная турбина ГТК-10-4 с подтвержденным ресурсом более 100 000 часов, в то же время имеет много проблем с регенератором пластинчатого типа, особенно при термоциклических нагрузках. Уже после нескольких десятков пусков и остановок ГПА регенератор начинает терять герметичность, что приводит к ухудшению параметров ГТУ в целом [3].

В связи с этим актуальной является задача контроля технического состояния регенератора, что требует использования методов и средств их диагностики.

3. Литературный обзор

Сегодня отсутствуют промышленные методы и средства диагностирования технического состояния регенераторов в процессе эксплуатации ГПА. В [2] отмечено, что техническое состояние регенератора косвенно можно контролировать путем измерения давления межтрубного пространства и измерения температуры среды на входе и на выходе межтрубного, и на входе трубного пространства. Для определения утечек воздуха из секций регенератора используются методы, которые, в основном, базируются на их опрессовке на неработающем агрегате.

В [4] предлагается в качестве диагностической термодинамической модели эффективности работы регенератора в эксплуатационных условиях использовать отношение параметров сходства Нуссельта (Иы) и Эйлера (Еы) от параметров Рейнольдса (Кг) и Прандтля (Рг), а также параметра Еы от параметров Кг и Рг (при соответствующих возможностях определения этих параметров при эксплуатации регенераторов). Проявление любых дефектов в процессе теплопередачи и гидравлических сопротивлений будет характеризоваться изменением этих зависимостей по параметру Кг (фактически относительно значения массовой скорости потока). Ухудшение состояния регенератора будет сопровождаться снижением его КПД при тех же режимных параметрах. В условиях эксплуатации для определения утечек воздуха из секций пластинчатых регенераторов используются следующие методы: по темпу падения давления в заглушенных секциях регенератора; по данным специально установленной мерной шайбы в режиме критического истечения воздуха через неплотности регенератора: за изменением параметров работы агрегатов, от которых осуществляется отбор воздуха на опрессовку. Кроме этого, в условиях эксплуатации на разных газопроводах делались попытки использовать и другие способы [4]:

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

- по приращению мощности ГТУ после заваривания щелей в секции регенератора.

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

4. Разработка математической диагностической модели состояния регенератора

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

Рассмотрим математическую диагностическую модель состояния характеристик материала регенератора на основе данных о распространении тепла в нем, которые определяются экспериментально.

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

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

Постановка задачи

Рассматривается модельная задача теплопроводности в двухмерной прямоугольной области

V = |(х,у),0< х<Ьх,;0< у<Ьу}. Основным уравнением является двухмерное уравнение теплопроводности:

дТ д ( , ч дТ \ д ( . ч дТ ^

ст Ы^'^дхГдуИ^су } (1)

Уравнения должно быть дополнено соответствующими граничными условиями:

AT (x, y, t)\av +B ( T — Tc )\sv = F (x, y, t),

(2)

A - имеет смысл коэффициента теплопроводности Я, В - коэффициент конвективного теплообмена, F(x, y, t) - закон изменения температуры. И начальными условиями:

T(x,y,t)U = T0 (x, y) . (3)

Согласно [6], при решении обратных задач, необходимо проверять не только выполнение уравнения (1), но и условий (2) и (3).

Пусть для исследуемой области известен закон распределения температуры T = T (x, y, t), полученный с использованием аппаратных средств. Возникает задача восстановления поля величины a(x, y, t) в

исследуемой области. При этом необходимо решить следующие задачи:

- построение регуляризующего алгоритма задачи восстановления поля a (x, y, t) по известным

распределениям температуры;

- построение регуляризующего алгоритма для задачи восстановления температурного поля по его значениям на определенном множестве точек.

Восстановление поля a (x, y, t). Известно, что в тех случаях, когда a = const,

Ä ср'

(4)

где А - теплопроводность материала, с - его удельная теплоемкость, р - плотность. Зная распределение а(х, у, t), можно фиксировать изменение указанных характеристик от начальных значений по известным характеристикам исходного материала и характеристикам материала с дефектом (вкрапления воздуха, газа и т. п.). Уравнения (1) записывается с использованием метода конечных разностей [7]:

rpn+1 rpn Tk ,i — Tk ,i

ar_ ax

— a,r

ar_ ax

ат_ ay

hx

arN ay' k ,i

a =

T

+

h

y

где fn = f (xk,yi,tn);r;hx,hy - шаг по координате соответственно t, x та y (рис. 1)

'ST л n dx

rj^n+1 _ rj^n

A,,+i«n++i + A+u<, --^L^-M. = Ak^• (9)

rpn rpn

-k+1,i — -k ,i

(s-jn

U J k ,i

К

rpn rpn

-k ,i+1 — -k ,i

(6) • (7)

Если выполняется условие

\Ak J[ > Кi|+| Ak+w|.

(10)

то итерационный процесс:

at! (s +1) =

A,

A,

k'+1 ak++1 (s) +

A

k+1,i n+1

A

i (s) —

, rr^n +1 _ rrin

1 - k ,i — - k,i

A,,. т

(11)

Уравнение (5) с условиями (6) и (7) позволяет определить коэффициент а"к1 во всех точках, кроме

А - значения в этой точке могут быть определены (например, путем линейной интерполяции), записывая (5) в виде:

/ \ n+1 (—1 "(ST j n+1 "

U J k ,i+1 an+1 + ak ,i+1 + ^дх J k+1,i

К К

Т'П+1 rpn

Ь; —-Ь i

ST дх

К

— j К

(8)

Вводя обозначения:

A,

A,

ST

¥

К

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

дх

К

дТj (дТ

А _ хдх Jk,i , 1дУ Jk,i At ,■ = : I

запишем (8) в виде:

является сходящимся.

Следует отметить [8], что условие (10) является достаточным условием сходимости итерационного процесса, то есть в некоторых случаях итерационный процесс (11) совпадает при условии, что (10) не выполняется. Для ускорения сходимости итерационной процедуры используется метод верхней релаксации при реализации которого делается коррекция процесса (11):

а;;1* (*+1) = а;;1* (*)+о (а-1 (,+1) - <'* (,)), (12)

где * - номер итерации, аЩ (* +1) - значение, вычисленное по (11); а;;1 (*) а;;1 (*)- значение, скорректированное по (12). В практических задачах 1 < о < 2, что соответствует методу верхней релаксации. Если итерационные методы не могут быть использованы, то, вводя невязку уравнения (9), где К и N - количество точек разбиения х и у, к = 0,...,К-1;

Sn+1 = A an+1 + A an+1 — uk ,i Ak ,i+1ak ,i+1 + Ak+1,iak+1,i

rpn+1 _ rpn

- k, i — - k, i

— Ak X+1

i = 0,...,N — 1. (13)

Величины аП+'i, могут быть найдены из условия:

12 ^ min (14)

k=0 i=0 "n+ Указанная задача имеет единственное решение, поскольку:

5 (< ) = Z

n+12 k i

(15)

является положительно определенной функцией, суммой квадратов линейных функций. Записывая необходимые условия экстремума

SS (аП+1)

дак+1

= 0

(16)

можно получить систему ^ +1) (К +1) -1 уравнений относительно N х К -1 неизвестных а£+1, к = 1,... , К; 1 = 1,... ,N. Изображая точки расположения невязок 5к; (рис. 2), получаем:

т

т

fin+1 — A un+1 + A u uk,i Ak,i+1uk,i+1 + Ak+1,iuk+1,i

rpn+1 rpn n+1 k ,i k, i * n+1.

- Ak,iak,i ■

T

rpn+1 rpn

fin+1 - A un+1 + A un+1 - k-1'i ~ k-1'i - A un+1 ■ °k-1,i — Ak-1,i+1ak-1,i + Ak,iuk,i Ak-1,iuk-1,i ■

T

n+1 n

fin+1 — a un+1 + A un+1 - k'i-1 kJ-1 - A un+1 °k,i-1 — Ak,iuk,i + Ak+1,i-1uk+1,i-1 Ak,i-1uk,i-1-

T

(17)

Рис. 2. Схематическое изображение распределения температурных коэффициентов

Таким образом, я^1, к 1.....К 1: 1=1,...,N-1

входит только в три слагаемые суммы (15). Величина аП+д входит только в одно слагаемое:

fin+1 — a un+1 + A

°s,N-1 — As,Nus,N + Ai

rpn+1 _ rpn

1 s, N-1 - 1 s, N-1

Л+1, N-1" s+1,N-1

- A un+1

As, N -1us, N -1

(23)

грН+\ грп

Яп+1 = А ап+1 + А ап+1 -100 0'0 - А ап+1

^0,0 ^Д^Д ^^1,01'Ч,0 0,0Ы0,0

т

Величина ап^\ входит в слагаемое:

(18)

fin+1 = Л п"+1 А- Л П"+1 -

^K-1,0 ^K-1,1" K-1,1 ^ ^K,0^K ,0

rpn+1 rpn

n+1 1 K-1,0 - 1K-1,0

— A un+1

K-1,0" K-1,0

(19)

Величина ап+м входит в слагаемое:

грН +1 rpn

Xn+l _ Л n+l I Л n+l 0,N-1 ~10,N-\ л n+\ (20)

Величины а п

уравнения:

c«+1 _ л n+1 I Д n+1

°0,i-1 — A0,iu0,i + A1,i-Aj-1-

i — 1,...,N-1 входять два

rpn +1 rpn

0,i-1 0,i-1 Д n+1 .

A0,i-1u0,i-1'

(21)

— A>,i+1a0>M1 + A1,iu1,+

Z'« + 1 rp]

n ; ~ -'n

n+1 +0,1 +0,1

— ■

Величины а"п\ , к = 1,...,К-1 входят в два уравнения:

Величины а^1 , г = 0,., N -1 входят в слага-

емое:

fin+1 — A u"+1 + A u"+1 -

UK-1,i AK- 1,i+1uK-1,i+1 + AK ,iuK ,i

rpn+1 rpn

1K- 1,i -1K- 1,i

- A u"+1

AK- 1,iuK- 1,i

(24)

С учетом (17)-(24) необходимые условия экстремума записываются в виде:

Для ак+1 , к = 1,...,К-1;г = 1,...,N-1

д

3S _

dun+1 — dun+1

duki duki - dfin+1

— k+ dun+1

№ )2+{ö:+i )2+(c-1 )2 )—

, dS?,. , ddt\

en+1 k-1,/ . en+1 k-1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+ АО, л .-— + AO, . ,

k-1,; n+1 k,t-1

du, ,.

dun

I rj-zn+1 rj-zn

— 2 I A un+1 + A un+1__—_— - A un+:

— 21 Ak,i+1uk,i+1 + Ak+1,iuk+1,i Ak,iuk,i

j(" ^^ , у

7m + 1 rpn

П 1

n+1 n+1 n+1

°k-1,0 — Ak-1,iuk-1,1 + Ak ,0uk ,0 '

rpn+1 rpn

— л n+i I л n+i °k,0 — Ak ,1uk ,1 + Ak+1,0uk+1,0

n+1

Ak-1,0uk-1,0'

-4,0<+01- (22)

Величины , 5 = 1,...,К-1 входят в одно

слагаемое:

I rpn+1 rpn ^

+2 ^ Ak-1,+1utl+1 + Ak .uz1 -Ak-1^+1 j Ak,i +

^ rpn+1 _ pn ^

+21 AkOZ1 + Ak+,t_-Jl^lZJl^L-Ak^_iuk+_1 I x

x Ak ,i — 0-

(25)

T

T

T

T

T

T

В уравнении 7 неизвестных величин (запишем их опуская верхний индекс «п+1»): аА(+

ак +1,/; ак ,1 ; ак-1,1+1; ак-1,1; ак+1,1 -1; ак ,1 -1-Для а^1:

dS _d« )2

da"+1 dan+1

f rpn+1 _ rpn Л

- 21 A an+1 + A an+1 - ~ ■ "" 1

- 2 I A01a°1 + A1,°a1°

1 , , , , z

Для a"K+1:

d(& )2

dS

daK ,0 ^aK ,0

2 I AK-1,1aK-1,1+ AKaK ,1

x AK,° -

1 TK-l,° TK-1,C

- A /7n+1 AK-1,°aK-1,°

(27)

Для a;+1, 5 - °,1,..., К-1:

dS

daTN

- 2 (As,

d(C-1 )2

da

n+1 -*s, N

Г' + A.

s+1,N-1" s+1,N-1

A°,°an;1 j(-A°,°)- °.(26)

rpn+1 _ rpn

Ts+1,N-1 - Ts+1,N-1

x As,N - °.

Для a"K+1, г -1,...,N-1:

- A an+1

As, N -1as,N-1

dS d(SK% )2

dan+1 dan+1 daK ,1 daK,i

(

2 | AK-1,i+1aK-1,i+1 + AK,iaK,1

(31)

rrin + 1 rj-in Л

n+1 TK-1,1 - TK-1,1 1

— A an

AK-1,iaK-1,1

В уравнениях (26) и (27) по три неизвестных в

(26):

a°,1; a1,°; a°,G ' в уравнении (27) aK-1,1, aK,° и aK_1,0.

Для a°n+N:

dS _ dfe-1 )2

dan+1 dan+1

(

j

rpn + 1 rpn

I A an+1 + A an+1 - T°,N-1 - T°,N-1 - a an+1 I

I °,N ^ Л1,N-1"1,N-1 °,N-1"°,N-1

^ , , , , Г )

x(-A°,n-1)- °.

В (28) - три неизвестные: a0N; aN_x; a0N_j.

Для a^1, i -1,..., N -1;

dS _d((c 1 )2+(c)2)_

(28)

da0n;1

dan

(

rpn+1 rj-zn ^

- 21A°,ia"o,i + A1,i-1a1ni-1 , " , A°,i-la<0,/-o1 J Л,1 +

f rrin+1 rr^n Л

„n+1 , Л „n+1 T°,i - T°,i л n+1

A°,ia°,i

+21 Ai+Ä + Axa+ -

:(-A°,i)- °,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для a°+01, к -1,...,K-1:

dS _d((c1,°)2 +(С )2)_

(29)

da;

k,° (

dal

Z'n+1 rpn и 1 л T 1,

- 21 Ak-1,1ak°^1,1 + Ak,° a"nß

f rpn+1 rpn

+21 A an+1 + A an+1 - '' +21 Ak,1ak,1 + Ak+1,° ak+1,°

. _ Л nn+1 I A _i_ Ak-1,° ak-1,° Ak,° +

Ak,° an,°

,J

x(-Ak,°)- (3Ö)

неизвестные величины: ak_x x;ak 0; ak_x 0;a,; ak+l 0;

xAk,I -

(32)

Неизвестные величины в (32 ): aK_lj+1; aK j; aK X j.

Таким образом, система линейных алгебраических уравнений (25)-(32) состоит из:

N * - (N -1)( K -1) +1 +1 +1 + N -1 +

+K --1 + K -1 + N -1 - NK - K -

- N +1 +1 +1 +1 + 2 N + 2K - 4 -

- NK + N + K - (N + 1)(K +1)-1 (33)

уравнений относительно (N +1)(K + 1)-1 неизвестных величин. Указанная система имеет единственное решение, если ее определитель не равен нулю. Система (25)-(32) позволяет восстановить поле a(x,y).

5. Выводы

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

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

- восстановление температурного поля площади поверхности регенератора по известным его значениям на некотором множестве точек поверхности;

- проведение экспериментальных исследований температурных полей поверхности регенератора;

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

Лггература

1. Jonsson, M. Humified gas turbines - a review of proposed and implemented cycles [Text] / M. Jonsson, J. Yan // Energy. - 2005. - Vol. 30, Issue 7. - P. 1013-1078. doi: 10.1016/j.energy.2004.08.005

2. Петрш, С. Я. Причини i фактори, що зобумовлю-ють змшу техшчного стану теплообмшних апаратш ГПА [Текст] / С. Я. Петрiв // Науюж вют! - 2012. - Вип. 2 (22).

3. Усеров, А. Г. Используемые аппараты и разработка новых эффективных регенераторов ГПА для повышения

т

неизвестные: as.N ; as+1N-1 ; as.N-1

неизвестными являются величины: a ; a ; a ;

a°,i+1; a1,i.

КПД и утилизация высокопотенциальной теплоты с целью снижения выбросов в атмосферу [Текст] / А. Г. Усеров, К. К. Шалбаев // Алмата, Вестн. НАН РК. - 2009. - № 2.

4. Козачанко, А. Н. Энергетика трубопроводного транспорта газов [Текст] / А. Н. Козачанко, В. И. Никишин, Б. П. Поршаков. - Гуп Издательство «Нефть и газ» РГУ нефти и газа им. И.М. Губкина. Москва, 2001. - С. 327-355.

5. Коздоба, Л. А. Методы решения нелинейных задач теплопроводности [Текст] / Л. А. Коздоба. - М.: Наука, 1975. - 170 с.

6. Тихонов, А. Н. Методы решения некоректных задач [Текст] / А. Н. Тихонов, В. Я. Арсенин. - М.: Наука, 1979. - 285 с.

7. Дьяченко, В. Ф. Основные понятия вычеслитель-ной математики [Текст] / В. Ф. Дьяченко. - М.: Наука, 1977. - 120 с.

8. Андерсон, Д. Вычеслительная гидромеханика теплообменника. Т. 1 [Текст] / Д. Андерсон, Д. Таннехел, Р. Плетчер. - М: Мир, 1990. - 384 с.

References

1. Jonsson, M., Yan, J. (2005). Humified gas turbines -a review of proposed and implemented cycles. Energy, 30 (7), 1013-1078. doi: 10.1016/j.energy.2004.08.005

2. Petriv, S. Y. (2012). Reasons and factors which contribute to changing the technical condition of heat exchangers compressor units. Science news, 2 (22).

3. Userov, A. G., Shalbaev, K. K. (2009). Used devices and development a new effective regenerators of compressor units for improve efficiency and utilization of high-grade heat in order to reduce emissions into the atmosphere. Almaty Herald. NAS RK, 2.

4. Kozachanko, A. N., Nikishin, V. I., Porshakov, B. P. (2001). Power engineering pipeline transport gas. PMU Publishing "Oil and Gas" Oil and gas RSU th. I. M. Gubkin. Moscow, 327-355.

5. Kozdoba, L. A. (1975). Methods for solving nonlinear heat conduction problems. Moscow: Science, 170.

6. Tikhonov, A. N., Arsenin, V. Y. (1979). Methods of solving incorrect problems. Moscow: Nauka, 285.

7. Diachenko, V. F. (1977). Basic concepts of computational mathematics. Moscow: Nauka, 120.

8. Anderson, D., Tannehel, D., Pletcher, R. (1990). Computational fluid dynamic of heat exchanger. Moscow: Mir, 1, 384.

9. Marchuk, G. I. (1989). Methods of computational mathematics. Moscow: Nauka, 608.

Дата надходження рукопису 17.03.2015

Замиховский Леонид Михайлович, доктор технических наук, профессор, заведующий кафедры, кафедра компьютерных технологий в системах управления и автоматики, Ивано-Франковский национальный технический университет нефти и газа, ул. Карпатская, 15, г. Ивано-Франковск, Украина, 76019 E-mail: leozam@ukr.net

Петрив Светлана Ярославовна, аспирант, кафедра компьютерных технологий в системах управления и автоматики, Ивано-Франковский национальный технический университет нефти и газа, ул. Карпатская, 15, г. Ивано-Франковск, Украина, 76019 E-mail: Svetlanapetriv@yandex.ru

УДК 004.896

DOI: 10.15587/2313-8416.2015.41213

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

© Л. М. Замиховский, Р. М. Матвиенко

Обоснована необходимость построения интеллектуальной системы поддержки принятия решений (ИС-ППР) при управлении технологическим процессом компримирования газа с включением в ее состав учебно-тренировочного модуля. Приведена структурная схема человеко-машинного взаимодействия в системе "Диспетчер - ИСППР - САУ ГПА - КЦ" и определены функции отдельных блоков интеллектуальной системы поддержки принятия решений при управлении процессом компримирования газа Ключевые слова: интеллектуальная система, система поддержки принятия решений, диспетчер, оператор, процесс компримирования газа

The necessity to construct the intelligent decision support systems (IDSS) in the control of the gas compression process with inclusion in structure of the training module is substantiated. А block diagram of human-computer interaction in the system "Dispatcher - IDSS - ASC GCU - CS" is shown, and the functions of the separate blocks of the intellectual decision support system in the control of the gas compression process are defined Keywords: intelligent system, decision support system, dispatcher, operator, gas compression process

1. Введение

Современная газотранспортная система имеет свои специфические особенности:

- большую территориальную протяженность;

- значительное число элементов, которые формируют систему, иерархическую структуру;

- централизованное управление технологическим процессом транспортировки целевого продукта;

- децентрализованное распределение целевого продукта;

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

i Надоели баннеры? Вы всегда можете отключить рекламу.