Электротехника
УДК 62-83:622(075)
АНАЛИТИЧЕСКОЕ И ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ СТАЦИОНАРНЫХ РЕЖИМОВ РАБОТЫ УСТАНОВОК ОХЛАЖДЕНИЯ ГАЗА КОМПРЕССОРНЫХ СТАНЦИЙ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ
А.М. Абакумов, С.В. Алимов, Л.А. Мигачева
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: [email protected]
На основе аналитических и экспериментальных исследований разработана модель для расчетов оптимальных стационарных режимов работы установок охлаждения газа. Проанализированы результаты, достигаемые при использовании различных алгоритмов управления частотой вращения двигателей вентиляторов аппаратов воздушного охлаждения.
Ключевые слова: установки охлаждения газа, оптимизация, электропривод вентиляторов.
В условиях дефицита электроэнергии важное место в повышении экономической эффективности предприятий Газпрома занимают мероприятия по энергосбережению на всех стадиях транспортировки газа.
Современная технология транспортировки газа по магистральным трубопроводам требует его охлаждения до определенной температуры. С этой целью на газотранспортных предприятиях - на установках комплексной подготовки газа, на компрессорных станциях магистральных газопроводов, дожимных компрессорных станциях, станциях подземного хранения газа и т.д. - широкое применение нашли аппараты воздушного охлаждения (АВО) газа. На компрессорных станциях магистральных газопроводов с газотурбинным приводом нагнетателей до 60...70% потребляемой электроэнергии расходуется на охлаждение газа, что определяет актуальность проблемы повышения энергоэффективности работы АВО.
Для разработки математической модели, ориентированной на оптимизацию алгоритмов управления АВО, целесообразно провести декомпозицию объекта управления, представив его в виде взаимосвязанных подсистем: процессов теплообмена, аэродинамических процессов и электромеханических процессов в электроприводах вентиляторов.
При проведении исследований рассматривалась схема АВО в виде двух последовательно включенных секций (частей теплообменника) 1 и 2 (см. рисунок) с возможностью независимого регулирования частоты вращения n1, n2 вентиляторов (В1,
Александр Михайлович Абакумов - д.т.н., профессор.
Сергей Викторович Алимов - аспирант.
Людмила Алексеевна Мигачева - к.т.н., доцент.
В2), приводимых в движение электроприводами (ЭП1, ЭП2). В стационарном режиме перепад температуры в частях теплообменника AT, является функцией следующих параметров: температуры газа Tin на входе, температуры наружного воздуха Toa, расхода газа Q, расхода или скорости потока V¡ охлаждающего воздуха i-й секции теплообменника, скорости потока воздуха смежной секции. Последнее объясняется тем обстоятельством, что во вторую часть теплообменника поступает газ, частично охлажденный в первой части теплообменника.
Общий перепад температуры на АВО
Tin ~ Твых = ДТ1 + ДТ2 = •
Тоа Q V Тоа Q Уг
Расчетная схема АВО
Известно, что при постоянном угле атаки лопастей ф потребная мощность на валу вентилятора зависит от относительной скорости в третьей степени
P = P
c ± c.nom
П
V nom J
(і)
где ппот - номинальная частота вращения, Рс.пот - номинальная мощность на валу вентилятора при номинальной частоте вращения.
Экспериментальные исследования, выполненные для аппаратов воздушного охлаждения газа типа 2АВГ-75, оборудованных двумя вентиляторами, приводимыми в движение электродвигателями типа ВАСО-16-14-24, подтверждают справедливость соотношения (1).
Управление аэродинамическими процессами осуществляется за счет изменения частоты вращения п вентиляторов с помощью частотно регулируемого привода. Результаты экспериментального исследования аэродинамических характеристик свидетельствует о возможности представления взаимосвязи скорости потока охлаждающего воздуха V и частоты вращения п (или частоты /напряжения на двигателе) в виде линейной зависимости
V = К" > (2)
где кьУ - коэффициент передачи вентилятора по частоте вращения двигателя.
При экспериментальных исследованиях тепловых процессов рассматривались обобщенные характеристики АВО, устанавливающие взаимосвязь между регули-
рующим воздействием - частотой напряжения, подаваемого на двигатель, и выходной переменной - перепадом температуры AT на АВО.
В процессе экспериментов исследовались зависимости перепада температур от регулирующего воздействия как при синхронном регулировании частоты вращения обоих вентиляторов, так и при отключении одного из вентиляторов и регулировании частоты другого.
На основе анализа полученных данных перепад температуры на АВО для управляющих воздействий f и V представлен в виде
\с + kff при f = 0 AT = \ 1 (3)
[C + kfifi + kf2/2, при f > 0.
ÍC + kLV], при V2 = 0 AT = j vl 1 ^ 2 (4)
[C + kvVi + kv2V2, при V2 > 0.
где C - перепад температуры при отключенных двигателях вентиляторов; k^ , k^ ,
k/2, k^j, kvi, kv2 - коэффициенты передачи для управляющих воздействий в виде
частоты f напряжения и скорости потока воздуха V первого и второго вентиляторов соответственно.
Результаты экспериментальных исследований показывают, что коэффициенты передачи для управляющих воздействий f и V¡ на первую часть теплообменника отличаются от соответствующих коэффициентов передачи для управляющих воздействий f2 и V2 на вторую часть теплообменника. При этом возникает задача отыскания оптимальных значений управляющих воздействий.
В типовой установке охлаждения газа (УОГ) содержится n параллельно включенных АВО. При допущении о равномерном распределении потока газа по отдельным АВО и однородности их тепловых характеристик общий температурный перепад на УОГ совпадает с температурным перепадом на отдельном аппарате воздушного охлаждения.
В соответствии с требованиями технологического задания перепад температуры на УОГ с учетом некоторой погрешности £ должен быть равен заданному ATzd значению
AT±£ = ATzd. (5)
При этом на частоту вращения вентиляторов и скорость потока воздуха наложены ограничения
n ■ < n < n * V ■ < V < V * 7—12 (6)
"min 1 ^ "max 1 ’ Y min 1 1 ^ Y max 1 ’ 1
При определенных значениях параметров T7n,Toa, Qg и номинальных значениях
частоты nnom7 вращения вентиляторов и, соответственно, номинальных значениях
скорости потока воздуха Vnom7 УОГ обеспечивает некоторое максимальное значение перепада температуры
C + kvlVnoml + kv2Vnom2 = ATmax ' (7)
С учетом соотношения (1) выражение для мощности на валу вентиляторов мож-
но записать в виде
PZ = PnomÁT^)3 + Р«от2(т^М3- (8)
V noml V nom2
Дополнительно будем считать номинальные мощности вентиляторов и номинальные значения скорости воздуха одинаковыми:
P = P = P •
noml nom2 nom ’
Vnomi Vnom2 Vnom *
(9)
Для получения обобщенных результатов перейдем к относительным единицам. За базовые значения примем: температуры - максимальное значение перепада температуры АТтса, мощности - Pnom; скорости воздуха - Vnom-Введем обозначения для относительных величин:
Ах =
AT
C
ATm
c = -
ATm
ki =
k
vl
AT„
k2 =
v2
ATm
kl =
k
vl
AL
V
Vo =■
; p = -
p
V V P 50 n
nom nom nom '1nom
* r *
; f = — ; n=-50
: f = n = V *
Ax = <
(10)
После преобразований соотношений (4), (5), (6), (8) в относительных единицах будем иметь
с + к/ у1, при V2 = 0 [с + &! V! + к2У2, при V2 > 0.
Ах = Ахтр ±4; (11)
(12) (13)
V ■ < V < V '
Fmin г — — — Fmaxг ’
с + + ^2 — 1;
рГ = ^3 + ^2 • (14)
При этом задача оптимизации стационарного режима работы УОГ может быть сформулирована как задача минимизации целевой функции (14) при выполнении требований технологического задания (11), уравнений связи (10), (13) и ограничений (12). Сформулированная задача представляет собой многопараметрическую задачу нелинейного программирования, и для её решения могут использоваться известные методы [1].
При принятой упрощенной постановке задача может быть решена аналитически. Рассмотрим общий случай, когда у1>0, у2>0.
Используя соотношения (11) и приняв ¿=0, найдем с учетом выражения (10) переменную
Введем обозначения:
a =
Axmp - С k1 V2 =-------------------1V1 .
k2 k2
Axmp ~ (ci + C2\ b_ k1
k2 ’
(15)
(16)
представим (15) в виде
\2 — а - Ъ\х • (17)
После подстановки выражения для у2 в соотношение для целевой функции (14) будем иметь
Pi = V3 + V3 = Vi3 + (a - bVi )3 = Vi3 + a3 - 3a2bVi + 3ab2V^ - b3V^ = = Vi3 (i - b )a3 - 3a2bVi + 3ab2Vi2
(18)
В результате исследования выражения (18) на экстремум найдены соотношения для оптимального значения переменных:
n
^ = а(Ъ3 ^ ; у2 арг = а - Ъ^1арг • (19)
Ъ — 1
Соотношения (19) справедливы для 6^1. Значение 6=1 соответствует частному случаю, когда к1 =к2
^\орЛ = У2 аР1 = а (20)
Из изложенного следует, что в случае к1фк2 оптимальное решение, соответствующее минимуму мощности рорЬ достигается при неравных значениях управляю-
щих воздействий. Как показывает анализ, при уменьшении отношения к1/к2<1 оптимальное значение v1opt уменьшается, а v2opt возрастает.
С использованием разработанной методики решены задачи оптимизации алгоритмов управления при наличии запаса тепловой мощности АВО.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Химмельблау Д. Прикладное нелинейное программирование / Д. Химмельблау. - М.: Мир, 1975. - 534 с.
Статья поступила в редакцию 3 сентября 2010 г.
UDC 62-83:622(075)
OPTIMIZATION OF COOLING UNITS STATIONARY OPERATING MODES OF MAIN GAS PIPELINES COMPRESSOR STATIONS
A.M. Abakumov, S.V. Alimov, L.A. Migacheva
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
On the basis of analytical and experimental researches the model is developed for calculations of optimum stationary operating modes of gas cooling units. The results were reached using of various rotation frequency managing algorithms of air cooling units fans engines, that were analyzed. The results reached at use of various rotation frequency management algorithms of air cooling units fans drives are analyzed.
Keywords: gas cooling units, optimization, fan electric drive.
Alexsander M. Abakumov - Doctor of Technical Sciences, Professor.
Sergey V. Alimov - Postgraduate student.
Luidmila A. Migacheva - Candidate of Technical Sciences, Associate professor.