Бесплатный автореферат и диссертация по наукам о земле на тему

Разработка метода локализации остаточных запасов нефти на поздних стадиях разработки

ВАК РФ 25.00.17, Разработка и эксплуатация нефтяных и газовых месторождений

Автореферат диссертации по теме "Разработка метода локализации остаточных запасов нефти на поздних стадиях разработки"

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

БАТАЛОВ ДМИТРИЙ АЛЕКСАНДРОВИЧ

РАЗРАБОТКА МЕТОДА ЛОКАЛИЗАЦИИ ОСТАТОЧНЫХ ЗАПАСОВ НЕФТИ НА ПОЗДНИХ СТАДИЯХ РАЗРАБОТКИ

Специальность 25.00.17 - Разработка и эксплуатация нефтяных и газовых месторождений

3 О С1М 23:5

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

Тюмень-2015

005562800

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего образования «Тюменский государственный нефтегазовый университет» (ТюмГНГУ) на кафедре «Разработка и эксплуатация нефтяных и газовых месторождений»

Научный руководитель - доктор технических наук

Стрекалов Александр Владимирович

Официальные оппоненты: - Мулявин Семен Федорович, доктор технических наук, Акционерное общество «Сибирский научно-исследовательский институт нефтяной промышленности», заведующий отделом проектирования и анализа разработки месторождений;

- Пуртова Инна Петровна, кандидат технических наук, Федеральное государственное унитарное предприятие «Западно-Сибирский научно-исследовательский институт геологии и геофизики», заместитель генерального директора по подсчету запасов, мониторингу и разработке месторождений углеводородов.

Ведущая организация - Общество с ограниченной ответственностью

«ТюменНИИгипрогаз» (ООО «ТюменНИИгипрогаз»)

Защита состоится 22 октября 2015 года в 9.00 часов на заседании диссертационного совета Д.212.273.01 на базе ТюмГНГУ по адресу: 625027, г. Тюмень, ул. 50 лет Октября, 38.

С диссертацией можно ознакомиться в библиотечно-издательском комплексе ТюмГНГУ по адресу: 625027, г. Тюмень, ул. Мельникайте, 72 а, каб. 32.

Автореферат разослан 22 сентября 2015 года.

Ученый секретарь

диссертационного совета,

кандидат технических наук, доцент

Аксенова Наталья Александровна

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

Актуальность работы

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

Степень разработанности темы исследования

Проблемам определения расположения слабодренируемых и застойных зон нефтяных залежей, поиска остаточных запасов и обоснование применения технологий вовлечения этих зон в процесс разработки посвящены работы многих отечественных и зарубежных специалистов: Крылова А.П., Абасова М.Т., Батурина Ю.Е., Боксермана A.A., Басниева К.С., Баишева Б.Т., Борисова Ю.П., Вахитова Г.Г., Горбунова А.Г., Гавуры В.Е., Гаттенбергера Ю.П., Давыдова A.B., Желтова Ю.В., Жданова С.А., Закирова С.Н., Леви Б.И., Максутова P.A., Максимова М.М., Мирзаджанзаде А.Х., Муслимова Р.Х.,

Малютиной Г.С., Маслянцева Ю.В., Михайлова H.H., Николаевского В.Н., Розенберга М.Д., Сургучева М.Л., Сазонова Б.Ф., Степановой Г.С., Стрижова И.Н., Симкина Э.М., Фурсова А.Я., Чижовой В.А., Цынковой О.Э., Черницкого A.B., Шахвердиева А.Х., Шавалиева A.M., Шалимова Б.В., Шарбатовой И.Н., Щелкачева В.П., Юрьева А. Н. и других.

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

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

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

Объект и предмет исследования

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

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

1. Анализ известных прямых и косвенных методов исследования продуктивных пластов на предмет нахождения расположения остаточных запасов нефти.

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

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

3. Апробация разработанного метода на гидродинамическом симуляторе и промысловых данных пласта БС]0 Савуйского месторождения.

Научная новизна выполненной работы

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

2. Доказано, что разработанный метод позволяет достичь приемлемую точность (выше 85 %) при условии накопления промыслового опыта процесса разработки более 60 % от предельного времени или 82 % от конечного коэффициента извлечения нефти.

Теоретическая значимость работы

