Моделирование заводнения нефтяного пласта с использованием термогелей тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

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

00461

5747

НУРУЛЛИН РУСТЕМ ФАРИТОВИЧ

МОДЕЛИРОВАНИЕ ЗАВОДНЕНИЯ НЕФТЯНОГО ПЛАСТА С ИСПОЛЬЗОВАНИЕ^ ТЕРМОГЕЛЕЙ

Специальность 01.02.05 — механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

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

2 5 ноя 2010

Тюмень - 2010

004613747

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

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

Никифоров Анатолий Иванович

Официальные оппоненты:

доктор физико-математических наук, профессор Федоров Константин Михайлович

доктор физико-математических наук, доцент Конюхов Владимир Михайлович

Ведущее предприятие: Татарский научно-исследовательский и

проектный институт нефти (ТатНИПИ) ОАО «Татнефть», г.Бугульма

Защита состоится 26 ноября 2010 года, в 1600 часов 'на заседании Диссертационного совета ДМ 212.274.09 в Тюменском государственном университете, по адресу: 625003, г.Тюмень, ул. Перекопская, д. 15А.

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета по адресу: 625003, г. Тюмень, ул. Семакова, 18.

Автореферат разослан 1с октября 2010 г.

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

к.ф.-м.н., доцент

Мусакаев Н. Г.

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

Актуальность темы

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

Среди физико-химических МУН пластов особое место занимают потокоотклоняющие технологии, в которых прирост КИН достигается за счет блокирования основных водопроводящих путей и вовлечения в активную разработку неподвижных и малоподвижных запасов нефти. К таким методам относятся технологии воздействия на пласты гелеобразующими и полимердисперсными системами.

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

Цель работы

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

Задачи и методы исследования

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

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

Моделирование процессов блокирования поровых каналов дисперсными примесями.

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

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

Научная новизна

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

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

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

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

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

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

Достоверность результатов

Достоверность результатов диссертационной работы обеспечивается:

1) математической моделью, основанной на общих законах и уравнениях механики сплошной среды;

2) тщательностью анализа физических процессов моделируемых явлений;

3) справедливостью используемых упрощений и приближений;

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

Практическая ценность

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

Апробация работы

Основные результаты диссертационной работы докладывались и обсуждались на V международной межвузовской научной конференции молодых ученых, аспирантов и студентов «Молодые - наукам о Земле», г.Москва 2010г.; конференции «Актуальные вопросы геолого-гидродинамического моделирования и переоценки нефтяных ресурсов Республики Татарстан», посвященной 15-летию ООО «ЦСМРнефть» при АН РТ и 75-летию академика Р.Х.Муслимова, г.Казань 2009г.; международной научно-практической конференции «Инновационные технологии в геологии и разработке углеводородов», г.Казань 2009г. .; международной научно-практической конференции «Инновации и технологии в разведке, добыче и переработке нефти и газа», г.Казань 2010г.; Всероссийской школе-конференции молодых исследователей и V Всероссийской конференции «Актуальные проблемы прикладной математики и механики», посвященной памяти академика А.Ф.Сидорова, Абрау-Дюрсо, 2010.

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

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

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

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

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

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

Первая глава посвящена литературному обзору материала по теме диссертации. Вопросам математического описания методов повышения нефтеотдачи пластов, в частности, моделированию потокоотклоняющих технологий посвящены работы JI.K. Алтуниной, П.Т. Зубкова, Ю.И. Капранова, Ф.И. Котяхова, Ф.Ф. Крейга, В.М. Конюхова, А.И. Никифорова, K.M. Федорова, P.C. Хисамова, А.Н. Чекалина и других исследователей.

Сказано о том, что идея использования полимеров в качестве потокоотклоняющей технологии возникла еще в конце 1970-х годов. Суть идеи состояла в том, чтобы полимер сам «нашел» место притока воды внутри скважины и заблокировал его, не мешая при этом течению нефти. Технология термогелевого заводнения основана на закачке гелеобразующих составов, полимеризация которых в пласте происходит за счет повышения температуры ГОС. На практике хорошо зарекомендовали себя составы, предложенные группой JI.K. Алтуниной. Найдены способы регулирования температуры и времени гелеобразования.

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

