Математическая модель скорости проходки для оптимального управления бурением скважин

Размещено на http://www.allbest.ru/

Математическая модель скорости проходки для оптимального управления бурением скважин

Модель Погарского А.А. [3] для механической скорости бурения хм позволяет неявно учитывать влияние расхода и давления промывочной жидкости и имеет вид:

хм = ,

где G - осевая нагрузка на долото,

n - скорость вращения долота,

Q - расход бурового раствора,

а - коэффициент буримости породы,

b и б - коэффициенты, задающие аналитическую форму модели.

Коэффициент b зависит от расхода раствора Q и гидравлической мощности, приложенной к долоту Nд и для максимальной скорости хм = max определяется как b = 1/. Математическая модель механической скорости проходки (1) доработана автором и отличается от модели Погарского тем, что в неё явно введён параметр Q, что является важным для системы управления процессом бурения.

,

где в - показатель степени, задающий форму функции хм = f(Q),

В практике бурения принято, что основным параметром управления является нагрузка на долото G, остальные параметры принимаются постоянными, их значения зависят от интервала проходимых пород и изменяются в проекте на бурении только при смене твёрдости породы. Поэтому модель (1) отражает принятые бурением принципы - управлять осевой нагрузкой G при фиксированных оборотах n и расходе Q.

Рисунок 1 Зависимости хм = f(Q) при G,n = const;

хм = f(n) при G,Q = const

Зависимость (1) пригодна для поиска оптимальных параметров бурения, однако она позволяет найти только локальный максимум функции скорости по параметру Gопт, т.к. он определяется для фиксированных оборотов долота n и расхода раствора Q.

Физически это объясняется тем, что если G = const, то при малых оборотах долота его зубья недостаточно сцепляются с забоем и скорость бурения мала. При больших n происходит проскальзывание долота по забою, что также снижает механическую скорость. Аналогично для расхода раствора Q: при небольших величинах расхода он не обеспечивает качественную очистку забоя, происходит повторное перемалывание шлама и скорость бурения невелика. При больших расходах напор бурового раствора препятствует вращению долота, т.к. ему приходится преодолевать дополнительное гидравлическое сопротивление в промывочных отверстиях долота и на забое.

Для поиска глобального максимума можно использовать в САУ математическую модель, которая имеет экстремумы для всех трёх параметров управления G, n и Q.

,

bn, bG, bQ -коэффициенты, определяющие аналитическую форму функции хм = f(G,n,Q)

Ограничения:

0 ? G ? Gmax;

0 ? n ? nmax;

0 ? Q ? Qmaх ,

где Gmax, nmax, Qmaх - максимально возможные по ресурсам бурового станка значения G, n, Q.

Модель получена автором в результате регрессионной обработки данных проводки 22 скважин в филиале "Краснодар бурение" ООО "Газпромбурение". Сводные данные по пробуренным скважинам приведены в таблице 1. В ней для определения функции хм = f(Q) были выделены однородные по буримости пород интервалы и фрагменты, в которых осевая нагрузка G и обороты ротора n поддерживались постоянными, изменялись только значения расхода Q и механической скорости хм - столбец хм (Q). Аналогично выделялись фрагменты с постоянными G и Q для получения зависимости хм = f(n)- столбец хм (n).

Полученная приближающая функция имеет вид:

Приближающие функции для хм = f(G), хм = f(n) и хм = f(Q) имеют вид:

,

,

Таблица 1

По уравнениям (3) - (6) рассчитаны данные таблицы 2. В ней приняты обозначения:

Kб, kn, kQ, kG, - размерные коэффициенты пропорциональности для определения регрессионных зависимостей хм = хм (G), хм = хм (n), хм = хм (Q) и хм = хм (GnQ).

bG, bn, bQ - коэффициенты формы кривых для соответствующих функций.

Графики опытных данных бурения и их приближающих функций приведены на рис. 2.

Таблица 2

Графики показывают хорошую точность аппроксимации кривых с опытными данными проводки скважин. Выявленный характер зависимостей известен в практике бурения, он совпадает с кривыми, приводимыми в трудах [1, 2, 3, 4 и др.].