1. Изложены элементы теории подземной гидромеханики в области решения её обратных задач для определения зон локализации остаточных запасов.

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

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

Практическая значимость работы

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

конец исследуемого периода разработки пласта применяются при проектировании ГТМ в ОАО «Сургутнефтегаз».

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

Методология и методы исследования

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

Положения, выносимые на защиту

1. Метод, основанный на поиске минимума целевой п-мерной функции относительных фазовых проницаемостей в зонах между скважинами, построенных в виде топологической схемы, включающей источники и приемники гидравлического импульса раздельно по каждой фазе.

2. Метод аппроксимации поверхностей нефтенасыщенности, основанный на отборе точек границами допуска и функции Бакли-Леверетга в качестве базовой для интерполяции.

Соответствие диссертации паспорту научной специальности

Область исследования соответствует паспорту специальности 25.00.17 - «Разработка и эксплуатация нефтяных и газовых месторождений», а именно: пункту 3 «Научные аспекты и средства обеспечения системного комплексного (мультидисциплинарного) проектирования и мониторинга процессов разработки месторождений углеводородов, эксплуатации подземных хранилищ газа, создаваемых в истощенных месторождениях и водонасыщенных пластах с целью рационального недропользования».

Степень достоверности результатов работы

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

Апробация результатов исследований

Результаты диссертационной работы и ее основные положения докладывались и обсуждались на: докладывались на XXVI-й научно-технической конференции молодых ученых и специалистов ОАО «Сургутнефтегаз» в 2006 г; XXVlI-й научно-технической конференции молодых ученых и специалистов ОАО «Сургутнефтегаз» в 2007 г.; VII-й конференции молодых специалистов организаций, осуществляющих виды деятельности, связанной с пользованием участками недр на территории Ханты-Мансийского автономного округа-Югры в 2007 г.; XXVIII-й научно-технической конференции молодых ученых и специалистов ОАО «Сургутнефтегаз» в 2008г; XXXIII-й научно-технической конференции молодых ученых и специалистов ОАО «Сургутнефтегаз», 2013 г.; XXXIV научно-технической конференции молодых ученых и специалистов ОАО

«Сургутнефтегаз» в 2014 г.; XIV конференции молодых специалистов, работающих в организациях, осуществляющих деятельность, связанную с использованием участков недр на территории Ханты-Мансийского автономного округа - Югры в 2014 г.

Публикации

Результаты выполненных исследований отражены в 8 печатных работах, в том числе 4 работы в изданиях, рекомендованных ВАК РФ.

Объем и структура работы

Диссертационная работа изложена на 165 страницах машинописного текста, содержит 7 таблиц, 87 рисунков. Состоит из введения, трех разделов, основных выводов и рекомендаций, списка использованных источников из 86 наименований.

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

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

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

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

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

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

Во втором разделе исследованы теоретические предпосылки вычислительной идентификации полей нефтенасыщенности.

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

1. координаты точек продуктивных забоев для скважин (возможно включая нагнетательные) -

2. результаты исследований кривой ОФП по нефти - к„ = <р„(сгв), для группы кернового материала, отобранного из продуктивного пласта;

3. мгновенный стационарный дебит нефти и воды - д" , где / - индекс скважины;

4. мгновенная депрессия, создаваемая в добывающих скважинах — .

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

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

Для скважины /', коэффициенты продуктивности по нефти и жидкости будут описываться согласно классическим отношениям:

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

где (х - величина динамической вязкости среды, Пас; к - фильтруемая (проницаемая) толщина пласта, м; г - приведенный или просто радиус скважины; к - величина абсолютной проницаемости, м2; Я - радиус контура дренирования/питания, м; 5 - коэффициент символизирующий «несовершенство» скважин или скин-эффект.

К =

2 лкИ

(2)

0,0

н

1

ж

0, м3/с

АР, Па

Ар

ИК по жидкости

Рисунок 1 - Пример построения и обработки индикаторной линии по одному замеренному режиму Соответственно, коэффициент продуктивности для нефти

и для жидкости

Величина вязкости среды в (4)

ч"

К'

2лк'И

Л

1п| 5. |+с

7+М-

(4)

(5)

ч'+ч' ч"+ч'

Приравняв (3) и (4) к рассчитанным по фактическим данным

продуктивностям получим

1 глк"^

1п| — 1 + 5,

1

п

2 пк?к

tg(a) Ар

(6)

(7)

Таким образом, фазовые проницаемости по нефти и среднюю по жидкости можно выразить из (6) и (7), чтобы получить для любой скважины ;

к" =■

1п| ~ I + 5,

М,

Ар2пИ%

(8)

Отсюда величина ОФП определяется отношением (например, для нефти)

(9)

4-Х-, к*

(10)

Так как функции ОФП известны, то возможно отыскать значение водонысыщенности (сг,), соотнесенное рассчитанному из (10) значению ОФП по нефти (рис. 2).

Поиск искомой водонасыщенности на основе известной функции ОФП нефти требует решения уравнения

(Н)

Отсюда, значение нефтенасыщенности о„= 1 -ов при незначительном газовом факторе. Разработанный метод позволяет вычислить зону текущей

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

к"

1.0

<РЛ<?.)

