Гидродинамическое моделирование

Содержание

Слайд 2

Цель Grid секции • Геометрия пласта и основные свойства породы •

Цель Grid секции

• Геометрия пласта и основные свойства породы
• Размеры и

глубины ячеек сетки, DX, DY, DZ, TOPS, или COORD и ZCORN
• Эти данные использованы для введения ячеек связанных с течением флюида.
• Пористость и песчанистость для каждой ячейки, PORO, NTG
• Проницаемости для каждой ячейки PERMX, PERMY, PERMZ, PERMR, PERMTHT, PERMZ
• Модификаторы проводимости – ключевые слова “MULT”
• Они использованы для создания массивов поровых объемов, глубин и проводимостей для каждой ячейки
• Выше, на рисунке минимальные требований для этой возможности.
• Численные и сеточные водоносные пласты должны быть определены в секции GRID
• Опция активации ячеек (ключевое слово ACTNUM)
• Опция контроля зон выклинивания и минимального порового объема
достаточного для активности ячейки.
• Ссылка на список ключевых слов доступных в GRID секции
Слайд 3

Для каждой ячейки сетки Цель Grid секции

Для каждой ячейки сетки

Цель Grid секции

Слайд 4

Цель Grid секции ECLIPSE не использует размеры ячеек и их свойства

Цель Grid секции

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

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

Цель Grid секции Поровый объем рассчитывается по формуле PV = DX

Цель Grid секции

Поровый объем рассчитывается по формуле
PV = DX × DY

× DZ × NTG × φ
Где NTG – песчанистось, φ - пористость. Размеры ячейки для блочно-центрированной геометрии по X, Y и Z передаются с помощью ключевых слов DX или DXV, DY или DYV, DZ, или COORD и ZCORN для геометрии угловой точки.
Величина пористости также должна быть задана в каждой ячейке при помощи ключевого слова PORO, песчанистость (ключевое слово NTG) для каждой ячейки имеет значение по умолчанию равное 1.
Слайд 6

Цель Grid секции Из уравнения Дарси Это уравнение применимо для горизонтального

Цель Grid секции

Из уравнения Дарси
Это уравнение применимо для горизонтального однофазного течения

в одномерном случае в поверхностных условиях. Член KA/L может быть определен как проводимость между ячейками расчетной сетки (т.е. величина, обратная сопротивлению для течения жидкости между ячейками сетки). ECLIPSE должен рассчитать проводимости в направлениях X, Y и Z для каждой ячейки для расчета потоков в трех направлениях, для чего необходимы проницаемости и размеры ячеек. Так как песчанистость также оказывает влияние на на величину потока между ячейками, то она должна учитываться при расчете проводимостей.
Слайд 7

Цель Grid секции С учетом песчанистости формула для расчета проводимости принимает

Цель Grid секции

С учетом песчанистости формула для расчета проводимости принимает вид:
Например,

для расчета проводимости в направлении X:
Слайд 8

Существует два способа задания геометрии сетки моделирования: блочно-центрированная геометрия и геометрия

Существует два способа задания геометрии сетки моделирования: блочно-центрированная геометрия и геометрия

угловой точки. При задании блочно-центрированной геометрии используются ключевые слова DX, DY, DZ и TOPS.
Ячейки в блочно-центрированной геометрии прямоугольные с горизонтальными верхней и нижней поверхностями и вертикальными сторонами. При задании геометрии угловой точки используются ключевые слова COORD и ZCORN. Ячейки в геометрии угловой точки могут иметь разнообразную форму. Поэтому, с помощью геометрии угловой точки можно довольно легко смоделировать сложные геологические структуры, такие как наклонные разломы, выклинивания и поверхности выветривания.

Цель Grid секции

Слайд 9

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

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

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

Цель Grid секции

Слайд 10

Декартовы координаты • Начало координат находится в верхнем левом углу •

Декартовы координаты
• Начало координат находится в верхнем левом углу
• Значения по

оси X увеличиваются слева направо; по оси Y сверху вниз
• В 3-D, блоки и узлы нумеруются слева на право, сверху вниз
• Начало сеточных координат не обязательно находится в углу ячейки (1, 1, 1)
• Оси X, Y и Z не обязательно параллельны направлениям I, J, K
Радиальные координаты
• Начало координат находится в центре сетки
• Значения увеличиваются в направлении R от центра.
• Угол θ измеряется по часовой стрелке.
• Z направлена вниз (т.е. как и в декартовой СК).

Соглашение на чтение данных

Слайд 11

Соглашение на чтение данных

Соглашение на чтение данных

Слайд 12

Соглашение на чтение данных Соглашение на чтение данных Для чтения данных

Соглашение на чтение данных

Соглашение на чтение данных
Для чтения данных в ячейках

применяется следующее соглашение:
• Первая ячейка в которую записывается или из которой считывается значение имеет координаты (1, 1, 1)
• В Декартовой системе координат ячейка (1, 1, 1) находится в верхнем левом углу графического дисплея
• В радиальной системе координат ячейка (1, 1, 1) находится в центре модели.
• В Декартовой системе координат данных считываются в следующем порядке: быстрее всего изменяется значение X, затем Y, затем Z
• В радиальной системе координат данных считываются в следующем порядке: быстрее всего изменяется значение R, затем THETA, затем Z
• Все симуляторы GEOQUEST используют это соглашение.
Для примера модель 20*5*10 как показано на рисунке, будет определена в секции RUNSPEC. ECLIPSE будет считывать данные о геометрии и свойствах 1000 ячеек для расчета поровых объемов и проводимостей в каждой ячейке. Первые 20 значений пористости будут считаны ECLIPSE начиная с ячейки (1, 1, 1) и заканчивая ячейкой (20, 1, 1), т.е. верхний дальний ряд модели. Следующие 20 значений пористости считаются из ячеек начиная с (1, 2, 1) по (20, 2, 1) включительно, которые составляют следующий после верхнего дальнего ряд модели.
Слайд 13

Соглашение на чтение данных Данные в моделях с радиальной и декартовой

Соглашение на чтение данных

Данные в моделях с радиальной и декартовой системами

координат считываются одинаково в последовательности I, J и K. Однако, в моделях с радиальной системой координат R заменяет X а θ заменяет Y. Z остается неизменной.
Слайд 14

Геометрическое представление • В ECLIPSE геометрия может быть описана двумя способами

Геометрическое представление

• В ECLIPSE геометрия может быть описана двумя способами
• Блочно-центрированная

геометрия (BC)
• Геометрия угловой точки (CP)
• У каждой геометрии есть свои преимущества и недостатки
• Выбор типа геометрии зависит от типа выбранной модели
• Блочно-центрированная геометрия и геометрия угловой точки не могут использоваться вместе, в одном файле данных
Слайд 15

Геометрическое представление Существуют два способа задания размеров и глубин ячеек. Блочно-центрированная

Геометрическое представление

Существуют два способа задания размеров и глубин ячеек.
Блочно-центрированная геометрия (BC)
Для

каждой ячейки требуется задать глубину верхней поверхности и размеры в направлениях X, Y и Z. Верхняя и нижняя поверхности ячеек горизонтальны, стороны вертикальны. Все ячейки имеют форму прямоугольников. Ключевые слова, используемые для задания блочно-центрированной геометрии в Декартовой системе координат, TOPS, DX (или DXV), DY (или DYV) и DZ. Ключевые слова, используемые для задания блочно-центрированной геометрии в радиальной системе координат, DR (или DRV), DTHETA (или DTHETAV) и DZ. Ключевые слова, заканчивающиеся буквой V это векторные аналоги стандартных ключевых слов.
В блочно-центрированной геометрии каждая ячейка описывается всего четырьмя вещественными числами, поэтому описание модели с помощью блочно-центрированной геометрии будет менее объемным, чем описание этой же модели с помощью геометрии угловой точки. Простые модели могут создаваться при помощи блочно-центрированной геометрии без использования программ подготовки данных, например такой как GRID.
Слайд 16

Геометрическое представление

Геометрическое представление

Слайд 17

Геометрическое представление Рассмотрим соседние ячейки в блочно-центрированной геометрии в наклонной структуре,

Геометрическое представление

Рассмотрим соседние ячейки в блочно-центрированной геометрии в наклонной структуре, когда

глубины верхних поверхностей ячеек различаются. Глубины верхних поверхностей ячеек по обе стороны разлома также будут разными. В блочно-центрированной геометрии невозможно отличить наклон от разлома.
Например, на разрезе, представленном на рисунке, нельзя определить без структурной карты, вертикальный разлом там, или наклонная структура. Блочно-центрированная геометрия не содержит информации достаточной для вычисления площади перекрытия между ячейками. Так как не известны глубины углов ячеек.
Обычно между соседними ячейками флюиды перетекают, поэтому соединение должно быть определено и следовательно проводимость должна быть рассчитана. ECLIPSE неявно предполагает, что ячейки имеющие соседние индексы – связаны и соединяются на одном уровне, хотя на самом деле, это не так: ячейки, которые кажутся имеющими общие поверхности в реальности не соединяются. Это результат того, что блочно-центрированная геометрия не учитывет различия между разломами и наклонными пластами. Таким образом, если моделируем наклонную структуру блочной геометрией, то имеем несоседние соединения ячеек (NNC), и проводимость вычисленная для этого «несоседнего» соединения оказывается заниженной. Это решение будет нефизично.
Слайд 18

Геометрическое представление Геометрия угловой точки (СР) Геометрия угловой точки основана на

Геометрическое представление

Геометрия угловой точки (СР)
Геометрия угловой точки основана на идее координатных

линий и глубинах углов ячеек. Координатные линии определяются для каждого края столбца ячеек. Координатные линии всегда прямые, но могут быть и не вертикальными. Каждую координатную линию определяют X, Y и Z координаты одной верхней и нижней точки. Ячейки затем определяются фиксированием их углов вдоль каждой выбранной координатной линии. Это позволяет ячейкам иметь любую физически возможную форму: поэтому могут быть корректно представлены наклонные поверхности, поверхности разломов, зоны выклинивания и поверхности выветривания.
Каждая ячейка определяются четырьмя координатными линиями и восемью глубинами, CP геометрия более объемная чем BC и почти все случаи требуют для конструирования сетки использования пре-процессора, такого как GRID.
CP геометрия содержит достаточную информацию для вычисления перекрытия между прилегающими ячейками так как глубины угловых точек известны. В случае как на рисунке только те ячейки которые находятся внутри сетки - связаны. Итак, жидкости в ячейках показанных стрелками не протекают через разлом.
Слайд 19

Геометрическое представление

Геометрическое представление

Слайд 20

Геометрическое представление

Геометрическое представление

Слайд 21

Пример геометрии центральных блоков • Модель 20*5*10 ячеек • Модель наклонная

Пример геометрии центральных блоков

• Модель 20*5*10 ячеек
• Модель наклонная в двух

направлениях из точки (1,1,1), которая самая
верхняя западная ячейка.
• Блоки величиной 300 ft в X направлении и 1000 ft в Y.
• Слои размерами 32, 22, 20, 4, 32, 4, 26, 26, 4, 28 ft сверху вниз.
Слайд 22

Пример геометрии центральных блоков Следующие ключевые слова содержат полное геометрическое описание

Пример геометрии центральных блоков

Следующие ключевые слова содержат полное геометрическое описание структуры:
TOPS
--Опредедяются

первые 20 TOPS (1, 1, 1) to (20, 1, 1)
6855.000 6865.000 6875.000 6885.000 6895.000
6905.000 6915.000 6925.000 6935.000 6945.000
7005.000 7015.000 7025.000 7035.000 7045.000
7055.000 7065.000 7075.000 7085.000 7095.000
--Cледующие 20 TOPS определяются (1, 2, 1) to (20, 2, 1)
6930.000 6940.000 6950.000 6960.000 6970.000
6980.000 6990.000 7000.000 7010.000 7020.000
7080.000 7090.000 7100.000 7110.000 7120.000
7130.000 7140.000 7150.000 7160.000 7170.000
-- Cледующие 20 TOPS определяются (1, 3, 1) to (20, 3, 1)
7030.000 7040.000 7050.000 7060.000 7070.000
7080.000 7090.000 7100.000 7110.000 7120.000
7180.000 7190.000 7200.000 7210.000 7220.000
7230.000 7240.000 7250.000 7260.000 7270.000
-- Cледующие 20 TOPS определяются (1, 4, 1) to (20, 4, 1)
7130.000 7140.000 7150.000 7160.000 7170.000
7180.000 7190.000 7200.000 7210.000 7220.000
7280.000 7290.000 7300.000 7310.000 7320.000
7330.000 7340.000 7350.000 7360.000 7370.000
-- Cледующие 20 TOPS определяются (1, 5, 1) to (20, 5, 1)
7205.000 7215.000 7225.000 7235.000 7245.000
7255.000 7265.000 7275.000 7285.000 7295.000
7355.000 7365.000 7375.000 7385.000 7395.000
7405.000 7415.000 7425.000 7435.000 7445.000
/ Это завершает TOPS для первого слоя
--Больше TOPS не нужно. Eclipse будет добавлять значение DZ к TOPS для расчета глубины следующего слоя
Слайд 23

DX --Все ячейки имеют DX=300. 1000*300 / DY --Все ячейки имеют

DX
--Все ячейки имеют DX=300.
1000*300 /
DY
--Все ячейки имеют DY=1000
1000*1000 /
EQUALS
--Устанавливается высота DZ

слоя за слоем
--Значения массива I1 I2 j1 j2 k1 k2
'DZ' 32 1 20 1 5 1 1 /
'DZ' 22 1 20 1 5 2 2 /
'DZ' 20 1 20 1 5 3 3 /
'DZ' 4 1 20 1 5 4 4 /
'DZ' 32 1 20 1 5 5 5 /
'DZ' 4 1 20 1 5 6 6 /
'DZ' 26 1 20 1 5 7 7 /
'DZ' 26 1 20 1 5 8 8 /
'DZ' 4 1 20 1 5 9 9 /
'DZ' 28 1 20 1 5 10 10 /
/

Пример геометрии центральных блоков

Слайд 24

Пример представления модели в геометрии угловой точки • Модель из предыдущего

Пример представления модели в геометрии угловой точки

• Модель из предыдущего примера

представлена теперь в CP версии.
• Координатные линии определяют вертикальные или наклонные стороны ячеек
• Соседние ячейки имеют общие координатные линии
• Эту информацию содержит ключевое слово COORD
• Восемь углов для каждой ячейки
• Эту информацию содержит ключевое слово ZCORN
Слайд 25

COORD --Определяет координатные линии -- X1 Y1 Z1 X2 Y2 Z2

COORD
--Определяет координатные линии
-- X1 Y1 Z1 X2 Y2 Z2
0. 0. 6825.000

0. 0. 7023.000
300.0000 0. 6835.000 300.0000 0. 7033.000
600.0000 0. 6845.000 600.0000 0. 7043.000
900.0000 0. 6855.000 900.0000 0. 7053.000
1200.000 0. 6865.000 1200.000 0. 7063.000
..........
..........
..........
/ Для модели размером 20 * 10, требуется 21*11 координатных линий, т.e. 231.
--Каждая линия определяется 6 числами, таким образом всего требуется 1386 чисел.

Пример представления модели в геометрии угловой точки

Слайд 26

ZCORN --Определяет глубины ячеек в порядке возрастания X (или R), затем

ZCORN
--Определяет глубины ячеек в порядке возрастания X (или R), затем Y

(или THETA), затем Z
1953.463745 1953.451538 1953.451538 1953.428467 1953.428467 1953.395020
1953.395020 1953.351440 1953.351440 1953.298706 1953.298706 1953.240845
1953.240845 1953.179932 1953.179932 1953.119385 1953.119385 1953.061646
1953.061646 1953.003540 1953.003540 1952.938965 1952.938965 1952.854248
1952.854248 2262.983154 2262.983154 2262.133057 2262.133057 2261.328613
2261.328613 2260.558350 2260.558350 2259.826172 2259.826172 2259.303955
2259.303955 2258.930664 2258.930664 2258.488770 2258.488770 2258.023926
2258.023926 2257.561768 2257.561768 2257.068359 2257.068359 2256.478516
2256.478516 2255.797363 2255.797363 2255.108154 2255.108154 2254.459961
2254.459961 2253.865479 2253.865479 2253.337646 2253.337646 2252.891113
2252.891113 2252.560791 2252.560791 2252.378906 2252.378906 2252.331055
2252.331055 2252.474121 2252.474121 2252.915771 2252.915771 2253.517822
2253.517822 2254.188232 2254.188232 2254.982666 2254.982666 2255.978271
2255.978271 2257.189941 2257.189941 2258.449463 2258.449463 2259.707275
2259.707275 2260.960938 2260.960938 2262.190674 2262.190674 1949.721802
1949.721802 1949.695190 1949.695190 1949.675781 1949.675781 1949.671509
1949.671509 1949.679565 1949.679565 1949.693481 1949.693481 1949.710693
1949.710693 1949.721069 1949.721069 1949.725098 1949.725098 1949.721558
............
............
...........
/ Для 1000 cells, требуется 8000 значений ZCORN.
--Здесь приведено только несколько

Пример представления модели в геометрии угловой точки

Слайд 27

• Параметры которые описывают размер и глубину каждой ячейки относятся к

• Параметры которые описывают размер и глубину каждой ячейки относятся к

геометрическим параметрам
• Ключевые слова используемые для задания геометрических параметров: TOPS, DX (or DXV), DY (or DYV), DZ для декартовой геометрии и TOPS, DR (or DRV), DTHETA (or DTHETAV) и DZ для радиальной геометрии, или для геометрии угловой точки - COORD и ZCORN.
• Параметры которые описывают пористость, проницаемость и т.д. относятся к свойствам.
• Ключевые слова используемые для задания свойств это: PORO (φ), PERMX (Kx), PERMY (Ky), PERMZ (Kz).
• Песчанистость задается используя либо NTG (коэффициент песчанистости) либо DZNET (толщина слоя).
• Ключевое слово используемое для указания активная ячейка или нет ACTNUM.
Равно 0 для неактивных ячеек, 1 для активных.

Свойства ячеек сетки

Слайд 28

Ключевые слова для задания свойств нуждаются в одном значении для каждой

Ключевые слова для задания свойств нуждаются в одном значении для каждой

ячейки. Есть много различных способов задания этих значений используя различные комбинации ключевых слов ECLIPSE. Значения свойств определяются в центрах ячеек и значение свойства распространяется на весь объем ячейки. Осреднение производится до введения данных в ECLIPSE, осредненное значение определяется пользователем. Заметим также что хотя некоторые ячейки исключены из моделирования (неактивные ячейки), должны быть предоставлены данные для расчета их поровых объемов и проводимостей. Это вызвано тем, что ECLIPSE содержит устройство для отключения ячеек из моделирования, основанное на расчете минимального порового объема ячеек. В ECLIPSE данные задаются только в виде явных значений. ECLIPSE не имеет устройства для введения данных в форме функций. Для примера пористость не может быть определена как функция проницаемости. Значения пористости
генерируются пре-процессором, таким как GRID, и данные для каждой ячейки импортируются в ECLIPSE.

Свойства ячеек сетки

Слайд 29

Как назначаются свойства ячейкам сетки

Как назначаются свойства ячейкам сетки

Слайд 30

Как назначаются свойства ячейкам сетки Как устанавливается одно значение свойства для

Как назначаются свойства ячейкам сетки

Как устанавливается одно значение свойства для каждой

ячейки сетки
Стандартный формат
KEYWORD
Value1 value2 value3……..value(NX*NY*NZ) /

К примеру,
PERMX
--K=1
10 10 10 10 10 10 10 10 10 10
10 10 10 10 10 10 10 10 10 10
--k=2
5 5 5 5 5 5 5 5 5 5
5 5 5 5 5 5 5 5 5 5
--k=3
100 100 100 100 100 100 100 100 100 100
100 100 100 100 100 100 100 100 100 100
--k=4
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
k=5
200 200 200 200 200 200 200 200 200 200
200 200 200 200 200 200 200 200 200 200
--k=6
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
--k=7
100 100 100 100 100 100 100 100 100 100
100 100 100 100 100 100 100 100 100 100
--k=8
50 50 50 50 50 50 50 50 50 50
50 50 50 50 50 50 50 50 50 50
--k=9
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
--k=10
50 50 50 50 50 50 50 50 50 50
50 50 50 50 50 50 50 50 50 50

Размерность модели 20*1*10

Слайд 31

Как назначаются свойства ячейкам сетки Значения для нескольких ячеек могут быть

Как назначаются свойства ячейкам сетки

Значения для нескольких ячеек могут быть определены

при помощи *,
к примеру
PERMX
20*10
20*5
20*100
20*2000
20*200
20*2000
20*100
20*50
20*2000
20*50
/

Размерность модели 20*1*10

Слайд 32

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

Как установить значения для ячеек сетки используя боксы

Вводимый бокс определяется значениями

I, J и K и относится к прямоугольному блоку ячеек. Бокс может быть участком, слоем, колонкой или рядом ячеек. В боксе могут быть установлены любые свойства ячеек сетки. Значение должно быть определено для всех ячеек бокса, включая находящиеся в на краях бокса. Бокс остается задействованным пока не начинается считывание нового бокса или не считано ключевое слово ENDBOX. Новый бокс закрывает предыдущий бокс и открывает следующий. По умолчанию, бокс включает все ячейки.
Например:
BOX
--I1 I2 J1 J2 K1 K2
1 20 1 1 1 10 /
PERMX
200*100/
ENDBOX
PERMX устанавливается равным 100 для ячеек от (1,1,1) до (20, 1, 10) т.е. для всей
модели.

Размерность модели 20*1*10

Слайд 33

Как установить значения для ячеек сетки используя боксы Меньшие боксы могут

Как установить значения для ячеек сетки используя боксы

Меньшие боксы могут использоваться

для установки проницаемости
послойно: Например, для слоя 3
BOX
--I1 I2 J1 J2 K1 K2
1 20 1 1 3 3 /
PERMX
20*200/
ENDBOX
Это изменяет значения PERMX установленные для всего слоя 3 на 200.

Размерность модели 20*1*10

Единицы измерения: METRIC: mD, FIELD: mD, LAB: mD, PVT-M: mD

Слайд 34

Как установить значения для ячеек сетки используя EQUALS Ключевое слово EQUALS

Как установить значения для ячеек сетки используя EQUALS

Ключевое слово EQUALS работает

со свойствами как с массивом и может служить альтернативой ключевому слову BOX. EQUALS можно использовать в последнем примере для установки PERMX в слое 3 :
BOX
--I1 I2 J1 J2 K1 K2
1 20 1 1 3 3 /
EQUALS
‘PERMX’ 100 /
ENDBOX

Размерность модели 20*1*10

Слайд 35

Как установить значения для ячеек сетки используя EQUALS EQUALS также может

Как установить значения для ячеек сетки используя EQUALS

EQUALS также может использоваться

для определения области без использования ключевого слова BOX.
Например:
EQUALS
--Arrayvalue I1 I2 J1 J2 K1 K2
‘PERMX’2000 / defaults to currently open box, i.e. entire reservoir
‘PERMX’10 1 20 1 5 1 1 /
‘PERMX’5 1 20 1 5 2 2 /
‘PERMX’100 1 20 1 5 3 3 /
‘PERMX’200 1 20 1 5 5 5 /
‘PERMX’100 1 20 1 5 7 7 /
‘PERMX’50 1 20 1 5 8 8 /
‘PERMX’50 1 20 1 5 10 10 /
/
определяет PERMX во всей секции. Этот бокс открывается ключевым словом EQUALS и закрывается вторым прямым слэшем (/) который служит концом ключевого слова EQUALS. Для уверенности в том, что значения определены для всех ячеек мы рекомендуем определить сначала значения для всей сетки, переприсваивая их затем выборочно используя BOX и EQUALS.

Размерность модели 20*5*10

Слайд 36

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

Как складывать, вычитать, умножать и делить данные по ячейкам

Обычный формат для

этих действий

‘KEYWORD’ Value I1 I2 J1 J2 K1 K2 /
‘KEYWORD’ Value I1 I2 J1 J2 K1 K2 /
‘KEYWORD’ Value I1 I2 J1 J2 K1 K2 /
/
Или
BOX
I1 I2 J1 J2 K1 K2 /

‘Keyword’ Value /
‘Keyword’ Value /
‘Keyword’ Value /
ENDBOX
Слайд 37

Может быть откорректировано любое количество массивов. Ограничение боксом не обязательно (опционально).

Может быть откорректировано любое количество массивов. Ограничение боксом не обязательно (опционально).

Если корректируется свойства всей модели, то в первом примере не указываются границы, во втором не применяются BOX и ENDBOX.
К примеру, проницаемость не всегда однородна и часто представляется как Kv/Kh =0.1 (вертикальная/горизонтальная). Это можно определить для всего поля следующим образом (сначала надо присвоить значение ‘PERMX’ ‘PERMZ’ ):
MULTIPLY
--Keyword Value I1 I2 J1 J2 K1 K2
‘PERMZ’ 0.1 / не выбираются пределы – исключены предыдущим открытым боксом
/ --этот последний слеш заканчивает ключевое слово MULTIPLY

Как складывать, вычитать, умножать и делить данные по ячейкам

Слайд 38

К примеру для установки соотношения Kv/Kh только для первого слоя используется:

К примеру для установки соотношения Kv/Kh только для первого слоя используется:
MULTIPLY
--Keyword

Value I1 I2 J1 J2 K1 K2
‘PERMZ’ 0.1 1 20 1 5 1 1 /
/
Или
BOX
--I1 I2 J1 J2 K1 K2
1 20 1 5 1 1 /
MULTIPLY
--массив значение
‘PERMZ’ 0.1 /
ENDBOX

Как складывать, вычитать, умножать и делить данные по ячейкам

Слайд 39

Для деления на какое-то число нужно умножить на обратную величину. Эти

Для деления на какое-то число нужно умножить на обратную величину. Эти

действия обратные. Таким образом, умножение на 10 после умножения на 2 приводит в результате к умножению на 20. Сложение происходит аналогично. К примеру, для увеличения пористости на 0,05 в той же области, используется:
ADD
--Ключевое слово Значение I1 I2 J1 J2 K1 K2
‘PORO’ 0.05 1 20 1 5 1 1 /
/ этот слеш завершает действие ключевого слова ADD
Или
BOX
--I1 I2 J1 J2 K1 K2
1 20 1 5 1 1 /
ADD
--Массив Значение
‘PORO’ 0.05 /
ENDBOX
Для вычитания добавьте отрицательное значение.

Как складывать, вычитать, умножать и делить данные по ячейкам

Слайд 40

Объем пор ячейки может быть увеличен с использованием MULTPV. Например, MULTPV

Объем пор ячейки может быть увеличен с использованием MULTPV.
Например,
MULTPV
200*1.01 /
увеличит

объем пор первых 200 ячеек на 1%. Это может быть необходимо, к примеру, для того чтобы удостовериться, что моделирование объема флюида в пластовых условиях (запасов, FIP – fluid in place) соотносится с оценкой из других источников.
Обратите внимание, чем больше объем пор, тем выше уровень поддержки давления в
ячейке (это происходит за счет сжимаемости флюида). Объем пор, таким образом, не может сильно варьироваться, а согласование изменений общего давления в пласте может проводиться только после окончания манипуляций с объемом пор. Изменение объема пор происходит в процессе согласования для достижения одной из или обоих целей.

Как умножать значения пористости ячейки, используя MULTPV

Слайд 41

Можно перенести данные из одной части сетки в другую, используя COPYBOX.

Можно перенести данные из одной части сетки в другую, используя COPYBOX.

Исходная и конечная области должны быть одного размера по координатам I, J,K.
Основная форма:
COPYBOX
Keyword Source box Destination box /
Ключевое слово Бокс источник Бокс получатель/
/
К примеру, для переноса значений пористости и проницаемости по координате х из слоя 1 в слой 2 используется:
COPYBOX
--Keyword I1S I2S J1S J2S K1S K2S
PORO 1 20 1 5 1 1
--I1D I2D J1D J2D K1D K2D
1 20 1 5 2 2 /
PERMX 1 20 1 5 1 1
1 20 1 5 2 2 /
/

Как копировать данные из одной части сетки в другую, используя COPYBOX

Слайд 42

Данные могут быть считаны прямо в файл данных ECLIPSE из внешнего

Данные могут быть считаны прямо в файл данных ECLIPSE из внешнего

файла, с использованием INCLUDE. Файл должен быть расположен в файловой системе и к нему можно обратиться по абсолютному или относительному пути. Не надо указывать путь, если файл расположен в той же директории, что и файл данных ECLIPSE. Все приводимые примеры имеют смысл при существовании файлов и путей к ним:
INCLUDE
COORD.GRDECL /
INCLUDE
ZCORN.GRDECL /
INCLUDE
‘../../GRID/PERMX.GRDECL’ /
INCLUDE
‘/tiny/user/hm2/griddata/outputs/PORO.GRDECL’ /
INCLUDE используется для ввода большого количества данных. Кавычки используются если путь содержит символы, которые ECLIPSE может неоднозначно интерпретировать, к примеру прямой слеш и точка. В Unix-е обратите внимание на прописные буквы в имени файла. Вводимые файлы могут быть вложенными, т.е. содержать внутри себя INCLUDE. INCLUDE не может применяться с другими ключевыми словами.

Как считывать данные из другого файла,
используя INCLUDE

Слайд 43

Любую ячейку можно сделать не активной, используя ACTNUM. ACTNUM в массиве

Любую ячейку можно сделать не активной, используя ACTNUM. ACTNUM в массиве

свойств ячеек сетки может использоваться во всех приведенных примерах. ACTNUM может принимать только два значения 0 для неактивной, и 1 для активной ячейки. Препроцессор, такой как GRID часто используется для определения ACTNUM для включения в ECLIPSE.
Неактивные ячейки исключаются из моделирования. ECLIPSE не рассчитывает потоки в этих ячейках. Графический вывод, однако, сильно не изменяется, т.к. расположение ячеек остается неизменным. ECLIPSE должен иметь достаточно информации, для расчета объема пор, мощности и проводимости неактивных ячеек (хотя они и неактивные).
Неактивные ячейки отображаются при выводе. Данных, связанных с этими ячейками очень мало, т.к. большинство параметров не определены; отсутствующие значения заменяются несколькими прочерками.

Как деактивировать ячейки,
используя ACTNUM

Слайд 44

Практическое занятие Исходные данные: Месторождение: Северо-Морское нефтяное Размерность модели 15х8х5 Ввод

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

Исходные данные:
Месторождение: Северо-Морское нефтяное
Размерность модели 15х8х5
Ввод в разработку 1 января

2012
Система единиц измерения российская
Геометрия угловой точки
(указать необходимые ключевые слова и количество значений для данной модели)
Задать песчанистость (70% для 1 и 5 слоев, 75% для 2 и 3 слоев и 65% для 4 слоя),
пористость (1,4,5 слои – 15%, 2 слой – 16%, 3 слой – 17%),
проницаемость в направлениях X и Y равна и составляет 15 мД для 1 слоя, 25 мД – для 2 и 3 слоев и 10мД для 4 и 5 слоев,
проницаемость по оси Z – 15% от проницаемости по горизонтали
Сделать:
секции Runspec и Grid по заданным параметрам
Слайд 45

Назначение секции RUNSPEC

Назначение секции RUNSPEC

Слайд 46

• Проводимость – свойство, определяемое для смежных ячеек. • Проводимость определяет

• Проводимость – свойство, определяемое для смежных ячеек.
• Проводимость определяет перетоки

флюидов между ячейками.
• Для правильного изменения проводимостей важно знать, как они определены.
• Проводимость в ECLIPSE между двумя ячейками рассчитывается в положительном направлении, т.е. в направлении вверх по потоку.
• Tx i,j,k между ячейками (I, J, K) и (I+1, J, K) в направлении I
• Tyi,j,k между ячейками (I, J, K) и (I, J+1, K) в направлении J
• Tzi,j,k между ячейками (I, J, K) и (I, J, K+1) в направлении K
• На рисунке приведен пример направлений соединений
• Направление соединения представляет собой направление расчета проводимостей между ячейками с соседними IJK индексами.
• Каждая ячейка имеет 6 направлений соединений. Ячейка (I, J, K) соединена с
(I-1, J, K), (I+1, J, K), (I, J-1, K), (I, J+1, K), I, J, K-1), (I, J, K+1).

Характеристики проводимости

Слайд 47

Характеристики проводимости

Характеристики проводимости

Слайд 48

Проводимость – свойство, определяемое для смежных ячеек, т.е. ячеек, между которыми

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

течение флюидов. Переток флюида расчитывается между
центрами ячеек. Величина этого перетока определяется проводимостью и
подвижностью флюида между соседними ячейками.
Рассмотрим переток между центрами двух соседних ячеек. Проводимость должна учитывать свойства каждой ячейки, т.е. должна быть своего рода средним от свойств обеих ячеек, также учитывая геометрию ячеек и площадь их совместной поверхности.
Варианты расчета проводимости, используемые ECLIPSE будут описаны в следующей лекции. Некоторые из них больше подходят для определенной геометрии ячеек, чем другие. Все расчеты проводимости, однако, ведутся в направлении, вверх по потоку т.е. на определение проводимостей ячейки (I, J, K) влияют потоки к ячейкам (I+1, J, K), (I, J+1, K) и (I, J, K+1).

Характеристики проводимости

Слайд 49

Распределение подвижности флюида между соседними ячейками не зависит от проводимости. В

Распределение подвижности флюида между соседними ячейками не зависит от проводимости. В

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

Характеристики проводимости

Слайд 50

Характеристики проводимости • Проводимость рассчитывается по-разному в геометрии блочно- центрированой и

Характеристики проводимости

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

точки
• Радиальная и декартовая проводимость также рассчитываются по-разному.
• Методы расчета запускаются с использованием ключевых слов OLDTRAN
(блочно-центрированная), NEWTRAN (угловая), OLDTRANR (блочно-центрированная)
• Если OLDTRAN, OLDTRANR или NEWTRAN не определены, ECLIPSE использует OLDTRAN для блочно-центрированной геометрии и NEWTRAN для геометрии угловой точки автоматически.
• OLDTRAN рассчитывается как (среднегармоническая проницаемость)
*(среднеарифметическую площадь)
• OLDTRANR - как среднегармоническое от (проницаемость*площадь)
• NEWTRAN – среднегармоническое проводимости половины ячейки
• Радиальная проводимость одинаково рассчитывается в блочно- центрированной и геометрии угловой точки
• Выражения основаны на радиальном потоке между радиусами с одинаковым давлением
Слайд 51

Моделирование глин • Каждый литологический слой может быть смоделирован как отдельный

Моделирование глин

• Каждый литологический слой может быть смоделирован как отдельный слой

сетки
• Глины могут объединены со слоями песка
• Слои глин могут быть смоделированы как разрывы между слоями сетки
Слайд 52

Моделирование глин Моделирование глин явно, как слоев сетки Этот способ способствует

Моделирование глин

Моделирование глин явно, как слоев сетки
Этот способ способствует контролю свойств

каждого слоя. Изменяя проницаемость, которая влияет на проводимость, можно контролировать поток, проходящий через глины. Хотя NTG глинистого барьера часто снижается до нуля, NTG не участвует в расчетах вертикальной проводимости. Однако, влияет на расчет горизонтальной проводимости, таким образом горизонтальный поток в глинах отсутствует.
Т.к. глинистые пласты очень тонкие, объем пор ячеек глинистого барьера сравнительно ниже, чем объем пор соседних ячеек.
Возможные трудности:
• Если модель использует ключевые слова PINCH и/или MINPV для деактивации
ячеек и установления связей через выклинивающиеся ячейки, ячейки, содержащие глины могут быть исключены из моделирования. Ячейки, деактивированные MINPV, становятся барьерами на пути потока. Ячейки тоньше пороговой величины исключаются из процесса моделирования и ячейки с разных сторон слоя взаимодействуют так, как будто глины отсутствуют. Таким образом, глины могут представляться в виде барьера на пути потока, или игнорироваться.
Слайд 53

Моделирование глин Моделирование глин явно, как слоев сетки Возможные трудности: •

Моделирование глин

Моделирование глин явно, как слоев сетки
Возможные трудности:
• Соседство ячеек с

большим и малым объемом пор может приводить к
проблемам схождения пропускной способности. Рассмотрим две соседние
ячейки с объемом пор 1 м3 и 1000 м3. ECLIPSE рассчитывает
нефтенасыщенность с точностью ±0.001. Если нефтенасыщенность
рассчитана как 0.5, то запасы нефти в пласте (OIP) в каждой ячейке составляют
0.5±0.001 м3 и 500±1 м3. OIP одной ячейки меньше чем погрешность
измерения OIP другой, т.е. наибольшая ячейка определяет пропускную
способность меньшей. Это может привести к проблемам схождения, чего
следует избегать.
• Если существует большое количество протяженных глинистых слоев, их
моделирование может потребовать введения неразумно большого количества
ячеек, что приведет к перерасходу времени.
Слайд 54

Моделирование глин Моделирование глин путем включения в более крупные песчаные ячейки

Моделирование глин

Моделирование глин путем включения в более крупные песчаные ячейки
Глины могут

объединятся с песчаными слоями. Для этого необходимо ввести чистую толщину песка, используя DZNET или NTG где

Отметим, что присваиваемая пористость ячейки должна быть равна пористости чистого песка, а не средней пористости всей ячейки.
Этот метод «размазывает» глину по всей ячейке, но распределение объема пор по высоте является некорректным. Также, вертикальная проводимость между геологическим слоем глин и соседним песчаным слоем больше не равна нулю. Горизонтальная проводимость в глине так же не нулевая. Горизонтальная проводимость в песке снижается. Все эти проводимости, в особенности в вертикальном направлении могут быть пересчитаны с использованием одного или нескольких ключевых слов “MULT”. Что позволяет решить проблему схождения пропускной способности и уменьшить количество ячеек в моделировании. Также снижаются или исчезают полностью сложности, связанные с тем, что неизвестно, какие именно ячейки в ECLIPSE будут деактивированы ключевыми словами MINPV или MINPVV.

Слайд 55

Моделирование глин Моделирование глин как промежутков между слоями песка Использование некоторых

Моделирование глин

Моделирование глин как промежутков между слоями песка
Использование некоторых препроцессоров для

подготовки сетки, таких как GRID, позволяет отделять один слой от другого. Пространство между слоями K=1 и K=2 не заполняется ячейками; оно пустое. Вертикальная проводимость между слоями K=1 и K=2, однако, не нулевая. Она рассчитывается так, как укажет пользователь, и разрыв не оказывает на это влияния. Таким образом глинистый слой исключается, а вертикальная проводимость между слоями может регулироваться MULTZ и/или MULTZ- увеличивая или уменьшая вертикальный поток. Горизонтальная проводимость и распределение порового объема в слоях 1 и 2 такие же, как если бы глины были явно включены в модель – т.е. не изменяется. Это позволяет решить проблему схождения пропускной способности и уменьшить количество моделируемых ячеек и изменять вертикальный поток между слоями вне зависимости от пористости и проницаемости глин.
Слайд 56

Изменение проводимости • Проводимость часто могут изменять (при адаптации модели) •

Изменение проводимости

• Проводимость часто могут изменять (при адаптации модели)
• Вверх по

ходу движения потока проводимость может быть явно установлена с использованием TRANX/TRANY/TRANZ, TRANR/TRANTHT (в разделе EDIT)
• Вниз по ходу движения потока проводимость может быть установлена явно используя TRANX-/TRANY-/TRANZ-, TRANR-/TRANTHT- (в разделе EDIT)
• Вверх по ходу движения потока проводимость может быть умножена с использованием MULTX/Y/R/THT/Z
• Вниз по ходу движения потока проводимость может быть умножена с использованием MULTX-/Y-/R-/THT-/Z-
• Коэффициент потоковой проводимости применяется при использовании GRIDOPTS в RUNSPEC
• Множитель проводимости не влияет на проницаемость
• Разлом может быть задан с использованием ключевого слова FAULTS
• Проводимость разлома вверх по ходу движения потока может быть умножена с
использованием ключевого слова MULTFLT
• Коэффициенты проводимости не накапливаются; они влияют на прямые и несоседние соединения
Слайд 57

Изменение проводимости Проводимость может быть изменена несколькими способами. Один из самых

Изменение проводимости

Проводимость может быть изменена несколькими способами. Один из самых простых

– изменение входной проницаемости. Это может быть сделано изменением проницаемости с использованием “PERM” или умножением проницаемости некоторый известный коэффициент, используя MULTIPLY. Однако, этот способ можно использовать не всегда.
К примеру, если глинистый барьер неопределенного размеров и проницаемости моделируется как разрыв между двумя пластами, а не как ячейки, должна быть задана проводимость через разрыв. Она может быть задана напрямую, через TRANZ и/или TRANZ- или умножением, используя MULTZ и/или MULTZ-.
Слайд 58

Действие ключевых слов изменения проводимости “TRAN” и “MULT” влияют на проводимость

Действие ключевых слов изменения проводимости

“TRAN” и “MULT” влияют на проводимость по

ходу течения потока,
“TRAN”- и “MULT”- в направлении против течения потока
Два последних могут использоваться только если установлен соответствующий параметр в GRIDOPTS в разделе RUNSPEC
MULTX, MULTY, MULTZ изменяют проводимость по ходу течения потока по направлениям X, Y и Z соответственно
MULTX-, MULTY-, MULTZ- изменяют проводимость против течения потока по направлениям X, Y и Z соответственно
MULTR, MULTTHT изменяют проводимость по ходу течения потока по радиусу и азимуту (по часовой стрелке)
MULTR-, MULTTHT- изменяют проводимость против течения потока по радиусу и азимуту (против часовой стрелки)
Взаимодействие между множителями проводимости может быть комплексным.
Слайд 59

Действие ключевых слов изменения проводимости Необходимо учитывать следующее: • Все ключевые

Действие ключевых слов изменения проводимости

Необходимо учитывать следующее:
• Все ключевые слова,

приведенные на предыдущем слайде, используются как в блочно-центрированной геометрии, так и в геометрии угловой точки
• MULTZ и MULTZ- используются в радиальной и в декартовой системах координат
• Множители “MULT” не накапливаются
• Если MULTIPLY используется для постоянного изменения проводимости, то оно накапливается
• Множители “MULT” влияют на прямые или нормальные соединения
• Множители “MULT” не влияют на несмежные соединения, созданные с использованием NNC (явные NNCs созданные пользователем)
• Множители “MULT” влияют на несмежные соединения, созданные ECLIPSE, например, соединения через разлом
• Если оба множителя “MULT”, по потоку и против потока определены между двумя ячейками, результат умножения относится к проводимости. К примеру, если MULTX задан для ячейки (I, J, K) а MULTX- для (I+1, J, K) рассчитанная проводимость между ячейками умножается на результат MULTX и MULTX-
• Множитель “MULT”, включенный в раздел GRID, влияет только на проводимость, рассчитанную в этом разделе
• Множитель “MULT”, включенный в раздел EDIT влияет на расчет проводимости, сделанный в разделе GRID
Слайд 60

Изменение проводимости

Изменение проводимости

Слайд 61

Несоседнее соединение (NNC) Обычно течение ожидается между соседними ячейками, и моделируется

Несоседнее соединение (NNC)

Обычно течение ожидается между соседними ячейками, и моделируется как

прямое соединение между ячейками с соседними индексами (I, J, K). Точное определение многих структур, однако, часто требует, чтобы ячейки с несмежными индексами (I, J, K) были бы расположены рядом друг с другом и тогда между ними существует поток. Примером является поток через разлом с очень большим сдвигом. Другие ситуации, требующие использования NNCs в геометрических описаниях, используемых ECLIPSE – локальные измельчения сетки (LGR) и подключение подошвенных вод.
Слайд 62

Несоседнее соединение (NNC) NNC допускает существование перетоков флюидов между ячейками с

Несоседнее соединение (NNC)

NNC допускает существование перетоков флюидов между ячейками с несоседними

индексами IJK
NNC необходимо если ожидаются перетоки между несоприкасающимися ячейками
Создание NNC возможно по умолчанию
Для отключения этой функции используйте NONNC в RUNSPEC
NNCs и проводимость между ними создаются автоматически в следующих случаях:
• При наличии разломов, для создания которых использовался NEWTRAN (геометрия угловой точки)
• При присутствии выклиниваний и поверхностей выветривания, созданных при помощи ключевых слов PINCH или PINCHOUT
• В локальных измельчениях и загрублениях сетки
• В моделях двойной пористости
В радиальных моделях с использованием ключевого слова COORDSYS явно задаются проводимости NNC «для замыкания окружности»
Для водоносного горизонта проводимости NNC должны быть определены явно
Кроме того, NNC могут быть созданы между любой парой ячеек сети в моделях, использующих ключевое слово NNC. Проводимость таких NNC должна задаваться явно
Слайд 63

Несоседнее соединение (NNC)

Несоседнее соединение (NNC)

Слайд 64

Создание несоседнего соединения через разлом • Разломные NNC создаются автоматически при

Создание несоседнего соединения через разлом

• Разломные NNC создаются автоматически при отсутствии

NONNC в разделе RUNSPEC
• Разломные NNC создаются при использовании NEWTRAN
• Значения проводимости вычисляются автоматически
Слайд 65

Создание несоседнего соединения через разлом Если NNC не запрещены ключевым словом

Создание несоседнего соединения через разлом

Если NNC не запрещены ключевым словом NONNC,

то при наличии разломов они создаются, если это необходимо и возможно.
В геометрии угловой точки, где NEWTRAN используется по умолчанию NNC создаются автоматически, без участия пользователя. Использование OLDTRAN в геометрии угловой точки не рекомендуется, так как это ключевое слово запрещает создание разломных NNC.
В блочно-центрированной геометрии, где OLDTRAN используется по умолчанию, ECLIPSE не хватает геометрической информации о том, связаны ли ячейки через разлом и разломные NNCs не создаются. Хотя в блочно- центрированной геометрии возможно использование NEWTRAN, это не дает дополнительной геометрической информации.
Слайд 66

Создание несоседних соединений через выклинивание • ECLIPSE позволяет слоям иметь нулевую

Создание несоседних соединений через выклинивание

• ECLIPSE позволяет слоям иметь нулевую толщину,

что позволяет моделировать зоны выклинивания и поверхности выветривания (когда слой нулевой толщины - имеем несоседнее соединение)
• PINCH создает NNC через слои с DZ меньше, чем заданное пороговое значение, если NONNC не применялся (т.е. если толщина ячейки меньше порогового значения имеем несоседнее соединение)
• PINCH содержит переключатели для создания NNC через слои с DZ больше пороговых значений, если ячейки деактивированы MINPV
Слайд 67

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

Создание несоседних соединений через выклинивание

Геометрия угловой точки позволяет ячейкам иметь непрямоугольную

форму, что является мощным средством для точного моделирования выклиниваний и поверхностей выветривания. Хотя эти структуры очень различаются геологически, с точки зрения геометрии сетки моделирования они определяются сходным образом. NNC создаются через неактивные ячейки с толщиной меньше заданного значения при использовании PINCH. К примеру,
PINCH
0.01 /
создаст NNC между ячейками расположенными выше и ниже ячейки с DZ (толщиной ячейки) меньше 0.01 фута или метра в зависимости от установленной системы исчисления.
Однако, определенное количество ячеек будут автоматически деактивированны ECLIPSE. Для того, чтобы избежать проблем сходимости, ECLIPSE деактивирует все ячейки с объемом пор меньше, чем 10-6 м3 или Rb (reservoir barrel). Это ограничение может быть увеличено, но не уменьшено, с использованием MINPV. Например,
MINPV
5000 /
будет деактивировать все ячейки с объемом пор меньше 5000 Rb в модели, использующей единицы FIELD. Эти ячейки не используются в моделировании и при решении уравнений. По умолчанию они не являются барьерами на пути потока и нет гарантии, что PINCH создаст NNC через них. Если второй пункт PINCH установлен на 'GAP', ECLIPSE создает NNC через ячейки,
деактивированные MINPV.
Слайд 68

Создание несоседних соединений в местах локальных измельчений сетки (LGR) • ECLIPSE

Создание несоседних соединений в местах локальных
измельчений сетки (LGR)

• ECLIPSE автоматически создает

NNC между основными ячейками и измельченными ячейками в местах локального измельчения сетки.
Слайд 69

Создание несоседних соединений в местах локальных измельчений сетки (LGR) Ячейки, входящие

Создание несоседних соединений в местах локальных
измельчений сетки (LGR)

Ячейки, входящие в локальное

измельчение сетки (LGR) используют индексы I, J и K локализованные в каждой LGR. Ячейки LGR редко будут иметь индексы смежные с соседними глобальными ячейками. При моделировании с использованием LGR, ECLIPSE сообщает о NNC между ячейками LGR и соседними основными ячейками. Это не совсем NNC; а скорее удобная форма сообщения, т.к. LGR рассчитывается отдельно от глобальной сети и не существует прямого потока между LGR и основной сеткой.
Пользователь не может изменить такие NNC.
Слайд 70

Несоседние соединения в моделях с двойной пористостью • ECLIPSE удваивает количество

Несоседние соединения в моделях с двойной пористостью

• ECLIPSE удваивает количество слоев

в моделях с двойной пористостью
• Верхняя половина – ячейки матрицы, нижняя половина – ячейки трещин
• Матричные и соответствующие ячейки разломов автоматически соединяются NNC
Слайд 71

Несоседние соединения в моделях с двойной пористостью В пластах с двойной

Несоседние соединения в моделях с двойной пористостью

В пластах с двойной пористостью

флюиды присутствуют в двух взаимосвязанных системах:
В матрице (теле) горной породы (rock matrix), которая обычно обеспечивает большую часть вмещающего объема пласта-коллектора
В высокопроницаемых трещинах (fractures), рассекающих тело горной породы
Если ячейки матрицы связаны только через систему трещин, то данная система называется системой двойной пористости - одинарной проницаемости (dual porosity single permeability system), в данном случае фильтрация флюида происходит только через систему трещин, а содержится флюид в основном в ячейках матрицы.
Если же существует и возможность перетока флюидов напрямую между ячейками матрицы, миную систему трещин, то данная система называется системой двойной пористости - двойной проницаемости (dual porosity dual permeability system)
Слайд 72

Несоседние соединения в моделях с двойной пористостью Пласты с двойной пористостью

Несоседние соединения в моделях с двойной пористостью

Пласты с двойной пористостью и

