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

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

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

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

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

АДИЕВ АЗАТ ЯВДАТОВИЧ

СОВЕРШЕНСТВОВАНИЕ СПОСОБОВ ОБРАБОТКИ ДАННЫХ ВЕРТИКАЛЬНОГО СЕЙСМИЧЕСКОГО ПРОФИЛИРОВАНИЯ ДЛЯ ИЗУЧЕНИЯ ОКОЛОСКВАЖИННОГО ПРОСТРАНСТВА

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

АВТОРЕФЕРАТ

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

2 9 МАР 1Ш

Уфа-2012

005012918

005012918

Работа выполнена в обществе с ограниченной ответственностью Научно-производственный центр «Геостра» Открытого акционерного общества «Башнефтегеофизика»

Научный руководитель - доктор геолого-минералогических наук, старший научный сотрудник Ленский Владимир Анатольевич

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

Валиуллин Рим Абдуллович - доктор технических наук, профессор, Башкирский государственный университет, кафедра «Геофизики», заведующий кафедрой;

Даниленко Виталий Никифорович - кандидат технических наук, ЗАО НПФ «ГИТАС», директор.

Ведущая организация - Открытое акционерное общество Научно-производственное предприятие «ВНИИГИС» (г.Октябрьский)

Защита состоится «13» апреля 2012 г. в 15-00 часов на заседании диссертационного совета Д 520.020.01 при открытом акционерном обществе «Научно-производственная фирма «Геофизика» (ОАО НПФ «Геофизика») по адресу: 450005, Республика Башкортостан, г. Уфа, ул. 8-ое Марта, 12.

С диссертацией можно ознакомиться в библиотеке ОАО НПФ «Геофизика»

Автореферат разослан «12» марта 2012 г.

Ученый секретарь '==^^г=^-=-__^Хисаева Дилара Ахатовна диссертационного совета

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

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

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

На начало постановки диссертационных исследований пакет программ "Оео8е1зРго" не обеспечивал полный цикл обработки данных ВСП и дополнялся пакетом программ обработки данных наземной сейсморазведки «ОеоС^ег-100». При этом использовались программные продукты, не в полной мере учитывающие специфику волновых полей на вертикальном профиле, что снижало качество результатов обработки и интерпретации. Кроме того, многократный переход из одного пакета в другой негативно влиял на технологичность обработки и, вследствие этого, на важнейшее преимущество скважинной сейсморазведки - экспрессность представления результатов.

Цель работы - совершенствование способов обработки и интерпретации данных ВСП в пакете "СсоЗспзРго" для повышения их экспрессности, информативности и достоверности.

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

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

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

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

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

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

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

5. Создание программ для обрабатывающего комплекса GeoSeis Pro на основе разработанных способов, их опробование и внедрение в процесс обработки и интерпретации данных методом ВСП в ОАО «Башнефтегеофизика».

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

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

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

2. Предложены приемы разделения волн по параметрам поляризации, повышающие качество выделения отраженных волн за счет подавления интен-

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

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

Основные защищаемые научные положения

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

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

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

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

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

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

Разработанные программы внедрены в технологический процесс обработки и интерпретации данных ВСП в ОАО «Башнефтегеофизика», что позволило существенно сократить сроки и одновременно повысить качество результатов обработки и интерпретации.

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

Апробация работы. Основные положения диссертационной работы докладывались на: 12-й международной геофизической конференции «Геомодель-2010» (Геленджик, 2010); научно технической конференции по разведочной геофизике в рамках нефтяного конгресса «Газ-Нефть» (Уфа, 2010); международном геофизическом симпозиуме «73rd EAGE Conference & Exhibition incorporating SPE EUROPEC 2011» (Вена, Австрия, 2011); XI ежегодной международной конференции и выставке «Гальперинские чтения - 2011» (Москва, 2011); «Продуктивные клиноморфные комплексы и возможности современной сейсморазведки» (Тюмень, 2011).

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

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

Диссертационная работа подготовлена в ООО НПЦ «Геостра» (ОАО «Башнефтегеофизика») под научным руководством доктора геолого-минералогических наук В.А. Ленского, которому автор искренне признателен.

Автор выражает благодарность кандидату технических наук Ю.Г. Антипину за постоянную методическую помощь в разработке алгоритмов. Автор считает приятным долгом выразить благодарность коллегам по работе, принимавшим участие в выполнении исследований и апробации разработанных программ: М.Т. Якупову, B.JI. Богдановой, Т.Н. Шаровой и другим.

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

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

В первой главе выполнены анализ и обобщение состояния обработки материалов вертикального сейсмического профилирования (ВСП). Приведен обзор отечественных и зарубежных программных комплексов обработки данных ВСП. В России наибольшее распространение получили программные комплексы ЗС-Interact (НПП «Гетек»), ВСП-ПОЛ (ОАО НПП «ВНИИГИС»), Univers (ООО «Геоверс»), Среди зарубежных известны комплексы VSProvess (Avalon Sciences Ltd), VISTA («GEDCO»), GeoCluster («CGG»), Q-BorSeis («Schlum-berger») и другие. Большинство этих комплексов ориентировано на использование определенных технологий. Они не позволяют в полной мере решать новые производственные задачи. Недостатком многих из них является неудобный пользовательский интерфейс, снижающий технологичность обработки. В связи с серьезными недостатками существующих обрабатывающих комплексов и учитывая большие объемы работ ВСП (около 60-70 скважин ежегодно) в ООО «Геостра» разрабатывается собственный комплекс программ обработки данных ВСП GeoSes Pro. Этот комплекс призван устранить ограничения существую-

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

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

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

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

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

