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

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

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

Ннкулышков Алексей Юрьевич

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

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

Автореферат диссертации на соискание ученой степени кандидата техшіческих наук

- 1 НОЯ 2012

Москва - 2012

005054143

005054143

Работа выполнена в Российском Государственном Геологоразведочном Университете имени Серго Орджоникидзе.

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

минералогических наук Ермолаева Галина Михайловна

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

кандидат геолого-минералогических наук, Жемчужников Евгений Глебович

Ведущая организация

доктор технических наук, Кириллов Сергей Александрович

ООО «Геофизическое партнерство»

Защита состоится 15 ноября 2012 года в 15 часов на заседании диссертационного совета Д 212.121.07 при Российском Государственном Геолого-разведочном университете им. Серго Орджоникидзе по адресу: 117997 Москва ул. Миклухо-Маклая д.23, аудитория 4-73.

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

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

ный секретарь

гертационного совета доктор А.Д. Каринский

ико-математических наук /

Актуальность

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

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

Объекты исследований

Объектами исследований являлись:

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

- данные сейсморазведки ЗЭ ряда площадей Западной Сибири и Латинской Америки.

Цель исследования

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

Задачи исследования

1. Анализ существующих алгоритмов вейвлет-преобразования.

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

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

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

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

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

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

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

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

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

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

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

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

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

Выполнено 2D моделирование клиноформного песчаного пласта с целью обоснования эффективности применения вейвлет-преобразования при прогнозе эффективной нефте- газонасыщенной мощности геологических объектов и других ФЕС.

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

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

Практическая значимость результатов

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

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

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

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

Реализация работы в производстве

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

Объем и структура работы

Диссертация состоит из введения, трех глав и заключения общим объемом 105 страниц, включая 39 рисунков и список литературы.

Личный вклад автора

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

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

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

Автор благодарит коллектив ООО «Геофизические Системы Данных» за интересные идеи, советы и поддержку при написании диссертационной работы.

СОДЕРЖАНИЕ Введение.

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

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

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

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

пиковой частоты сейсмического сигнала от временной мощности тонкого пласта

8

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

1. Алгоритмы вейвлет-преобразования.

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

Грубую классификацию алгоритмов вейвлет-анализа можно сделать, выделив непрерывное (CWT — Continuous Wavelet Transform) и дискретное (DWT — Discrete Wavelet Transform) вейвлет-преобразования. Получить набор вейвлет-коэффициентов в случае дискретного преобразования быстрее, и оно дает достаточно точное представление о сигнале при меньшем объеме получаемых в результате данных. Непрерывное преобразование требует больших вычислительных затрат, но, вместе с этим, позволяет детальнее рассмотреть структуру сигнала. Вследствие бурного развития вычислительной техники большое количество вычислений при выполнении математических алгоритмов не является серьезным недостатком. Поэтому при выборе алгоритмов спектрального анализа нестационарных сигналов (в первую очередь сейсмической записи), исходят из их эффективности. Основное отличие вейвлет-преобразования от преобразования Фурье заключается в том, что оно обеспечивает двумерную развертку одномерного сигнала. При этом частота и

9

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

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

1 . ¿о

/<». г) = — / Ги (<г. т)Мош)е-ш- -уг.

7-х от-

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

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

О 50 100 1»

Рис. 1. Пример частотно-временного амплитудного спектра сейсмической трассы.

2. Анализ и интерпретация результатов вейвлет-преобразования.

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

Существует несколько основных способов интерпретации результатов вейвлет-преобразования. Первый способ заключается в трансформации амплитуд в частотный диапазон (другими словами, результатом являются кубы частот, где по направлению Z отображены амплитуда или фаза в диапазоне от 1 до 100 Гц). Такой объем называется «тюнинг кубом». По отдельности эти результаты трудно интерпретируемы, особенно количественно. Анализ тюнинг-куба оптимизируется с помощью цветового смешивания (RGB), что предоставляет возможность интерпретатору оценить резервуар надежнее, в некоторых случаях давать количественную оценку его характеристикам (рис. 2). Так же, этот подход позволяет выделять структуры при наличии на территории исследований тонких песчаных газо- и нефтенасыщенных линз и каналов.

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

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

