Бесплатный автореферат и диссертация по наукам о земле на тему

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

ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

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

ИИ46И6879

На правах рукописи

ЕКИМЕНКО Антон Валерьевич

ПРОГНОЗ свойств ГЕОЛОГИЧЕСКОГО РАЗРЕЗА СРЕДНЕДЕВОНСКО-НИЖНЕФРАНСКОГО

НЕФТЕГАЗОНОСНОГО КОМПЛЕКСА ТИМАНО-ПЕЧОРСКОЙ НЕФТЕГАЗОНОСНОЙ ПРОВИНЦИИ ПО СЕЙСМИЧЕСКИМ ДАННЫМ

Специальность 25.00.10 - Геофизика, геофизические

методы поисков полезных ископаемых

Автореферат диссертации на соискание ученой степени кандидата геолого-минералогических наук

- 1 ШЛ 7070

САНКТ-ПЕТЕРБУРГ 2010

004606879

Работа выполнена в государственном образовательном учреждении высшего профессионального образования Санкт-Петербургском государственном горном институте имени Г.В.Плеханова (техническом университете).

Научный руководитель -

доктор геолого-минералогических наук

Егоров Алексей Сергеевич

Официальные оппоненты:

доктор геолого-минералогических наук, профессор

Кашубин Сергей Николаевич,

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

Киселёв Юрий Васильевич

Ведущая организация - ФГУНПП «Севморгео»

Защита диссертации состоится 16 июня 2010 г. в 15 ч на заседании диссертационного совета Д 212.224.01 при Санкт-Петербургском государственном горном институте имени Г.В.Плеханова (техническом университете) по адресу: 199106 Санкт-Петербург, 21-я линия, д.2, ауд.4312.

С диссертацией можно ознакомиться в библиотеке Санкт-Петербургского государственного горного института.

Автореферат разослан 14 мая 2010 г.

УЧЕНЫЙ СЕКРЕТАРЬ диссертационного совета кандидат геолого-минералогических наук, доцент И.Г.КИРБЯКОВА

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Необходимость развития методики комплексного анализа кинематических и динамических характеристик сейсмических волн определяется повышением требований к качеству геологической интерпретации данных сейсморазведки МОВ-ОГТ 2Б и ЗО в условиях увеличения объёмов детальных геофизических работ в пределах Тимано-Печорской нефтегазоносной провинции (ТПНГП).

Исследованиями потенциальных геофизических полей в регионе занимались Гафаров Р. А, Запорожцева И. В. Вопросы нефтегазогеологического районирования освещены в работах Аминова Л. 3., Белонина М. Д., Прищепы О. М. Геологическому изучению региона посвящены работы Никонова Н.И. Малышева Н. А., Тимонина Н. И, Богацкого В. И., Лебедева Б. А., Дедеева В. А и других.

Сейсморазведочные работы в северных районах ТПНГП проводятся со второй половины 60-х годов. К началу 70-х годов поисковыми работами Печорской геофизической экспедиции был исследован Колвинский мегавал. В 1975 г. к работе приступила Ленинградская сейсморазведочная экспедиция (ныне ОАО "Нарьян-Марсейсморазведка"), а с 1981 г. ПГО "Печорагеофизика" (ныне ОАО "Севергеофизика")- В настоящее время большая часть работ в рассматриваемом регионе проводится по методике пространственной сейсморазведки. Целью этих работ является не только структурная интерпретация, но и прогноз свойств пласта, картирование зон развития литолого-фациальных неоднородностей.

При относительно высокой геолого-геофизической изученности региона значительное число скважин оказываются непродуктивными (успешность бурения в регионе составляет в среднем 60%), что свидетельствует о необходимости разработки более надежных критериев прогноза свойств геологического разреза на основе детального анализа всей совокупности характеристик сейсмической записи. Наиболее перспективным направлением развития методики интерпретации сейсмических данных представляется прогноз параметров разреза на основе комплексной интерпретации кинематических и динамических характеристик волнового поля.

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

терригенных отложений среднедевонско-нижнефранского нефтегазоносного комплекса (НТК) в пределах Колвинского мегавала. Этот комплекс является одним из наиболее ресурсонасы щенных - залежи углеводородов (УВ) выявлены здесь на Харьягинском, Ошском, Лекхарьягинском и Инзырейском месторождениях.

Результаты исследований особенностей волнового поля и разработанные методические приёмы интерпретации сейсмической записи могут быть использованы при решении задач структурной и динамической интерпретации данных, как в пределах Колвинского мегавала, так и в других районах ТПНГП.

Цель работы. Повышение достоверности прогноза свойств геологического разреза на основе комплексной интерпретации динамических и кинематических характеристик волнового поля.

Идея работы. Геологическую интерпретацию сейсмических данных следует выполнять с привлечением всех характеристик сейсмической записи: динамических и кинематических, с учётом результатов моделирования волнового поля.

Основные задачи исследований:

• анализ данных сейсморазведки МОВ-ОГТ 30 и данных исследования скважин с целью определения геометрических, упругих и коллекторских свойств среднедевонско-нижнефранского НГК, создание сейсмогеологической модели среднедевонско-нижнефранского НГК территории Колвинского мегавала;

• расчёт и анализ синтетических волновых полей, установление связи между параметрами моделей и особенностями волнового поля, выработка методических приёмов геологической интерпретации сейсмических данных;

• уточнение геологического строения исследуемых площадей Колвинского мегавала по результатам комплексной интерпретации кинематических и динамических характеристик сейсмических данных МОВ-ОГТ 30.

Научная новизна:

• обоснованы геометрические, упругие параметры и коллекторские свойства сейсмогеологической модели среднедевонско-нижнефранского НГК;

• получены сейсмические образы (волновые поля), характеризующие изменения геометрических, упругих, коллекторских свойств разреза; установлены зависимости между геометрическими, упругими и коллекторскими свойствами терригенных отложений среднедевонско-нижнефранского НГК и динамическими и кинематическими характеристиками волнового поля

Защищаемые положения

1. Созданная на основе сбора, обработки и анализа геолого-геофизнческой информации сейсмогеологнческая модель отражает все наиболее значимые свойства геологического разреза средиедевонско-нижнефранского НГК территории Колвинского мегавала, а именно геометрические свойства, упругие характеристики пластов, коллекторские свойства.