геометрического расхождения фронта используется информация о скоростных характеристиках разреза v(t). По функции v(t) рассчитывают радиус фронта волны r(t) для каждого момента времени t, функция компенсации спада энергии за счет геометрического расхождения фронта имеет вид: c'(t) = r(t)=v(t)-t.

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

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

c(t)=t2 или c(t)=t15, предложенные Клаербоутом, или экспоненциальная функция

c-'(t) = v(t)-t-10bt,

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

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

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

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

Важнейшим вопросом при обработке данных ВСП является повышение временной разрешенности сейсмической записи. С целью сокращения длительности сигналов используется деконволюция. Решением теоретических и методических вопросов деконволюции сейсмических данных занимались Е.А.Робинсон, Г.Кюнец, Р.Райс, Е.А. Козлов, Ф.М. Гольцман, С.В.Гольдин, И.К. Кондратьев, А.А.Зенов, С.А. Кац, А.К. Яновский, W.R. Burus и др. Для реализации процедуры деконволюции в качестве входной информации необходимы данные о форме сейсмической волны, распространяющейся в геологической среде. Ввиду сложности получения такой информации из данных наземной сейсморазведки, вместо реального сигнала часто используется условный импульс, рассчитываемый по сейсмической трассе.

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

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

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

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

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

Во второй главе приведены исследования по разработке двух способов компенсации затухания сейсмических сигналов при ВСП. С целью ослабления влияния помех выполнена оптимизация затухания в способе, использующем непосредственно сейсмическую запись. В этом способе для расчета функции c(t), компенсирующей затухание сейсмической записи s0(ti), используется непосредственно реальная запись колебательного процесса s(t). Вычисление последовательности c(tK) выполняется по формуле:

i+wind/2

c(tk) = N-s0(t;)-wind/ 2 |s0(ti+k)|, (1)

i-wind/2

где wind- окно регулирования.

Из-за влияния помех поведение функции c(t) не соответствует истинному затуханию целевых волн и её применение приведет к искажению их динамических параметров. Для устранения искажений предложен следующий способ (рис.1).

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

А = 2>2/п,

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

где X; - значения трассы; п -число точек в интервале анализа. Затем каждое значение трассы приводится к заданному уровню усиления №

X; = N *Х, / А.

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

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

О = 3\У/4ЛЬ2= ШОЭ2А2, где XV - полная энергия колебаний; А - амплитуда волны; со - частота колебаний; Ь - радиус расхождения волнового фронта. Из приведенного соотношения видно, что амплитуда волны затухает как 1/Ь.

С целью учета преломления на сейсмических границах при компенсации геометрического расхождения вместо упрощенной функции Ь=г(Д)= автором предлагается рассчитывать реальную длину пути волн Ь на основе решения прямой кинематической задачи по заданной скоростной модели в виде пластов мощностью Ц и соответствующих им скоростей ^ (рис. 2).

Источник _*_

\ Ь \ V Ы. VI \

¡\ Приемник ./ \ 7

ЬЗ, УЗ \|/

Рис.2. Схема распространения сейсмического луча

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

этОи) / V! = этОг) / у2 =.. .= зт(1;) / ^ = р, где 1 - угол между лучом и нормалью к границе с номером]. Длина пути пробега волны в каждом пласте 1] равна: 1;= Ц / зш(^), где зш(^) = у>-р Полная длина пути Ь от источника до точки приема для каждого момента времени равна Ь = 1. Для вычисления параметра р использован способ, основанный на представлении длины пути Ь в виде ряда Тейлора ограниченной степени щ, предложенный Ю.Г. Антипиным:

Х-Аор-А1р3-...-Атр2т+1=0,

1-3-...- (2т-1) „ ,

где А0 = I Ь, ъ , Ат = - £( ^ у>+|),

2-4-2ш

X - проекция Ь на горизонтальную плоскость, равная проекции расстояния от источника до приемника.

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

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

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

В первом способе разделение волн выполняется во временной области по их поляризации. Трехкомпонентная запись преобразуется из первоначальной системы координат Х,У,2 в новую систему координат XI,У 1,71, которая повернута вокруг оси Ъ на угол а относительно первоначальной. Угол а определяется непосредственно по записям компонент X и У таким образом, чтобы координатная ось XI совпала с проекций полного вектора колебаний подавляемой волны на плоскость ХУ:

а = агсов( X / бяП (X2 + У2)).

Преобразование поворота системы координат Х,У,г в систему XI,У 1,21 осуществляются по известным формулам:

XI = X • соэ (а) + У • зт(а), (2)

У1 =-Х-БШ (а) + У-соз(а).

При этом траектория распространения волны при подходе к точке регистрации будет располагаться в вертикальной плоскости XI Ъ.

Далее производится сканирование плоскости XI Ъ посредством изменения угла ф наклона оси Ъ от 90 до 0 градусов с заданным шагом Лф и последующим вычислением компоненты волнового поля, сориентированной на соответствующий угол ф = 90 - ,)Дф. Вычисления производятся по записям волнового поля в каждой точке регистрации с использованием формул, аналогичных выше приведенным:

Х^Х^соБОЛф^г^БтуДф), (3)

Х^зт( ,|Лф) + ]Дф ).

Для каждой компоненты из полученного набора компонент вычисляется энергия в заданной временной области, содержащей волну-помеху, и выбирается та компонента, на которой энергия волны-помехи наименьшая. При этом соответствующая компонента Х1з будет сориентирована на направление, при котором выделяемые волны регистрируются с наибольшим отношением сигнал/помеха. На рис.4 приведены результаты обработки описанным способом для выделения отраженных продольных РР и обменных РБ волн с одновременным подавлением продольных и поперечных падающих волн соответствен-

Рис. 4. Выделение полей отраженных продольных РР(а) и обменных РБ (б) волн.

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

¡=п-1

Б^а) = 1/п X з,(1) [собГол) + 8т(оЛ)], ¡=0

где Б^со) - спектр одной трассы - б^); со - угловая частота; п - число отсчетов в трассе. Полученный результат в виде двумерной матрицы в координатах «частота - глубина» транспонируется и для каждой строки матрицы снова вычисляется одномерное преобразование Фурье. Окончательный результат вычислений показан на рис.5а в виде двумерного спектра волнового поля. На этой картине интерпретатором выделяется область помех, от которых необходимо очистить исходное волновое поле. Для примера на изображении двумерного спектра выделена область, соответствующая обменной отраженной РБ-волне, которая в данном случае принята за помеху. Далее в автоматическом режиме осуществляется перебор углов подхода волны и последующее преобразование волновых полей по формулам (3), каждый раз вычисляется двумерный спектр и его энергия в выделенной области. Из полученных последовательностей значений энергии и соответствующих этим значениям углов выбирается тот угол, при котором энергия в заданной области спектра минимальна (рис.5б). Затем, путем обратного двумерного преобразования Фурье, выбранная реализация волнового поля преобразуется во временную область.

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

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

точке регистрации ср^ который вычисляется путем решения прямой кинематической задачи по заданной скоростной модели разреза (рис.6, ось У перпендикулярна плоскости рисунка).

21 ^ г I, ^ 1

\ ---- \ ———■_____ V \ Ь х ' / \ 15 Г ш

\ 13 11 /! \ / ' \ / * \ / И

ТО

Рис. 6. Лучевая схема к расчету угла подхода отраженной волны

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

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

Преобразования осуществляются по формулам аналогично (3): XI = Х-соз(90°- Фз) + г-5'т(90°- <р,), Ъ\ =-Х-5т(90°- <рз) + г-соз(90° - (д ).

При этом учитывается, что угол ф; в одной и той же точке регистрации изменятся со временем ^ вследствие изменения глубины отражения.

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

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

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

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

ния Ц). Тогда сигнал s¡ на выходе системы определится известной формулой свертки последовательностей:

N

= Х Ц _; * и; .

¡=0

В частотной области спектр 5(со) результата свертки последовательностей Ц и и, выражается произведением соответствующих им спектров К(<м) и и(ю):

8(со) = К(<о) и(ш). (4)

При деконволюции ставится задача рассчитать импульсную характеристику системы Ц - то есть последовательность коэффициентов отражения, если заданы сейсмическая трасса б, и форма сейсмического импульса иг Из (4) определим частотную характеристику К(са) как:

К(со) = 5(со) • 1/и(ш).

Здесь выражение 1Ш(оз) определяет идеальную частотную характеристику обратного фильтра ц"'. Реальная импульсная характеристика обратного фильтра и;"' определяется как результат обратного преобразования Фурье от его частотной характеристики и"'(со):

и,"' = X и"'(ш) ' [а^сй) + шп(ш1)], где и-'(сй) = 1/и(со) = [а(со) - Ь(ш)] / [а2(со) + Ь2(со) + Я.2]; а(ш) и Ь(со) соответственно реальная и мнимая составляющие комплексного спектра и(ю); X - «отбеливающий» параметр, определяемый как доля от максимального значения в амплитудном спектре А(ш):

А(ш) = 5цП(а2(со) + Ь2(со)).

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

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

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

5 я ж № 0! я Я а » я т $т»шаятятят

-*-=-

4_н)

ГН 1

а б

Рис. 7. Данные расчета импульсной реакции обратного фильтра с нулевой (а) и оптимальной (б) задержкой 0- исходный сигнал; 1- импульсная реакция обратного фильтра; 2- результат свертки исходного сигнала с импульсной реакцией обратного фильтра.

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

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

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

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

Разработанные по теме диссертационной работы способы, рассмотренные в главах 2-4, реализованы в виде обрабатывающих программ в рамках

Рис.8. Исходное волновое поле (слева) и результат деконволюции (справа)

комплекса GeoSeis Pro. Программы прошли производственное опробование при обработке материалов ВСП и НВСП, полученных в различных регионах России более чем в 50 скважинах. Дополнение комплекса разработанными программами позволило с 2011 г. выполнять почти полный цикл обработки данных ВСП и НВСП в комплексе GeoSeis Pro, что привело к значительному повышению производительности и качества обработки. В результате за последние 4 года сроки подготовки отчета о результатах работ ВСП-НВСП были сокращены вдвое - с 8 недель до 4 недель при практически двойном увеличении объема получаемой и представляемой геофизической информации. Одновременно возросла точность результатов геологической интерпретации данных ВСП и НВСП. В качестве примера рассмотрены результаты применения усовершенствованной технологии обработки данных ВСП с использованием комплекса программ GeoSeis Pro в производственном режиме в двух скважинах, иллюстрирующие высокую информативность скважинной сейсморазведки.

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

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

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

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

