На правах рукописи

ЕНИКЕЕВ РУСЛАН РИНАТОВИЧ ПОМЕХОУСТОЙЧИВЫЕ АЛГОРИТМЫ ОБРАБОТКИ ДАННЫХ ПРОМЫСЛОВЫХ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН Специальность 25.00.17 - Разработка и эксплуатация

нефтяных и газовых месторождений

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Уфа - 2004

Работа выполнена в Уфимском филиале ООО Нефтеюганский научно исследовательский и проектный институт нефти.

Научный консультант: доктор физико-математических наук, профессор Булгакова Гузель Талгатовна

Официальные оппоненты: доктор технических наук, профессор Пономарев Александр Иосифович кандидат технических наук, доцент Левченко Владимир Сидорович Ведущая организация Научно - производственное объединение Нефтегазтехнология

Защита состоится л22 октября 2004 года в 11.30 на заседании диссертационного совета Д 212.289.04 при Уфимском государственном нефтяном техническом университете по адресу: 450062, Республика Башкортостан, г.Уфа, ул. Космонавтов, 1.

С диссертацией можно ознакомиться в библиотеке Уфимского государственного нефтяного технического университета.

Автореферат разослан л 2004 года

Ученый секретарь диссертационного совета Матвеев Ю. Г.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Проектирование и осуществление контроля над процессом разработки нефтяных месторождений неотъемлемо связаны с определением коллекторских свойств и изучением фильтрационных параметров продуктивных коллекторов.

Расчет фильтрационных параметров коллектора лучше всего проводить по данным гидродинамических исследований скважин и пластов (ГДИ), которые отражают непосредственный процесс фильтрации жидкости в пластовых условиях и позволяют получить усредненную информацию о значительной части пласта.

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

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

В результате научно-технического прогресса в последние годы появились приборы нового поколения, характеризующиеся высокой чувствительностью (по давлению до 0.001 МПа, предел измерения 100 МПа, по температуре до 1000С), с электронной памятью, малогабаритные, с возможностью снятия кривых длительностью от 30 суток и более.

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

В частности, неустойчивость проявляется в условиях малых выборок, когда в координатах метода удается спрямить только небольшой участок КВД.

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

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

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

Цель работы. Создание экспертной системы интерпретации результатов промысловых исследований скважин на основе модификации классических методик с учетом априорной информации об исследуемом объекте (скважине, пласте) методами теории некорректно поставленных задач.

Усовершенствование методики интерпретации и технологии проведения промысловых гидродинамических исследований скважин и пластов.

Основные задачи

исследования 1. Разработка помехоустойчивого алгоритма интерпретации результатов гидродинамических исследований скважин и пластов на неустановившихся режимах фильтрации и его апробация.

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

3. Создание экспертной системы интерпретации промысловых гидродинамических исследований скважин и пластов.

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

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

Научная новизна 1. Разработан помехоустойчивый метод интерпретации результатов гидродинамических исследований скважин на неустановившемся режиме фильтрации с использованием теории А.М.Тихонова (теории некорректных задач).

2. Предложена методика обработки недовосстановленных кривых восстановления давления, основанная на регуляризации производной давления.

3. Разработана методика интерпретации результатов исследований скважин методом гидропрослушивания осложненных совмещением с КВД.

4. Предложен метод учета переходных режимов скважины для интерпретации результатов гидродинамических исследований скважин на установившихся режимах.

На основе разработанных подходов и методик построены соответствующие алгоритмы и создана экспертная система интерпретации исследований скважин (Свидетельство № 2001611144 об официальной регистрации программы "Обработка данных гидродинамических исследований скважин").

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

Предложена методика интерпретации недовосстановленных кривых при помощи производной давления.

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

Предложена методика учета переходных процессов, происходящих при проведении гидродинамических исследований скважин на установившемся режиме фильтрации. Учет переходных процессов позволяет оценить дополнительные фильтрационные параметры пласта (скин-фактор, пьезопроводность, проницаемость), которые согласуются с результатами других методов.