2. Расчёт синтетических волновых полей позволяет выявить особенности разреза, существенно влияющие на волновое поле. Установлено влияние зоны выклинивания на создание скоростной модели разреза за счет образования здесь обменной волны; определена геолого-структурная природа AVO-аномалии; определён характер влияния коэффициента пористости пласта — коллектора на волновое поле.

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

Методика исследования. Для решения поставленных задач проведены обобщение и анализ результатов геофизических исследований скважин и данных сейсморазведки МОВ-ОГТ 3D, обеспечившие обоснование параметров сейсмогеологической модели среднедевонско-нижнефранского НГК. Синтетические волновые поля рассчитаны с использованием программы конечно-разностного моделирования волнового поля TESSERAL-2D (Tesseral Technologies Inc.). Обработка фактических материалов выполнена с использованием программ FOCUS, GEODEPTH (Paradigm

Geophysical).

Достоверность. Достоверность построения эффективной сейсмогеологической модели среднедевонско-нижнефранского НГК определяется использованием результатов геофизических исследований 13 скважин и данных по 200 км2 сейсмической съёмки МОВ-ОГТ 3D.

Эффективность предлагаемых приёмов интерпретации сейсмических данных определяется сопоставлением результатов сейсмогеологического моделирования с результатами обработки данных натурных наблюдений МОВ-ОГТ 3D.

Практическая значимость. Разработанные методические приёмы позволяют провести эффективную интерпретацию сейсмического материала, с целью прогноза свойств разреза среднедевонско-нижнефранского НГК и, таким образом, повысить надежность локализации участков недр наиболее благоприятных для обнаружения УВ.

Реализация результатов работы. Результаты, полученные в настоящей работе, использованы для интерпретации сейсмических данных компанией ОАО «Нарьян-Марсейсморазведка».

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

Апробация работы. Основные результаты, полученные автором, докладывались на международных конференциях: «Санкт-Петербург 2006», «Геомодель 2007», «Геомодель 2008»; «Геофизика 2007» и «Геофизика 2009», международном молодежном форуме «Проблемы недропользования» (СГПТИ (ТУ) 2007 г).

Публикации. По теме диссертации опубликовано 6 работ, в том числе 2 статьи в изданиях рекомендованных ВАК РФ.

Автор выражает искреннюю благодарность научному руководителю д. г.-м. н., профессору Егорову Алексею Сергеевичу. Автор признателен д. ф.-м. н., профессору Караеву Назиму Алигейдаровичу за помощь и консультации на протяжении всего

периода обучения в СПбГУ и аспирантуре Горного института.

Автор благодарен коллективу кафедры геофизических и геохимических методов поисков и разведки за участие в обсуждении основных результатов и советы по рассматриваемым в диссертации вопросам. Решение многих задач поставленных в работе велось при активном участии старшего преподавателя кафедры геофизики СПбГУ Кашкевича В. И., автор благодарен ему за ценные советы.

Автор признателен коллективу компании ОАО «Нарьян-Марсейсморазведка». Автор выражает благодарность руководителю ОАО «Нарьян-Марсейсморазведка» Гоцу И.А. и главному геологу Вискуновой К.Г. за созданные благоприятные условия для написания диссертации.

Структура и объем диссертационной работы.

Диссертационная работа состоит из введения, 4 глав и заключения, изложенных на 125 страницах машинописного текста, содержит 47 иллюстраций, 10 таблиц, библиографический список из 91 наименования.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Введение содержит общую характеристику работы. Первая глава посвящена геологическому обзору исследуемой территории, формулировке основных проблем связанных с изучением среднедевонско-нижнефранского НГК. Во второй главе, на основании обобщения данных ГИС, результатов исследования керна и данных сейсморазведки МОВ-ОГТ 3D обоснованы геометрические, упругие и коллекторские параметры сейсмогеологической модели геологического разреза среднедевонско-нижнефранского НГК. В главе обосновывается первое защищаемое положение. Третья глава содержит результаты расчета и анализа синтетических волновых полей, отражающих проявление геометрических, упругих и коллекторских свойств. В главе доказывается второе защищаемое положение. Четвёртая глава содержит результаты интерпретации реальных сейсмических данных, полученных на Ошской площади Колвинского мегавала. В главе рассмотрены вопросы калибровки сейсмических и скважинных данных, структурная и динамическая интерпретация. В главе доказывается третье защищаемое положение. Заключение содержит общие выводы по результатам исследования, основные

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

ОБОСНОВАНИЕ ЗАЩИЩАЕМЫХ ПОЛОЖЕНИЙ

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

Актуальность создания сейсмогеологической модели разреза среднедевонско-нижнефранского НГК обусловлена необходимостью иметь представление о свойствах разреза при сейсмических работах на площадях в разной степени освещенных результатами геологических и геофизических исследований.

Создание сейсмогеологической модели разреза предполагает обоснование следующих параметров необходимых при проведении сейсмических исследований: геометрические свойства разреза, упругие характеристики пластов, коллекторские свойства. Эти параметры определяют волновую картину, регистрируемую при проведении сейсмических наблюдений.

Геометрические параметры модели включают в себя глубины залегания пластов, мощности пластов, углы падения, тип границ между пластами. Глубины залегания и мощности определены на основе анализа стратиграфических отметок по скважинам, углы падения определялись по результатам сейсмических работ (анализ глубинных разрезов и карт изогипс). Мощности пластов изучаемого НГК колеблются от 50 до 100 метров. Углы падения, определенные по структурным картам, не превышают 3 градусов. Предположение о типе границ высказывается на основании анализа историко-геологических представлений. Разрез изучаемого НГК характеризуется наличием поверхности регионального размыва, что позволяет высказывать предположение о наличии в среде незеркальной (шероховатой) отражающей границы.

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

Помимо самих параметров были определены их соотношения, например соотношение скоростей продольных и поперечных волн. Это имеет большое значение при изучении площадей региона, недостаточно обеспеченных данными ГИС.

Установлено, что дифференциация геологического разреза по скоростям распространения волн значительна: пласты глин характеризуются низкими значениями скоростей (3400 м/с), песчаные пласты высокими значениями (4800 м/с). Значения плотности менее дифференцированы: пласты глин имеют плотность в среднем 2.2 г/см3, песчаные пласты характеризуются в среднем плотностью 2.4 г/см3. Для данного региона автором были установлены соотношения скоростей продольных и поперечных волн для пород разного состава.