5. Разработанные способы реализованы в виде алгоритмов и соответствующих обрабатывающих программ, опробованы на материалах ВСП и НВСП, полученных в различных регионах России более чем в 50 скважинах, и включены в состав комплекса обработки данных ВСП GeoSeis Pro, что позволило выполнять почти полный цикл обработки данных ВСП в этом комплексе и обеспечило значительное повышение производительности и качества обработки и, в конечном итоге, повышение информативности скважинной сейсморазведки. На основе полученных материалов подготовлены и внедрены в ОАО «Башнеф-тегеофизика» методические рекомендации по обработке данных ВСП и НВСП.

Основное содержание диссертации опубликовано в следующих научных трудах:

в изданиях, рекомендованных ВАК Минобрнауки РФ:

1. Ленский В .А. Патент № 2433426 Россия, МПК G01V 1/00. Способ определения азимутального направления трещиноватости пород / В.А. Ленский, А.Я. Адиев, P.A. Ахтямов, Е.В. Ленская // ООО НПЦ «Геостра»,-№2010112035/28; заявлено 29.03.2010; опубл. 10.11.2011; Бюл. №31.

2. Якупов М.Т. Патент № 2437124 Россия, МПК G01V 1/40. Способ сейсморазведки в криволинейных скважинах / М.Т. Якупов, P.A. Ахтямов, А.Я. Адиев, Ю.Г. Антипин // ООО НПЦ «Геостра»,- №2010122567/28; заявлено 02.06.2010; опубл. 20.12.2011; Бюл. №35.

3. Адиев А.Я. Развитие методики поляризационной обработки данных скважинной сейсморазведки /АЛ. Адиев, М.Т. Якупов // НТЖ «Оборудование и технологии для нефтегазового комплекса». - Москва: ВНИИОЭНГ. - 2011. -№ 6. -С.4-10.

в других изданиях:

4. Адиев АЛ. Возможность определения азимутального направления трещиноватости пород-коллекторов с целью выбора оптимального направления горизонтальных скважин на месторождениях Оренбургской области / А.Я. Адиев, A.B. Спиридонова, Е.В. Ленская // Сборник докладов второй научно-практической конференции «Разведочная геофизика в XXI веке». Выпуск 2. -Уфа: НПФ «Геофизика». - 2010. - С.37-40.

5. Адиев А.Я. Возможности архивной скважинной сейсморазведки при планировании геолого-технологических мероприятий / АЛ. Адиев, Д.А. Иванов, Е.В. Ленская // Сборник докладов второй научно-практической конференции «Разведочная геофизика в XXI веке». Выпуск 2. - Уфа: НПФ «Геофизика». -2010. -С.59-63.

6. Адиев АЛ. Анализ результативности работ НВСП в условиях Западной Сибири / АЛ. Адиев, В.А. Ленский, Е.В. Ленская и др. // Тезисы докладов 12-ой международной научно-практической конференции «Геомодель-2010»: Геленджик, 13-17 сентября 2010 [Электронный ресурс]: многопредмет. науч.

журн. /Тезисы и публикации конференций EAGE, Geomodel 2010. - Режим доступа к сборн.: http://www.earthdoc.org.

7. Адиев А.Я. Современное состояние и некоторые возможности повышения эффективности метода вертикального сейсмического профилирования / А.Я. Адиев, Ю.Г. Антипин // НТЖ «Нефть. Газ. Новации». - Самара: Агни. -2010. -№ 3. - С.26-31.

8. Иркабаев Д.Р. О применении метода ВСП, НВСП в районах с различной степенью изученности / Д.Р. Иркабаев, А.Я. Адиев, Т.Н. Шарова // Сборник докладов симпозиума «Тюмень-2011» «Продуктивные клиноформные комплексы и возможности современной сейсморазведки»: Тюмень, 28-29 марта 2011 [Электронный ресурс]: многопредмет. науч. журн. /Тезисы и публикации конференций EAGE, Turnen 2011. - Режим доступа к сборн.: http://www.earthdoc.ore.

9. Adiev A.Ya. Geological effectiveness of offset vertical seismic profiling at Orenburg oil fields / A.Ya. Adiev, A.A. Tikhonov // «73rd EAGE Conference & Exhibition incorporating SPE EUROPEC 2011». Vienna, 2011 23 - 26 May [Электронный ресурс]: многопредмет. науч. журн. /Тезисы и публикации конференций EAGE, Vienna 2011. - Режим доступа к сборн.: http://www.earthdoc.org.

10. Адиев А.Я. Разработка технологии компенсации затухания сейсмических сигналов / А.Я. Адиев, М.Т. Якупов // «Гальперинские чтения-2011». Материалы XI ежегодной международной конференции и выставки: Москва, 24-28 октября 2011 [Электронный ресурс]: многопредмет. науч. журн. /Тезисы и публикации конференций, Москва 2011. - Режим доступа к сборн.: http://www.geovers.ru.