Разработанная программа интерпретации результатов гидродинамических исследований скважин с использованием помехоустойчивых алгоритмов (Свидетельство № 2001611144 об официальной регистрации программы "Обработка данных гидродинамических исследований скважин" (КВД)) используется в ОАО Самаранефтегаз, ОАО Юганскнефтегаз, ОАО Томскнефть-ВНК и ООО Западно-Малобалыкское при составлении проектов разработки нефтяных месторождений.

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

Апробация работы. Основные положения и результаты диссертацион ной работы докладывались на 47-й научно-технической конференции студентов, аспирантов и молодых ученных Уфимского государственного нефтяного технического университета (г.Уфа, 1996г.), на VI-й научно технической конференции молодых ученых и специалистов ДООО БашНИПИнефть (Уфа, 2001), на IX Всероссийской школе-коллоквиуме по стохастическим методам (Ростов-на-Дону, 2002), III-м Всероссийском симпозиуме по промышленной и прикладной математике (Ростов-на-Дону, 2002).

Публикации. По материалам диссертации автором опубликовано печатных работ, в том числе 7 статей, тезисы трех докладов и получено свидетельство об официальной регистрации программы.

ичный вклад автора в совместных публикациях:

1. [1].

Еникеев Р. М. - промысловые исследования скважин.

Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Калиновский Ю. В. - постановка задачи, обсуждение результатов.

Каримов М. У-Г. - промысловые исследования скважин.

2. [3].

Еникеев Р. М. - промысловые исследования скважин.

Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Калиновский Ю. В. - постановка задачи, обсуждение результатов.

Каримов М. У-Г. - промысловые исследования скважин.

3. [4].

Галлеев Р. Р. - постановка задачи, обсуждение результатов.

Ганиев Э. Р. - постановка задачи, обсуждение результатов.

Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Каримов М. У-Г. - промысловые исследования скважин.

4. [5].

Еникеев Р. М. Баширов И. Р., Сайфутдинов Ф. Х., Шавалеев А. М.

- промысловые исследования скважин, обсуждение результатов.

Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Каримов М. У-Г. - проведение расчетов 5. [7].

Гарифуллин Р. Н. - графическое представление результатов Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Хасанов М. М. - постановка задачи, обсуждение результатов.

6. [10].

Булгакова Г. Т. - постановка задачи, обсуждение результатов.

Еникеев Р. Р. - постановка задачи, разработка алгоритма, получение и анализ результатов.

Структура и объем работы Диссертационная работа состоит из введения, пяти глав, выводов, списка литературы, состоящего из 129 наименований. Текст работы изложен на 165 страницах, включая 68 рисунков и 10 таблиц.

Автор выражает глубокую благодарность и искреннюю признательность научному руководителю доктору физико-математических наук Булгаковой Г.Т.

за постоянное внимание и помощь в работе.

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

Автор искренне благодарит заведующего лабораторией ГДИ ДООО БашНИПИнефть кандидата технических наук Еникеева Р.М., доцента кафедры РЭГ и ГКМ УГНТУ, кандидата технических наук Калиновского Ю.В.

за поддержку и сотрудничество при выполнении настоящей работы, сотрудников лабораторий разработки и моделирования УФ ООО ЮганскНИПИнефть, ОАО Юганскнефтегаз, Самаранефтегаз, Томскнефть за поддержку и квалифицированные советы при создании и внедрении комплекса программ.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе приведен литературный обзор работ по определению фильтрационных параметров однородных пластов, таких авторов, как Аметов И.М., Баренблатт Г.И., Блинов А.Ф., Богачев Б.А., Борисов Ю.П., Бузинов С.Н., Бурдэ Д., Васильевский В.Н., Горбатова А.Н., Дияшев Р.Н., Желтов Е.П., Исайчев В.В., Каменецкий С.Г., Кульпин Л.Г., Крылов А.П., Ли Юн Шан, Литвинов А.А,, Максимов В.А., Муравьев И.М., Мирзаджанзаде А.Х., Мясников Ю.А., Пикуза В.И., Требин Ф.А., Умрихин И.Д., Хуань Коужень, Чарный И.А., Чекалюк Э.Б., Чернов В.С., Шагиев Р.Г., Щелкачев В.Н. и др.

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