двойной проницаемостью моделируются в ECLIPSE путем разделения матрицы и трещин на отдельные группы ячеек, соединенных автоматически создаваемыми NNC.
Пользователь может создать сетку пласта как обычную. Затем использовать ключевое слово DPGRID для того, чтобы дать команду ECLIPSE удвоить количество слоев в модели и копировать свойства каждого дополнительного слоя с уже существующего. Верхние слои назовем ячейками матрицы, а нижние – ячейками трещин. Пользователь должен ввести данные по свойствам ячеек трещин (пористость), коэффициент перетока, связывающий матрицу и трещины и функции насыщения, как для матричной ячейки, так и для трещин.
В моделях с двойной пористостью не существует перетока между ячейками матрицы, поэтому перфорация скважины должна попадать в нижнюю часть модели (т.е. быть связанной с трещиной). В моделях с двойной проницаемостью поток матрица-матрица имеет важное значение, а перфорации скважин должны указываться в ячейках трещин и матрицы.
Слайд 73

Создание несоседних соединений для подошвенных вод • Водоносные горизонты (aquifer) подсоединяют

Создание несоседних соединений
для подошвенных вод

• Водоносные горизонты (aquifer) подсоединяют к

пласту коллектору при помощи NNC
• Местоположение данных NNC должно быть задано явно
• Проводимость NNC может быть изменена
Слайд 74

Создание несоседних соединений в радиальных моделях • Ячейки, имеющие границы с

Создание несоседних соединений
в радиальных моделях

• Ячейки, имеющие границы с θ

= 0° и 360° несмежные
• Для создания NNC необходимо включить несоседнее соединение и «замкнуть
окружность»
Слайд 75

Создание несоседних соединений в радиальных моделях Как правило, поток в окрестностях

Создание несоседних соединений
в радиальных моделях

Как правило, поток в окрестностях скважины

не обязательно направлен радиально к стволу скважины. Для полного замыкания тангенсальной составляющей потока, должно быть создано NNCs между ячейками с гранями 0° и 360°. Это относится к блочно-центрированной радиальной модели и к радиальной модели угловой точки. Создание NNC между этими ячейками включается ключевым словом COORDSYS. Для замыкания окружности в верхних четырех слоях шестислойной модели используется
COORDSYS
--K1 K2 Complete?
1 4 COMP /
5 6 INCOMP /
/
При использовании ключевого слова NONNC в RUNSPEC данное действие будет запрещено.
Слайд 76

Радиальные модели • Для включения радиальной функции используйте ключевое слово RADIAL

Радиальные модели

• Для включения радиальной функции используйте ключевое слово RADIAL в
RUNSPEC
•

Значения I,J,K относятся к радиальному, тангенсальному и вертикальному
направлениям
• Тангенсальные расстояния измеряются по часовой стрелке в градусах во всех
системах единиц (metric, field и т.д.)
• Внутренний радиус задается ключевым словом INRAD
• Количество ячеек в радиальном направлении задается NR в RUNSPEC
• Радиальные размеры ячейки устанавливаются с помощью DR (или DRV) или
рассчитываются ECLIPSE если определен внешний радиус (OUTRAD)
• Тангенсальные и вертикальные параметры устанавливаются с помощью
DTHETA ( или DTHETAV) и DZ
Слайд 77

Радиальные модели Геометрия радиальной модели

Радиальные модели

Геометрия радиальной модели

Слайд 78

Радиальные модели По умолчанию используется Декартова геометрия модели. Для создания радиальной

Радиальные модели

По умолчанию используется Декартова геометрия модели. Для создания радиальной модели

используйте ключевое слово RADIAL в разделе RUNPSEC . Размеры модели определяются ключевым словом DIMENS они определяются числом ячеек в радиальном, тангенсальном и вертикальном направлениях (NR, NTHETA and NZ) вместо NX, NY и NZ.
Многие ключевые слова используемые для определения в радиальной геометрии и свойств ячеек сети отличаются от ключевых слов используемых в Декартовой версии. В общем X меняется на R, Y на THETA, а Z остается неизменной. В таблице приведены основные отличия между радиальной и Декартовой геометриями:
Слайд 79

