Научная статья на тему 'ИДЕНТИФИКАЦИЯ ГИДРОДИНАМИЧЕСКОЙ МОДЕЛИ СКВАЖИНЫ С ЭЛЕКТРОЦЕНТРОБЕЖНЫМ НАСОСОМ ПО ДАННЫМ КОНТРОЛЯ ВОЗМУЩЁННЫХ РЕЖИМОВ ЭКСПЛУАТАЦИИ'

ИДЕНТИФИКАЦИЯ ГИДРОДИНАМИЧЕСКОЙ МОДЕЛИ СКВАЖИНЫ С ЭЛЕКТРОЦЕНТРОБЕЖНЫМ НАСОСОМ ПО ДАННЫМ КОНТРОЛЯ ВОЗМУЩЁННЫХ РЕЖИМОВ ЭКСПЛУАТАЦИИ Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
106
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СКВАЖИНА / ЭЛЕКТРОЦЕНТРОБЕЖНЫЙ НАСОС / ОСЛОЖНЯЮЩИЕ ФАКТОРЫ / ДЕКОМПОЗИЦИЯ / РЕГУЛЯРИЗАЦИЯ / ОЦЕНКИ / ИДЕНТИФИКАЦИЯ / МЕТОД НАИМЕНЬШИХ КВАДРАТОВ / НАБЛЮДАТЕЛИ / ФИЛЬТРЫ СОСТОЯНИЯ / ПРОБНЫЕ ВОЗМУЩЕНИЯ

Аннотация научной статьи по энергетике и рациональному природопользованию, автор научной работы — Соловьев Илья Георгиевич, Говорков Денис Александрович, Цибульский Владимир Романович

Актуальность исследования связана с проблемами создания и сопровождения динамических моделей скважин с электроцентробежным насосом в контурах оперативного контроля и управления состоянием систем, в условиях дрейфа параметров подъёмника и притока. Оценка первопричин нарушения предписанного технологического режима работы скважины и прогноз ближайших негативных последствий имеют важное значение для принятия взвешенных (оптимальных) решений по коррекции режимных состояний и обеспечению стабильных условий эксплуатации. Цель: разработка алгоритма оценивания дрейфующих параметров подъёмника, включая электроцентробежный насос, по методу наименьших квадратов и притока динамической модели скважины и приёмов стабилизации оценок в условиях реальной эксплуатации по данным зашумлённых измерений. Методы: численное интегрирование нелинейных дифференциальных уравнений, методы регрессионного анализа по схеме наименьших квадратов, статистических вычислительных испытаний, пробных возмущений, приёмов регуляризации алгоритмов оценивания. Результаты. Разработан двухэтапный алгоритм оценивания дрейфующих параметров подъёмника (износы, засорения) и притока (параметры продуктивности и подпора пласта) для упрощенного представления комплексной динамической модели скважины класса: «приток-подъёмник- электроцентробежный насос-устье» по данным двухтемповой шкалы контроля режимных состояний системы с реализуемым наблюдателем динамики притока в условиях пробных частотных возмущений. Выводы. Предложенные правила стабилизации оценок параметров модели на основе генерации пробных возмущений, снижения размерности за счет декомпозиции комплексной модели на составные блоки, фильтрации зашумлённых сигналов измерений посредством аппроксимации целыми функциями, реализации динамического наблюдателя притока формируют основу конструкции алгоритма идентификации с надёжным оцениванием основной группы параметров деформации напорной характеристики насоса и притока. Несмотря на низкие показатели точности идентификации гидросопротивления подъёмника, ошибка воспроизведения динамики переменных состояний системы возмущённых режимов эксплуатации с выходов наблюдателя не превышает 2 % что говорит о низкой чувствительности штатных режимов эксплуатации подъёмника в рамках используемых модельных решений к искомой настройке параметра. Требуемая точность идентификации гидросопротивления трубы должна соотноситься с условиями регламента периодической очистки подъёмника от сорбентов.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по энергетике и рациональному природопользованию , автор научной работы — Соловьев Илья Георгиевич, Говорков Денис Александрович, Цибульский Владимир Романович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

IDENTIFICATION OF A WELL HYDRODYNAMIC MODEL WITH AN ELECTRIC SUBMERSIBLE PUMP USING THE PERTURBED OPERATION MODES CONTROL DATA

The relevance of the research is associated with the problems of creating and maintaining dynamic models of electrical submerged pump-equipped well as a part of operational monitoring and control systems under conditions of the elevator and inflow parameters deviation. Assessment of the root causes of well desired operating mode violation and the forecast of the immediate negative consequences are important for balanced (optimal) decisions on operational modes correction and ensuring stable operating conditions. The main aim of the research is to develop the least-squares algorithm - estimating the drift parameters of the elevator including the electrical submerged pump and the inflow of a dynamic well model and methods for stabilizing these estimates in real-time operations based on noisy measurements data. Methods: numerical integration of nonlinear differential equations, regression least-squares analysis methods, statistical computational tests, test perturbations, methods of regularizing estimation algorithms. Results. The authors have developed the two-stage algorithm for estimating the drift parameters of the elevator (wear, impurities) and inflow (productivity and reservoir parameters) to simplify the presentation of a complex dynamic model «inflow-elevator-electrical submerged pump-wellhead» according to two speed scales for monitoring the operating conditions of the system with inflow dynamics observer under conditions of test frequency disturbances. Summary. The proposed rules for stabilizing the model parameters estimates based on the generation of test perturbations, reducing dimension by decomposing the complex model into composite blocks, filtering noisy measurement signals by approximation, implementing a dynamic inflow observer form the basis for identification algorithm design with reliable estimation of the main group of deformation parameters of pump pressure-head characteristics and inflow. Despite the low accuracy of the elevator hydraulic resistance identification, the error in reproducing the dynamics of the variable states of the system in disturbed operating modes from the observer outputs does not exceed 2 %, which indicates the low sensitivity of the standard operating modes of the elevator to the desired parameter setting within the framework of the model solutions used. The required accuracy of the pipe hydraulic resistance identification should be consistent with the conditions for the periodic cleaning of the lift from sorbents.

Текст научной работы на тему «ИДЕНТИФИКАЦИЯ ГИДРОДИНАМИЧЕСКОЙ МОДЕЛИ СКВАЖИНЫ С ЭЛЕКТРОЦЕНТРОБЕЖНЫМ НАСОСОМ ПО ДАННЫМ КОНТРОЛЯ ВОЗМУЩЁННЫХ РЕЖИМОВ ЭКСПЛУАТАЦИИ»

УДК 622.279:681.5

ИДЕНТИФИКАЦИЯ ГИДРОДИНАМИЧЕСКОЙ МОДЕЛИ СКВАЖИНЫ С ЭЛЕКТРОЦЕНТРОБЕЖНЫМ НАСОСОМ ПО ДАННЫМ КОНТРОЛЯ ВОЗМУЩЁННЫХ РЕЖИМОВ ЭКСПЛУАТАЦИИ

Соловьев Илья Георгиевич12,

solovyev@ikz.ru

Говорков Денис Александрович12,

dagovorkov@mail.ru

Цибульский Владимир Романович1,

v-tsib@yandex.ru

1 Федеральный исследовательский центр Тюменский научный центр СО РАН, Россия, 625026, г. Тюмень, ул. Малыгина, 86.

2 Тюменский индустриальный университет, Россия, 625000, г. Тюмень, ул. Володарского, 38.

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

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

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