К

..........] [Ч

0,0

к' 1.0

ЭГ

1.0 а.

Рисунок 2 - Графическое решение для поиска текущих нефте- и водонасыщенности на основе кривой ОФП Аналогично (3) и (4) получено, что

Ар'

к-=-

2лк'М, = д.

Отсюда депрессия м'

■а-

2лк"И1 — 2лк;И,

и соотношение фазовых проницаемостей нефти и воды

= ^ г = 0»М, = 3

(12)

(13)

(14)

(15)

к" к' к' <2вц. Вместе с тем через кривые ОФП фазовые проницаемости для /-й скважины

к" = •«'.к), К =к'к, = к, ■<?.{&,) (16)

также позволяют отразить соотношение (15) через функции ОФП

Построив вспомогательную функцию У(о"в)и решив уравнение

Г(ст.) = д, (18)

относительно а„ задача также решается без использования вязкости жидкости.

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

Изменение давления в любой точке пласта можно рассчитать по известному автомодельному решению в условиях упруго-водонапорного режима

= (19)

где Ар„(А/) - изменение пластового давления в точке через время д<; е — величина гидропроводности, м2/с; г{ - дистанция от /-ой точки импульса до точки наблюдения - т; Дд, - изменение притока жидкости за время ы\х ~ пьезопроводность, отражающая скорость передачи давления в пористой среде через жидкость; Л^-количество скважин.

Величина пьезопроводности определяется к"

Х" = ^[тРи+(\-т)рХ (20)

где Д,- коэффициент сжимаемости нефти; /?„- коэффициент сжимаемости

скелета; т - коэффициент открытой пористости.

Для учета изменений в создаваемых скважинами импульсах введем дискретную функцию притока жидкости от времени для любой скважины / -<?,(')> причем <7,(')о. когда происходит отбор флюида, а (')>(), когда идет нагнетание флюида в исследуемый продуктивный пласт. Сам импульс за условный квант времени - Д/ (минимально прослеживаемое время, например 1 сутки = 86400 с) будет

!,(/) = (0-9,С-ДО, если / =0, то /,(0 = <?,('). (21)

Время квантуется величиной Аt от нуля (начало исследуемого периода) и до конца исследуемого интервала.

Любой гидравлический импульс прослеживается в точке наблюдения, в которой будет скважина -ус функцией забойного давления - p,(t). Импульс будет являться изменением давления к определенному и дискретному времени t

UJ(I) = PjO)-PJOO), (22)

где t0 - время начала исследуемого интервала, с.

В результате для любой скважины - у", в которой прослеживается импульс от скважины-источника - i изменение, действующее на забойное давление скважины -у" и относящееся по физическому смыслу только к скважине-источнику - / будет

(23)

где I, - количество импульсов скважины - г,к - индекс режима к е [О, Ь -1]; - гидропроводность на линии от скважины - /' к скважине - у; г0^л -расстояние между скважинами - / и

Так как в зоне интерферирующего гидродинамического взаимодействия скважины -j находится более одной скважины, включая и ее саму, то согласно суперпозиции следует суммировать импульсы от всех скважин принадлежащих множеству -Я}, соответствующему скважине - у , в которой наблюдается изменение давления согласно условию: скважина - / принадлежит множеству К, ~ если г(,„у)<гтах, где гтах - максимальный радиус влияния, за

пределами которого интерференция не учитывается. Скважина у" безусловно входит в множество Л, первой.

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

й

(24)

Предполагая, что выражение (24) будет относится к каждой фазе (флюиду) отдельно, введем обозначение для воды - (в) и нефти - (н), тогда "парциальные" изменения давления будут:

и

Г с>=1£

/.'"с.)

(к,

£7

1«к, к, =0 4

(25)

(26)

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

УМ,+(!-")

¿,-1 /(»¡Г/ Л,/"'

"у (о=иум+иум-

4**£ЛА

Е1

+

У (27)

И, + к

где А = -'-у-^ - средняя эффективная толщина пласта в зоне между скважинами,

м (берется по геологической модели); - фазовая проницаемость по воде в зоне между 7-ой и ¿-ой скважинами, м2.

Расстояние до самой скважины - у принимается равным радиусу

скважины - гс.

Поиск неизвестных фазовых проницаемостей и в количестве

определенным множеством Л/ опишем классической задачой оптимизации -минимизации функции для каждой ]-й скважины, в виде суммы относительных отклонений между расчетным и фактическим изменением давлений в скважине

к-О

О /«.<•) Л.(«) £<•) Ь(») Ь<») 1гМ V

(я,2-м>'

^ -л' -л' ->л'—^ т1П > (28)

где ¿у — количество режимов у-ой скважины; и®(/) - функция фактического изменения давления на забое скважины -у.

Принятые в работе постулаты:

1) за некоторый достаточно малый интервал времени - Т распределение насыщенности не претерпевает значительных изменений;

2) интервал времени Т достаточен для прохождения волны изменения давления между скважинами;

3) фильтрация на участках между скважинами в отдаленности более 20% расстояния между скважинами является плоскопараллельной и гидродинамически объединяет две фазы (вода+нефть) в упругоемкости системы в соответствии с текущими насыщенностями.

Это позволило после решения задачи (28) через полученные фазовые проницаемости определить текущее значение насыщенности по всем выделенным на схеме согласно множеству Л/ линиям — (/'-> у) через функцию -Г {т.).

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

Для более достоверного отражения полей насыщенности в работе предложен метод аппроксимации, основанный на прототипе «Гептоид».

Известны - п опорных точек в трехмерном пространстве - у^ о>). Требуется вычислить координату <%,) точек поверхности нефтенасыщенности лежащих в области определения [хтт, дгт(Ц], [утш, утах] согласно шагам сетки

= + = (29)

Ох-\

= У, = Ау-г + Ут1п,1 = 0..СУ -1, (30)

и у — 1

где г — индекс точки по оси У;у — индекс точки по оси X.

Каждая из опорных точек (**, Ук, о*) проецируются на искомую - (*,, ур сг^,)), определим расстояние между ними в плоскости Х-У