Анализ данных по скважинам, расположенным в разных частях Колвинского мегавала, позволил установить, что значения скоростей сейсмических волн и плотностей горных пород, входящих в одни и те же литолого-стратиграфические комплексы, на разных площадях отличаются. В северной части ТПНГП значения плотностей и скоростей волн превышают значения, установленные в более южных регионах. Это явление объясняется большими глубинами залегания пород изучаемого комплекса на севере (более 4000 м).

Коллекторские свойства изучены на основе анализа фондовой литературы ОАО «Нарьян-Марсейсморазведка» и Тимано-Печорского научно исследовательского центра.

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

Обобщение этого комплекса информации позволяет обосновать параметры сейсмогеологической модели среднедевонско-нижнефранского НГК (рис. 1). Такая модель может быть использована при сейсмических исследованиях разных площадей Колвинского мегавала, поскольку включает весь комплекс необходимых параметров разреза: вариации глубин, мощностей и

Рис. 1. Сейсмогеологическая модель среднедевонско-нижнефранского НГК. 1 - переслаивание аргиллитов и алевролитов; 2 - аргиллиты; 3 - песчаники; 4 - переслаивание алевролитов и песчаников (джъерские отложения); 5 - доломиты; 6 - эрозионное срезание (шероховатая граница)

углов падения пластов; тип границ; коллекторские свойства; скорости распространения волн и плотности пород.

2. Расчёт синтетических волновых полей позволяет выявить особенности разреза, существенно влияющие на волновое поле. Установлено влияние зоны выклпнпвапия па создание скоростной модели разреза за счет образования здесь обменной волны; определена геолого-структурная природа AVO-аномалпи; определён характер влияния коэффициента пористости пласта - коллектора на волновое поле.

Для выявления' особенностей проявления целевого интервала в волновом поле в рамках данной работы были созданы четыре группы моделей. Каждая из групп, по сути, детализирует элементы общей сейсмогеологической модели (см. рис. 1): модель границы раздела двух сред (для определения влияния на волновое поле контакта разных пород); модель шероховатой границы (поверхность несогласия в геологическом разрезе); модель выклинивающегося пласта (выход под эрозионный срез отдельных пачек исследуемого НГК), модель изменения коллекторских свойств (изменение пористости «четвертой пачки» песчаников D2St-IV).

Расчет волнового поля проводился с использованием программы Tessera! 2D, которая рештизует расчет волнового поля методом конечных разностей. Анализировались кинематические и динамические характеристики сейсмической записи.

Моделирование границы раздела двух сред. В исследуемом регионе пласты песчаников среднедевонско-нижнефранского НГК перекрываются либо глинами (например, глинистые пачки Ошской площади), либо песчаниками (например, контакт четвертой пачки песчаников D2st-IV и песчаников джъерского горизонта Djdzr). Это обстоятельство обуславливает различия в волновом поле и необходимость изучения такого типа моделей.

С учётом того, что в производственной практике сейсморазведчикам часто приходится работать в условиях отсутствия полного комплекса ГИС, автором рассчитаны волновые поля для моделей, соответствующих разному уровню изученности методами ГИС: 1 - все упругие характеристики пласта известны достоверно по скважинным данным; 2 - отсутствует информация о

скоростях поперечных волн; 3 - отсутствует информация о скоростях поперечных волн и плотности.

Выполненный анализ амплитуд отраженных волн свидетельствует о том, что отсутствие достоверных данных об упругих параметрах, присущих конкретному региону, ведёт к ошибке в определении амплитудного отклика и делает невозможным применение методических приёмов основанных на интерпретации амплитуд волн.

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

Полноволновое моделирование позволило опробовать технологию AVO - анализа. По сейсмограммам были рассчитаны специальные атрибуты, которые характеризуют поведение амплитуд сейсмических волн в зависимости от удаления источник-приёмник. Анализ кроссплотов этих атрибутов позволил определить характерную форму кроссплотов для случая выклинивания пластов. Как видно из приведенного кроссплота (рис. 2 б), наблюдается разделение двух облаков точек - это связано с тем, что на разрезе имеет место интерференция отражений от кровли и подошвы клина. В случаях когда имеет место зашумленная волновая картина такое поведение точек может быть воспринято как аномалия, связанная с нефте- и газо-насыщением.

Моделирование шероховатой отражающей границы. Для исследования волновых полей, образуемых на шероховатой границе, автором были построены шесть моделей. Латеральный размер единичной неоднородности (L) задавался в пределах от 0 (гладкая граница) до 300 метров; вертикальный размер неоднородности задавался равным 1/4L. На основании анализа синтетических

волновых полей сделаны следующие выводы:

- в граф обработки следует обязательно включать процедуру миграции до суммирования;

- наличие неоднородностей размером до 0.5 долей длины волны не приводит к значительному изменению амплитуды отраженной волны и процедура миграции компенсирует это падение;

- неоднородности размером 0.5-1.6 длины волны влияют на амплитуды отраженных волн сильнее, чем неоднородности большего или меньшего масштаба (рис. 2 в); геометрия таких границ четко восстанавливается в процессе обработки, но падение амплитуд достигает 15 процентов;

- влияние неоднородностей размера большего 1.6 длины волны в процессе обработки полностью компенсируется, и такие неоднородности ведут себя как гладкие наклонные границы.

Таким образом, на основании анализа синтетических волновых полей были изучены особенности волнового поля, обусловленные элементами модели: показано влияние неполноты скважинной информации на динамическую интерпретацию; для модели шероховатой границы определено влияние размера единичной неоднородности на амплитуду отраженной волны, показано искажение скоростного вертикального спектра при наличии в среде шероховатых границ; для модели клина определены кинематические искажения вызванные образованием обменной волны, определено влияние геолого-структурного фактора на амплитуды (АУО аномалии).

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

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

а)

100

200

400

600

800

б)

в)

100-.....

^80 ь §

4 ?

5 &

I 160

< а.

40

..... *- ......-Ж.....

"......•----- -......

..

■X

>

ч

•Амш итуда *а нем 1грнра ванно.' 1 разр*

*Амшг !| 1а миг мфова^ном р азрезе

0 25а 256 50 75 100 150 300

Размер неоднородности границы (м.)

Рис. 2. Сейсмограмма ОПВ (а) и лучевая схема для модели выклинивания; кроссплот атрибутов интерсепт и градиент для модели выклинивания (б); график амплитуд волны отраженной от шероховатой границы (в).

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