Результаты. Разработан двухэтапный алгоритм оценивания дрейфующих параметров подъёмника (износы, засорения) и притока (параметры продуктивности и подпора пласта) для упрощенного представления комплексной динамической модели скважины класса: «приток-подъёмник- электроцентробежный насос-устье» по данным двухтемповой шкалы контроля режимных состояний системы с реализуемым наблюдателем динамики притока в условиях пробных частотных возмущений. Выводы. Предложенные правила стабилизации оценок параметров модели на основе генерации пробных возмущений, снижения размерности за счет декомпозиции комплексной модели на составные блоки, фильтрации зашумлённых сигналов измерений посредством аппроксимации целыми функциями, реализации динамического наблюдателя притока формируют основу конструкции алгоритма идентификации с надёжным оцениванием основной группы параметров деформации напорной характеристики насоса и притока. Несмотря на низкие показатели точности идентификации гидросопротивления подъёмника, ошибка воспроизведения динамики переменных состояний системы возмущённых режимов эксплуатации с выходов наблюдателя не превышает 2 % что говорит о низкой чувствительности штатных режимов эксплуатации подъёмника в рамках используемых модельных решений к искомой настройке параметра. Требуемая точность идентификации гидросопротивления трубы должна соотноситься с условиями регламента периодической очистки подъёмника от сорбентов.

Ключевые слова:

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

В работе рассматривается развитие задач по построению и моделированию динамики переходных процессов в скважине с регулируемым электроцентробежным насосом (ЭЦН) [1-5]. Речь пойдет о конструировании алгоритма идентификации модели, приведенной ранее в статье [6], по обновляемым данным контроля ее режимных состояний. Заметим, что даже в условиях регламентированных испытаний по оценке параметров притока [7] проблема надёжности и качества оценивания [8, 9] всегда остаётся первостепенной. В рассматриваемом случае проблема корректности ещё более усугубляется, т. к. размерность вектора оценок удваивается и помимо притока одно-

временно идентифицируются параметры подъёмника, которые в процессе эксплуатации могут дрейфовать. Заметим сразу, что проблема идентификации нестационарных систем в работе не исследуется, так как скорости дрейфов параметров модели много меньше гидродинамики переходных процессов скважины. По сути, нижеизложенное есть обоснование регламента оценивания [10, 11] параметров модели скважины по данным непрерывного контроля «слабо» возмущённых режимов нормальной эксплуатации. Как и в классическом анализе, качество оценок напрямую зависит от точности первичного контроля её состояний [12, 13].

DOI 10.18799/24131830/2020/5/2649

181

Отметим ряд важных моментов.

1. Исходная гидродинамическая модель, ориентированная на задачи оперативного контроля и управления режимами эксплуатации скважин, учитывая взаимодействие главных компонент и действия осложняющих факторов [14-16], по математической конструкции связей (корреляций) между ключевыми точками изначально распределенного подъемника остается упрощенной.

2. Редуцированный (упрощенный) характер описания с минимальным количеством настраиваемых параметров позволяет рассчитывать на создание корректных процедур параметрического оценивания, реализуемых по месту организации служб диспетчерского контроля и сбора первичных данных.

3. Идентификация в классе упрощенных моделей, в нашем случае - это, прежде всего, допущение о постоянстве удельного веса газо-жидкостной смеси (ГЖС) среднего по стволу подъёмника, по сути, означает аппроксимацию реальных наблюдаемых процессов классом заданных уравнений. Возможное снижение точности моделирования, сужающее горизонты прогноза и планирования, может компенсироваться регламентом регулярной подстройки параметров в обязательном порядке при смене режимных состояний системы.

4. Программе производственных испытаний предшествует этап анализа проблемных задач оценивания на модельных экспериментах с «зашумлён-ными измерениями». Результаты данного этапа обсуждаются в заключительной части работы. Поверка работоспособности алгоритма по данным глубинного контроля осуществляется в части идентификации параметров притока по незавершённым [17] пуско-остановочным режимам работы скважины.

Функциональная схема скважины, приведенная к

вертикали, изображенная на рис. 1, иллюстрирует

принятые обозначения переменных состояний (к -

координата дискретного времени) и параметров:

• Pr, p(2,k), p(1 ,k), p(3,k), p(4,k), p(8,k), pL(k) - распределение давлений [МПа] вдоль линии тока от пласта, призабойной зоны (ПЗ), забоя на уровне HR, подвески насоса на уровне HN, затрубья до перепускного клапана, выход насосно-компрессорной трубы (НКТ), до линейного нефтесборника;

• q(k), qN(k), q^k), qL(k) - объемные расходы [м3/сут] ГЖС с обводненностью в и средним удельным весом у по стволу скважины [МПа/м] от забоя q(k)=qN(k)+qI{k) по каналам НКТ, затрубья (ЗТ) и на устье;

• гои, Т rN, rK - известные по данным конструкции гидросопротивления, линеаризованной модели потерь напора на трение на устьевом штуцере, перепускном клапане и в сегментах подъемника: затрубье, НКТ, низ обсадной колонны [сут/м2];

• кусочно-линейная нормированная напорная характеристика насоса задается таблицей своих коэффициентов (Ло(/)Л(0) (рис. 2 в [6]), в общем

случае настраиваемых по каталогам изделии

[18, 19], и параметрами типоразмера силовоИ

/; 0 0\

установки {п ) по номинальным значениям напора [м] и расхода [м3/сут]; переменные состояния по входу: и(А)={1,0} -идентификатор режима работы по типу {ВКЛ, ВЫКЛ}, ю(к) - относительная частота вращения вала двигателя, в случае частотно-регулируемого привода пропорциональная частоте питающего напряжения (в номинале а>(к)= 1).

Рис. 1. Функциональная схема скважины с погружным электронасосом (отдельно выделены измеряемые переменные состояния и идентифицируемые параметры) Fig. 1. Functional diagram of a well with a submersible electric pump (measured state variables and identified parameters are separately highlighted)

Параметры модели, подлежащие оцениванию, делятся на две части.

Первая часть (дрейфующие характеристики подъемника):

• rS - дополнительное гидросопротивление НКТ, обусловленное сорбированием парафинов, солей и прочих поллютантов (изначально rS=0);

• vh, vq - коэффициенты деформации (снижения) напорной характеристики насоса по напору (износ) и расходу (засорение) (изначальное vh=vq=1). Вторая часть (дрейфующие параметры притока):

• pR - давление подпора пласта на контуре питания ПЗ;

• r(2) - гидросопротивление перехода «пласт-ПЗ»;

• /(1) - гидросопротивление перехода «ПЗ-забой» или их обратные величины м(2)=1//(2), ^(1)=1/г(1), именуемые гидропроводностями переходов, которые в совокупности дают оценку коэффициенту продуктивности линейной модели

^(1)^(2)

притока скважины ж, =-;

* w(1) + w(2)

• Г(2) - постоянная времени ПЗ, связанная с параметром ее объемной гидроупругости: Г(2)=/(2)г(2) [20], в первом приближении оценивается по длительности периода окончания переходного процесса на забое скважины - Тквд, например по кривой восстановления давления (КВД) [7, 14]

Т (2) = ГКВД/3. (1)

Обратим внимание на ряд особенностей организации работы процедур параметрического оценивания.

1. Выделенные на рис. 1 контролируемые переменные состояния имеют разномасштабные периоды измерений:

• состояния по давлениям, равно как и по регулируемой частоте

Ь ={Рь (к), Р(4, к), р(8, к), р(3, к), оо(к) к еК) , (2)

считываемые в моменты времени 4 , имеют сравнительно высокую частоту опроса 4+1=4+А/р, А/р~30 с, кеК={0,1,2,...} - идентификатор дискретного времени по шкале «плотных» измерений;

• замеры расходов на устье

=(д1 (к), к е Кч с к) (3)

осуществляются намного реже А/?~30 мин и более, К - шкала «редких» измерений, синхронизированная по времени со шкалой К

2. Контроль давления у приема насоса -р(3,к) - позволяет разбить описание модели и, как следствие, алгоритм оценивания на два блока в соответствии с указанным ранее разделением идентифицируемых параметров. Реализуемая таким образом декомпозиция задачи идентификации шестимерного вектора на две задачи меньшей размерности с векторами

cn = [ri ^ vh f, cr = [PR r(1) r(2)f

(4)

где повышает уровень обусловлен-

ности оценок по методу наименьших квадратов

(МНК).

3. Условия идентифицируемости параметров (4) модели скважины напрямую зависят от уровня информативности (вариабельности) данных контроля (2), (3), что чаще не обеспечивается стационарными режимами нормальной эксплуатации системы. Типовые методики анализа притока по КВД [7, 20] не экономичны в технологиях реального времени, не информативны по состоянию подъемника и могут быть не корректными по нелинейным проявлениям, связанным с наличием свободного газа на забое. Должный уровень информативности данных первичного контроля

обеспечивается в нашем случае режимом пробных частотных возмущений о(к), кеК работы ЭЦН (без останова скважины).

В отличие от классических методов анализа целых решений, например по КВД, где метод пробных возмущений для произвольных о(к) не применим, в рассматриваемой схеме анализа переходных процессов метод пробных возмущений [21] оказывается не только возможным, но и необходимым инструментом повышения корректности решения обратной задачи.

Отмеченное выше разделение модели от состояния р(3,к) на составные блоки сверху и снизу включает:

Канал «насос-НКТ», описываемый в условиях ф,(к)>0 [6, 16] равенством

"и (к)ю(к)2укк0Л0 (I) -

p(3, k) + r

V h0

-u(k )o(k) ^(i)qN (k)

vq

= p(8, k) + r( HN + rqN (k)),

(5)

где I - номер участка напорной характеристики (рис. 2 в [6]), определяемый условием

i = arg

qN(k) o (k) ■ q°

Vq ]qN(i +1), qN(i)]

Возможный режим противотока в НКТ (д^(к)<0, Н(8,к)>0, 1?(к)=1 [6]) в данном описании игнорируется. Канал «затрубье», описываемый равенством

р(3, к) = р(4, к) + +у(нм - Н(4, к) - (1 -14 (к)) гТ (д(к) - дм (к))), (6)

где динамический уровень раздела фаз Н(4,к) с идентификатором отсутствия подачи через перепускной клапан из затрубья (14(к)=1), с площадью кольцевого сечения - 8Т описывается уравнением

Бт^-Н (4, г) = 1М%(<) - 4(1)). са