Помимо стандартных методик интерпретации также рассмотрена методика интерпретации КВД при помощи производной давления.

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

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

Наиболее простым способом обработки кривых восстановления давления (КВД) является метод Миллера, Дайеса, Хатчинсона (касательной), в соответствии с которым зависимость P = f (t) перестраивается в координатах P = f (ln(t)), спрямляется конечный участок и по уравнению прямой определяются фильтрационные характеристики пласта.

На рис.1 изображена фактическая КВД. Как видно из рисунка, фактическая кривая не является гладкой, а имеет некоторую погрешность.

22. 22. 22. 22. 10 20 30 40 50 60 70 80 t, час Рис. 1. Кривая восстановления давления Для интерпретации исходная кривая перестраивается в координатах P = ln(t). При попытке спрямления конечного участка (отрезок AB см рис.2) возникает задача минимизации функционала n I(a,b) = (a * xi + b - yi )2, (1) n i= который представляет собой квадрат среднеквадратичного отклонения прямой y = a x + b от точек спрямляемого участка. Коэффициенты уравнения находятся методом наименьших квадратов. Решение системы a*,b* имеет вид:

xy - x y x2 y - xxy a* =, b* =, (2) 2 x2 - x x2 - x P, МПа знаменатели этих формул являются квадратом среднеквадратичного отклонения xi, величина этого отклонения мала. Из вышесказанного будет следовать, что при малом возмущении исходных данных отклонения a и b от их невозмущенного значения могут быть большими. В этом проявляется неустойчивость данной задачи. Это легко видеть из рис.3, на котором представлен разброс спрямляемых кривых при разбросе начальных данных в пределах статистической погрешности (в данном случае порядка 1%). Разброс a и b вызывает такой же разброс в результатах интерпретации КВД.

0, 0, 0, 0, ln t Рис. 2. Результаты эксперимента по наложению шумов (линия 1 - спрямление исходных данных, линия 2 и 3 соответственно крайние нижние и верхние прямые спрямления при наложении шумов на исходные данные в 1%) Из вышесказанного можно сделать вывод о том, что данную задачу можно отнести к некорректно поставленным задачам, которые характеризуются:

Х неединственностью решения;

Х неустойчивостью решения;

Х отсутствием классического решения.

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

Н. Тихоновым, М. М. Лаврентьевым и многими другими.

Появляется необходимость регуляризации данной задачи. Для этого предлагается использовать некоторую априорную информацию. В данном случае такой информацией может служить проницаемость пласта, определяемая геофизическими исследованиями. Это приводит к уменьшению разброса величины а, то есть к новому функционалу (a,b) = (a - a0 )2 где а0 - определяется из проницаемости пласта, найденной из геофизических исследований. В общем случае величина а и а0 не совпадают. Следовательно, P, МПа задача минимизации и функционала I(a,b) и функционала (a,b) не имеет решения, она является двухкритериальной задачей. Ее предлагается свести к однокритериальной задаче минимизации функционала M (a,b) = I(a,b) +(a,b) > 0. (3) Ее решения находятся из системы M a = (4) M = 0, b решения которой имеют вид:

xy + a0 - x y (x2 + )y - x(xy + a0 ) a** =, b** =. (5) 2 + x2 - x + x2 - x В этих формулах знаменатель становится порядка единицы при не очень малых, а решения задачи становится устойчивой. В задаче остался произвол в выборе величины. Эта величина в данной работе выбиралась из следующих соображений: вычислялось значение I(a*,b*), в качестве бралась максимальная величина, при которой I(a**(),b**()) отличается от I(a*,b*) не больше, чем среднестатистическая погрешность задания yi. В большинстве случаев величина оказывается больше 0.5, следовательно, при данной регуляризации задача становится устойчивой. Это легко видеть из рис.3, на котором аналогично рис.2 приведен эксперимент по наложению шумов в начальных данных.

A B 0. 0. 0. Интервал 0. ln t Рис. 3. Результаты применения регуляризации В третьей главе рассматривается задача интерпретации промысловых исследований с помощью производной давления.