Калибровка амплитуд сейсмической записи по скважинным данным. На современном этапе исследований ставится задача калибровки сейсмической записи путем сопоставления амплитуд на реальной сейсмограмме и амплитуд синтетической сейсмограммы, рассчитанной по скважинным данным (рис. 3 а). Эти амплитудные кривые совпадают при том условии, что все скважинные данные учтены, и обработка сейсмических данных проведена с сохранением амплитуд.

Анализ амплитудных кривых (для Ошской площади) демонстрирует согласованное поведение графиков амплитуд наблюденной сейсмограммы (кривая 1) с теоретической сейсмограммой, при расчете которой учтены обоснованные значения упругих параметров (кривая 2) (См. рис. 3 а). Использование при расчетах амплитудных характеристик волнового поля, полученных без учета региональных соотношений, приводит к значительным ошибкам в калибровке амплитуд (кривая 3) (См. рис. 3 а).

Таким образом, отсутствие достоверных данных ГИС или неучёт соотношений между упругими параметрами, характерных для данного региона, делает невозможным использование амплитуд сейсмических волн для геологической интерпретации сейсмических материалов. На основании анализа данных ГИС, полученных на территории Колвинского мегавала, установления соотношений между скоростями продольных и поперечных волн и плотностью горных пород, удалось точно выполнить калибровку сейсмических данных.

Обоснование структурной природы АУО-аномалии. На

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

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

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

В рассматриваемом случае, на основе анализа результатов инверсии, полученных на Ошской площади, было установлено, что высокие значения акустической жесткости приурочены к зоне выклинивания. Зоны повышенных значений акустической жесткости коррелируются с изолинией мощности клина 0.6 (X, - длина волны). Таким образом, стало возможно отделить область, где прогноз свойств по результатам инверсии достоверен от зоны неоднозначности.

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

а)

90

1)

¿. 70 гз

Ч

?> 50

с?

С

I

б)

---з

В)

«

1 ^ ,1 о

В" о н

о 3 250

X Щ о о га 240

со ь

1 « 230 о

^ 220 о о

п £ 210

С« о 200

«

ь 8 190 5

"То 15 20 25 30" Угол падения (градусы)

* Точки, полученные по синтетическим трассам .23 Точки, полученные по полевым данным ( над точкой номер скважины)

интерсепт (у.е.)

15 17

19

интерсепт (у.е.)

>

7 9 11 13 Кп <%)

Рис.3. Графики амплитуд отраженной волны, полученные при калибровке сейсмических данных

(а), кроссплоты атрибутов интерсепт и градиент для синтетических (вверху) и реальных (внизу) данных

(б), соотношение значений коэффициента пористости и амплитуды сейсмической записи (в).

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

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

ЗАКЛЮЧЕНИЕ

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

Проведенные в рамках данной работы обобщение и анализ геолого-геофизических данных по площадям центральной части Колвинского мегавала позволили создать эффективную модель строения среднедевонско-нижнефранского НГК. Эта модель может быть использована на других (в том числе неизученных бурением) площадях Колвинского мегавала.

Анализ результатов моделирования позволяет заключить следующее:

1) установлено, что недостаточная изученность площадей методами ГИС приводит к погрешностям, величина которых соизмерима с аномалиями, возникающими в реальных условиях в результате изменения литологии пород, их коллекторских свойств и характера насыщения. Для динамического анализа наблюденного сейсмического поля необходимы измеренные в скважине упругие параметры или использование зависимостей характерных для

данного региона;

2) определено влияние шероховатости границ на амплитуды отраженных волн;

3) установлены особенности проявления обменной волны, формирующейся в пласте - коллекторе среднего девона, в зоне выклинивания;

4) определена природа АУО - аномалий, полученных для модели выклинивающегося пласта; идентифицированы АУО - аномалии в реальном волновом поле среднедевонского интервала разреза;

Для Ошской площади обоснован следующий подход к определению акустических свойств пород по сейсмическим данным:

1) картирование линии выклинивания отложений перекрывающих основной пласт коллектор (БгвИУ);

2) картирование зоны изменения мощности перекрывающих отложений от 0 до 0.61 (X - длина волны);

3) в зоне изменения мощностей от 0 до 0.6Х для определения акустических свойств пород используются только данные ГИС и интерполяцию между скважинами;

4) в зоне мощностей более 0.6Х. используются статистические зависимости между сейсмическими и петрофизическими параметрами.

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

Результаты исследования могут быть использованы при интерпретации сейсмических материалов, направленных на изучение не только объектов в среднедевонско-нижнефранском НГК, но и объектов ордовикско-нижнедевонского НГК, обширные по площади зоны выклинивания которого расположены на восточном борту Хорейверской впадины. К зонам выклинивания нижнедевонских коллекторов приурочены значительные по запасам нефти месторождения им. Р.Требса и А.Титова, которые в настоящее время относятся к государственному резерву, но

планируются к введению в освоение уже в ближайшем будущем.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Екименко A.B. Использование динамической интерпретации сейсмических данных с целью прогноза коллекторских свойств в условиях карбонатного разреза. // Записки Горного института/ СПб. СПГГИ (ТУ) им. Г.В.Плеханова. 2009. Т. 182 С. 16-18.

2. Екименко A.B. Прогнозирование коллекторских свойств с использованием куба акустического импеданса. // Записки Горного института/ СПб. СПГГИ (ТУ) им. Г.В.Плеханова. 2009. Т. 183 С.235-237.

3. Екименко А. В. Использование динамических параметров отраженных волн при изучении среднедевонско-нижнефранского НГК Тимано-Печорской нефтегазоносной провинции. I Международная научно-практическая конференция молодых ученых и специалистов, посвященная памяти академика А. П. Карпинского, 24-27 февраля 2009 г. Тезисы докладов. СПб.: Изд-во ВСЕГЕИ, 2009. 534 с. Электронное издание на компакт-диске (CD).

4. Екименко A.B., Вазаева А. А. Эффективность выделения ловушек УВ сейсморазведкой MOB в северной части Тимано-Печорской нефтегазоносной провинции. Сборник тезисов докладов 9-й международной научно-практической конференции и выставки «Геомодель 2007», Геленджик 2007, С. 151. Электронное издание на компакт-диске (CD).

5. Вазаева А. А. Екименко A.B., Кузнецова И. Ф. Прогноз коллекторских свойств карбонатных пород по результатам динамической интерпретации в северной части Тимано-Печорской нефтегазоносной провинции. Геофизические методы исследования Земли и недр: Материалы 6-й международной научно-практической конференции «Геофизика-2007» - СПб.:С.-Петерб.ун-т, 2008. С. 5658.

