САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
На правах рукописи
КИСЕЛЕВА Мария Алексеевна
УСТОЙЧИВОСТЬ И КОЛЕБАНИЯ БУРОВЫХ
УСТАНОВОК
01.02.01 – Теоретическая механика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата физико-математических наук
Санкт-Петербург 2012
Работа выполнена на кафедре прикладной кибернетики математикомеханического факультета Санкт-Петербургского государственного университета.
Научный руководитель: член-корреспондент РАН, доктор физико-математических наук, профессор ЛЕОНОВ Геннадий Алексеевич
Официальные оппоненты: доктор физико-математических наук, профессор БЕЛЯЕВ Александр Константинович (Институт проблем машиноведения Российской академии наук (ИПМаш РАН), заместитель директора по научной работе) доктор физико-математических наук, профессор ЮШКОВ Михаил Петрович (Санкт-Петербургский государственный университет, профессор)
Ведущая организация: Институт прикладной физики Российской академии наук (г. Нижний Новгород)
Защита состоится 1 ноября 2012 г. в 14.00 часов на заседании совета Д 212.232.30 по защите докторских и кандидатских диссертаций при СанктПетербургском государственном университете по адресу: 198504, СанктПетербург, Петродворец, Университетский пр., 28, математико-механический факультет, ауд. 405.
С диссертацией можно ознакомиться в Научной библиотеке им.
М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская наб., 7/9.
Автореферат разослан “” _ 2012 г.
Ученый секретарь диссертационного совета Кустова Е.В.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
В диссертации рассмотрены электромеханические модели буровых установок использующих асинхронный двигатель в качестве привода.
В случае, когда бур является абсолютно твердым, с помощью аналитических и численных методов решена задача о предельной нагрузке.
Для двухмассовой электромеханической модели буровой установки с буром, испытывающим крутильные колебания, проведено компьютерное моделирование.
Актуальность темы. Выход из строя бурового оборудования является частой проблемой в нефтегазодобывающей промышленности.
Поэтому задача исследования переходных процессов возникающих в буровых установках при бурении является актуальной.
Основное внимание в настоящей работе уделено изучению динамики буровых установок с асинхронным электродвигателем.
Разработанный критерий устойчивости для простейшей модели буровой установки позволяет получить допустимую нагрузку на бур при смене среды бурения. В случае двухмассовой модели буровой установки найдены скрытые колебания, т.е. колебания, которые не устанавливаются после переходного процесса из окрестностей стационарных состояний.
Таким образом, поломки бурового оборудования могут быть обусловлены наличием данных колебаний.
Цель работы. Целью работы является разработка математических моделей буровых установок, использующих асинхронный двигатель в качестве привода, и изучение влияния различных нагрузок на данные модели с использованием аналитических и численных методов исследования динамических систем, современных вычислительных средств и специализированных математических пакетов.
Методы исследования. Методы исследования включают в себя аналитические (метод функций Ляпунова, методы исследования дифференциальных уравнений с разрывной правой частью) и численные методы исследования устойчивости динамических систем.
Результаты, выносимые на защиту
.
• Разработана математическая модель буровой установки с абсолютно твердым буром, приводимой в движение асинхронным двигателем.
Решена задача о предельной нагрузке.
несимметричного сухого трения и показано, что при определенных условиях предельно допустимая резкопеременная нагрузка определяется значением максимальной постоянной нагрузки, при которой система имеет стационарный режим.
• Описана двухмассовая математическая модель буровой установки, приводимой в движение асинхронным двигателем. Проведены исследования данной модели в пакете Matlab и найдены скрытые колебания.
Достоверность результатов. Все полученные аналитические результаты математически строго доказаны.
Научная новизна. Все основные результаты, представленные в диссертации, являются новыми.
Теоретическая и практическая ценность. Разработанные в диссертации модели позволяют производить более эффективный анализ работы буровых установок.
Апробация работы. Результаты данной работы докладывались на международных конференциях “4th IEEE International Conference on Nonlinear Science and Complexity” (Будапешт, Венгрия – 2012), “Устойчивость и колебания нелинейных систем управления”, конференция Пятницкого (Москва, Россия – 2012), “Шестые поляховские чтения” (Санкт-Петербург, Россия – 2012), “МКПУ-2011” (Дивноморское, Россия – 2011), конференция памяти В.Я. Ривкинда (Финляндия, Ювяскюля – 2010) и на семинарах кафедры прикладной кибернетики (2009 – 2012).
Публикация результатов. Основные результаты диссертации опубликованы в 6 печатных работах, в том числе в 2 журнальных статьях.
Статьи [1,2] опубликованы в изданиях, рекомендованных ВАК РФ.
В работах [1,3-4,6] соавторам принадлежит общая постановка допустимой резкопеременной нагрузки для модели буровой установки моделирование. В работе [3] Киселевой М.А. введена двухмассовая математическая модель буровой установки, приводимой в движение асинхронным двигателем, и проведен локальный анализ полученной системы. В работах [4,6] Киселева М.А. провела анализ систем буровых установок и компьютерное моделирование, в результате которого были обнаружены скрытые колебания.
Объем и структура диссертации. Диссертационная работа состоит из введения, двух глав, разбитых на параграфы, трех приложений, списка литературы, включающего 103 наименования, изложена на страницах машинописного текста и содержит 35 рисунков.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
двигателей и различных моделей буровых установок и представлен обзор литературы, посвященной изучению буровых установок. Также во введении обосновываются актуальность и научная новизна диссертации, формулируются задачи, решаемые в диссертации.
электромеханической модели буровой установки с асинхронным приводом и решению задачи о предельной нагрузке для данной модели.
В работах H. Nijmeijer, N. van de Wouw, N. Mihailovi и др.
была рассмотрена модель буровой установки, состоящая из верхнего диска, приводимого в движение двигателем, упругого стального стержня, нижнего диска и тормозного механизма. Верхний и нижний диски соединены стальным стержнем и оба диска могут вращаться относительно своей оси (см. Рис. 1). Тормозной механизм используется для моделирования силы трения, действующей на нижний диск.
Дифференциальные уравнения вращения верхнего и нижнего дисков запишутся следующим образом:
Juu (t) + k (u(t) l (t)) + b(u(t) l (t)) + Tf u(u) km u = 0, Здесь u и l - угловые смещения верхнего и нижнего дисков, Ju и Jl - моменты инерции, k, b, km - неотрицательные коэффициенты, u постоянное напряжение, Tf u(u) и Tf l (l ) - моменты трения, действующие на верхний и нижний диски.
Данная двухмассовая модель удобна для проведения анализа поведения модели буровой установки. Однако, недостатком предложенной модели является то, что в ней не учитывается в полной мере динамика двигателя, который приводит в действие верхний диск.
Предположим теперь, что бур является абсолютно твердым телом, жестко связанным с ротором, который вращается под действием магнитного поля, создаваемого статором асинхронного двигателя.
Величина взаимодействия бура с породой определяется величиной момента сопротивления, который возникает в процессе бурения. Такая система испытывает резкопеременные нагрузки при вхождении бура в твердую породу, поэтому необходимо исследовать поведение асинхронного двигателя при скачке нагрузки, т.е. при резком изменении момента сопротивления, действующего на бур.
В качестве уравнений электромеханической модели буровой установки возьмем уравнения асинхронного двигателя, дополненные моментом силы сопротивления бурению Mf :
Здесь – угол поворота бура относительно создаваемого статором магнитного поля, вращающегося с постоянной угловой скоростью mf, i1, i2 – токи в обмотках ротора, R – сопротивление обмоток, L – индуктивность обмоток, SB – магнитный поток сквозь ротор, I – момент инерции бура, – коэффициент пропорциональности, = + mf – угловая скорость вращения бура относительно неподвижной системы координат. Примем, что момент силы сопротивления Mf – кулоновского типа. Здесь, в отличие от классического закона кулоновского трения с симметричной разрывной характеристикой момент силы сопротивления Mf имеет несимметричную характеристику, изображенную на рис. 2.:
M, T0 > 0, число M – принимается достаточно большим. Это адекватно тому, что бурение (или сверление) происходит только при > 0. В реальных системах на переходных режимах такая характеристика “запрещает” переход от положительных к отрицательным. В этом случае система может лишь “застрять” на значениях = 0 на некотором промежутке времени. Эти эффекты часто наблюдаются в процессе бурения.
Заменой переменных система (2) сводится к следующей системе:
электрические величины в обмотках ротора, а переменная s определяет скольжение ротора. Для (s, y) справедливо следующее доопределение по Филиппову:
Введем параметр max =. Локальный анализ состояний равновесия системы уравнений (3) показывает, что при 0 < max система (3) имеет единственное асимптотически устойчивое состояние равновесия.
Действительно, в случае = 0 система (3) имеет единственное асимптотически устойчивое состояние равновесия s = 0, y = 0, x = 0, которое соответствует вращению бура с постоянной угловой скоростью, совпадающей с угловой скоростью вращения магнитного поля (работа на холостом ходу).
При (0, max) система (3) имеет одно состояние равновесия:
где s0 - меньший корень уравнения В этом случае бур вращается в том же направлении, что и магнитное поле, но с меньшей угловой скоростью s0 < mf.
Пусть в момент времени t = происходит скачок нагрузки с 0 до 1, где 0 < 0 < 1 < max. Такая ситуация возникает при переходе в более твердую среду бурения. При = 0 система имеет единственное. Важно, чтобы в новом переходном режиме решение s(t), x(t), y(t) В работе доказаны следующие утверждения:
Теорема 1. Пусть выполнены условия Тогда решение системы (3) с = 1 и начальными данными s( ) = состоянию равновесия этой системы.
Следствие 1. Пусть выполнены условия Тогда решение системы (3) с mf = c, = 1 и начальными данными s( ) = t + к состоянию равновесия этой системы.
Следствие 2. Пусть M – достаточно большое положительное число, mf = c, 0 = 0 и Тогда решение системы (3) с = 1 и начальными данными s( ) = 0, y( ) = 0, x( ) = 0 стремится при t + к состоянию равновесия этой системы.
Рис. 3. 1 область допустимой нагрузки в силу условий теоремы, область допустимой нагрузки, полученная моделированием системы.
В случае 2c2 < для значений 1 2c2, (т.е. условие (10) не выполняется) было проведено компьютерное моделирование системы (3) (область 2 на рис. 3), которое показало, что сохраняется утверждение следствия.
Во второй главе рассматривается электромеханическая модель буровой установки с асинхронным приводом, учитывающая деформацию бура при кручении.
Дополним рассмотренную ранее двухмассовую модель буровой установки уравнениями асинхронного двигателя:
Li1 + Ri1 = SB(sinu)u, Li2 + Ri2 = SB(cosu)u, Здесь u, l угловые смещения ротора и нижнего диска относительно вращающегося магнитного поля, mf – скорость вращения магнитного поля, Tf l (mf + l ) – момент силы трения.
Сделав замену переменных получим систему 5-го порядка Рассмотрим случай, когда сила трения имеет несимметричную характеристику, рассмотренную ранее:
где l = mf u. Здесь по-прежнему M, > 0.
Тогда может быть проведен локальный анализ системы, и справедлива следующая теорема:
имеет одно асимптотически устойчивое состояние равновесия:
где s0 – это наименьший корень уравнения В ходе численного анализа системы удалось обнаружить устойчивые режимы работы буровой установки и режимы, при которых бур останавливается.
Рассмотрим более сложную модель трения – срывное трение.
Пусть теперь момент силы трения имеет следующий вид (см. рис. 4) где Tsl, Tf l, sl, sl и bl – неотрицательные коэффициенты.
Проведение качественного анализа системы подобного типа является сложной задачей из-за наличия срывного трения.
В работе было проведено компьютерное моделирование данной системы.
При определенных параметрах в системе возникают так называемые скрытые колебания – колебания, область притяжения которых не содержит Рис. 5. Скрытые колебания в математической модели буровой установки, приводимой в движение асинхронным двигателем - проекция на {, s, u}.
Рис. 6. Скрытые колебания в математической модели буровой установки, приводимой в движение асинхронным двигателем - проекция на {x, s, u}.
окрестность состояния равновесия. На рис. 5 - 7 изображены найденные в системе устойчивое состояние равновесия и устойчивый предельный цикл, что свидетельствует о наличии скрытых колебаний.
Этот результат показывает, что такой сложный эффект как скрытые колебания возникает даже в достаточно простых моделях. Локализация Рис. 7. Скрытые колебания в математической модели буровой установки, приводимой в движение асинхронным двигателем - проекция на {y, s, u}.
скрытых колебаний является трудной задачей т.к. часто не удается найти скрытые колебания при моделировании системы со случайными данными мала в силу малости области притяжения. Поэтому необходимо развитие аналитических методов исследования подобных систем.
Приложения.
В Приложении 1 представлен компьютерный код алгоритма построения траекторий трехмерной системы буровой установки.
В Приложении 2 представлен компьютерный код алгоритма вычисления и локального исследования состояний равновесия для пятимерной системы буровой установки со срывным трением.
В Приложении 3 представлен компьютерный код построения устойчивого состояния равновесия и скрытых колебаний (устойчивого предельного цикла) для пятимерной системы буровой установки со срывным трением.
Публикации автора по теме диссертации.
2012, Том 444, Вып. 2, с. 160–164.
установок приводимых в движение асинхронным двигателем // Вестник С.-Петерб. ун-та, 2012, сер. 1, Вып. 3, с. 39–41.
3. M. Kiseleva, N. Kuznetsov, G. Leonov, P. Neittaanmki. Mathematical Problem for Drilling System // Abstracts of the International Workshop “Mathematical and Numerical Modelling in Science and Technology”, 2010, Jyvskyl, Finland, p. 16.
4. Киселева М.А., Леонов Г.А. Задача устойчивости буровой установки при резкопеременных нагрузках // Международная научная конференция по механике “Шестые поляховские чтения”, Тезисы докладов, 2012, Санкт-Петербург, Россия, с. 45.
5. Киселева М.А. Математическая модель буровой установки с асинхронным двигателем // XII Международная конференция “Устойчивость и колебания нелинейных систем управления” (конференция Пятницкого), Тезисы докладов, 2012, Москва, Россия, с. 171.
6. Kiseleva M., Kuznetsov N., Leonov G., Neittaanmki P. Drilling Sysa tems Failures and Hidden Oscillations // Proceedings of the 4th IEEE International Conference on Nonlinear Science and Complexity (NSC-2012), 2012, Budapest, Hungary, p. 109.