модели базируются на кинетических соотношениях вида —■ = f(m,...). К

примеру, таким образом моделируется процесс заводнения с применением термогеля в работах К.М. Федорова и соавторов. Вместе с тем отмечается, что предпринимались попытки использования теории коагуляции Смолуховского. В частности в работе S. Hubbard, L.J. Roberts, K.S. Sorbie данная теория коагуляции используется для моделирования кинетики ксантан/хромной гелевой системы. В отличие от настоящей работы ядро коагуляции в ней считается не зависящей от температуры и концентрации реагента.

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

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

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

at ôt

— (mRS„. + a) + div(RU,) = 0, et

где U0,UW - скорости фильтрации нефти и воды; So, Sw - соответственно нефтенасыщенность и водонасыщенность, 0 s S, s 1, ( i = о, w ), причем S0 + Sw = 1 ; R - концентрация примеси; m - пористость; a = mRSwr - изотерма сорбции Генри.

Пренебрегая влиянием гравитационных и капиллярных сил, уравнения движения фаз взяты в виде обобщенного закона Дарси:

и о =-K-grad(P), и, =-^grad(P), M. M.

где Mi - динамическая вязкость; Р - давление в фазах; К, = - фазовые проницаемость, к - абсолютная проницаемость пласта.

В работах C.B. Избаша, Ю.М. Шехтмана показано, что скорость движения частиц размером менее 10 мкм не отличается от скорости движения самой жидкости, на основании чего принято, что собственно частицы движутся со скоростью воды, не изменяя ее вязкости.

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

ОмД Т=Т5

Рисунок 1. Область расчета поля температур

Для вычисления температуры пластовой системы использована расширенная в направлениях от кровли и подошвы область (рис.1). Во всей области поле температур описывалась уравнением:

5СТ

д1

- йЫ

Т ^ с,р,и, = Я <И\ Т,

где С = т + (1 - 5„,)с0р0) + (1 - т)сррр - объемная теплоемкость пористой

среды; с, - удельная теплоемкость фаз; р, - плотность фазы; Л - коэффициент теплопроводности. На кровле АР и подошве ВЕ пласта выполняется условия

равенства потоков тепла и температур: Я

дТ дп;,

Х8Т дп

-Т |г.. Полагается,

что на входе в пласт (на границе АВ) температура равна температуре нагнетаемой воды Т=Т\ а на удаленных границах СБ и ОН температура неизменна и может отличаться по значению. На остальной части границы (ВС, Эв и НА) задано условие равенства нулю градиента температур по нормали к границе. Во внешних областях от кровли и подошвы пласта предполагается, что происходит только кондуктивный перенос тепла. На удаленных границах задается постоянная температура.

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

д

записывается так называемое "уравнение сплошности": ^ (",

)+и= 0, где

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

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

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

где Р (/? > 0) - коэффициент формы частицы; ч„ - средняя скорость движения жидкости в поровом канале, которая связана с суммарной скоростью фильтрации и соотношением:

где - коэффициент извилистости капилляра.

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

функцией распределения частиц по объемам у/ = , для которой

<> » дш

выполняется уравнение сплошности : —- = и., где и( - интенсивностью

изменения количества частиц одного типоразмера в потоке, которая

иц = талр

\и'~ г

и,

представлена в виде двух слагаемых: и( = и" +ис(, где первое слагаемое отвечает за конвективный перенос частиц: 1

и™ =--—и „Чу/,

т5.(\ + Г)

второе - за объединение агрегатов между собой, которое представлено в виде уравнения Смолуховского:

1 г

о

где 0 = С, 1п(1 + Г-Г„ж1)(/гх.,)£- - ядро коагуляции, коэффициенты С „С.

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

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

|ггфс1г <1г

о /о

г— /'„„

к = к° \г*(р<Ь \г\° ¿г

о

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

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

12

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

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

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

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

В пункте 3.3 представлена модель, основанная на применении кинетических соотношений для проницаемости и соотношения Козени-Кармана для пористости. На базе примера пункта 3.2 показано, что можно

13

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

Рисунок 3. Коэффициент извлечения нефти

В пункте 3.4 анализировалось влияние сроков закачки полимера и вязкости нефти на эффективность применения технологии. Расчеты проводились на фрагменте вертикального сечения трехслойного пласта. Проницаемости слоев полагались равными: в верхнем и нижнем пропластке -20 мД, в среднем - 200 мД. Рассматривалось три случая:

1) вязкость нефти 5 мПа-с, закачка полимера при обводненности 90%;

2) вязкость нефти равна 20 мПа с, закачка полимера при обводненности 90%;

3) вязкость нефти 20 мПа с, закачка полимера при 95%. Во всех случаях гель закачивался в течении суток при заданных значениях забойного давления. Результаты расчетов приведены на рисунках 4-6.