11. Адиев А.Я. Развитие методики поляризационной обработки данных скважинной сейсморазведки / А.Я. Адиев, М.Т. Якупов // «Гальперинские чте-ния-2011». Материалы XI ежегодной международной конференции и выставки: Москва, 24-28 октября 2011 [Электронный ресурс]: многопредмет. науч. журн. /Тезисы и публикации конференций, Москва 2011. - Режим доступа к сборн.: http://www.geovers.ru.

Подписано к печати 07.03.2012. Формат 60x84 1/16. Бумага писчая. Гарнитура «Тайме». Усл. печ.л. 1,45. Уч-изд.л. 1,75. Тираж 100 экз. Заказ №7.

Отпечатано с готовых авторских оригиналов на ризографе в салоне полиграфии «СТЭФА» 450009, г.Уфа, пр.Октября, 12 тел: (347) 275-23-17

Содержание диссертации, кандидата технических наук, Адиев, Азат Явдатович

Введение.

1. Анализ состояния обработки материалов вертикального сейсмического профилирования.

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

1.2. Обзор исследований в области компенсации затухания сейсмических сигналов при ВСП.

1.3. Современное состояние поляризационной обработки данных ВСП.

1.4. Обзор способов деконволюции сейсмических записей при ВСП.

Выводы.

2. Разработка способов компенсации затухания сейсмических сигналов при ВСП.

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

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

2.3. Опробование разработанных способов компенсации затухания на материалах НВСП.

Выводы.

3. Развитие поляризационной обработки данных ВСП.

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

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

3.3. Выделение оптимальных компонент отраженных волн на основе пластовой скоростной модели разреза.

Выводы.

4. Разработка способа переменной во времени деконволюции по заданному набору сигналов.

4.1. Способ расчета обратного фильтра с оптимальной задержкой.

4.2. Способ переменной во времени деконволюции по заданному набору сигналов.

4.3. Опробования способа переменной во времени деконволюции на реальных материалах.

Выводы.

5. Результаты применения усовершенствованной обработки данных ВСП в производственных условиях.

5.1. Результаты ВСП в скважине 426-Бузулукская.

5.2. Результаты ВСП в скважине 1-Лагинская.

Выводы.

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

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

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

На начало постановки диссертационных исследований пакет программ «GeoSeis Pro» не обеспечивал полный цикл обработки данных ВСП и дополнялся пакетом программ обработки данных наземной сейсморазведки «GeoCluster-ЮО». При этом использовались программные продукты, не в полной мере учитывающие специфику волновых полей на вертикальном профиле, что снижало качество результатов обработки и интерпретации. Кроме того, многократный переход из одного пакета в другой негативно влиял на технологичность обработки и, вследствие этого, на важнейшее преимущество скважинной сейсморазведки - экспрессность представления результатов.

Цель работы - совершенствование способов обработки и интерпретации данных ВСП в пакете «GeoSeis Pro» для повышения их экспрессности, информативности и достоверности.

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

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

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

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

122 Выводы

Способы обработки, предложенные в разделах 2-4, реализованы в виде обрабатывающих программ, применяющихся в рамках комплекса обработки данных ВСП GeoSeis Pro при выполнении этапов поляризационной обработки, повышения разрешенности записей и корректировки (восстановления) амплитуд. Дополнение комплекса разработанными программами позволило ОАО «Башнефтегеофизика» с 2011 года выполнять почти полный цикл обработки данных ВСП и НВСП в комплексе GeoSeis Pro и усовершенствовать технологию обработки, за счет чего достигнуто значительное повышение производительности и качества обработки и в конечном итоге повышение качества интерпретации данных скважинной сейсморазведки.

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

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

123

ЗАКЛЮЧЕНИЕ

В связи с неудовлетворенностью существующими обрабатывающими комплексами и учитывая большие объемы работ ВСП (около 60-70 скважин ежегодно), в ОАО «Башнефтегеофизика» выполняется разработка собственного комплекса программ обработки данных ВСП GeoSeis Pro. Этот комплекс призван устранить ограничения существующих комплексов по технологичности обработки и интерпретации данных, недоучету специфики сейсмических наблюдений в скважинах и возможности решения новых геологических задач с использованием динамических параметров волн.

В процессе исследований по теме диссертационной работы автором получены следующие основные результаты:

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

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

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

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

5. Разработанные способы реализованы в виде алгоритмов и соответствующих обрабатывающих программ, опробованы на материалах ВСП и НВСП, полученных в различных регионах России более чем в 50 скважинах, и включены в состав комплекса обработки данных ВСП GeoSeis Pro, что позволило выполнять почти полный цикл обработки данных ВСП в этом комплексе и обеспечило значительное повышение производительности и качества обработки и, в конечном итоге, повышение информативности скважинной сейсморазведки. На основе полученных материалов подготовлены и внедрены в ОАО «Башнефтегеофизика» методические рекомендации по обработке данных ВСП и НВСП.

Библиография Диссертация по наукам о земле, кандидата технических наук, Адиев, Азат Явдатович, Уфа

1. Адиев А.Я. Современное состояние и некоторые возможности повышения эффективности метода вертикального сейсмического профилирования / А.Я. Адиев, Ю.Г. Антипин // НТЖ «Нефть. Газ. Новации». Самара: Агни. - 2010. - № 3. - С. 26-31.