Рис. 2. Анализа результатов вейвлет-преобразования. а) амплитуды на частоте 20 Гц; б) амплитуды на частоте 40 Гц; в) амплитуды на частоте 60 Гц; г) результат

ЯСВ-смешивания.

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

2.1. Сейсмическое моделирование характеристик тонкого пласта.

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

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

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

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

Разница в частоте отражения между амплитудным спектром, рассчитанным в большом (рис. 3) и коротком окне (рис. 4) очевидна.

Коэффициенты Сейсмическая

отражения * Импульс + Помеха .—у с ческая

трэссэ

г £

х 5

« га

II

Ш

V

Амплитуда

V

Амплитуда Амплитуда

5 (О

5 О

-*

>1

г «

у 1 р 1 V

>

V

Амплитуда

всей длине. (РагШса 1999)

Коэффициенты отражения ^

„ ^Сейсмическая

Импульс + Помеха с=ф трасса

§ ¡5

т о

о а

ж с П) из о. о са

Сп&ктр имупьса

Рис. 4. Модель сейсмической трассы. Переход в частотную область с использованием короткого временного окна (Рагйка 1999).

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

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

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

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

Свойство наложения отражений от тонких пластов может быть проиллюстрировано с помощью простой модели клина (рис. 5). Временное отражение состоит из двух пиков одинаковой, но противоположной амплитуды. Кровля клина отмечена отрицательным коэффициентом отражения, а подошва положительным. Мощность клина изменяется от 0 до 50 метров. При свертке геологической модели с импульсом (8-10-40-50 Гц) происходит эффект наложения, вызванный изменением в мощности. Отражения от кровли и подошвы разделяются при большой мощности и накладываются в одно отражение при истончении клина.

Амплитудный спектр (рис. 6) был рассчитан в коротком окне. Временная мощность клина определяет период максимумов амплитудного спектра в зависимости от частоты. 1/Р, где Р (Гц) - период максимумов амплитудного спектра от частоты, и I (с) - мощность пласта.

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

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

Модель коэффициента отражения , 0_10

Временная мощность (с) 20_30

40

і 100'

200'

Модель свернутая с сейсмическим импульсом Временная мощность (с)

10

50

1 100' &

200-

' І М I і 11 > И М і • і М і I ! 11 і І 111 [(І