200 300 400

Жидкость, тыс.мЗ

Рисунок 4. Зависимость добытой нефти от отобранной жидкости.

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

-без закачки ГОС -с закачкой ГОС

- без закачки ГОС

- закачка ГОС при обводненности 90% ■ закачка ГОС при обводненности 95%

Рисунок 5. Обводненность продукции.

Время, годы

Рисунок 6. Коэффициент извлечения нефти.

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

Рисунок 7. Изменение проницаемости пласта после закачки термогеля для разных случаев неоднородности (мД).

В пункте 3.6 изучается влияние учета зависимости вязкости флюидов от температуры на результаты расчетов. Рассмотрено два варианта: в первом случае предполагается, что вязкости флюидов не зависят от температуры; во втором - используются следующие зависимости вязкости нефти и воды от температуры

0.001787 (1 +0.03377"+ 0.000221Г3) ' 17

для воды формула Пуазейля: //, =

для нефти зависимость из литературы (рис.8): I1,, —

Рпс. 99. Зависимость ннзкт.тп исфти ц от температуры Т

Рисунок 8. Зависимость вязкости нефти от температуры.

О "'1............1..... 1м*ми;|- '..........."............—О дпп

Рисунок 9. Изменение проницаемости: слева - для постоянных вязкостен, справа-для пересчитываемых.

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

Рисунок 10. Коэффициент извлечения нефти и обводненность продуктии

бремя, гады

-постоянная вязкость

— — переменная вязкость

Жидкость, тыс.мЗ

Рисунок 11. Кривые вытеснения.

В пункте 3.7 анализировалось влияние температуры реагирования термогеля на параметры разработки пласта. Рассматривались варианты закачки в пласт термогеля при прочих равных параметрах с разной температурой реагирования: 30 °С, 50 °С и 78 "С. Результат приведен на рисунке 12. Выявлено, что корректный подбор параметров термогеля является одним кз определяющих факторов эффективного применения технологии.

Рисунок 12. Проницаемость пласта после закачки термогеля с температурой реагирования, °С: а) 30, б) 50, в) 78.

Рисунок 13. Кривая вытеснения и коэффициент извлечения нефти.

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

Основные результаты и выводы, полученные в диссертации

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

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

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

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

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

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

Список опубликованных работ по теме диссертации

1. Нуруллин Р. Ф„ Никифоров А.И. Об эффективности применения термогелей при заводнении нефтяных пластов // НТЖ Нефтяное хозяйство. -2010. - №6. - С.65-67.

2. Нуруллин Р. Ф„ Никифоров A.M. Анализ эффективности применения тёрмогелей при заводнении нефтяных пластов // НТЖ Нефть. Газ. Новации. -2010. - № 1. - С.64-68.

3. Нуруллин Р.Ф., Никифоров А.И. Некоторые результаты моделирования заводнения нефтяных пластов с применением термогелей // Сборник трудов Седьмой международной научно-практической конференции "Исследование, разработка и применение высоких технологий в промышленности". - Санкт-Петербург, 2009. - С.380

4. Нуруллин Р.Ф., Никифоров А.И. Об эффективности применения термогелей при заводнении нефтяных пластов. // Сборник трудов Международной научно-практической конференции «Инновационные технологии в геологии и разработке углеводородов»,- Казань: Изд-во НПО «Репер»,2009. С.340-345.