6. Ekimenko А. V. Seismic Image of Earth Crust in the Northern Tatarstan. // Expanded absracts of EAGE&SEG Saint-Petersburg Conference, 16-19 October 2006, Saint-Petersburg, Russia. Электронное издание на компакт-диске (CD)

РИЦ СПГГИ. 07.05.2010. 3.253 Т. 100 экз. 199106 Санкт-Петербург, 21-я линия, д.2

Содержание диссертации, кандидата геолого-минералогических наук, Екименко, Антон Валерьевич

Введение

Глава 1. Геологическое строение Печоро-Колвинского авлакогена

Тимано-Печорской нефтегазоносной провинции (ТПНГП)

1.1 Структурно-тектоническое положение Печоро-Колвинского авлакогена в пределах ТПНГП

1.2. Литолого - стратиграфическая характеристика осадочного чехла

1.3. Нефтегазоносность Печоро-Колвинского авлакогена

Глава 2. Обоснование параметров сейсмогеологической модели разреза среднедевонско-нижнефранского НТК

2.1 Фактологическая база исследований

2.2. Обоснование структурно-вещественных параметров сейсмогеологической модели среднедевонско-нижнефранского НГК

2.2.1 Геометрические характеристики среднедевонско-нижнефранского НГК

2.2.2 Упругие параметры среднедевонско-нижнефранского НГК

2.2.3 Коллекторские свойства среднедевонско-нижнефранского НГК

Глава 3. Моделирование волновых полей и анализ их кинематических и динамических параметров

3.1 Моделирование раздела двух сред

3.2 Моделирование шероховатой отражающей границы

Кинематические особенности волнового поля шероховатой границы

Динамические особенности волнового поля

Некоторые особенности сейсмической записи от шероховатой границы

3.3. Модель выклинивающегося пласта

Кинематические особенности волнового поля

Динамические особенности волнового поля

Глава 4. Интерпретация сейсмических данных с использованием выявленных особенностей волнового поля

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

4.2 Определение положения линии выклинивания с привлечением амплитудных характеристик горизонта

4.3. Обоснование структурной природы AVO аномалий в зоне выклинивания

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

4.4 Определение коллекторских свойств целевого пласта

Введение Диссертация по наукам о земле, на тему "Прогноз свойств геологического разреза среднедевонско-нижнефранского нефтегазоносного комплекса Тимано-Печорской нефтегазоносной провинции по сейсмическим данным"

Актуальность темы

Актуальность комплексного анализа сейсмических данных МОВ-ОГТ 2D и 3D на основе детальных исследований особенностей распространения упругих волн определяется активным развитием поисковых, разведочных и добычных работ в пределах Тимано-Печорской нефтегазоносной провинции (ТПНГП), эффективность которых во многом зависит от детальности и качества их геологической интерпретации.

Провинция имеет более чем 150-летнюю историю геологического изучения. В XIX веке исследования носили преимущественно геолого-географический характер. Начиная с 20-х годов XX века Северной промысловой экспедицией проводилось изучение Тимана, Пай-Хоя, Вайгача и Новой Земли. С 1933 по 1940 годы в Арктическом регионе работали экспедиции Арктического института и Горно-геологического управления Главсевморпути. Начиная с 1947 года, в регионе проводятся планомерные комплексные геологосъёмочно-поисковые работы. Большая часть территории Малоземельной и Большеземельской тундр с начала 70-х до конца 80-х г.г. была изучена геологической съёмкой масштаба 1:50000. Поисковые сейсморазведочные работы в северных районах ТПНГП проводились со второй половины 60-х г.г. К началу 70-х г.г. поисково-детальными работами Печорской геофизической экспедиции был исследован Колвинский мегавал. В 1975 г. к работе в НАО приступила Ленинградская сейсморазведочная экспедиция (ныне ОАО "Нарьян-Марсейсморазведка"), а с 1981 г. ПГО "Печорагеофизика" (ныне ОАО "Севергеофизика").

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

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

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

В качестве объекта исследований выбран ареал развития терригенных отложений среднедевонско-нижнефранского нефтегазоносного комплекса (НТК) (в пределах Печоро-Колвинского авлакогена) - одного из наиболее ресурсонасыщенных в ТПНГП [5]. Залежи УВ в этих отложениях выявлены на крупном по запасам Харьягинском месторождении (6 залежей), Ошском (2 залежи), Лекхарьягинском и Инзырейском месторождениях.

Выполненный комплекс исследований на теоретических моделях и фактическом материале позволил: решить задачи структурной интерпретации, определить природу появления аномалий AVO (Amplitude variation with offset -изменение амплитуд с удалением) на сейсмограммах, определить ограничения сейсмического метода при проведении интерпретации динамических характеристик сейсмозаписи.

Цель работы. Повышение достоверности прогноза свойств геологического разреза на основе комплексной интерпретации динамических и кинематических характеристик волнового поля.

Идея работы. Геологическую интерпретацию сейсмических данных следует выполнять с привлечением всех характеристик сейсмической записи: динамических и кинематических, с учётом результатов моделирования волнового поля.

Основные задачи исследований:

• анализ данных сейсморазведки МОВ-ОГТ 3D и данных исследования скважин с целью определения геометрических, упругих и коллекторских свойств среднедевонско-нижнефранского НТК, создание сейсмогеологической модели среднедевонско-нижнефранского НГК территории Колвинского мегавала;

•расчёт и анализ синтетических волновых полей, установление связи между параметрами моделей и особенностями волнового поля, выработка методических приёмов геологической интерпретации сейсмических данных;

•уточнение геологического строения исследуемых площадей Колвинского мегавала по результатам комплексной интерпретации кинематических и динамических характеристик сейсмических данных МОВ-ОГТ 3D.

Научная новизна:

• обоснованы геометрические, упругие параметры и коллекторские свойства сейсмогеологической модели среднедевонско-нижнефранского НГК;

• получены сейсмические образы (волновые поля), характеризующие изменения геометрических, упругих, коллекторских свойств разреза; установлены зависимости между геометрическими, упругими и коллекторскими свойствами терригенных отложений среднедевонско-нижнефранского НГК и динамическими и кинематическими характеристиками волнового поля

Защищаемые положения

1. Созданная на основе сбора, обработки и анализа геолого-геофизической информации сейсмогеологическая модель отражает все наиболее значимые свойства геологического разреза среднедевонско-нижнефранского НГК территории Колвинского мегавала, а именно геометрические свойства, упругие характеристики пластов, коллекторские свойства.

