Прямые поиски залежей углеводородов в последнее время приобретают всевозрастающее значение в связи с тем, что объектом направленного поиска с помощью геофизических методов становятся средние и малые месторождения на значительных глубинах. Проверка с помощью бурения перспективных ловушек такого типа становится все менее рентабельной.
Изучение зависимости амплитуд отраженных волн от угла наклона лучевых траекторий (или от расстояния источник-приемник) проводится уже давно [35,37,39,41,44]. Был разработан метод AVO (Amplitude Versus Offset), что можно перевести как «зависимость амплитуды от расстояния источник-приемник». Далее фразу «расстояние источник-приемник» будем для краткости обозначать «РИП». Этот метод позволяет с достаточной точностью определять тип порового флюида коллектора по характеру изменения амплитуд с изменением угла падения волн на отражающие границы в продуктивной части разреза (как правило, на подошву и кровлю пласта-коллектора). Однако никто не изучал эту зависимость для карбонатов и доломитов, насыщенных газом, нефтью или водой. Все авторы ограничивались рассмотрением только терригенных пород. Между тем, приводимая здесь табл. 1 доказывает, что на долю карбонатных коллекторов (даже без учета смешанных коллекторов) приходится более трети всего углеводородного сырья. Поэтому задача распространения уже разработанного метода AVO на карбонатные коллекторы представляется.
Антологические Нефть Газ типы пород ж—* ~0 л s л л Табл. 1. Соотношение ресурсов.
1 ерригенные Zo.4 о4.1 углеводородов в коллекторах Карбонатные 44.0 28.1 различных литологических типов,.
Смешанные 27.6 7.8% [11] нам важной и актуальной.
Научная новизна данной работы заключается в том, что в ней изучаются карбонатные коллекторы, отличающиеся значительной неоднородностью по сравнению с терригенными, а также анизотропией скорости, что существенно усложняет решение поставленных задач. Полученные результаты могут быть использованы для проведения прямых поисков залежей УВ в карбонатных коллекторах, так как позволяют делать заключения о типе порового флюида (газ или вода), насыщающего горные породы, только с помощью анализа амплитудных кривых в ОГТ-выборках.
Цель работы — обосновать возможность применения для прямых поисков залежей углеводородов в карбонатных коллекторах сейсмического метода ОГТ (в модификации АЛ/О), основанном на детальном изучении динамических особенностей отраженных волн.
Задачи диссертации формулируются следующим образом.
1. Провести исследования зависимости амплитуд отраженных волн от РИП для конкретного карбонатного коллектора.
2. Разработать компьютерную программу для моделирования А/0~ эффектов в анизотропных упругих толстослоистых средах. Под толстыми слоями будем понимать такие слоя, мощность которых превышает длину сейсмической волны.
3. Проверить с помощью математического моделирования правильность проведения амплитудного анализа полевого материала.
4. Сделать выводы о перспективности поисков залежей УВ в карбонатных коллекторах с помощью А/0-анализа.
Для решения первой задачи была взята сейсмическая информация, полученная в процессе проведения полевых работ (АО «Волгограднефтегеофизика») по методу ОГТ на профиле 860 529 Романовской площади (район Доно-Медведицких поднятий). Выбор профиля обусловлен наличием на сейсмическом временном разрезе перспективной антиклинальной структуры по продуктивному горизонту трещиноватые известняки воронежского горизонта, D3fr), а также тем, что скважина 20-Нижнекоробковская, пробуренная вблизи этого профиля, дала воду с примесью нефти из пласта-коллектора. На сейсмограммах О. ГТ изучались амплитуды осей синфазности, связанных с волнами, отраженными от карбонатного коллектора. Было взято несколько сейсмограмм, чьи средние точки находились как над предполагаемой залежью, так и вдали от нее. Полученные в результате соотношения показали, что характер зависимости амплитуды от РИП имеет разную форму над предполагаемой залежью и вдали. Было сделано предположение, что это связано в большей степени с изменением в коллекторе типа порового флюида (с газа на воду) вдоль профиля, хотя допускается также и более слабое влияние смены литологии и коэффициэнта пористости.
Вторая задача — составление компьютерной программы для двумерного математического моделирования сейсмических эффектов в анизотропных средах — была решена с помощью системы визуального программирования Delphi 3.0. Программа написана на языке Object Pascal 7.0. Она обладает возможностью рассчитывать с точностью нулевого приближения лучевого метода кинематические и динамические характеристики отраженных волн, распространяющихся в упругой среде, состоящей из небольшого (10−20) количества «толстых» анизотропных (трансверсально-изотропных) слоев. При этом параметры слоев могут изменяться вдоль профиля.
Для решения третьей задачи была построена двумерная модель среды по ПР 860 529 Романовской площади. При этом учитывались как возможная анизотропия слоев, так и изменение их упругих характеристик по латерали. Исходной информацией для этого послужили: временной разрез по профилю, данные сейсмического и акустического каротажа в скважине 20-Нижнекоробковская.
Результатом моделирования стали вычисленные зависимости амплитуды отраженных волн на синтетических сейсмограммах ОГНГ в тех точках профиля, которые соответствуют моделируемой залежи и вдали от нее.
Решение первых трех задач позволило решить и последнюю. При сравнении экспериментальных и теоретических А/0~ зависимостей обнаружилось их значительное сходство. Был сделан вывод о том, что проведение исследований амплитуд волн, отраженных от границ карбонатного коллектора имеет большое практическое значение, так как позволяет делать предположения о типе флюида, заполняющего поры карбонатного коллектора. Это, в свою очередь, позволит проводить прямые поиски залежей УВ.
Таким образом, основными защищаемыми положениями диссертации являются:
1. Доказана эффективность применения метода АЛ/О для обработки сейсмических данных ОГТ с целью прямых поисков залежей углеводородов в карбонатных коллекторах.
2. Разработана компьютерная программа для моделирования полей отраженных волн в анизотропных горизонтально-неоднородных карбонатных коллекторах.
3. Доказана возможность применения метода А/0 для определения оптимальных параметров суммирования по ОГТ при исследовании карбонатных коллекторов.
Работа состоит из введения, пяти глав и заключения. В 1-й главе «Обзор работ, связанных с проблемой изучения зависимости амплитуд отраженных волн от расстояния источник-приемник» проводится обзор исследований по проблеме прямых поисков залежей углеводородов с помощью АЛ^О-анализа. Как можно увидеть, все исследования проводились только для терригенных коллекторов.
В процессе решения поставленных задач мы натолкнулись на определенные трудности, связанные с петрофизическими особенностями карбонатных пород. Как показано во П-й главе «Выбор адекватной математической модели карбонатного разреза», известняки и доломиты очень часто характеризуются существенной изменчивостью физических свойств (для нас особенно интересны скорости продольных и поперечных волн). Карбонаты нередко имеют различные скорости распространения сейсмических волн в матрице породы по различным направлениям, то есть обладают свойством анизотропии скоростей. Это вызывается закономерными изменениями физических свойств матрицы породы (в результате выщелачивания, доломитизации, каль-цитизации и др.), раздроблением породы системами трещин. Предлагается использовать в качестве математической модели карбонатного разреза совокупность «толстых» трансверсально-изотропных однородных слоев. С помощью такой модели можно учитывать не только тонкую слоистость разреза, но и существенную неоднородность карбонатов.
В III-й главе «Теория и алгоритмы компьютерной программы математического моделирования волновых процессов в анизотропных средах» для математического моделирования кинематических и динамических особенностей отраженных волн в анизотропных средахприводятся основные теоретические положения и алгоритмы программы, разработанной автором и использованной им для исследований.
В ГУ-Й главе «Применение А0-анализа к полевым материалам сейсмического метода ОГТ на примере Нижнекоробковского поднятия (Доно-Медведицкий вал)» дается описание конкретного сейсмического полевого материала метода ОГТ. Производится обработка этой информации по методу А0 и выявляются особенности поведения графиков зависимости амплитуды отраженных волн от РИП.
— я глава «Проверка проведенных исследований с помощью математического моделирования» посвящена проверке и подтверждению исследований, описанных в IV главе. Для этого используется математическое моделирование с помощью разработанной автором программы. Делаются выводы относительно влияния типа насыщающего флюида на характер А/0-зависимости.
5.3. ВЫВОДЫ К ГЛАВЕ V.
Изучение амплитудных зависимостей с помощью математического моделирования позволило подтвердить правильность проведения А0-анализа над Нижнекоробковской структурой. При этом производился учет анизотропии скорости в карбонатных породах, а также вариаций физических свойств вдоль профиля вследствие изменения коэффициента пористости и типа насыщающего флюида. Значит имеется возможность проводить прямые поиски залежей углеводородов с помощью метода АЛ/О. Доказательством этого является факт совпадения реальных сейсмических материалов с синтетическими.
ЗАКЛЮЧЕНИЕ
.
Данная работа посвящена одному из методов прямых поисков залежей углеводородов по сейсмическим данным ОГТ. Этот метод в иностранной печати получил название АУО-анализ, суть которого заключается в изучении амплитуд отраженных волн в выборках ОГТ до суммирования. Автор показал, что А0-анализ является мощным инструментом для прямых поисков залежей У В. С его помощью можно по форме амплитудных кривых в выборках ОГТ различать коллекторы, которые содержат воду, от коллекторов, содержащих газ. Кроме того метод позволяет разделять аномалии, обусловленные газосодержащими песчаниками, и аномалии, чьим источником являются высокоскоростные пласты магматических пород. Прямой поиск нефтенасыщенных коллекторов также возможен, хотя и требует более тонкого анализа, так как амплитудные зависимости у таких коллекторов близки по форме и уровню с водонасыщенными.
АЮ-анализ позволяет выбирать правильные параметры суммирования для таких типов коллекторов, у которых наблюдается смена знака амплитуды отраженной волны в выборке ОГТ. При суммировании подобных выборок происходит взаимное уничтожение сигналов на разных удалениях. Поэтому на сейсмическом разрезе такие коллекторы имеют плохую амплитудную выраженность. А0-анализ позволяет определить максимальное расстояние источник-прибор, при котором деструктивный эффект во время суммирования происходить не будет.
В диссертации обосновывается, что наиболее полно петрофизические свойства карбонатного коллектора для математического моделирования могут быть описаны только при учете изменчивости свойств вдоль слоя, а также анизотропии скорости. Следовательно программа, которую можно использовать для математического моделирования должна обладать возможностью представлять среду в виде толстых (под толстым понимаем такой слой, чья мощность превышает длину регистрируемых сейсмических волн) трансверсально-изотропных неоднородных слоев. С помощью такого подхода можно учитывать не только тонкую слоистость разреза, но и существенную неоднородность карбонатов.
Существует немало возможностей для применения программы, рассчитывающей волновые поля в анозотропных средах. Например, для трещиноватых сред кинематические и, особенно, динамические особенности сейсмических волн могут быть смоделированы только с учетом анизотропных свойств пород. Разработанная автором программа обладает возможностью учитывать анизотропию скорости и изменчивость свойств слоев вдоль профиля. Вычисления производятся в рамках нулевого приближения лучевого метода, что делает программу эффективной и быстродействующей. Таким образом, предлагаемая автором программа может применяться для изучения АЛ/О-зависимостей для карбонатных коллекторов.
Результаты, полученные с помощью АЛ/О-анализа по данным метода ОГТ, проведенного на Романовской площади показали, что наблюдается существенное различие между формой амплитудных зависимостей для волн, отраженных от границ карбонатного коллектора (франские трещиноватые известняки) в сводовой части Нижнеко-робковской структуры и на некотором удалении от нее. Этот факт связывается автором со сменой порового флюида, насыщающего коллектор с газа в сводовой части на воду вне структуры.
Изучение амплитудных зависимостей с помощью математического моделирования позволило подтвердить правильность проведения А^/О-анализа над Нижнекоробковской структурой. При этом производился учет анизотропии скорости в карбонатных породах, а также вариаций физических свойств вдоль профиля вследствие изменения коэффициента пористости и типа насыщающего флюида. Значит имеется возможность проводить прямые поиски залежей углеводородов с помощью метода А/0. Доказательством этого является факт совпадения реальных сейсмических материалов с синтетическими.