Радиальные модели Центр радиальной модели обычно скважина (пустое место в центре

Радиальные модели

Центр радиальной модели обычно скважина (пустое место в центре диаграммы

с радиусом INRAD). Это необходимо для того, чтобы скважина не занимала каких-либо соседних ячеек и чтобы учитывался объем пор в пространстве вокруг ствола скважины (при несовпадении INRAD с радиусом скважины выводится предупреждение).
Радиальная геометрия модели может быть задана 3 способами; наиболее популярен третий:
Определение промежуточных радиусов по умолчанию
Определение промежуточных радиусов пользователем
Неполное определение радиусов пользователем
Слайд 80

Радиальные модели Определение промежуточных радиусов по умолчанию Внешний радиус может быть

Радиальные модели

Определение промежуточных радиусов по умолчанию
Внешний радиус может быть задан используя

OUTRAD. В таком случае количество ячеек между INRAD и OUTRAD равняется NR, определенному в RUNPSEC , а величина их радиусов соответствует логарифмическому распределению, основанному на зависимости:

где
Ri+1 - радиус i+1ой ячейки
Ri - радиус iой ячейки
RNR-1 - внутренний радиус крайнего блока
RU - последний радиус, определенный пользователем
NR - количество ячеек в радиальном направлении

Слайд 81

Радиальные модели Определение промежуточных радиусов пользователем Внешний радиус может быть задан

Радиальные модели

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

Средние радиусы ячеек могут быть установлены с помощью DR или DRV, также с использованием INRAD.
Неполное определение радиусов пользователем
По умолчанию средние радиусы рассчитывается ECLIPSE на основании логарифмической прогрессии, что может привести к слишком малым размерам ячейки в пространстве около ствола скважины, что в свою очередь может привести к трудностям при вычислении. По распространенной практике радиусам около скважины присваиваются определенные пользователем значения с помощью DR или DRV.
Слайд 82

Выходной контроль GRIDFILE По умолчанию ECLIPSE выводит файл .GRID содержащий геометрическое

Выходной контроль

GRIDFILE
По умолчанию ECLIPSE выводит файл .GRID содержащий геометрическое описание

только активных ячеек. Для включения геометрии всех ячеек используйте
GRIDFILE
2 /
Файл .GRID содержит только геометрию пласта в закодированной форме. По умолчанию он двоичен.
NOGGF
Ключевое слово NOGGF запрещает вывод файла геометрии сетки.
INIT
Требует вывод файла .INIT, содержащего статические свойства резервуара. К ним относятся:
• Все значения определенные в разделе GRID на каждую ячейку, такие как пористость, три проницаемости, коэффициент песчанистости, проводимости и коэффициенты проводимости.
• Функции насыщения (капиллярное давление и относительная проницаемость)
• Таблица PVT
• Области пласта, которые определены в REGION секции.
• Области течения
• Информация о NNC для построения в GRAF.
У INIT нет параметров. Файл .INIT по умолчанию не выводится. Этот файл кодированный и по умолчанию бинарный. Его используют такие пост-процессоры как GRAF.
Слайд 83

Выходной контроль RPTGRID ECLIPSE может быть настроена для вывода печатных отчетов

Выходной контроль

RPTGRID
ECLIPSE может быть настроена для вывода печатных отчетов о характеристиках

сетки в PRT файлы при помощи RPTGRID. Это ключевое слово зависит от определенного количества аргументов, в среднем около пятидесяти. В отличие от большинства других ключевых слов, аргументы должны быть в кавычках. К примеру
RPTGRID
'PORV' 'TRANX' 'PERMZ' 'ACTNUM' /
выдаст в PRT отчет об объеме пор ячеек, проводимостях в направлении X, проницаемости в направлении Z и номера активных ячеек.
BOUNDARY
Распечатка всего раздела GRID может быть весьма объемной. Если пользователя интересует определенный участок сетки, BOUNDARY используется для ограничения вывода RPTGRID определенной областью ячеек. К примеру,
BOUNDARY
--I1 I2 J1 J2 K1 K2
1 10 1 5 1 10 / ограничивает вывод выбранной областью.
Слайд 84

Выходной контроль

Выходной контроль

Слайд 85

Ключевые слова секции Grid Декартова геометрии Радиальная геометрии Свойства ячеек сетки

Ключевые слова секции Grid

Декартова геометрии
Радиальная геометрии
Свойства ячеек сетки для всех геометрий
Контроль

выклинивания и отключения ячеек
Задание проводимостей
Распределение проводимостей и модификаторов
Разломы
Численные водоносные пласты
Операторы (умножение, копирование и т.п.)
Двойная пористость / проницаемость
Опция граничных потоков
Термальная опция
Опция вертикального уравновешивания
Выходной контроль и другие
Слайд 86

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 87

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 88

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 89

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 90

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 91

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 92

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 93

Ключевые слова секции Grid

Ключевые слова секции Grid

Слайд 94

Лекция 6 СЕКЦИЯ PROPS Гидродинамическое моделирование

Лекция 6 СЕКЦИЯ PROPS

Гидродинамическое моделирование

Слайд 95

Цель секции PROPS • Этот раздел ограничивается заданием PVT свойств флюида

Цель секции PROPS

• Этот раздел ограничивается заданием PVT свойств флюида и

сжимаемости горных пород
• Для оценки PVT свойств всех флюидов в любой момент времени необходимо ввести следующие данные:
объемный коэффициент флюида в пласте,
вязкость,
соотношение нефть/газ (GOR) и/или соотношение газ/нефть (OGR)
• Необходимо задать условия перерасчета свойств флюидов из пластовых условий в поверхностные
• Также может быть явно задан начальный газовый фактор
• Если свойства внутри резервуара изменяются, то может быть указано более одной таблицы заданных свойств
Слайд 96

Цель секции PROPS

Цель секции PROPS

Слайд 97

Цель секции PROPS ECLIPSE представляет расчет материального баланса на каждом этапе

Цель секции PROPS

ECLIPSE представляет расчет материального баланса на каждом этапе моделирования.

Для этого необходимо рассчитать плотность каждой фазы. Т.к. плотность каждой фазы зависит от давления и количества дополнительных растворенных компонентов в каждой фазе, фазовые свойства PVT являются идеальным способом определения этих флюидных характеристик. В противном случае они выводятся из сочетания лабораторных экспериментов и полевых исследований. Этот тип данных требует определенных условий, т.е. поверхностной плотности каждого компонента.
ECLIPSE строит изотермические модели; температурные изменения не рассматриваются и совершенно очевидно, что все компоненты модели находятся при одной температуре.
Слайд 98

Фазовое поведение углеводородной системы Нефтяные и газовые флюиды в пласте представляют

Фазовое поведение углеводородной системы

Нефтяные и газовые флюиды в пласте представляют собой

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

Фазовое поведение углеводородной системы Если рассмотреть постепенное снижение давления (линия 1-2-3

Фазовое поведение углеводородной системы

Если рассмотреть постепенное снижение давления (линия 1-2-3 на

рисунке) при температуре между критической и крикодентермой. При снижении давления от точки 1 точка росы достигается и конденсируется жидкость, в результате в точке 2 система состоит из 5% жидкости и 95% газа, таким образом снижение давления приводит к появлению в системе жидкости .
Данные эффект называется ретроградной конденсацией.
При переходе от точки 2 к точке 3 количество жидкости в системе снижается и флюид вновь пересекает точку росы, становясь при этом 100% газообразным.
Слайд 100

Фазовое поведение углеводородной системы Распространенные УВ флюиды: Low-shrinkage oil (heavy oil

Фазовое поведение углеводородной системы

Распространенные УВ флюиды:
Low-shrinkage oil (heavy oil - black

oil) – слабосжимаемая нефть
High-shrinkage oil (volatile oil) – высокосжимаемая нефть
Retrograde condensate gas - газоконденсат
Wet gas – «жирный» газ
Dry Gas – «сухой» газ
Слайд 101

Фазовое поведение углеводородной системы

Фазовое поведение углеводородной системы

Слайд 102

Обзор Black oil • По определению, в поверхностных условиях «мертвая» нефть

Обзор Black oil

• По определению, в поверхностных условиях «мертвая» нефть не

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

Обзор Black oil Линии AA-HH представляют собой траектории по которым фазовое

Обзор Black oil

Линии AA-HH представляют собой траектории по которым фазовое состояние

нефть и/или газ может меняться в процессе разработки залежи.
Линия AA представляет нефть выше точки насыщения газом. Давление в
резервуаре изотермически снижается и вязкость и сжимаемость флюида могут
изменяться, при этом флюид остается 100% жидкостью. Концентрация растворенного газа не изменяется. Свободный газ, однако, добывается, т.к. газ выделяется из раствора в стволе скважины. В терминологии ECLIPSE это определяется как «мертвая» нефть.
Линия BB представляет нефть первоначально выше точки насыщения газом. При
уменьшении давления вязкость и сжимаемость нефти могут изменяться, оставаясь выше точки насыщения газом. Газовый фактор (GOR) постоянен. Составляющие остаются без изменений, пока не достигнута точка насыщения газом, затем газ начинает выделятся из раствора. Если это происходит в резервуаре, может формироваться газовая шапка. Если это происходит на забое скважины – свободный газ поступает в ствол скважины. В любом случае, т.к. давление снижается изотермически вдоль линии BB траектория флюидов пересекает качественную линию внутри фазовой огибающей и из раствора выделяется больше газа. GOR добытой нефти сравнительно ниже, чем нефти в залежи. В терминологии ECLIPSE это определяется как «живая» нефть.
Слайд 104

Обзор Black oil

Обзор Black oil

Слайд 105

Обзор Black oil Линия CC представляет двухфазную смесь. Выше газонефтяного контакта

Обзор Black oil

Линия CC представляет двухфазную смесь. Выше газонефтяного контакта (GOC)

давление ниже точки насыщения газом и газ существует в свободной фазе (газовая шапка). Ниже контакта газ-нефть давление выше точки насыщения; нефтяная фаза содержит растворенную газовую составляющую. Т.к. давление снижается изотермически вдоль линии CC смесь пересекает качественную линию внутри фазовой огибающей и из раствора выделяется больше газа В терминологии ECLIPSE это определяется как «живая» нефть.
Линия DD представляет первоначально околокритический флюид. Флюид, проходя через критическую точку (Pc, Tc) вдоль линии DD, превращается в двухфазную смесь, но это не будет явно заметно. При соответствующем сочетании пластового давления и давления на забое скважины (BHP) в скважину поступает однофазный флюид, который переходит в двухфазное состояние внутри скважинного ствола.
Линия EE представляет первоначально однофазный газ. Этот газ содержит
испаренные нефтяные компоненты, но при снижении давления вдоль линии EE
он не может достичь точки росы, так как находится вне фазовой огибающей газа
В терминологии ECLIPSE это определяется как сухой газ.
Слайд 106

Обзор Black oil

Обзор Black oil

Слайд 107

Обзор Black oil Линия FF это также двухфазная смесь, первоначально внутри

Обзор Black oil

Линия FF это также двухфазная смесь, первоначально внутри фазовой

огибающей. При снижении давления вдоль линии FF жидкая фаза постепенно испаряется, до пересечения точки росы, когда жидкости не остается. В терминологии ECLIPSE этот флюид определяется как жирный газ.
Линия GG представляет летучую нефть, первоначально выше точки росы. При повышении давления летучая фаза конденсируется. В терминологии ECLIPSE это также определяется как жирный газ.
Линия HH представляет нефть, которая изменяется при адиабатическом
переходе от однофазной нефти до двухфазной газонефтяной смеси. Это может
произойти в сепараторе.
Считается, что пластовые флюиды находятся в пласте, сепараторе и в поверхностных условиях при фиксированных температурах, хотя эти температуры могут быть различными. К примеру, температура при нормальных условиях обычно отличается от пластовой температуры. В таких случаях, данные PVT нефти могут быть представлены в таблицах свойств как функция давления, при определенной температуре.
Слайд 108

Обзор Black oil Однофазные системы, рассматривающие только газ, также хорошо приспособлены

Обзор Black oil

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

моделирования с использованием ECLIPSE. Как и в случае с черной нефтью, флюид находится далеко от критической точки и не пересекает линию точки росы таким образом концентрация летучей нефти (если она имеется) остается постоянной В терминологии ECLIPSE это определяется как сухой газ.
Программы моделирования темной нефти, такие как ECLIPSE не могут точно смоделировать переход компонентов из одной фазы в другую. Если необходимо смоделировать такой процесс, как высвобождение газа или выпадение конденсата, это должно происходить не напрямую, а путем изменения раствора GOR (Rs) и пара OGR (Rv). Флюиды, обладающими такими свойствами относятся к «живой» нефти и
жирному газу, соответственно.
Подобные процессы корректней моделировать при помощи Eclipse 300.
Слайд 109

Обзор Black oil Для того чтобы флюиды с использованием модели черной

Обзор Black oil

Для того чтобы флюиды с использованием модели черной нефти

могли быть смоделированы достаточно точно должны выполняться следующие условия:
• Количество выпавшего конденсата или высвободившегося газа должно составлять небольшую часть от общего количества углеводородов
• Оставшиеся углеводородные соединения не должны сильно меняться при высвобождении газа или выпадении конденсата
• Траектория флюидов должна быть удалена от критической точки
• Процесс должен быть изотермическим
Слайд 110

Модель черной нефти по сравнению с композиционной • Большую часть времени

Модель черной нефти по сравнению с композиционной

• Большую часть времени расчета

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

Модель черной нефти по сравнению с композиционной Моделирование темной нефти возможно

Модель черной нефти по сравнению с композиционной

Моделирование темной нефти возможно в

случаях, когда композиционные изменения незначительны и флюид удален от критической точки. В этом случае свойства флюида относительно стабильны и могут быть представлены в виде ряда таблиц свойств PVT как функции давления. Временные затраты при расчете PVT свойств с использованием таблиц очень малы. Основная часть машинного времени при моделировании черной нефти тратится на решение связанной системы уравнений потока для каждой ячейки, плюс ввод скважин и добыча флюидов.
Композиционное моделирование применяется в случаях, когда модель черной нефти работает не корректно. Программы композиционного моделирования, такие как ECLIPSE 300 представляют такие же расчеты потока флюидов, что и программы моделирующие темную нефть, но требуют больше компьютерных ресурсов.
Пластовый флюид описывается в терминах псевдокомпонентов. Каждый псевдокомпонент представляет группу отдельных углеводородных (и не только) соединений. После расчета потоков флюидов каждый псевдокомпонент должен быть приведен в равновесное состояние. Это процесс пошагового приближения. Затем, для каждого компонента решается кубическое уравнение состояния (EoS). Это также процесс пошагового приближения. Как и расчет равновесия, так и решение EoS выполняется для каждого момента времени во всех ячейках сетки.
На практике расчет потоков занимает меньше 50% машинного времени композиционного моделирования; расчет равновесия и решение EoS занимает остальное время.
Слайд 112

Уравнение состояния нефти • Все значения в таблицах зависят только от

Уравнение состояния нефти

• Все значения в таблицах зависят только от давления.
•

(s) относится к поверхностным условиям, (r) – к пластовым.
• В нефтяной фазе содержится растворенная газовая составляющая Vg(r)
• Уравнение не решается, оно заносится в таблицы, которые интерполируются и экстраполируются

Уравнение состояния для модели «черной» нефти

Слайд 113

Уравнение состояния нефти Уравнение состояния черной нефти (EoS) рассматривает нефть как

Уравнение состояния нефти

Уравнение состояния черной нефти (EoS) рассматривает нефть как однофазный

флюид, чьи свойства зависят только от давления. В терминах PVT черная нефть удалена от критической точки. Состав черной нефти представляется неизменным; другие свойства, такие как плотность и вязкость могут быть описаны достаточно легко, так как их поведение вдали от критической точки несложно. Черная нефть берется при фиксированной температуре, температуры пласта, скважины и поверхности также постоянны. В таких случаях нефтяные PVT данные могут быть представлены в таблицах свойств как функции только давления.
На практике модель черной нефти часто применяется ниже точки насыщения. «живая» нефть в терминологии ECLIPSE это нефть, которая может опускаться ниже точки насыщения, тогда как «мертвая» нефть не может этого делать. Если резервуар содержит газовую шапку, точка насыщения находится на газонефтяном контакте. Давление выше GOC, в газовой фазе, должно быть ниже точки насыщения. И «живая», и «мертвая» нефть содержат газовые составляющие, растворенные в нефтяной фазе. Это различие между фазами и компонентами является центральным в понимании модели черной нефти. Назначение элемента Rs в уравнении состояния черной нефти состоит в описании влияния составляющих растворенного газа на свойства нефтяной фазы. Для «мертвой» нефти он фиксирован; для «живой» нефти Rs может изменятся.
ВНИМАНИЕ «мертвая» нефть содержит газ. Он не выделяется из раствора в пластах, но растворенная газовая составляющая может быть высвобождена в любом месте, где давление ниже давления точки насыщения, например в стволе скважины или на поверхности.
Слайд 114

Уравнение состояния нефти Если пласт содержит только «мертвую» нефть с одним

Уравнение состояния нефти

Если пласт содержит только «мертвую» нефть с одним GOR

и точкой насыщения, единичное значение Rs может быть определено в точке насыщения с использованием RSCONST Если, с другой стороны, некоторое количество «мертвой» нефти полностью отделено (например блоковыми разломами или стратиграфическими ловушками), нефть с различными GOR и точкой насыщения может быть определена в различных частях пласта, используя ключевое слово PVTNUM вместе с некоторым количеством PVT таблиц для каждой нефти. Ключевое слово RSCONSTT должно быть использовано для определения точки насыщения GOR в каждой области PVT. Модель черной нефти может быть использована для моделирования «живой» нефти при некоторых условиях, встречающихся до определенного уровня, т.е PVT модель удовлетворяет следующим условиям:
• Нефтяная составляющая не изменяется, когда газ покидает раствор. В идеале, это значит, что распределение углеводородных молекул по весу одинаково для нефти и газа. Поскольку это невозможно, следует довольствоваться тем, что при испарении газа из раствора, свойства нефти изменяются слабо.
• Количество газа, испаряющегося из раствора, составляет небольшую долю от общего объема углеводородов. Так как выделяющийся газ изменяет состав остающейся нефти, это следует из предыдущего условия.
• Выделение газа не подводит смесь слишком близко к критической точке.
В моделях «живой» нефти начальный Rs определяется как GOC для начальной модели и как функция глубины, с использованием RSVD или PBVD.
Моделирование газа с использованием модели темной нефти происходит аналогично. Сухой газ соответствует «мертвой» нефти в терминологии ECLIPSE. Сухой газ может содержать компоненты нефтяных паров, описываемых соотношением нефть-газ Rv.
Слайд 115

Ввод данных PVT для ’мертвой’ нефти с помощью PVDO • Ключевое

Ввод данных PVT для ’мертвой’ нефти с помощью PVDO

• Ключевое слово

PVDO используется для определения свойств нефти выше точки насыщения (недонасыщенная нефть)
• Это таблицы зависимости коэффициента пластового объема и вязкости от давления
• Сводные таблицы PVDO могут повторять единичное PVDO если в пласте заведомо находится «мертвая» нефть
• Таблицы PVDO привязываются к определенным ячейкам или группам ячеек.
• Моделирование прерывается, если давление в любой ячейке сети опускается ниже точки насыщения.
Слайд 116

Ввод данных PVT для ’мертвой’ нефти с помощью PVDO Формат ключевого

Ввод данных PVT для ’мертвой’ нефти с помощью PVDO

Формат ключевого слова

включает параметры давления, коэффициента пластового объема нефти и вязкости, располагаемые слева направо. По умолчанию введенные таблицы интерполируются ECLIPSE, хотя давление не может быть задано по умолчанию. Таблица заканчивается одиночным прямым слэшем.
Значение давления монотонно увеличиваться вниз по таблице. Коэффициент пластового объема должен монотонно убывать, с увеличением давления. ECLIPSE не выдает предупреждения или сообщения об ошибке, если это не указано особо. Т.к. пластовое давление не может опускаться ниже точки насыщения, не имеет смысла включать данные ниже Pb и минимальным давлением в таблице должно быть давление в точке насыщения. Если давление в какой-нибудь ячейке достигает точки насыщения, моделирование прерывается. Нефтяная фаза содержит фиксированное количество растворенной газовой составляющей. Это определяется с использованием RSCONST. Если представлено более одной нефти, RSCONSTT используется для задания точек насыщения GOR для каждой нефти. Обратите внимание, что единицы измерения FIELD в RSCONST и RSCONSTT Mscf/stb, где М обозначает 1000.
Слайд 117

Ввод данных PVT для ’мертвой’ нефти с помощью PVCDO • Ключевое

Ввод данных PVT для ’мертвой’ нефти с помощью PVCDO

• Ключевое слово

PVСDO используется для определения свойств нефти выше точки насыщения (недонасыщенная нефть)
• Вместо определения зависимости Bo и μo от давления, эффективнее определять Bo и μo и наклон прямой в точке насыщения
• Максимальное давление определяется с использованием PMAX
• Несколько таблиц PVCDO могут быть указаны после ключевого слова PVCDO если в пласте находится несколько «мертвых» нефтей
• Таблицы PVCDO привязываются к определенным ячейкам или группам ячеек
• Моделирование прерывается, если давление в какой-либо ячейке сетки опускается ниже точки насыщения
Слайд 118

Ввод данных PVT для ’мертвой’ нефти с помощью PVCDO Формат ключевого

Ввод данных PVT для ’мертвой’ нефти с помощью PVCDO

Формат ключевого слова

включает параметры давления, коэффициента пластового объема нефти, вязкости и сжимаемости вязкости, располагаемые слева направо. Таблица заканчивается одиночным прямым слэшем. Сжимаемость определяется как:

или

Вместо определения кривых Bo и μo в зависимости от давления, определяются их значения и наклон в точке насыщения. Эта альтернатива PVDO применима, если недонасыщенные кривые являются прямыми линиями. В принципе, ECLIPSE может экстраполировать любое давление выше точки насыщения. Для предотвращения этого должно быть введено ключевое слово PMAX. Если давление в какой-нибудь ячейке достигает точки насыщения, моделирование прерывается.
Нефтяная фаза содержит фиксированное количество растворенной газовой составляющей. Это определяется с использованием RSCONST. Если представлено более одной нефти, RSCONSTT используется для задания точек насыщения GOR для каждой нефти. Обратите внимание, что единицы измерения FIELD в RSCONST и RSCONSTT Mscf/stb, где М обозначает 1000.

Слайд 119

Ввод данных PVT для ‘живой’ нефти с помощью PVTO • Ключевое

Ввод данных PVT для ‘живой’ нефти с помощью PVTO

• Ключевое слово

PVTO используется для определения свойств нефти выше (недонасыщенная нефть) и ниже (насыщенная нефть) точки насыщения
• Это таблицы зависимости давления, объемного коэффициента и вязкости от точки насыщения отношения газ/нефть
• Объемный коэффициент и вязкость для недонасыщенной нефти должны быть определены для максимального давления
• Несколько таблиц PVTO могут быть после ключевого слова PVTO если в пласте находится несколько нефтей.
• Таблицы PVTO привязываются к определенным ячейкам или группам ячеек
Слайд 120

Ввод данных PVT для ‘живой’ нефти с помощью PVTO

Ввод данных PVT для ‘живой’ нефти с помощью PVTO

Слайд 121

Ввод данных PVT для ‘живой’ нефти с помощью PVTO Формат ключевого

Ввод данных PVT для ‘живой’ нефти с помощью PVTO

Формат ключевого слова

для насыщенных областей включает параметры GOR при насыщенном давлении, давление насыщения, объемный коэффициент и вязкость, в порядке слева направо. Для недонасыщенных областей значения Rs опускаются, пока GOR фиксируется выше давления насыщения. По умолчанию, введенные данные интерполируются ECLIPSE, хотя значения Rs не могут быть заданы по определению. Таблица заканчивается одиночным прямым слэшем Пластовое давление может падать ниже точки насыщения и данные давления не обязательно будут выше наибольшего давления точки насыщения. Предположим давление в водо-нефтегазовом пласте на GOC составляет 3814.7, наибольшее давление точки насыщения в таблицах PVTO будет занимать предыдущую строку.
Требуется знать данные PVT нефти, при максимальном давлении подошвенной воды, т.к. ECLIPSE рассчитывает Rs в каждой ячейке. Все это также относится и к газовой шапке.
Слайд 122

Ввод данных PVT для ‘живой’ нефти с помощью PVTO При падении

Ввод данных PVT для ‘живой’ нефти с помощью PVTO

При падении давления

на GOC оно падает ниже давления насыщения и газ высвобождается. Rs нефти снижается по линии насыщения. Если в пласте поддерживается постоянное давление, в ячейке, скажем, 2614.7 содержащей свободный газ, будет происходить ре-поглощение высвободившегося газа, вдоль по кривым насыщения Rs от Pb. Если к тому моменту газ мигрирует выше и для репоглощения ничего не останется, ECLIPSE интерполирует кривую недонасыщения Rs при 0.465 между другими кривыми недонасыщения при Rs=0.241 и 0.77.
Слайд 123

Ввод данных PVT для ‘живой’ нефти с помощью PVCO • Ключевое

Ввод данных PVT для ‘живой’ нефти с помощью PVCO

• Ключевое PVCO

используется для определения свойств нефти выше (недонасыщенная) и ниже (насыщенная) точки насыщения
• Вместо определения зависимости Bo и μo от давления, определяется Bo и μo и наклон прямой в точке насыщения
• Несколько таблиц PVCO могут быть после ключевого слова PVCO если в пласте находится несколько нефтей.
• Таблицы PVCO привязываются к определенным ячейкам или группам ячеек.
Слайд 124

Ввод данных PVT для ‘живой’ нефти с помощью PVCO

Ввод данных PVT для ‘живой’ нефти с помощью PVCO

Слайд 125

Ввод данных PVT для ‘живой’ нефти с помощью PVCO Формат ключевого

Ввод данных PVT для ‘живой’ нефти с помощью PVCO

Формат ключевого слова

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

Уравнение состояния газа • Все значения зависят только от давления •

Уравнение состояния газа

• Все значения зависят только от давления
• (s) относится

к поверхностным условиям, (r) – к пластовым
• В газовой фразе содержится испаренная нефтяная составляющая Vg(r)
• Уравнение не решается, оно заносится в таблицы, которые интерполируются и экстраполируются
• Для сухого газа Rv фиксировано, а давление всегда выше точки росы
• Для жирного газа Rv получается при давлении ниже точки росы

Уравнение состояния газа для модели «черной» нефти

Слайд 127

Уравнение состояния газа Уравнение состояния черной нефти (EoS) для газа рассматривает

Уравнение состояния газа

Уравнение состояния черной нефти (EoS) для газа рассматривает газ

как одну фазу, чьи свойства зависят только от давления. В терминах PVT черная нефть удалена от критической точки. Состав газа представляется неизменным; другие свойства, такие как плотность и вязкость могут быть описаны достаточно легко, так как их поведение вдали от критической точки несложно. Газ берется при фиксированной температуре, температуры пласта, скважины и поверхности также постоянны. В таких случаях газовые PVT данные могут быть представлены в таблицах свойств как функции
только давления.
На практике модель черной нефти для газа часто применяется в двухфазных областях. Жирный газ в терминологии ECLIPSE это газ, который может пересекать кривую точки росы. Если резервуар содержит газовую шапку, точка росы находится на газонефтяном контакте. Давление выше GOC, в газовой фазе, должно быть ниже точки росы или выше, если точка росы ниже GOC. И сухой и жирный газ содержат испаренную нефтяную составляющую в газовой фазе. Это различие между фазами и компонентами является центральным в понимании модели черной нефти.
Назначение элемента Rs в уравнении состояния газа состоит в описании влияния составляющих испаренной нефти на свойства газовой фазы. Для сухого газа он фиксирован; для жирного газа Rv может изменяться.
Сухой газ может содержать нефть. Она не конденсируется в пласте или в стволе скважины, но испаренные компоненты конденсата высвобождаются в поверхностных условиях.
Слайд 128

Уравнение состояния газа Если пласт содержит только сухой газ с одним

Уравнение состояния газа

Если пласт содержит только сухой газ с одним OGR

и точкой росы, единичное значение Rv может быть определено в точке росы с использованием RVCONST.
Если, с другой стороны, некоторое количество сухого газа полностью отделено (например блоковыми разломами или стратиграфическими ловушками), газы с различными OGR и точками росы могут быть определены в различных частях пласта, используя ключевое слово PVTNUM вместе с некоторым количеством PVT таблиц для каждого газа. Ключевое слово RSCONSTT должно быть
использовано для определения точки росы OGR в каждой области PVT.
Модель черной нефти может быть использована для моделирования жирного газа при некоторых условиях:
• Состав газа не изменяется, когда нефть конденсируется из фазы
• Распределение углеводородных молекул по весу одинаково для нефти и газа
• Поскольку это невозможно, следует довольствоваться тем, что при конденсации нефти свойства газа изменяются слабо
• Количество нефти, выпадающей из газовой фазы, составляет небольшую долю от общего объема углеводородов. Т.к. конденсирующаяся нефть изменяет состав остающегося газа, это следует из предыдущего условия.
• Конденсация нефти не подводит смесь слишком близко к критической точке.
Слайд 129

Уравнение состояния газа Изменения составляющих не могут быть смоделированы напрямую, с

Уравнение состояния газа

Изменения составляющих не могут быть смоделированы напрямую, с использованием

ECLIPSE или любой другой модели темной нефти. Вместо этого, отношение Rs для газонефтяного пара используется для имитации эффекта изменения составляющих, когда нефть попадает или покидает газовую фазу. Для этого, Rs должен изменяться и в моделях жирного газа Rs заносится в таблицы как функция давления. В моделях жирного газа начальный Rs определяется на GOC при инициализации модели и как функция глубины с использованием ключевых слов RSVD или PBVD.
Слайд 130

Ввод PVT свойств сухого газа с помощью PVZG • Ключевое слово

Ввод PVT свойств сухого газа с помощью PVZG

• Ключевое слово PVZG

используется для определения свойств газа ниже точки росы или далеко от критической точки.
• Первая запись – это единичный пункт, определяющий относительную температуру для остальных данных.
• Следующие записи содержат данные о давлении, z-факторе и вязкости в приведенной последовательности.
• Несколько таблиц PVZG могут быть после ключевого слова PVZG если в пласте находится несколько газов.
• Таблицы PVZG привязываются к определенным ячейкам или группам ячеек.
• Моделирование прерывается, если давление в любой ячейке сети достигает точки росы.
Слайд 131

Ввод PVT свойств сухого газа с помощью PVZG Формат ключевого слова

Ввод PVT свойств сухого газа с помощью PVZG

Формат ключевого слова включает

относительное давление в первой записи с последующими параметрами давления, газового z-фактора и вязкости, располагаемыми слева направо. По умолчанию введенные таблицы интерполируются ECLIPSE, хотя давление не может быть задано по умолчанию.
Таблица заканчивается одиночным прямым слэшем.
Значение давления монотонно увеличиваться вниз по таблице. Z-фактор должен монотонно увеличиваться, с увеличением давления. ECLIPSE не выдает предупреждения или сообщения об ошибке, если это не указано особо. Т.к. пластовое давление не может опускаться ниже точки росы, не имеет смысла включать данные ниже Pd и минимальным давлением в таблице должно быть давление в точке росы. Если давление в какой-нибудь ячейке достигает точки росы, моделирование прерывается.
Газовая фаза содержит фиксированное количество испаренных компонентов. Это определяется с использованием RVCONST. Если представлено более одного газа, RVCONSTT используется для задания точки росы GOR для каждого газа. Обратите внимание, что единицы измерения FIELD в RSCONST и RSCONSTT Mscf/stb, где М обозначает 1000.
Слайд 132

Ввод PVT свойств сухого газа с помощью PVDG • Ключевое слово

Ввод PVT свойств сухого газа с помощью PVDG

• Ключевое слово PVDG

используется для определения свойств газа ниже точки росы или значительно ниже критической точки
• Это таблица зависимости объемного коэффициента и вязкости от давления
• Несколько таблиц PVDG могут быть после ключевого слова PVDG если в пласте находится несколько газов.
• Таблицы PVDG привязываются к определенным ячейкам или группам ячеек.
• Моделирование прерывается, если давление в любой ячейке сети достигает точки росы.
• PVDG является альтернативой PVZG
Слайд 133

Ввод PVT свойств сухого газа с помощью PVDG Формат ключевого слова

Ввод PVT свойств сухого газа с помощью PVDG

Формат ключевого слова включает

параметры давления, коэффициента пластового объема газа и вязкости, располагаемые слева направо. По умолчанию введенные таблицы интерполируются ECLIPSE, хотя давление не может быть задано по умолчанию. Таблица заканчивается одиночным прямым слэшем.
Значение давления монотонно увеличиваться вниз по таблице. Коэффициент пластового объема должен монотонно увеличиваться, с увеличением давления. ECLIPSE не выдает предупреждения или сообщения об ошибке, если это не
указано особо.
Т.к. пластовое давление не может опускаться ниже точки росы, не имеет смысла включать данные ниже Pd и минимальным давлением в таблице должно быть давление в точке росы. Если давление в какой-нибудь ячейке достигает точки росы, моделирование прерывается. Газовая фаза содержит фиксированное количество испаренных компонентов. Это определяется с использованием RVCONST. Если представлено более одного газа, RVCONSTT используется для задания GOR в точке росы для каждого газа. Обратите внимание, что единицы измерения FIELD в RSCONST и RSCONSTT Mscf/stb, где М обозначает 1000.
Слайд 134

Ввод PVT данных жирного газа с использованием PVTG • Ключевое слово

Ввод PVT данных жирного газа с использованием PVTG

• Ключевое слово PVTG

используется для определения свойств газа выше (недонасыщенный) и ниже (насыщенный) точки росы
• Это таблицы зависимости Rv, Bg и вязкости от давления
• Недонасыщенные Rv, Bg и μg должны быть определены для максимального Pg.
• Несколько таблиц PVTG могут быть после ключевого слова PVTG если в пласте находится несколько жирных газов.
• Таблицы PVTG привязываются к определенным ячейкам или группам ячеек.
Слайд 135

Ввод PVT данных жирного газа с использованием PVTG Формат ключевого слова

Ввод PVT данных жирного газа с использованием PVTG

Формат ключевого слова для

насыщенных областей включает параметры давления, отношения нефть/газ, коэффициент сжимаемости газа и вязкость, в порядке слева направо. Для недонасыщенных областей значения Pg допускаются в дополнительных строчках данных. Введенные по умолчанию данные интерполируются ECLIPSE, хотя значения Pg не могут быть заданы по умолчанию. Таблица заканчивается одиночным прямым слэшем.
Зависимость коэффициента пластвого объема и вязкости от Rv. Прямые линии, пересекающие левую вертикальную ось представляют пересечение фазовой огибающей при фиксированном давлении.
Слайд 136

Уравнение состояния воды • Все значения зависят только от давления •

Уравнение состояния воды

• Все значения зависят только от давления
• (s) относится

к поверхностным условиям, (r) – к пластовым
• Водная фаза содержит только воду.
• Нефть и газ не могут растворятся в воде и наоборот
• Для каждой области PVT необходима одна строка данных в PVTW

Уравнение состояния воды

Слайд 137

Уравнение состояния воды Уравнение состояния модели черной нефти для воды (EoS)

Уравнение состояния воды

Уравнение состояния модели черной нефти для воды (EoS) рассматривает

ее как однофазный флюид, чьи свойства зависят только от давления. PVT свойства воды задаются с использованием ключевого слова PVTW. К примеру
PVTW
--Pref Bw Cw μw Cwμ
4000 1.03 3.0E-6 040 0.0 /
где
Pref Давление для которого заданы все параметры
Bw Объемный коэффициент воды при относительном давлении
Cw Сжимаемость воды

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

Слайд 138

Относительные плотности • В пласте жидкая углеводородная фаза это нефть обычно

Относительные плотности

• В пласте жидкая углеводородная фаза это нефть обычно с

растворенным в ней газом
• В пласте газообразная углеводородная фаза это газ, возможно с растворенной в нем нефтью
• Характеристики нефтяного и газового компонентов замеряются на линии сепараторов.
• Плотности нефтяной и газовой составляющей задаются с использованием DENSITY или GRAVITY.
Слайд 139

Относительные плотности В модели черной нефти все PVT свойства являются функциями

Относительные плотности

В модели черной нефти все PVT свойства являются функциями давления.
Поверхностные

свойства каждого компонента также зависят от давления и температуры. Стандартные давление и температура не определяются в ECLIPSE; они могут принимать любые значения на поверхности. Плотности в стандартных условиях могут быть заданы через плотности или через удельный вес нефти в единицах API. К примеру:
GRAVITY
-- Плотность нефти Отн. плотность Отн. плотность
-- API воды газа
32 1.050 0.700 /
или
DENSITY
--Плотность Плотность Плотность
-- Нефти Воды Газа
-- Kg/m3 Kg/m3 Kg/m3
865 1050 0.9051 /
Где gravity – удельный вес жидкости, а API измеряется в градусах. Удельный вес воды это отношение плотности воды к плотности чистой воды. Удельный вес газа – это отношение плотности газа к плотности воздуха. Этот и другие методы пересчета единиц измерения можно найти в ECLIPSE 100 TECHNICAL DESCRIPTION.
Слайд 140

Фазы модели черной нефти • Могут быть одна, две или три

Фазы модели черной нефти

• Могут быть одна, две или три активных

фазы: OIL, WATER и GAS.
• Фаза OIL может содержать растворенную газовую компоненту, определенную DISGAS
• Фаза GAS может содержать испаренную нефтяную компоненту, определенную VAPOIL
• Одиночная фаза не может содержать компонентов, которые образуют другие фазы в процессе моделирования.
Слайд 141

Фазы модели черной нефти

Фазы модели черной нефти

Слайд 142

Фазы модели черной нефти Трехфазное моделирование В трехфазном моделировании OIL, WATER

Фазы модели черной нефти

Трехфазное моделирование
В трехфазном моделировании OIL, WATER и GAS

представлены и определены в разделе RUNSPEC. В этом случае должны присутствовать «живая» нефть или жирный газ. В результате газ выделяется из раствора, когда нефть пересекает точку насыщения, или нефть конденсируется из паров, когда газ пересекает точку росы, в зависимости от того, находится ли фазовая смесь слева или справа от критической точки.
• В трехфазной модели слева от критической точки система должна насыщаться нефтью ниже критической температуры и давления. Газ сухой, но нефть содержит газовую составляющую, которая может выделяться из раствора. Ключевые слова раздела RUNPSEC - OIL, WATER, GAS и DISGAS. Также могут быть некоторые изменения в отношении газ/нефть в зависимости от глубины, которые первоначально устанавливаются с использованием ключевого слова RSVD в разделе SOLUTION. Или же изменения точки насыщения в зависимости от глубины, которые устанавливаются с использованием PBVD в разделе SOLUTION. Значение Rs может изменяться в ходе процесса.
Слайд 143

Фазы модели черной нефти Трехфазное моделирование • В трехфазной модели справа

Фазы модели черной нефти

Трехфазное моделирование
• В трехфазной модели справа от критической

точки система должна быть
газоконденсатной. Нефть «мертвая», но газ содержит испаренную нефтяную
компоненту, которая может конденсироваться. Ключевые слова раздела
RUNPSEC - OIL, WATER, GAS и VAPOIL. Также могут быть некоторые изменения в
отношении нефть/газ в зависимости от глубины, которые первоначально
устанавливаются с использованием ключевого слова RVVD в разделе SOLUTION,
или изменения точки насыщения в зависимости от глубины, которые
устанавливаются с использованием PDVD в разделе SOLUTION. Значение Rv
может изменяться в ходе процесса.
Слайд 144

Фазы модели черной нефти Двухфазное моделирование В двухфазном моделировании возможны сочетания

Фазы модели черной нефти

Двухфазное моделирование
В двухфазном моделировании возможны сочетания OIL и

WATER, WATER и GAS,
GAS и OIL.
• В системе WATER и OIL газовая фаза не определена, таким образом свободный
газ отсутствует. Нефть всегда выше критической точки. Постоянное значение Rs нефти может быть задано в каждой области PVT с использованием RSCONST или RSCONSTT. RSVD можно не использовать.
• В системе GAS и WATER нефтяная фаза не определена, таким образом жидкая нефть отсутствует. Газ не пересекает точку росы. Постоянное значение Rv газа может быть задано для каждой области PVT с использованием RVCONST или RVCONSTT. RVVD можно не использовать.
• В системе OIL и GAS допускается наличие газа, растворенного в нефти и паров
нефти в газе. Могут, но не обязательно, использоваться ключевые слова VAPOIL и DISGAS вместе или по отдельности. Ключевые слова RSVD и/или RVVD могут применяться при необходимости.
Слайд 145

Фазы модели черной нефти Однофазное моделирование В модели используется одна фаза.

Фазы модели черной нефти

Однофазное моделирование
В модели используется одна фаза. Это может

быть OIL, WATER или GAS.
• При моделировании нефти газовая фаза не определена, таким образом нефть
всегда выше критической точки. Постоянное значение Rs нефти может быть
задано в каждой области PVT с использованием RSCONST или RSCONSTT. RSVD
можно не использовать
• При моделировании газа нефтяная фаза не определена, таким образом газ не
пересекает точку росы. Постоянное значение Rv газа может быть задано для
каждой области PVT с использованием RVCONST или RVCONSTT. RVVD можно не
использовать.
• При моделировании воды представлена только водяная фаза.
Слайд 146

Описание различных PVT свойств с помощью PVT регионов • По умолчанию

Описание различных PVT свойств
с помощью PVT регионов

• По умолчанию моделирование

разных типов нефтей в отдельных областях состоит в задании всех PVT свойств для каждой отдельной области.
• Каждой ячейке в каждой области назначается значение с использованием
PVTNUM в разделе REGIONS
• PVTNUM это целое, которое указывает, какие таблицы PVT должны быть
использованы в данной ячейке
• Сводные таблицы PVT вводятся под каждым соответствующим ключевым
словом
• Таблицы нумеруются в порядке ввода
• Одинаковые таблицы могут быть устранены. Это может быть, к примеру, если
PVT нефти изменяется, а свойства остальных флюидов - нет.
• Флюид, находящийся в ячейке, обладает PVT свойствами этой ячейки
Слайд 147

Описание различных PVT свойств с помощью PVT регионов

Описание различных PVT свойств
с помощью PVT регионов

Слайд 148

Описание различных PVT свойств с помощью PVT регионов Если пласт может

Описание различных PVT свойств
с помощью PVT регионов

Если пласт может быть

разделен на несколько областей, PVT свойства которых
значительно отличаются, простейший метод описания свойств состоит в
составлении таблиц PVT для каждой области. Место для таблиц должно быть
также зарезервировано в разделе RUNSPEC.
Как определить разные PVT свойства с помощью PVT регионов
В качестве примера рассмотрим модель «мертвой» нефти и воды 20*5*10,
включающую нефть с различными точками насыщения и поверхностной
плотностью на участках (1-10, 1-5, 1-10) и (11-20, 1-5, 1-10). Соответствующими
ключевыми словами PVT могут быть
PVTW
4000 1.03 3.0E-6 0.40 0.00 / PVT таблица 1 для воды
--PVT таблица 2 для воды аналогична таблице 1
/
Слайд 149

Описание различных PVT свойств с помощью PVT регионов PVDO 2500 1.260

Описание различных PVT свойств
с помощью PVT регионов

PVDO
2500 1.260 0.5
3000 1.257

0.5
3500 1.254 0.5
4000 1.251 0.5
4500 1.248 0.5
/ PVT таблица 1 для «мертвой» нефти
2550 1.191 0.5
3050 1.198 0.5
3550 1.205 0.5
4050 1.213 0.5
4550 1.220 0.5
/ PVT таблица 2 для «мертвой» нефти
GRAVITY
32 1.050 0.700 / Таблица 1 плотности флюида в поверхн. условиях.
33.5 1.050 0.700 / Таблица 2 плотности флюида в поверхн. условиях.
RSCONSTT
0.656 2500 / Точка насыщения GOR для нефти в таблице 1
0.670 2550 / Точка насыщения GOR для нефти в таблице 2
Любые таблицы PVT для газа должны выглядеть таким же образом.
Слайд 150

Описание различных PVT свойств с помощью PVT регионов Дополнительные ключевые слова

Описание различных PVT свойств
с помощью PVT регионов

Дополнительные ключевые слова раздела

RUNSPEC могут быть следующими
TABDIMS
--NTSFUN NTPVT
1* 2 /
где NTSFUN количество вводимых таблиц функций насыщения, а NTPVT количество введенных таблиц PVT.
Каждая область PVT должна быть определена с использованием PVTNUM в разделе REGIONS. Здесь соответствующими ключевыми словами будут REGIONS
EQUALS
'PVTNUM' 1 /
'EQLNUM' 1 /
'PVTNUM' 2 11 20 /
'EQLNUM' 2 11 20 /
/
Области EQLNUM соответствующие областям PVTNUM должны быть определены, т.к. ECLIPSE требует, чтобы каждая область PVT была сбалансирована отдельно. Существует дополнительное требование по определению количества сбалансированных областей в разделе RUNSPEC с использованием ключевого слова EQLNUM.
Слайд 151

Описание различных PVT свойств с помощью PVT регионов Хотя этот метод

Описание различных PVT свойств
с помощью PVT регионов

Хотя этот метод достаточно

прост в применении, существуют следующие неудобства:
• Бессмысленно предполагать, что PVT свойства флюида в реальности зависят от того, в какой ячейке он оказался.
• Могут возникнуть проблемы несходимости. К примеру, рассмотрим нефть с GOR 0.5, попадающую в ячейку с максимальным GOR 0.4. Для сохранения массы должен быть высвобожден газ. Если нет соответствующих данных в PVT таблице газа в ключевом слове PVDG, ECLIPSE будет экстраполировать газовые таблицы PVT по необходимости. При этом экстраполяция может быть некорректной из-за неполной информации.
Слайд 152

Ключевые слова для API трассировки и начальное API нефти • С

Ключевые слова для API трассировки
и начальное API нефти

• С помощью

API трассировки начальная плотность API устанавливается в
модели. PVT свойства воды и газа остаются неизменными.
• ECLIPSE отслеживает плотность нефти в каждой ячейке
• Поверхностная плотность нефти берется из таблиц PVT
• Соотношение в ячейке нефти с различной плотностью может быть рассчитано из потока
• Плотность на поверхности смеси это плотность в поверхностных условиях нефтяных компонент, взвешенная в соответствии с их концентрацией
• Таблица PVT смешанной нефти интерполируется в каждой ячейке из нефтяных компонент, пропорционально их концентрации.
• Плотность смеси в условиях пласта рассчитывается из таблиц PVT
• Для активизации этой опции используется ключевое слово API в разделе RUNSPEC
• Начальная плотность API нефти определяется с помощью OILAPI или APIVD. Не разделяйте пласт на разные PVT и уравновешивающие регионы
Слайд 153

Ключевые слова для API трассировки и начальное API нефти

Ключевые слова для API трассировки
и начальное API нефти

Слайд 154

Ключевые слова для API трассировки и начальное API нефти Описание различных

Ключевые слова для API трассировки
и начальное API нефти

Описание различных PVT

свойств с помощью API трассировки
Опция API трассировки используется для расчета PVT свойств смеси различных нефтей.
Как выполнить API трассировку
• Ключевое слово API должно присутствовать в разделе RUNSPEC.
• Эта опция не требует деления пласта на PVT и уравновешивающие регионы.
• В разделе SOLUTION начальная API плотность нефти задается напрямую в ячейке с использованием OILAPI или как функция глубины, с использованием APIVD.
• Должны быть заданы хотя бы два комплекта таблиц PVT. Они должны охватывать все значения первоначальной API плотности нефти.
• Память для таблиц PVT резервируется с использованием TABDIMS в разделе RUNPSEC.
• Так как API отслеживание рассчитывает только PVT свойства нефти, дублирующие таблицы PVT для газа и воды могут быть опущены.
Слайд 155

Ключевые слова для API трассировки и начальное API нефти В процессе

Ключевые слова для API трассировки
и начальное API нефти

В процессе моделирования:
•

ECLIPSE рассчитывает количество различной API нефти, перетекающей из ячейки в ячейку.
• Рассчитывается плотность в поверхностных условиях каждой нефти
• Средняя плотность линейно интерполируется на основе пропорций каждой представленной нефти.
• Таким же образом рассчитываются средние PVT.
• PVT затем подгоняются под условия пласта.
К примеру, половину ячейки с нефтью 30° API занимает нефть с 20° API. Это поверхностные плотности. ECLIPSE рассчитывает поверхностную плотность нефтяной смеси как 25° API. Это используется в слагаемом ρo(s) уравнения состояния. Новые PVT свойства нефти рассчитываются как среднее арифметическое от таблиц PVT для нефти с 20° и 30°.
Слайд 156

Ключевые слова для API трассировки и начальное API нефти Этот метод

Ключевые слова для API трассировки
и начальное API нефти

Этот метод изменения

свойств PVT нефтяной смеси имеет следующие преимущества:
• Могут приниматься во внимание усредненные нефтяные смеси
• PVT свойства нефти не изменяются при перетекании из одной ячейки в другую
Принципиальный недостаток в следующем -
• ECLIPSE линейно интерполирует таблицы PVT пропорционально массе каждой нефти, но зависимость свойств нефти от API нелинейна. Количество таблиц PVT должно быть дополнено PVT таблицами для промежуточных плотностей нефти, для уверенности в адекватности результата интерполяции. Чем больше это количество, тем, скорее всего, будет лучше результат
Слайд 157

Ключевые слова для API трассировки и начальное API нефти Пример данных

Ключевые слова для API трассировки
и начальное API нефти

Пример данных API

трассировки
В качестве примера рассмотрим нефть, со значениями от 40° на 10800 футах до 24° на 11920 футах. PVT свойства газа и воды в пласте не изменяются. API плотность измерена на восьми глубинах. Раздел PROPS для свойств флюида может содержать следующее:
PVDG
-- PVT таблица No 1 для сухого газа.
15 235 .008
500 6.711 .0109
1000 3.145 .0139
1500 1.987 .0168
2000 1.431 .0197
2500 1.117 .0226
3000 .923 .0256
3500 .798 .0285
4000 .716 .0314
4500 .663 .0344
5000 .630 .0373
5500 .612 .0402
5600 .609 .0408
6000 .599 .0431
6500 .596 .0460
/
--Все остальные PVT таблицы сходны с таблицей No.1,
--потому что API трассировка влияет на свойства нефти. /
Слайд 158

Ключевые слова для API трассировки и начальное API нефти Пример данных

Ключевые слова для API трассировки
и начальное API нефти

Пример данных API

трассировки
PVTO
-- PVT таблица No.1 для «живой» нефти
.200 600 1.185 .60 /
.400 1400 1.285 .60 /
.622 2200 1.395 .60 /
.686 2420 1.430 .60 /
.746 2650 1.460 .60 /
.804 2830 1.485 .60 /
.900 3130 1.530 .60 /
1.050 3590 1.600 .60 /
1.358 4430 1.750 .60
6500 1.738 .60 /
1.882 5600 2.020 .60
6500 2.000 .60 /
2.285 6500 2.230 .60
7000 2.218 .60 /
/ PVT таблица No.2 для «живой» нефти
.200 700 1.190 .60 /
.400 1500 1.295 .60 /
.622 2350 1.400 .60 /
.686 2570 1.430 .60 /
.746 2800 1.465 .60 /
.804 2990 1.490 .60 /
.900 3290 1.535 .60 /
1.050 3750 1.605 .60 /
1.358 4600 1.745 .60
6500 1.700 .60 /
1.882 5970 1.970 .60 /
2.285 7020 2.140 .60
7500 2.126 .60 /
Слайд 159

Ключевые слова для API трассировки и начальное API нефти Пример данных

Ключевые слова для API трассировки
и начальное API нефти

Пример данных API

трассировки
/ PVT таблица No.3 для «живой» нефти
.200 770 1.190 .60 /
.400 1610 1.300 .60 /
.622 2460 1.405 .60 /
.686 2700 1.435 .60 /
.746 2920 1.470 .60 /
.804 3110 1.495 .60 /
.900 3430 1.535 .60 /
1.050 3900 1.595 .60
6500 1.530 .60 /
1.358 4850 1.720 .60 /
1.882 6470 1.930 .60 /
2.285 7720 2.090 .60
8000 2.081 .60 /
/ PVT таблица No.4 для живой нефти
.200 830 1.195 .60 /
.400 1700 1.305 .60 /
.622 2600 1.410 .60 /
.686 2840 1.440 .60 /
.746 3060 1.470 .60 /
.804 3270 1.495 .60 /
.900 3600 1.535 .60
6500 1.470 .60 /
1.050 4120 1.595 .60 /
1.358 5130 1.715 .60 /
1.882 6850 1.920 .60 /
2.285 8170 2.070 .60
8500 2.060 .60 /
/ PVT таблицы 5, 6, 7 и 8 аналогичны таблице 4
/
/
/
/
Слайд 160

Ключевые слова для API трассировки и начальное API нефти PVTW --

Ключевые слова для API трассировки
и начальное API нефти

PVTW
-- PVT таблица

No 1 для воды
5000 1.0000 3.0E-6 0.3 0 /
-- PVT таблицы 2-8 аналогичны таблице No. 1
-- потому что API трассировка влияет только на нефтяную фазу/
/
/
/
/
/
/
ROCK
-- Таблица сжимаемости породы № 1
50003.5E-06 /
-- Таблицы сжимаемости 2-8 аналогичны таблице 1
/
/
/
/
/
/
/
Слайд 161

Ключевые слова для API трассировки и начальное API нефти GRAVITY --

Ключевые слова для API трассировки
и начальное API нефти
GRAVITY
-- 8 таблиц

поверхностной API плотности
40.0 1.16 0.824 /
38.7 1.16 0.824 /
36.4 1.16 0.824 /
34.5 1.16 0.824 /
31.9 1.16 0.824 /
29.7 1.16 0.824 /
26.8 1.16 0.824 /
24.0 1.16 0.824 /
SOLUTION
APIVD
--8 таблиц зависимости API плотности от глубины
10800 40.0
11100 38.7
11350 36.4
11500 34.5
11670 31.9
11740 29.7
11810 26.8
11920 24.0 /
Слайд 162

Сжимаемость горной породы • Следует определять сжимаемость горной породы, т.к. объем

Сжимаемость горной породы

• Следует определять сжимаемость горной породы, т.к. объем пор

изменяется при изменении давления
• Если сжимаемость породы обратима и одинакова используйте ключевое слово ROCK. Данные RUNPSEC в этом случае не нужны.
• В противном случае используйте опцию сжатия породы
• Установите опции сжимаемости и гистерезиса с использованием ROCKCOMP в RUNSPEC
• Используйте ROCKTAB или ROCKTABH для определения данных о сжатии горной породы, включая модификаторы объема пор и проводимости в зависимости от давления
• ROCKTAB используется для определения обратимой сжимаемости
• ROCKTABH используется для определения необратимой сжимаемости, т.е. гистерезиса сжатия.
• Используйте ROCKNUM для обозначения областей различных типов пород вместе с различными таблицами ROCKTAB или ROCKTABH
• Ключевое слово OVERBURD может быть использовано в любом случае для определения таблиц перегрузки пород
Слайд 163

Сжимаемость горной породы

Сжимаемость горной породы

Слайд 164

Сжимаемость горной породы Где индексы R и S относятся, соответственно, к

Сжимаемость горной породы

Где индексы R и S относятся, соответственно, к пластовому

и поверхностному давлению, а C – сжимаемость породы. ECLIPSE использует приблизительную форму.
ROCK используется для задания однородной сжимаемости в каждой области PVT, т.е. количество таблиц в ключевом слове ROCK должно быть равно количеству используемых таблиц PVT которое задается параметром NTPVT в RUNSPEC
ключевым словом TABDIMS.
Если установлен ROCKOPTS переключатель 3=’SATNUM’, таблицы сжимаемости пород связываются с регионами насыщенности (SATNUM), а не с PVT регионами (PVTNUM).
ROCKTAB и ROCKTABH используются для установки обратимой и гистерезисной сжимаемости породы, соответственно. Введенное ключевое слово следует за таблицей NTROCC как установлено в RUNSPEC ключевым словом ROCKCOMP. Каждая таблица содержит колонки давления, множителя объема пор и соответствующего множителя проводимости.
Слайд 165

Функции насыщенности и масштабирование концевых точек • Минимальные требуемые моделью данные

Функции насыщенности и
масштабирование концевых точек

• Минимальные требуемые моделью данные –

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

Функции насыщенности и масштабирование концевых точек Функция насыщенности предназначена для: •

Функции насыщенности и масштабирование концевых точек

Функция насыщенности предназначена для:
• Определения верхних

и нижних границ насыщенности каждой активной фазы. Что в свою очередь используется для определения начальной насыщенности каждой фазы в газовой, нефтяной и водной зонах.
• Для определения капиллярных давлений так, чтобы можно было рассчитать первоначальную переходную зону для каждой фазы.
• Представить данные относительной проницаемости для вычисления мобильности флюидов и решения уравнений потока между ячейками, из ячейки в скважину и наоборот.
Группа 1
SWOF Krw, Krow, Pcow vs. Sw @ (связанный газ – не активный)
SGOF Krg, Kgog, Pcog vs. Sg @ связанная вода
SLGOF Krg, Kgog, Pcog vs. Sl @ связанная вода
Группа 2
SWFN Krw, Pcw vs Sw (в 3-фазных или водонефтяных системах при связанном газе)
SGFN Krg, Pcg vs. Sg (в 3-фазных или газонефтяных системах при связанной воде)
SOF3 Kro, Krg vs So
SOF2 Kro vs So
SOF32D Kro vs Sw and Sg
Ключевые слова из разных групп нельзя смешивать
Слайд 167

Функции насыщенности и масштабирование концевых точек Функции насыщенности это таблицы относительной

Функции насыщенности и масштабирование концевых точек

Функции насыщенности это таблицы относительной проницаемости

и капиллярного давления в зависимости от насыщенности. Они должны быть представлены для всех фаз.
Правила создания и оформления таблиц:
• Каждая таблица содержит множество колонок данных
• Во всех колонках любой таблицы должно быть одинаковое количество значений.
• Количество рядов данных в каждой таблице должно быть не меньше двух и не больше параметра NSSFUN ключевого слова TABDIMS в разделе RUNPSEC.
• При использовании нескольких таблиц их количество должно быть определено при помощи NTSFUN в разделе RUNPSEC.
• Фазовая насыщенность и относительная проницаемость замещающей фазы должны находиться в пределах между 0 и 1 и монотонно возрастать вниз по колонке, за исключением ключевых слов SOF2, SOF3 и SOF32D.
Слайд 168

Функции насыщенности и масштабирование концевых точек Правила создания и оформления таблиц:

Функции насыщенности и масштабирование концевых точек
Правила создания и оформления таблиц:
• Относительная

проницаемость замещаемой фазы должна находиться в пределах между 0 и 1 и монотонно убывать вниз по колонке, за исключением ключевых слов SOF2, SOF3 и SOF32D.
• Капиллярное давление должно быть больше или равно нулю, оставаться неизменным или убывать вниз по таблице.
• Значения насыщенности не могут быть заданы по умолчанию
• Заданные по умолчанию значения относительной фазовой проницаемости будут интерполированы ECLIPSE.
• Каждая таблица заканчивается прямым слэшем (/) и каждое ключевое слово может содержать несколько таблиц.
• Все последующие таблицы может быть определены по умолчанию, с
использованием прямого слэша, после первой.
• Всегда So+Sw+Sg=1.
Слайд 169

Функции насыщенности и масштабирование концевых точек Первый набор ключевых слов функций

Функции насыщенности и масштабирование концевых точек

Первый набор ключевых слов функций насыщенности


SWOF
Это ключевое слово может использоваться в моделях, содержащих как воду, так и нефть. Если газ также находится в активной фазе, функции насыщения нефть/газ могут быть введены с использованием SGOF или SLGOF.
SWOF
--Sw krw krow Pcow
0.2 0.00 0.90 4.0
0.3 0.05 0.75 2.0
0.4 1* 0.55 1.0
0.5 2* 0.5
0.6 3*
0.7 0.40 0.00 0.10 /
т.е. колонки водонасыщенность, относительная проницаемость воды, относительная проницаемость нефти в воде, когда в системе представлены только две фазы нефть и вода, а также водо-нефтяное капиллярное давление. Относительные проницаемости и капиллярное давление – функции от водонасыщенности.
Данная таблица интерпретируется как относительная проницаемость нефть/вода при начальной насыщенности газа (начальная насыщенность газа обычно равна нулю). Если присутствует свободный газ, то относительная проницаемость нефть/вода при максимальной насыщенности нефти должна быть такой же, как относительная проницаемость нефть/газ при минимальном газонасыщении.
Слайд 170

Функции насыщенности и масштабирование концевых точек Первый набор ключевых слов функций

Функции насыщенности и масштабирование концевых точек

Первый набор ключевых слов функций насыщенности


SGOF
Это ключевое слово может использоваться в моделях, содержащих как газ, так и нефть как активные фазы . Если вода тоже активна, функции насыщения вода/нефть должны быть введены с использованием SWOF.
SGOF
--Sg krg krog Pcog
т.е. колонки газонасыщенности, относительной проницаемости газа и относительной проницаемости нефти в газе при водонасыщенности на уровне связанной воды (connate water saturation). Относительные проницаемости и капиллярное давление – функции от газонасыщенности.
Эта таблица интерпретируется, как двухфазная относительная проницаемость нефть/газ при наличии в системе только связанной воды. При наличии свободной воды, относительная проницаемость нефть/газ при начальной насыщенности газа должна равняться относительной проницаемости нефть/вода при максимальной нефтенасыщенности.
Слайд 171

Функции насыщенности и масштабирование концевых точек Первый набор ключевых слов функций

Функции насыщенности и масштабирование концевых точек

Первый набор ключевых слов функций насыщенности


SLGOF
Это ключевое слово может использоваться в моделях, содержащих газ и нефть в активных фазах. Если вода тоже активна, функции насыщения вода/нефть должны быть введены с использованием SWOF
SLGOF
--Sl krg krog Pcog
т.е. колонки насыщенность жидкостью, относительная проницаемость газа, относительная проницаемость нефти в газе при водонасыщенности на уровне связанной воды и капиллярное давление нефть/газ. Относительные проницаемости и капиллярное давление – функции от насыщенности жидкостью.
Второй набор ключевых слов функций насыщенности
Остается на самостоятельное изучение. Он реже используется (например, заводнение со смешиванием (Miscible flood) или использование функции растворитель (Solvent option), а также в случае 2-хфазного запуска без нефти (водогазового расчета))
Слайд 172

Функции насыщенности и масштабирование концевых точек • SWL – насыщенность связанной

Функции насыщенности и масштабирование концевых точек

• SWL – насыщенность связанной воды,

обозначаемая Swco. Это наименьшее значение насыщенности воды для любой рассматриваемой функции насыщенности воды и часто определяется как неснижаемая водяная насыщенность.
• SWCR – критическая насыщенность воды, обозначаемая Swcr. Это наибольшее значение насыщенности воды во всех таблицах для неподвижной воды (krw =0).
• SWU – максимальная насыщенность воды, обозначаемая Swu. Это наибольшее значение насыщенности воды для любой рассматриваемой функции насыщенности воды.
Слайд 173

Функции насыщенности и масштабирование концевых точек • SGL - начальная насыщенность

Функции насыщенности и масштабирование концевых точек

• SGL - начальная насыщенность газа,

обозначаемая Sgco. Это наименьшее значение насыщенности газа для любой взятой функции насыщенности газа.
• SGCR – критическая насыщенность газа, обозначаемая Sgcr. Это наибольшая насыщенность газа во всех таблицах для неподвижного газа (krg =0).
• SGU – максимальная насыщенность газа, обозначаемая Sgu. Это наибольшее значение насыщенности газа для любой взятой функции насыщенности газа.
• SOWCR – критическая насыщенность нефти в воде, обозначаемая Sowcr. Это наибольшее значение насыщенности нефти, при котором нефть неподвижна в воде.
• SOGCR – критическая насыщенность нефти в газе, обозначаемое Sogcr. Это наибольшее значение насыщенности нефти, при котором нефть неподвижна в газе.
Слайд 174

Функции насыщенности и масштабирование концевых точек • SGL - начальная насыщенность

Функции насыщенности и масштабирование концевых точек

• SGL - начальная насыщенность газа,

обозначаемая Sgco. Это наименьшее значение насыщенности газа для любой взятой функции насыщенности газа.
• SGCR – критическая насыщенность газа, обозначаемая Sgcr. Это наибольшая насыщенность газа во всех таблицах для неподвижного газа (krg =0).
• SGU – максимальная насыщенность газа, обозначаемая Sgu. Это наибольшее значение насыщенности газа для любой взятой функции насыщенности газа.
• SOWCR – критическая насыщенность нефти в воде, обозначаемая Sowcr. Это наибольшее значение насыщенности нефти, при котором нефть неподвижна в воде.
• SOGCR – критическая насыщенность нефти в газе, обозначаемое Sogcr. Это наибольшее значение насыщенности нефти, при котором нефть неподвижна в газе.
Слайд 175

Функции насыщенности и масштабирование концевых точек • Целью масштабирования концевых точек

Функции насыщенности и масштабирование концевых точек

• Целью масштабирования концевых точек является

обеспечение относительно небольшого количества общих функций насыщенности
• Функции насыщенности изменяются и адаптируются к существующим типам пород
• Существуют три основных способа масштабирования:
1) Масштабирование по оси насыщенности с изменением относительной проницаемости и капиллярного давления (“масштабирование конечной точки”)
2) Масштабирование относительной проницаемости (“вертикальное масштабирование”)
3) Масштабирование капиллярного давления
• Все функции насыщенности – это своего рода псевдофункции, т.к. все они зависят от масштабирования
• Масштабирование функций насыщенности – это самый простой метод создания псевдофункций.
Слайд 176

