Бесплатный автореферат и диссертация по наукам о земле на тему
Прогнозирование напряженного состояния коллекторов и флюидоупоров нефтегазовых залежей
ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Автореферат диссертации по теме "Прогнозирование напряженного состояния коллекторов и флюидоупоров нефтегазовых залежей"

На правах рукописи УДК 552 082:550.832834(571 1)

ИБРАЕВ ВАЛЕРИЙ ИВАНОВИЧ

ПРОГНОЗИРОВАНИЕ НАПРЯЖЕННОГО СОСТОЯНИЯ КОЛЛЕКТОРОВ И ФЛЮИДОУПОРОВ НЕФТЕГАЗОВЫХ ЗАЛЕЖЕЙ (на примере Западно-Сибирской нефтегазоносной провинции)

Специальность 25.00.10 - «Геофизика, геофизические методы поисков полезных ископаемых»

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора технических наук

Москва - 2006

003068005

Работа выполнена в Российском Государственном Университете нефти и газа им. И М. Губкина на кафедре разведочной геофизики и компьютерных систем

Научный консультант'

доктор геолого-минералогических наук, профессор А К. Урупов Официальные оппоненты

доктор технических наук Е.А. Козлов

доктор технических наук, профессор М.Б Рапопорт

доктор геолого-минералогических наук Ф.З Хафизов

Ведущее предприятие: Тюменское отделение «СургутНИПИнефть»

Защита состоится« 20 » марта 2007 г в 15°° часов, ауд 523

на заседании диссертационного совета Д 212.200.05 при Российском Государственном Университете нефти и газа им. И.М. Губкина по адресу:

Москва, В-296 ГСП-1,119991, Ленинский проспект, д. 65.

С диссертацией можно ознакомиться в библиотеке РГУ нефти и газа им. И М Губкина

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

Автореферат разослан« 15 » декабря 2006 г

Учены» секретарь диссертационного совета, доцеш

6 - Л э

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

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

Разработка методики и технологии площадного прогноза аномальных давлений на основе данных сейсморазведки и геофизических исследований скважин (ГИС), вертикального сейсмического профилирования (ВСП) применительно к условиям Западной Сибири базировалась на исследованиях следующих авторов: А Г. Авербуха, Б.Л. Александрова, К.А. Аникиева, П.М. Бондаренко, Н.В. Вассоевича, Е.И. Гальперина, Г.Н. Гого-ненкова, А.И. Гальченко, В.М. Добрынина, A.A. Карцева, O.K. Кондратьева, П.Н. Кропоткина, O.JI. Кузнецова, Ю.А. Курьянова, Е.В. Кучерука, A.C. Мелик-Пашаева, И И. Нестерова, О.М. Нелепченко, И.К. Овчинникова, А.А Орлова, С Н. Птецова, М.Б. Рапопорта, Ф.К. Салманова, В.А. Серебрякова, А.К. Урупова, В.Д. Швецова, Л.П. Шендерей и ряда зарубежных авторов:У. Фертля, К. Магары, М. Хабберта, В. Ил-линга, Д Девиса, К. Хотгмана, Р. Чемпена и других.

Актуальность работы

Прогнозу напряженного состояния коллекторов и флюидоупоров (глинистых покрышек) в условиях Западной Сибири отводится важная роль в общем комплексе геофизических исследований Это связано с необходимостью увеличения глубинности исследований и со сложностью горно-геологических условий перспективных на нефть и газ объектов Как показывает практика бурения глубоких скважин на севере Западной Сибири в зонах аномальных пластовых (в коллекторах) и поровых (в глинах) давлений

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

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

Цель работы

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

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

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

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

3. Создание методики обобщенного профиля для оценки достоверности определения кинематических характеристик и структурных построений

4. Разработка технологии и методики геосейсмического моделирования.

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

Методика исследований включает:

- анализ и использование широко применяющихся в настоящее время методов и способов прогноза аномально высоких пластовых и поровых давлений (АВПД, АВПоД);

анализ геологических результатов математического, геосейсмического, тек-тонофизического моделирований;

прогноз напряженного состояния нефтегазоносных коллекторов и глинистых покрышек (АВПД и АВПоД) на месторождениях севера Западной Сибири, с использованием разработанных автором и под его руководством методик, программно-алгоритмических средств;

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

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

2. Установлено влияние термобарических условий на физические свойства терриген-ных отложений.

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

4. Разработана методика и технология геосейсмического моделирования

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

6 Установлена возможность использования прогнозных параметров напряженного состояния горных пород по геофизическим данным в качестве интерпретационных признаков для выделения углеводородосодержащих объектов и геотектонически -активных зон.

7 Показана перспективность использования разработок как основы для создания региональной базы напряженного состояния осадочного чехла Западной Сибири

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

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

лекторов) и поровых (для глинистых покрышек) давлений до стадии бурения скважин

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

3. Методика и технология площадного прогноза оценки свойств флюидоупоров (глинистых покрышек) нефтегазоносных коллекторов

4 Методика и технология геосейсмического моделирования

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

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

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

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

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

4. Прогнозное распределение АВПоД по площади для глинистых покрышек до стадии бурения позволяет оценить сохранность углеводородов в выявленных ловушках и одновременно может являться поисковым признаком.

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

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

Реализация разработки и её апробация

Результаты, полученные при выполнении настоящей работы, апробированы на материалах сейсморазведочных съемок и скважинных исследований по Бованенков-ской группе месторождений, Самбургскому, Ем-Еговскому, Маслиховскому и Ай-Пимскому месторождениям. Отдельные элементы разработки применены при оперативном прогнозе по Гыданскому району на площадях: Лунная, Минховская, Тото-Яхинская, Приречная, Глубинная, Сопочная, Овражная Разработанная методика и технология прогноза напряженных свойств флюидоупоров (глинистых покрышек) апробирована на Пулытьинской-2, Кандырской площадях и Семивидовско-Толумском месторождении Оценка достоверности определения интервальных скоростей и структурных построений выполнена по Бованенковской группе месторождений (семь месторождений), Ново-Надымском, Яунлорском, Восточно-Сургутском, Чапровском, Рогожников-ском месторождениях.

Основные результаты исследований докладывались на:

семинарах и конференциях, проводимых ЦГЭ, посвященных проблемам прогнозирования геологического разреза (ПГР),

всесоюзной научно - практической конференции в г. Ленинграде, 1990г; международной геофизической конференции, организованной европейским геофизическим сообществом^ Ев), БЕв / Москва - 92;

международной геофизической конференции, организованной европейским геофизическим сообществом (БЕО) и евроазиатским геофизическим обществом (ЕАГО), БЕС - ЕАГО/Москва-93;

конференциях «Геомодель-2002», «Геомодель-2003», «Геомодель-2004»; ежегодной научно-практической конференции «Пути реализации нефтегазового потенциала ХМАО» в 2002,2003,2005 гг. в г. Ханты-Мансийске.

Личный вклад

Основу диссертации составили систематизация и обобщение многолетних разработок автора, а также результаты:

научно - исследовательских, тематических, опытно-методических работ выполненных автором в институте ФГУП ЗапСибНИИГГ (ЗапСибНИИГеофизика), связанные с разработкой технологии прогноза напряженного состояния терри-генных отложений в условиях севера Западной Сибири;

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

рамках деятельности ОАО «Тюменнефтегеофизика» (трест «Тюменнефтегеофи-зика», управление «ЗапСибнефтегеофизика»).

Исследования автора проводились в период 1984-2006 гг. Завершение работы выполнено в период прикрепления автора соискателем ученой степени доктора технических наук приказом № 7-А от 17 февраля 2006 г. по Российскому государственному университету нефти и газа им. И.М Губкина. Публикации

Основные положения диссертации представлены более чем в 40 научно-исследовательских, тематических работах и 19 публикациях, в том числе 12 работ опубликованы в монографии и статьях в ведущих рецензируемых научных журналах и изданиях согласно перечню ВАК Объем работы

Работа состоит из введения, 6 глав и заключения. Содержит 180 страниц машинописных листов, 73 рисунка, 11 таблиц, одно приложение, библиография включает 96 наименований Благодарности

Автор выражает искреннюю признательность за сотрудничество, консультации и ценные советы при работе над методиками и технологиями моим коллегам: д.г. - м.н. Ахиярову В.Х., к т.н. Анпенову С В , к.г. - м.н. Белкину Н.М. д г.-м.н. профессору Бондаренко П.М., д ф. - м н. Голошубину Г.М., к г - м н Гальченко А И , д.т н Ждановичу В В., к ф. - м.н Кравченко Г.Г., д.т.н. Курьянову Ю.А., к.г.-м.н Теплоухо-ву В.М , к г.-м н Хабарову В.В

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

Автор выражает признательность коллегам по работе Ковалеву В.Г., Кондратенко М.А., Коровкину М А., Ширшову А С., Кузнецовой Л.Г , Романову Д.В

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

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

ГЛАВА 1. Теоретические и физические предпосылки прогноза напряженного состояния коллекторов и флюидоупоров в геологических условиях севера Западной Сибири.

