CC BY
136
УДК 622.276: 519.6
ОПРЕДЕЛЕНИЕ ГРАНИЦ УЧАСТКОВ ДИНАМОГРАММЫ ПРИ ЕЕ ОБРАБОТКЕ
В.Б. Садов
SPECIFICATION OF A DYNAMOMETER CHART BOUNDS AT ITS PROCESSING
V.B. Sadov
Рассматривается задача определения границ участков динамограммы на основе обработки вейвлетных преобразований сигналов.
Проводится анализ различных вариантов базисных функций и их пригодность для локализации границ участков динамограмм.
Сформулированы принципы нового метода определения границ характерных участков динамограммы.
Ключевые слова: динамограмма, обработка.
The article considers a problem of specification of a dynamometer chart bounds on the basis of wavelet signal transformation processing. Different variants of basic functions and their availability to localize dynamometer charts bounds are analyzed. Principles of a new method of specification of a dynamometer chart bounds are formulated.
Keywords: dynamometer chart, processing.
Введение. Одной из ключевых задач при скважинной добыче нефти является диагностика состояния скважины и добывающего оборудования. Самым распространенным способом диагностики является анализ динамограмм, т.е. зависимостей нагрузки на полированном штоке от перемещения штока при его полном двойном ходе. Внутри динамограммы есть участки, характеризующие каждую фазу движения полированного штока. Возможные варианты динамограмм с разными дефектами оборудования приведены, например в [1].
Автоматизированная обработка динамограммы может производиться на основе двух различных подходов к решению этой задачи. При первом подходе на основе всей динамограммы целиком получаем какое-то значение критерия, по которому судим о дефектах скважины и оборудования.
Такой подход применялся, например в [2]. При втором подходе вся динамограмма разбивается на отдельные участки, и дефекты определяются по поведению динамограммы на этих участках. Такой подход описан в [3]. Автору представляется более интересным второй подход, так как обработка всей динамограммы целиком может представлять большую вычислительную и алгоритмическую задачу.
При формализации этого подхода главной задачей будет являться процедура определения точек
Садов Виктор Борисович - канд. техн. наук, доцент кафедры систем управления, Южно-Уральский государственный университет; 8У 2005@inbox.ru____________
начала/конца характерных участков динамограммы. На границах участков из-за изломов (или переходов) меняется спектр сигнала. Это можно использовать для локализации участков динамограммы. Для этой цели можно применить вейвлетное преобразование сигналов, так как сама динамограмма представляет собой периодическую кривую.
Любую функцию s(t), принадлежащую гиль-бертовому пространству, можно разложить в ряд по базису ^mk (0 [4,5]:
= X -тк V mk(t), (1)
m,k=-^
где 5'тк - проекции сигнала на ортогональный базис функций. Они могут быть определены скалярным произведением по формуле
•тк = Vтк ^^ = | • Утк « ' ^ . (2)
Вейвлет-функция у^) должна удовлетворять следующим требованиям [6]: ее интеграл
| у(^) • А = 0; сама функция принадлежит пространству Гильберта 1?[Щ ( | у2^) • А конечен).
Victor Borisovich Sadov - Candidate of Science (Engineering), associate professor of Control Systems Department of South Ural State University; sv 2005@inbox.ru
Вейвлетный базис функционального пространства образуется путем масштабных преобразований и сдвигов порождающего вейвлета у(г):
у(а, Ь, ґ) = |а|
т2
-1/2
у[(ґ -Ь)/ а], а,Ь є Я, а Ф 0,
¥(ґ) є Ьл (Я). (3)
Анализ согласно (3) начинаем со значения а = 1 в сторону увеличения его значения. Для нашего случая ищем проекцию только порождающего вейвлета для всех значений сдвигов Ь . Смотрим поведение этой функции. Вблизи точек перегиба эта функция резко меняет свое значение. Это и можно использовать в качестве индикатора смены участка динамограммы.
Варианты выбора порождающих вейвлетов и результаты моделирования
В качестве примера рассмотрим развернутую динамограмму, изображенную на рис. 1 (величины по осям условные). Ось перемещений дана в нормированном виде в диапазоне от 0 до 1 (0-0,5 - движение полированного штока вверх, 0,5-1 - вниз).
Рассмотрим варианты выбора порождающего вейвлета, так как от него зависят алгоритмы определения точек перегиба динамограммы. Автором было рассмотрено большое количество таких функций, некоторые из которых приведены в данной статье.
В качестве порождающего вейвлета примем функцию
у (г) = е 151 г| • ^(62,8 • г), (4)
где выбор коэффициентов функции обусловлен следующим: условное суммарное время всех участков принимаем равным единице; нас будут интересовать составляющие на порядок большие, чем частота движения полированного штока; затухание порождающей функции на границах ее периода задаем большим 3. График этой функции показан на рис. 2.
На рис. 3 приведены графики развернутой динамограммы и функции А'ь (масштабированные значения).
Видно, что график а'ь по форме практически повторяет график динамограммы и на нем заметны колебания, характерные для порождающей функции. Это говорит о слабой пригодности порождающего вейвлета класса (4) для определения границ участков.
Рассмотрим порождающий вейвлет:
0, г < 0;
1, 0 < г < 0,02;
-1/3, 0,02 < г < 0,08;
0, г > 0,08.
График этой функции приведен на рис. 4, а на рис. 5 приведены масштабированные графики исходной динамограммы и функции а'ь .
Из рис. 5 видно, что вейвлет (5) можно ис-
¥(ґ) =
(5)
Рис. 1. Динамограмма развернутая
Рис. 2. График порождающей функции
Рис. 3. Значения коэффициента разложения и исходная динамограмма
пользовать для определения точек перегиба динамограммы. Подобными свойствами обладает и порождающий вейвлет:
0, г < 0;
50 • г, 0 < г < 0,02;
у (г) = -¡1, 0,02 < г < 0,04; (6)
-0,75, 0,04 < г < 0,08;
0, г > 0,08.
График этой функции приведен на рис. 6, а на рис. 7 приведены масштабированные графики исходной динамограммы и функции а'ь .
Для сравнения на рис. 8 приведены данные по расчету с использованием линейной модели. Коэффициенты модели получались с использованием метода наименьших квадратов. При этом для расчета коэффициентов брался интервал по оси х [Ь - 0,025;Ь + 0,025] с шагом 0,005. Далее вводился аналогично вейвлет-преобразованию сдвиг по всей оси от 0 до 1.
Из всех этих зависимостей видно, что точки начала участков здесь можно определить как точки локальных максимумов/минимумов коэффициентов вейвлет-разложения или точки пересечения оси абсцисс с графиком а1(Ь) после прохождения локального максимума/минимума. Интересным вариантом является комбинация использования
типа у = ^0 + а1 • х по вейвлету (5) и исходная динамограмма
линейной модели и вейвлет-аппроксимации. На рис. 9 приведены графики исходной динамограммы и коэффициента а1 модели у = а0 + а1 • х, полученной в результате обработки коэффициентов вейвлет-разложения типа (5).
Из графиков видно, что началом интервала здесь является точка пересечения оси абсцисс и зависимости а1(Ь) после локального максимума/минимума.
Таким образом, задача свелась лишь к выбору зоны нечувствительности при фиксации локальных максимумов/минимумов.
Общий алгоритм получения границ участков динамограммы
Можно вывести общий алгоритм получения границ интервалов.
1. Всю развернутую диаграмму заполняем по горизонтали точками Ь с каким-либо шагом (можно равномерной сеткой).
2. Для каждой точки исходной динамограммы, используя порождающий вейвлет у(г) согласно (5), вычисляем значения коэффициента разложения в ряд •Ь по формуле (2).
3. Для каждой точки исходной динамограммы, используя модель вида у = а0 + а1 • х, вычисляем зависимость а1(Ь) по методу, например, наименьших квадратов.
4. Для каждой точки коэффициента , используя модель вида у = а0 + а1 • х, вычисляем зависимость а13 (Ь) по методу, например, наименьших квадратов.
5. Путем анализа поведения функций •Ь ,
а1(Ь) и а15 (Ь) определяем точки начала участков.
Примечание: так как вейвлет-функция достаточно быстро затухает, вместо интеграла на бесконечном отрезке из (2) можно использовать конечный отрезок.
Для анализа поведения функции •Ь можно
предложить следующую последовательность действий:
1. Вычисляем точки максимума -тах и минимума 5^п функции •Ь на диапазоне Ь е [0;1].
2. Поскольку -тах и 5^п имеют разные знаки, то в качестве зоны нечувствительности можно использовать зону [§5^п; §£тах], 8е(0;1). Смотрим точки пересечения кривой а'ь с прямыми 8 •тах и 8 5^п. Между ними должна находиться точка локального максимума/минимума. Она является претендентом на точку начала участка.
3. По значению абсциссы и ординаты претендента и значению этих координат ближайших претендентов делается вывод о правильности определения точки начала участка. Для этого несложно сформировать какой-либо алгоритм фиксации точки начала каждого участка с учетом поведения и значений точек исходной динамограммы при этом значении абсциссы Ь .
4. Если алгоритм фиксации точек дает неверный результат, то меняем зону нечувствительности (меняем значение 8 ) и продолжаем алгоритм с действия 2.
Для анализа поведения коэффициентов разложения функций а1(Ь) и а13 (Ь) можно предложить следующую последовательность действий.
1. Вычисляем точки максимума атах и минимума а^п функции а1(Ь) на диапазоне Ь е [0;1].
2. Поскольку атах и amin имеют разные знаки, то в качестве зоны нечувствительности можно использовать зону [8 amin; 8 атах ], 8е (0;1). Смотрим точки пересечения кривой а1(Ь) с прямыми 8 атах и 8 а^. Между ними должна находиться точка локального максимума/минимума. Точка пересечения с абсциссой после локального максимума/минимума - претендент на точку начала участка. Надо также следить, пересекают ли линии зависимости а1(Ь) слева и справа от точки локального максимума/минимума ось абсцисс.
3. По значению абсциссы и ординаты претендента и значению этих координат ближайших претендентов делается вывод о правильности определения точки начала участка.
4. Если алгоритм фиксации точек дает неверный результат, то меняем зону нечувствительности (меняем значение 8 ) и продолжаем алгоритм с действия 2.
Проверка одного алгоритма на модельных примерах
Рассмотрим работу алгоритма при использовании функции а1(Ь), полученной на основе обработки результатов линейной аппроксимации функции •Ь , использующей порождающий вейвлет (5). Примерный вид такой обработки графически показан на рис.10.
Рис. 10. Результаты обработки динамограммы (определение границ участков)
Здесь штриховкой показаны участки выхода графика а1(Ь) за границы зоны нечувствительности, а линиями - точки-претенденты на начало участка (потенциально существенные точки перегиба). Из графика видно, что эти точки получаются с некоторой (хотя и вполне допустимой) погрешностью. Это связано с конечностью импульса порождающей функции (рис. 10) и сглаживанию благодаря используемому методу наименьших квадратов (для вычисления коэффициента а1 использовался участок [Ь - 0,025;Ь + 0,025]).
Вообще говоря, точки развернутой динамограммы снимает какое-либо цифровое устройство (обычно микроконтроллер), поэтому можно использовать быстрые цифровые алгоритмы, но в данном случае это несущественно, так как динамограммы не обрабатываются на контроллере в реальном масштабе времени.
Рассмотрим развернутую динамограмму, изображенную на рис. 11 (величины по осям условные). Ось перемещений дана в нормированном виде в диапазоне от 0 до 1. Эта динамограмма не является стандартной и составлена для случая удара плунжера внизу и большого давления газа у приема насоса. Форма динамограммы была немного искажена, имитируя реальные условия съема динамограммы.
На рис. 12 приведены масштабированные графики исходной динамограммы и функции , полученной с использованием вейвлета (5).
На рис. 13 приведены графики исходной динамограммы и коэффициента а1 модели у = а0 + а1 • х, полученной в результате обработки коэффициентов вейвлет-разложения типа (5).
Результаты обработки динамограммы по описанному ранее алгоритму приведены на рис. 14.
Из рис. 14 видно, что точка, соответствующая
обратному выбегу кривой динамограммы (третья слева вертикальная линия), может быть получена как точка участка. Здесь надо учитывать:
- предполагаемый номер участка;
- значение усилия на штоке в начале и конце предполагаемого участка;
- общий диапазон изменения усилия на предполагаемом участке;
- общий диапазон усилий динамограммы.
Несложно сформировать набор таких правил,
характерных для каждого участка.
Примечание: так как возможны выбросы точек динамограммы из-за сбоев аппаратуры и технических причин (легкое заедание механизмов и пр.), то для логического исследования полученных точек лучше всего использовать не исходную динамограмму, а ее усредненную кривую. Она отражает все особенности исходной динамограммы, но является немного сглаженным ее вариантом. На рис. 15 приведены графики исходной и усредненной динамограммы, каждая точка которой есть среднее арифметическое точек исходной динамограммы О(-) на интервале [• - 0,025; • + 0,025]. При снятии динамограммы, например, микроконтроллером с АЦП, исходная динамограмма, в отличие от показанной на рис. 15, будет еще и зашумлена. Усреднение уменьшит величину шума.
°0 0,2 0,4 0,6 0,8 1
Рис. 11. Динамограмма развернутая
Рис. 13. Значения коэффициента а1 разложения типа у = а0 + а1 • х по вейвлету (5) и исходная динамограмма
Рис. 12. Значения коэффициента разложения и исходная динамограмма
..........
‘4 <///
0 0,2 0,4 0,6 0,8 1
Рис. 14. Результаты обработки динамограммы (определение границ участков)
°0 0,2 0,4 0,6 0,8 1
Рис. 15. Исходная и усредненная динамограммы
Литература
1. Мищенко, И. Т. Скважинная добыча нефти: учебн. пособие для вузов / И.Т. Мищенко. - М.: ФГУП Изд-во «Нефть и газ» РГУ нефти и газа им. И.М. Губкина, 2003. - 816 с.
2. Тагирова, К.Ф. Автоматизация управления технологическим процессом добычи нефти из ма-
лодебитных скважин на основе динамических моделей: автореф. дис. ... д-ра техн. наук /К.Ф. Тагирова. - Уфа: ГОУ ВПО Уфимский гос. авиационный техн. ун-т, 2008. - 32 с.
3. Светлакова, С. В. Информационно-
измерительная система динамометрирования скважин, оборудованных штанговыми глубинными насосами: автореф. дис. . канд. техн. наук / С.В. Светлакова. - Уфа: Уфимский государственный нефтяной технический университет, 2008. -16 с.
4. Новиков, Л. В. Основы вейвлет-анализа сигналов: учеб. пособие /Л.В. Новиков. - СПб.: ИАнП РАН, 1999. - 152 с.
5. Петухов, А.П. Введение в теорию базисов всплесков: учеб. пособие / А.П. Петухов. - СПб: Изд-во СПбГТУ, 1999. - 132 с.
6. Давыдов, А.В. Лекции. Тема 19. Основы вейвлет-преобразования сигналов /А.В. Давыдов. -М!М!М!.ргойт>.пагой. ги.
Поступила в редакцию 5 сентября 2012 г.
