Аннотация: Разработана волновая модель крутильных автоколебаний однородной бурильной колонны в форме торсионного маятника. Выполнено компьютерное моделирование нестационарных крутильных движений долота бурильной колонны при вращении верха колонны с ускорениями, которые моделируют её разгон или торможение.
Ключевые слова: Бурильная колонна, долото, торсионные автоколебания.
Технические науки
УДК 539.3
Глушакова Ольга Владимировна
кандидат физико-математических наук, доцент,
доцент кафедры высшей математики
Национальный транспортный университет
Glushakova Olga
PhD, associate professor
National Transport University
ПЕРЕХОДНЫЕ ПРОЦЕССЫ ДИНАМИЧЕСКОГО КРУЧЕНИЯ КОЛОНН ГЛУБОКОГО БУРЕНИЯ ПРИ ВРАЩЕНИИ
TRANSIENT PROCESSES OF DYNAMIC TORSION OF ROTATING DRILL STRINGS IN DEEP BOREHOLES
Аннотация: Разработана волновая модель крутильных автоколебаний однородной бурильной колонны в форме торсионного маятника. Выполнено компьютерное моделирование нестационарных крутильных движений долота бурильной колонны при вращении верха колонны с ускорениями, которые моделируют её разгон или торможение.
Ключевые слова: Бурильная колонна, долото, торсионные автоколебания.
Summary: A wave model of torsional autovibrations of homogeneous drill string is elaborated in the form of a torsional pendulum. The computer simulation of non-stationary torsion motions of the drill string bits is performed for the cases of the drill string rotations with accelerations. They simulate the effects of the drill string accelerating.
Key words: Drill string, bit, torsional autovibrations.
Постановка задачи нестационарной волновой динамики кручения бурильных колонн. В работах [2-4, 6-8] на основе волновой модели торсионного маятника рассмотрены задачи о самовозбуждении упругих периодических крутильных колебаний долот бурильных колонн при постоянной угловой скорости вращения их верхнего конца.
Однако наряду с явлением самовозбуждения периодических колебаний долота при стационарном вращении верхнего конца бурильной колонны значительный практический интерес представляет также задача моделирования переходных процессов торсионных колебаний долота при изменении величины в состояниях начала вращательного движения, его разгона, перехода к стационарному вращению и торможению. Формулировка этой задачи имеет некоторую специфику, обусловленную необходимостью выбора системы отсчета, обеспечивающую наибольшую простоту разрешающих уравнений. В отличие от случаев стационарного вращения долота или его установившихся колебаний при постоянной, которые удобно рассматривать в системе координат, вращающейся вместе с верхним концом колонны, задачу о нестационарном движении нужно формулировать в инерциальной системе отсчета, поскольку иначе приходится учитывать, что элементы колонны осуществляют сложное вращение с переменной скоростью. Поэтому ниже для вывода разрешающих уравнений использован подход, базирующийся на рассматривании движения системы относительно недвижимой системы отсчета.
Примем, что бурильная колонна длиной своим верхним концом жестко связана с приводным механизмом, который обеспечивает ее вращение с заданной переменной угловой скоростью. К нижнему концу колонны прикреплено долото (рис. 1). Вращение и упругие волны кручения в колонне будем рассматривать в неподвижной системе координат,начало которой совпадает с центром долота, осьнаправлена вверх вдоль оси колонны. При вращении колонны на долото действует крутящий момент, вызванный фрикционным взаимодействием долота с разрушаемой скальной породой. Кинематическое возмущение, действующее на верхний конец БК, и динамическое воздействие на ее нижний конец приводят к формированию в теле БК бегущих упругих волн кручения, определяемых углом закручивания. Эти волны описываются уравнением
(1)
Здесь ― плотность материала, ― его модуль упругости при сдвиге, ― полярный момент инерции площади поперечного сечения трубы БК.
Уравнение (1) имеет решение как в форме синфазных колебаний с модами гармоничных стоячих волн, так и в форме бегущих волн с профилями произвольного вида. Для анализа торсионных колебаний элементов БК будем применять волновые решения уравнения (1). Введя обозначение, где ― скорость распространения поперечной (торсионной) волны, приведем уравнение (1) к стандартной форме
(2)
Здесь точкой обозначено дифференцирования по времени, штрихом ― по независимой переменной.
Можно убедиться [1], что уравнение (2) имеет решение в форме (решение Д'Аламбера).
(3)
где, произвольные непрерывные, не обязательно дифференцируемые функции, первая из которых определяет волну, распространяющуюся в положительном направлении оси, вторая ― в противоположном направлении.
На краю используется уравнение баланса моментов сил инерции, сил трения и сил упругости, которые действуют на долото, условно отделенное от колонны,
(4)
Момент сил инерции подсчитывается по формуле
(5)
где ─ момент инерции долота относительно оси, ─ угловое ускорение долота относительно инерциальной системы координат.
Фрикционный момент является функцией угловой скорости долота. Ее вид определяется твердостью и прочностью разрушаемой породы, диаметром режущей части долота и его конструкцией, а также силой, с которой долото прижимается к дну буровой скважины. Этот момент можно представить в виде зависимости
(6)
для случая вязкого трения и в виде
(7)
для случая проявления трения Кулона. Здесь.
Крутящий момент сил упругости, действующий на долото, вычисляется по формуле
(8)
Если подставить правые части равенств (5) – (8) в уравнение (4), получим уравнение с частными производными относительно. Его можно привести к обыкновенному дифференциальному уравнению, выразив через с помощью фазовых переменных.
После ряда преобразований приводим уравнение (4) к виду
(9)
Это обыкновенное дифференциальное уравнение второго порядка содержит одну искомую переменную ― функцию. Уравнение (9) является уравнением с запаздывающим аргументом.
Интегрирование уравнения осуществляется методом Рунге – Кутта с включением в алгоритмические процедуры искомых функций, зависящих от их значений в моменты времени. Компьютерная реализация такого алгоритма позволила проанализировать процессы нестационарных крутильных движений как долота, так и всех элементов бурильной колонны, и проследить за переходными процессами динамического кручения системы при ее движении вне диапазона автоколебаний системы, а также при ее входе в этот диапазон, переходе через него и выходе из него. Описанный алгоритм дал возможность исследовать влияние величин ускорения и замедления вращения БК и времени прохождения диапазона самовозбуждения автоколебаний. Шаг интегрирования разрешающего уравнения (9) выбирался из условий сходимости вычислительного процесса. Он оказался равным.
Результаты анализа крутильных автоколебаний бурильной колонны при ее ускоренном вращении. С использованием изложенной выше методики выполнены исследования нестационарного динамического поведения долота и БК при различных значениях времени прохождения переменной через диапазон.
Рассмотрены случаи равноускоренного изменения по закону, исследовалось влияние значения на интенсивность нестационарных автоколебаний стальной БК длиной с размерами поперечного сечения .
На рис. 2 – 5 представлены диаграммы, отображающие сценарий самовозбуждения, развития и выхода из режима автоколебаний при (рис. 2, 3) и при (рис. 4, 5).
Литература: