MODELING THE DYNAMICS OF THE CIRCULATION SUB OF THE DRILL STRING
- Authors: Kobityansky A.E1, Mulkov A.N1, Beloborodov V.S1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: Vol 21, No 4 (2019)
- Pages: 73-80
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mm/article/view/2987
- DOI: https://doi.org/10.15593/2224-9877/2019.4.09
- Cite item
Abstract
The dynamics of the circulation sub used to strengthen the walls of the well during drilling is studied. The design and functioning principle of the sub in the process of their operation is considered. The main stages of the work of the sub are highlighted. The description of each stage is given and the corresponding design schemes are presented, which take into account the geometric, power and operational parameters of the system. Taking into account the accepted assumptions the mathematical models presented in the form of a set of relations characterizing mechanical and hydrodynamic processes of functioning of the translator are constructed. The conditions that impose operational, technological and structural restrictions on the modes of each of the stages of its operation are introduced in the calculation. Under the accepted assumptions, the value of the critical pressure required to push the ball through the seat of the bushing structure is determined. Relations are established between this pressure and the geometric characteristics of the ball. The procedure for calculating the circulation mode of the plugging fluid or drilling fluid when opening the side holes is given. The method of obtaining ratios for determining the pressure on the riser required to ensure the required flow rate of the drilling pump is presented. For each stage of the circulation sub, the characteristic of the algorithms and the order of numerical simulation of the dynamics of the considered design are given. A software module has been formed to implement the solution of the obtained equations. According to the results of mathematical modeling the features of dynamic phenomena in the process of circulating sub are revealed. Numerical values of the system characteristics that affect its dynamics by varying the design and technological parameters of the sub are obtained. The calculation results are formed in the form of tables and graphs. A number of them are presented in the article in the form of corresponding diagrams.
Full Text
В процессе бурения возникает необходимость укрепления стенок скважины [1-3] и исключения поглощения бурового раствора пластами породы. Это достигается применением специальных устройств, среди которых следует выделить конструкции циркуляционных переводников [4-8]. Основная функция циркуляционных переводников, относящихся к клапанным системам, заключается в переключении потока бурового раствора между центральным каналом и затрубным пространством колонны труб, а также в подаче тампонирующих материалов [8]. Проблема выбора конструктивных, режимных и технологических параметров, обеспечивающих эффективное функционирование этих устройств, - важная и актуальная задача, теоретическим исследованиям которой посвящены данные материалы. Работа типовой конструкции циркуляционного переводника состоит из четырех этапов (рис. 1). Когда клапан не активирован, буровой раствор проходит через центральный канал устройства. На первом этапе для активации переводника в бурильную колонну сбрасывается полимерный шар. Этот шар прокачивается буровым раствором до посадки в седло втулки переводника. Течение бурового раствора приостанавливается (см. рис. 1, а). Давление жидкости возрастает. Втулка переводника с двумя боковыми отверстиями перемещается, сжимая пружину. Второй этап работы конструкции начинается после достижения втулкой конечного положения (упора). Открываются циркуляционные порты и буровой раствор поступает в затрубное пространство в обход элементов компоновки низа бурильной колонны (КНБК) (см. рис. 1, б). Третий этап происходит после окончания закачки бурового (тампонирующего) раствора. Буровые насосы отключают. Перепад давления между центральным каналом и затрубным пространством достигает нулевого значения. Сила упругости, создаваемая сжатой пружиной, сдвигает втулку с седлом вверх, преодолевая силы трения. Циркуляционные порты перекрываются (см. рис. 1, в). На четвертом этапе при завершении цикла работы переводника в бурильную колонну сбрасывают два стальных шара. Включаются насосы. Втулка снова сдвигается вниз, открывая порты. За счет давления раствора, подаваемого насосом, шары фиксируются в отверстиях портов (см. рис. 1, г). Циркуляция останавливается. Давление на насосах увеличивают, и полимерный шар продавливается через посадочное седло. Пружина возвращает втулку в исходное положение. Буровой раствор вновь поступает через центральный канал. При закрытых боковых портах все шары системы попадают в шароуловитель. Таким образом, реализация всех четырех этапов возможна при определенных значениях конструктивных, эксплуатационных и технологических параметров системы. Расчетные схемы основных этапов работы переводника (см. рис. 1) составлены с учетом ряда следующих допущений: - на всех участках движения втулки трение в уплотнениях постоянно; - вязкое трение при движении втулки пропорционально ее скорости движения. Обозначения и параметры расчетных схем представлены в табл. 1. а б в г Рис. 1. Расчетные схемы: а - первый этап; б - второй этап; в - третий этап; г - четвертый этап Таблица 1 Параметры расчетных схем и их обозначение Характеристика Обозначение Единицы измерения Давление жидкости p МПа Сила трения в уплотнениях Fтр Н Сила упругости Fупр Н Жесткость пружины с Н/м Внутренний диаметр втулки d м Масса полимерного шара m кг Масса втулки М кг Площадь внутреннего сечения гильзы А м2 Коэффициент вязкого сопротивления h Н·с·м-1 Диаметр уплотнений D м Коэффициент трения в уплотнениях q - Сила давления жидкости Fд = рА Н Диаметр отверстия седла dс м Коэффициент трения между шаром и стенками втулки k1 - Диаметр шара Dш м Модуль Юнга Еш Па Площадь сечения шара Аш м2 На первом этапе работы переводника шар в процессе движения должен находиться в седле. Для исключения продавливания шара и его вылета в накопитель определяется «критическое» значение давления жидкости, при котором шар остается в седле. На рис. 2 представлены расчетные схемы до продавливания шара и в момент его движения с последующим вылетом из седла. а б Рис. 2. 3D-модель и расчетная схема: а - положение шара во втулке до начала продавливания; б - случай продавливания шара При этом приняты следующие допущения: - форма шара при продавливании через втулку - эллипсоид; - относительное движение шар-втулка характеризуется сухим трением. Математическая модель первого этапа работы циркуляционного переводника представляет уравнение совместного движения втулки с шаром, а также условие для давления жидкости, обеспечивающего наличие шара в седле [9-11]: (1) На втором этапе работы переводника определяются потери давления всей циркуляционной системы скважины [12, 13]. Полные потери в скважине - суммарные потери напора во всех элементах циркуляционной системы. С учетом схемы рис. 1, б потери давления в долоте и кольцевом пространстве не учитываются. Потери в бурильных трубах рассчитываются на основе пластической модели Бингама, так как в большинстве своем буровой раствор (или тампонирующая жидкость) относится к неньютоновским жидкостям [14-16]. Последовательность расчета представлена в табл. 2. Таблица 2 Алгоритм определения потерь давления в циркуляционной системе № п/п Параметр Зависимость 1 Средняя скорость жидкости в трубах, м/с 2 Критическая скорость жидкости в трубах, м/с 3 Потери давления в трубах, МПа Турбулентный поток: Ламинарный поток: 4 Перепад давления на портах, МПа 5 Расход жидкости через порты, л/с Примечание: Q - расход, м3/с; D - внутренний диаметр трубы, м; YР - предельное напряжение сдвига, Па; L - длина трубы, м; рст - потери давления на стояке; µ - коэффициент расхода в портах; f - суммарная площадь сечения портов; ρт - плотность тампонирующего раствора, кг/л; PV - пластическая вязкость, 10-3 Па·с. Чтобы шар не продавливался через седло, перепад давления на циркуляционных портах рп не должен превышать значения критического давления ркр (см. формулу (1)). С учетом требуемого расхода через циркуляционные порты решена обратная задача. По перепаду давления на портах определяется необходимое давление на буровых насосах. Решение этой задачи осуществлено в процессе компьютерного моделирования. Уравнение совместного движения гильзы с шаром на третьем этапе работы переводника аналогично соотношению (1) и принятым обозначениям (см. табл. 1). Отличие лишь в том, что сила давления жидкости Fд равна нулю: (2) Расчет ведется до достижения втулкой верхнего положения (Х = 0). Во время четвертого этапа работы переводника втулка возвращается в начальное положение. При этом действуют соотношения, описанные в третьем этапе функционирования переводника. Расчет основных режимных и конструкторско-технологических параметров процесса функционирования циркуляционного переводника осуществлялся с помощью математического моделирования в среде MATLAB методом Рунге - Кутты четвертого порядка [17-19]. В процессе математического моделирования определяется минимальное давление, необходимое для смещения втулки на полный ход, и критическое давление, которое продавливает шар через седло (первый этап функционирования). Расчет второго этапа работы переводника реализован по схеме обратной задачи. Первоначально определяется суммарная площадь сечения портов, площадь бурильных труб, расход жидкости через порты и критическая скорость жидкости в трубах. Также определяются потери в обвязке наземного оборудования. Затем путем сравнения полученного расхода (скорости) жидкости с критическим определяются потери давления в трубах. На завершающем этапе рассчитываются потери давления на стояке как сумма потерь давления в трубах, в обвязке и перепада давления на циркуляционных портах. При расчете третьего этапа значение давления обнуляется, а в начальных условиях вводится соответствующее значение координаты. Результаты моделирования формируются в виде таблиц и графиков. Для оценки влияния конструктивных элементов на динамические характеристики проведено варьирование рядом параметров переводника при фиксированном значении одного из них. Такими параметрами для первого и третьего этапов выделены: количество уплотнений, коэффициент жесткости пружины, ход втулки. Для этих этапов в качестве примера на рис. 3 приведены результаты расчетов, определяющих значения перемещения Х и скорости V втулки при различных коэффициентах жесткости пружины и фиксированных значениях количества уплотнений n = 4 и ходе втулки S = 0,08 м. Из диаграмм расчетов видно, что особенно при движении втулки на третьем этапе работы возможно появление эффекта удара из-за скачка скорости. На первом этапе это воздействие менее значимо по сравнению с третьим этапом, т.е. при достижении втулкой крайних точек хода необходимо сгладить этот эффект. На основе методик [10, 20-22] проведен расчет напряжений, возникающих во втулке в момент удара. С учетом особенностей конструкции устройства напряжения в момент удара при жесткости пружины, равной 7 кН/м, составляют σ = 250 МПа < [σ] = 830 МПа, т.е. не превышают допускаемых значений. В других случаях, например при жесткости пружины в 10 кН/м, σ = 498 МПа < [σ] = 830 МПа. В непредвиденных случаях ударные нагрузки могут увеличиться вследствие неправильно заданных режимов тампонирования или промывки, скачка давления на бурильных насосах и т.п. Для предотвращения удара в конструкцию переводника введено разрезное кольцо, один из вариантов которого представлен на рис. 4. Оно уменьшает скорость подвижного элемента за счет создания дополнительной силы трения, рассчитываемой по методикам [23-25]. Рис. 4. Разрезное кольцо Рис. 3. Результаты влияния жесткости пружины: а - первый этап работы переводника; б - третий этап Схема расположения кольца в механизме показана на рис. 5. Рис. 5. Вариант расположения разрезного кольца в конструкции переводника: 1 - корпус; 2 - втулка; 3 - разрезное кольцо; 4 - упор Для второго этапа варьируемыми параметрами являются перепад давления на портах, диаметр боковых портов и плотность жидкости. Изменение диаметра портов от 0,02 до 0,04 м, увеличение перепада давления на портах в диапазоне 5-15 МПа и варьирование плотностью бурового раствора от 1,1 до 1,5 кг/м3 приводит к изменению расхода жидкости на портах от 2,8 до 11,3 л/с, а также давления на стояке в диапазоне 5,14-15,87 МПа. По результатам проведенных исследований можно сделать следующие выводы: - разработана методика исследования динамики циркуляционного переводника на основных этапах его функционирования; - сформированы расчетная схема и математическая модель динамических процессов переводника с учетом взаимосвязи конструкторско-технологических и эксплуатационных параметров; - проведено математическое моделирование динамики переводника и получены численные значения основных параметров конструкции, обеспечивающие рациональные режимы его функционирования; - полученные результаты позволяют перейти к направленному поиску оптимальных конструктивных и режимных параметров функционирования циркуляционных переводников.About the authors
A. E Kobityansky
Perm National Research Polytechnic University
A. N Mulkov
Perm National Research Polytechnic University
V. S Beloborodov
Perm National Research Polytechnic University
References
- Ваталева М.В. Методы и модели управления рецептурой при производстве тампонажных смесей для крепления скважин с учетом особенностей нефтяных и газовых месторождений: дис. … канд. техн. наук: 05.13.06 / Перм. нац. исслед. политехн. ун-т. - Пермь, 2019. - 192 с.
- Булатов А.И., Проселков Ю.М., Шаманов С.А. Техника и технология бурения нефтяных и газовых скважин: учеб. для вузов / ООО «Недра-Бизнесцентр». - М., 2003. - 1007 с.
- Ивачев Л.М. Промывка и тампонирование геологоразведочных скважин: справ. пособие. - М.: Недра, 1989. - 247 с.
- Циркуляционный переводник: пат. 2658851 Рос. Федерация: МПК Е 21В 34/14, Е 21В 21/10 / Фуфачев О.И., Осипов Д.А. - Заявка № 2017128977; заявл. 14.08.2017; опубл. 25.06.2018, Бюл. № 18.
- Циркуляционный переводник бурильной колонны: пат. 155225 Рос. Федерация: МПК Е 21В 21/10 / Сударев Р.П., Мялицин Н.Ю., Караваев В.М. - Заявка № 2014144890; заявл. 06.11.2014; опубл. 27.09.2015, Бюл. № 27.
- Циркуляционный переводник [Электронный ресурс]. - URL: http://тмс-групп.рф/netcat_files/9/4/ Tsirkulyatsionnyy_ perevodnik.pdf (дата обращения: 14.07.2019).
- Циркуляционный переводник бурильной колонны ООО «ВНИИБТ-Буровой инструмент» / Р.П. Сударев, Ю.А. Коротаев, Н.Ю. Мялицин, В.М. Караваев // Бурение и нефть. - 2018. - № 6. - С. 45-47.
- Циркуляционные переводники PBL [Электронный ресурс]. - URL: http://sbdr.ru/catalog/2 (дата обращения: 14.07.2019).
- Экспериментально-расчетные исследования эффективности различных способов демпфирования колебаний деталей газотурбинных двигателей / Б.Ф. Шорр, Н.Н. Серебряков, А.Н. Стадников, Д.В. Шадрин, Е.С. Руденок, А.В. Каначкин, А.Д. Бортников // Вестник Самар. гос. аэрокосм. ун-та. - 2015. - № 3. - С. 171-182.
- Пановко Я.Г. Основы прикладной теории колебаний и удара. - 3-е изд., перераб. и доп. - Л.: Машиностроение, 1976. - 320 с.
- Писаренко Г.С., Яковлев А.П., Матвеев В.В. Справочник по сопротивлению материалов / под ред. Г.С. Писаренко. - 2-е изд., перераб. и доп. - Киев: Наукова думка, 1988. - 736 с.
- Банникова О.Ю. Совершенствование технологии приготовления и применения буровых растворов на основе сухих полимерных смесей: дис. … канд. техн. наук: 25.00.15. - Уфа, 2015. - 210 с.
- Осложнения и аварии при строительстве нефтяных и газовых скважин: учеб. пособие / С.В. Каменских [и др.]; Ухтин. гос. техн. ун-т. - Ухта, 2014. - 231 с.
- Слезкин Н.А. Динамика несжимаемой жидкости / под ред. Е.И. Обросковой. - М.: Гос. изд-во техн.-теор. лит., 1955. - 520 с.
- Рабиа Х. Технология бурения нефтяных скважин / пер. с анг. В.Г. Григулецкого и Ю.М. Кисельмана; под ред. В.Г. Григулецкого. - М.: Недра, 1989. - 413 с.
- Осипов П.Ф. Гидравлические и гидродинамические расчеты при бурении скважин: учеб. пособие / Ухтин. гос. техн. ун-т. - Ухта, 2004. - 71 с.
- Дьяконов В.П. MATLAB. Полный самоучитель. - М.: ДМК-Пресс, 2012. - 768 с.
- Shampine L.F., Reichelt M.W. The MATLAB ODE suite // SIAM J. on Sci. Comp. - 1997. - No. 18. - Р. 1-22.
- Мироновский Л.А., Петрова К.Ю. Введение в MATLAB: учеб. пособие / СПбГУАП. - 2005. -122 с.
- Биргер И.А., Шорр Б.Ф., Иосилевич Г.Б. Расчет на прочность деталей машин: справ. - 4-е изд., перераб. и доп. - М.: Машиностроение, 1993. - 640 с.
- Одельский Э.Х. Гидравлический расчет трубопроводов разного назначения. - Минск: Вышэйшая школа, 1967. - 103 с.
- Балденко Ф.Д. Расчеты бурового оборудования / РГУ нефти и газа им. И.М. Губкина. - М., 2012. - 428 с.
- Анурьев В.И. Справочник конструктора-машиностроителя: в 3 т. / под ред. И.Н. Жестковой. - 9-е изд., перераб. и доп. - М.: Машиностроение, 2006. - Т. 1. - 928 с.
- Анурьев В.И. Справочник конструктора-машиностроителя: в 3 т. / под ред. И.Н. Жестковой. - 9-е изд., перераб. и доп. - М.: Машиностроение, 2006. - Т. 2. - 960 с.
- Амосов А.А., Дубинский Ю.А., Копченова Н.В. Вычислительные методы для инженеров: учеб. пособие. - М.: Высш. шк., 1994. - 544 с.
Statistics
Views
Abstract - 53
PDF (Russian) - 25
Refbacks
- There are currently no refbacks.