- •1.Представление дифференциальных уравнений параболического типа в разностной форме
- •Рассмотрим представление в конечно – разностной форме одномерного дифференциального уравнения параболического типа
- •2.Понятие устойчивости и сходимости разностных уравнений. Условие устойчивости явной и неявной разностной схемы.
- •4.Постановка задачи о притоке реального газа в круговом пласте к скважине. Уравнение, граничные и начальные условия.
- •7. Особенности постановки задач двумерной однофазной фильтрации.
- •9.Метод переменных направлений.
- •10. Локально-одномерная схема Самарского а.А.
- •11.Схематизация залежи в случае однофазной фильтрации. Способ задания реальных скважин на модели.
- •13.Правило фаз Гиббса. Сведение многокомпонентных смесей к системе бинарных и тройных смесей.
- •16.Примеры решения площадных и профильных задач двухфазной фильтрации.
- •17.Об активном воздействии на водонапорный режим при разработке газовых месторождений
- •18. Методика решения задачи в случае трехмерной, трехфазной фильтрации (sip-метод)
- •19. Особенности задания начальных и граничных условий в случае пространственных задач фильтрации
- •20.Задание начальных и граничных условий в двумерном случае.
- •21. Результаты решения трехмерных трехфазных задач фильтрации применительно к разным технологиям разработки.
- •22. Влияние отдельных факторов на характер поведения кривых давления и насыщенности при решении задач фильтрации трёхфазных смесей.
- •23. Результаты прогнозирования показателей разработки реальной залежи.
- •24. Моделирование залежи. Получение замыкающих соотношений для решения трехфазных задач. Моделирование залежи: схематизация, размещение скважин по площади.
- •Получение замыкающих соотношений для решения трехфазных задач. Изменение плотности газа рассчитывалось по уравнению состояния для реального газа
- •25. Выбор исходных геолого-промысловых параметров при моделировании залежи.
- •Оценка запасов газа, конденсата и нефти по третьему объекту
- •26. Оценка снижения коэффициентов нефтеотдачи оторочек подстилающего типа за счет опережающей разработки газовой части пласта.
- •27. Расчет начальных и граничных условий при моделировании реальной залежи.
- •28. Разработка математической модели.
- •29. Модельные исследования.
- •Часть 2!!!!!!!!!!!!!!
- •Понятие о моделировании.
- •2. Моделирование фильтрации на макро и микро уровне
- •3. Геологические и гидродинамические модели
- •4. Разномасштабные модели фильтрации.
- •5. Существующие лицензированные программные пакеты. Понятие программного продукта.
- •6. Структура типового программного пакета для реализации задач моделирования нефтегазовых месторождений.
- •7. Возможности пре- и пост - процессора лицензированных программ.
- •8. Программное обеспечение для построения геологической модели (Petrel)
- •9. Исходные данные для гидродинамического моделирования процессов разработки нефтегазовых месторождений
- •10. Схематизация пласта и выбор расчетной модели.
- •11. Выбор модели фильтрации.
- •12. Постоянно- действующие модели. Методы и цель создания.
- •13. Основные этапы создания гидродинамической модели месторождения.
- •Воспроизведение истории разработки. Основные методы.
- •Уравнение сохранения массы при многофазной многокомпонентной фильтрации.
- •Абсолютная проницаемость. Методы получения. Способ задания.
- •Относительные фазовые проницаемости при двухфазной фильтрации. Методы получения. Способ задания.
- •Относительные фазовые проницаемости при трехфазной фильтрации. Методы получения. Способ задания.
- •Модель трехфазной фильтрации (Black oil).
- •Свойства флюидов и породы, учитываемые в модели трехфазной фильтрации (Black oil). Характерный вид зависимости.
- •Пористость и емкостные свойства пласта. Источники информации. Способ задания.
- •22. Моделирование трещиновато-порового пласта.
- •23. Данные о насыщенности и капиллярном давлении. Источники информации и способ задания.
- •24. Факторы, влияющие на определение размеров расчетных блоков.
- •26. Переход от геологической модели к гидродинамической. Процедура масштабирования (upscaling).
26. Оценка снижения коэффициентов нефтеотдачи оторочек подстилающего типа за счет опережающей разработки газовой части пласта.
Рассмотрим количественные оценки потерь нефти, которая уходит в газовую зону из-за снижения в ней давления на примере БУ10,11
Запасы газа, конденсата и нефти
№ УКПГ |
Газ, млрд. куб. м |
Конденсат, млн. т. |
Нефть, млн. т. |
1,2,5,8 |
448,5 |
74,0 |
481,7 |
Всего было просчитано 6 вариантов разработки залежи до 1995 г. /11 лет разработки/.
Первый вариант предполагает отбор в течении 8 лет по 2% от общего объема газа с каждого УКПГ. Темп отбора нефти:
- с 3(28скв) – 11(146скв)годы- 8 УКПГ
- с 5(28скв) – 11(200скв)годы- 2 УКПГ
- с 7(28скв) – 11(172скв)годы- 1 УКПГ
от 20т/сут – до12 т/сут
общее число СКВ.530 шт
Второй вариант - отбор газа в течение первых двух лет, а затем прекращение его отбора или проведение полного сайклинга, когда количество добытого газа и нагнетаемого равны.
Третий вариант - отбор только нефти
Четвертый вариант - характер отбора газа и нефти как в первом варианте, но отношение горизонтальной проницаемости k гор к вертикальной k вер на порядок выше, т.е. k гор / k вер = 1000.
/ Ухудшение вертикальной проницаемости имитирует увеличение глинистых пропластков в разрезе/.
Пятый вариант - все как в четвертом варианте, но отбор газа в первые три года по фактическим данным.
Шестой вариант - пятый вариант до 1990 г., затем прекращение отбора /имитация сайклннга/ на УКПГ- 1В, УКПГ- 2В
Потери нефти /млн. т/ во времени для различных вариантов
Варианты/ Годы |
1 |
2 |
3 |
4 |
5 |
6 |
1 |
4,3 |
4,3 |
- |
3,3 |
1,5 |
1,5 |
11 |
55,26 |
13,63 |
- |
48,483 |
46,86 |
38,00 |
Таким образом, рассмотрение рассчитанных вариантов показывает, что:
- разработка месторождений подобного типа с опережающим отбором газа приводит к расформированию нефтяной оторочки;
- консервация же запасов газовой шапки /проведение сайклинга в ней/ позволяет снизить или даже предотвратить потери нефти в нефтяной оторочке.
27. Расчет начальных и граничных условий при моделировании реальной залежи.
Если границы пласта непроницаемы, что имеет место на кровле и подошве пласта, то поток через такую границу равен нулю. Через боковые границы пласта может поступать вода. При этом величина притока флюида через внешнюю границу моделируется с помощью источников, расположенных в граничных узлах пласта. Плотность такого источника определяется по формуле
где
-
плотность источника в узле (i,j,k)
; к'-
коэффициент проницаемости в данном
узле; k3
,
- коэффициент фазовой проницаемости,
плотность, доля воды, коэффициент
динамической вязкости, давление в
водяной фазе, соответственно;
-
начальное пластовое, давление;
-
размеры ячейки пласта;
1
- коэффициент, учитывающий удаленность
контура питания пласта, от залежи
/например, при
= 10 в нашем случае контур питания находится
на расстоянии в 35 км/.
При задании начальных условий необходимо учитывать, что первоначально флюиды в пласте находятся в условиях капиллярно-гравитационного равновесия.
На рис.1 приведено примерное распределение по толщине пласта флюидов, находящихся в условиях капиллярно-гравитационного равновесия. Из рассмотрения рис.1 следует, что весь пласт может быть разделен на пять зон;
- Зону газовую, где нефть и вода содержатся как остаточные.
- Переходную зону между газом и нефтью, где насыщенности газа и нефти подчиняются капиллярно-гравитационному равновесию, а водонасыщенность равна остаточной.
- Нефтяную зону, где газ и вода присутствуют как остаточные, а нефтенасыщенностъ максимальна.
Рис 1.
Начальное распределение флюидов по
толщине пласта.
И, наконец, пятую зону - водяную, где насыщенность по воде максимальна, а газ и нефть присутствуют как остаточные.
С помощью капиллярных кривых определяются величины капиллярного давления, соответствующие величинам остаточной водо- и газонасыщенности, а также величинам максимальной насыщенности по воде.
Исходя из принятых для модели данных, величина капиллярного давления при остаточной газонасыщенности S1 = 0,060 равняется Рк1 = 0,0001 /все величины безразмерные/; величина капиллярного давления Рк2 = 0,001466 при максимальной насыщенности по газу S1 = 0,88; величина капиллярного давления Рк3 = 0,00008 при максимальной водонасыщенности S3 = 0,90; величина капиллярного давления Рк4 = 0,0007198 при остаточной водонасыщенности S3 = 0, 10.
Для определения начальных значений давлений и насыщенности флюидов в пласте сначала рассчитаем границы зон, которые выделяются на рис1.
Будем рассматривать способ задания начальных условий сверху вниз.
Учитывая, что толщина переходной зоны равняется примерно 7,5 м, а газонефтяной контакт /ГНК/ находится на отметке 2780м, определяем значение = 1,0306 (2772,5/2690) .
Величина определяется из условий капиллярно-гравитационного равновесия по формуле:
,
=1.0336
/здесь
и далее все величины безразмерные/,
-
плотность газовой фазы;
-
плотность нефтяной фазы;
G
=
; g
– ускорение силы тяжести;
=ZD0=
- 2690;
,
Р
0
- характерные значения плотности и
давления.
Давление в газовой фазе на границе переходной зоны определяется по формуле:
,
где - давление в газовой фазе в точке , - давление в газовой фазе на кровле пласта, т.е. при .
Давление в нефтяной фазе на границе переходной зоны, т.е. при определяется так:
Толщина нефтяной оторочки по промысловым данным равна R1= 10/2690 и поэтому величина
= + R1, =1.0373
Значение , а, следовательно, и толщина переходной зоны между нефтью и водой при условии капиллярно-гравитационного равновесия определяется по формуле:
,
=1.0393
Значение давления в нефтяной фазе на границе с нефтяной зоной определяется так:
Значение давления в нефтяной фазе на границе с переходной зоной между нефтью и водой, т.е. при , определяется по формуле:
Давление в этой же точке в водяной фазе определяется следующим образом:
Величина давления в водяной фазе между водяной зоной и переходной определяется по формуле:
После определения границ зон и давления на них, глубины узлов разностной сетки, моделирующей пласт, сравниваются с этими границами и рассчитываются фазовые давления в них.
При , т.е. в газовой зоне фазовые давления определяются по формуле:
Если , то узел находится в переходной зоне между газом и нефтью и фазовые давления определяются следующим образом:
В нефтяной зоне при
При узел разностной сетки расположен в переходной зоне между нефтью и водой и фазовые давления, определяются так:
В водяной зоне и фазовые давления определяются следующим образом:
