Прогнозирование энергоэффективности УЭЦН
Аннотация
электроцентробежный насос полезный трансформатор
Авторы представляют опыт разработки методики проведения расчетов КПД элементов типичных УЭЦН и энергоэффективности УЭЦН в целом.
Annotation
The authors present development experience of efficiency computing methods of elements of typical electric-centrifugal pumps' rigs and their energy efficiency as a whole.
Одним из критериев эффективности добычи нефти установками электроцентробежных насосов (УЭЦН) является величина удельного расхода энергии, например, на тонну добытой продукции. Поскольку добыча пластовой жидкости заключается в передаче ей энергии для транспортировки из пласта на поверхность, то энергоэффективностью добычи является отношение энергии необходимой для подъема жидкости к затраченной энергии, т.е. КПД УЭЦН.
В основном энергия расходуется на работу насоса (ЭЦН) и погружного электродвигателя (ПЭД), хотя потери энергии в предвключенных устройствах, а также в кабеле, трансформаторе и др., также значимы и их необходимо учитывать в общем балансе энергии. Кроме того, имеет место взаимное влияние элементов УЭЦН. Так при увеличении потребления энергии, например, ЭЦН - растет потребление ПЭД, потери в кабеле, трансформаторе и др.
Современная УЭЦН представляет собой линейную цепь последовательно соединенных элементов: станции управления, повышающего трансформатора, кабельной линии, ПЭД, предвключенных устройств (газосепаратора, диспергатора, мультифазного насоса, сепаратора мехпримесей и др.) и ЭЦН. Обозначим через N0 мощность, поступающую на вход первого элемента УЭЦН (станции управления), а через N1 мощность, передаваемую от первого элемента - второму. Через Nn обозначим мощность на выходе с последнего элемента (ЭЦН). Тогда выражение для КПД УЭЦН можно записать следующим образом:
з=NnN0=N1N2?N2N3?Nn?1Nn?NnN0=з1?з2?зn?1?зn, (1)
где зi - КПД i-го элемента УЭЦН. Из (1) следует, что по критерию энергоэффективности добычи, лучшей будет УЭЦН в которой максимально произведение КПД отдельные ее элементов. Поскольку КПД элементов, см. выше, могут быть взаимосвязаны, нужна методика расчета, учитывающая как условия работы УЭЦН, так и взаимосвязь КПД элементов. Разработка такой методики и проведение ряда расчетов энергоэффективности типичных УЭЦН, стали целью данной работы.
Методика расчета КПД элементов УЭЦН
КПД ЭЦН и предвключенных устройств
Мощность, передаваемая от насоса к жидкости Nout=сgHQ+ДP?Q+ДNTP, где ДP - заданное избыточное давление, создаваемое на устье скважины (его величина порядка 20 атм), ДNTP - потери на трение в НКТ. Поскольку обычно сgHQ?ДNTP, то КПД ЭЦН вычисляется следующим образом
зESP=сgH(Q)?Q+ДP?QNin(Q),(2)
где Nin - мощность, потребляемая насосом, Q - подача, H - напор, создаваемый насосом, с - плотность жидкости, g - ускорение свободного падения.
Напор и мощность, потребляемые насосом, являются функциями Q, n и свойств рабочей жидкости. Обычно в каталогах фирм, производителей ЭЦН приводятся зависимости H и Nin от подачи Q, полученные в испытаниях на воде при n=3000 или 3500 об/мин, и от частоты вращения n, вычисленные из данных этих испытаний в предположении, что Q ~n, H ~n2, Nin ~n3.
В [1] показано, что использование зависимости H ~n2 приводит к относительной ошибке ±1 - 2%, а зависимости Nin ~n3 к ошибке ±5 - 7%.
Влияние вязкости скважинной жидкости обычно учитывают по корреляционным зависимостям [2 - 6]. Согласно [7] это приводит к относительной ошибке ±10 - 25%.
Потери энергии в предвключенных устройствах обычно не превышают 1 - 2 кВт, что много меньше типичного значения мощности, потребляемой ЭЦН. Поэтому, обычно эти потери можно не учитывать. Ели же требуется высокая точность прогноза энергоэффективности УЭЦН, необходимо знать величину этих потерь, которые в настоящее время приводят далеко не все компании, их изготавливающие.
КПД погружных электродвигателей
По определению, КПД равен отношению механической мощности на валу ПЭД к потребляемой из электрической сети:
зM=N2MN2M+ДNLOST(N2M),(3)
где N2M - мощность на валу двигателя, NLOST - суммарные потери в двигателе. Обычно по данным, приводимым производителями ПЭД, можно установить зависимость NLOST от мощности, потребляемой из сети, а значит, и от N2M.
При расчете необходимо учитывать, что мощность на валу двигателя N2M равна суммарной мощности, потребляемой всеми устройствами, приводимыми в движение: гидрозащитой, ЭЦН, предвключенными устройствами (газосепаратором, диспергатором, мультифазным насосом, сепаратором механических примесей и др.).
КПД станции управления
В настоящее время КПД станции управления определяется производителями как постоянная величина, указанная для номинального режима работы станции. Поскольку режим работы станции управления выбирается таким образом, чтобы не превышать максимальное значение тока при любых допустимых значениях загрузки двигателя, то потери энергии в станции управления будут ограничены и не превысят справочного значения.
КПД трансформатора
Потери в трансформаторе характеризует КПД трансформатора, равный отношению мощности, передаваемой из трансформатора в кабельную линию N1K к мощности, поступающей на его первичную обмотку N1K+ДNT:
зT=N1KN1K+ДNT,(4)
где ДNT - мощность, теряемая в трансформаторе. Если выразить ДNT через параметры трансформатора, см. Приложение 1, получим:
зT=N1KN1K+NXX(UHM+ДUK(N1M)U2H)2+NK3(IM(N1M)I1HU2HU1H)2,(5)
где N1K - мощность потребляемая на входе кабеля, NXX - потери холостого хода трансформатора, N3 - потери трансформатора в опыте короткого замыкания, UHM - номинальное напряжение двигателя, ДUK - падение напряжения на кабеле (зависит от потребляемой двигателем мощности - см. Приложение 2), U2H - напряжение отпайки, IM - ток потребляемый двигателем (зависимость от потребляемой мощности обычно приводится производителем), I1H, U1H - номинальные ток и напряжение в первичной обмотке, U2H - напряжение отпайки.
Поскольку в основе этой формулы лежит эмпирическая зависимость, то ошибка в расчетах будет определяться точностью измерительных приборов и точностью аппроксимации экспериментальной зависимости, и будет составлять не более 1% [8].
Потери в кабеле
При определении КПД кабеля полезной является мощность, поступающая на вход ПЭД, т.е. N1M=N2M+ДNLOST, см. (3). Обозначив потери в кабеле через ДNK, получим следующее выражение для КПД кабеля:
зK=N1MN1M+ДNK.(6)
Если выразить потери в кабеле ДNK через его характеристики (см. Приложение 2), то
зK=N1MN1M+ncI2M(N1M)?с?LKSK(1+бT(TK?T0)),(7)
где с - удельное сопротивление меди при температуре T0=20°С, составляет 0,0175 Ом·мм2/м, LK - длина кабеля, м, SK - площадь сечения жилы кабельной линии, мм2, бT - температурный коэффициент сопротивления меди, равный 0,0041/°С, TK - абсолютное значение температуры жилы, °С, nc- число жил кабеля, обычно равняется трем, IM - ток потребляемый двигателем (зависимость от потребляемой мощности обычно приводится производителем).
Поскольку в основе этой формулы лежит температурная зависимость сопротивления меди, то в диапазоне температур от -50°С до +200°С эта зависимость почти линейная и погрешность метода не превышает 1% [8].
Примеры расчетов КПД УЭЦН различной комплектации
Будем считать заданными напор и подачу насоса. Зафиксируем частоту вращения вала (n=2910 об/мин) и будем считать, что рабочая жидкость - вода. Пусть кабель имеет длину 950 м и применяется трансформатор ТМПН-250 фирмы BIRZST, отпайка 2168 В, потери холостого хода 0,58 кВт, потери короткого замыкания 3,9 кВт.
Проведем расчет энергоэффективности, т.е. вычислим КПД, для УЭЦН различной комплектации. Используются ЭЦН? имеющие одинаковую номинальную подачу и суммарный напор, но разный КПД: ВНН5А-100 (кол-во ступеней 125) и ВНН5А-100Э (кол-во ступеней 228), асинхронный ПЭДН63-117 и вентильный ПВЭДН63-117. Все компоненты производства ЗАО «Новомет-Пермь». Также будем использовать некоторое гипотетическое предвключенное устройство (ПУ), потребляющее 1 кВт, чтобы оценить влияние таких устройств на общий КПД УЭЦН. В расчете используется КПД станция управления СУ-160 имеющая КПД, равный 98%.
Табл. Результаты расчета КПД при различной комплектации УЭЦН
Из результатов расчетов, приведенных в табл., видно, что КПД установки с вентильным двигателем выше на 31%, чем УЭЦН с асинхронным ПЭД, если используется обычный насос. В случае применения энергоэффективного насоса КПД установки при использовании вентильного двигателя выше на величину 32% по сравнению с асинхронным. Использование вентильного двигателя в совокупности с энергоэффективным насосом дает улучшение энергоэффективности на 36% по сравнению с обычной комплектацией (асинхронный двигатель + обычный насос).
Кроме того, видно, что в некоторых случаях (например, АД + ВНН з = 0,39 и АД+ ВНН+ПУз = 0,4) использование предвключенных устройств дает увеличение значения общего КПД. Это объясняется следующим: предвключенное устройство не забирает энергию из потока передаваемой жидкости, а увеличивает нагрузку на валу двигателя. А поскольку для выбранного оборудования рабочая точка двигателя была выбрана левее максимума соответствующей кривой, то дополнительная нагрузка сдвинула рабочую точку вправо, что привело к увеличению КПД двигателя и увеличению КПД УЭЦН. Этот пример показывает важность оценки КПД в комплексе, с учетом влияния отдельных компонентов друг на друга.
Выводы
1. Поскольку добыча пластовой жидкости заключается в передаче ей энергии для транспортировки из пласта на поверхность, то энергоэффективностью добычи естественно считать отношение энергии, необходимой для подъема жидкости, к затраченной энергии, т. е. КПД УЭЦН.
2. Предложены формулы для расчета КПД отдельных элементов УЭЦН: насоса, предвключенных устройств, ПЭД, кабеля, трансформатора, станции управления. Оценена точность расчета, определяемая не столько приближениями, сделанными при выводе формул, сколько отсутствием данных о характеристиках элементов УЭЦН в технической документации предприятий - изготовителей этих элементов.
3. Приведены примеры расчета КПД УЭЦН различной комплектации, позволяющие увидеть вклад отдельных элементов УЭЦН в энергоэффективность установки в целом.
Литература
1. Нефтяное хозяйство. 2013. №6.
2. Hydraulic Institute, Standard for Effects of Liquid Viscosity on Rotodynamic (Centrifugal and Vertical) Pump Performance, 1983, USA.
3. Gulich J.F. Centrifugal pumps. Second Edition, Springer, 2010.
4. Li V.G. Experimental Research of Technical Characteristic centrifugal pumps // World Pumps. 2002. №26.
5. Gilmar Amaral, Valdir Estevam, Petroleo Brasileiro and Fernando A. Franco, Influence of Viscosity on ESP Performance, 2007 SPE Annual Technical Coonference and Exhibition, Anaheim, California, 11 - 14 November.
6. Ляпков П.Д. О влиянии вязкости на характеристику погружных центробежных насосов. Труды ВНИИ, вып. 41. М.: Недра, 1964.
7. Абахри С.Д., Перельман М.О., Пещеренко С.Н., Рабинович А.И. Влияние вязкости на рабочие характеристики центробежных насосов // Бурение и нефть. 2012. №3. С. 22 - 26.
8. Методы и средства измерения тепловых параметров. [Электронный ресурс] URL: http://xeon128.narod.ru/lection/6.htm (дата обращения 1.08. 2013).
9. Теоретические основы электротехники. [Электронный ресурс] URL: http://toe-kgeu.ru/toe1/259-toe1 (дата обращения 2.08.2013).
10. Писаревский Ю.В., Ген Ж.А. Проблемы разработки высокоскоростных электрических двигателей для стоматологии // Электромеханические преобразователи. 2011. №4. С. 46 - 51.
11. Забудский Е.И. Электрические машины. Ч.1. Трансформаторы. Учебное пособие для вузов. М.: МГАУ, 2002. 167 с.
12. ТНК ВР. Методика оценки мероприятий, направленных на экономию электроэнергии и выбору УЭЦН по критерию энергоэффективности.
13. Сипайлов В.А. Оптимизация режимов работы установок электроцентробежных насосов механизированной добычи нефти: дисс. на соискание ученой степени к.т.н. Томск: ТПУ, 2009. 196 с.
Приложение 1. Расчет потерь энергии в трансформаторе
Величина ЭДС, индуцируемой в первичной обмотке, определяется скоростью изменения магнитного потока. Результирующий магнитный поток запишем в виде ЦC(t)=Цmcosщt=Цmcos2рft, где Цm - амплитуда этого потока, щ - циклическая частота, f - частота. Тогда получим:
e1(t)=W1dЦC(t)dt=W1ddtЦmcos2рft=2рfW1Цmsin2рft=E1msinщt,(П1)
где W1 - число витков в первичной обмотке, E1m=2рfW1Цm - амплитуда ЭДС, индуцируемая в первичной обмотке, тогда действующее значение ЭДС
E1=E1m2v=2рfW1Цm2v.(П2)
Магнитный поток Цm можно выразить через индукцию Bm
Цm=BmSC,(П3)
где SC - площадь поперечного сечения сердечника.
Тогда величина ЭДС первичной обмотки
E1=2рfW12vBmSC.(П4)