Исследован механизм формирования геотектонических напряжений в сплошных и пористых средах. Дана физическая характеристика аномально высоким пластовым и поровым давлениям коллекторов и флюидоупоров (глинистых покрышек) Рассмотрены представления о генезисе аномально высоких пластовых давлений. Выполнено обобщение исследований обзорного характера по рассматриваемой тематике Охарактеризована геологическая обстановка разрезов рассматриваемой территории Представлено районирование территории по месторождениям с АВПД и типам разрезов. Определены направления исследований

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

Доминирующее положение в арсенале методов оценки и прогнозирования АВПД и АВПоД занимают методики, основанные на изучении закономерностей изменения свойств пород в результате их уплотнения (в процессе осадконакопления). Для прогноза АВПД и АВПоД по данным геофизических исследований скважин (В.М.Добрынин, В. А. Серебряков, Б Л Александров, У.Фертль и другие) в качестве критерия оценки используются свойства глин, зависящие от степени их уплотнения при аномальных поровых давлениях.

Необходимость комплексного решения задачи по выделению зон аномальных пластовых и поровых давлений, прогнозирования давлений отмечается многими исследователями. Доказана возможность и целесообразность проведения этих работ на основе сейсмических исследований до постановки бурения. Примеры использования сейсмических данных приводятся У. Фертлем, коллективом исследователей ГАНГ (РГУ нефти и газа им И.М Губкина) В.М. Добрыниным, М Б. Рапопортом, В А. Серебряковым, Г.А. Карапетовым, А.Л. Малкиным, В.Г. Мартыновым, В.Ю. Матусеевичем, с использованием способа, основанного на расчете спектров скоростей.

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

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

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

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

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

Объектом исследований диссертационной работы является платформенный чехол Западно-Сибирской плиты Геология рассматриваемой территории представлена в работах В С Бочкарева, Ю Н. Карогодина, А.Э Конторовича, С В. Крылова, Н Я. Куни-на, И И. Нестерова, Н.Н Ростовцева, Ф К Салманова, A.A. Трофимука, Н.А Туезо-вой, Ф.З. Хафизова и многих других исследователей

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

необходимость комплексной интерпретации данных сейсморазведки и ГИС с целью прогноза аномальных давлений в терригенных отложениях до стадии бурения. ВЫВОДЫ К ГЛАВЕ

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

ТИПЫ РАЗРЕЗОВ:

1. Ямальский

2. Надымский

3. Уренгойский

4. Южно-Русский

5. Пурпейский

6. Приуральский

7. Сургутский

8. Нижневартовский

9. Александровский

Рис.1. Районирование территории по типу разреза.

флюидоупоров (глинистых покрышек) нефтегазоносных залежей в Западно -Сибирской нефтегазоносной провинции

2 Рассмотрено формирование геотектонических напряжений в сплошных и пористых средах.

3 Дана физическая характеристика аномально высоким пластовым и поровым давлениям Рассмотрены представления о генезисе аномально высоких пластовых давлений

4 Выполнено обобщение исследований обзорного характера по рассматриваемой тематике

5 Дана краткая характеристика геологических условий Западной Сибири

6 Выполнено районирование территории Западной Сибири по месторождениям с АВПД и типам разрезов изменения пористости глин с глубиной.

ГЛАВА 2. Методика и технология площадного прогноза напряженного состояния терригенных отложений

Во второй главе рассмотрена деформационная модель формирования аномальных давлений. Исследованы влияющие факторы на кинематические и динамические характеристики упругих волн Представлены аналитические распределения скорости, плотности и пористости с глубиной для глин и песчаников в условиях Западной Сибири Проведен анализ результатов лабораторных исследований кернового материала для определения влияния термобарических условий на скоростные свойства пород, необходимый для решения задач площадного прогноза напряженных свойств (аномальных давлений) коллекторов и флюидоупоров (глинистых покрышек) Исследованы методы оценки аномальных давлений и способы определения напряженных свойств горных пород на основе данных ГИС, предложенных Добрыниным В.М. и Лимбергером Ю.А. (методика кривых нормального уплотнения, 1971г.), Александровым Б.Л (методика эквивалентных глубин, 1973 г.), Гальченко А И. (методика коэффициентов увеличения давлений, 1977г.), Добрыниным В М , Серебряковым В А. и Сребродольским А Д (методика компрессионной кривой, 1982г.). Методика эквивалентных глубин - графический способ, методики линии нормального уплотнения, коэффициентов увеличения давления, компрессионной кривой - аналитические способы.

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

12

напряжений реконструируемой поверхности на определенное геологическое время. Реализация расчетов производится с использованием сеточной функции F(x,y) - карта выбранной поверхности по формуле: Нг = Hi-(H-(Fmm+K)), где Нг и Hi- соответственно рассчитанная и исходная глубины реконструируемой поверхности; Н - глубина поверхности приведения; Fmin - минимум функции поверхности приведения; К - поправочный коэффициент, определяемый по результатам стратиграфической привязки

Для рассматриваемого типа задач реконструкция производится по поверхности фундамента и в расчетах используется поверхность F(x,y), определенная в узлах регулярной прямоугольной сети, на то геологическое время с которого восстановленный структурный план носит унаследованный характер

Карта распределения зон максимальных напряжений представляет суммарную поверхность, полученную в результате дифференцирования по направлениям X и У с заданным шагом и по секторам в а°. В окончательном виде имеем функцию :

п

F = XFak/n k=l

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

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

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

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

- обоснование оптимальной корреляционной зависимости интервальной скорости от коэффициента аномальности пластовых или поровых давлений.

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

дольных волн в нормально уплотняющихся глинах (без аномально высоких поровых давлений - АВПоД) Следовательно, необходимым условием является наличие зависимости скорости в нормально уплотненных глинах (Унг) от глубины Зависимость Ут=Р(Н) может быть получена двумя способами: по данным акустического каротажа (АК) в глинистых интервалах без АВПоД и аналитическим способом, включая интервалы с аномальными поровыми давлениями, по данным площадной сейсмической съемки

Первый способ предполагает наличие глинистых отложений без АВПоД во всем диапазоне глубин изучаемого разреза На практике это условие не всегда выполняется В частности на Бованенковской площади на глубине более 2300 м такие отложения отсутствуют, что не позволяет охарактеризовать основные аномально проявляющиеся отложения в интервалах реперных горизонтов М, Во, Б и глубже. Поэтому автором реализован второй подход, в котором решается прямая задача нахождения УИНт в нормально уплотняющихся глинах по известным его величинам; коэффициентам аномальности поровых давлений и с учетом типа насыщения песчаников

Изложенная методика реализована в виде автоматизированной технологии прогноза напряженного состояния терригенных отложений для условий Западной Сибири. Автоматизированная технология состоит из отдельных программных средств и процедур обработки данных сейсморазведки и ГИС, объединенных в единый пакет прикладных программ

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

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

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

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

-расчет, коэффициента аномальности и порового давления по методике эквивалентных глубин,

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

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

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

-выделения однородных по геофизическим характеристикам интервалов и расчета средних значений параметров;

-определения индекса литологии и насыщения;

-вычисления удельного электрического сопротивления глинистых интервалов; -построения линии нормального уплотнения;

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

Результаты площадного прогноза напряженных свойств коллекторов и флюидо-упоров (глинистых покрышек) в базе СОИ позволяет объединять данные с целью создания региональной базы напряженного состояния платформенного чехла Западной Сибири

ВЫВОДЫ К ГЛАВЕ

1. Рассмотрена деформационная модель формирования давлений

2. Исследовано влияние термобарических условий на физические свойства тер-ригенных отложений

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

4. Теоретический обоснована и разработана технология прогноза напряженных свойств коллекторов и флюидоупоров (АВПД, АВПоД) нефтегазоносных залежей до стадии бурения, на основе методов компрессионной кривой и эквивалентных глубин.

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

6. Разработана система обмена информации.

ГЛАВА 3. Результаты применения методики и технологии прогноза напряженного состояния коллекторов и флюндоуноров нефтегазоносных залежей севера Западной Сибири

Рассмотрено применение технологии прогноза напряженных свойств коллекторов и флюидоупоров нефтегазоносных залежей Западной Сибири.

Анализ материалов региональных работ, результатов районирования Западной Сибири по изменению пористости терригенных отложений с глубиной, региональная прослеживаемость опорных отражающих горизонтов, относительная выдержанность интервалов по мощности в пределах площадей сейсмических съемок свидетельствуют о благоприятных условиях использования данных сейсморазведки для прогноза АВПД и АВПоД

На первом этапе опробования технологии по Бованенковской группе месторождений (7 месторождений) построен обобщенный вертикальный профиль распределения интервальных скоростей от среднеинтервальных глубин по площади в интервалах опорных горизонтов' Г, Гг, М, М, В0, Б, Т^, А