Лекция 7 СЕКЦИЯ REGIONS Гидродинамическое моделирование

Лекция 7 СЕКЦИЯ REGIONS

Гидродинамическое моделирование

Слайд 177

Назначение раздела REGIONS • Раздел REGIONS необязателен (опционален) • Раздел REGIONS

Назначение раздела REGIONS

• Раздел REGIONS необязателен (опционален)
• Раздел REGIONS делит пласт

на несколько областей в зависимости от характеристик пласта или целей расчета
• Большинство ключевых слов раздела REGIONS имеют форму XXXNUM
• Для каждого типа области ячейке присваивается номер области от 1 до максимально возможного номера области.
• Несколько ключевых слов XXXNUM не относятся к разделу REGIONS
• Ключевое слово REGDIMS в RUNSPEC используется для установления верхнего предела количества областей.
Слайд 178

Назначение раздела REGIONS Раздел REGIONS используется: • Для назначения особых свойств

Назначение раздела REGIONS

Раздел REGIONS используется:
• Для назначения особых свойств или характеристик

ячейкам или группам ячеек. Свойства задаются таблицами в файлах данных ECLIPSE.
• Для составления отчета о начальных запасах флюидов в определенных частях пласта.
К примеру, если пласт содержит три различные нефти, будут представлены три набора таблиц PVT, по одной на область. Предполагается, что таблицы пронумерованы в порядке введения в файл данных. Каждой ячейке сетки присваивается номер от 1 до 3, в зависимости от того, какую таблицу мы используем в этой ячейке. По умолчанию номер области для каждой ячейки равен 1. Этот номер относится как к номеру области, так и к номеру PVT (или PVTNUM – ключевому слову, используемому для задания этих значений в каждой ячейке). Ячейка, которой присвоен номер 1 будет использовать первую таблицу PVT. Тем, которым присвоен номер 2, будут использовать вторую таблицу, и тем, которым присвоен номер 3 - третью. Если существует только одна таблица PVT, не нужно обозначать области PVT, т.к. 1 по умолчанию присваивается каждой ячейке. В этом случае раздел REGIONS не нужен и может быть исключен. Различные типы областей, назначенные для каждой ячейки сети, прикрепляют эту ячейку к различным типам таблиц.
Слайд 179