-(у,-Л?. (31)

Для каждой точки (х,, ур можно составить множество М=[йо, Лк,... •/?„_]], в котором помещается расстояние /?* от каждой опорной точки до искомой в плоскости Х-У. После сортировки по росту расстояния из М выделяется множество Е, взятия первых N элементов. По-сути N - количество точек влияния. Для предлагаемого метода Стрекалова-Баталова - N=4.

Максимальное расстояние от опорных точек, соответствующих выборке -Е, до искомой будет величина

E^=max(E)=EN^. (32)

Нормируя элементы множества Е по максимальному расстоянию Е^ получим множество нормированных расстояний

(33)

шах

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

(34)

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

Коэффициент влияния опорной точки

(35)

где J - ,1-функция Бакли-Леверетта, которая используется в качестве базовой функции интерполяции.

Суммированием всех Рк получим коэффициент сглаживания

(36)

/•о

на основании которого нормируем полученные коэффициенты (35)

Искомая координата вычисляется, как сумма произведений -Р'к

на, соответствующие им координаты ак опорных Лоточек:

N-1

Последовательность вычислений выглядит следующим образом.

1. Для каждой искомой точки (у) отсеивается множество, состоящее из N опорных точек с минимальными расстояниями в плоскости Х-¥.

2. Находятся коэффициенты, выражающие степень влияния каждой опорной точки на искомую, в зависимости от расстояния (33-35).

3. Вычисляется значение насыщенности в искомой точке (38). Перебирая величины Х01 получим искомую поверхность г](х,у).

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

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

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

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

(38)

привязки к геометрии вскрытия и определенным элементам системы разработки.