Оценка достоверности и соответствия распределения Уинт. общей закономерности изменения скоростей по ВСП производится на основе скорректированных трафиков скоростей. Полученные интервальные скорости приводятся к условиям глинистого разреза и используются для расчета коэффициента аномальности по формуле:

Р V -V

Ка= 1+ ( —^ -1)(—--),

4 Р лУ -V

ГНИ * НГ

где Рг и Рщд-горное и гидростатическое давления соответственно; Унт - интервальная скорость нормально уплотненных глин, Угл - интервальная скорость, приведенная к условиям глинистого разреза; V - скорость в точке пересечения линии нормального уплотнения с осью ОХ т.е. при Н=0 является постоянной величиной для площади исследований. Для Бованенковской площади скорость приведения (при Н=0) равна 1450 м/с.

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

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

В технологии реализовано два варианта прогноза: с использованием единой зависимости нормального уплотнения глин и поинтервальной оценкой положения "линии глин". Имеющиеся статистические данные сопоставления коэффициентов аномальности, определенных по данным ГИС и рассчитанных с использованием программных средств пакета не позволяют однозначно отдать предпочтение первому или второму варианту.

Результаты интерпретации материалов ГИС достаточно апробированы и многократно подтверждены испытанием скважин, а хорошая сходимость сопоставления расчетных значений по ГИС и сейсморазведке является доказательством высокой эффективности площадного прогноза. Из общего количества прогнозных значений 92 % удовлетворяют требованиям точности оценки технологических параметров бурения в зонах АВПД (рис. 26) Дополнительная оценка коэффициента аномальности по данным вертикального сейсмического профилирования и результатам интерпретации ГИС, как более детальная, характеризует достоверность выполненных расчетов на основе кинематической интерпретации.

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

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

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

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

Принятое условно превышение пластового или норового давления в 1,3 раза над гидростатическим давлением считается аномальным, а сама величииа носит название коэффициента аномальности (Ка). На АЛ - Пимеком месторождении граничное значение Ка принято равной 1,2.

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

Рис. 2, Площадной прогноз Ка в интервалах разреза Бованепконской группы месторождений.

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

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

Возможности разработанной технологии площадного прогноза напряженного состояния терригенных отложений с использованием кинематических и динамических параметров рассмотрены применительно к Пойменной площади Самбургского месторождения в интервалах усть - тазовской серии и баженовской свиты. Коэффициент аномальности в исследуемых интервалах "М-В2","В2-В4" характеризуется относительно равномерным распределением по площади. Повышенные значения порового давления отмечаются в интервале «В4-Б» с глубины 3500 метров, а по площади в этом интервале градиентное возрастание коэффициента аномальности (Ка) имеет северную и северозападную направленность. Расчеты Ка с привлечением динамических характеристик выполнены по профилю 19 на пикете 126 5, где вместо Уинт использованы акустические жесткости. Оценка достоверности прогноза проведена сопоставлением расчетных значений, полученных на основе динамических характеристик и определений Ка по ГИС, и подтверждена результатами испытаний в скважине 700 Самбургского месторождения. Достоверность прогноза Ка высокая. Погрешность площадного прогноза АВПД и АВПоД в среднем составляет ± 3 МПа.

Практическое применение элементов разработанной технологии проведено при оперативном прогнозе АВПД на месторождениях: Лунное, Минховское, Тотояхинское, Приречное, Глубинное, Сопочное, Овражное

ВЫВОДЫ К ГЛАВЕ

1. Рассмотрена оценка напряженных свойств поверхностей сложенных коренными горными породами на основе обработки результатов площадной палеоре-конструкции Бованенковской группы месторождений

2. Представлены результаты исследований по прогнозу АВПД и АВПоД на 17 месторождениях Западной Сибири, в том числе по десяти площадям проведен площадной прогноз, по семи площадям оперативный прогноз.

ГЛАВА 4. Методика обобщенного вертикального профиля интервальных скоростей.

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

19

разреза исследуемой территории по данным вертикального сейсмического профилирования, сейсмического каротажа (СК) или акустического каротажа (АК) и плотностные свойства по данным гамма-гамма каротажа плотностного ГГК (П).

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

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

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

Предлагаемая автором методика позволяет решить одновременно две задачи.

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

Вторая задача, решаемая предлагаемой методикой это оценка точности структурных построений.

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

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

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

Количество точек для одного интервала для сейсмической съемки модификации 2Д может достигать 20000 и больше, а для модификации ЗД более 500000. Поэтому допускается равномерное прореживание, т.е. использовать не все точки сетки, а каждую 5, 10 или 100. Полученные данные по всем интервалам площади сводятся в один массив Уинт в другой массив Ней Сводные массивы переносятся в пакет программ для построения вертикального профиля.

Обобщенный вертикальный профиль интервальных скоростей был построен по Бованенковской группе месторождений. Количество профилей на площади 160, интервалов 7, точек на обобщенном вертикальном профиле больше 5500. Часть точек при равных величинах интервальной скорости и средней интервальной глубины накладываются, что в целом не влияет на получаемое распределение. Интервалы Гг-М1 и М-Во имеют отклонение от общей закономерности изменения интервальной скорости с глубиной, рассчитанной по данным ВСП. После уточнения корреляции определены величины поправок которые составили для горизонта М1 +128м, а для горизонта М +110 м. Результаты повторного построения после редакции свидетельствуют о корректности определения скоростей по интервалам, которые могут быть использованы для дальнейших расчетов для прогноза напряженных свойств коллекторов и флюидоупоров (глинистых покрышек). Для сравнительной оценки распределения интервальных скоростей на обобщенном вертикальном профиле и рассчитанных по материалам ВСП (СК, АК) необходимо представлять совместно.

Для решения второй задачи оценки точности структурных построений по предлагаемой методике были выполнены построения по Ново-Надымской, Яунлорской, Восточно-Сургутской, Чапровской, Рогожниковской площадям.

ВЫВОДЫ К ГЛАВЕ

1. Обоснована и реализована возможность оценки достоверности определения ин-

тервальных скоростей по методике обобщенного профиля.

2. Установлена возможность уточнения и коррекции привязки сейсмических отра-

жений к геологическим объектам.

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

строений на основе методики обобщенного профиля по 6 месторождениям.

4. Методика обобщенного профиля является экспресс методом уточнения интер-

вальных скоростей и оценки достоверности структурных построений.

ГЛАВА 5. Методика площадной оценки качества флюидоупора (глинистой покрышки) нефтегазоносных коллекторов.

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

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

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

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

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

22

Рпл. - пластовое давление - нагрузка, испытываемая поровым пространством пласта (толщи) коллектора.

Связь скорости продольной волны и порового давления определяются формулой:

у2 ГРг РпорМо^ ^рп рп ) Кп