Типы ключевых слов раздела Regions Определение регионов XXXNUM e.g. PVTNUM Вывод

Типы ключевых слов раздела Regions

Определение регионов XXXNUM e.g. PVTNUM
Вывод в отчет

начальных запасов по региону FIPNUM
Дополнительные области FIPXXXXX
Ориентированные ключевые слова KRNUMX/X-/Y/Y-Z/Z-
Операторы EQUALS, ADD, COPY, BOX,ENDBOX
Контроль вывода RPTREGS, BOUNDARY, INIT, RPTSCHED
Исключения FLUXNUM, RESVNUM, NINENUM, PINCHNUM, которые указываются в секции GRID

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

Слайд 180

Типы ключевых слов раздела Regions Ключевые слова, определяющие области Наиболее часто

Типы ключевых слов раздела Regions

Ключевые слова, определяющие области
Наиболее часто используются EQLNUM,

PVTNUM и SATNUM.
EQLNUM связывает балансирующие регионы с ячейками или группами ячеек.
PVTNUM используется для связи таблиц PVT с определенными ячейками.
Т.к. каждая область PVT содержит флюиды различной плотности, следовательно они должны быть сбалансированы отдельно. Таким образом каждая область PVT должна соотносится с балансирующей областью, даже если все области PVT имеют одинаковый набор межфлюидных контактов. Однако наоборот, один набор фаз может иметь различные контакты. Каждый из множества межфлюидных контактов определяется отдельной записью в EQUIL, отделенной прямым слэшем (/). Двойные записи данных межфлюидных контактов могут быть определены по умолчанию после первой записи в EQUIL.
Слайд 181