2. Расчёт синтетических волновых полей позволяет выявить особенности разреза, существенно влияющие на волновое поле. Установлено влияние зоны выклинивания на создание скоростной модели разреза за счет образования здесь обменной волны; определена геолого-структурная природа AVO-аномалии; определён характер влияния коэффициента пористости пласта - коллектора на волновое поле.

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

Методика исследования. Для решения поставленных задач проведены обобщение и анализ результатов геофизических исследований скважин и данных сейсморазведки МОВ-ОГТ 3D, обеспечившие обоснование параметров сейсмогеологической модели среднедевонско-нижнефранского НГК. Синтетические волновые поля рассчитаны с использованием программы конечно-разностного моделирования волнового поля TESSERAL-2D (Tesseral Technologies Inc.). Обработка фактических материалов выполнена с использованием программ FOCUS, GEODEPTH (Paradigm Geophysical).

Достоверность. Достоверность построения эффективной сейсмогеологической модели среднедевонско-нижнефранского НГК определяется использованием результатов геофизических исследований 13

У скважин и данных по 200 км" сейсмической съёмки МОВ-ОГТ 3D.

Эффективность предлагаемых приёмов интерпретации сейсмических данных определяется сопоставлением результатов сейсмогеологического моделирования с результатами обработки данных натурных наблюдений МОВ-ОГТ 3D.

Практическая значимость. Разработанные методические приёмы позволяют провести эффективную интерпретацию сейсмического материала, с целью прогноза свойств разреза среднедевонско-нижнефранского НГК и, таким образом, повысить надежность локализации участков недр наиболее благоприятных для обнаружения УВ.

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Екименко, Антон Валерьевич

Заключение

Проведенные в рамках данной работы обобщение и анализ геолого-геофизических данных по площадям центральной части Колвинского мегавала позволили создать эффективную модель строения среднедевонско-нижнефранского НГК. Эта модель может быть использована на неизученных бурением площадях Колвинского мегавала.

С использованием математического аппарата, позволяющего рассчитывать распространение упругих волн в неоднородных средах и учитывать конкретные системы сейсмических наблюдений, используемые в настоящее время на практике, было проведено моделирование волновых полей набора теоретических моделей, имитирующих вариацию структурных и вещественных параметров целевого НГК.

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

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

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

1) Установлено, что в случае отсутствия измеренных по стволу скважин параметров Vs и р, использование их расчетных величин (без учёта характерных для региона соотношений) приводит к погрешностям, величина которых соизмерима с аномалиями, возникающими в реальных условиях в результате изменения литологии пород, их коллекторских свойств и характера насыщения. Этот факт указывает на необходимость ограниченного использования пересчетных параметров Vs и р при обработке и интерпретации современных сейсмических данных. Использование пересчетных параметров Vs и р допустимо только для кинематической интерпретации, на стадии проведения процедуры привязки волнового сейсмического поля к разрезу скважин; для динамического анализа наблюденного сейсмического поля необходимы измеренные в скважине упругие параметры или использование зависимостей характерных для данного региона.

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

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

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

3) Установлено наличие обменной волны, формирующейся в пласте-коллекторе среднего девона, и выявлены геологические условия, при которых волна оказывает искажающее влияние на наблюденное волновое поле.

4) Определена природа AVO — аномалий, полученных для рассчитанных моделей выклинивающихся пластов-коллекторов, идентифицированы AVO — аномалии в реальном волновом поле среднедевонского интервала разреза.

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

6) На примере среднедевонских терригенных отложений определено влияние изменения значений пористости коллекторов на параметры волнового поля и сделан расчет пористости для коллекторов на площади конкретного месторождения Колвинского мегавала.

7) В зонах выклинивания, являющихся в настоящее время основным объектом поисковых работ, обоснован следующий подход к определению акустических свойств пород по сейсмическим данным: i. картирование линии выклинивания; й. картирование зоны изменения мощности от 0 до 0.6 длины волны; iii. в зоне изменения мощностей от 0 до 0.6 длины волны для определения акустических свойств пород используются только данные ГИС и интерполяцию между скважинами; iv. в зоне мощностей более 0.6 длины волны используются статистические зависимости между сейсмическими и петрофизическими параметрами.

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

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

Выявленные особенности волнового поля могут быть использованы при интерпретации сейсмических материалов, направленных на изучение не только объектов в среднедевонско-нижнефранском НГК, но и объектов ордовикско-нижнедевонекого НГК, обширные по площади зоны выклинивания которого расположены на восточном борту Хорейверской впадины. К зонам выклинивания нижнедевонских коллекторов приурочены значительные по запасам нефти месторождения им.Р.Требса и А.Титова, которые в настоящее время относятся к государственному резерву, но планируются к введению в освоение уже в ближайшем будущем.

Библиография Диссертация по наукам о земле, кандидата геолого-минералогических наук, Екименко, Антон Валерьевич, Санкт-Петербург

1. Авербух А.Г. Изучение состава и свойств горных пород при сейсморазведке. М.: Недра, 1982, 232 с.

2. Аки К. и Ричарде П., Количественная сейсмология, М, Изд. «Мир», 1983.

3. Аминов JI.3., Баннов В.И., Боровинских А.П., Дмитриева Т.Е. и др. Минерально-сырьевой комплекс Республики Коми: проблемы и перспективы развития./ Сыктывкар.- 1999

4. Аминов JI.3., Белонин М.Д., Богацкий В.И., Боровинских А.П. и др. Методология и практика геолого-экономической оценки краевых систем древних платформ (на примере Тимано-Печорской нефтегазоносной провинции)/.- СПб.: ВНИГРИ, 2000

5. Белонин М.Д., Прищепа О.М., Теплов E.JL, Буданов Г.Ф., Данилевский С.А., Тимано-Печорская провинция: геологическое строение, нефтегазоносность и перспективы освоения. С-Пб, 2004, с. 119-132.

6. Белонин М.Д. Прищепа О.М. О стратегии воспроизводства запасов нефти и газа в Тимано-Печорской провинции и ее акваториальном продолжении // Минеральные ресурсы России. Экономика и управление. 2005. - № 1 . - С. 27-34.

7. Белонин М.Д., Буданов Г.Ф., Головань А.С. и др. Поиски залежей углеводородов на больших глубинах Тимано-Печорской нефтегазоносной провинции. Сов.геол. №3.- М.- 1990, С. 159-164.