Канал «низ колоны»

р(1, к) = р(3, к) + у (НК - Ны + Гкд(к)). (7) Блок притока

(Т(2)В +1) р(2, к) = рК - г (2)д(к),

q(k) = w(1)

D = d/dt, p(1, k) (p(2, k) - p(1, k)) P(1, k) + Hr (P GO P(1, k))+

„ _(1 -P)G ■ p0

HR =

(8)

Pgo - Po

где G - газовый фактор нефти [22], pGO - давление насыщения, ниже которого инициируются процессы разгазирования нефти, p0~0,1 МПа - давление атмосферное.

Техника параметрического оценивания по МНК [23, 24], связанная с приведением модели объекта идентификации к регрессионному виду

y(k) = cT v(k),

(9)

основана на построении и последующем решении системы линеиных алгебраических уравнении

F • С = b, F = X M(k)v(k)v(k)T

k eK0

b ) v(k) y(k)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(10)

относительно вектора оценок искомых параметров. Здесь /(к)>0 - весовые коэффициенты ^ /(к) = 1,

кеК

ранжирующие значимость данных первичного контроля (2), (3) с использованием которых формируются вторичные данные 1ц={у(к),у(к) кеК0сК), предназначенные для реализации алгоритма (10). В примерах /(к)=1/М0, где М0 - мощность множества К0 (число слагаемых в сумме, охватывающих значимую часть возмущённого режима работы системы).

Преобразование модели (5) к регрессионному виду

(9)

Ум (к) = С/ УN (к) устанавливает правило формирования пары

Ум(к) = - (р(3, к) - р(8, к))-Нм

Г

vN (k) =

qN(k) u(k)m(k)qN(k) — Xl(i) -u(k)rn(k)2 к" A0(i)

q

в условиях, что оценка Г(2) (1) априорно известна, связано со следующими преобразованиями уравнении (8).

Используя алгебру линейных дифференциальных операторов [25, 26], перепишем модель динамики ПЗ (8) в виде:

р(2,к) = рК-г(2)д(к\ (15)

где ц(к) воспроизводится на выходе реализуемого динамического фильтра состояний [23]

[Т{2)0 + \)д{к) = ч{к). (16)

Подставляя (15) в выражение для притока в забой скважины (8), после приведения к общему знаменателю и деления результата на м>(1)р(1,к), приходим к искомому представлению пары (14) в виде

уК{к) = р{\к\

(17)

УЛ(*) = [1 Sp(\,k)q(k) -q(k)}

где

Sp (1, k) = 1 + -1

(11)

и, согласно (10), алгоритма оценивания неизвестных (дрейфующих) параметров подъемника

Р • См = Ьм . (12)

Реализация алгоритма возможна, если пара (11) доступна формированию по данным 1Р и 1Ч первичного контроля (2), (3). Если выход регрессионной модели у^(к) (11) формируется, согласно (5), по «плотной» шкале кеК измерений 1Р, то восстановление компонент вектора связано с измерениями расходов $у(к). Заметим, что устьевые замеры расходов qL(k) осуществляются по «редкой» временной шкале 1ф кеКщ (3), причем равенство qL(k)=qN(k) имеет место в условиях 14(к)=1, когда подача qт(k) через перепускной клапан отсутствует (6).

Если одновременно с НКТ наблюдается подача через клапан 14(к)=0, Я(4,к)=0, то выполнено qL(k)=q(k) и на основании (6) можно записать

qN (к) = qL (к) -1(1 - 14(к)) • Ут (к),

гт

Ут(к) = — (Р(3, к)-р(4,к))-Нм, (13) 7

что свидетельствует о возможности формирования переменной qN(к) по шкале «редких» замеров - Кч, следовательно, и о реализуемости алгоритма (12).

Построение регрессионной модели для неизвестных параметров (4) притока (8) в виде

.p(1, k)

и алгоритму оценивания по правилу (10)

FR • cд = bR. (18)

Как и в случае с (12), реализация алгоритма (18) возможна, если пара (17) доступна формированию по данным контроля (2), (3).

Согласно (17), формирование пары (yR(k),vR(k)) основывается на расчетах давления -р(1,к), переменной q(k) и ее динамического аналога - д(к), восстановленного с помощью фильтра (16). В отличие от схемы формирования данных (11) для идентификации параметров подъемника (12), которые, по сути, собираются по шкале «редких» измерений - Kq, корректное воспроизведение выхода фильтра на возмущенных режимах эксплуатации динамический системы требует построения пары (17) по более «плотной» временной шкале с периодами квантования At0 менее Г2/30~10...15 мин, что в условиях формирования выборки Iq (3) может быть недостижимым.

Рассмотрим схему построения оценок q(k), qN (k), кеК процессов q(k), qN(k) на основе реализации настраиваемого наблюдателя [27-30].

В условиях реализации первичной процедуры МНК оценивания параметров подъемника: С = [J Sv vh ]T, динамика расхода в НКТ - qN(k),

согласно (5) и (12), может быть оценена расчетным путем на основе данных IP с «плотной» временной шкалой по выражению

qN(k) =

u(k)а(k)2 Укк0Л0 (i) + yN(k)

„ к0 r, + u (k)a(k)S v —r A (i)

q

k eK, (19)

yR (k) = СRT vR (k)

(14)

где учтена блокировка возможного противотока qN(k)<0 герметичным обратным клапаном на выкиде

ЭЦН (знак «+» на позиции индекса) при выключении (м(к)^0) насоса.

В режиме 14(к)=0, согласно (13), искомая оценка притока д(к) рассчитывается с использованием дК(к) по правилу

д(к) = ды (к) + -(1 - 14(к)) • Ут (к).

гт

Но в условиях 14(к)=1, когда динамика движения уровня Н(4,к), равно как и расхода дт(к), остаётся скрытой, приведённое соотношение оказывается неполным для оценивания притока д(к). Для решения возникшей задачи воспользуемся моделью «затрубъя» (6), из которой, согласно (13), в условиях 14(к)=1, следует

д(к) = дм(к)-8тН(А,к),

Н (4, к) = Ны-1 (р(3, к) - р(4, к)) = - ут (к).

У

Отсюда, наблюдатель притока для любых состояний 14(к) принимает вид

д(к) = д„(к) +14(к)Бт ут(к) + -(1 -1 ¿к)) ■ ут(к). (20)

гт

Построение оценки переменной р(1, к) по «плотной» временной шкале ке К осуществляется согласно (7) по выражению

р(1, к) = р(3, к) + у (НЕ - Ны + Гк4(к)), (21)

где также используется ранее восстановленные (20) значения 4(к).

Полученный ранее алгоритм идентификации параметров притока (17), (18) реконструируется с использованием оценок переменных (р(1, к), 4(к), д (к)) (20), (21) к виду

FR ■CR ЬR

(22)

где вместо (16) переменная д (к) воспроизводится по правилу

(Т (2)В +1) 4 (к) = 4(к) (23)

с «плотной» временной шкалой кеКд.

Вышеизложенное является доказательством следующего результата.

Утверждение.

Пусть гидродинамическая модель скважины с частотно-регулируемым ЭЦН (рис. 1), заданная системой уравнений «верха» (5), (6) с неопределенным вектором дрейфующих параметров см и «низа» (7), (8) с неопределенными параметрами притока ск (4), представлена данными контроля 1р и 1Ч (2), (3), с разномасштабными интервалами («плотностью») синхронизированных измерений своих состояний.

Тогда двухэтапная процедура идентификации неизвестных параметров см, с^ модели по данным контроля состояний 1р и а также режима 14(к) подачи из затрубъя на основе МНК с одновременным

оцениванием динамики притока {р(Ц, к), д(к, подачи дИ(к) и уровня H(4,k)=-yT(k) посредством

реализации настраиваемого по оценкам CN наблюдателя (19)-(23) с применением режима пробных частотных возмущений rn(k) ЭЦН задается нижеследующей инструкцией: НАЧАЛО

1. Ввод исходных данных о параметрах конструкции, ЭЦН и флюида.

2. Синхронизация данных выборки IP и Iq с высокой (2) keK и низкой (3) keKq частотами опроса.

3. Формирование данных входо-выходных состояний (11), (12) модели регрессии подъёмника по шкале keKq.

4. Расчёт МНК оценки (12) параметров подъёмника CN.

5. Реализации настраиваемого по CN наблюдателя (19)-(21), (23), с последующей расчётной оценкой графиков динамики притока

^p(1, k), q(kи входо-выходных состояний

его модели регрессии (17) по шкале keK.

6. Расчёт МНК оценки (22) параметров притока CR.

КОНЕЦ.

Тестирование представленного алгоритма идентификации осуществлялось по данным моделирования частотно o(k) возмущённых режимов работы системы с настройками из табл. 1. Данные измерений формировались на основе модельных состояний с добавлением искусственной помехи измерений.

Результаты тестирования свидетельствуют, что по показателю относительной ошибки оценивания:

Sci =(|c -c,\lci)-100% реализуемый по модельным

измерениям алгоритм идентификации доставляет практически точные оценки искомых параметров. Наблюдаемая малая погрешность оценки притока обусловлена ошибкой разностного вычисления производной сигнала qTk) (20). Однако наложение помехи измерений (стандартная функция «rand» ПО MatLab [31]) с нулевым средним и амплитудами: 1 МПа для p(3,k), 10 м3/сут для qL(k) и 0,25 МПа для p(4,k) и p(4,k), уже искажает результат оценивания, что свидетельствует о плохой обусловленности рассматриваемых задач (cond(F.)~106) и указывает как минимум на два момента:

1) необходимость предварительной фильтрации измеряемых сигналов;

2) необходимость применения процедур регуляризации алгоритмов МНК с внесением дополнительной априорной информации об искомых моделях. Учитывая необходимость расчёта производной

(20), решение первой задачи осуществлялось посредством аппроксимации зашумлённых сигналов контроля [32-34] на интервалах постоянства частотных возмущений гладкими кривыми, предоставляемыми инструментальной средой MatLab, модуль Curve Fitting Toolbox [31]. Сравнительный анализ результатов фильтрации позволил выбрать вариант базиса в виде дробно рациональный функции.

Таблица 1. Параметры модели скважины Table 1. Well model parameters

Наименование, [ед. изм.] Name, [units] Значение Value Обозначение Notation

Среднепластовое давление, [МПа] Average reservoir pressure, [MPa] 22 Pr

Глубина забоя по вертикали, [м] Well depth, [m] 3000 Hr

Коэффициент продуктивности скважины, [м3/(МПа-сут)] Well productivity, [m3/(MPaday)] 8 Wr

Коэффициент гидропроводности притока из пласта в ПЗ скважины, [м3/(МПа-сут)] Reservoir productivity coefficient, [m3/(MPaday)] 10 w(1)

Коэффициент гидропроводности притока из ПЗ в забой скважины, [м3/МПа-сут] Inflow productivity coefficient, [m3/(MPaday)] 40 w(2)

Объёмная доля воды, [ед.] Water cut, [units] 0,3 ß

Удельный вес, [МПа/м] Specific gravity, [MPa/m]

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

• нефти/oil 0,68 10-2 Yo

• воды/water 1,1 10-2 Yw

• жидкости/liquid 0,810-2 Y

• газа/gas 0,810-5 YG0

Весовая характеристика газа, [МПа/м] Gas weight characteristic, [MPa/m] 5,6-Ю-4 Yg

Газовый фактор, [ед.] Gas factor, [units] 100 G

Давление в линии нефтесбора, [МПа] Oil collector line pressure, [MPa] 1,5 Pl

Давление насыщения, [МПа] Saturation pressure, [MPa] 18 Pg0

Гидросопротивления, [сут/м2] Hydroresistances, [day/m2]

• обсадной колоны/casing 1 Гк

• НКТ/tubing 1,6 rN

• устьевого штуцера/wellhead choke 1 ru

• затрубного пространства/annulus 0,8 rZT

• перепускного клапана из ЗТ wellhead annulus choke 2 Г GU

• потерь на засорение в НКТ pressure loss due to clogging 1 rs

Параметры насоса Pump parameters

• уровень подвески, [M]/pump depth, [m] 2200 HN

• номинальный напор, [м] nominal head, [m] 1900 h0

• номинальный дебит, [м3/сут] nominal rate, [m3/day] 60 q0

Площади сечений, [м2] cross-sectional areas, [m2]

• НКТ/tubbing 0,002 Sn

• ЗТ/annulus 0,01 St

Постоянная времени ПЗ, [сут] Inflow zone response time, [day] 0,4 T(2)

ная регрессионная модель объекта идентификации (9) имеет вид

у(А-) = с^ (А-) + с2\>2 (А-) + с3у3 (А-), А- е А^ IIК2, (24)

где К1={1,2,_#1} - дискретные моменты времени стационарного режима работы, а К2={Л1+1Д1+2,...Л2} -последующий период возмущений. В условиях зашумлённости данных измерений равновесное состояние системы на периоде ке К1 может быть представлено моделью

у = су— + ^ + ^ с вычисляемыми осреднёнными состояниями

У =ЕУ(к)/м1, V = £V.(к)/м-, . е{1,2,3}.

кеК1 / кеК /

В условиях У1(к)>0, кеК\иК2 справедливо преобразование с переходом от (24) к модели и алгоритму меньшей размерности вида

у2 (к) = 4г(к), Рг

• cz = bz ,

где

C1 = (У - C2V2 + C3V3)/ VV CZ = [C2 cîi-

Fz = X M(k)z(k)z(k)T,

k eK2

bz = X К)z(k)yz(k), К) = VN -

к eK2

N ' 2 N1

а выход и регрессоры укороченной модели удовлетворяют условиям

"T

yz (k) = - У, z(k) =

Vi(k) vi

v2(k )

v3(k )

По аналогии с типовыми подходами ГДИС [7, 14] реконструкция (регуляризация) процедур оценивания (12), (22), сформированных по данным частотно возмущенных режимов работы, осуществлялась введением дополнительного уравнения связи, восстановленного для осреднённых состояний стационарных режимов эксплуатации до периода пробных возмущений. Правило реконструкции с понижением размерности матрицы МНК в следующем. Пусть исход-

у—(к) V— у—(к) V

Именно данная схема редукции была использования для реконструкции процедур оценивания подъёмника (12) и притока (22).

На рис. 2 вынесены типовые графики стационарного и последующих частотно возмущенных режимов работы системы с наложенной помехой измерений для ключевых состояний р(3,к) - давления у приёма и q^,(k) - объёмного расхода на устье. Отдельно на рис. 2 вынесены укрупнённые фрагменты исходных графиков, дополненные для сравнения кривыми результатов выходов наблюдателя, настроенного по результатам идентификации.

В серии из двадцати статистических реализаций наложенных помех измерений отдельно для трёх вариантов реализаций гидродинамики фильтра (23) с параметрами Т(2)={0,3;0,4;0,5} (истинное значение Т(2)=0,4) сигнальная ошибка наблюдения на возмущённых режимах работы модели в среднем по всем экспериментам и сигналам не превышала 2 %. Осред-ненные данные результатов работы алгоритмов оценивания параметров подъёмника и притока сведены в табл. 2 и иллюстрируются диаграммами на рис. 3.

Сравнительный анализ результатов параметрического оценивания свидетельствует, что относительная ошибка идентификации параметров деформации напорной характеристики насоса 8у и уп в среднем не превышает 4 и 6 %, в то время как ошибка идентифи-

кации суммарного гидросопротивления подъёмника Т\ может доходить до 100 %. При этом точность воспроизведения сигнальных процессов (рис. 2) на выходе наблюдателя не выходит за уровень в 2 %. Подобная неравномерность ошибки оценивания имеет место и для параметров притока. Так, для давления

пласта - рд и гидросопротивления - т(1) ближней зоны притока ошибки идентификации во всех модельных экспериментах не превышают 2 и 5 %. В то время как ошибка по т(2) доходит до 30 %, что, опять-таки, не нарушает хорошую сходимость воспроизведения гидродинамики притока.

Рис. 2. Графики возмущённых переменных состояний моделирования скважины и тестовые «измерения» после

наложения помехи максимальной амплитуды Fig. 2. Well model perturbed variable states and test «measurements» with maximum amplitude noise

Таблица 2. Статистика результатов работы алгоритма идентификации параметров скважины по данным модельных экспериментов

Table 2. Statistics of the well parameters identification algorithm results on model experiments data

Схема идентификации подъемника Lift identification scheme Схема идентификации притока Inflow identification scheme

Параметр Parameter Ошибка оценивания, % Estimation error Параметр Parameter Ошибка оценивания, % Estimation error

Описание эксперимента Experiment definition

ri Sv Vh Г1 Sv Vh Pr r(1) r(2) Pr r(1) r(2)

Исходные параметры Initial parameters 2,84 1,125 0,9 - - - 22 0,1 0,025 - - -

Идентификация без помехи и фильтра Identification without noise and filter 2,839 1,125 0,899 0,002 0,005 0,001 21,96 0,098 0,025 0,15 1,2 2,05

Идентификация без помехи с фильтром Identification without noise and with filter 2,45 1,121 0,892 13,7 0,4 0,82 22,96 0,107 0,023 0,17 1,3 5

Идентификация c помехой и фильтром (среднее значение при T(2)=0,3) Identification with noise and filter (average for T(2)=0,3) 1,28 1,11 0,87 76,1 3,5 4,9 21,87 0,10 0,03 0,7 3,0 16,3

Идентификация c помехой и фильтром (среднее значение при T(2)=0,4) Identification with noise and filter (average for T(2)=0,4) 0,93 1,10 0,86 80,3 2,8 4,4 21,99 0,10 0,03 0,8 2,4 17,7

Идентификация c помехой и фильтром (среднее значение при T(2)=0,5) Identification with noise and filter (average for T(2)=0,5) 1,39 1,11 0,87 99,7 4,0 5,1 22,12 0,10 0,02 1,0 3,5 16,8

/'1 <5r Vh " pr r(l) r( 2)

Рис. 3. Диаграммы средних ошибок оценивания параметров подъемника (A) и притока (B) Fig. 3. Diagrams of average errors in estimating elevator (A) and inflow (B) model parameters

Для проверки работы системы по натурным измерениям использовалась выгрузка данных глубинного контроля давления у приёма со станции управления ЭЦН на незавершённом пуско-остановочном режиме работы скважины одного из месторождений Ноябрьского региона. Исходная информация по скважине и устьевым измерениям сведена в табл. 3 и представлена графиками на рис. 4. Нормированная напорная характеристика насоса, используемая при расчете наблюдателя в схеме идентификации притока, соответствует типовой из [6].

Таблица 3. Параметры модели, результаты ГДИС и идентификации Table 3. Model parameters, test and identification results

Конструкция скважины/Well structure

Параметры флюида/Fluid parameters

Параметр/Parameter

Hr

SN

I

Hn I h0

Yo

I

Yw

в I G

Единицы измерения/Units

м m

м m

0,0042 I 0,056

1569 I

м /сут m3/day

МПа/м MPa/m

I

ед. units

I

Значение/Value

2102

1670

200

0,81 10-

1,09-10"'

0,07

13

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Результаты ГДИС Inflow test results

Расчетные параметры Calculated parapeters

Результаты идентификации Identification results_

I w(2) |

Параметр/Parameter

Ж.

WR

T(2)

_£r

W(1)

Wr

Единицы измерения/Units

МПа

MPa

м3/(МПа-сут) m3/(MPa-day)

сут day

сут/м2 day/m2

МПа

MPa

м3/(МПа-сут) m3/(MPa-day) 38,5 I 108 I 28,4

Значение/Value

13,3

31

13,09

5 4.5 Время, сут

Рис. 4. Графики измерения давлений у приёма насоса при останове и повторном запуске, их аппроксимации и результаты моделирования, а также среднесуточные данные о равновесных устьевых давлениях и расходах перед ГДИС и анализируемым переходным процессом Fig. 4. Pump intake pressure charts during stop and consecutive start, their approximations and model results and daily average data of steady intake pressures and rates before well test and analyzed transient process

q

м

2

m

2

1

Как видно из графика, шумы измерений незначительны, поэтому аппроксимация графика точечных измерений давления у приёма насоса непрерывной функцией использовалась для аналитического расчёта производной сигнала p(3,k) для оценки q (k) (20), с постоянным среднесуточным уровнем давления p(4). Результат ЫНК идентификации по правилу (19)—(23) указывают, что ошибка параметрического оценивания относительно данных ГДИС (пластовое давление и коэффициент продуктивности) лежит в пределах 5 %, а сигнальная ошибка по графику p(3,k) не превышает 6 %, что приемлемо для задач диспетчерского контроля режимных условий эксплуатации системы.

Выводы

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

2. Корректная работа алгоритма в условиях непрерывной эксплуатации скважины во многом обеспечивается следующими обстоятельствами:

а наличием точки глубинного контроля давления p(3,k) у приёма насоса, позволяющей декомпозировать комплексную задачу идентификации с плохо различимыми вариациями контролируемых устьевых состояний на две задачи меньшей размерности с информативным сигналом p(3,k) и более высокими факторами надёжности ЫНК оценок; а введением процедуры пробных частотных возмущений, обеспечивающих повышенный уровень информативности данных контроля состояний на переходных режимах, без останова ЭЦН. Точечный контроль на интервалах пробных возмущений позволяет отслеживать дрейфы параметров модели с заданной периодичностью или по запросу оператора, что более практично в задачах реального времени, чем в классических пуско-остановочных регламентах контроля [35]; а предварительной фильтрацией зашумленных данных первичного контроля путём аппроксимации их гладкими кривыми, что позволяет корректно восстанавливать исходные сигналы и их производные, сокращая уровни относительных ошибок оценивания;

СПИСОК ЛИТЕРАТУРЫ

1. Proactive utilization of ESP performance monitoring to enhance productivity / Hashim Al-Sadah, Mohannad Abo Khamseen, Abdullah Al-Ghamdi, Ali Fardan // SPE middle east oil and gas show and conference. - Bahrain, 2019. URL: https://doi.org/10.211S/194925-MS

2. Rapid and comprehensive artificial lift systems performance analysis through data analytics, diagnostics and solution evaluation / Lichi Deng, Ehsan Davani, Hamed Darabi, Vural Sander Suicmez, David Castineira // SPE middle east artificial conference and

• приемами регуляризации на основе осреднения данных стационарных состояний регрес-соров, восстановленных по равновесному состоянию системы до периода пробных возмущений.

3. Проведенные исследования, выявленные проблемные точки конструирования и корректного оценивания свидетельствуют, что рассмотренный вариант решения задачи идентификации, наблюдаемая неравномерная точность оценивания параметров подъёмника и притока никак не закрывают проблему автоматического сопровождения модели скважин с ЭЦН в условиях реальной эксплуатации. Однако требование к точности идентификации параметра r1 подъёмника следует соотносить с условиями инициализации мероприятий по очистке НКТ от сорбентов, когда наблюдаемая производительность системы qL(k) заметно снижается на десять и более процентов. Причины таких последствий система идентификации опознаёт более надёжно.

Типовые правила ГДИС оценивают два главных показателя: подпор пласта на контуре питания - pR и суммарное гидросопротивление притока -rR=r(1)+r(2), связанное с коэффициентом продуктивности скважины условием wR=1/rR. Согласно статистике машинных испытаний (табл. 2), средняя относительная ошибка оценивания параметра rR в нашем случае не превышает 4 %.

Выводы по результатам натурных измерений

1. Динамика давлений для скважины с ЭЦН заметно отличается на пуске и останове, что связано с наличием в системе нелинейного звена - насоса с особым режимом импульсной откачки из затрубья в момент пуска.

2. Рассмотренная схема оценивания учитывает данную особенность и позволяет решать задачу идентификации притока в условиях пусков и остановов и при скачкообразных возмущениях частоты.

3. Предпринятая аппроксимация динамики приза-бойной зоны моделью первого порядка оказывается грубой для «хорошего» воспроизведения реальных переходных режимов притока, что на незавершённых процессах [17] может приводить к смещению оценок пары (pR,wR). Относительно данных ГДИС ошибка оценивания по указанным параметрам в приведённом анализе лежит в пределах 5 %, что приемлемо для задач оперативного контроля причин отклонений технологического режима от плана, но критично для задач с горизонтами длительного прогнозирования.

exhibition, - Bahrain, 2018. URL: https://doi.org/10.2118/192460-MS (дата обращения 10.12.2019).

3. Mriganga Das, Zadil Hayat Bin Naharuddin, Dhuha Al-Khalid. Optimizing production and opération through automated well and network modeling // SPE Kuwait oil and gas show and conférence. - Kuwait, 2015. URL: https://doi.org/10.2118/175205-MS (дата обращения 10.12.2019).

4. Reda El Mahbes, Wilfried Manfoumbi, Benedicte Kadio-Morokro. Real-time remote monitoring to enhance electrical submersible

pump system run life and maximize production // SPE middle east artificial lift conference and exhibition. - Bahrain, 2018. URL: https://doi.org/10.2118/192470-MS (дата обращения 10.12.2019).

5. Gabor Takacs. Three inventions shaping the future of ESP technology // Journal of Petroleum Science and Engineering. -2019. - V. 182. URL: https ://doi.org/10.1016/j.petrol.2019.106330 (дата обращения 10.12.2019).

6. Модель переходных процессов в скважине с частотно-регулируемым электроцентробежным насосом / И.Г. Соловьёв, Д.А. Говорков, С.С. Белашевский, Ю.А. Ведерникова // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2019. - Т. 330. - № 1. - С. 110-120. URL: https://doi.Org/10.18799/24131830/2019/1/55 (дата обращения 10.12.2019).

7. A review paper on the application of deconvolution technique in well test analysis: Tal Block Pakistan case study / Muhammad Tauqeer, Safwan Arshad, Zaheer Ahmed, Jawad Ahmed // SPE/PAPG Pakistan section annual technical conference and exhibition. - Pakistan, 2017. https://doi.org/10.2118/191279-MS

8. Optimal control of field development in a closed loop / E.S. Zakirov, S.N. Zakirov, I.M. Indrupskiy, O.V. Lyubimova, D.P. Anikeev, I.M. Shiryaev, M.N. Baganova // SPE Russian Petroleum Technology Conference. - Russia, 2015. URL: https://doi.org/10.2118/176642-MS (дата обращения 10.12.2019).

9. ESP-DST well testing in a complex reservoir in the Barents Sea: establishing new methodologies and lessons learned / De Leonardis D., Y. Shumakov, K.L. Morton, S. Sarac // SPE Bergen one day seminar. - Bergen, 2017. URL: https://doi.org/10.2118/185892-MS (дата обращения 10.12.2019).

10. Альберт А. Регрессия, псевдоинверсия и рекуррентное оценивание. - М.: Наука, 1977. - 224 с.

11. Robust identification and feedback design: an active noise control case study / R.S. Sánchez Peña, M.A., Cugueró, A. Masip, J. Quevedo, V. Puig // Control Engineering Practice. - 2008. -V. 16. - Iss. 11. - P. 1265-1274.

12. Valuable improvements of data preparation and controlling operational regime practices at mature field as a basic part of a smart field concept usage case study / D. Kyrnaev, V. Fomin, A. Maslanov, V. Karpov, A. Amirov, D. Sleptsov, I. Batilov, K. Ratanov, A. Kozhin, A. Klemba, A. Rybalko, M. Guletzskiy // Offshore technology conference Asia. - Kuala Lumpur, 2018. URL: https://doi.org/10.4043/28293-MS (дата обращения 10.12.2019).

13. Автоматизация сбора и подготовки данных (в)ТМС для проведения гидродинамических исследований скважин с использованием «виртуального расходомера» / А.А. Пашали, М.А. Александров, А.Г. Климентьев, И.Р. Ямалов, А.С. То-польников, А.В. Жонин, А.В. Колонских, В.Г. Михайлов // Нефтяное хозяйство. - 2016. - № 11. - С. 60-63.

14. Ипатов А.И., Кременецкий М.И. Геофизический и гидродинамический контроль разработки месторождений углеводородов. - М.: НИЦ «Регулярная и хаотичная динамика», 2006. -780 с.

15. Соловьёв И.Г. Структурно-функциональный анализ системы оптимального управления для осложнённого фонда скважин с ЭЦН // НТЖ «Автоматизация, телемеханизация и связь в нефтяной промышленности». - 2019. - № 8. - С. 14-19.

16. Вычислительный анализ режимов эксплуатации скважины с ЭЦН в условиях действия газового фактора / С.С. Белашев-ский, Д.А. Говорков, Д.В. Житов, И.Г. Соловьев // Автомати-

зация, телемеханизация и связь в нефтяной промышленности. -2018. - № 8. - С. 26-30.

17. Сергеев В.Л., Нгуен Т.Х.Ф. Модели и алгоритмы адаптивной интерпретации результатов комбинированных газогидродинамических исследований интеллектуальных скважин // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2018. - Т. 329. - № 10. - C. 67-75.

18. Каталог погружного оборудования УЭЦН TPS-Line (REDA). https://www.slb.ru/upload/iblock/8c4/tps_line_esp-booklet_rus.pdf

19. Технический каталог. Погружное оборудование для добычи нефти и комплексный сервис. URL: https://www.rimera.com/upload/catalog/pdf/Alnas-rus.pdf (дата обращения: 03.10.2019).

20. Подземная гидравлика / К.С. Басниев, А.М. Власов, И.Н. Кочина, В.М. Максимов. - М.: Недра, 1986. - 303 с.

21. Растригин Л.А. Системы экстремального управления. - М.: Наука. Физматлит, 1974. - 632 с.

22. Брилл Дж., Мукерджи Х. Многофазный поток в скважинах. -М.; Ижевск: Институт компьютерных исследований, 2006. -384 с.

23. Рао С.Р. Линейные статистические методы и их применения. -М.: Наука. Физматлит, 1968. - 548 с.

24. Соловьев И.Г., Говорков Д.А. Факторы устойчивости МНК-оценок параметров модели притока вертикальной скважины // НТЖ Автоматизация, телемеханизация и связь в нефтяной промышленности. - 2009. - № 9. - С. 31-36.

25. Андреев Ю.Н. Управление линейными конечномерными объектами. - М.: Наука. Физматлит, 1976. - 424 с.

26. Павлов Б.В., Соловьёв И.Г. Системы прямого адаптивного управления. - М.: Наука, 1989. - 136 с.

27. Luenberger D.G. An introduction to observers // IEEE Transactions on Automatic Control. - 1971. - V. 16 - P. 596-602.

28. Коровин С.К., Фомичёв В.В. Наблюдатели состояния для линейных систем с неопределённостью. - М.: Наука. Физмат-лит, 2007. - 224 с.

29. Fomin V., Solovyev I. Adaptive observer design for hydrodynamic states of an oil-well // Proceeding of CSCS. - Bucharest, 2011. -P. 286-289.

30. Nikolov N.N., Lukov V.Y., Alexandrova M.I. Discrete adaptive real-time state observer development using least-squares method // XXVI International Scientific Conference Electronics. - Sozopol, 2017. - P. 1-4.

31. MATLAB - высокоуровневый язык технических расчетов. URL: https://matlab.ru/products/matlab (дата обращения: 03.10.2019)

32. Катковник В.Я. Непараметрическая идентификация и сглаживание данных. - М.: Наука. Физматлит, 1985. - 336 с.

33. Залоруйко Н.Г. Прикладные методы анализа данных и знаний. - Новосибирск: Изд-во Ин-та математики, 1999. - 270 с.

34. Schoukens J., Rolain Y., Pintelon R. On the use of parametric and non-parametric noise-models in time- and frequency domain system identification // 49th IEEE Conference on Decision and Control (CDC). - Atlanta, 2010. - P. 316-321.

35. Сергеев В.Л., Донг Ван Хоанг, Фам Динь Ан. Адаптивная интерпретация результатов гидродинамических исследований горизонтальных скважин на прогнозирующих моделях // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2019. - Т. 330. - № 1. - C. 165-172.

Поступила 05.05.2020 г.

Информация об авторах

Соловьев И.Г., кандидат технических наук, доцент, старший научный сотрудник, Федеральный исследовательский центр Тюменский научный центр СО РАН; доцент, Тюменский индустриальный университет. Говорков Д.А., кандидат технических наук, старший научный сотрудник, Федеральный исследовательский центр Тюменский научный центр СО РАН; доцент, Тюменский индустриальный университет. Цибульский В.Р., доктор технических наук, профессор, ведущий научный сотрудник, Федеральный исследовательский центр Тюменский научный центр СО РАН.

UDC 622.279:681.5

IDENTIFICATION OF A WELL HYDRODYNAMIC MODEL WITH AN ELECTRIC SUBMERSIBLE PUMP USING THE PERTURBED OPERATION MODES CONTROL DATA

Ilya G. Solovyev1,2,

solovyev@ikz.ru

Denis A. Govorkov12,

dagovorkov@mail.ru

Vladimir R. Tsibulsky1,

v-tsib@yandex.ru

1 Tyumen Scientific Centre SB RAS,

86, Malygin street, Tyumen, 625000, Russia.

2 Tyumen Industrial University,

38, Volodarsky street, Tyumen, 625000, Russia.

The relevance of the research is associated with the problems of creating and maintaining dynamic models of electrical submerged pump-equipped well as a part of operational monitoring and control systems under conditions of the elevator and inflow parameters deviation. Assessment of the root causes of well desired operating mode violation and the forecast of the immediate negative consequences are important for balanced (optimal) decisions on operational modes correction and ensuring stable operating conditions. The main aim of the research is to develop the least-squares algorithm - estimating the drift parameters of the elevator including the electrical submerged pump and the inflow of a dynamic well model and methods for stabilizing these estimates in real-time operations based on noisy measurements data.

Methods: numerical integration of nonlinear differential equations, regression least-squares analysis methods, statistical computational tests, test perturbations, methods of regularizing estimation algorithms.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Results. The authors have developed the two-stage algorithm for estimating the drift parameters of the elevator (wear, impurities) and inflow (productivity and reservoir parameters) to simplify the presentation of a complex dynamic model «inflow-elevator-electrical submerged pump-wellhead» according to two speed scales for monitoring the operating conditions of the system with inflow dynamics observer under conditions of test frequency disturbances.

Summary. The proposed rules for stabilizing the model parameters estimates based on the generation of test perturbations, reducing dimension by decomposing the complex model into composite blocks, filtering noisy measurement signals by approximation, implementing a dynamic inflow observer form the basis for identification algorithm design with reliable estimation of the main group of deformation parameters of pump pressure-head characteristics and inflow. Despite the low accuracy of the elevator hydraulic resistance identification, the error in reproducing the dynamics of the variable states of the system in disturbed operating modes from the observer outputs does not exceed 2 %, which indicates the low sensitivity of the standard operating modes of the elevator to the desired parameter setting within the framework of the model solutions used. The required accuracy of the pipe hydraulic resistance identification should be consistent with the conditions for the periodic cleaning of the lift from sorbents.

Key words:

Well, electric centrifugal pump, complicating factors, decomposition, regularization, assessment, identification, least squares method, observers, state filters, test disturbances.

REFERENCES 5. Gabor Takacs. Three inventions shaping the future of ESP tech-

i tt i ■ aicju Ait. a a Vt ill ... .. nology. Journal of Petroleum Science and Engineering, 2019, vol.

1. Hashim Al-Sadah, Mohannad Abo Khamseen, Abdullah Al- J . /lnlnl(:/'? . ,-fñimn«™

Ghamdi, Ali Fardan. Proactive utilization of ESP performance J81 ^'l^ at:u https:;'(doi.org/10.1016/j.petrol.2019.106330

■. ■ . , . .■ ■. mi7 . ., j (accessed 10 December 2019).

monitoring to enhance productivity. SPE middle east oil and gas ^ , , „ . t^ . . . ^tt. ■.

, , r d i. ■ omo i ■] y , 6. Solovyev i.G., Govorkov D.A., Belashevsky S.S., Vedernikova

show and conference. Bahrain, 2019. Available at: TT . i, . ,,r„-,r . .

https://doi.org/10.2118/194925-MS (accessed 10 December 2019). UA TT^ ^TT ^ , T W frequency-regulated

iT-V™ a T-. -u icjc- ESP. Bulletin of the Tomsk Polytechnic University. Geo Assets

2. Lichi Deng, Ehsan Davani, Hamed Darabi, Vural Sander Suicmez, „ „„■{„ , „„„ , ™ » ... t,. ■. ^ . , ,f lift . Engineering. 2019. vol. 330. No 1. pp. 110-120. in Rus. Available David Castineira. Rapid and comprehensive artificial lift systems 5 f. ,,nlimn,ii5isSmin/ii(w jmr.

_f i ■ . , ■ .■ .■ , , at: https://doi.org/10.18799/24131830/2019/1/55 (accessed 10 De-

performance analysis through data analytics, diagnostics and solution evaluation. SPE middle east artificial conference and exhibi- „,- >,,,7, ,T .

tion. Bahrain, 2018. Available at: https://doi.org/10.2118/192460- 7 Mu T* Tauqeer, f hf Í

MS (accessed 10 December 2019) Ahmed. A review paper on the application of deconvolution tech-

3. MngÜngrDas, ZadiT^Ha^at Bin Naharuddin, Dhuha Al-Khalid. ^, test analysis: Tal ,Bl°,ck P^T case st"dy

. .■ . .■ , . . . .. . SPE/PAPG Pakistan section annual technical conference and ex-

Optimizing production and operation through automated well and . . „„,„ > ■, ,,

. . . .■ QnT7 „ .. ., , , , r hibition. Pakistan, 2017. Available at:

network modeling. SPE Kuwait oil and gas show and conference. . ... . .

Kuwait, 2015. Available at: https://doi.org/10.2118/175205-MS „ https://doi.org/10.2118/191279-MS (accessed 10 December 2019).

8. Zakirov E.S., Zakirov S.N., indrupskiy , i.M. Lyubimova O.V.,

4. Re^H MahbeTWilfned Manfoumbi, Benedicte Kadio-Morokro. Anike?v Df, Shiryaev ^ Baganova MN Optimal COntr0l of

Real-time remote monitoring to enhance electrical submersible pump F d Developmentin a ^^ Loop. PR T .Pf,roleum

. ir- . ■ . .■ m¡7 , . , Technology Conference. Moscow, 2015. Available at:

system run life and maximize production. SPE middle east artificial , ,,, ■ „„„,,„„.,^,,„,,0, . .

f.„ r , ■ omo . ■, . / . https://doi.org/10.2118/176642-MS (accessed 10 December 2019).

lift conference and exhibition. Bahrain, 2018. Available at: „ ^ rT „ , T „ „ „„„ „ '

J J De Leonardis D., Shumakov Y., Morton K.L., Sarac S. ESP-DST

https://doi.org/10.2118/192470-MS (accessed 10 December 2019).

well testing in a complex reservoir in the Barents Sea: establishing

new methodologies and lessons learned. SPE Bergen one day seminar. Bergen, Norway, 5 April 2017. Available at: https://doi.org/10.2118/185892-MS (accessed 10 December 2019).

10. Albert A. Regressiya, psevdoinversiya i rekurrentnoe otsenivanie [Regression, pseudo-inversion and recurrence estimation]. Moscow, Nauka Publ., 1977. 224 p.

11. Sánchez Peña R.S., Cugueró M.A., Masip A., Quevedo J., Puig V. Robust identification and feedback design: An active noise control case study. Control Engineering Practice. 2008, vol. 16, Iss. 11, pp. 1265-1274.

12. Kyrnaev D., Fomin V., Maslanov A., Karpov V., Amirov A., Sleptsov D., Batilov I., Ratanov K., Kozhin A., Klemba A., Rybalko A., Guletzskiy M. Valuable improvements of data preparation and controlling operational regime practices at mature field as a basic part of a smart field concept usage case study. Offshore technology conference Asia. Kuala Lumpur, 2018. Available at: https://doi.org/10.4043/28293-MS (accessed 10 December 2019).

13. Pashali A.A., Aleksandrov M.A., Klimentev A.G., Yamalov I.R., Topolnikov A.S., Zhonin A.V., Kolonskih A.V., Mikhaylov V.G. Automatization of collecting and preparation of telemetry data for well testing using «virtual flowmeter». Oil industry, 2016, no. 11, pp. 60-63. In Rus.

14. Ipatov A.I., Kremenetskiy M.I. Geofizicheskiy i gidrodinamich-eskiy kontrol razrabotki mestorozhdeniy uglevodorodov [Geophysical and hydrodynamic control of hydrocarbon field development]. Moscow, Regulyarnaya i khaotichnaya dinamika Publ., 2006. 780 p.

15. Solovyev I.G. Structural-functional analysis of the optimal control system for complicated fund of wells with ESP. Scientific-Technical Journal «Automation, telemechanization and communication in oil industry, 2019, no. 8, pp. 14-19. In Rus.

16. Belashevskiy S.S., Govorkov D.A., Solovyev I.G. Periodical service regulation of complicated wells. Scientific-Technical Journal «Automation, telemechanization and communication in oil industry, 2018, no. 1, pp. 21-26. In Rus.

17. Sergeev V.L., Phuong T.H. Nguyen. Models and algorithms for adaptive interpretation of combined well test data of intelligent wells. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2018, vol. 329, no. 10, pp. 67-75. In Rus.

18. Katalog pogruzhnogo oborudovaniya UECN TPS-Line (REDA) [Catalog of submersible equipment of ESP]. Available at: https://www.slb.ru/upload/iblock/8c4/tps_line_esp-booklet_rus.pdf (accessed 10 December 2019).

19. Tekhnicheskiy katalog. Pogruzhnoe oborudovanie dlya dobychi nefti i kompleksny servis [Technical catalog. Submersible equipment for oil production and comprehensive service.] Available at: https://www.rimera.com/upload/catalog/pdf/Alnas-rus.pdf (accessed 3 October 2019).

20. Basniev K.S., Vlasov A.M., Kochina I.N., Maksimov V.M. Pod-zemnaya gidravlika [Underground hydraulics]. Moscow, Nedra Publ., 1986. 303 p.

21. Rastrigin L.A. Sistemy ekstremalnogo upravleniya [Extreme control systems]. Moscow, Nauka Fizmatlit Publ., 1974. 632 p.

22. Brill Dzh., Mukerdzhi Kh. Mnogofazny potok v skvazhinakh [Multiphase flow in wells]. Moscow; Izhevsk, Computer research institute Publ., 2006. 384 p.

23. Rao S.R. Lineynye statisticheskie metody i ikhprimeneniya [Linear statistical methods and their applications]. Moscow, Nauka. Fiz-matlit. Publ., 1968. 548 p.

24. Solovyev I.G., Govorkov D.A. Least-sqaure estimations of oil well inflow characteristics: factors of stability. Scientific-Technical Journal«Automation, telemechanization and communication in oil industry, 2009, no. 9, pp. 31-36. In Rus.

25. Andreev Yu.N. Upravlenie lineynymi konechnomernymi obyektami [Management of linear finite-dimensional object]. Moscow, Nauka. Fizmatlit. Publ., 1976. 424 p.

26. Pavlov B.V., Solovyov I.G. Sistemypryamogo adaptivnogo uprav-leniya [Direct adaptive control systems]. Moscow, Nauka Publ., 1989. 136 p.

27. Luenberger D.G. An introduction to observers. IEEE Transactions on Automatic Control, 1971, Vol. 16, pp. 596-602.

28. Korovin S.K., Fomichov V.V. Nablyudateli sostoyaniya dlya lineynkh sistem s neopredelennostyu [State observer for linear systems with uncertainty]. Moscow, Nauka. Fizmatlit. Publ., 2007. 224 p.

29. Fomin V., Solovyev I. Adaptive observer design for hydrodynamic states of an oil-well. Proceeding of CSCS -18. Bucharest, 2011. pp. 286-289.

30. Nikolov N.N., Lukov V.Y., Alexandrova M.I. Discrete adaptive real-time state observer development using least-squares method.

XXVI International Scientific Conference Electronics. Sozopol, 2017. pp. 1-4.

31. MATLAB. Available at: https://www.mathworks.com/products/ matlab.html (accessed 3 October 2019).

32. Katkovnik V.Ya. Neparametricheskaya identifikatsiya i sgla-zhivanie dannykh [Nonparametric identification and data smoothing.]. Moscow, Nauka. Fizmatlit. Publ., 1985. 336 p.

33. Zaloruyko N.G. Prikladnye metody analiza dannykh i znaniy [Applied methods of data and knowledge analysis.]. Novosibirsk, Mathematics Institute Publ., 1999. 270 p.

34. Schoukens J., Rolain Y., Pintelon R. On the use of parametric and non-parametric noise-models in time- and frequency domain system identification. 49h IEEE Conference on Decision and Control (CDC). Atlanta, 2010. pp. 316-321.

35. Sergeev V.L., Dong Van Hoang, Pham Dinh An. Adaptive interpretation of the results of horizontal well production testing using forecasting models. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2019, vol. 330, no. 1, pp. 165-172.

Received: 5 May 2020.

Information about the authors

Ilya G. Solovyev, senior researcher, Tyumen Scientific Centre SB RAS; associate professor, Tyumen Industrial University.

Denis A. Govorkov, senior researcher, Tyumen Scientific Centre SB RAS; associate professor, Tyumen Industrial University.

Vladimir R. Tsibulsky, lead researcher Tyumen Scientific Centre SB RAS.

i Надоели баннеры? Вы всегда можете отключить рекламу.