Из построенной карты текущей нефтенасыщенности через 1900 сут (рис. 3) следует, что вытеснение нефти водой происходит по разному: в районе 0019, 0012, 005, 0022, 0018, 0015 вытеснение равномерно с небольшим островком нефти, а по линии 009, 006, 004, ООП, 002, 000 крайне неравномерно с образованием множества "целиков нефти". Установлена неподвижность остаточных запасов нефти. Наличие приливов на прямоугольной границе ГДМ вызвано оттеснением запасов нефти законтурной водой, поступающей согласно принятым граничным условиям.

500 00 1000.00 ' ■1 < ч 000 500.00 1000 00 1500 00

Рисунок 3 - Поле нефтенасыщенности ГДМ через 1900 сут с начала разработки

Доказательством работоспособности разработанного метода является физика волнового процесса вытеснения. Период волн насыщенности существенно больше волн давления (рис. 4), что позволяет косвенно через кривую давления исследовать среднюю насыщенность за исследуемый период. Период относительно слабо изменяющейся насыщенности согласно их динамике по данным ГДМ составляет 1-3 месяцев. При этом кривая давления за этот период успевает изменить свою производную не менее 5-10 раз. Это позволяет ее полноценно использовать для решения обратной задачи (раздел 2).

Для анализа выбран интервала времени: выбран 10.08.2002-06.02.2003. Дата начала моделирования ГДМ 01.01.2000. Данный интервал времени разделен на три периода по 60 сут.

(пх= 16;п\=25;яЬ=0; пласт:ВС10)

Рисунок | - Динамика некоторых ФЕС в зоне около 009 скв

Первым этапом идентификации текущих насыщенностей является поиск насыщенностей в стационарной постановке на основе индикаторных линий и кривых ОФП.

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

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

На рис. 6 показана кривая зависимости погрешности от относительного времени опыта разработки (Т, д.ед.).

1 ^у» 11 21 ¿1 41 я

¡4-

—*— факт -•— Теория 3

оГ

КнО, Кн1, Кн2, КнЗ. Кн4, Кн5, КвО, Кв1, Кв2, КвЗ, Кв4, Кв5,

мД мД мД мД мД мД мД мД мД мД мД мД

49.521 - 335.04 - - - 319.51 772.12 - 52.627 - 41.7

Рисунок 5 - Сопоставление расчетной и фактической кривой изменения забойного давления (период I) 60 сут от начала анализа (10.08.2002)

Рисунок 6 - Зависимость относительной погрешности метода локализации от относительного времени опыта разработки

Погрешность результатов по методу Стрекалова-Баталова экспоненциально падает по мере накопления опыта (в виде режимов работы скважин) разработки. Зона минимально необходимой достоверности (погрешность менее 25 %) достигается после 46 % проектного времени разработки и 82 % выработки извлекаемых запасов.

На рис. 7 приведена автоматически построенная схема интерференционных связей. Анализируется взаимодействие скважин, поэтому число связей составляет 2312. Не работающие скважины и скважины без изменения режима условно назовем простаивающими, т.е. не имеющими импульсов. Рассмотрим область пласта БСю вокруг скважины 621 (рис. 8). В зоне скважин 619-616-224 существует довольно большой целик нефти, при том, что эти скважины были остановлены, из-за большой обводненности. Однако такого целика нет в радиусе 1500 м от 712-й скважины. Такой же по размерам целик находится между скважинами 415 и 615. Причем на более раннем периоде данный целик еще не был сформирован (рис. 9) также как и зоне скважин 728-411.

Рисунок 7 - Интерференционная схема анализа БС10 Савуйского месторождения (копия из ПО «Немезида») на 01.06.2011

Рисунок 8 - Фрагмент карты нефтенасыщенности пласта БС10 Савуйского месторождения на 01.04.2012 пятиточечное заполнение интерференционных связей (с топологической схемой исследованной интерференции)

Рисунок 9 - Фрагмент карты нефтенасыщенности пласта БС10 Савуйского месторождения на 01.09.2011 пятиточечное заполнение интерференционных связей (с топологической схемой исследованной интерференции)

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

Использование метода для построения карт остаточной нефтенасыщенности объекта БСю Савуйского месторождения позволило разработать комплексную программу выработки остаточных запасов нефти:

1) обосновано строительство боковых стволов на 4 скважинах;

2) использование на 30 скважинах химических методов увеличения нефтеотдачи;

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

Ее внедрение позволило обеспечить дополнительную добычу нефти в размере более 30 тыс.т нефти с продолжающимся эффектом.

ОСНОВНЫЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ

1. Научно обосновано применение метода локализации остаточных запасов нефти на поздних стадиях освоения месторождения, разработка которого основана на принципах гидродинамических исследований и решения обратных задач подземной гидромеханики с применением прототипа, предложенного И.П. Пуртовой в 2007 г.

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

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

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

Стрекалова-Баталова для локализации насыщенностей выявил его работоспособность с погрешностью в пределах 8-15 % на период разработки более 60 % от всего времени освоения месторождения или 82 % от конечного коэффициента извлечения нефти.

5.Выявлены ограничения применения метода. Причиной неработоспособности (вероятность менее 7 %) метода, является отсутствие выраженных (с амплитудой более 6-10 %) колебаний дебитов жидкости. В случае несущественного отличия в сжимаемости и вязкости нефти и воды разделение импульсов по фазам невозможно.

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

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

В изданиях, рекомендованных ВАК РФ

1. Баталов Д.А. Технология интенсификации добычи нефти на водоплавающих залежах месторождений Западной Сибири / Д.А. Баталов, А.Т. Хусаинов // Оборудование и технологии для нефтегазового комплекса. - 2013.-№4. -С. 56-58.

2. Баталов Д.А. Методика прогнозирования доизвлечения остаточных запасов на водоплавающих нефтяных залежах находящихся на поздних стадиях разработки / Д.А. Баталов, А.Т. Хусаинов // Электронный научный журнал «Нефтегазовое дело». - 2013. - № 2. -С. 161-167. URL:http://www.ogbus. ru/authors/ BatalovDA/BatalovDA _l.pdf.

3. Стрекапов A.B. Проблемы формирования геологических моделей / A.B. Стрекалов, Д.А. Баталов, А.Т. Хусаинов // Электронный научный журнал «Нефтегазовое дело». - 2014. - № 1. -С. 65-93. URL:http://www.ogbus.ru/ authors/BatalovDA/BatalovDA_2.pdf.

4. Стрекалов A.B. Общие проблемы теории фильтрации при совместном движении двух флюидов / A.B. Стрекалов, Д.А. Баталов // Электронный научный журнал «Нефтегазовое дело». - 2014. - №5. - С. 91-103. http://ogbus.ru/article/obshhie- problemy-teorii-filtracii-pri-sovmestnom-dvizhenii-dvux-flyuidov.

В других изданиях

5. Баталов Д.А. Оптимизация системы разработки пласта ЮС2 Русскинского месторождения // Сб. докл. победителей XXVI научн.-техн. конф. молодых ученых и специалистов ОАО «Сургутнефтегаз», 2006 г. - С. 55 - 58.

6. Баталов Д.А. Увеличение эффективности разработки водоплавающих залежей на примере пласта БС12 Родникового месторождения // Сб. докл. победителей XXVII научн.-техн. конф. молодых ученых и специалистов ОАО «Сургутнефтегаз», 2007 г. - С. 31 - 35.

7. Баталов Д.А. Увеличение эффективности разработки высокопродуктивных пластов на поздней стадии методами деформационно-гидродинамического воздействия на примере пласта БС10 Савуйского месторождения // Сб. докл. победителей XXVIII научн.-техн. конф. молодых ученых и специалистов ОАО «Сургутнефтегаз», 2008г. - С. 32 - 37.

8. Баталов Д.А. Повышение рентабельности разработки нефтяных месторождений на поздней стадии // Комитет по молодежной политике ХМАО, сборник докладов победителей окружного молодежного конкурса «Золотое Будущее Югры 2008 г.» - С. 27 - 32.

Соискатель

Д.А. Баталов

Издательство «Вектор Бук» Подписано в печать 20.08.2015 г. Формат 60x84/16. Бумага офсетная. Печать Riso. Усл. печ. л. 1,44. Тираж 100 экз. Заказ215. Отпечатано с готового набора в типографии издательства «Вектор Бук». 625004, г. Тюмень, ул. Володарского, 45. Тел. (3452) 46-54-04,46-90-03.