Области с начальными запасами флюидов Например, для модели «мертвая» нефть с

Области с начальными запасами флюидов

Например,
для модели «мертвая» нефть с водой:
PROPS
PVDO

1> /
/
/
/
PVTW
/
/ table 2 по умолчанию равна table 1
/
RSCONSTT
/
/ table 2 по умолчанию равна table 1
/
REGIONS
PVTNUM
/
/
COPY
‘PVTNUM’ ‘EQLNUM’ /
/
SOLUTION
EQUIL
/
/ table 2 по умолчанию равна table 1
/

Количество областей PVTNUM и EQLNUM (в данном случае две) устанавливается с использованием NTPVT в TABDIMS и NTEQUL в EQLDIMS (секция RUNSPEC), соответственно. SATNUM используется для связи функций насыщения с ячейками или группами ячеек. Количество областей SATNUM устанавливается с использованием NTSFUN в TABDIMS.

Слайд 182

Области с начальными запасами флюидов Области с начальными запасами флюидов могут

Области с начальными запасами флюидов

Области с начальными запасами флюидов могут определяться

с использованием FIPNUM, которое задает целые значения для каждой ячейки текущей области. Дополнительные области с начальными запасами флюидов могут задаваться с использованием FIPXXXXX.
Области FIPNUM и FIPXXXXX могут частично пересекаться. Общее количество этих областей должно быть задано NTFIP в REGDIMS раздела RUNPSEC.
К примеру, рассмотрим модель 20*5*10. Области от I=1 до 10 и от I=11 до 20 являются отдельными блоками, ограниченными сбросами, о начальных запасах флюидов которых необходимо составить отдельные отчеты, также необходимо составить отчеты о начальных запасах флюидов в каждом слое. Эти области будут определены так:

В разделе RUNPSEC
REGDIMS
--NTFIP NMFIPR
10 2 /
И в разделе REGIONS
EQUALS
--Array Value I1 I2
'FIPNUM' 1 /
'FIPNUM' 2 11 20 /
/
FIPLAYER
--Fiplayer - FIPXXXXX
--имеет значение от 1 до 10
--FIPLAYER=1 for layer 1
100*1
--FIPLAYER=2 for layer 2
100*2
--FIPLAYER=3 for layer 3
100*3
--FIPLAYER=4 for layer 4
100*4
--FIPLAYER=5 for layer 5
100*5
--FIPLAYER=6 for layer 6
100*6
--FIPLAYER=7 for layer 7
100*7
--FIPLAYER=8 for layer 8
100*8
--FIPLAYER=9 for layer 9
100*9
--FIPLAYER=10 for layer 10
100*10 /

Слайд 183

Лекция 8 СЕКЦИЯ SOLUTION Гидродинамическое моделирование

Лекция 8 СЕКЦИЯ SOLUTION

Гидродинамическое моделирование

Слайд 184

Области с начальными запасами флюидов • Раздел SOLUTION определяет начальные условия

Области с начальными запасами флюидов

• Раздел SOLUTION определяет начальные условия моделирования
•

Начальные условия могут быть определены путем балансировки, перечисления или
перезапуска (restart)
• Аналитические подошвенные воды задаются в разделе SOLUTION
• Данные раздела SOLUTION могут быть введены в файлах PRT и файлах перезапуска

Переходная газо-нефтяная зона

Переходная водо-нефтяная зона

Нефть

Слайд 185

Назначение раздела SOLUTION Раздел SOLUTION определяет начальные условия моделирования. Такие как:

Назначение раздела SOLUTION

Раздел SOLUTION определяет начальные условия моделирования. Такие как:
• начальное

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

Балансировка • При балансировке определяются контакты и данные давления и глубины

Балансировка

• При балансировке определяются контакты и данные давления и глубины
• Предполагается

гидростатическая балансировка
• Балансировка подходит для изначально неразбуренного пласта с плоскими контактами
• Там может быть входной уровень капиллярного давления
• EOS темной нефти используется для расчета гидростатического давления каждой фазы
• Фазовое насыщение в каждой зоне берется из функции насыщения в разделе PROPS
• В переходных зонах Sw и Sg устанавливается путем обратного просмотра кривых капиллярного давления, представленных как часть функций насыщения.
Слайд 187

Балансировка Опция балансировки для установки начальных давлений и насыщений основана на

Балансировка

Опция балансировки для установки начальных давлений и насыщений основана на функциях

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

Здесь (r) и (s) определяют пластовые и поверхностные условия, соответственно,
ρo – плотность нефти, ρg – плотность газа. Задав контакты флюидов, можно решить уравнение состояния для газа и воды схожим образом, определив начальное гидростатическое давление каждой фазы в любом месте пласта на момент добычи. На практике, ECLIPSE рассчитывает фазовое давление на 100 глубинах, равномерно распределенных по пласту. Это может быть изменено, путем перенастройки NDPRVD в ключевом слове EQLDIMS раздела RUNPSEC

и гидростатическое давление нефтяной фазы

Слайд 188

Начальная фаза насыщения Если начальное давление рассчитано, ECLIPSE задает фазовое насыщение

Начальная фаза насыщения

Если начальное давление рассчитано, ECLIPSE задает фазовое насыщение в

каждой зоне. Это не относится к переходным зонам, которые рассматриваются ниже. Фазовое насыщение берется из граничного насыщения, определенного функциями насыщения в разделе PROPS, а нефтяное насыщение рассчитывается на основе насыщения воды и газа. Поэтому:
• В газовой зоне Sg максимально. Это SGU, наибольшее газонасыщение во вводимой функции газового насыщения. Sw минимально, SWL, наименьшее значение водонасыщенности во вводимой таблице функций насыщения воды, т.е. начальное или минимальное водяное насыщение Swco. Нефтяное насыщение: So=1-SGU-SWL.
• В водяной зоне Sw максимально. Это SWU, наибольшее водонасыщение во вводимой функции водяного насыщения. Часто это значение равно 1. Sg минимально, SGL, наименьшее значение газового насыщения во вводимой таблице функций насыщения газа, т.е. начальное газовое насыщение Sgco. Оно обычно нулевое. Нефтяное насыщение So=1-SWU-SGL. Обычно So=1-Swco.
• В нефтяной зоне, значения воды и газа минимальны. Sg=SGL=Sgco и Sw=SWL=Swco. Нефтяное насыщение So=1-SWL-SGL. Обычно, So=1-Swco.
Слайд 189

Начальное фазовое насыщение в переходной зоне В переходной зоне фазовое насыщение

Начальное фазовое насыщение в переходной зоне

В переходной зоне фазовое насыщение зависит

от капиллярного давления.

и

Которое может быть легко подсчитано на любой глубине, как разница между гидростатическим давлением фаз. ECLIPSE затем использует обратный просмотр кривых капиллярного давления в водимых функциях насыщения для определения переходных зон насыщения газа и воды. Теоретически, кривая капиллярного давления воды при Pcow→∞ стремится к Sw=0. На практике, кривая капилярного давления должна заканчиваться в Swco для уверенности, что у нее нет разрыва в Sw при наименьшем водяном насыщении в функции насыщения.
Рассмотрим наклонную ячейку, пересекающуюся с водо-нефтяным контактом.
Насыщение ячейки задается путем расчета глубины центра ячейки и задания
начальных Soi, Swi и Sgi для центра ячейки, как показано на рисунке. Этот метод
позволяет достичь стабильного начального решения, в котором фазовое
насыщение согласуется с гидростатическим давлением. Этот метод известен как
центроблоковая балансировка.

Слайд 190

Начальное фазовое насыщение в переходной зоне

Начальное фазовое насыщение в переходной зоне

Слайд 191

Начальное фазовое насыщение в переходной зоне • В двухфазной модели вода-нефть

Начальное фазовое насыщение в переходной зоне

• В двухфазной модели вода-нефть или

газ-вода газо-нефтяной контакт (GOC) игнорируется
• В двухфазной модели нефть-газ водо-нефтяной контакт (OWC) игнорируется
• В конденсатной модели, изначально выше точки росы, OWC должно быть установлено на той же глубине, что и GOC
• Если имеется только свободный газ с испаренной нефтью, GOC может находиться ниже подошвы месторождения
• При отсутствии подвижной воды изначально OWC может находиться ниже подошвы месторождения
Слайд 192

Перечисление • Начальные условия могут быть заданы явно • Начальное давление

Перечисление

• Начальные условия могут быть заданы явно
• Начальное давление должно быть

установлено для каждой ячейки
• Начальное насыщение каждой фазы, кроме нефтяной, должно быть установлено для каждой ячейки
• Начальные соотношения растворов (Rs и/или Rv) должны быть определены в каждой ячейке
• Капиллярное давление и фазовое насыщение должны быть согласованы для обеспечения стабильных начальных условий
• Все остальные существенные количественные признаки растворов, такие как концентрация трассеров, должны быть также определены
• Перечисление может применяться в пластах с изначально наклонными контактами или несбалансированных
Слайд 193

Начальное соотношение растворов • Концентрация растворенного в нефти газа, Rs или

Начальное соотношение растворов

• Концентрация растворенного в нефти газа, Rs или ее

изменения с глубиной, т.е. RS или RSVD
• Концентрация испаренной нефти, Rv или ее изменения с глубиной, т.е. RV или RVVD
• Изменения с глубиной точки насыщения и/или точки росы, т.е.. PBVD и/или PDVD
- Предыдущая
Птицы и звери зимой
Следующая -
Материалы для рисования. Колорит и перспектива

Похожие презентации

Мемориальный комплекс Три штыка: подвиг народа бессмертен
Отрасли животноводства
9 мая
на 29.08.22
Элементы зонной теории
20131017_zhivotnye._predlog_iz-za
Механизмы и процессы. Промышленные роботы фанук и процесс сварки
Итак 90-е. Это было
Зелень Шееле. Чистые субботы Список приёма сырья
20140803_veselye_gnomiki
Какие СМИ я смотрю/слушаю/читаю
маршруты для Диадока
Порядок и примеры заполнения документов валютного контроля
Дизайн и моделирование одежды
Разходи за туристически пътувания на населението
Синяя полоса
Музыка и живопись (музыкально-педагогический аспект)
Виды замеса макаронного теста, особенности технологических расчетов
Ag Artificial Gravity oxword-Hogvarts-Moscow ge geneo
ПЕДСОВЕТ 19.01.2022 (ЗАВ. ОТДЕЛЕНИЕ)
Кабель витая пара. Лекция 5
Модернизация агрегатного участка СТО Автомир
Однофазный однополупериодный выпрямитель. Однополупериодный инвертор. Преобразователи частоты в электроприводе
Балет
Понятие здания, сооружения. Внутренее пространство зданий
Патология наследственности
Синквейны Венцловский
Занятие мои помощники правила
Вы можете прислать нам Вашу презентацию и мы опубликуем ее на сайте.
Обратная связь
Для правообладателей
Email: Нажмите что бы посмотреть