Для того, чтобы продифференцировать кривую (рис.1), а именно, выборку значений давления ui в разные моменты времени, проинтерполируем ее в классе ломаных кривых по точкам (ui,ti ), то есть в качестве функции возьмем функцию вида ui+1 - ui u(t) = (t - ti ) + ui, ti t < ti+1. (6) ti+1 - ti P, МПа В качестве производной КВД возьмем значения производной этой ti+1 + ti ~ ломаной в точках ti =. Рассмотрим эту кривую на рис. 4. График производной находится справа от самой КВД.

n Теперь в выборку значений {ui}i=1 внесем относительную погрешность в n ~ 0.05 %, то есть рассмотрим новую выборку значений {ui}i=1, где ~ ui = ui (1+ i * 0.0005), а i - случайные числа от -1 до 1. Как и прежде для выборки n n- ~ ~ {ui}i=1 найдем значения {ui'}i=1. Графики обоих кривых представлены на рис.5.

А) 22, В) 0, 22, 22, 22, 0, 10 30 50 70 22, 10 30 50 70 t, час -0, t, час Рис. 4. Кривая восстановления давления (кривая А) и ее производная (кривая В) При относительной погрешности на кривую восстановления давления в n- ~ 0.05 % относительная погрешность при получении значений {ui'}i=1 была в среднем 100 - 300 %. Пример показывает неустойчивость операции дифференцирования, что приводит к некорректности задачи нахождения производной КВД.

В частности, неустойчивость проявляется в условиях малых выборок, когда в координатах метода удается спрямить только небольшой участок КВД.

Очень часто приходится обрабатывать так называемые недовосстановленные КВД, полученные в исследованиях, прерванных по каким - либо причинам. На характер кривой восстановления давления влияют многие факторы. Одним из основных является немгновенное прекращение притока нефти к забою. Для учета переменного притока q(t) необходим учет изменения дебита скважины на забое, а именно, учет изменения дебита от установившегося значения до остановки (нуля). Немгновенное изменение дебита при смене режимов работы показано на рис.6.

P, МПа P', МПа / час В) А) 22, 0, 22, 22, 0, 10 30 50 70 22, 22, 10 30 50 70 -0, t, час t, час Рис. 5. Результаты математического эксперимента по наложению шума на кривую восстановления давления в 0.05 % Рис. 6. Изменение дебита скважины при смене режимов работы скважины Зная Q1 - установившийся дебит до остановки, значения дебита в момент времени i - qi, а также зная депрессию P(i ), в преобразованных координатах i - t1 i -i- [Q1 ln(i - t1) + q1 ln +... + qi-1 ln - qi ln(i -i-1)] i -1 i -i- = x Q1 - qi (7) = P(i ) y Q1 - qi можно построить экстраполирующую прямую и, таким образом, найти значение углового коэффициента. Проводить экстраполирующую прямую 4kh = i + A можно с помощью регуляризирующего метода наименьших y x квадратов (глава 2). Регуляризация поиска такой прямой ведется как и в разделе 2.

На практике непосредственно измерить темп затухания притока с необходимой точностью чрезвычайно трудно, а в малодебитных скважинах P, МПа P', МПа / час достоверный учет притока из пласта технически невозможен. По этой причине приток определяется расчетным путем, который связан с выполнением операций дифференцирования экспериментальных значений давления.

Для расчета производной давления используется регуляризирующий метод дифференцирования функции Тихонова А.М., повышающий устойчивость этой процедуры.

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

Технология проведения исследований в большинстве случаев такова, что скважина останавливается для спуска глубинного манометра. С точки зрения полноты исследования, если скважина оборудована глубинным манометром, то выгодно не только осуществить гидропрослушивание, но и снять в этой скважине кривую восстановления давления (КВД). Понятно, что при совмещении этих двух видов исследований становится сложной обработка КВД, если вначале проводилось гидропрослушивание, и наоборот, усложняется обработка КР, если перед гидропрослушиванием снималась КВД. Однако, на наш взгляд, и в этом случае существует удовлетворительная технология исследования, которая заключается в следующем. В простаивающую эксплуатационную скважину спускается глубинный манометр;

если скважина работающая, то она останавливается, в нее спускается манометр и она выдерживается до стабилизации забойного давления. После этого скважина пускается в работу с некоторым постоянным дебитом Q, при этом снимается кривая падения забойного давления (КПД). Затем в некоторый момент времени T1 скважина останавливается и в ней снимается КВД до тех пор, пока прирост забойного давления станет незначительным, т.е. давление практически стабилизируется. Затем в некоторый момент времени Т2 с дебитом Q пускается в работу (или резко меняет режим) возмущающая скважина (как правило, нагнетательная), а в эксплуатационной скважине снимается КР.

Схематично изменение давления в эксплуатационной скважине показано на рис. 7.

Полагая, что процесс фильтрации в пласте описывается линейным уравнением пьезопроводности и используя принцип суперпозиции, легко записать зависимость изменения прироста давления во времени в реагирующей скважине для времени t > T2 :

Q r Q c P(t) = P (t)- P (0) = - - E - + c c i 4 k h 4 t 4 k h (8) 2 r R c - E - - E -.

i i 4 (t - T ) + Q 4 k h 4 (t - T ).

1 T1 T2 TN - 1 000 2 000 3 000 4 000 5 000 6 t, мин Рис.7. Изменение давления в эксплуатационной скважине при гидропрослушивании При записи (8) учтено, что дебит скважины-стока имеет знак минус, скважины-источника - плюс;

R - расстояние между нагнетательной и реагирующей скважиной.

Первые две интегральных экспоненты с большой точностью могут быть заменены известным асимптотическим представлением для малых значении аргумента:

- Ei(- x) = -C - ln(x), x<1, (9) C = 0.577216 - константа Эйлера, тогда после элементарных преобразовании получим t Q -T1 - R P(t)= ln + Q Ei - 4 (t -T 4 k h t ), t>T2. (10) Записав выражение (10) для двух моментов времени t1 и t2 (для двух замеров), разделим обе части этих выражении друг на друга и после преобразований получим t - T R 2 P(t ) ln + Q E Q - i 1 t 4 (t - T ) - P(t ) 2 2 (11) t - T R 1 ln + Q E Q - i t 4 (t - T ) = 0.

1 1 Как видно, (11) есть трансцендентное уравнение относительно пьезопроводности, легко решаемое любым численным методом. В процессе решения для интегральной экспоненты можно использовать различные разложения и асимптотические представления, в зависимости от значения аргумента. После нахождения из соотношения (10), записанного для t1 и t2, Рзаб, МПа k h можно найти соответствующие значения гидропроводностей. Наконец, k h для всех использованных пар точек значения и можно оcреднить.

Моменты времени t1 и t2 нельзя выбирать близкими, т.к. в этом случае на результаты обработки может повлиять погрешность замеров. Обычно t выбирается из первой половины КР, t2 - из второй. Практика обработки КР показала, что и при таком выборе точек при погрешности замеров 0.05 МПа k h некоторые текущие значения и сильно отличаются от средних значении. Поэтому для уменьшения влияния погрешностей предлагается другой способ обработки, который можно назвать интегральным.

Проинтегрируем равенство (10) один раз в пределах от T2 до TN, другой раз в пределах от Т2 до ТС (ТС < ТN), затем почленно разделим полученные выражения и после элементарных преобразований получим Tc -T1 T2 -T1 Tc -T IN ln -T2 ln -T1 ln + Q IEc Q Tc Tc T2 T2 -T (12) TN -T1 T2 -T1 TN -T Ic ln -T2 ln -T1 ln + Q IEN = Q TN TN T2 T2 -T где TN TC IN = dt IC = dt P P T2 T (13) TN TC R2 R IEN = Ei IEC = Ei - - 4 (t -T2) dt - - 4 (t -T2 ) T2 T Уравнение (12) также есть трансцендентное уравнение относительно и легко решается, после чего, используя проинтегрированное выражение (10), k h можно найти. Интегралы (13) подсчитываются любым методом численного интегрирования.

Что же касается КВД, то она описывается выражением (10), в котором отброшена интегральная экспонента (Т1 < t < T2), и для ее обработки можно использовать обычный метод Хорнера.

В пятой главе рассматривается методика обработки результатов гидродинамических исследований скважин с учетом переходных процессов.

Идентификационный подход к описанию переходных режимов с одного установившегося состояния в другое был предложен А. Х. Мирзаджанзаде. В данной работе этот подход получил дальнейшее развитие и обоснование.

При исследовании скважины методом установившихся отборов, при гидродинамических исследованиях на нефтегазовых месторождениях, строится индикаторная диаграмма в координатах Q = f (РЗ ).

Основное назначение индикаторной диаграммы состоит в том, чтобы по данным небольшого числа исследований предсказать добывные возможности скважины при изменении перепада давлений РПЛ - РЗ на забое скважины. В рамках модели установившейся фильтрации однофазной жидкости к центральной скважине, дренирующей однородный круговой пласт, индикаторная диаграмма аппроксимируется линейной зависимостью, которую можно классифицировать как дедуктивную (детерминированную) модель.

Q = -KPЗ + A = K(РПЛ - РЗ ), (14) где Q - дебит скважины, РПЛ - пластовое давление, РЗ - забойное давление, K - коэффициент продуктивности скважины, характеризующий добывные возможности скважины. Коэффициенты A и K определяются методом наименьших квадратов из условия минимизации невязки I = (A - KPЗi - Qi), (15) i где Qi, PЗi - экспериментальные значения дебитов и соответствующих забойных давлений.

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

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

Исследования переходных режимов проводятся на скважинах, оборудованных центробежными насосами (УЭЦН) с датчиками давлений на приеме. По фактическому давлению на приеме насоса Рпр определяется забойное давление и устанавливается зависимость Q = f(Pзаб), позволяющая построить индикаторную кривую. При снижении забойного давления ниже давления насыщения можно воспользоваться поправкой Вогеля для учета разгазирования.

При пластовых давлениях выше давления насыщения, а забойного давления ниже давления насыщения дебит скважины определяется по так называемой композитной кривой (рис.8).

Q = K (P - P ), пр пл з P, при P P, з з нас (16) P = P P P з нас з з 4 + + 4, при P < P, з нас 9 P P нас нас где Qнас - дебит скважины при забойном давлении, равном давлению насыщения.

Рис. 8. Композитная интегральная кривая Вогеля То есть дебит скважины при забойных давлениях ниже давления насыщения можно рассчитать по стандартной формуле Дюпюи, если использовать эффективное забойное давление.

Учет переходных процессов позволяет оценить дополнительные фильтрационные параметры пласта (скин-фактор, пьезопроводность, проницаемость), которые согласуются с результатами других методов.

Для анализа переходных режимов скважины используется идентификационная модель. Рассмотрим элемент системы разработки, состоящий из контура нагнетания и добывающей скважины. Характеристики пласта от контура нагнетания до контура питания эксплуатационной скважины Rк, на котором определяется пластовое давление Рпл, в окрестности эксплуатационной скважины описываются уравнением Q = Kпр (Pн - Pпл ), (17) где Кпр - коэффициент приемистости нагнетательных скважин, Рн - давление на контуре нагнетания. Характеристики эксплуатационной скважины задаются уравнением Q = Kпр (Pпл - P3). (18) Графики этих уравнений в координатах Рпл, Q представлены на рис. 9.

При известном забойном давлении Р31 дебит скважины Q1 определяется по точке пересечения характеристик пласта и скважины. В идеальной безынерционной системе изменение дебита вызывает мгновенное установление забойного давления Р32 (участок АС).

Реальные процессы установления забойного давления являются нестационарными (рабочая точка переходит в состояние С через промежуточную точку В, описывая нелинейную траекторию). Изменение забойного давления от значения Р31 до значения Р32 должно описываться нестационарным (идентификационным) уравнением вида dPЗ Q T + PЗ = Pпл -, (19) dt Kпр Pз(t = 0)= P31.

Рпл Рн B A РЗ C РЗ Q Q Рис.9. Характеристики скважины и пласта Q Rk Время установления процесса Т порядка, где - коэффициент пьезопроводности. Уравнение (19) можно записать в виде Pпл - P3(t) dP3 = T /Q +, (20) Q dt Kпр т.е. в виде прямой линии в координатах х, у :

Pпл - P3(t) dP = ;

= T / Q (21) у х Q dt с угловым коэффициентом i и свободным членом А, равными i = T ;

A =. (22) Kпр Учитывая, что преобразованные кривые имеют некоторый разброс точек, прямая аппроксимируется методом наименьших квадратов. Величина пластового давления в уравнении (19) определяется из индикаторных диаграмм.

По вычисленным значениям i и А определяется коэффициент продуктивности Кпр и время установления стационарного состояния в пласте Т.

При изменении забойного давления изменяются размеры зоны влияния вокруг скважины, за пределами которой возмущение, вызываемое нарушением стационарного притока к скважине, практически не сказывается. Радиус депрессионной воронки выражается через время установления Т в виде Rк = 2 T. (23) Зная Т и Rк, определяем коэффициент пьезопроводности, через который при известных значениях вязкости нефти н и упругоёмкости пласта * определяется коэффициент проницаемости пласта k k = *. (24) Из рассчитанного значения коэффициента продуктивности определяется скин-фактор 2 k h Rк S = - ln. (25) Kпр rc На рис.10 показана корреляционная зависимость между значениями коэффициента проницаемости, рассчитанными по данным нормальной эксплуатации скважин (ДНЭ) одного из месторождений ОАО Юганскнефтегаз и с учетом переходных режимов. Коэффициент корреляции составляет 0.85.

На рис.11 представлены результаты расчета скин - фактора по формуле (25) и по данным ДНЭ. Коэффициент корреляции составляет 0.87.

Корреляция 0. 02468 10 Проницаемость с учетом переходных режимов, 10-3 мкм Рис.10. Корреляционная зависимость проницаемости Корректировка коллекторских свойств в окрестности скважин позволяет согласовать расчетные и фактические параметры работы скважин и может рассматриваться как предварительная процедура по этапу адаптации истории разработки при создании гидродинамической модели месторождения.

- Проницаемость по ДНЭ, мкм Корреляция 0. - -5 0 5 10 15 20 25 Скин-эффект с учетом переходных режимов Рис.11. Корреляционная зависимость скин-эффекта

ОСНОВНЫЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ 1. Обратная задача интерпретации кривых восстановления давления (КВД) является некорректной задачей, поэтому интерпретация КВД не должна быть ограничена простым расчетом фильтрационных параметров по готовым формулам. Показано, что обязательным этапом обработки КВД является проверка полученных результатов на устойчивость. Для повышения устойчивости данных гидродинамических исследований скважин предложен регуляризирующий (повышающий устойчивость) алгоритм, заключающийся в привлечении априорной информации о пределах, в которых могут лежать искомые параметры. Предложенный алгоритм устойчивой оценки фильтрационных параметров, основанный на регуляризации А.Н.Тихонова, можно назвать регуляризирующим методом наименьших квадратов.

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

3. Предложены сглаживающие алгоритмы, основанные на теории некорректно поставленных задач, существенно снижающие влияние шумов на производную давления. Рассмотрена возможность применения методики интерпретации некачественных (недовосстановленных) кривых восстановления давления при помощи производной с использованием регуляризирующих алгоритмов. Регуляризирующий метод наименьших квадратов позволяет Скин эффект по ДНЭ интерпретировать КВД, представленную в модифицированных координатах, учитывающих переменный дебит притока после остановки скважины.

4. Разработан метод интерпретации результатов ГДИ, основанный на совмещении кривой реагирования и КВД. Методика интерпретации является интегральным методом, сглаживающим погрешности замеров. В результате применения данного метода повышается степень достоверности определения фильтрационных характеристик коллектора по результатам гидродинамических исследований скважин.

5. Проведен анализ переходных режимов при изменении режима эксплуатации скважины. Предложена методика интерпретации результатов промысловых исследований скважин с учетом переходных режимов. При исследовании скважин на установившихся режимах моделирование переходных процессов позволяет дополнительно определить фильтрационные характеристики пласта, такие как проницаемость, скин эффект, коэффициент пьезопроводности. Корреляция значений фильтрационных характеристик пласта, рассчитанных по предлагаемой методике с данными нормальной эксплуатации скважин и ГИС, составляет порядка 90%.

6. На основе разработанного алгоритма и методик создана экспертная система, прошедшая успешное опробование, широко используемая для решения конкретных задач на месторождениях НК Юкос и получившая свидетельство об официальной регистрации (Свидетельство № об официальной регистрации программы "Обработка данных гидродинамических исследований скважин").

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ 1. Еникеев Р. М., Еникеев Р.Р., Калиновский Ю.В., Каримов М. У-Г.

Обработка кривых реагирования после снятия КВД в реагирующей скважине // Системный анализ процессов разработки нефтяных месторождений и транспорта нефти и нефтепродуктов: Тезисы докладов.

- БашНИПИнефть, 1996. - С.19.

2. Еникеев Р.Р. Новое оборудование для проведения гидродинамических исследований скважин // Материалы XXXXVII-й научно-технической конференции студентов, аспирантов и молодых учен. - Уфа: УГНТУ, 1996. -Т.1. -С.34-35.

3. Еникеев Р. М., Еникеев Р.Р., Калиновский Ю.В., Каримов М. У-Г.

Обработка осложненных кривых реагирования при гидропрослушивании // Межвузовский сборник научных статей. Вып. 1. Нефть и газ. - Уфа, 1997. - С.64-66.

4. Галеев Р. Р., Ганиев Э. Р., Еникеев Р.Р., Каримов М. У-Г. Использование приборов нового поколения для гидродинамических исследований скважин // Разработка и совершенствование методов увеличения нефтеотдачи трудноизвлекаемых запасов. Проблемы и решения.-Уфа:

Гилем, 1998.

5. Баширов И. Р., Еникеев Р. М., Еникеев Р.Р., Каримов М. У-Г., Сайфутдинов Ф. Х., Шавалеев А. М. Оценка границ распространения коллектора гидродинамическим методом // Современные инструменталь ные физико-химические и гидродинамические методы исследований пластовых флюидов, пород и продуктивных пластов. -Уфа:

БашНИПИнефть, 1999. -Вып. 97.

6. Еникеев Р. Р. Опыт применения гидродинамических исследований скважин для оценки границ распространения коллектора // Нефтепромысловое дело.- 2001.- № 5. - С.29-30.

7. Гарифуллин Р. Н., Еникеев Р. Р., Хасанов М. М. Новый подход к интерпретации кривых восстановления давления // Вестник инжинирин гового центра Юкос. - 2001.- № 2. C.13-16.

8. Еникеев Р. Р. Интерпретация результатов исследований скважин на неустановившемся режиме как некорректно поставленная задача // Молодые ученые - решению важнейших проблем нефтедобычи АНК Башнефть: Материалы VI научно-технической конференции молодых ученых и специалистов ДООО БашНИПИнефть (Уфа,, 30 мая 2001г.).

- Уфа: Башнипинефть, 2001.

9. Еникеев Р. Р. Об одном способе обработки результатов гидро прослушивания // Научно технический вестник УЮКОСФ.- 2003. -№7. - С.28-30.

10. Булгакова Г.Т., Еникеев Р.Р. Рациональная схема исследования параметров пласта и призабойной зоны // Интервал.- 2003. - № 6-7. C.75 78.

11. Свидетельство № 2001611144 об официальной регистрации программы Обработка данных гидродинамических исследований скважин (КВД)/ М.М.Хасанов, Г.Т.Булгакова, Р.Р.Еникеев и др. //Реестр программ для ЭВМ, 06.09.2001.

бесплатно Еникеев Руслан Ринатович ПОМЕХОУСТОЙЧИВЫЕ АЛГОРИТМЫ ОБРАБОТКИ ДАННЫХ ПРОМЫСЛОВЫХ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Распечатано 06 июля 2004г.

Формат 60х84/ Бумага офсетная. Печать плоская. Усл. - изд. л. 1.01. тираж 100. Заказ № Центр анализа и прогнозирования Юкос ЭП    Книги, научные публикации