На роль основного уравнения модели механического разрушения породы долотом при гидравлической промывке забоя для детерминировано-адаптивной системы управления процессом бурения могут претендовать зависимости (1) и (2). В практике бурения принято, что управление процессом выполняется, в основном, регулированием осевой нагрузки G при фиксированных n и Q, заданных проектом, однако предпочтение следует отдать модели (2), т.к. она позволяет определять глобальный максимум механической скорости по трём параметрам управления, а не по одному.

Рисунок 2 Графики опытных данных бурения и их приближающих функций для хм = f(G), хм = f(n), хм = f(Q) и хм = f(G,n,Q)

Детерминированность системы обусловлена использованием для управления модели скорости проходки как определённого (детерминированного, не стохастического) уравнения, а адаптивность - тем, что модель подстраивается под условия в скважине изменением её параметрических коэффициентов и показателей степени. Модель используются для постоянного решения на буровой обратной задачи - вычисления параметров модели kб, bn, bG и bQ в ходе процесса бурения с целью оперативного определения управляющих параметров G, n и Q по их текущим значениям в ходе проводки скважины.

Принцип адаптации модели к забойным условиям заключается в следующем: при поступлении новых значений хм, G, n и Q через заданный интервал опроса датчиков САУ, например, через каждые 10 мин. параметры kб, bn, bG, bQ пересчитываются, т.е. модель постоянно уточняется для текущих условий бурения. При каждом пересчёте по заданному критерию, например, "максимум рейсовой скорости" находятся оптимальные управляющие параметры Gопт, nопт и Qопт и передаются оператору бурения (буровому мастеру) для их установки на пульте управления бурением. В течение 10 мин. бурение проводится на этих оптимальных в данное время параметрах, затем полученные реальные значения хм, G, n и Q вновь снимаются с датчиков и вводятся в управляющую программу САУ, модель снова уточняется, находятся оптимальные параметры и т.д. Контроль достижения оптимума производится по минимуму вибрации бурильной колонны. буровой скважина месторождение долото

Немаловажным является и тот факт, что постоянная адаптация модели к условиям на забое снижает требования к точности описания моделью реального процесса, т.к. оперативная подстройка параметрических коэффициентов модели нивелирует все погрешности модели и увеличивает её достоверность.

Поскольку модель постоянно адаптируется к забойным условиям, то в принципе не существенно, какой вид имеет функция хм = f(G,n,Q). Кроме детерминированной модели ТП бурения, в качестве базовой можно использовать стандартный многочлен, например, второго порядка в виде

хм = а0 + а1G + а2 n + а3Q + а4Gn+а5GQ + а6nQ+a7 G2 + а8 n2 + а9Q2 (7)

и корректировать его коэффициенты а0 ,…, а9 как описано выше. При этом происходит кусочно-степенная аппроксимация многочленом (7) реальной функции изменения механической скорости в ходе бурения. При этом возможно получение новых видов моделей для описания ТП бурения.

В принципе для адаптивного управления возможно применение и многочлена первой степени

хм = а0 + а1G + а2 n + а3Q

с кусочно-линейной аппроксимацией функции хм = f(G, n,Q) по коэффициентам а0 ,…, а3, однако точность отработки такой линейной модели в процессе управления бурением может оказаться ниже. Использование многочлена третьей степени не целесообразно, т.к. ведёт к излишней подробности описания процесса бурения и к переусложнению системы пересчёта адаптивных коэффициентов (появятся ещё 9 переменных а10 ,…, а18 при третьих степенях параметров G, n, Q).

Литература

1. Бревдо Г.Д. Проектирование режима бурения. - М.: Недра, 1988. -200 с.: ил.

2. Бревдо Г.Д., Гериш К. Оптимизация параметров режима бурения. Обзорная информация, сер. "Бурение", М., ВНИИОЭНГ, 1980. - 59 с.

3. Григулецкий В.Г. Оптимальное управление при бурении скважин. - М.: Недра, 1988, 229 с.: ил.

4. Справочник инженера по бурению. Под ред. В.И. Мищевича и Н.А. Сидорова. Т. 1, М.: Недра, 1973, 516 с.