где Рг - горное давление (кГ/м2), Рпор- поровое давление(кГ/м2), g - ускорение свободного падения (м/с2), 103 - множитель согласования единиц измерения, Кп- пористость (в долях единиц 0<Кп<1), рп- плотность((кГ/м3)

Из формулы связи скорости продольной волны и порового давления следует:

рпор^р^^рп

Ю3ё

Нагрузка на пласт (толщу), определяемая геостатическим давлением, распределяется

между скелетом породы и поровым пространством т е ■

Рг = Рэфф. + Рпор и соответственно Рпор. = Рг - Рэфф

Подставляя правые части равенств в формулу Рпор., получим: Рэфф.= ^

1036

Геостатическое давление для средней интервальной глубины определяется как' Рг = Ней рп,

где Нси-глубина до средины интервала в текущей точке,

рп - плотность породы, рассчитывается по ГГК (П) или принимается для терригенных отложений 2.31 г/см3, градиент геостатического давления при этом составляет 0,0231 МПа/м.

Предложенная автором методика опробована при выполнении интерпретационных работ и площадного прогноза поровых давлений глинистых покрышек на Пулытьин-ской -2, Кандырской площадях и Семивидовско-Толумском месторождении

ВЫВОДЫ К ГЛАВЕ

1. Установлена взаимосвязь между поровым давлением и интервальной скоростью, эта взаимосвязь является основой предложенной автором методики.

2. Представлено описание методики для практической реализации.

3. Приведены примеры использования методики на Пулытьинской-2, Кандырской площадях и Семивидовско-Толумском месторождении

ГЛАВА 6. Технология и методика геосейсмического моделирования.

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

1. Первичная обработка результатов измерения поля,

2. Трансформация экспериментальных значений поля;

3. Решение обратной задачи;

4. Интерпретация полученных результатов.

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

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

При методическом руководстве и непосредственном участии автора совместно с Томским отделением ЗапСибНИИГеофизика разработана технология математического

моделирования в виде автономного пакета программ Технологическая реализация пакета программ геосейсмического моделирования предусматривает:

- формирование, коррекцию, хранение, копирование моделей;

- формирование и коррекцию инициирующего импульса,

- задание геометрии наблюдения;

- осреднение каротажных данных,

- расчет синтетической трассы, временных разрезов и сейсмограмм;

-графическое представление исходных и конечных данных на экране дисплея

или на принтере;

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

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

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

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

дополнительным сервисом, копировать всю базу моделей или отдельную модель, визуализировать.

Блок подготовки и коррекции импульса является независимым модулем и позволяет формировать различные типы инициирующих импульсов (Берлаге, Риккера, Пузыре-ва и Синусоидальный) с изменением параметров в диалоговом режиме и визуальный просмотр графика импульса

Математический аппарат позволяет методом подбора и последовательных приближений уточнять параметры объектов содержащих углеводороды и объектов, подверженных литологической изменчивости по латерали

Основным преимуществом рассмотренной технологии, по сравнению с существовавшими раньше, является.

- диалоговый режим работы;

- отсутствие ограничений по толщине слоев,

- оперативный доступ к информации;

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

Рассмотрены примеры подбора модели среды при решении задач по определению влияния на формирование волнового поля:

- типа насыщения и литологического замещения (глинизации);

- выклинивания тонкослоистого пласта коллектора;

- изменение физических свойств пород коры выветривания (трещинные коллектора) и свойств вмещающих пород.

Как уже отмечалось, одномерное моделирование позволяет решать ряд сложных задач, но не дает полной картины изменения волнового поля по профилю. А это, в свою очередь, осложняет стратиграфическую корреляцию ОВ, связанных с широким разнообразием тонкослоистых геологических объектов, клиноформных напластований и границами их выклинивания. Для решения обозначенных задач целесообразно привлечение двумерного моделирования, динамического анализа и сейсмофациального анализа Рассмотрена практическая реализация технологии двумерного моделирования при решении задач - выделения границ фациального замещения коллекторов вогулкинской толщи абалакской свиты в приуральской зоне Западной Сибири на одном из месторождений Шаимского нефтегазоносного района с привлечением отмеченных выше видов анализа В результате выполнено районирование территории по условиям осадконакопления и по просьбе Заказчика автором был представлен прогноз на наличие коллектора вогулкинской толщи абалакской свиты в разрезе проектной скважины 10006П Семивидовско-Толумского месторождения (рис. 3) Прогнозная

эффективная толщина составляла 3-5 метров, а вскрытая бурением 3,8 метра. Результаты прогноза полностью подтвердились В скважине из песчаников вогулкинской толщи абалакскои спиты, получен фонтанный приток нефти.

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

Двумерное геосенсмическоё моделирование для установления влияния латеральной изменчивости физических свойств пород коры выветривания и свойств вмещающих пород на формирование волнового поля рассмотрен на примере Даниловского месторождения. Результат геосейсмического моделирования свидетельствуют о доминирующем влиянии на формирование отраженной воины «А» изменение физических свойств подстилающих толщ (слоев). Влияние изменения физических свойств покры-

Рие.З. Точка заложения проектной скважины а - фрагмент сейсмического разреза; б - геологическая модель; в - геосеисмичсскам модель.

вающей толщи менее значимо. Результаты моделирования доказывают возможность картирования коры выветривания по площади, а широкий диапазон изменения плотности пород коры выветривания (от 2.13 г/см3 до 2.52 г/см3) предполагает взаимосвязи динамических характеристик волнового поля с геологическими параметрами.

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

Кроме рассмотренных выше средств моделирования для оценки влияния флюидонасыщения пласта коллектора на волновое поле, исходя из основ теории поля, такую оценку можно провести на основе сопоставления волновых полей разных типов волн. Продольная волна характеризуется распространением в среде упругой деформации (сжатия и растяжения), а поперечная волна распространением сдвиговой деформации. Под воздействием горного давления флюид в поровом пространстве породы находится под давлением, отличающимся от напряжения в скелете породы. При замкнутой системе, препятствующей оттоку флюида из порового пространства, деформация осадочных пород сопровождается, полной передачей нагрузки на флюид и в этом случае давление противодействия боковому распору равно поровому g = Рпор, а через горное давление напряжение противодействия боковому распору g = Рг 6/(1-5), где Рг - горное (геостатическое) давление;

5 - коэффициент Пуассона Приравнивая правые части уравнений и преобразуя, получим.

8=Рпор/(Рпор+ Рг)

Из этой формулы следует, что коэффициент Пуассона является деформационной характеристикой осадочных пород, которая рассматривается как система " скелет-флюид", а при Рпор = Рг достигает значения 0.5.

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

Динамические условия равновесия после замены в них напряжений деформацией переходят в уравнения упругих волн распространения деформаций в среде. В жидкостях и газах наблюдается упругость объема, а не формы, поэтому возникают только продольные волны Коэффициент Пуассона через соотношение скорости поперечной волны Vs к скорости продольной волны Vp определяется соотношением:

5= (1-2(У8/Ур)2) / 2(1-(Уз/Ур)2).

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

Рассмотрено влияние флюидонасыщения пластов БВ7 и БВв на волновое поле продольных (РР) и обменных (РБ) волн, по результатам не продольного вертикального сейсмического профилирования в скважине 399.

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

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

Основы тектонофизического моделирования были заложены Гзовским В.М. и продолжены Бондаренко П.М. в ИГиГ СО РАН, г. Новосибирск. Неоднократные встречи автора с Бондаренко П М. по вопросам применимости элементов тектонофизического моделирования при интерпретации материалов сейсморазведки позволили выработать подходы для структурно-динамической интерпретации с позиций тектонофизического моделирования.

Автором представлены результаты применения методики на материалах сейсморазведки на Ем-Еговской и Ай-Пимской площадях.

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

часть разреза в зону дифференциальных движений на времени 1.3 -1.7 с. Особый интерес в нефтегазоносном отношении представляет прифундаментная часть разреза, так как именно в этом интервале сосредоточены продуктивные объекты. В интервале времени 1.85-2.2 с имеет место множество тектонических нарушений, вызванных блоковым строением фундамента В центре антиклинали по оси поднятия направленность тектонических нарушений носит субвертикальный характер, в краевых частях разломы ориентированы больше в южном направлении. Напряжения, вызванные блоковым поднятием фундамента, распространяются, как вертикально, так и по латерали, образуя зоны деструкций и дислокаций Последние могут контролировать литологические ловушки и являются основополагающими факторами при формировании трещиноватых коллекторов в коре выветривания. Насыщением этих коллекторов и объясняется повышенное поглощение в интервале горизонтов "К-А"

ВЫВОДЫ К ГЛАВЕ

1. Разработана технология геосейсмического моделирования.

2. Представлены результаты опробования технологии для:

стратиграфической привязки,

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

выделения границ фациального замещения коллекторов вогулкинской толщи абалакской свиты,

создания модели выклинивания тонкослоистого пласта коллектора; моделирования волнового поля по типу насыщения и литологического замещения (глинизации);

изменение физических свойств пород коры выветривания (трещинные коллектора) и свойств вмещающих пород.

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

ЗАКЛЮЧЕНИЕ

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

- эффективной автоматизированной технологии площадного прогноза аномальных давлений;

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

- методики площадного прогноза АВПоД для оценки свойств флюидоупоров (глинистых покрышек) нефтегазоносных залежей;

- технологии геосейсмического моделирования;

- системы обмена и хранения геолого-геофизической информации.

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

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

Показана возможность повышения точности прогноза поровых давлений и обеспечения их оценок с погрешностью ±ЗМПа при решении задачи площадного прогноза АВПД на основе сейсмических данных. С привлечением результатов интерпретации ГИС на стадии поисково-разведочных работ может быть достигнута меньшая погрешность.

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

тверждаемость прогноза напряженного состояния коллекторов и флюидоупоров составляет 92%.

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

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Ибраев В И. Обработка и интерпретация материалов пространственной сейсморазведки программными средствами СЦС-3 ПГР и ИНТЕРСЕЙС. Геология нефти и газа, № 7, М., 1990, (совместно со Ждановичем В.В., Щекиным С.Н.).

2. Ибраев В.И. Формирование волновых полей реальных геологических сред методами математического моделирования. - Сборник научных трудов ЗапСибНИГНИ, Тюмень, 1989. С. 84-89, (совместно с Пустовитом В.И.).

3. Ибраев В.И. Прогнозирование АВПД по данным ГИС - сейсморазведка в геологических условиях Западной Сибири,- Сборник ВНИГРИ: АВПД и нефтеносность, г. Ленинград, 1992 (Совместно с Анпеновым С.В , Гальченко А.И., Теплоуховым В М.).

4. Ибраев В.И. Прогноз аномальных давлений в терригенных отложениях Западной Сибири на основе комплексирования геофизических методов. - Геология нефти и газа, № 11, М.,1993.

5. Ибраев В И. Монография «Интерпретация материалов объемной сейсморазведки». - Разведочная геофизика, обзорная информация, выпуск 9 М., 1991, (совместно со Ждановичем В.В., Потаповым O.A., Кузнецовым В И.) - 64 с.

6. Ибраев В.И. Новые возможности изучения скоростных разрезов в сейсморазведке. - SEG/Москва, 1992, (совместно с Епинатьевой A.M.,

Соловьевым B.C., Голошубиным Г.М., Рябиновой Л.Г.).

7. Ибраев В.И. Прогноз напряженного состояния терригенных отложений на основе комплексной интерпретации материалов ГИС и сейсморазведки - SEG EAGO/Москва, 93, (совместно с Анпеновым С.В , Гальченко А.И., Теплоуховым В М.).

8. Ибраев В.И Геосейсмическое моделирование при решении задач выделения тонкослоистых коллекторов// Геофизика «Технологии сейсморазведки- 2». Специальный выпуск - М.: ЕАГО. - 2003. С. 22-26, (совместно с. Белкиным Н.М)

9. Ибраев В.И. Методическое руководство по реализации прогноза напряженного состояния терригенных отложений (аномально высоких пластовых и поровых давлений АВПД, АВПоД) на основе комплексной интерпретации сейсмических данных и материалов геофизических исследований скважин Тюмень. 2002, (совместно с Белкиным Н.М. Стариковой Т.И.), 24с.

10. Ибраев В.И. Прогноз напряженного состояния верхнеюрских отложений Сургутского свода по комплексу ГИС - сейсморазведка. Геология нефти и газа 2002/5. (совместно с Белкиным Н.М., Косом И.М.).- С. 42-45.

11. Ибраев В.И. Частотно-резонансный анализ сейсмических данных. В кн. Пути реализации нефтегазового потенциала ХМАО (5-ая научно-практическая конференция) Ханты-Мансийск. 2002 Том 1. С. 310-313, (совместно с. Белкиным Н.М.). 12 Ибраев В И Влияние литологической изменчивости и типа флюидонасыгцения пласта - коллектора на волновое поле В кн. Пути реализации нефтегазового потенциала ХМАО (8-ая научно-практическая конференция) Ханты-Мансийск. 2005. С. 77-86, (совместно с Белкиным Н М.).

13. Краткий обзор использования данных сейсморазведки, ВСП (СК), ГИС для площадного прогноза напряженного состояния терригенных отложений до стадии бурения// Геофизика «Технологии сейсморазведки». Специальный выпуск - М.: ЕАГО -2004 - С. 22-26, (совместно с. Белкиным Н.М.).

14. Ибраев В И Разработка технологии прогноза напряженного состояния терригенных отложений на основе комплексной интерпретации данных сейсморазведки и ГИС в условиях Западной Сибири Диссертация на соискание ученой степени кандидата технических наук

15. Ибраев В И. Динамический анализ и моделирование для выделения границ фа-циального замещения коллекторов абалакской свиты. Геология нефти и газа, № 4, М., 2006,-С 32-38.

16 Ибраев В.И Методика площадной оценки качества флюидоупора (глинистой покрышки). Геофизика № 3, М , 2006. -С. 27-29

17. Ибраев В.И. Влияние термобарических условий на физические свойства терригенных отложений. - Известия высших учебных заведений «Нефть и газ», № 3, 2006. С. 4-10

18. Ибраев В И Методика структурно-динамической интерпретации временных сейсмических разрезов с позиций тектонофизического моделирования. Геофизика № 4, М, 2006.

19. Ибраев В.И Методика обобщенного вертикального профиля интервальных скоростей. - Известия высших учебных заведений «Нефть и газ», № 4,2006.- С. 6-10.

Подписано к печати 27.06.2006 г Объем 1,8 п.л. Тираж 120 экз Заказ № 73

ОАО «Тюменнефтегеофизика» 625023, Тюмень, ул. Республики, 173

Содержание диссертации, доктора технических наук, Ибраев, Валерий Иванович

Введение.

Глава 1. Теоретические и физические предпосылки прогноза напряженного состояния коллекторов и флюидоупоров в геологических условиях севера Западной Сибири.

1.1. Механизм формирования геотектонических напряжений в сплошных средах.

1.2. Формирование напряжений пористых сред.

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

1.4. Представления о генезисе аномальных пластовых давлений.

1.5. Изученность проблемы прогноза напряжённого состояния коллекторов (АВПД) и флюидоупоров (АВПоД) нефтегазовых залежей с привлечением данных сейсморазведки.

1.6. Геологические условия Западной Сибири.

1.7. Районирование территории Западной Сибири по месторждениям с АВПД и типам разрезов изменения пористости глин с глубиной.

1.8. Выводы к главе.

Глава 2. Методика и технология площадного прогноза напряженного состояния терригенных отложений.

2.1. Деформационная модель формирования давлений.

2.2. Физические свойства терригенных отложений.

2.3. Методы оценки аномальных давлений и способы определения напряженного состояния горных пород на основе данных ГИС.

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

2.5. Система обмена информации (СОИ).

2.6. Программы обработки сейсмических данных.

2.7. Программы обработки данных ГИС.

2.8. Выводы к главе.

Глава 3. Результаты применения методики и технологии прогноза напряженного состояния коллекторов и флюидоупоров нефтегазовых залежей севера Западной Сибири.

3.1. Бованенковская группа месторождений.

3.2. Пойменная площадь.

3.3. Гыданский район.

3.4. Утлорская площадь, Маслиховское месторождение.

3.5. Ай-Пимское месторождение.

3.6. Выводы к главе.

Глава 4. Методика обобщенного вертикального профиля интервальных скоростей.

4.1. Актуальность проблемы.

4.2. Принцип построения обобщенного профиля.

4.3. Практическая реализация.

4.4. Выводы к главе.

Глава 5. Методика площадной оценки качества флюидоупора (глинистой покрышки) нефтегазовых коллекторов.

5.1. Постановка задачи.

5.2. Методы и способы решения.

5.3. Результаты опробования авторской методики.

5.4. Выводы к главе.

Глава 6. Технология и методика геосейсмического моделирования.

6.1. Разработка технологии геосейсмического моделирования.

6.2. Методика геосейсмического моделирования и практические решаемые задачи.

6.2.1. Подбор модели среды.

6.2.2. Моделирование и динамический анализ при решении задач выделения границ фациального замещения коллекторов абалакской свиты.

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

6.3. Методика структурно-динамической интерпретации с позиций тектонофизического моделирования

6.4. Выводы к главе.

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

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

Разработка методики и технологии площадного прогноза аномальных давлений на основе данных сейсморазведки и геофизических исследований скважин (ГИС), вертикального сейсмического профилирования (ВСП) применительно к условиям Западной Сибири базировалась на исследованиях следующих авторов: А.Г. Авербуха, Б.Л. Александрова, К.А. Аникиева, П.М. Бондарешсо, Н.В. Вассоевича, Е.И. Гальперина, Г.Н. Гогоненкова, А.И. Гальченко, В.М. Добрынина, А.А. Карцева, O.K. Кондратьева, П.Н. Кропоткина, O.J1. Кузнецова, Ю.А. Курьянова, Е.В. Кучерука, А.С. Мелик-Пашаева, И.И. Нестерова, О.М. Нелепченко, И.К. Овчинникова, А.А. Орлова, С.Н. Птецова, М.Б. Рапопорта, Ф.К. Салманова, В.А. Серебрякова, А.К. Урупова, В.Д. Швецова, Л.П. Шендерей и ряда зарубежных авторов:У. Фертля, К. Магары, М. Хабберта, В. Иллинга, Д. Девиса, К. Хотгмана, Р. Чемпена и других.

Актуальность работы

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

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

Цель работы

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

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

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

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

3. Создание методики обобщенного профиля для оценки достоверности определения кинематических характеристик и структурных построений.

4. Разработка технологии и методики геосейсмического моделирования.

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

Методика исследований включает: анализ и использование широко применяющихся в настоящее время методов и способов прогноза аномально высоких пластовых и поровых давлений (АВПД, АВПоД); анализ геологических результатов математического, геосейсмического, тектонофизического моделирований; прогноз напряженного состояния нефтегазовых коллекторов и глинистых покрышек (АВПД и АВПоД) на месторождениях севера Западной Сибири, с использованием разработанных при участии автора программно-алгоритмических средств и методик; изучение и анализ многочисленных публикаций по проблеме исследований.

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

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

2. Установлено влияние термобаричееких условий на физические свойства терригенных отложений.

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

4. Разработана методика и технология геосейсмического моделирования.

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

6. Установлена возможность использования прогнозных параметров напряженного состояния горных пород по геофизическим данным в качестве интерпретационных признаков для выделения углеводородосодержащих объектов и геотектонически -активных зон.

7. Показана перспективность использования разработок как основы для создания региональной базы напряженного состояния осадочного чехла Западной Сибири.

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

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

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

3. Методика и технология площадного прогноза оценки свойств флюидоупоров (глинистых покрышек) нефтегазовых коллекторов.

4. Технология и методика геосейсмического моделирования.

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

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

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

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

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

4. Прогнозное распределение АВПоД по площади для глинистых покрышек до стадии бурения позволяет оценить сохранность углеводородов в выявленных ловушках и одновременно может являться поисковым признаком.

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

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

Реализация разработки и ее апробация

Результаты, полученные при выполнении настоящей работы, апробированы на материалах сейсморазведочных съемок и скважинных исследований по Бованенковской группе месторождений, Самбургскому, Ем-Еговскому, Маслиховскому и Ай-Пимскому месторождениям. Отдельные элементы разработки применены при оперативном прогнозе по Гыданскому району на площадях: Лунная, Минховская, Тото-Яхинская, Приречная, Глубинная, Сопочная, Овражная. Разработанная методика и технология прогноза напряженных свойств флюидоупоров (глинистых покрышек) апробирована на Пулытьинской-2, Кондырской площадях и Семивидовско-Толумском месторождении. Оценка точности определения интервальных скоростей и достоверности структурных построений выполнена по Бованенковской группе месторождений (семь месторождений), Ново-Надымском, Яунлорском, Восточно-Сургутском, Чапровском, Рогожниковском месторождениях.

Основные результаты исследований докладывались на:

Всесоюзной научно - практической конференции в г. Ленинграде, 1990г; на семинарах и конференциях, проводимых ЦГЭ, посвященных проблемам прогнозирования геологического разреза (ПГР); на международной геофизической конференции SEG / Москва - 92; на международной геофизической конференции SEG - ЕАГО/Москва-93; на конференции «Геомодель-2002, Геомодель-2003, Геомодель-2004»; на 5,6, 8 научно-практических конференциях «Пути реализации нефтегазового потенциала ХМАО» в 2002,2003,2005 гг. в г. Ханты-Мансийске.

Личный вклад

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

Исследования автора проводились в период 1984-2006 гг. Завершение работы выполнено в период прикрепления автора соискателем ученой степени доктора технических наук приказом № 7-А от 17 февраля 2006 г. по Российскому государственному университету нефти и газа им. И.М. Губкина.

Публикации

Основные положения диссертации представлены более чем в 40 научно-исследовательских, тематических работах и 21 публикации, в том числе 14 работ опубликованы в двух монографиях и статьях в ведущих рецензируемых научных журналах и изданиях согласно перечню ВАК.

Объем работы

Работа состоит из введения, 6 глав и заключения. Содержит 181 страницу машинописных листов, 73 рисунка, 11 таблиц, одно приложение, библиография включает 109 наименований.

Благодарности

Автор выражает искреннюю признательность за сотрудничество, консультации и ценные советы при работе над технологиями и методиками моим коллегам: д.г. - м.н. Ахиярову В.Х., к.т.н. Анпенову С.В., к.г. - м.н. Белкину Н.М., д.г.-м.н. профессору Бондаренко П.М., д.ф. - м.н. Голошубину Г.М., к.г. - м.н. Гальченко А.И. д.т.н. Ждановичу В.В., к.ф. - м.н. Кравченко Г.Г., д.т.н. Курьянову Ю.А., к.г.-м.н. Теплоухову В.М.

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

Автор благодарен коллегам по работе Ковалеву В.Г., Кондратенко М.А., Ширшову А.С., Кузнецовой Л.Г., Романову Д.В. за помощь, оказанную в оформлении графики, Коровкину М.А. за реализацию алгоритмов при создании системы обмена информации.

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

Выход

1 8 1

•to •-JM* kvl •to 14- И7-Ю tol 1»* 14- J-« •-11-51 kwj jum IMVI*

•to )-♦! 11-»- s itoeJ *rr JUm К-Н-» tin »- 4 — U-11-54