8. Берзон И.С., Динамические характеристики сейсмических волн в реальных средах, Изд. АН СССР, 1963

9. Боганик Г.Н., Гурвич И.И. Сейсморазведка. Тверь: АИС, 2006, 744 с. Ю.Бреховских JI. М. Волны в слоистых средах, Изд. Наука, 1973, 343 с. П.Ведерников Г.В. Методика и технологии сейсморазведочных работ.

10. Новосибирск-Томск-Нортхэмтон, STT, 2006.

11. Воскресенский Ю.Н., Иноземцев А.Н., Резванов Р.А., 2005, Атрибутный AVO-анализ для разделения угольных и УВ-аномалий. Технологии сейсморазведки; 1/2005, издательство ЕАГО.

12. З.Воскресенский Ю. Н. Изучение изменений амплитуд сейсмических ^ отражений для поисков и разведки залежей углеводородов // М. РГУ нефти и газа, 2001.

13. Воскресенский Ю.Н., 2002, Состояние и перспективы развития методов анализа амплитуд сейсмических отражений для прогнозирования залежей углеводородов: Обзорная информация: М., ООО "Геоинформцентр"

14. Гогоненков Г.Н., Кириллов С.А., Ларичев В.А., Максимов Г.А. Решение обратной динамической задачи восстановления свойств тонкослоистых пластов. Тезисы Международной научно-практической конференции «Санкт-Петербург 2006», Санкт-Петербург 2006.

15. Гогоненков Г.Н. Прогнозирование геологического разреза по сейсмическим данным. Геология нефти и газа, 1980, № 1

16. Гогоненков Г.Н., Захаров Е.Т., Эльманович С.С. Восстановление детального акустического разреза по данным сейсморазведки.- В кн.: Прикладная геофизика. М., 1980, № 97

17. Гогоненков Г. Н. Изучение детального строения осадочных толщ сейсморазведкой. Москва. Недра. 1987

18. Дедеев В.А., Малышев Н.А., Аминов Л.З. и др. Научные основы развития поисково-разведочных работ на нефть и газ в Тимано-Печорской нефтегазоносной провинции, Сыктывкар.- 1987

19. Дедеев В.А., Аминов Л.З., Гецен В.Г и др Тектонические критерии прогноза нефтегазоносности Печорской плиты,- Л.: Наука.- 1986.

20. Екименко А.В. Использование динамической интерпретации сейсмических данных с целью прогноза коллекторских свойств в условиях карбонатного разреза. Записки Горного института Т. 182 С. 16.

21. Екименко А.В. Прогнозирование коллекторских свойств с использованием куба акустического импеданса. Записки Горного института Т. 183С.235

22. Ехлаков А.Ю., Белоконь Т.В., Корбух Ю.А., Звягин Г.А., Сиротенко Л.В., Коблова А.З. Новые данные по геологии нижнего девона и силура по результатам бурения Колвинской параметрической скважины // Советская геология, 1991. № 8. С.86-95.

23. Караев Н.А., Лукашин Ю.П., Прокатор О.М., Семенов В.Н., Козлов Е.А. Физическое моделирование сейсмических волн в поротрещинных средах. Тезисы Международной научно-практической конференции «Санкт -Петербург 2008», Санкт-Петербург 2008

24. Караев Н.А., Рабинович Т.Я. Рудная сейсморазведка. М.: Геоинформарк, 2000.

25. Караев Н. А. Сейсмический метод отраженных волн в рудных районах, Л. Недра, Ленинградское отделение 1982.

26. Козлов Е.А., Модели среды в разведочной сейсмологии. Тверь: ГЕРС, 2006 г.

27. Лаврик С.А., Логинов Д.В. Оценка значимости-сейсмических атрибутов на основе Многослойных нейронных сетей. Тезисы Международнойнаучно-практической конференции «Санкт-Петербург 2008», Санкт -Петербург 2008.

28. Ларионова З.В., Арасланова P.M. Новые данные о строении эйфельских отложений в Юго-Восточном Притиманье / / Изв.вузов Геол.и разв. -1988.-№ 8. с.3-14.

29. Методика поисков и разведки нефтегазоносных объектов нетрадиционного типа, отв.редактор А.Г. Алексин, М. Наука 1990

30. Мушин И.А., Бродов Л.Ю., Козлов Е.А. и др., Структурно-формационная интерпретация сейсмических данных. М.: Недра, 1990.

31. Отчёт о проведении работ по переобработке и интерпретации сейсмического материала на Нирмалинско-Антоновском объекте, Отв. исп.Мерщий Р.Ф. Нарьян-Мар, 2004, 159 с.

32. Отчёт по обобщению геолого-геофизических материалов в центральной части Колвинского мегавала, Нарьян-Мар, Мерщий Р.Ф. и др. 2004, 233 с.

33. Отчёт "Переобработка и переинтерпретация сейсморазведочных материалов по Ошской площади", Ухта, 2007.

34. Пакет геолого-геофизической информации по участку недр "Ошское месторождение". Нарьян-Мар, 2006.

35. Перродон А. Формирование и-размещение месторождений нефти и газа. Пер. с франц.- М.: Недра.- 1991.- 359 с.

36. Природные резервуары в терригенных формациях Печорского нефтегазоносного бассейна / Коми научный центр УрО Российской академии наук. Сыктывкар ,1993. -154с.

37. Прищепа О. М., Аминов JI. 3., Белонин М. Д., Богацкий В. И. и др. Нефтегазовый потенциал и геолого-экономические показатели подготовки и освоения углеводородного сырья в Тимано-Печорской провинции, СПб., 1995

38. Прищепа О. М., Аминов JI. 3., Белонин М. Д., Боровинских А. П. и др. Нефтеперспективные объекты Республики Коми (геолого-экономический анализ), Ухта, ТП НИЦ, 1999

39. Прищепа О. М. Условия формирования зон нефтегазонакопления Печоро-Колвинского авлакогена, Автореф. Дисс. канд. геол.-мин. н. ВНИГРИ, Ленинград, 1991.

40. Прищепа О. М. Стратегия сбалансированного воспроизводства запасов нефти игаза северо-западного региона России, Автореф. Дисс. докт. геол-мин. н. ВНИГРИ, СПб., 2007.

41. Прищепа О.М. Прогноз и изучение зон нефтегазонакопления в последовательном геологоразведочном процессе, Электронное издание Всероссийского нефтяного научно-исследовательского геологоразведочного института (ВНИГРИ)