(і > >>.>>>> к>>>> ьі>.іХ>.> ''Д'.АЛ* 'ЛКк-V' А ККЬ?ЛЛ

50

Временная мощность

Рис. 5. Сейсмическая модель клиноформного пласта. (Рагііка 1999) 1

Временная мощность 50 мс

Рис. 6. Амплитудно-частотный спектр сейсмической модели.

20 40 60 80 100 Частота (Гц)

2.2. Визуализация и анализ результатов спектральной декомпозиции с помощью RGB-смешивания.

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

їмо

1050 1100 uso 1200 С 12»

що

1350 1400 1450

Рис. 7. Частотно-временные амплитудные спектры сейсмической записи.

На следующем этапе полученное представление сейсмической трассы трансформируется в так называемый «тюнинг куб». Это название подразумевает сейсмический куб амплитуд, включающий в качестве четвертого измерения -частоту. «Четырёхмерное» представление сейсмических данных после спектральной декомпозиции позволяет выполнять детальный анализ полученных результатов. На рис. 8. показаны сечения «тюнинг куба» вдоль интерпретируемого горизонта. Как видно из данной иллюстрации, в центральной части исследуемой площади распределение амплитуд на разных частотах варьируется. Тем более, картина совершенно иная по сравнению с абсолютными амплитудами (рис 8А). Это вызвано тем, что в отличие от изображений, полученных с помощью вейвлет-преобразования, здесь в значения амплитуд

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

Рис. 8. Значения амплитуд куба данных 3D: А - абсолютные (исходные) амплитуды; Б - амплитуды после вейвлет-преобразования на 60 Гц; В -амплитуды после вейвлет-преобразования на 100 Гц; Г - амплитуды после вейвлет-преобразования на 140 Гц.

Для удобства визуализации и дальнейшей интерпретации результатов

вейвлет-преобразования используется алгоритм RGB-смешивания, позволяющий

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

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

рассчитывается значение цвета пикселя в красно- зелено- синем спектре.

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

конкретной точке (рис. 9). Как видно из данной иллюстрации, на результате RGB

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

18

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

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

Рис. 9. Результат RGB-смешивания амплитуд сейсмической записи после вейвлет-преобразования.

вейвлет-преобразования.

2.3. Типизация геологического разреза с помощью атрибутного анализа.

По данным многих исследователей туронские продуктивные отложения имеют региональное распространение, однако их изучению уделялось недостаточно много внимания. К наиболее изученным из них относятся газовые залежи газсалинской пачки Южно-Русского, Заполярного, Харампурского, Новочасельского, Тэрельского и Ленского месторождений.

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

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

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

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

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

Рис.

11. Классификация точек в трехмерном пространстве атрибутов.

mi« Numb« Po>№ Piopcrtion Distortion Porosity Imp

None 175 00875 1231552 0232035 4S84315

None 345 01725 1 021534 0 210456 «19.491

Ncne 391 0 1955 1.050W5 0199509 4950582

None 515 02575 033684 0220366 4326.771

Nont 574 0287 0855328 0203012 4823.772

Amp

0.576068 | ЛуЧШИЙ

0425901 коллектор

0.228564 0235649 0186362

3. Статистический анализ атрибутов сейсмической записи.

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

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

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

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

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

х2к

где п - число наблюдений; к - число свойств.

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

23

-4 1

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

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

V = + Й.ДГ. -... + -г Ь,

где XI, х2, ..., хк- свойства; аь а.2, ...,аь ¿-постоянные коэффициенты.

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

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

На рис. 12 приведен пример анализа данных вейвлет-преобразования на различных частотах и атрибута пористости по методу главных компонент. Анализ матриц факторных нагрузок и корреляционных связей позволяет говорить, что наиболее информативным является первый фактор, в который больший вклад делают амплитуды на 50 и 55 Гц. Это свидетельствует о том, что амплитуды на данных частотах имеют максимальную корреляцию с данными пористости пласта. Такая тесная связь говорит о фактическом отображении литологических свойств пласта на атрибутах амплитуды 50 и 55 Гц.

45 Гц 50 Гц 55 Гц 60 Гц

РСА1 0.436255 0.515114 0.515249 0.482423

РСА2 0.59491 037344 -0.359717 -0.614189

РСАЗ -0.453769 0366075 0.587594 -0.561084

РСМ -0.451366 0679109 -0.50976 0.27427

45 Гц 50 Гц 55 Гц 60 Гц

РСА1 0510912 0.96497 0965222 030373

РСА2 0.406795 0-255355 -0.24597 ■0.419977

РСАЗ -00665381 0.0536786 0.0861605 ■0.0822739

РСМ -0.0181609 0.0273205 -00205099 00110333

РСА1 VI АтрЬикЬ АшрМцбс ДтрМшЯс БрмР ЗО-ТОбу- РСА1 л - АтрИиНс 1ресО 30 70Сц1-

1 |4—^

1. ■¿и уИ1

45 Гц -2 0 2 4 6 РСМ 60 Гц -2 0 2 46 РСА1

РСА1 ¿трМийс Дгт,Е1....1. Дпяр>4и<1е 30-70*5^ РСА1 ДлфЫивс Длпри>к>е ■ ДтрМиск 5рееО 30 Т0ДО

1П рЬЬ З.СИ0

а л

50 Гц

55 Гц

Рис. 12. Таблица факторных нагрузок (а) и коэффициентов корреляции (в) компонент анализа, а так же изображение анализируемых данных в оси первой

компоненты.

Заключение

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

1. Частотно-временная модификация вейвлет-преобразования является наиболее надежным способом спектрального разложения сейсмической записи.

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

3. Технология 1ШВ-смешивания энергии единичных частот позволяет в значительной степени повысить надёжность определения положения стратиграфических элементов геологической среды в пространстве.

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

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

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

По теме диссертации опубликованы следующие работы:

1. Никульников А.Ю., Применение спектральной декомпозиции для повышения эффективности совместной интерпретации данных сейсморазведки 3D и ГИС: Международная научная конференция «Молодые - наукам о Земле» -М, 2011

2. Никульников А.Ю., Применение атрибутного анализа сейсмических данных с использованием результатов спектральной декомпозиции: Геоинформатика, M., 1., 2012

3. Никульников А.Ю., Применение спектральной декомпозиции сейсмической записи с целью картирования тонких стратиграфических элементов осадочных пород: Геология, геофизика и разработка нефтяных и газовых месторождений, М., 11., 2012

Другие опубликованные работы автора:

1. Никульников А.Ю., Определение скоростей распространения упругих продольных волн по данным детальной сейсморазведки 2D с учетом скважинных данных ГИС: Международная научная конференция «Молодые -наукам о Земле» - М., 2008

2. Никульников А.Ю., Геологическая эффективность применения миграции сейсмограмм до суммирования на примере месторождений Западной Сибири: Международная научная конференция «Геоперспектива - 2009» - М., 2009

3. Зимирева B.JI., Славкин B.C., Никульников А.Ю. Влияние учёта неоднородностей В 4P на динамические характеристики сейсмической записи в условиях Западно Сибирского НГБ: Технологии сейсморазведки,М.,2., 2010

4. Лобанов A.M., Филатов В.Г., Петров A.B., Овсепян M.JL, Гласко Ю.В., Беспрозванный П.А., Зиновкин C.B., Никульников А.Ю., Филатов Г.В., Интропродолжение и эпигенетическое магнитоминералообразование в нефтегазоразведке: М., 2009

Подписано в печать 12 октября 2012 г.

Формат 60x90/16

Объём 1,5 п.л.

Тираж 100 экз.

Заказ № 121012445

Оттиражировано на ризографе в ООО «УниверПринт»

ИНН/КПП 7728572912\772801001

Адрес: г. Москва, улица Ивана Бабушкина, д. 19/1.

Тел. 740-76-47, 989-15-83.

http://wmv.univerprint.ru

Содержание диссертации, кандидата технических наук, Никульников, Алексей Юрьевич

Оглавление.

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

Список иллюстраций.

Список сокращений.

Введение.

1. Теоретические основы вейвлет-анализа.

1.1. История возникновения и обзор существующих технологий вейвлет-преобразования.

1.2. Преобразование Фурье в фиксированном окне (вТРТ).

1.3. Частотно-временное непрерывное вейвлет-преобразование (ТРС\¥Т).

1.4. Выводы.

2. Анализ и интерпретация результатов вейвлет-преобразования.

2.1. Расчет временной мощности тонкого пласта.

2.2. Визуализация и анализ результатов спектральной декомпозиции с помощью ЫСВ-смешивания.

2.3. Типизация коллекторов с использованием результатов вейвлет-преобразования.

2.4. Выводы.

3. Статистический анализ атрибутов сейсмической записи.

3.1. Одномерный статистический анализ и его применение.

3.2. Двумерный статистический анализ.

3.3. Многомерный статистический анализ.

3.4. Выводы.

Введение Диссертация по наукам о земле, на тему "Технология интерпретации результатов вейвлет-преобразования сейсмической записи"

В настоящее время запасы углеводородов уникальных и крупных месторождений на территории Российской Федерации и стран СНГ выработаны уже более чем на 50 %. В этих условиях все большую роль приобретает освоение небольших залежей с трудно извлекаемыми запасами.

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

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

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

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

3.4. Выводы.

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

Заключение.

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

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

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

3. Технология ЯОВ-смешивания энергии единичных частот позволяет в значительной степени повысить надёжность определения положения стратиграфических элементов геологической среды в пространстве.

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

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

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

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

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

1. Авербух А. Г. Методика интерпретации данных сейсморазведки при интегрированном изучении нефтегазовых резервуаров. Геофизика, №1, ЕАГО, М., 1998, с.13-19.

2. Ампилов. Ю. П. Сейсмическая интерпретация: опыт и проблемы. М., Геоинформмарк, 2004. 278 с.

3. Ампилов. Ю. П. От сейсмической интерпретации к моделированию и оценке месторождений нефти и газа. М., СПЕКТР, 2008, 384 с.

4. Айвазян С.А., Бежаева З.И., Староверов О.В. многомерных наблюдений. М, Статистика 1974.

5. Боганик Г. Н., Гурвич И. И. Сейсморазведка. 2006, Тверь, АИС,743 с.

6. Голикова Г. В., Санников К. Ю., Мочалов А.П. Некоторые элементы структуры резервуара, определяемые по полям интерференционных волн, Технологии сейсморазведки № 4, с. 35 416 2010.

7. Дунаева К.А., Сагайдачная О.М., Сальников A.C. Вопросы применения вейвлет-преобразования сейсмических полей с реконструкцией во временную область, Технологии сейсморазведки, 2011, № 3, с. 71-83

8. Дьяконов В.П. Вейвлеты. От теории к практике. М.: СОЛОН- Р,2002.

9. Добеши И. Десять лекций по вейвлетам. Пер. с англ. Мищенко Е.В.-М.:РХД,2001

10. Дэвис Дж.С. Статистический анализ данных в геологии. Пер. с англ. В 2 кн./Пер. В.А. Голубевой: Под ред. Д.А. Родионова. Кн. 1. М.: Недра, 1990. -319 с.

11. ЗемцоваД.П., Никитин A.A., Пискун П.В. Вейвлет-волнового поля при решении детализационных задач сейсморазведки. Тезисы докладов УП-ой международной научно-практической конференции Геомодель-2005, Геленджик, 11-17 сентября 2005г., стр.68

12. Каплан С. А., Зотова Е.В., Шелавина Е.Ю.и др. Возможности прогноза интегральных емкостных свойств коллектора по данным ГИС и сейсморазведки, Геофизика, 1997, № 5, с. 9-13.

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

14. Кондратьев И. К., Рыжков В. И., Бондаренко М. Т.,Лапина Е. В. Эффективность прогнозирования коллекторов способами динамической интерпретации в Восточной Сибири, Технологии сейсморазведки № 4, с. 2634., 2010.

15. Кондратьев И.К., Бондаренко М.Т., Каменев С.П. Динамическая интерпретация данных сейсморазведки при решении задач нефтегазовой геологии, Геофизика, 1996, № 56, с. 41-47.

16. Короновский A.A., Храмов А.Е. Непрерывный вейвлегный анализ и его приложения. М.: Физматлит, 2003.

17. Крылов Д.Н. К оценке определения литологии и коллекторских свойств по данным сейсморазведки. Геология нефти и газа, №3, М., Недра, 1992, с.27-32

18. Кунин Н.Я., Кучерук Е.В. Сейсмостратиграфия в решении проблемы поиска и разведки месторождений нефти и газа. — В кн.: Месторождения горючих полезных ископаемых, т. 13, Итоги науки и техники, ВИНИТИ АН СССР, М., 1984, 196с.

19. Малярова Т.Н., Иванова H.A., Современные методы сейсмофациального анализа на реальных примерах /Тезисы докладов VIII международной научно-практической конференции Геомодель-2006, стр.136.

20. Михальцев A.B., Мушин И.А., Погожев В.Н. Обработка динамических параметров сейсморазведки. М., Недра, 1990, с.250.

21. Мушин И.А., Белоусов Г.А., Городков А.Б. СВАН-сейсморазведка. (спектрально-временной анализ в технологиях сейсморазведки). Геофизика 5/2005, 3-9 с.

22. Никитин А.А. Теоретические основы обработки информации. М.: Недра, 1986.

23. Поротов Г.С., Математические методы моделирования в геологии: Учебник, Санкт-Петербургский государственный горный институт (технический университет). СПб, 2006. 223 с.

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

25. Сагайдачная О.М., Дунаева К.А., Сальников А.С. Декомпозиция и анализ сейсмических полей на основе слоев вейвлет-разложения. Геофизика 5/2010.

26. Славкин B.C., Копилевич Е.А. Моделирование природных резервуаров нефти и газа на основе структкрно-литологической интерпретации данных сейсморазведки и бурения. — М., ВНИОЭНГ, 1995,167с

27. Трапезникова Н.А. Методика спектральных вариаций для прогнозирования свойств геологического разреза, Геофизика ,1997, № 2, с. 12-16.

28. Якимов И.Е. Разработка и исследование методов и технологий освоения трудноизвлекаемых запасов газа : диссертация . кандидата технических наук : 25.00.17.- Тюмень, 2008.- 175 е.:

29. Chakraborty, A., and D. Okaya, Frequency-time decomposition ofseismic data using wavelet-based methods: Geophysics, 60, 1906 -1916. , 1995

30. Castagna, J. P., S. Sun, and R. Seigfried, 2003, Instantaneous spectral analysis: Detection of low-frequency shadows associated with hydrocarbons: The Leading Edge, 22, 120-132.

31. Castagna, J., M. Burnett, E. Méndez, G. Rodríguez, L. García, J. Martínez, M. Avilés, and R. Villaseñor, 2003, Application of spectral decomposition to gas basins in Mexico: The Leading Edge, 22, 1130-1134.

32. Castagna, J., and S. Sun, 2006, Comparison of spectral decomposition methods: First Break, 24, 75-79.

33. Doyen, P., Seismic reservoir characterization: EAGE. , 2007

34. Giroldi, L., and F. Alegria, Using spectral decomposition to identify and characterize glacial valleys and fluvial channels within the Carboniferous section in Bolivia: The Leading Edge, 24, 1152-1159, 2005.

35. Grana, D. and E. Delia Rossa, 2010, Probabilistic petrophysical prop-erties estimation integrating statistical rock physics with seismic in-version: Geophysics, 75, no. 3, 021-037, doi: 10.1190/1.3386676.

36. Jianlei, L., Spectral decomposition and its application in mapping stratigraphy and hydrocarbons, Faculty of the Department of Geosciences University of Houston, A Dissertation, 2006

37. Mallat, S., and Z. Zhang, 1993, Matching pursuit with time frequency dictionaries: IEEE Transactions on Signal Processing, 41, 3397-3415.

38. Maklad, M. S., and J. K. Dirstein, 2007, Spectral detection of attenuation and lithology: CSPG CSEG Convention, 501-505.

39. Miao, X., and S. Cheadle, 1998, High resolution seismic data analysis by wavelet transform and matching pursuit decomposition: Geo-Triad, CSEG, CSPG and CWLS Joint convention, 31-32.

40. Partyka, G., J. Gridley, and J. Lopez, 1999, Interpretational applications of spectral decomposition in reservoir characterization: The Leading Edge, v. 18, p. 353-360

41. Partyka, G., L. Peyton, and R. Bottjer, 1998, Interpretation of incised valleys using new 3-D seismic techniques: A case history using spectral decomposition and coherency: The Leading Edge, 17, 1294-1298.

42. Peyton, L., R. Bottjer, and G. Partyka, Interpretation of incised valleys using new 3-D seismic techniques: A case history using spectral decomposition and coherency: The Leading Edge, 17, 1294-1298. , 1998

43. Sierra, J., H. Campos, M. Bonilla, W. Marin, D. Joseph, and S. Cardinez, 2008, Seismic multi-attribute integration for well prospect generation: South Main Soldado Field, Gulf of Paria, Trinidad & Tobago: XIV Congreso Venezolano de Geofísica.

44. Sinha, S. K., P. S. Routh, P. D. Anno, and J. P. Castagna, Time-frequency attribute of seismic data using continuous wavelet transform: 73rd Annual International Meeting, SEG, Expanded Abstracts, 1481-1484. , 2003

45. Steeghs, P., and G. Drijkoningen, 2001, Seismic sequence analysis and attribute extraction using quadratic time-frequency representations: Geophysics, 66, 1947-1959.

46. Stockwell, R. G., L. Mansinha, and R. P. Lowe, Localization of complex spectrum: The S-transform: IEEE Transactions on Signal Processing, 998-1001. , 1996