П ' ЙН4И1Г> НЩМ

Г4 - Рвщлктвг Ktmic пни

Ft - Mt«t*»mMnrt П • Chvih HIIM fit * itiH me hw mimm Каш - И* «им Itiwr* mum* fai* |H4I tkuJ «to It- MI 1-И-И dK) •to 14- Ml •-14-31 vtol •to 14- Hi ИМ» tot •to 14- 3-H Ml -«

Ik»] •I» u-n-i»

-<Шв •to т-« U-«- 3

Cat**tr» m«4)

ШШЯШМй ь t< I i - T «Cvri ■ n-f > I N«f «И •< f * W

I ■ M"«ftM) И

It - «at 44«^mn4»Hiint»« >

- a» ihlHtw

- Щчмьц* **»*< rt»4l>w) If

11

-'JJ

ПП<я«в ElC-biKH

Рис. 6.1.1. Технологическая схема сейсмического моделирования

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

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

Свертка импульса с коэффициентом отражения предназначена для расчета синтетических трасс одномерных моделей. Временной разрез для двумерных моделей представляет собой набор синтетических трасс.

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

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

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

Модуль настройки предназначен для настройки системы и указания размеров динамически распределяемой памяти для массивов. Данные между блоками передаются в форме файлов определенной структуры.

