CC BY
151
УДК 621.65:681.586 ББК 33.361.3:32.96-044.3
А. И. Исмагулова, Н. Г. Романенко
МОДЕЛИРОВАНИЕ ПРОЦЕССА ИЗМЕРЕНИЯ ДИНАМИЧЕСКОГО УРОВНЯ ЖИДКОСТИ В ШТАНГОВОЙ ГЛУБИННО-НАСОСНОЙ УСТАНОВКЕ НА ОСНОВЕ СРЕДСТВ ДИНАМОМЕТРИРОВАНИЯ
A. I. Ismagulova, N. G. Romanenko
MODELING OF THE PROCESS OF MEASURING THE DYNAMIC FLUID LEVEL IN ROD DEEP-WELL PUMPING UNIT ON THE BASIS OF DYNANOMETRY DEVICES
Приводится математическая модель построения динамограммы уровня жидкости в штанговой глубинно-насосной установке. Полученная динамограмма используется для регулирования уровня жидкости в скважине с целью повышения эффективности и надежности ее работы. Полученные математические соотношения формализованы в виде алгоритмической процедуры расчета характеристик динамограммы. Приведен численный пример, обосновывающий работоспособность полученных соотношений и алгоритмов.
Ключевые слова: штанговые глубинно-насосные установки, динамограмма, уровень жидкости в скважине, математическая модель.
The mathematical model of construction of dynamometer card of fluid level in rod deep-well pumping unit is presented. The obtained dynamometer card is used to control the fluid level in the well in order to increase the efficiency and reliability of its work. The resulting mathematical relationships are formalized in an algorithmic procedure for calculation of the characteristics of dynamogram.
A numerical example substantiating the performance of these correlations and algorithms is given.
Key words: rod deep-well pumping units, dynamometer card, fluid level in the well, mathematical model.
Введение
Надежность и производительность штанговой глубинно-насосной установки (ШГНУ) в значительной степени определяется эффективностью работы системы автоматического регулирования динамического уровня жидкости в скважине, что, в свою очередь, целиком и полностью зависит от качества измерения текущих значений требуемых параметров установки. В отсутствие стационарных средств прямого измерения динамического уровня, устанавливаемых на скважине, для непрерывного контроля уровня жидкости в скважине используют косвенные методы измерения, основанные на информации, доступной для контроля средствами системы управления ШГНУ
- динамограммах или ваттметрограммах. При оснащении ШГНУ стационарными средствами ди-намометрирования датчиками усилий и хода устьевого штока, которые позволяют осуществлять непрерывный контроль за работой установки, в качестве исходной информации для оценки динамического уровня вполне естественно выбрать данные динамометрирования. Основываясь на модели усилий в подвеске устьевого штока, проведем исследование зависимости составляющих усилий от характеристик скважины и значения динамического уровня.
Математическая модель изменения уровня жидкости в ШНГУ
В практике эксплуатации ШГНУ для оценки состояния и качества работы оборудования широко используется динамометрирование установок [1, 2], т. е. снятие диаграмм нагрузки на устьевой шток в зависимости от его хода. Процедура практического использования динамограмм работы штангового насоса состоит в их сравнении с так называемыми «теоретическими» динамограммами, которые находятся заранее на основе аналитических методов и обычно учитывают только статические нагрузки и эффект растяжения колонны штанг и насосно-компрес-
сорной трубы (НКТ) под действием периодически перекладываемого веса жидкости. Теоретические динамограммы могут быть использованы также и при оценке и верификации результатов моделирования процессов в ШГНУ. Сходные вопросы рассмотрены в [3].
Приведем основные соотношения для расчета параметров «теоретической» динамограммы, используя известные положения [1, 4]. При работе ШГНУ штанги постоянно находятся в жидкой среде, поэтому для вертикально установленной скважины статические нагрузки в точке подвеса колонны штанг обусловлены весом штанг в жидкости —Ш ж и весом поднимаемого
столба жидкости —Ж . При ходе вверх нагнетательный клапан закрыт, поэтому статическая нагрузка на штанги равна —..ф) = —Ш ж + —Ж ; при ходе вниз нагнетательный клапан открывается, нагрузка —Ж снимается со штанг и передается на трубы, т. к. связанный с ними всасывающий клапан закрыт, и поэтому статическая нагрузка на штанги при ходе вниз равна -^(в) = —Ш ж . Вес штанг в жидкости с учетом действия архимедовой выталкивающей силы
—Ш.Ж —Ш ^Ш ■ (Р Т Р 0) —ш
1 -
^Ш ■ (Р Т Р 0)
-Ш
-ШЖ • Ь , (1)
Ш.Ж арх’ 4 ;
где —Ш - вес штанг в воздухе; 5”Ш • (р Т-р 0) - архимедова сила; 5”Ш - площадь сечения штанг; р т - давление жидкости в трубах над плунжером, действующее на нижний торец штанг; р{) -атмосферное давление, действующее на верхний торец штанг; Ьарх - коэффициент, учитывающий архимедову силу или потерю веса штанг в жидкости (коэффициент плавучести штанг):
Ь = 1 - ^Ш •(рт-ро)
арх — '
Ш
Давление в трубах р т определяется суммой гидростатического давления столба жидкости в трубах р 1, потерь давления на трение жидкости в трубах, устьевого давления (на выкиде
скважины) за вычетом давления разгрузки в результате газлифтного эффекта (выполнения подъемной работы энергией расширения выделяющегося из нефти газа). Поскольку гидростатическое давление столба жидкости р 1 значительно больше суммы всех остальных составляющих давления р т , то, пренебрегая также величиной р0, получаем
Ь = 1 _ (р т-р р) = 1 _ ' ¿ырж8 = 1 - Рж (2)
арх -Ш ¿ыРс^ш Рс ’
где р 1= ^Ш • ¿ыРж8; Ьы - глубина спуска насоса (длина колонны штанг); рЖ - плотность
жидкости; 8 - ускорение свободного падения; рс - плотность материала штанг (стали).
Вес штанг в воздухе —Ш = ¿ыРс; однако, с учетом секционирования колонны штанг, ее вес будет определяться по общему выражению
N
-Ш = Е ^ш/ш,'Рс8 , (3)
1=1
где ^ш г- и ¿ш ; - площадь сечения и длина 1-й ступени штанг.
Нагрузка —Ж обусловлена разницей давления жидкости над и под плунжером насоса, т. е.
РЖ = ^ (РТ - Р ), (4)
Ж Т г вс.ц' ’ ' '
где - площадь поперечного сечения плунжера (цилиндра) насоса; рвс ц - давление под плунжером (всасывания цилиндра):
Рвс.ц = Рпр - ДРКЛ.Б . (5)
Здесь Рпр - давление на приеме насоса:
Р = И • р 2 + р , (6)
г пр г затрд г затр ’ 4 '
И - глубина погружения насоса под динамический уровень; рз - средняя плотность жидкости в затрубном пространстве; Рз - давление газа в затрубном пространстве на уровне жидкости; АРКЛ б - перепад давления во всасывающем клапане (местное сопротивление).
На основе вышесказанного можно записать следующее выражение для приближенного нахождения нагрузки РЖ (Р Т@ Р1; аРкл б @ о; Рзатр = 0; рзатр = рЖ):
Рж @ Бп •(ХиРж2 - ИРж2) = Бп •(¿и — И)рж2 . (7)
Статическая нагрузка на НКТ, исходя из сказанного выше, будет иметь следующие значения при ходе штока вверх и вниз соответственно:
Р ТУИЧ = и Р т,ц, = -Р -щ- + -Рт- , (8)
ст.Т(Н) Т.Ж ст.Т(Н) Т.Ж Ж ’ ' '
где Рр ж - вес НКТ с учетом частичного (под динамический уровень) погружения труб в жидкость:
ИД + (ЬН — ИД)
РТЖ = РТ Д Н—^ь ; (9)
Т.Ж Т т арх ' '
ЬН
Рт - вес колонны НКТ в воздухе:
Рт = ^т Ья р с 2 . (10)
Здесь - площадь сечения металла труб, определяемая по выражению
где Р>Т и 5т - наружный диаметр и толщина стенки НКТ.
Таким образом, выражения (1)-( 11) позволяют получить значения статических нагрузок в точке подвеса колонны штанг и НКТ при ходе штока вверх и вниз. Гидростатическая нагрузка
—ж , попеременно действуя на штанги или трубы, вызывает их упругие деформации 1ш и 1т
в соответствии с законом Гука. С учетом секционирования колонны штанг имеем
— NC I
1Ш = —Ж^Ы / (Еу^Ш) = Е Е ^ и ^Т = —ж^ы / (Еу^Т), (12)
у 1 =1 Ш.1
где Еу - модуль упругости (Юнга); Ш г- и $ш г- - длина и сечение штанг 1-й секции.
В результате действия нагрузки Рж перемещение плунжера вверх относительно цилиндра насоса начинается только после того, как точка подвеса штанг своим перемещением вверх скомпенсирует деформацию (удлинение) штанг и деформацию (укорочение) труб. Естественно,
при этом длина хода плунжера 5^ по сравнению с длиной хода устьевого штока уменьшается на длину общей деформации
Р
1 = 1Ш + 1Т = —Ж-
Ш Т Е
у
(13)
^ПЛ равно ^ПЛ 1.
Полученные соотношения позволяют записать координаты вершин параллелограмма, представляющего собой общий вид «теоретической» динамограммы в осях 5-— :
Впё = [(51 ^, —0 .ж ),(I о + 1 —0 .ж + —ж );(5 а 6, —0 .ж + —ж ),(5 а о — 1 —0 .ж ^
где ^нМт , 5ВМт - значение позиции подвески устьевого штока соответственно в нижней и верхней мертвых точках кривошипно-коромыслового механизма.
Приняв координату 5НМТ за начало отсчета хода штока по оси 5, получим:
Впё = [(0,—ш.Ж ) , (1,—Ш.Ж + —Ж ) ; (5ШТ ,—Ш.Ж + —Ж ) , (5ШТ - 1,—Ш.Ж )] . (14)
Построенная динамограмма описывает возможные границы изменения уровня жидкости в ШНГУ в текущий момент времени. Покажем на численном примере процедуру использования полученных соотношений для измерения уровня жидкости в ШНГУ.
Пример анализа динамического уровня на построенной модели
Проведем анализ степени влияния на значения статических усилий в подвеске устьевого штока, а соответственно, и на расчетное значение динамического уровня по (6), изменяющихся во времени параметров скважины. К таким параметрам относятся плотность жидкой смеси в НКТ - рж и затрубное давление рз . На рис. 1-5 приведены расчетные зависимости составляющих усилий от указанных параметров. Расчеты проведены в среде МаШСа^ Все зависимости получены для скважины, имеющей следующие характеристики оборудования:
- насос: глубина спуска насоса ЬН = 1401 м, диаметр плунжера насоса ВН = 44 мм;
— колонна штанг: длина штанги 1 ступени /ш1 = 500 м, диаметр штанги 1 ступени Вш1 = 22 мм; длина штанги 2 ступени /ш2 = 900 м, диаметр штанги 2 ступени Вш2 = 19 мм, плотность материала штанг рс = 7,8 г/см3 .
Вес штанг в жидкости, как следует из (4), линейно зависит от плотности жидкой смеси, степень влияния которой может быть оценена возможным диапазоном изменения коэффициента
плавучести штанг Ьарх . В реальном диапазоне изменения плотности жидкости от 0,5 до 2,0 г/см3 значение коэффициента Ьарх будет изменяться не более чем на 12 % относительно средней величины. Для наглядности на рис. 1 приведен график, рассчитанный по (1)-(4), демонстрирующий влияние плотности жидкой смеси на вес штанг в жидкости, из которого следует, что статическое усилие на нижней ветви динамограммы может незначительно изменяться.
Плотность жидкости, г/см3
Рис. 1. Зависимость веса штанг в жидкости от плотности жидкости
Алгоритм расчета полностью соответствует описанной выше процедуре вычислений. Приведем его основные этапы.
1. Выбираем величину кванта времени At, величину Т& базисного времени и полагаем t = 0.
2. Вводим исходные значения параметров, перечисленные в начале раздела.
3. На основе последовательного применения формул (1), (2), (3), (6), (5), (4) находим соответственно ж , 6арх, , рпр , рвс ц и для данного текущего момента времени.
4. По формуле (13) находим X и на основе полученных значений формируем динамограмму динамического уровня жидкости Нд для текущего момента времени.
5. Формируем текущие управляющие решения на основе анализа полученных динамограмм.
6. Увеличиваем текущее время на квант At. Если полученное значение больше Гб, то
процедура прекращается, в противном случае переходим к шагу 3.
Отметим, что вес штанг с жидкостью, как следует из (1)-(3), зависит не только от динамического уровня, но и от плотности жидкости и затрубного давления. Графики на рис. 2 подтверждают этот вывод и показывают, что статическое усилие на верхней ветви динамограммы существенно зависит от динамического уровня и в меньшей степени - от плотности жидкости.
Рис. 2. Зависимость веса штанг с жидкостью от плотности жидкости и динамического уровня: Р3 = 2 атм
Зависимость этого усилия, назовем его «верхним», от затрубного давления намного слабее, чем от плотности жидкости, что достаточно наглядно демонстрируют графики на рис. 3. Заметим, что влияние плотности жидкости на верхнее усилие в наименьшей степени проявляется при низких динамических уровнях, что вполне понятно, т. к. в этой области относительная доля составляющей веса жидкости сравнительно невелика, и, следовательно, ее отклонения не приводят к существенным изменениям результирующего усилия. При этом степень влияния затрубного давления на верхнее усилие от величины динамического уровня не зависит.
. НПО
35000
Рэ. *ги
« О в Т.* 1.Т 2
Плотность жидкости, г/см3
а
ОЛ 1.1 1.4 1.7 1
Плотность жидкости, г/см3
б
Рис. 3. Зависимость веса штанг с жидкостью от плотности жидкости и затрубного давления:
а - Нд = 400 м; б - Нд = 1 000 м
Наконец, для оценки возможного диапазона изменения веса жидкости и сопоставления его с составляющей компенсирующей силы от затрубного давления на рис. 4 и 5 приведены соответствующие зависимости, рассчитанные по (5) и (6). Как следует из анализа этих зависимо -стей, сила затрубного давления может быть на порядок меньше веса жидкости.
Р3, атм
Рис. 4. Зависимость силы затрубного давления от значения давления
Рис. 5. Зависимость веса столба жидкости в скважине от плотности жидкости и затрубного давления:
а - Нд = 400 м; б - Нд = 1 000 м
Приведенные диаграммы иллюстрируют возможности и процедуры использования метода динамограмм в процессе контроля уровня жидкости в скважине.
Анализ адекватности построенной модели
Для проверки адекватности построенной модели целесообразно было бы сравнить значения перечисленных выше характеристик, вычисленные на основе приведенных соотношений, с реальными их значениями, полученными в результате их измерения в реальных условиях эксплуатации. Ввиду отвутствия указанной возможности для анализа адекватности было выполнено математическое моделирование динамического поведения характеристик для двух видов регуляторов - П-регуляторов и ПИ-регуляторов. Известны некоторые особенности, отличающие эти регуляторы, и в процессе моделирования предполагалось выявить эти отличия.
Для моделирования были взяты те же исходные данные, что и в предыдущем разделе при построении численного примера. В качестве системы основных базисных величин для сравнения взяты Тб = Т - базисное время; Нб = Нтах - базисный уровень, соответствующий максимально
возможному динамическому уровню. Производные базисные величины: Рб - давление, Qб - объем подачи и Сб - объем незаполнения рассчитываются через основные величины в соответствии с вышеприведенными уравнениями динамической модели объекта управления.
Учитывая достаточно большую кратность отношения постоянной времени скважины к рекомендуемому периоду квантования, с целью удобства представления результатов основные расчеты выполним, считая широтно-импульсивный модулятор (ШИМ) усилительным звеном. Вначале рассмотрим процессы регулирования с П-регулятором динамического уровня. Рассмотрим расчет звеньев со следующими параметрами: Тс = 10, Кр = 10, Кс = 1. Результаты расчетов для процессов накопления и откачки жидкости в скважине при изменении динамического уровня, а также компенсации скачка пластового давления представлены на рис. 6.
/>а. О, о.е.
1
V
г~
*, о.е.
О 10 20 30 40 50 60 70 80 90 100
•— Динамический уровень ■■ Подача насоса — Заданный уровень
Рис. 6. Процессы регулирования динамического уровня с П-регулятором
Из приведенных диаграмм видно, что в системе, как это известно из опыта практического использования П-регуляторов (см. [1, 2]), присутствует статическая ошибка по уровню, однако процесс при соответствующей настройке регулятора не имеет перерегулирования. При увеличении коэффициента регулятора статическая ошибка уменьшается. При достаточно большом Кр (Кр > 100) система переходит в автоколебательный режим на частоте модуляции ШИМ, характеризующийся чередованием крайних значений периодов включенного и отключенного состояний ШГНУ.
Применение ПИ-регулятора уровня позволяет устранить статическую ошибку в системе (рис. 7; Кр = 10 и Тр = 1), что также известно из опыта практического использования ПИ-регуляторов.
Динамический уровень — Подача насоса -Заданный уровень
Рис. 7. Процессы регулирования динамического уровня с ПИ-регулятором
При больших коэффициентах усиления в данном варианте построения системы возможны субгармонические режимы, которые также наблюдаются в практическом опыте [1]. Таким образом, полученные соотношения позволяют объяснить ряд наблюдаемых на практике свойств поведения динамического уровня жидкости в скважине, что можно считать косвенным обоснованием их адекватности.
Заключение
Полученные в работе соотношения для динамограммы позволяют оценить пределы колебаний текущего уровня жидкости в трубе скважины и на основе этой информации определить наиболее благоприятный режим работы ШНГУ, обеспечивающий наиболее приемлемое сочетание требований как по продуктивности скважины, так и по обеспечению приемлемых условий ее работы. Описанная процедура может быть использована при создании системы автоматического управления процессом работы скважины.
СПИСОК ЛИТЕРАТУРЫ
1. Хакимзянов М. И., Светлакова С. В., Сидоров М. Е. Опыт диагностирования состояния глубиннонасосного оборудования системой динамометрирования ДДС-04 // Обеспечение промышленной безопасности на предприятиях нефтяной и газовой отрасли: сб. науч. тр. - Уфа: Изд-во УГНТУ, 2007. - С. 115-120.
2. Белоусов В. С., Смирнов Д. Н. Математическое моделирование индикаторных диаграмм поршневых насосов сверхвысокого давления // Науч. вестн. Новосибир. гос. техн. ун-та. - 2008. - № 2. - С. 59-70.
3. Имплозионный гидрогенератор давления многократного действия: пат. 2318985 Россия, МПК Е 21 В 37/00 (2006.01), F 16 F 5/00 (2006.01) / Совпель В. В. № 2006125246/03; заявл. 13.07.2006; опубл. 10.03.2008.
4. Добыча при помощи погружного электронасоса, размещенного в водоотделяющей колонне подводной скважины / D. Cocciolone, M. Parker, T. Pitts, M. Ohl // Нефтегазовые технологии. - 2007. - № 12. - С. 55-61.
REFERENCES
1. Khakimzianov M. I., Svetlakova S. V., Sidorov M. E. Opyt diagnostirovaniia sostoianiia glubinno-nasosnogo oborudovaniia sistemoi dinamometrirovaniia DDS-04 [Experience of diagnosting of the condition deep-well pumping units using the dynamometry system]. Obespechenie promyshlennoi bezopasnosti na predpriiatiiakh neftianoi i gazovoi otrasli: cb. nauch. tr. Ufa, UGNTU Publ., 2007, pp. 115-120.
2. Belousov V. S., Smirnov D. N. Matematicheskoe modelirovanie indikatornykh diagramm porshnevykh na-sosov sverkhvysokogo davleniia [Mathematical modeling of indicating diagrams of piston pumps of superhigh pressure]. Nauchnyi vestnikNovosibirskogo gosudarstvennogo tekhnicheskogo universiteta, 2008, no. 2, pp. 59-70.
3. Sovpel' V. V. Implozionnyi gidrogenerator davleniia mnogokratnogo deistviia [Implosion hydro generator of pressure of multiple usage]. Patent RF, no. 2318985, 2008.
4. Cocciolone D., Parker M., Pitts T., Ohl M. Dobycha pri pomoshchi pogruzhnogo elektronasosa, raz-meshchennogo v vodootdeliaiushchei kolonne podvodnoi skvazhiny [Mining with the use of submersible electric pump placed in riser of underwater well]. Neftegazovye tekhnologii, 2007, no. 12, pp. 55-61.
Статья поступила в редакцию 16.01.2013
ИНФОРМАЦИЯ ОБ АВТОРАХ
Исмагулова Агиба Ишимовна — Атырауский институт нефти и газа, Республика Казахстан; доцент кафедры «Электроэнергетика»; umd_aing@mail.ru.
IsmaguIovaAgiba Ishlmovna -Atyrau Institute of Oil and Gas, Republic of Kazakhstan; Assistant Professor of the Department "Electricity"; umd_aing@mail.ru.
Романенко Николай Геннадьевич — Астраханский государственный технический университет; канд. техн. наук, доцент; зав. кафедрой «Электрооборудование и автоматика судов»; nikolayrom@yandex.ru.
Romanenko Nickolai Gennadievich - Astrakhan Technical State University; Candidate of Technical Sciences, Assistant Professor; Head of Department "Electrical Equipment and Automation of Ships"; nikolayrom@yandex.ru.
