Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва




Скачать 97.79 Kb.
НазваниеЮлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва
Дата публикации26.03.2013
Размер97.79 Kb.
ТипДокументы
shkolnie.ru > Математика > Документы
Михаил Владимирович Дмитриевский -

кандидат физико-математических наук,

ЗАО «ИНТЕРА»
Юлия Александровна Плиткина -

аспирант кафедры математического моделирования ТюмГУ
Математическое моделирование полей давлений с целью определения направления трещины гидроразрыва

УДК 519.63+519.65




АННОТАЦИЯ. В данной работе предлагается комбинированный метод восстановления полей давлений, основанный на двух классических подходах: интерполяции «пластовых» давлений и математическом моделировании процессов фильтрации. Предлагаемый метод имеет вариационную постановку (минимизация последовательности функционалов), и позволяет учитывать зачастую противоречивые значения дебитов и забойных давлений, заданных в скважинах. Такая особенность метода дает возможность по выбранному критерию определить примерное направление трещины гидроразрыва.
ANNOTATION. This paper presents the combined method of pressure fields reconstruction, based on two classical approaches: interpolation of “formation” pressures and mathematical modeling of filtration processes. This method has variation statement (consecutive minimizing of functionals), and allows us to take into account values of rates and bottom-hole pressures given in the wells, which are often inconsistent. This feature of the method enables to define an approximate direction of a fracture by the chosen criterion.



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

  1. снижение или увеличение пластовых давлений в отдельных зонах пласта;

  2. зоны отбора нефти, зоны влияния нагнетательных скважин;

  3. средневзвешенные пластовые давления, рассчитываемые по зонам отбора;

  4. направление движения контура нефтеносности.


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

  1. замеры забойных давлений - замеры, произведенные на забоях скважин во время их работы;

  2. замеры “пластовых” давлений - замеры, произведенные на забоях скважин через некоторое время после их остановки, необходимое для того, чтобы давление в ''достаточной'' степени стабилизировалось;

  3. дебиты / приемистости вертикальных, наклонных и горизонтальных скважин;

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

  5. геометрические и фильтрационные параметры трещин гидроразрыва пласта (ГРП);

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

  7. фильтрационные параметры пласта – коэффициенты проницаемости, вязкости, пьезопроводности и т.д.


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

  1. ^ Интерполяция “пластовых” давлений. В качестве исходной информации для подобных методов служат замеры “пластовых” давлений в скважинах, а также значения давлений, заданные на контуре питания. Заметим, что подобная информация не может быть получена на определенную дату одновременно по всем скважинам - для этого пришлось бы временно остановить разработку всего месторождения, поэтому замеры “пластовых” давлений, полученные на предыдущую дату, пересчитываются на дату построения карты. В результате этого пересчета в данные вносится дополнительная погрешность. Для получения карты полученные исходные данные пересчитываются в узлы равномерной сетки одним из интерполяционных методов (крайгинг, сплайн-интерполяция, метод Шепарда и др.). Достоинствами методики являются простота программирования алгоритмов и использования полученных программ пользователями. К недостаткам можно отнести тот факт, что получаемая картина часто физически противоречива. Для достижения правдоподобно выглядящего результата обычно приходится прибегать к добавлению в расчет новых, искусственно создаваемых исследователем данных - как правило это контура, с заданными на них значениями пластовых давлений.

  2. ^ Математическое моделирование процессов фильтрации. Основу данного подхода составляет создание математической модели, описывающей процесс фильтрации жидкости в пласте [3,4]. Обычно такая модель описывается системой дифференциальных уравнений с заданными начальными и граничными условиями. В зависимости от типа выбранной модели варьируется вид и количество уравнений. Существующие математические модели позволяют учитывать массу разнородной информации, что можно рассматривать и как преимущество и как недостаток: часто приходится адаптировать модель к исходным данным, варьируя ее параметрами. Этот процесс очень трудоемкий и занимает основное время построения фильтрационной модели.


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

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

(1)

где ,

,

,

- координаты центра первой ячейки. Задача построения поля давлений сводится к определению неизвестных коэффициентов , которые, в виду выбора “ступенчатой” базисной функции, определяют значение давления в ячейке.

Относительно ячеек будем считать истинными следующие гипотезы:

  1. давление в ячейке постоянно;

  2. вязкость и плотность жидкости в ячейке постоянна;

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

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

, ,

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

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


  • ^ Учет информации о давлениях.

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

, (2)

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

. (3)

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

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

, (4)

где - элемент контура, проходящий через ячейку .

  • Учет информации о дебитах вертикальных, наклонных и горизонтальных скважин, зонах отсутствия фильтрации, линиях ГРП.

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

,



(5)

, (6)

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

. (7)

Дифференцирование выражения (7) по переменной приводит к закону сохранения массы, записанного для ячейки :