Математический аппарат позволяет методом подбора и последовательных приближений уточнять параметры объектов содержащих углеводороды и объектов подверженных литологической изменчивости по латерали.

Основным преимуществом рассмотренной технологии, по сравнению с существовавшими раньше, является:

- диалоговый режим работы;

- отсутствие ограничений по толщине слоев;

- оперативный доступ к информации;

- неограниченное количество вариантов модели.

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

6.2. Методика геосейсмического моделирования и практические решаемые задачи 6.2.1. Подбор модели среды

Рассмотрим несколько примеров влияния на формирование волнового поля:

- типа насыщения и литологического замещения (глинизации);

- выклинивания тонкослоистого пласта коллектора;

- изменение физических свойств пород коры выветривания (трещинные коллектора) и свойств вмещающих пород.

На рисунке 6.2.1.1 представлен результат подбора скоростных и плотностных параметров в разрезе скважины 12 по профилю 6 ПК-104.25 до максимального приближения реальной и синтетической трасс на основе одномерного геосейсмического моделирования с учетом и оценкой влияния каждой границы на формирование волнового поля. Оценка влияния производится последовательным исключением границ. Графики параметров на рисунке 6.2.1.1 и их значения в таблице 6.2.1.1 представлены для исследуемого интервала с 39 по 64 границы.

Результаты подбора по 26 границам (синтетический разрез) представлены на рисунке 6.2.1.2. а. В процессе подбора производится уточнение значений пластовой скорости (Упл.) и плотности (а). Порог уточнения для скоростей составляет 11%, а для плотностей 8%. Окончательный коэффициент корреляции между подобранной синтетической и реальной трассами составляет 0.931, что свидетельствует о высокой степени достоверности уточнения параметров скорости и плотности.

Сопоставление фрагмента временного (рис. 6.2.1.2 б) и синтетического (рис. 6.2.1.2 а) разрезов совместно с литологией, полученной по результатам ГИС в скважине 12, позволяет установить принадлежность литологических границ 46,47 и 50,51 (кровля и подошва исследуемых объектов-коллекторов ачимовской толщи) отраженной волне на времени 2,2 секунды. Интерференционный характер ОВ на этом времени объясняется в большей степени малой толщиной этих пластов (БС16-18 - 6м, БС19.20 - 8м) и в меньшей типом флюидонасыщения.

Резкий перепад акустической жесткости за счет характера насыщения, особенно а) б)

Рис. 6.2.1.1. Подбор параметров пластов геологического разреза ПР-6, ПК-104,25, скв.12 графики исходных параметров геосейсмической модели; - [ рафики подобранных параметров геосейсмичской модели

Упл. - пластовая скорость: ST - синтетическая трасса; а - плотность; RT - реальная трасса; G акустическая жесткость. б)