2. Адиев А.Я. Развитие методики поляризационной обработки данных скважинной сейсморазведки / А.Я. Адиев, М.Т. Якупов // НТЖ «Оборудование и технологии для нефтегазового комплекса». Москва: ВНИИОЭНГ. -2011.-№ 6.-С. 4-10.

3. Адиев Р.Я. Программное обеспечение обработки и интерпретации данных скважинной сейсморазведки GeoSeis Pro / Р.Я. Адиев, Ю.Г. Антипин, В.В. Лесников // Свидетельство об официальной регистрации программы для ЭВМ №2003610472.

4. Агаев Х.Б. Алгоритм поляризационного анализа трехкомпонентных записей в сейсморазведке / Х.Б. Агаев, Е.И. Гальперин // Экспресс-информ., ВИЭМС. Сер. «Региональная, разведочная и промысловая геофизика», 1982.-Вып. 18.-С. 24-28.

5. Александров С.И. Методы поляризационной обработки сейсмических сигналов / С.И. Александров // Экспресс-информ., ВИЭМС. Сер. «Региональная, разведочная и промысловая геофизика», 1979. - Вып. 22. - 58 С.

6. Антипин Ю.Г. Система обработки данных вертикального сейсмического профилирования на ЭВМ (ОС ВСП) / Ю.Г. Антипин, Н. Р. Сивков // Библ. программ для обработки геофизических данных на ЭВМ. М.: ВНИИ-Геофизики. - 1975. - Вып. 32. - С. 54-62.

7. Антипин Ю.Г. Расчет профиля ВСП с учетом преломления волн на границах / Ю.Г. Антипин, P.C. Исмагилова // Разведочная геофизика. Вып. 102.-М.: Недра, 1986.-С. 88-93.

8. Антипин Ю.Г. Фазовый спектр сейсмического сигнала как индикатор уровня шума скважинного сейсмического прибора / Ю.Г. Антипин, С.Ю. Антипин // Геофизический вестник: Ежемес. информ. бюл. ЕАГО, 2003. № 3-4.-С. 37-41.

9. Бельфер И.К. Подсистема цифровой обработки материалов поляризационных сейсмических наблюдений (подсистема ПОЛЯР системы СЕИС-ПАК) / И.К. Бельфер, И.В. Сапункова. М.: ГФАП, 1979. - 248 С.

10. Бреховских Л.М. Волны в слоистых средах / Л.М. Бреховских. М., 1957.-342 С.

11. Быков И.А. Алгоритмы цифровой обработки скважинных азимутальных наблюдений / И.А. Быков, И.М. Тихонова // Вопросы динамической теории распространения сейсмических волн. Л.: Наука, 1980. - Вып. 20. -С.135-155.

12. Быков И.А. Прослеживание сейсмических волн при трехкомпо-нентных скважинных наблюдениях / И.А. Быков // Вопросы динамической теории распространения сейсмических волн. JL: Наука, 1982. - Вып. 22. -С. 150-162.

13. Быков И.А. Цифровая обработка трехкомпонентных скважинных сейсмических наблюдений / И.А. Быков, И.М. Тихонова // Скважинные сейсмические и акустические исследования в рудных районах. Д.: Наука, 1982.-С. 19-34.

14. Бондарев В.И. Сейсморазведка. Учебник для студентов высших учебных заведений / В.И. Бондарев. Информационно-издательский центр. -Екатеринбург, 2008. - 703 С.

15. Ван Юйчао Метод обработки данных 3D ВСП / Юйчао Ван // Галь-перинские чтения 2007, VIII Ежегодная международная конференция «ВСП и трехмерные системы наблюдений в сейсморазведке». - С. 71-72.

16. Гальперин Е.И. Вертикальное сейсмическое профилирование. 2-е изд., доп. и перераб. / Е.И. Гальперин. - М.: Недра, 1982. - 344 С.

17. Гальперин Е.И. Результаты применения поляризационного метода вертикального сейсмического профилирования (ПМ ВСП) на территории

18. Краснодарского края / Е.И. Гальперин, Ю.Д. Мирзоян, В.Я. Ойфа // Экспресс-информация. Регион, развед. и промысловая геофизика. М.: ВИЭМС, 1980.-Вып. 16.-С. 6-15.

19. Гальперин Е.И. Поляризационный метод сейсмических исследований / Е.И. Гальперин. М.: Недра, 1977. - С. 276.

20. Гальперин Е.И. Об одном свойстве поляризационных колебаний в сейсмических волнах / Е.И. Гальперин, A.B. Фролова, P.M. Гальперина // ДАН СССР, 1982. Т. 264. - № 3. - С. 596-600.

21. Гальперин Е.И. Поляризационный метод основной метод сейсморазведки месторождений твердых полезных ископаемых / Е.И. Гальперин, JI.A. Певзнер // XXX Международный геофизический симпозиум. - М., 1985. - (1). - С. 80-89.

22. Гальперин Е.И. Цифровая обработка материалов поляризационного метода сейсморазведки / Е.И. Гальперин // Нефтяная геофизика: Сб. докладов. -М.: СЭВ, 1982. Т. 1 (Сейсморазведка). - С. 248-253.

