CC BY
79
УДК 621.3
Н. А. Автушенко, Г. С. Леневский, канд. техн. наук, доц.
ПОСТРОЕНИЕ СИСТЕМ УПРАВЛЕНИЯ ПЕРЕХОДНЫМИ ПРОЦЕССАМИ В МАГИСТРАЛЬНЫХ ТРУБОПРОВОДАХ ГОРЯЧЕГО ВОДОСНАБЖЕНИЯ
В статье рассмотрена методика расчета регулятора давления в системах горячего водоснабжения, разработана математическая модель в среде МЛТЬЛБ. Расчет основан на преобразовании гиперболических передаточных функций трубопровода в инерционное звено второго порядка. Полученное решение позволяет контролировать выходной параметр в любой точке магистрального трубопровода горячего водоснабжения. Такое решение направлено на исключение гидравлических ударов, повышение надежности систем и, как следствие, продление эксплуатационного срока.
Теоретические исследования в области систем автоматического управления регулируемых электроприводов турбомеханизмов, работающих на длинные трубопроводы, вызваны ростом числа энергосберегающих технологий, а также необходимостью автоматизации работы насосных станций и распределительных систем трубопроводов. Между тем, применение частотных приводов насосов на ряде объектов требует специального подхода по причине возникновения гидравлических ударов и колебаний давления в трубопроводе при регулировании частоты вращения турбомеханизма.
Возникновение колебаний обусловлено тем, что длинные трубопроводы являются объектами с распределенными параметрами. Изменение давления в любой точке вызывает волну давления, распространяющуюся по трубопроводу со скоростью, определяемой скоростью звука в среде трубопровода. Колебания давления возникают также при наличии в системе трубопроводов замкнутых контуров или объектов, отражающих волну давления.
Принимается, что давление в обратном сетевом трубопроводе постоянно и равно 0,2 МПа, давление на входе сетевых насосов второго подъема формируется насосами первого подъема, изменение
давления в обратном сетевом трубопроводе близко к нулю и, как следствие, не учитывается в дальнейшем расчете. Таким образом, контроль перепада давления вырождается в контроль давления в конкретной точке прямого магистрального трубопровода. Напряжение задания всех частотных приводов насосов остается одинаковым.
В качестве исследуемого объекта используется участок магистральной тепловой сети Могилевских тепловых сетей от ТЭЦ-2 до павильона П2 общей протяженностью 4173,6 м. Система имеет следующие параметры: количество работающих насосов второго подъема - 4, тип насоса СЭ-2500-140
(Он = 2500 м3, рН = 1,37 МПа), суммарный расход в прямом коллекторе ТЭЦ-2 - 9087 м , расход исследуемой магистрали - 4479 м , номинальное давление в прямом коллекторе - 1,2 МПа, условный диаметр трубопровода - 1 м, температура наружного воздуха - 0 оС, температура теплоносителя - 98 оС.
Для расчета регулятора в конце трубопровода необходимо рассчитать передаточные функции насосного агрегата и магистрального трубопровода.
Расчет передаточной функции насосного агрегата имеет ряд трудностей, т. к. содержит квадрат входного воздей-
ствия - расхода воды (скорости вращения крыльчатки): Н = a■Q2 + Ь-() + с; (Н = 0,0000037^2 - 0,0023330^ + + 169,1289177). Для получения искомой передаточной функции частично используется расчет А. К. Аракеляна и А. В. Шепе-лина [1].
(р) = ^ = *н 1 + ^
ю(р)
1 + Хкп • е-2тр
.13
е-2т - 1 _ 2тр + 2т2р2 - ^3^, (1)
где т - время распространения волны в один конец трубопровода:
I
т = ■
(2)
кн - коэффициент линеаризации кривой напора насоса от скорости:
кн = -
Рн.а.ном Рн.а.нач
ю - ю
н.а.ном н.а.нач
(3)
Далее в расчетах используется модуль значения (3).
Коэффициент степени колебательности системы X для инженерных расчетов принимается равным 0,6 и определяется по формуле
X = 1 -
2кт • А
рс
(4)
где А - площадь сечения трубопровода, Ь = жё2/4; кнх - коэффициент жесткости характеристики в области рабочей точки.
Учитывая вышестоящую формулу, может быть найден коэффициент жесткости характеристики в области рабочей точки
= (1 - х )рс нх 2 • А
(5)
Коэффициент отражения волны от конца трубопровода
кА -1
к =-РС
ко = А
(6)
к--------+1
Рс
Коэффициент, учитывающий степень закрытия задвижки в конце трубопровода:
, Рх-1
к . (7)
а
X=1
Передаточная функция насосного агрегата примет вид:
ш =
гг Н
7750,16 0,13 р +1
(8)
Передаточная функция магистрального трубопровода, полученная из уравнений Навье-Стокса и уравнения неразрывности среды Ж^, /), и преобразованная передаточная функция в виде, удобном для моделирования в среде МЛТЬЛВ, Ж]^, /) определяются как
1
2 ВЛ (5)
22 (5)
5^[0( 5), I ] + с^[9( 5), I ]
Ш1(5,1) = ■
с^[9( 5)/] + -
к____
2ВЛ (5)
^( 5)
5к[0(5)1]
где 21(£), 22(£), ХвАъ) - концевые операторные сопротивления; К - поправочный коэффициент, К = _/(/).
При оценке поведения динамической составляющей давления теплоносителя трубопровода было замечено сходство в получении передаточных функций [2] с квадратичными полиномами п-го порядка. Наряду с использованием регулятора давления с помощью гиперболических функций [3] искомую передаточную функцию можно представить в полиномиальном виде. Передаточная функция магистрального трубопровода имеет вид:
с
.(9)
Для исследуемого участка протяженностью 4173,6 м передаточная функция примет вид:
Шр 2 р1 (5 ) =
1
(4,5452 +1)1,8552 +1)1,1152 +1)
Анализ распределения резонансных частот выполняется с помощью логарифмических амплитудно-частотных характеристик (ЛАЧХ). При 5 = ]ш ЛАЧХ пе-
редаточных функций Ж(я,/), Ж](я,/) и Жр2р1{^), Ж](?) имеют вид:
1(ф) - 20^ [0&, /)] ;
Ь1(т) = 20(§- [( ja>, I)] ;
Ь2{а>) = 20/^\Жр 2 р1Ца)) ] .
Анализ распределения резонансных частот с помощью ЛАЧХ представлен на рис. 1.
0,01 0,1 Гц 1
ш
---------------->-
2п
Рис. 1. ЛАЧХ магистрального трубопровода по давлению
Учитывая рекомендации по выбору постоянных времени [2], используем только вторую резонансную частоту как базовую. Реально в регуляторе оператор Лапласа со степенью выше второй влияния на систему не оказывает - оставляем множитель знаменателя от базовой резонансной частоты юРЕЗ = 1Д/1,85 . Кроме того, если систему оставить в виде
1,85 52 + 1
-
то очевидно,
р 2 р14173,6 (5)
что это передаточная функция для незату-
хающих колебаний, хотя в реальности колебания имеют затухающий характер. Для этого по аналогии с [2] вводится демпфирующая составляющая.
В общем виде передаточная функция (9) примет вид:
1
(р 52 + 2 • у • ТТР +1)
(10)
где ТТР - постоянная времени трубопровода ?2рЕЗ = VЮ2РЕЗ .
Для исследуемого участка функция примет вид:
Wp2pl (s) —
1
(l,85s2 + 2 .у VTT85 +1)
(11)
где у - параметр аппроксимирующей модели объекта с распределенными параметрами (РП - объекта) [2].
Обратная связь по давлению может быть рассчитана следующим образом:
КДР(р) —
-L ГГ^
(12)
Расчет параметров выполняется исходя из рекомендаций [2] по принципу подчиненного регулирования. Малая не-компенсируемая постоянная времени разрабатываемой системы автоматического управления (САУ) давлением в трубопроводе выбирается исходя из поиска наименьшей постоянной времени элементов системы. Таковой является постоянная времени преобразователя частоты ТТП. Малая некомпенсируемая постоянная САУ примет вид:
T — 2. T
1. ^ 1ТП ■
(13)
Передаточная функция регулятора, настроенного на модульный оптимум (МО):
ш (5) = (ТР5 + 2 • У • ТТР + Мдс (тН / А5 + 1) (14)
Р 2 Т К К К К 5 V/
/л ДР Г2^Н / АЛ
Для получения оптимального переходного процесса устанавливается входной фильтр, а САУ настраивается на симметричный оптимум. Тогда передаточная функция регулятора и фильтра определится как
Wф (S) —
1
Т* s+1
(15)
Передаточная функция регулятора давления, настроенного на симметричный оптимум (СО), примет вид:
WPEГ ( S ) —
abT. s +1 abT. s
(тТрs2 + 2У-Ттр +1)Кдс (Тн/АS + 1)
2 • Т.КДР КТЕЭКГ2Кн /А S
. (16)
Для моделирования в среде МЛТЬЛВ регулятор преобразовывается таким образом, чтобы показатель степени оператора Лапласа знаменателя был не менее показателя степени оператора Лапласа числителя (реально каждое из слагаемых третьего множителя знаменателя с учетом его порядка малости, за исключением 1, вырождается в 0 и не оказывает влияния на переходной процесс).
^^^ег (S) — ■
abTu s +1 ___№____.
abTs
, КДС (Тн / А S + 1)
2 • Т.КДР КTES КГ2 Кн / А S
ТТр s 2 + 2 • у • ТТР +1 0,00001s2 + 0,0001s +1.
(17)
Передаточная функция фильтра для симметричного оптимума
WФСО(s) —
1
a№ s+1
М
(18)
Pегулятор давления, настроенный на СО (блок REG), представлен на рис. 2, математическая модель САУ давлением теплоносителя в трубопроводе тепловых сетей в среде MATLAB представлена на рис. 3.
Анализ параметров регулятора. Pассматривается параметр а в диапазоне от 1,2 до 8. Графики переходных процессов давления в павильоне П2 исследуемого трубопровода для различных а представлены на рис. 4 и 5.
Из полученных графиков следует, что изменение параметра а даже в большом диапазоне, когда система обладает довольно значительной инерционностью, на переходной процесс влияния практически не оказывает.
X
X
1З2
Pna 2. Математическая модель блока REG (регулятора давления) электропривода и системы управления поведением теплоносителя в магистральных трубопроводах тепловых сетей в среде MATLAB
Так как время моделирования (^ПП = 0,95^ст) принимается по аналогии с переходным процессом, описанным с использованием гиперболических функций, -30 с, то предварительно ТФ = ПП/3,4.
Выполняется моделирование для оценки влияния фильтра на входе в САУ. Результаты моделирования представлены на рис. 5.
График переходного процесса (рис. 6) используется для сравнения с переходным процессом, полученным для гиперболических функций, откуда видно, что оптимальным является переходной процесс с параметром фильтра ТФ = 4хТТР.
Оценка качества переходного процесса выполняется исходя из условия попадания в 95 % от установившегося значения давления в исследуемой точке при времени переходного процесса 30 с.
По результатам моделирования видно, что более предпочтительным является постоянная времени фильтра ТФ = 4хТТР = = 7,4 с.
Моделирование выполняется по схеме, представленной на рис. 3, система регулирования выполнена по схеме, представленной на рис. 2 с настройкой регу-
лятора давления на МО (отсутствуют блоки ге§р1, ге§р3).
Моделирование с настройкой регулятора давления на СО выполняется по схеме, представленной на рис. 3. Результаты моделирования представлены на рис. 7.
При использовании регулятора по давлению, рассчитанного по формулам (15), (17), (18), колебательность заметно снижается в сравнении с нерегулируемой схемой.
Далее выполняется проверка полученной системы регулирования для различных сигналов задания. Данные соответствия сигналов заданий и требуемых давлений для исследуемой системы представлены в табл. 1. Пересчет выполняется по (12).
Результаты моделирования переходных процессов для СУЭП в магистральных трубопроводах с внешней обратной связью и регулятором, настроенным на СО для заданий из табл. 1, представлены на рис. 7.
Результаты моделирования совпадают с заданием с требуемой точностью в 1 %.
Электротехника
Рис. 3. Математическая модель САУ давлением теплоносителя в трубопроводе тепловых сетей в среде МАТЬАВ
Вестник Белорусско-Российского университета. 2009. № 3(24)
t------------------>•
Рис. 4. Графики переходных процессов давления в трубопроводе в павильоне П2 при изменении параметра а в диапазоне от 1,2 до 8
0 50 100 с 150
t---------->►
Рис. 5. Графики переходных процессов давления в трубопроводе в павильоне П2 при пуске четырех насосов при изменении постоянной времени фильтра (ТФ = 1хТТР ... 5хТТР )
10
105-Па
8
Р р2 гед 3,7
/ Р р2 гед 5,5
Р р2 по геа' Р р2 гед 7,4
10
15
20
t ■
30
Рис. 6. График переходного процесса давления в П2 для нерегулируемой и регулируемой схем с
Тф = 2хТтр = 3,7 с, Тф = ЗхТтр = 5,5 с, Тф = 4хТтр = 7,4 с
9
105-Па
7
20 40 60 80 100 120 140 с 180
t ■
Рис. 7. Графики переходных процессов давления в павильоне П2: 1 - одновременное включение четырех насосов по регулируемой схеме в внешней ОС по давлению с настройкой на СО; 2 - веерное включение четырех насосов по нерегулируемой схеме
6
4
2
0
5
с
5
3
1
0
Табл. 1. Соответствие сигналов заданий и давлений в П2
Сигнал задания напряжения САУ, В Давление в П2, Па
10 845590
9 762712
8 677966
7 593220
Для исследуемой схемы оценивается адекватность системы регулирования при моделировании штатной ситуации - изменение сигнала задания с 9 по 9,37 В при увеличении давление с 0,76 до 0,8 МПа. Моделирование выполняется при мгновенном изменении задания (нерекомен-
дуемый режим задания) (рис. 8) и плавном задании за 10 с (рис. 9). Отклик системы запаздывает на 4,2 с ввиду значительной протяженности трубопровода. На основании полученных результатов можно судить об адекватности полученного регулятора.
9
105-Па
7
1
і і !и=1СК/ ;
с II ¥
І І І и=8и і
С <
І і І І І І
0 20 40 60 80 100 120 140 с 180
г---------->-
Рис. 8. Семейство кривых давления в павильоне П2 при использовании системы регулирования с настройкой на СО
5
3
Выводы
Полученная система позволяет контролировать ОН-параметры в любой точке трубопровода, исключая нестационарные процессы в системе.
Система может быть использована для прогнозирования результатов конкретного режима, а также для исключения гидравлических ударов и колебания напора в коммуникационной сети путем формирования переходных процес-
сов пуска и останова системы электропривода группы насосных агрегатов. Представление длинного трубопровода
как динамического объекта позволит находить оптимальные законы управления системой.
г--------->-
Рис. 9. График переходного процесса давления в павильоне П2 при одновременном пуске насосных агрегатов при плавном увеличении сигнала задания с 9 до 9,37 В за 10 с
СПИСОК ЛИТЕРАТУРЫ
1. Аракелян, А. К. Способы построения систем автоматического управления электроприводами насосов, работающих на длинные трубопроводы / А. К. Аракелян, А. В. Шепелин // Электричество. - 2001. - № 2. - 212 с. : ил.
2. Киселев, Н. В. Электроприводы с распределенными параметрами / Н. В. Киселев, В. Н.
Мядзель, Л. Н. Рассудов. - Л. : Судостроение, 1985. - 220 с. : ил.
3. Автушенко, Н. А. Построение математической модели системы горячего водоснабжения магистральных трубопроводов в среде МЛТЬЛБ / Н. А. Автушенко, Г. С. Леневский // Вестн. Белорус.-Рос. ун-та. - 2007. - № 3. -212 с : ил.
Белорусско-Российский университет Материал поступил 27.09.2008
N. A. Avtushenko, G. S. Lenevsky Construction of control systems by transients in the main hot waterway pipelines
The design procedure of a regulator of pressure in the main hot waterway pipelines is shown in the scientific paper. A mathematical model in MATLAB has been developed. Account is based on transformation of hyperbolic transfer functions of the pipeline to a second order relaxation circuit. The received decision allows supervising output parameter in any point of the main hot waterway pipelines. It is directed on exception of hydraulic impacts, reliability of systems and operational term improvement.