IBHBBBS г (I • г < ■ i f i i 'i '(((ii'i.'(' 1 f: 111 f 11 f i f 111 in fit ттШШШВЗ^т »i "И Щ HI fM П'.Ш.ПИИ III 111 111 . 11 lii'iiT' номера границ i К * '? . 'й " ™ * V !! V мттттттмтжжмъ a i

Рис. 6.2.1.2. Влияние типа флюидонасыщения и литологической изменчивости пластов BCi6-ib, БС19-20 на волновое поле продольных волн а) волновое попе при нефтзнасыщэнии дя глина б) фрагмент реального сейсмического разреза через скв 12 i I песчаник в) волновое попе при глинизации ^Н аргиллит г) волновое поле при водонэсыщении алевролит

Заключение

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

- автоматизированной технологии площадного прогноза аномальных давлений;

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

- методики площадного прогноза АВПоД для оценки свойств флюидоупоров (глинистых покрышек) нефтегазовых залежей;

- технологии геосейсмического моделирования;

- системы обмена и хранения геолого-геофизической информации.

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

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

Показана возможность повышения точности прогноза поровых давлений и обеспечения их оценок с погрешностью ±ЗМПа при решении задачи площадного

Показана возможность повышения точности прогноза поровых давлений и обеспечения их оценок с погрешностью ±ЗМПа при решении задачи площадного прогноза АВПД на основе сейсмических данных. С привлечением результатов интерпретации ГИС на стадии поисково-разведочных работ может быть достигнута меньшая погрешность.

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

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

Библиография Диссертация по наукам о земле, доктора технических наук, Ибраев, Валерий Иванович, Москва

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

2. Александров Б.Л. Аномально-высокие пластовые давления в нефтегазоносных бассейнах.-М.: Недра, 1987. 216 с.

3. Александров Б.Л., Афанасьев B.C. Определение и учет факторов, влияющих на достоверность количественной оценки АВПД по электрометрическим и акустическим данным.-М: Изд. АН СССР МНП ИГиРГИ, 1980. -41с.

4. Аникиев К.А., Введенский А.Я., Шендерей Л.П. Классификация гипотез происхождения аномально высоких пластовых давлений (АВПД). М.: Геология и разведка, 1985.-С. 93-99.

5. Берзон И.С, Сейсмическая разведка вертикально-слоистых сред фундамента. М.: Недра, 1977.-320 с.

6. Булатов В.В. Глубинная геомеханика. -М.: Недра, 1990. 264 с.

7. Буряковский Л.А., Джеваншир Р.Д., Алияров Р.Ю. Геофизические методы изучения геофлюидальных давлений.-Баку: ЭЛМ, 1986. 146 с.

8. Введение в сейсмическую интерпретацию: Пер. с англ. / Р. Мак-Куиллин, М. Бекон, У. Барклай. М,: Недра, 1985. - 308 с.

9. Вычислительные математика и техника в разведочной геофизике: Справочник геофизика /Под ред. В.И. Дмитриева. М.: Недра,1990. - 498 с.

10. Войтенко B.C. Управление горным давлением при бурении скважин. М.: Недра, 1985. - 178 с.

11. Гальперин Е.А. Вертикальное сейсмическое профилирование. М: Недра, 1971. -344 с.

12. Гальченко А.И. Способы оценки и прогноза геофлюидальных давлений на основе данных ГИС в нефтегазовых районах Западной Сибири. Диссертация на соискание ученой степени кандидата геолого-минералогических наук. Тюмень, ЗапСибНИИГеофизика, 1991. - 216 с.

13. Геология нефти и газа Западной Сибири.// А.Э. Конторович, И.И. Нестеров, Ф.К. Салманов и др.,-М.: Недра, 1975. 680 с.

14. Геология нефти и газа Сибирской платформы // Под ред. А.Э. Конторовича, B.C. Суркова, А.А. Трофимука. М.: Недра, 1981. - 509 с.

15. Гельчинский Б.Я., Белозеров А.А., Бердинникова Н.И. и др. Метод эффективной сейсмической модели. Л., Изд-во Ленинградского ун-та, 1975.

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

17. Гогоненков Г.Н. Прогнозирование геологического разреза по сейсмическим данным .- М.: Геология нефти и газа, 1981, № 1. С. 48-55.

18. Гогоненков Г.Н. Расчет и применение синтетических сейсмограмм, М.: Недра, 1972.-С. 54-73.

19. Гогоненков Г.Н., Авербух А.Г., Птецов С.Н. и др. Программно-методический комплекс для решения задач прогнозирования геологического разреза (СЦС-3-ПГР).-М.: ЦГЭ, 1984.

20. Гольдберг В.М., Скворцов Н.П. Проницаемость и фильтрация в глинах. М.: Недра, 1986.- 160 с.

21. Гольдин С.В. Линейные преобразования сейсмических сигналов. М.: Недра, 1974.-352 с.

22. Гельфапд В.А. Уточнение модели среды с помощью синтетических сейсмограмм. Нефтегазовая геология и геофизика, 1977 г. № 5. С. 32-36.

23. Гурвич И.И., Боганик Г.Н. Сейсмическая разведка М.: Недра, 1980.

24. Дахнов В.Н. Интерпретация результатов геофизических исследований разрезов скважин. И.: Недра, 1982.-448 с.

25. Добрынин В.М., Лимбергер Ю.А. Определение аномальности высоких пластовых давлений по электрометрии скважин, 1971, № 10. С. 69-74.

26. Добрынин В.М., Серебряков В.А. Методы прогнозирования аномально высоких пластовых давлений. М.: Недра, 1978. - 232 с.

27. Добрынин В.М., Серебряков В.А. Геолого-геофизические методы прогнозирования аномальных пластовых давлений. М.: Недра, 1989. - 287 с.

28. Дмитриев В.И., Морозов, В.А. Жданов М.С. и др. Вычислительные математика и техника в разведочной геофизике. М.: Недра, 1990.

29. Добрынин В.М. Деформация и изменение физических свойств коллекторов нефти и газа. М.: Недра, 1970. - 239 с.

30. Добрынин В.М., Вендельштейн Б.Ю., Кожевников Д.А. Петрофизика. М.: Недра, 1991.-368 с.

31. Епинатьева A.M., Голошубин Г.М., Литвин А.Л. и др. Метод преломленных волн.-М.: Недра, 1990.-297 с.

32. Жданович В.В. и др. Интегрированная система комплексной интерпретации данных сейсморазведки и ГИС на ЕС ЭВМ.- Прага, тр. 33 Международного геоф. конгресса, 1988, т.1. С 139-147.

33. Желтов Ю.П. Деформация горных пород. М.: Недра, 1966. - 197 с.

34. Желтов Ю.П. Механика нефтегазоносного пласта. -М.: Недра, 1975. 216 с.

35. Ибраев В.И., Пустовит В.И. Формирование волновых полей методами математического моделирования. Тюмень: изд. ЗапСибНИИГеофизика, 1989. -С. 84-89.

36. Ибраев В.И. Обработка и интерпретация материалов пространственной сейсморазведки программными средствами СЦС-3 ПГР и ИНТЕРСЕЙС. Геология нефти и газа, № 7, М., 1990. С. 27-31, (совместно со Ждановичем В.В., Щекиным С.Н.).

37. Ибраев В.И. Новые возможности изучения скоростных разрезов в сейсморазведке. SEG/Москва, 1992. - С. 408, (совместно с Епинатьевой A.M., Соловьевым B.C., Голошубиным Г.М., Рябиновой Л.Г.).

38. Ибраев В.И. Прогноз аномальных давлений в терригенных отложениях Западной Сибири на основе комплексирования геофизических методов. М.: Недра, «Геология нефти и газа», № 11,1993. - С. -23-26 .

39. Ибраев В.И., Теплоухов В.М. Отчет "Результаты экспериментальных опытно-методических работ по изучению низкочастотной составляющей волнового поля и скоростной анизотропии (дог. 37). Тюмень, 1991.

40. Ибраев В.И. Геосейсмическое моделирование при решении задач выделения тонкослоистых коллекторов// Специальный выпуск журнала Геофизика «Технологии сейсморазведки-2». -М.: ЕАГО, 2003. С. 22-26, (совместно с. Белкиным Н.М.).

41. Ибраев В.И. Интерпретация материалов объемной сейсморазведки (совместно со Ждановичем В.В., Потаповым О.А., Кузнецовым В.И.) Обзор Разведочная геофизика. - М: ВИЭМС 1991. - 65 с.

42. Ибраев В.И. Прогнозирование АВПД по данным ГИС-сейсморазведка в геологических условиях Западной Сибири, (совместно с Анпеновым С.В., Гальченко A.M., Теплоуховым В.М.) Сборник ВНИГРИ: АВПД и нефтегазоносность, г. Ленинград, 1992.

43. Ибраев В.И. Прогноз напряженного состояния верхнеюрских отложений Сургутского свода по комплексу ГИС сейсморазведка. Геология нефти и газа 2002/5, (совместно с Белкиным Н.М., Косом И.М.). - С. 42-45.

44. Ибраев В.И. Частотно-резонансный анализ сейсмических данных. В кн. Пути реализации нефтегазового потенциала ХМАО (5-ая научно-практическая конференция). Ханты-Мансийск. 2002. Том 1. С. 310-313, (совместно с. Белкиным Н.М.).

45. Ибраев В.И. Разработка технологии прогноза напряженного состояния терригенных отложений на основе комплексной интерпретации данных сейсморазведки и ГИС в условиях Западной Сибири. Диссертация на соискание ученой степени кандидата технических наук.

46. Ибраев В.И. Динамический анализ и моделирование для выделения границ фациального замещения коллекторов абалакской свиты. Геология нефти и газа, №4, М., 2006. С. 32-38.

47. Ибраев В.И. Методика площадной оценки качества флюидоупора (глинистой покрышки). Геофизика № 3, М., 2006. С. 27-29.

48. Ибраев В.И. Влияние термобарических условий на физические свойства терригенных отложений. Известия высших учебных заведений «Нефть и газ», №3,2006.-С. 4-10.

49. Ибраев В.И. Методика структурно-динамической интерпретации временных сейсмических разрезов с позиций тектонофизического моделирования. Геофизика №4, М., 2006. С. 9,10,30,31.

50. Ибраев В.И. Методика обобщенного вертикального профиля интервальных скоростей. Известия высших учебных заведений «Нефть и газ», № 4, 2006. -С. 6-10.

51. Ибраев В.И. Монография «Становление интерпретации материалов объемной сейсморазведки в Западной Сибири», Тюмень: ОАО «Тюменский дом печати», (совместно со Ждановичем В.В., Потаповым О.А., Кузнецовым В.И.), 2006. 82 с.

52. Ибраев В.И. Монография «Прогнозирование напряженного состояния коллекторов и флюидоупоров нефтегазовых залежей в Западной Сибири», Тюмень: ОАО «Тюменский дом печати», 2006. 208 с.

53. Интерпретация результатов геофизических исследований нефтяных и газовых скважин: Справочник./Под ред. В.М. Добрынина. М.: Недра, 1988. - 476 с.

54. Карцев А.А. Гидрогеологические условия проявления сверхгидростатических давлений в нефтегазоносных бассейнах. М.: Геология нефти и газа, 1980, № 4. -С. 40-43.

55. Кондратьев O.K. Сейсмические волны в поглощающих средах. М.: Недра,1988. -С. 40-53.

56. Комплектование геофизических методов при решении геологических задач/Под ред. В.Е. Никитского, В.В. Бродового. М.: Недра, 1987. - 471 с.

57. Корж М.В., Филина С.И., Зонн М.С. Применение палеогеографических методов для поисков нефтегазовых залежей в ловушках литологического и комбинированного типов. М., ВНИИ0ЭНГ, 1982.

58. Котельников Д.Д., Конюхов А.И. Глинистые минералы осадочных пород. М.: Недра, 1986.-247 с.

59. Котяхов Ф.И. Физика нефтяных и газовых коллекторов. М.: Недра, 1977.- 287с.

60. Кузнецов O.JL, Чиркин И.А., Курьянов Ю.А., Роговицкий Г.В., Дыбленко В.П. Экспериментальные исследования. М.: Государственный научный центр Российской федерации - ВНИИгеосистем, 2004. - Том 2. - С. 90-166.

61. Кунин Н.Я., Кучерук Е.В. Сейсмостратиграфия в решении проблем поиска и разведки месторождений нефти и газа / Итоги науки и техники. ВИНИТИ, том 13. М., 1984.

62. Кунин К.Я., Усенко С.В., Черкасов Г.Р. Области некомпенсированного прогибания Западной Сибири и их сейсмостратиграфическая характеристика, -М„ ВИНИТИ, 1985.-22 с.

63. Кучерук Е.В., Шендерей Л.П. Современные представления о природе аномально высоких пластовых давлений. М., ВИНИТИ 1975. 165 с.

64. Кучерук Е.В., Люстих Т.Е. Прогнозирование и оценка аномальных пластовых давлений по материалам геофизических исследований. М., ВИНИТИ, том 7,1986. -С. 70-115.

65. Ляховицкий Ф.М. Сейсмические волны в гетерогенных средах.// Междуведомственный геофизический комитет при Президиуме АН СССР. М., 1988.- 149 с.

66. Мармоштейн Л.М. Петрофизические свойства осадочных пород при высоких давлениях и температурах. М.: Недра, 1985. - 189 с.

67. Методические указания по прогнозу и оценке аномально высоких пластовых давлений (АВПД)/Л.,ВНИГРИ-1987. С. 15-30.

68. Магара К. Уплотнение пород и миграция флюидов// Перев. с англ. -М.: Недра, 1982.-295 с.

69. Мелик-Пашаев B.C. Халимов Э.М. Серегина В.Н. Аномально высокие пластовые давления на нефтяных и газовых месторождениях. М.: Недра, 1983.- 161 168 с.

70. Месчян С.Р. Механические свойства грунтов и лабораторные методы их определения. М.: Недра, 1974. - 192 с.

71. Мешбей В.И. Сейсморазведка методом общей глубинной точки. М.: Недра, 1973.- 152 с.

72. Мкртчян О.М., Гребнева И.Л., Игошкин В.П. и др. Сейсмогеологическое изучение клиноформных отложений Среднего Приобья. М.: Наука, 1990. - 108 с.

73. Муше Ш.П., Митчелл А. Аномальные пластовые давления в процессе бурения: Происхождение-прогнозирование-выявление-оценка: Техн. руководство: Пер. с англ.-М.: Недра, 1991. 287 с.

74. Нестеров И.И., Потеряева В.В., Салманов Ф.К. Закономерности распределения крупных месторождений нефти и газа в земной коре. М.: Недра, 1975. - 276 с.

75. Нестеров И.И. Уплотнение глинистых пород. М. Советская геология, 1965, № 12.- С. 13-18.

76. Никитин А.А. Теоретические основы обработки геофизической информации. М.: Недра, 1986.-342 с.

77. Николаевский В.Н. Механика пористых и трещиноватых сред. М.: Недра, 1984.229 с.

78. Овчинников И.К. Теория поля. М.: Недра, 1971.-312 с.

79. Орлов А.А. Роль покрышек в формировании и сохранении сверхгидростатических пластовых давлений в нефтяных и газовых залежах. Киев, Наукова думка, АН УССР, Сер. Геология, вып. 63,1984. - С. 1-6.

80. Орлов А.А. Тектонические условия формирования сверхгидростатических пластовых давлений в газо,-нефте,-водоносных структурах. -М., Геология и разведка, № 9, 1985. С. 116-119.

81. Петрофизическая характеристика осадочного покрова нефтегазоносных провинций СССР: Справочник // Под ред. Г.М. Авчяна, M.J1. Озерской. М.: Недра, 1985.- 192 с.

82. Программирование на Фортране 77 :Пер. с англ. / Дж. Ашкрофт, Р. Элдридж Р. Полсон, Г. Уилсон, М. :Радио и связь, 1990.-272 с.

83. Программно-методический комплекс для решения задач прогнозирования геологического разреза (СЦС-З-ПГР)/ Г.Н. Гогоненков, А.Г. Авербух, С.Н. Птецов и др.-М., ЦГЭ, 1984.

84. Птецов С.Н. Анализ волновых полей для прогнозирования геологического разреза.-М.: Недра, 1989. 135 с.

85. Рапопорт М.Б. Автоматическая обработка записей колебаний в сейсморазведке. -М.: Недра, 1973.

86. Рудкевич М.Я., Бочкарев В.С, Максимов Е.М. и др. Объяснительная записка и атлас структурных, палеонтологических карт. ЗапСибНИГНИ, вып. 28, Тюмень, 1970.- 174 с.

87. Сейсмогеологический анализ нефтегазоносных отложений Западной Сибири/ О.М. Мкртчян, J1.J1. Трусов, Н.М. Белкин, В.А. Дегтев. М.: Наука, 1987. - 126 с.

88. Старцев А.А., Симонов В.И. Состояние и пути совершенствования технологии глубокого бурения в Западной Сибири./Обзор ВИЭМС, М., 1987. - 45 с.

89. Статистика и анализ геологических данных: Пер. с англ. / Дж. Девис. М.: Мир, 1977.-569 с.

90. Сейсмическая стратиграфия. Использование при поисках нефти и газа // Под ред. Ч. Пейтона. Пер. с англ., в 2 ч.-М.:Мир, 1982. 839 с.

91. Сейсморазведка: Справочник геофизика. В двух книгах /Под ред .В.П., Номоконова.-М.: Недра, 1990. 736 с.

92. Селли Р.К. Введение в седименталогию. М.: Недра, 1981. - 363 с.

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

94. Урупов А.К., Лапин С.И. Скорости сейсмических волн в анизотропных слоисто-трещиноватых средах/Прикладная геофизика; вып.67-М.: Недра, 1972. С. 3-16.

95. Ушатипский И.Н., Зарипов О.Г. Минералогические и геохимические показатели нефтегазоносности мезозойских отложений Западно-Сибирской плиты. Свердловск: Средне Уральское книжное издательство, 1978. - 208 с.

96. Фертль У.Х. Аномальные пластовые давления: Пер. с англ. М .- Недра, 1980.-Пер. изд. - 398 с.

97. Чемпен Р.Е. и др. Геология и вода./Пер. с англ. Л.: Недра, 1983. - 159 с.

98. Шаблинская Н.В. Разломы по геофизическим данным и их роль в формировании структур (Западно-Сибирской и Тимано-Печорской плит.//Разломы земной коры. -М.: Наука, 1977.-С. 160-166.

99. Эринчек П.Т. Опробование залежей с АВПД на площадях севера Западной Сибири. Изучение геологического разреза и прогнозирование АВПД.//Сб.науч.трудов, Л., ВНИГРИ, 1987. С. 138-142.

100. Rafipour B.J. Seismic Response for Reservoir-Fluid Evalution, SPE Formation Solution, March, 1989. p. 45-48.

101. Stein Natan Seismic data used to predict formation pressures//0il and Gas,-1992.-Vol.90, № 48. p. 55-58.

102. Hans Rische Kombinirte P-S-Seismik.//Zeitschrifl fur Angewandte Geologic, Leipzig,Band 35/ 1989/Heft 12.

103. F.Christian Hagner Das Bruchspannungs-/Reibungsspannungs-verhaltnis und Energieumsatz seismischer Ereignisse im Carnal lit.// Z. angew.Geol., 36/ 1990/5.