5. Нуруллин Р.Ф., Никифоров А.И. Математическое моделирование процессов довытеснения запасов нефти с использованием термогелей II Сборник трудов 5 международной научной конференции студентов, аспирантов и молодых ученых «Молодые - наукам о Земле». - Москва, 2010. - С.240.

6. Никифоров А.И., Нуруллин Р. Ф. О моделировании заводнения нефтяных пластов с применением термогеля. II Тезисы докладов Всероссийской школы-конференции молодых исследователей и V Всероссийской конференции «Актуальные проблемы прикладной математики и механики», посвященной памяти академика А.Ф.Сидорова. Екатеринбург: УрО РАН, 2010. С.62-63.

Отпечатано в типографии «АВС-Компания» тел. 8-91-96-46-70-13

Заказ № 61 от 21.04.2010 г. Формат 60x84 1/16. Усл. печ. л. 1,5. Бумага 80 г. Печать ризографическая. Тираж 120 экз.

ВВЕДЕНИЕ.

ГЛАВА 1. СОСТОЯНИЕ ПРОБЛЕМЫ ПРОГНОЗИРОВАНИЯ РЕЗУЛЬТАТОВ ПРИМЕНЕНИЯ ПОТОКО ОТКЛОНЯЮЩИХ МУН.

ГЛАВА 2. ПОСТАНОВКА ЗАДАЧИ ДВУХФАЗНОЙ ФИЛЬТРАЦИИ С ВЗАИМОДЕЙСТВУЮЩИМИ ПРИМЕСЯМИ И ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ.

2.1 Уравнения движения и сохранения масс фаз и примеси.

2.2 Распределение температуры в пласте.

2.3 Функция распределения пор по размерам.

2.4 Функция распределения частиц по размерам.

2.5 Проницаемость.

2.6 Пористость.

2.6 Блокирование поровых каналов.

2.7 Состав частиц в потоке.

2.8 Численная реализация.

2.9 Анализ сходимости и сравнение с аналитическим решением

ГЛАВА 3. ЗАВОДНЕНИЕ С ПРИМЕНЕНИЕМ ТЕРМОГЕЛЕЙ.

3.1 Изменение структуры порового пространства.

3.2 Вытеснение нефти в послойно неоднородной пластовой системе.

3.3 Об одном подходе к моделированию изменения фильтрационных свойств пласта при закачке полимера.

3.4 Анализ влияния сроков закачки полимера и вязкости нефти на эффективность применения технологии.

3.5 Элемент пятиточечной системы заводнения.

3.6 Учет зависимости изменения вязкости нефти и воды от температуры.

3.7 Влияние температуры реагирования.

3.8 Анализ влияние параметров кровли и подошвы пласта на результат прогнозирования эффективности технологии.

Актуальность темы.

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

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

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

Цель работы

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

Задачи и методы исследования.

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

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

Моделирование процессов блокирования поровых каналов дисперсными примесями.

Моделирование изменения пористости и проницаемости нефтесодержащего пласта, вызванного изменением структуры порового пространства.

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

Научная новизна

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

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

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

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

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

Достоверность результатов

Достоверность результатов диссертационной работы обеспечивается: 1) математической моделью, основанной на общих законах и уравнениях механики сплошной среды;

2) тщательностью анализа физических процессов моделируемых явлений;

3) справедливостью используемых упрощений и приближений;

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

Практическая ценность

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

Апробация работы

Основные результаты диссертационной работы докладывались и обсуждались на V международной межвузовской научной конференции молодых ученых, аспирантов и студентов «Молодые - наукам о Земле», г.Москва 2010г.; конференции «Актуальные вопросы геолого-гидродинамического моделирования и переоценки нефтяных ресурсов Республики Татарстан», посвященной 15-летию ООО «ЦСМРнефть» при АН РТ и 75-летию академика Р.Х.Муслимова, г.Казань 2009г.; международной научно-практической конференции «Инновационные технологии в геологии и разработке углеводородов», г.Казань 2009г. .; международной научно-практической конференции «Инновации и технологии в разведке, добыче и переработке нефти и газа», г.Казань 2010г.; Всероссийской школе-конференции молодых исследователей и V Всероссийской конференции «Актуальные проблемы прикладной математики и механики», посвященной памяти академика А.Ф.Сидорова, Абрау-Дюрсо, 2010.

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

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

