4.1.2 Построение структурных моделей залежей

Построение структурной модели проводилось в три этапа.

На первом этапе была создана опорная модель горизонтов. При этом использовались структурные поверхности отражающих горизонтов, построенные по результатам интерпретации 3D сейсморазведки. Для верейского и башкирского объекта использовались поверхности по отражающим горизонтам ОГ IIв (кровля отложений верейского горизонта) и ОГ IIб (кровля отложений башкирского яруса). При построении модели объектов нижнего карбона были использованы поверхности ОГ II (кровля тульских терригенных отложений), ОГ IIп (кровля черепетского горизонта турнейского яруса), а также поверхности по кровле бобриковских и малевско-упинских отложений.

На втором этапе на основе опорной модели создана модель изохор для всех моделируемых пластов. В качестве исходных данных использованы результаты интерпретации материалов ГИС. Следующим шагом явилось создание детальной модели горизонтов. В качестве базовой послужила созданная ранее опорная модель горизонтов, в качестве исходных данных –модель изохор и отбивки пластов по скважинам. На рисунках 4.1 и 4.2 приведены структурные модели отложений среднего и нижнего карбона соответственно.

Рисунок 4.1 –Структурная модель отложений среднего карбона

Рисунок 4.2 –Структурная модель отложений визейского объекта

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

Литологическое моделирование проводилось на основе данных интерпретации ГИС по зонам между геологическими границами кровли и подошвы пластов, выделенными в процессе корреляции. Все породы были разделены на литотипы: «коллектор» и «неколлектор». На первом этапе литологического и петрофизического моделирования произведено осреднение скважинных данных (Block wells). Скважинные данные содержат следующие кривые, необходимые для построения 3D моделей: дискретную кривую литологии, непрерывные кривые пористости и нефтенасыщенности. Кривая литологии создана таким образом, что интервалам, соответствующим литотипу «коллектор», присвоено значение 1, а «неколлектор» –0.

Второй этап литологического моделирования заключался в распределении пород на сетку модели. Эта операция выполнена с применением инструмента RMS Parameter Interpolation. Для всех моделируемых пластов использовался детерминистический метод интерполяции непрерывной кривой литологии, с дальнейшей дискретизацией по граничному значению 0,5 на коллектор-неколлектор. Результат распространения коллекторов представлен на рисунках 4.3 и 4.4.

Рисунок 4.3 –Вертикальный срез куба литологии (верейский+башкирский объект)

Рисунок 4.4 –Вертикальный срез куба литологии (визейский объект)

Петрофизическое моделирование проводилось только в пределах коллекторов. Пористость распространялась посредством Parameter Interpolation. Исходной информацией для моделирования пористости по всем пластам явились кривые результатов интерпретации ГИС по скважинам, осредненные на ячейки. В качестве примера приведено распространение пористости коллекторов по разрезу визейского объекта (рис. 4.5).

Рисунок 4.5 –Вертикальный срез куба пористости (визейский объект)

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

Таблица 4.2 – Петрофизические зависимости для расчета проницаемости

Объект	Петрофизическая зависимость Кпр = f (Кп)
B-II	Кпр=0,0119*exp(48.213*Кп)
B-III	Кпр=0,0119*exp(48.213*Кп)
A4	Кпр=6,868*105*Кп4,767
C-I+II ... C-VI	Кпр=8.144*(106)*(Кп6.777)
Ct-II+III, Ct-IV	Кпр=0,0082*exp(52.567*Кп)

Пример распространения проницаемости по разрезу представлен на рисунке 4.6.

Рисунок 4.6 –Вертикальный срез куба проницаемости (визейский объект)