Аннотация

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ СКОРОСТИ ПРОХОДКИ ДЛЯ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ БУРЕНИЕМ СКВАЖИН

Цуприков Александр Александрович к.т.н., доцент РИНЦ SPIN-cod=6454-3658

Кубанский государственный технологический университет, Краснодар, Россия

Известные математические модели механической скорости проходки являются функциями двух основных параметров управления - осевой нагрузки на долото и скорости вращения долота. Третий параметр - расход бурового раствора, от которого зависит очистка скважины от выбуренной породы, содержится в явном или неявном виде только в модели А.А. Погарского. На буровой управление процессом проводки скважины выполняется регулированием только одного параметра - нагрузки на долото, при этом остальные параметры поддерживаются постоянными согласно проекту на бурение. Однако, практика бурения свидетельствует, что функция скорости бурения имеет экстремумы и для скорости вращения долота, и для расхода раствора. В работе выполнен регрессионный анализ экспериментальных данных бурения для получения математической модели механической скорости проходки как функции трёх параметров управления, при этом получены модели скорости бурения, зависящие отдельно от каждого из параметров при поддержании двух остальных постоянными. Построены графики приближающих функций с точками данных бурения из рапортов буровых мастеров.

Описана методика адаптивного оптимального управления процессом с помощью полученной модели трёх переменных: по текущим данным бурения через заданный интервал времени с помощью метода наименьших квадратов постоянно пересчитываются коэффициенты модели бурения - тем самым модель постоянно подстраивается под условия на забое скважины. По адаптированной модели определяются оптимальные параметры управления бурением, они устанавливаются на буровом станке и производится разбуривание забоя в оптимальном режиме. При смене породы коэффициенты модели снова подстраиваются под забойные условия, определяются оптимальные для данной породы параметры управления и т.д. Поскольку модель постоянно адаптируется к забою, в качестве модели можно взять стандартный полином второй степени и пересчитывать его коэффициенты. Это позволит также определять новые виды моделей для управления процессом бурения

Ключевые слова: МОДЕЛЬ МЕХАНИЧЕСКОЙ СКОРОСТИ БУРЕНИЯ. АДАПТИВНОЕ ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ

MATHEMATICAL MODEL OF SPEED OF THE DRIVING FOR OPTIMUM CONTROL OF DRILLING OF WELLS

Tsоuprikov Aleksandr Aleksandrovich Cand.Tech.Sci., associate professor RSCI SPIN-code=6454-3658

Kuban state technology university, Krasnodar, Russia

Known mathematical models of mechanical speed of a driving are functions of two key parameters of management - axial load of a chisel and speeds of rotation of a chisel. The third parameter - a consumption of boring solution on which cleaning of a well of destroyed breed depends exists only in the A.A. Pogarsky model in an obvious or implicit way. The drilling management of the process of conducting of the well is carried out by regulation only of one parameter - loads of a chisel, thus other parameters are supported by constants according to the project on drilling. However, practice of drilling testifies that function of speed of drilling has extremum and for the speed of rotation of a chisel, and for a solution consumption. In work the regression analysis of experimental data of drilling for receiving mathematical model of mechanical speed of a driving is made as functions of three parameters of management, thus are received the models of speed of drilling depending separately on each of parameters at maintenance of two others constants. Schedules of the approaching functions with points of data of drilling from official reports of drilling foremen are constructed. The technique of adaptive optimum control of process by means of the received model of three variables is described: according to the current data of drilling through the set time interval by means of a method of the smallest squares drilling model coefficients constantly are recalculated - thereby the model is constantly arranged under conditions on a well face. Optimum parameters of management of drilling are determined by the adapted model, they are established on the drilling rig and the face drilling in the optimum mode is made. At breed change coefficients of model are again arranged under bottomhole conditions, optimum parameters of management for this breed, etc. are defined. As the model constantly adapts to a face, as a model it is possible to take a standard polynomial of the second degree and to recalculate its coefficients. It will also allow defining new types of models for management of drilling process

Keywords: MODEL OF MECHANICAL SPEED OF DRILLING. ADAPTIVE OPTIMUM CONTROL

Размещено на Allbest.ru