23. Гамбурцев Г.А. О возможности выделения поперечных отраженных волн при сейсморазведке по методу отражений / Г.А. Гамбурцев. Изв. АН СССР. - Сер. геоф. № 2, 1939. - С. 155-173.

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

25. Гольцман Ф.М. Статистические алгоритмы восстановления сигналов и полей по их линейным преобразованиям / Ф.М. Гольцман М.: Изв. АН СССР. - Сер. «Физика Земли», № 4. - 432 С.

26. Горелик Г.С. Колебания и волны / Г.С. Горелик. 2 изд. — М.: Наука, 1959.-281 С.

27. Гурвич И.И. Сейсмическая разведка / И.И. Гурвич. М.: Гостоп-техиздат, 1960. - 504 С.

28. Зенов A.A. Самонастраивающийся цифровой фильтр для повышения разрешенности сигналов на сейсмической трассе / A.A. Зенов // Изв. высш. учеб. заведений. Геология и разведка. -М., 1973. № 12. - С. 67-74.

29. Кац С.А. Синтез оптимальных обратных фильтров для сигналов, искаженных аддитивными и мультипликативными помехами / С.А. Кац // Интерпретация и обнаружение сейсмических волн в неоднородных средах. -М.: Недра, 1971. С. 27-33.

30. Кац С.А. Корректирующий фильтр для автоматической увязки сейсмограмм / С.А. Кац // Сейсмические волны в тонкослоистых срезах.

31. АН СССР, Институт физики Земли им. О.Ю. Шмидта. М.: Наука, 1973. - С. 42-47.

32. Канасевич Э.Р. Анализ временных последовательностей в геофизике / Э.Р. Канасевич. М.: Недра, 1985.-401 С.

33. Клаербоут Д.Ф. Сейсмическое изображение земных недр. Пер. с англ., редактор перевода Потапов O.A. М.: Недра, 1989. - 407 С.

34. Козлов Е.А. Цифровая обработка сейсмических данных / Е.А. Козлов, Г.Н. Гогоненков, Б.Л. Лернер и др. М.: Недра, 1973. - 312 С.

35. Козлов Е.А. Об обратной фильтрации записей MOB / Е.А. Козлов // Прикладная геофизика. Вып. 51. - М.: Недра, 1968. - С. 45-52.

36. Козлов Е.А. Об эффективности обратной фильтрации записей MOB / Е.А. Козлов, P.A. Калинина, Т.С. Кедрова и др. // Прикладная геофизика. -Вып. 30. -М.: Недра, 1968. С. 37^4.

37. Кондратьев И.К. Корректность и помехоустойчивость способов обратной фильтрации сейсмических трасс / И.К. Кондратьев // Прикладная геофизика». Вып. 50. - М.: Недра, 1967. - С. 3-22.

38. Кондратьев И.К. К теории цифровой обратной фильтрации волновых процессов / И.К. Кондратьев // Изв. АН СССР. Сер. «Физика Земли», № 4. -М., 1968.-С. 22^5.

39. Кондратьев И.К. Использование обратной и веерной фильтрации в методе ОГТ / И.К. Кондратьев, A.A. Зенов, Н.И. Крылова и др. // Прикладная геофизика. Вып. 68. -М.: Недра, 1972. - С. 80-85.

40. Кондратьев И.К. О выборе формы сигнала для обратной фильтрации сейсмических записей / И.К. Кондратьев // Изв. АН СССР. Сер. «Физика Земли», № 10. - М., 1969. - С. 47-60.

41. Кондратьев И.К. Оценка искажений результатов обратной фильтрации сейсмических трасс в реальных условиях / И.К. Кондратьев, Л.И. Мор-жина // Прикладная геофизика. Вып. 60. - М.: Недра, 1970. - С. 58-79.

42. Кравчук A.A. Исследование некоторых алгоритмов деконволюции сейсмических материалов / A.A. Кравчук // Сейсмические методы поиска и разведки полезных ископаемых. Киев: Техника, 1974. - С. 38^49.

43. Ленский В.А. Пакет программ обработки и интерпретации данных ВСП для полевых аппаратурно-методических комплексов / В.А. Ленский, А.Л. Ишбулатова // НТВ «Каротажник». Тверь: АИС, 1997. - Вып. 30. - С. 56-60.

44. Патент № 2433426 Россия, МПК G01V 1/00. Способ определения азимутального направления трещиноватости пород / В.А. Ленский, А.Я. Адиев, P.A. Ахтямов, Е.В. Ленская; ООО «НПЦ «Геостра». № 2010112035/28; заявлено 29.03.2010; опубл. 10.11.2011. Бюл. №31.

45. Лобусев М.А. Методика построения скоростной модели среды в методе ВСП на основе использования данных о поляризации сейсмических волн. Дис. канд. техн. наук. Москва, 2006 г. 119 С.

46. Малкин А.Л. Способы оценивания временной разрешенности сейсмической записи / А.Л. Малкин, Д.Б. Фиников // Геология и геофизика, 1986. -№ 12. С. 54-59.

47. Мирзоян Ю.Д. Опыт и результаты применения ПМ ВСП в сложно-построенных средах / Ю.Д. Мирзоян, С.Э. Камбарли, В.Я. Ойфа, Э.Я. Куцен-ко // Геология, геофизика и разработка нефтяных и газовых месторождений. -2003.-№ 11.-С. 24-38.