Диссертация состоит из введения, 3 глав, заключения, списка использованной литературы из 93 наименований, содержит 117 страниц машинописного текста и 60 рисунков.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

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

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

1. Avery, M. R. Field evaluation of production well treatments in Kansas using crosslinked cationic polymer gel Text. / M. R. Avery, J. A. Sutphen // Proc. 8th Tertiary Oil Recovery Conf., Wichita, Kansas, 1991, v. 11. P. 31-42.

2. Chauveteau, G. Controlling in-situ gelation of Polyacrylamide by zirconium for water shutoff Text. / G. Chauveteau, R. Tabary, M. Renard, A. Omari // Symp. on oilfield chemistry, Houston, TX, 16-19 Feb. 1999. SPE 50752. p. 1-10.

3. Coen, E. M. Modeling particle size distributions and secondary particle formation n emulsion polymerization Text. / E. M. Coen, R. G. Gilbert, B. R. Morrison, H. Leue, S. Reach // Polymer, 1998, 39, № 26. P. 70997112.

4. Ding, Y. A Control Volume for Flow Simulation in Well Vicinity for Arbitrary Well Configurations Text. / Y. Ding, P. Lemonnier, T. Estebenet, J-F. Magras // Society of Petroleum Engineers, 1998, № 48854, pp. 1-11.

5. Gruesbeck, C. Entrainment and Deposition of Fine Particles in Porous Media Text. / C. Gruesbeck, R. E. Collins // Society of Petroleum Engineers Journal, December 1982. P. 847-856.

6. Hubbard, S. Experimental and theoretical investigation of time-setting polymer gels in porous media Text. / Hubbard S., Roberts L.J., Sorbie K.S. // SPE Reservoir Engineering, April 1988. P.1257-1267.

7. Patil, D. P. An analytical solution to continuous population balance model describing floe coalescence and breakage A special case Text. / D. P. Patil, J. R. G. Andrews // Chemical Engineering Science, 1998, v. 53, № 3. -P. 599-601.

8. Sahimi, M. Hydrodynamics of Particulate Motion in Porous Media Text. / M. Sahimi, A. O. Imdakm // Phys. Rev. Let., March 1991, 66, № 9. P. 1169-1172.

9. Saiers, J. E. First- and second-order kinetics approaches for modeling the transport of colloidal particles in porous media Text. / J. E. Saiers // Water Resources Research, 1994, 30, № 9. P. 2499-2506.

10. Sharma, M. M. Transport of Particulate Suspensions in Porous Media: Model Formulation Text. / M. M. Sharma, Y. C. Yortsos // AIChE Journal, October 1987, 33, №10. P. 1636-1643.

11. Todd, B. J. A Mathematical Model of In-Situ Gelation of Polyacrylamide by a Redox Process Text. / B.J. Todd, G. P. Willhite, D. W. Green // SPE Reservoir Engineering, February 1993. P. 51-58.

12. Азиз, X. Математическое моделирование пластовых систем Текст. : [пер. с англ.] / X. Азиз, Э. Сеттари. М.: Недра, 1982. - 407 с.

13. Алтунина, Л. К. Увеличение нефтеотдачи пластов композициями ПАВ Текст. / Л. К. Алтунина, В. А. Кувшинов. .Новосибирск: Наука, 1995 .-198с.

14. Алтунина, Л.К. Регулирование кинетических и реологических характеристик гелеобразующих систем для увеличения нефтеотдачи

15. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидромеханика: Учебник для вузов.-М.: Недра, 1993, 416 с.

16. Булыгин, В. Я. Гидромеханика нефтяного пласта Текст. / В. Я. Булыгин М.:Недра, 1974. - 232 с.

17. Булыгин, В. Я. Одновременная фильтрация двух жидкостей в нефтяном пласте Текст. / В. Я. Булыгин // Уч. зап. Казанского ун-та. Вопросы подземной гидромеханики / Казань: изд-во КГУ. 1966. -Т.125, KH.8.-C. 3-36.