. (8)
Функционал (7) позволяет учитывать дебиты/приемистости вертикальных и горизонтальных скважин, а также влияние трещин ГРП на поле пластовых давлений. Массовая интенсивность внешних источников и стоков для ячейки складывается из дебита и приемистости скважин, попадающих в эту ячейку. Для моделирования наклонных и горизонтальных скважин, а также трещин ГРП, их массовая интенсивность распределяется по какому-либо закону (например, равномерно) по ячейкам, в которых проходит ствол скважины или трещина.
Границы зон отсутствия фильтрации, описанные при помощи кривых, пересчитываются в зоны, заданы с помощью граней ячеек (для этого очень хорошо подходит алгоритм Брезенхема [5]). Тем самым у некоторых ячеек появляются грани, через которые не может происходить фильтрация жидкости. В этом случае функционал должен быть модифицирован так, чтобы соответствующая скорость фильтрации в формуле (5) через непроницаемую грань равнялась нулю.
^ Ранжирование данных по степени их достоверности осуществляется за счет комбинирования функционалов (2) и (7) для каждой из ячеек сетки. Рассмотрим возможные ситуации:

  1. ^ В ячейке абсолютно точно известно давление. В этом случае для ячейки минимизируется функционал (2).

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

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

, (9)

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


Рис. 2. Учет противоречивой информации
Всего было построено три поля давлений:

  1. Первая картинка - в центре круга была задана скважина с забойным давлением 200 атм, при этом записывался функционал вида (2).

  2. Вторая картинка - в центре круга была задана скважина с дебитом жидкости, равным , при этом записывался функционал вида (7).

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


Как видно из рис. 2, поле давлений, изображенное на третьей картинке, является "промежуточным" между первыми двумя.

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


ЛИТЕРАТУРА





  1. Данилов В.Л. Вариационный принцип наименьшей скорости рассеяния энергии при фильтрации жидкостей в пористой среде и его приложения. Москва-Ижевск: Институт компьютерных исследований, 2003. 108 с.

  2. Михлин С.Г. Вариационные методы в математической физике. Издательство ''Наука'', Главная редакция физико-математической литературы. Москва. 1970.

  3. Каневская Р.Д. Математическое моделирование гидродинамических процессов разработки месторождений углеводородов. Москва-Ижевск: Институт компьютерных исследований, 2003. 128 с.

  4. Костюченко С.В. Математическое моделирование полей давлений в нефтяных резервуарах с произвольными системами скважин различных профилей // Нефтяное хозяйство. 2000. №10. С.70-77.

  5. Порев В.Н. Компьютерная графика. Спб.: БХВ-Петербург, 2002. 432 с., илл.

Похожие:

Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconВопросы к экзамену по дисциплине «Математическое моделирование в...
Понятие математического моделирования. Требования, предъявляемые к математическим моделям
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconМатематическое моделирование электорального поведениЯ населения в...
Целью настоящей работы является попытка применения математического аппарата динамических систем на сетях для анализа некоторых вопросов...
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconМетод ускорения расчета самосогласованной задачи моделирования инжекционного лазера
Ключевые слова: инжекционный лазер, многослойные волноводы, математическое моделирование
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва icon"компьютерное моделирование электронных приборов" Методические указания...
Изучение принципов построения математических моделей физических процессов, лежащих в основе работы электронных приборов, принципов...
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconОтчет по курсу: «Математическое моделирование инженерно геодезических...
«Построение модели ситуации и рельефа по результатам имитационного моделирования на ЭВМ карты масштаба 1: 500»
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconМоделирование переноса взвешенных веществ на океаническом шельфе
Овос может быть решена только с использованием современных методов математического моделирования, поскольку в каждом конкретном случае...
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconМатематическое моделирование взаимодействия мощного лазерного фемтосекундного...
Математическое моделирование взаимодействия мощного лазерного фемтосекундного импульса с диэлектрическими материалами
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва icon"Моделирование колебаний математического маятника"
Смоделировать колебания математического маятника, изображение нити которого на экране менялась бы в размерах
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconТехнический отчет математическое моделирование инженерно-геодезических...
Построение модели ситуации и рельефа по результатам имитационного моделирования на ЭВМ топографической карты масштаба 1: 500 с сечением...
Юлия Александровна Плиткина аспирант кафедры математического моделирования Тюмгу математическое моделирование полей давлений с целью определения направления трещины гидроразрыва iconОмер Юсиф Мохамед Математическое моделирование и алгоритмы оценки...
Математическое моделирование и алгоритмы оценки эффективности системы администрирования локальной компьютерной сети
Вы можете разместить ссылку на наш сайт:
Школьные материалы


При копировании материала укажите ссылку © 2014
shkolnie.ru
Главная страница