48. Михальцев A.B. Применение обратной фильтрации при цифровой обработке данных сейсморазведки ОГГ / A.B. Михальцев, И.А. Липовецкий // Экспресс-информ. ВИЭМС. Сер. IX «Per. разв. и промысл, геофизика». -Вып. 10, 1974.-С. 43-45.

49. Редекоп В.А. Применение программного комплекса VimSeis ЗБ для обработки данных ВСП / В.А. Редекоп, В.В. Помазанов, Д.Г. Риле // Материалы научно-практической конференции «ВСП и трехмерные системы наблюдений в сейсморазведке». М., 2004. - С. 86-87.

50. Рыжков В.И. Пакет интерактивной обработки трехкомпонентных многоволновых данных ВСП ЗС-ШТЕЯАСТ / В.И. Рыжков, В.М. Кузнецов // Научно-производственное предприятие ГЕТЭК. 118 С.

51. Робинсон Э.А. Статистические методы сжатия сейсмического импульса / Э.А. Робинсон. ТИИЭР, 1984. - Т. 72,- №10.-214 С.

52. Сильвия М.Т. Обратная фильтрация геофизических временных рядов при разведке на нефть и газ / М.Т. Сильвия, Э.А. Робинсон. М.: Недра, 1973.-413 С.

53. Табаков Г.Г. Система обработки данных сейсморазведки для ЭВМ ЕС-СЦС-3 / Г.Г. Табаков, В.И. Афанасенко, В.Б. Маневич // Методика и техника сейсморазведки. М.: Недра. - 1979. - С. 82-85.

54. Тимошин Ю.В. Совместное применение обратной фильтрации и интерференционных систем / Ю.В. Тимошин // Сейсмич. методы поиска и разведки полезных ископаемых. Киев, 1973. - С. 33-38.

55. Худзинский Л.Л. Комплексная программа цифровой обработки данных скважинной сейсморазведки МОВ / Л.Л. Худзинский, Г.Е. Руденко // Изв. вузов. Геология и разведка. М., 1976. - № 6. - С. 119-123.

56. Хургин Я.Н. Методы теории целых функций в радиофизике, теории связи и оптике / Я.Н. Хургин, В.П. Яковлев. М.: Физматгиз, 1962. - 220 С.

57. Шевченко A.A. Скважинная сейсморазведка. / A.A. Шевченко. M.: РГУ нефти и газа, 2002. - 129 с.

58. Шехтман Г.А. Изучение структуры околоскважинного пространства по разностным и дифференциальным вертикальным годографам / Г.А. Шехтман, В.В. Попов // Технология сейсморазведки, 3. 2007. - С. 39-43.

59. Патент № 2437124 Россия, МПК G01V 1/40. Способ сейсморазведки в криволинейных скважинах / М.Т. Якупов, Р.А. Ахтямов, А.Я. Адиев, Ю.Г. Антипин; ООО «НПЦ «Геостра». № 2010122567/28; заявлено 02.06.2010; опубл. 20.12.2011. - Бюл. № 35.

60. Bois P. Analyse sequentielle / P. Bois // Geophys.Prosp., v. 20, № 3, 1972.-97 P.

61. Bois P. Determination de l'impulsion sismique (wave let impulses) / P. Bois // Geophys. Prosp, v. XVI, № 1,1968. 514 P.

62. Burns W.R. A statistically optimized deconvolution. / W.R. Burns // Geophysics, v. 33, №2,1968.-232 P.

63. Hardage B.A. Vertical seismic profiling principles (2nd ed.) / B.A. Hardage. - Amsterdam, Geophysical Press, 1985. - 509 P.

64. Deregowski S.M. Optimum digitale filtering and inverse filtering in the frequensy domain. / S.M. Deregowski // Geophys. Prosp., v. 19, № 4,1971. 391 P.

65. D'Oeraene J. A self-adoptiny deconvolution and deringing program / J. D'Oeraene, Wlodarczak R. Protee // Geophysical Prosp., v. 16, № 3,1968. 207 P.

66. Dillon W.G. The effect of noise on iterative VSP inversion / W.G. Dillon, T.W. Spencer // 54th Annual International Meeting of SEG. Exp. abstracts, 1984.-P. 564-569.

67. Geocluster, Seismic Processing Software, Release 4.1 May 2006. Compagnie Generale de Geophysigue 1, Rue Leon Migaux 91341 Massy Cedex-France. -918 P.

68. Kunetz G. Essai d'analyse da traees sismiques / G. Kunetz // Geophys. Prosp., v. IX., № 3, 1961. 201 P.

69. Kunetz G. Efficient deconvolution of marine seismic records / G. Kunetz, J. Fourmann // Geophys. Prosp., v. 33, 1968. 194 P.

70. Robinson E.A. Predictive deconvolution of seismic traces / E.A. Robinson // Geophysies, v. 22, № 4, 1957. 178 P.

71. Robinson E.A. Predictive deconvolutien of time series with application to seismic exploration / E.A. Robinson // Geophysics, № 3, 1967. 151 P.

72. Rice R. Invers convolution filters / R. Rice // Geophysics, v. 27, № 2, 1962.- 196 P.

73. Shrwood J. Minimum-phase and related properties of the response of a horizontally stratified absorptive earth to plane acoustic waves (Deconvolution transform) / J. Shrwood // Geophusice, № 2, 1965. 134 P.