УДК 622.692
СНИЖЕНИЕ ПУЛЬСАЦИЙ ДАВЛЕНИЯ В МАГИСТРАЛЬНЫХ ТРУБОПРОВОДАХ ГАЗОСНАБЖЕНИЯ
Е.К. Синиченко, Ф.В. Рекач
Российский университет дружбы народов ул. Орджоникидзе, 3, Москва, Россия, 115419
Рассматривается гаситель давления, принцип действия которого основан на упруго-демпфирующем действии на поток газа.
Ключевые слова: пульсации давления, газоснабжение.
Область применения компрессоров очень широка. Мощные компрессоры применяются в нефтеперерабатывающей и нефтехимической, химической и газовой промышленности, энергетике. В газовой промышленности компрессорные машины являются основным технологическим оборудованием. Все это предъявляет чрезвычайно высокие требования к эффективности эксплуатации компрессорных машин. Часто компрессоры эксплуатируются с большим перерасходом электроэнергии, чему в немалой степени способствуют вредные пульсации давления и расхода, имеющие место на выходе из компрессора и затем передающиеся в линию. Кроме того, наличие вредных пульсаций значительно снижает КПД компрессорных машин, создает неудобства при наладке и эксплуатации установок, часто является причиной преждевременного износа оборудования. Учитывая вышесказанное, можно утверждать, что работы, связанные с разработкой и применением в компрессорном хозяйстве гасителей вредных пульсаций давления, имеют важное экономическое значение.
Проектирование рациональных систем компрессорных установок, их эффективная и надежная эксплуатация невозможны без учета режимов работы трубопроводных систем. Большое значение имеют вопросы, связанные с рациональным проектированием и эксплуатацией трубопроводных систем, находящихся под действием пульсирующих потоков газа, пара и т.д.
Теоретические основы расчета колебаний потоков и вибраций систем, вызванных возвратно-поступательным движением поршней, рассматриваются в фундаментальных исследованиях С.П. Тимошенко [1] и Д.В. Рэлея [2]. Основы неустановившегося движения сжимаемой жидкости в напорных трубопроводах были заложены в работе Н.Е. Жуковского [3] и в лекциях по гидродинамике, читавшихся им в Московском университете.
В настоящее время предложено и используется много различных методов борьбы с вредными пульсациями давления. Например, снижение уровня пульсаций давления в трубопроводных системах энергетических установок на практике достигается путем изменения конструкции трубопроводов, уменьшения возмущающих сил за счет уравновешивания ротора и компрессора, увеличения толщи-
ны стенок, виброизоляции трубопроводов от источников возбуждения, рациональной укладки трубопроводов и опор. Как правило, такие способы борьбы с вредными пульсациями давления оказываются весьма дорогостоящими. В связи с этим возникла необходимость создания специальных устройств — гасителей пульсаций давления, которые при включении в систему трубопровода существенно снижали бы в нем уровень вредных пульсаций давления. Различным типам, конструкциям и методам расчета стабилизаторов давления в трубопроводах также посвящена весьма обширная литература.
Одним из распространенных типов гасителей пульсаций давления являются различные резонансные звукопоглощающие системы. Изучением таких систем успешно занимался С.Н. Ржевкин [4]. Резонатор представляет собой замкнутую полость с жесткими стенками, сообщающуюся через узкий канал с трубопроводом, в котором необходимо устранить вредные пульсации давления. В настоящее время предложено большое число различных конструкций гасителей указанного типа, общие принципы управления такими системами изложены в монографии А.Г. Бутковского «Методы управления системами с распределенными параметрами» [5].
Статистические данные последних лет эксплуатации магистральных паро-и газопроводов показывают, что применяемые методы защиты от вредных пульсаций давления являются недостаточно эффективными. Поэтому задача разработки стабилизаторов давления новой конструкции является весьма актуальной.
Авторами разработана конструкция более эффективных стабилизаторов давления, изготовлен опытный образец и проведены экспериментальные его исследования на универсальном комплексе. Полученные результаты показали высокую эффективность работы новых стабилизаторов давления. Однако с целью обоснования параметров нового стабилизатора давления необходимо выполнить теоретические исследования и сравнить их с результатами эксперимента.
Принципиальная схема нового стабилизатора давления (опытного образца) приведена на рис. 1.
Рис. 1. Расчетная схема стабилизатора диссипативного типа с одним слоем перфорации: 1 — труба трубопровода; 2 — камера стабилизатора
Для математического описания движения пара (газа) через стабилизатор используются уравнения акустики, написанные в цилиндрической системе координат. Рассматривается случай короткого трубопровода. В этом случае перепад давления на входе и выходе из трубопровода намного меньше пульсаций давления. Следовательно, сопротивлением движению пара (газа) в трубопроводе можно пренебречь. Влияние силы тяжести на движение пара (газа) также не учитывается. При расчетах предполагается, что движение пара (газа) происходит при постоянной температуре, т.е. рассматривается изэнтропическое движение идеального газа. В этом случае идеальный газ можно рассматривать как ба-ротропную среду. Цилиндрическая система координат г, ф, 2 выбирается таким образом, чтобы ось 2 совпадала с продольной осью трубы паро-, газопровода. Уравнения движения и неразрывности такой среды при отсутствии внешних массовых сил имеют вид [6]
2
диг диг иф диг диг иф 1 др + иг I + и2 — ,
дt дг г дф дг г р дг
диф диф иФ диф диф ифиг — 1 др
1 ■ и„--1 + и 1
dt dr r 5ф dz r p • r 5ф duz duz иф duz duz _ 1 dp
(1)
dt r dr r дф z dz p dz dp + d(p • ur • r) + + r d(p^ uz) =
dt dr дф dz
Стационарное движение газа (в дальнейшем значок 0) в одномерной системе газопровод—стабилизатор без пульсаций давления описывается уравнениями движения Эйлера [6]:
duz 1 йр й (ри2)
z dz p dz dz Из (2) получаем: puz _ q _ const,
_ 0. (2)
duz dp d r 2 .... pu^—- + [pu2 + p] _ 0.
dz dz dz
т.е. pu;? + p _ const.
На входе в трубопровод имеем:
uz _ u° _ uBX _ const; p_p0 _ q / uBX _ const; p _ p0 _ p(p0) _ const.
Следовательно,
pu2 + p _pbxuB2x + Pвх, puz _ pbXuBX _ q,
p _ p(p0).
Полученная система имеет решение:
0 0 0 ы2 = ы2, р = р , р = р .
Между перфорированной трубой и корпусом стабилизатора пар (газ) при отсутствии пульсаций давления неподвижен. Следовательно,
и2 = 0, р = р0, р = р0.
Нестационарное движение (значок ') пара (газа) в трубе паро- и газопровода описывается уравнениями акустики [7], полученными из уравнений (1), при этом возмущения всех величин считаются малыми:
0 ди'г ди'г др' 0 диф диф др' р0—L + q—- = ——, рг—- + qr—- = ——, дt дг дг дt дг дф
р0 К + qK = _&, р'= (р0)р',
дt д2 д2 й р
др ' 0 д(и>) 0 диф 0 ди' 0 др ' Л г-*- + р0 4 г ' + р0—^ + р0г-^ + ги°° — = 0. дt дг дф дг дг
(3)
(4)
В стабилизаторе имеем:
0 сК = _р 0 ди'ф = др'
дt дг ' дt дф'
р0 ди^ = _<р ,= йр^(р0)р, дt дг йр
г др;+р0 Мт»+р0 Ч+р0 г =0.
дt дг дф дг
В трубе газопровода иг = иг, иф = иф, иг = и° + иг, р = р0 +р ' , р = р0 + р',
в стабилизаторе иг = иг, иф= иф, иг = и'г, р = р0 +р ' , р = р0 + р'.
Осредним уравнения (3) и (4) по поперечному сечению газопровода и поперечному сечению стабилизатора по формулам
1 2п г0
— II г/(г, ф, г, t)йгйф =/(г, t) (труба), (5)
^тр 0 0
1 2п г1
-II г/(г, ф, г, t)йгйф = /ст (г, t) (стабилизатор), (6)
^ст 0 г0
где 5'тр — площадь поперечного сечения трубы газопровода, 5'1 — площадь поперечного сечения стабилизатора (см. рис. 1).
Для газопровода (в дальнейшем значок ' опускается) получим:
0 дй7 дй7 др 2„ р0^Т + = *, р = с р,
дt дг дг
^ + р0^ + и^ = 0. (7)
дt дг дг
где с2 — йр/ йр — квадрат скорости звука в газе; V — 5трЬ и /тс (г, t) — поток массы через боковую поверхность трубы газопровода в сечении г, (кг/сек).
В поперечных сечениях со стабилизатором давления третье уравнение в (7) будет иметь вид
др 0 дй 0 др 1 ч
—+ р° — + и0 — —--/тс (г, t), (7а)
& дг 2 дг 2л^р
где V — £Ь и 1тс (г, t) — поток массы через боковую поверхность трубы газопровода в сечении г, (кг/сек):
2п
4 (z, 0 — 2пЬР0 | [Ч ]г —г0 йф.
0
Положим
ртр (г, t) — р(г, t) + диг (г, t), (8)
где р1р (г, ^ — давление в трубе без учета давления д^ (г, t) от переносного ускорения стационарного члена д — и°р° нестационарной скоростью Й2 (г, t).
Подставляя (8) в (7) и учитывая, что р — -у р, получим:
с
дрТр , 0 0 дртр , 0Г 2 , 0x2-, диг
+ р [с - (иг ) — 0,
дt дг дг
р0
ди дртр
(9)
дt дг
Исключая из уравнений (9) Й2 (г, t), получим:
д2 ртр д2 р1р д2 ртр
^др- + * ^ - [с2 - (и0)2]^ — 0. (10)
Выполняя аналогичные выкладки для уравнения (7 а), получим:
д2 ртр д2 ртр д2 п1? с2 дР
+ 2и0 ^ - [с2 - (и°0)2 ]^ — --. (10а)
дt2 г д^с^ 1 г дг2 2п^гр й
Возмущения скорости, плотности и давления в стабилизаторе будем обозначать соответственно и2 (г, t), рст, рст (г, t).
Для стабилизатора с постоянным объемом Уст приток (отток) газа массы Аш сопровождается увеличением (уменьшением) осредненной плотности газа в стабилизаторе на величину Арст; Аш — Арст Уст. Для рассматриваемой баротропной среды
Арст
Арст — ■
2
с
где Арст (г, t) — осредненное приращение давления в стабилизаторе.
Следовательно,
v ст
Арс
Ат =АР_ Vt . (11)
с
Приток газа массы т равен
т = J /тс (z, t)dt. (12)
Подставляя (12) в (11), получим
2
рсТ (Z, t) = V- J /тс (Z, t)dt. (13)
Перетекание газа через отверстия. Будем считать, что длина стабилизатора
Ь << X,
где X — длина волны основного тона пульсаций массового расхода.
Поэтому давление по длине стабилизатора в момент времени ^ считается постоянным. Систему перфорированных отверстий, равномерно распределенных по длине стабилизатора заменим одним отверстием в сечении 2 = 2П. Общий поток газа через это отверстие будет равен
Ь
1п (0 = { /Тс (2, 1)й2. (14)
0
Проведенные к настоящему моменту исследования показали, что зависимость скорости перетекания упер газа через отверстие от разности давлений в магистрали и стабилизаторе можно аппроксимировать формулой (8):
рТр(Гц, 0 - Р0"(211, 0 = 4^пер + В^ ^пер|, (15)
где А1 и В1 — постоянные коэффициенты, полученные экспериментально.
Уравнение (15) можно линеаризовать, полагая В1 = 0 при упер << Vтр.
Обозначим
4
A = ■
р • F
гтс тр
где ртс — плотность газа, перетекающего из трубы в стабилизатор, ртс « р0. Тогда формула (15) примет вид
р^ (Z11, t) - рст (Z11, t) = A • In (t), (16)
^11 (t) = Упер Ф- Fотв ,
где FоТв — площадь отверстия в сечении zn, равная суммарной площади перфорированных отверстий. Коэффициент A характеризует зависимость проводимости отверстий от формы и геометрических размеров.
Движение газа в трубопроводе. Разобьем магистраль трубопровода на два участка:
I участок (г = 1): г е [-Х1; г11],
II участок (г = 2): г е [г11; Ь2].
Будем искать решение уравнения (10) для г-го участка трубопровода в виде ряда Фурье:
р,т(г, %) = X<(г)со8К%) + Ь?(г)ып(шк%), (17)
к=1
где шк = 2пк/Т, Т — главный период пульсаций давления.
Подставляя (17) в уравнение (10) (коэффициент г опускается), получим: "с2 - (М0)21- 2шк„о^ + ш2акр(г) = 0,
(18)
й 1атк - 2шки° йЪ?
йг 2 йг
й 2Ъ? + 2шки0 йа1^
йг1 йг
с2 -(и0)2 ^ + 2шА0^ + ш2кЪ?(1) = 0.
В дальнейшем будем рассматривать случай м° / с = 8 << 1 и пренебрежем чле-
йЪ^ йа^к _ _
нами, содержащими —— и ——. Отбрасывая указанные члены и члены поряд-
йг йг
2
ка 8 , получим следующую систему уравнений:
й 2атр
^От + 14 акР( г) = 0,
2 (19)
^ + ЪД г) = 0,
йг
шк 2пк
где ^к = — =
с сТ
Система (19) распадается на два отдельных одинаковых (с точностью до обозначений) уравнения для функций акр( г) и ЪД г). Решение системы имеет вид
ак = С1 С08(Цкг) + С2к вт^), (20)
3 4 (20)
Ък = Ск С0й(Цкг) + С4 81П0кг). Тогда давление в магистрали для г-го участка примет вид
да .
рк (г%) = Е1Гс1 со^кг) + С181п(Цкг)] соэ(шк%) +
к=Л (21)
Сгк
Из второго соотношения (9) имеем:
Р тр д% тр дг
Массовый расход газа
G - =р0 /ТД =-^Тр J^ dt. (22)
Подставляя (21) в (22), получим
F <»
Gk(z, t) = £{[Ck kZ) - C3 sin(^z) 1 cos(mkt) +
c k=1lL 1 (23)
\_-Cl cos(^kz) + Cik sln(^kz)] sln(mkt)}.
Краевые условия
1. Пусть в некотором сечении газопровода z = - L1 заданы пульсации массового расхода
G (-L1, t) = Go (t). (24)
Разложим Go(t) в ряд Фурье:
Go (t) = Z gk ( z ) cos( mkt) + hk ( z) sin( mktX (25)
k=1
2. При z = L2 примем давление в магистрали равным нулю, т.е.
ртр (L, t) = 0. (26)
3. Условие равенства массовых расходов газа в узле z = z11 :
G1P (zn, t) + In (t) - Gf (zn, t) = 0. (27)
4. Давление при z = z11 - 0 и z = z11 + 0 должны быть равны, т.е.
ртр (zn, t) = р2? (zn, t). (28)
5. Условие (16) перетекания газа через перфорационное отверстие (через систему отверстий) при z = z11:
ртр (zn, t) - рст (zn, t) = A • I11 (t). (29)
Исходные уравнения. Массовый расход газа через отверстие в сечении z = z11 представим в виде ряда Фурье:
I11(t) = Zak cos(mk0 + Pk sin(mkt). (30)
k=1
Подставляя (14) в (13) и используя (30), получим
c2 г ~ c2 ( а В Л
рст (zn,t) = —J/11 (t)dt = —Z — sin(mkt) cos(mkt) . (31)
V J v . . mt mt
ст к=1V " 1к "1к У
Подставляя выражения для рг1р (2, ?), рс (211, ?), О^ (2, ?), /11 (I) в краевые условия (24), (26)—(29) при фиксированном к (к в уравнениях опускается) и группи-
руя выражения при ео8() и 8т(), получим систему из десяти уравнений относительно С, С2, С3, с4, С, С22, С23, С24, а, в:
-С- ( sin(L + Cj4 cos(^Xj)) = g,
Ftl (-CÍ sin(|jL2) - C22 cos(|oL2)) = A,
C22 cos(^L2) + C^ sin(^L2) = 0, C23 cos(^L2) + C4 sin(^L2) = 0, F
(-C23 sin(|jz22) + Cj4 cos(^z22) + C23 sin(^z22) - C^ cos(^z22) ) + а = 0,
—тр( sin(|jz22) + CÍ2 cos(^z22) -C2 sin(^z22) + C22 cos(^z22)) + в = 0,
Cj cos(^z22) + C22 sin(^z22) - C2 cos(^z22) - C2 sin(^z22) = 0, C23 cos(^z22) + Cj4 sin(^z22) - C 23 cos(^z22) - C4 sin(^z22) = 0,
(32)
Cj2 cos(^z2 2) + C22 sin(^z2 2) -
mV,'
C23 cos(^z2 2) + C24 sin(^z22) -
c2a
mV
= Aa,
= Ap.
Решая систему (32) для каждого к и суммируя решения для n членов разложения G0(t), получим выражения для pjp (z, t), рст (z22, t), GJP (z, t), I22 (t) согласно формулам (22), (32), (23), (30). Алгоритм реализован в среде MAPLE (GasTigrV). Пример
Исходные данные:
L2 = 20 м, L = 5 м, L2 = 25 м, Гтр = 0,7854 м2 (диаметр трубы 2 м), FCT = 39,27 м2. Массовый расход G0(t) = 200 • sin (rat) кг/с, где ra = 4п 2/c (T = 0,5 c). c = 500 м/c. Коэффициент проводимости отверстий A = 30 /мс.
На рис. 2 изображены графики давления в сечении z = 7,5 м без стабилизатора (кривая 2) и со стабилизатором (кривая 2).
Рис. 2. Давление в магистрали, МПа
Выводы. Компактных размеров стабилизатора давления можно добиться при относительно малых периодах пульсаций давления (Т < 0,1 с). В противном случае необходимо устанавливать упругие камеры.
При увеличении частоты массового расхода (и (или) уменьшении скорости звука) при расчетах возникают ограничения по длине рассматриваемой магистрали. Для более протяженных систем необходимо введение дополнительных участков магистрали.
ЛИТЕРАТУРА
[1] Тимошенко С.П. Колебания в инженерном деле. — М.: Наука, 1967.
[2] Рэлей Д.В. Теория звука: в 2 т. — М.: Гостехиздат, 1955.
[3] Жуковский Н.Е. О гидравлическом ударе в водопроводных трубах. — М.: Л.; Гостех-издат, 1949.
[4] Ржевкин С.Н. Курс лекций по теории звука. — М.: Изд-во МГУ, 1960.
[5] Бутковский А.Г. Методы управления системами с распределенными параметрами. — М.: Наука, 1975.
[6] Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. — Ч. 1. — М.: Физмат-гиз, 1963.
[8] Мясников М.П., Низамова Г.Х. Расчет стабилизатора давления диссипативного типа // Отчет № 4286 Института механики МГУ, 1993.
[9] Рахматулин Х.А. Обтекание проницаемого тела // Вестник МГУ. Серия физико-математических и естественных наук. — 1950. — С. 41—45.
REDUCTION OF PRESSURE OSCILLATIONS IN MAGISRAL GAS TUBES
E.K. Sinichenko, F.V. Rekach
Peoples' Friendship University of Russia Ordzhonikidze, 3, Moscow, Russia, 115419
The absorber pressure, the principle of action is based on the elastic-damping effect on the flow of gas have been considered.
Key words: pressure oscillations, movement of gas current.