18. Волощук, В. А. Кинетическая теория коагуляции Текст. / В. А. Волощук. — Л.: Гидрометеоиздат, 1984. 284 с.

19. Газизов, А. Ш. Воздействие полимер дисперсных систем на «языки» прорыва воды Текст. / А. Ш. Газизов [и др.] // Численные методы решения задач фильтрации и оптимизации нефтедобычи. Казань: Физ.-техн. ин-т Казан, фил. АН СССР, 1990. - С. 25-30.

20. Газизов, А. Ш. О математической модели фильтрации нефти в неоднородном по толщине пласте при воздействии полимердисперсными средами Текст. / А. Ш. Газизов [и др.] //

21. Вопросы математического моделирования процессов фильтрации и рациональной разработки нефтяных месторождений. Казань: Физ.-техн. ин-т Казан, фил. АН СССР, 1989. - С. 5-9.

22. Газизов, А. Ш. Оценка эффективности технологии применения полимеро-дисперсной системы по результатам промысловых исследований Текст. / А. Ш. Газизов, Р. X. Низамов // Нефтяное хозяйство, 1990, № 7. С. 49-52.

23. Галкин В.А. Уравнения Смолуховского. М.: Физматлит, 2001. - 336 с.

24. Де Грот С., Мазур П. Неравновесная термодинамика. М.:Мир, 1964. - 456с.

25. Ентов, В. М. Гидродинамика процессов повышения нефтеотдачи Текст. / В. М. Ентов, А. Ф. Зазовский. М.: Недра, 1989. - 232 с.

26. Желтов Ю.П. Механика нефтегазоносного пласта. М., «Недра», 1975, 216 с.

27. Зак, С. А. Оценка снижения проницаемости в процессах кольматации пористой среды Текст. / С. А. Зак, Е. Д. Миневич, Э. Чен-Син // Интенсификация добычи нефти. М.: 1983, вып.85. - С. 42-50.

28. Закиров, Э. С. Трехмерные многофазные задачи прогнозирования, анализа и регулирования разработки месторождений нефти и газа Текст. / Э. С. Закиров. М.: Изд-во «Грааль», 2001. - 303 с.

29. Закиров, Р. X. Экспериментальные исследования движения дисперсных частиц в пористой среде Текст. / Р. X. Закиров, Т. А. Захарченко // Фильтрация многофазных систем (материалы Всесоюзного семинара). Новосибирск: 1991, с.24-26.

30. Золотухин, А. Б. Моделирование процессов извлечения нефти из пластов с использованием методов нефтеотдачи Текст. / А. Б. Золотухин. М.: МИНГ, 1990. - 267 с.

31. Зубарев, А. Ю. Модель фильтрации суспензии с учетом кольматации и суффозии Текст. / А. Ю. Зубарев, Б. Хужаеров // Инженерно-физический журнал, 1988, т. 55, № 3. С. 442-447.

32. Зубков, П. Т. Влияние гелевых барьеров на течение воды и нефти в неоднородном пористом пласте Текст. / П. Т. Зубков, К. М. Федоров // Изв. РАН. МЖГ, 1995, № 2. С. 99-107.

33. Зубков, П. Т. Механизм формирования высоковязких барьеров в неоднородных нефтяных пластах Текст. / П. Т. Зубков, К. М. Федоров // Изв. РАН. МЖГ, 1994, № 1. С. 98-103.

34. Избаш, С. В. Фильтрационные деформации грунта Текст. / С. В. Избаш // Изв. ВНИИГ, 1933, № 10. С. 5-21.

35. Капранов, Ю. И. Изменения поровой структуры в потоке монодисперсной взвеси Текст. / Ю. И. Капранов // ПМТФ, т.41, №2, 2000.-С. 113-121.

36. Капранов, Ю. И. О фильтрации взвеси твердых частиц Текст. / Ю. И. Капранов // Прикладная математика и механика, т. 63, вып. 4, 1999.-С. 620-628.

37. Капранов, Ю. И. Структурная модель процесса механической кольматации пористой среды Текст. / Ю. И. Капранов // Динамика сплошной среды. АН СССР. Сиб. Отд-ние. Ин-т гидродинамики, 1989. Вып. 90.-С. 27-39.