42. Прищепа О.М. Зоны нефтегазонакопления методические подходы к их выделению, обеспечивающие современное решение задач отрасли, Электронное издание Всероссийского нефтяного научно-исследовательского геологоразведочного института (ВНИГРИ)

43. Прищепа О.М., Орлова JT.A Состояние сырьевой базы углеводородов и перспективы ее освоения на Северо-Западе России, Электронное издание Всероссийского нефтяного научно-исследовательского геологоразведочного института (ВНИГРИ)

44. Пузырев Н.Н. и др., Сейсморазведка методом обменных и поперечных волн, Издательство «Недра», 1985.

45. Пузырев Н.Н. Методы и объекты сейсмических исследований // Введение в общую сейсмологию. Новосибирск: Изд-во СО РАН, 1997.

46. Сейсмическая стратиграфия, под ред. Н.Я. Кунина, Г.Н. Гогоненкова, М. Мир 1982.

47. Тимано-Печорский седиментационный бассейн. Атлас геологических карт (литолого-фациальных, структурных и палеогеологических) // Никонов Н.И., Богацкий В.И., Мартынов А.В. и др. Ухта, Республика Коми, 2000.

48. Уайт Дж.Э. Возбуждение и распространение сейсмических волн. Пер. с англ. М: Недра, 1986. 262 с.

49. Урупов А.К., Левин А.Н., Определение и интерпретация скоростей в методе отраженных волн. М.: Изд. «Недра», 1985.

50. Урупов А. К. Сейсмические модели и парметры геологических сред. 1985

51. Хаттон Л., Уэрдингтон М., Мейкин Дж. Обработка сейсмических данных. Теория и практика. Пер. с англ. -М.: Мир, 1989.

52. Шерриф Р. и Гелдарт Л., Сейсморазведка: М, Изд. «Мир», 1987.

53. Adriansyah and George A. MeMechan. Effects of attenuation and scattering on AVO measurements. Geophysics, vol. 63, no. 6 (november-december 1998)

54. Avseth P., Mukerji T. and Mavko G. Quantitative Seismic Interpretation Cambridge University Press 2005

55. Bacon M., Simm R. and Redshaw Т., 3-D seismic interpretation, University press, Cambridge2003

56. Bruce S. Gibson, Alan R. Levander, Modelling and processing of scattered waves in seismic reflection surveys, Geophysics, vol.53,N0.4 (APRIL1988)

57. Bruce S. Gibson, Alan R. Levander, Apparent layering in common-midpoint stacked images of two-dimensionally heterogeneous targets, Geophysics, vol.55,NO. 11 (NOVEMBER 1990)

58. Castagna, J.P., and Backus, M.M., 1993, Offset-Dependent reflectivity -Theory and Practice of AVO Analysis: Investigation in geophysics №8

59. Castagna J.P., Bazle M.L, Eastwood, R.L, 1985, Relationships between compressional wave and shearwave velocities in clastic silicate rocks: Geophysics, 50, p. 571 -581.

60. Castagna, J.P., and Smith, S.W., 1994, A comparison of AVO indicators: A modeling study: Geophysics, 59, 1849-1855.

61. Castagna, J.P. and Swan, H.W., 1997, Principles of AVO crossplotting: The Leading Edge.

62. Castagna, J. P., Swan, H.W., and Foster D., 1998, Framework for AVO gradient and intercept interpretation: Geophysics, 63, 948-956.

63. Castagna, J. P., 1991, Seismic lithology overview: 61st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1051-1053.

64. Connolly, P., 1999, Elastic impedance: The Leading Edge, vol. 18, 438-452.

65. Florence Delprat-Jannaud, Patrick Lailly. A fundamental limitation for the reconstruction of impedance profiles from seismic data. Geophysics, vol. 70, no. 1 (january-february 2005)

66. Guillaume Cambois, Can P-wave AVO be quantitative? THE LEADING EDGE NOV 2000

67. Knapp R. W., Fresnel zones in the light of broadband data. Geophysics, vol. 56. no. 3 (march 1991): p. 354-35975.Knapp R. W

68. Marmalyevskyy N., Roganov Y., Kostyukevych A., Gazarian Z. AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack. EAGE 68th Conference & Exhibition, 2006, Vienna, Austria

69. Marmalyevskyy N., Roganov Y., Kostyukevych A., Gazarian Z. Frequency Depending AVO for a Gas-Saturated Periodical Thin-Layered Stack. 76th SEG Conference, New Orleans, USA, 2006

70. Richard H. Groshong, Jr. 3-D Structural Geology Springer Berlin Heidelberg New York 2006

71. Robert H. Tatham* and Paul 1. Stoffa. Vp/Vs-a potential hydrocarbon indicator. Geophysics, vol. 41, no. 5 (october 1976). p. 837-849

72. Ruben D. Martinez. Wave propagation effects on amplitude variation with offset measurements: A modeling study. Geophysics, vol. 58, no. 4 (april); p. 534-543

73. Rutherford, S. R., and Williams, R. H., 1989, Amplitude versus offset variations in gas sands: Geophysics, 54, 680-688.

74. Sheriff, R.E. 2002. Encyclopedic dictionary of exploration geophysics.SEG, Tulsa, 429.

75. Shree К. Nayar , Katsushi Ikeuchi , Takeo Kanade, Surface Reflection: Physical and Geometrical Perspectives, IEEE Computer Society Washington, DC, USA, 1991

76. Shuey, R. Т., 1985, A simplification of the Zoeppritz equations: Geophysics, 50, 609-614.

77. Steven R. Rutherford and Robert H. Williams. Amplitude-versus-offset variations in gas sands Geophysics, vol. 54. no. 6 p. 688.

78. M. B. Widess. How thin is a thin bed? Geophysics vol.38, (1973). p.l 176 (1973).

79. William L. Sorokan, Thomas J. Fitchz, Kirk H. Van Sickle and Philip D. North§Successful production application of 3-D amplitude variation with offset: The lessons learned, Geophysics, VOL. 67, NO. 2 (MARCH-APRIL 2002)

80. Yilmaz Ozdogan, Seismic data processing, vol. l,vol. 2, SEG 1987

81. Yinbin LiuG and Douglas R. Schmittz. Amplitude and AVO responses of a single thin bed. Geophysics, vol. 68, no. 4 (july-august 2003); p. 1161-1168