Full Text

Использование штанговых глубинных насосных установок является одним из самых распространенных способов добычи нефти. Выявление факторов, влияющих на динамическое поведение штанговой колонны и развитие методов расчета, способных эти факторы учесть, может помочь в повышении надежности этих установок. Модели, описывающие поведение штанговых колонн, начали разрабатываться с 50-х годов прошлого века [1-3]. Однако в этих моделях расчеты носят в основном квазистатический характер, не учитывается сильно нелинейная природа действующих на колонну сил. Более полными являются подходы, основанные на решении волнового уравнения для штанги [4-7]. Они эффективны для скважин, близких к вертикальным, где на боковую поверхность насосных штанг действует только гидродинамическое трение, линейно зависящее от скорости, и штанговая колонна сама по себе (без насоса и привода) является линейной системой. В наклонных и особенно криволинейных скважинах в дело вступает кулоново трение, действующее на боковую поверхность штанг, и штанговая колонна становится сама по себе нелинейной системой. Для этого случая вышеупомянутый подход неприменим. Работа клапанов насоса и кулоново трение между штангами и насосно-компрессорной трубой приводят к многозначным зависимостям силы от скорости. Возникающие при этом переключения типа механических связей с силового на кинематический и обратно вызывают довольно интенсивные вибрации с частотой собственных колебаний колонны. Для систем, включающих многозначные нелинейные зависимости необходимы специальные методы расчета, например, такие, которые используют формулировку задачи в виде вариационного неравенства [8-10]. Этот подход часто используется при исследовании механических систем с кулоновым трением и динамических задач пластичности [10-13]. Силы, возникающие в насосных системах, в том что касается их зависимости от скорости, имеют некоторое сходство с силами кулонова трения, поэтому данный подход был использован авторами в более ранних работах [14-17]. В этих работах скорость нижнего конца колонны в фазе покоя (когда все клапаны закрыты) предполагается равной нулю. Это соответствует случаю, когда жидкость в цилиндре насоса несжимаема. Однако в реальных условиях жидкость обладает некоторой сжимаемостью, обусловленной в основном наличием газа в газожидкостной смеси. Как следствие, имеет место некоторое движение нижнего конца штанги при закрытых клапанах насоса. В настоящей работе разработана модель, являющаяся развитием моделей, предложенных в [14-17], и учитывающая сильную сжимаемость газожидкостной смеси в цилиндре насоса. В этой модели газовый компонент подчиняется закону Бойля, а сжимаемость жидкостного компонента, растворение газа в жидкости и выделение газа из жидкости, потери давления в клапанах и инерция срабатывания клапанов не учитываются. На базе предложенной модели рассмотрено влияние содержания газа в смеси, заполняющей цилиндр насоса, на динамограмму штанговой колонны в верхнем сечении. Кроме того, рассмотрено различие поведения штанговой колонны в двух случаях: когда газ равномерно распределен (но не растворен) в цилиндре насоса и когда в результате недостаточного заполнения газ в цилиндре насоса занимает локализованный объем. 1. Постановка задачи и решение С точки зрения механики деформируемого твердого тела штанговая колонна представляет собой упругий ступенчатый стержень, расположенный в упругой трубе. Верхний конец стержня совершает кинематически заданное приводом движение, а нижний находится под действием силы, зависящей от направления его движения относительно соответствующего конца трубы, давления на входе насоса и давления на выходе насоса . Величина этой силы (1) где - давление на выходе насоса; - давление в цилиндре насоса; - рабочая площадь плунжера; - площадь сечения штанги возле насоса; - величина силы трения между плунжером и цилиндром насоса; - скорость перемещения нижнего конца стержня и связанного с ним плунжера насоса. В качестве положительного направления перемещения взято направление, соответствующее движению штанговой колонны вниз. Для того чтобы вычислить давление в цилиндре насоса, рассмотрим поведение газожидкостной смеси в насосе. Рабочий цикл состоит из четырёх стадий, которые представлены на рис. 1 как сегменты кривой перемещение плунжера-давление в цилиндре насоса. Точки перехода от одной стадии к другой помечены цифрами от 1 до 4. Соответствующие позиции плунжера показаны на рис. 2. Первая стадия (1-2). Впускной и выпускной клапаны закрыты. Масса газожидкостной смеси и объем жидкостного компонента не меняются. При движении плунжера вниз происходит сжатие газового компонента газожидкостной смеси. Давление в цилиндре насоса повышается от до . Вторая стадия (2-3). Когда давление в насосе достигает , выпускной клапан открывается и жидкость из насоса начинает поступать в насосно-компрессорную трубу. Плунжер движется вниз, Рис. 1. Рабочий цикл насоса Fig. 1. Working cycle of the pump впускной клапан остается закрытым. Давление в цилиндре насоса остается постоянным и равным . Стадия длится до тех пор, пока плунжер не примет крайнее нижнее положение и не начнется его движение вверх. Рис. 2. Положение плунжера в ключевых точках рабочего цикла Fig. 2. Plunger position in the key points of the working cycle of the pump Третья стадия (3-4). С началом движения плунжера вверх выпускной клапан закрывается. Как и на первой стадии, впускной и выпускной клапаны закрыты, а масса газожидкостной смеси и объем жидкостного компонента остаются постоянными. При движении плунжера вверх происходит расширение газового компонента газожидкостной смеси. Давление в цилиндре насоса понижается от до . Четвертая стадия (4-1). Когда давление в насосе достигает , впускной клапан открывается и жидкость из скважины поступает в цилиндр насоса. Плунжер движется вверх, выпускной клапан остается закрытым. Давление в цилиндре насоса остается постоянным и равным . Стадия длится до тех пор, пока плунжер не примет крайнее верхнее положение и не начнется его движение вниз. Тогда впускной клапан закрывается и весь процесс повторяется с самого начала. Рассматриваются два случая присутствия газа в цилиндре насоса: газ распределенный (но не растворенный) в жидкости, и газ, занимающий локализованный объем в цилиндре насоса в результате его недостаточного заполнения (рис. 3). Начнем с ситуации, когда газ в насосе распределен в откачиваемой смеси. В этом случае масса газа неизменна только на тех стадиях, когда оба клапана закрыты. Когда Рис. 3. Два типа присутствия газа в цилиндре насоса: а - газ распределен в смеси; б - газ занимает локализованный объем в цилиндре насоса Fig. 3. Two types of the gas content in the pump: (a) the gas is distributed in the mixture; (b) thegas is localized in some parts of the pump впускной (выпускной) клапан открыт, газ поступает (вытекает) в насос в составе газожидкостной смеси. Поэтому масса газа, а следовательно, и сжимаемость газожидкостной смеси различны на стадии сжатия и расширения. Пусть задано - объемное содержание газа при . При изменении давления благодаря сжатию или расширению газового компонента меняется и текущее объемное газосодержание. Если жидкостный компонент несжимаем, а газовый компонет подчиняется закону Бойля, то зависимость текущего объемного газосодержания от давления может быть записана как (2) На первой и третьей стадии, когда впускной и выпускной клапаны закрыты, давление в цилиндре насоса изменяется согласно закону Бойля: где - объем газа в цилиндре насоса; - объем жидкой составляющей в цилиндре насоса, остающийся постоянным на этих стадиях; и - перемещение плунжера и давление в цилиндре насоса в начале стадии процесса; - расстояние между дном цилиндра и позицией плунжера, соответствующей , т.е. ненапряженному состоянию штанги. Таким образом, получаем выражение для текущего давления в цилиндре насоса (3) и производной давления по перемещению плунжера u (4) где На первой стадии процесса , , а на третьей , . Во второй ситуации, соответствующей неполному заполнению цилиндра, газ в цилиндре насоса занимает локализованный объем, и масса газа неизменна на всех стадиях процесса. В этом случае текущее давление в цилиндре насоса и производная давления по перемещению плунжера u вместо (3) , (4) записываются, как (5) (6) где - давление в начале стадии процесса; - объемы газа и жидкости в начале стадии процесса, отнесенные к единице площади плунжера. Перемещение нижнего конца стержня относительно нижнего конца трубы во всех фазах движения может быть представлено как сумма упругой составляющей и перемещения скольжения: (7) При закрытых клапанах в цилиндре происходит упругая деформация газожидкостной смеси. Нижний конец колонны в этом случае совершает некоторое обратимое движение ue. Давление в насосе (8) Скорость скольжения равна нулю в этой фазе: (9) Это условие может быть представлено в виде вариационного неравенства . (10) Здесь - выпуклый суперпотенциал, определенный, как , где . Численный анализ рассматриваемой задачи выполняется на основе принципа возможных перемещений [14]: (11) где - эквивалентная сила, учитывающая давление, действующее на боковую поверхность насосной штанги [14]; - плотность материала штанг; - коэффициент гидродинамического сопротивления; - пробная функция, равная нулю в точках, где перемещение задано. Исходя из принципа возможных перемещений и субдифференциальных граничных условий, представляя пробную функцию в виде , можно сформулировать следующую задачу: найти поле перемещений u(t), удовлетворяющее динамическому квазивариационному неравенству (12) начальным условиям для u и и такое, что Здесь где - соответствует кулонову трению, действующему на боковую поверхность штанговой колонны. Для дискретизации по времени мы используем схему Ньюмарка, где выражаем перемещения и ускорения в момент времени через скорости в данный момент времени, а также через скорости перемещения и ускорения, вычисленные на предыдущем шаге. Учитывая свойства билинейной формы , скалярного произведения и линейность функционала , малый размер временного шага, где давление можно считать изменяющимся линейно: (13) мы сводим решение динамического квазивариационного неравенства к решению на каждом временном шаге следующей задачи: найти , удовлетворяющее квазивариационному неравенству (14) где . (15) . (16) Здесь и - вспомогательные величины, зависящие от перемещений, скоростей и ускорений, вычисленных на предыдущем временном шаге; и - вычисленные на предыдущем шаге значения давления в цилиндре и производной давления по перемещению. В формулах (15), (16) по сравнению с аналогичными формулами, представленными в [14, 17], появляются дополнительные члены, содержащие и , отвечающие за влияние сжимаемости газа в цилиндре насоса. Как показано в работах [11-13], решение квазивариационного неравенства может быть найдено путем последовательного решения вариационных неравенств. Если размер временного шага мал, то для достижения необходимой точности достаточно одного шага. Решение вариационного неравенства эквивалентно нахождению минимума выпуклого негладкого функционала. С помощью метода конечных элементов [18] задача сводится к минимизации конечномерного выпуклого негладкого функционала, которая может быть реализована путем сведения к последовательности гладких задач нелинейного программирования с помощью метода множителей [19] или алгоритма Удзавы [9]. Детали конечноэлементной реализации и решения задачи с помощью метода множителей можно найти в [14]. Разработанный алгоритм был использован для исследования поведения колонны насосных штанг в условиях, близких к реальным. 2. Примеры решения Рассматриваемая колонна состояла из трех частей (снизу вверх): 1 - диаметр 19 мм, длина 272 м; 2 - диаметр 22 мм, длина 392 м; 3 - диаметр 25 мм, длина 328 м. Максимальный зенитный угол 34 угловых градуса на глубине 560 м. Давление на входе насоса 4,815 МПа, давление на выходе насоса 8,563 МПа, диаметр плунжера 32 мм. Расстояние между дном цилиндра и нулевой позицией плунжера равно 2 м. Частота возмущения 6 ходов в минуту. Длина хода верхнего конца колонны равна 3 м. В первом примере исследовано влияние газосодержания на динамическое поведение штанговой колонны. Представлен случай, когда газ распределен в откачиваемой смеси. Рассматривались три значения газосодержания. На рис. 4 показаны кривые перемещение-сила для верхнего и нижнего концов колонны. Рис. 4. Кривые перемещение-сила (динамограммы) в случае, когда газ распределен в смеси для: (а) верхнего конца колонны; (б) нижнего конца колонны Fig. 4. Displacement-force curves (dynagraphs) for the case when the gas is distributed in the mixture (a) is for the upper end; (b) is for the lower end Видно, что сегмент кривой в фазе сжатия имеет более слабый наклон, чем в фазе расширения. Это объясняется тем, что масса газа, а следовательно, и сжимаемость смеси в фазе сжатия больше, чем в фазе расширения. Этот факт также объясняет ослабление свободных колебаний, которое гораздо более значительно для той фазы, когда плунжер движется вниз, чем для фазы, когда плунжер движется вверх. Чем выше газосодержание, тем меньше наклон сегмента кривой в фазе расширения и особенно в фазе сжатия, тем слабее свободные колебания. Рис. 5. Кривые перемещение-сила (динамограммы) в случае, когда газ занимает локализованный объем в цилиндре насоса для: (а) верхнего конца колонны; (б) нижнего конца колонны Fig. 5. Displacement-force curve (dynagraphs)for the case when the gas is localized in some parts of the pump cylinder (a) is for the upper end; (b) is for the lower end Во втором примере рассматривается ситуация, когда газ занимает локализованный объем в цилиндре насоса. Газосодержание было взято равным 0,5. Результат показан на рис. 5. В этом случае масса газа постоянна во всех фазах процесса. Поэтому, как видно из рис. 5, фрагмент кривой в фазе сжатия имеет ту же форму, что и в фазе растяжения. Повышенного ослабления вибраций при движении плунжера вниз в этом случае не наблюдается. Заключение Разработанная авторами модель продольных колебаний штанговой глубинной насосной установки была усовершенствована с целью учета сжимаемости газожидкостной смеси в цилиндре насоса. Показано значительное влияние газосодержания в смеси, заполняющей насос, на динамическое поведение штанговой колонны. Рассмотрены различия в поведении штанговой колонны для случаев, когда газ распределен в смеси и когда он занимает локализованный объем в цилиндре насоса. В первом случае поведение штанговой колонны значительно различается в фазе сжатия и в фазе расширения, во втором случае этих различий не наблюдается.

About the authors

I N Vasserman

Institute of Continuous Media Mechanics

I N Shardakov

Institute of Continuous Media Mechanics

References

Statistics

Views

Abstract - 182

PDF (Russian) - 54