38. Киркпатрик, С. Перколяция и проводимость Текст. / С. Киркпатрик // Теория и свойства неупорядоченных материалов. Сер. Новости физики твердого тела. М.: Мир, 1977. - С. 94-115.

39. Кондратьев, В. Н. Фильтрация и механическая суффозия в несвязанных грунтах Текст. / В. Н. Кондратьев. Симферополь: Крымиздат, 1958. - 76 с.

40. Коннор, Дж. Метод конечных элементов в механике жидкости. Пер. с англ. Текст. / Дж. Коннор, К. Бреббиа. JL: Судостроение, 1979. -264 с.

41. Коновалов, А. Н. Задачи фильтрации многофазной несжимаемой жидкости Текст. / А. Н. Коновалов. Новосибирск: Наука. Сиб. Отд-ние, 1988.-166 с.

42. Конюхов, В. М. Образование и фильтрация оторочки загустителя в слоистых неоднородных пластах Текст. / В. М. Конюхов,

43. A. В. Костерин, А. Н. Чекалин // Изв. РАЕН, серия МММИУ, 1997, т. 1,№ 1.-С. 84-109.

44. Конюхов, В. М. Образование, динамика и эффективность оторочек загустителя при заводнении пластов сложного строения Текст. /

45. Коробов, К. Я. Изучение фильтрации жидкостей различной вязкости при малых градиентах Текст. / К. Я. Коробов, В. А. Трефилов, Ю. В. Антипин // Тр. УНИ. Уфа: Башкнигоиздат, 1975, вып. 8. С. 84-90.

46. Котяхов, Ф. И. Основы физики нефтяного пласта Текст. / Ф. И. Котяхов. М.:Гостоптехиздат, 1956. — 362с.

47. Крейг, Ф. Ф. Разработка нефтяных месторождений при заводнении Текст./ Ф. Ф. Крейг. М.:Недра, 1974. - 192 с.

48. Левин, В. Г. Физико-химическая гидродинамика Текст. / В. Г. Левич.- М.:Физматгиз, 1959. 699 с.

49. Лозин, Е. В. Результаты разработки и принципы доразработки объектов с заводнением Текст. / Е. В. Лозин // Нефтяное хозяйство, 1987, №12. С.36-39.

50. Лойцанский, Л. Г. Механика жидкости и газа • Текст., / Л. Г. Лойцанский. М.: Наука, 1973. - 848 с.

51. Минц, Д. М. Кинетика фильтрации малоконцентрированных суспензий на водоочистных фильтрах Текст. / Д. М. Минц // Докл. АН СССР, 1951, т. 72.

52. Михайлов, Н. Н. Изменение физических свойств горных пород в околоскважинных зонах Текст. / Н. Н. Михайлов. М.: Недра, 1987.- 152 с.

53. Муслимов, Р. X. Геология, разрабока и эксплуатация Ромашкинского нефтяного месторождения Текст. В 2-х томах. Т. 2. / Р. X. Муслимов, А. М. Шавалиев, Р. Б. Хисамов, И. Г. Юсупов,- М.: ВНИИОЭНГ, 1995.-286 с.

54. Нигматулин, Р. И. Динамика многофазных сред Текст. В 2 ч. 4.2. / Р. И. Нигматулин. -М.: Наука. Гл. Ред. Физ.-мат. Лит. 1987, 360 с.

55. Никифоров, А. И. О моделировании вытеснения нефти водой с гелеобразующими добавками Текст. / А. И. Никифоров, С. В. Анохин // Математическое моделирование. 2002. - Т. 14.-12. - С. 117-127

56. Никифоров, А. И. Перенос частиц двухфазным фильтрационным потоком Текст. / А. И. Никифоров, Д. П. Никаныиин // Математическое моделирование. 1998. — Т. 10. - 6. - С. 42-52.

57. Нуруллин Р.Ф., Никифоров А.И. Об эффективности применения термогелей при заводнении нефтяных пластов // НТЖ Нефть. Газ. Новации. -2010. №2.

58. Нуруллин Р.Ф., Никифоров А.И. Об эффективности применения термогелей при заводнении нефтяных пластов // НТЖ Нефтяное хозяйство. 2010.

59. Островский, Г.М. Прикладная механика неоднородных сред Текст. / Г. М. Островский. СПб.: Наука, 2000. - 359 с.

60. Резник, Г. Д. О поглощении частиц в фильтрах водопроводных очистных сооружений Текст. / Г. Д. Резник. Инженерно-физический журнал, 1999, т.72, № 1. - С. 38-42.

61. Саламатин, А. Н. Математические модели дисперсных . потоков Текст. / А. Н. Саламатин. Казань: Изд-во Казанского ун-та, 1987. -172 с.

62. Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Проектирование разработки Текст. / Ш. К. Гиматудинов [и др.] М.:Недра, 1983. - 463 с.

63. Сургучев, М. JI. Выступления. М. JI. Сургучев Текст. / Регулирование процессов эксплуатации нефтяных залежей. М.: Наука, 1977.-С. 215-220.

64. Технология ограничения водопритоков на основе алюмосиликата и математическое моделирование ее применения в продуктивных пластах Текст. / Ю. А. Котенев [и др.] // Разработка и эксплуатация нефтяных месторождений, 2004, № 4. — С. 60-63

65. Тодес, О. М. К теории коагуляции и укрупнения частиц в золях. Кинетика укрупнения частиц при "перегонке" вещества через гомогенную фазу Текст. / О. М. Тодес // Журн. Физ. Химия, 1946, Т. 20, Вып. 7, с. 629-644.

66. Тронов, В. П. Фильтрационные процессы и разработка нефтяных месторождений Текст. / В. П. Тронов. — Казань: Изд-во «Фэн» Академии наук РТ, 2004. 584 с.

67. Федоров, К. М. Аналитические исследования процесса гелеобразования в призабойной зоне скважин Текст. / К. М. Федоров // Изв. РАН. МЖГ, 1997, №4. С. 80-87.

68. Федоров, К. М. Геолого-математическое моделирование термогелеобразующего воздействия на пласт Текст. / Химия нефти и газа: Материалы IV международной конференции "Химия нефти и газа". В 2-х т. Томск: "8ТТ", 2000. Т. 2. - С.134-139.

69. Флетчер, Р. Численные методы на основе метода Галеркина: Пер. с англ. Текст. / Р. Флетчер. М.: Мир, 1988. - 352 с.

70. Хазе Р. Термодинамика необратимых процессов. М.:Мир, 1967. -544с.

71. Хейфец, Л. И. Многофазные процессы в пористых средах Текст. / Л. И. Хейфец, А. В. Неймарк. М.: Химия, 1982. - 320 с.

72. Хужаеров, Б. Модель многокомпонентной кольматационносуффози-онной фильтрации Текст. / Б. Хужаеров // Инженерно-физический журнал, 1994, т. 66, №4. С. 427-433.

73. Хужаеров, Б. Модель фильтрации суспензии с учетом кольматации и суффозии Текст. / Б. Хужаеров // Инженерно-физический журнал, т. 63, № 1, 1992.-С. 72-79.

74. Чекалин, А. Н. Исследование двух- и трехкомпонентной фильтрации в нефтяных пластах Тест. / А. Н. Чекалин, Г. В. Кудрявцев, В. В. Михайлов. Казань: Из-во Казанского университета, 1990. - 147 с.

75. Численные методы в динамике жидкости. Пер. с англ. Текст. / Э. Джеймсон [и др.]; под ред. Г. Вирца, Ж. Смолдерна. М.: Мир, 1981. - 407 с.

76. Шехтман, Ю. М. Фильтрация малоконцентрированных суспензий Текст. / Ю. М. Шехтман. М.:Недра, 1961. - 161 с.

77. Щелкачев В.Н., Лапук Б.Б. Подземная гидравлика Текст. / В.Н. Щелкачев, Б.Б. Лапук Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001, - 736 с.

78. Энгельгардт, В. Поровое пространство осадочных пород. Перевод с немецкого Текст. / В. Энгельгардт. М.: Недра, 1964. - 232 с.