ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ...................................................3
ГЛАВА 1.
ОБЗОР РАБОТ НО МЕТОДАМ ИНТЕРПРЕТАЦИИ
ПОВЕРХНОСТНЫХ ГОДОГРАФОВ ОТРАЖЕННЫХ
ВОЛН И ТРЕХМЕРНОЙ МИГРАЦИИ................................6
1.1. Определение скорости но поверхностным и линейным годографам отраженных волн.............................7
1.2. Построение отражающих поверхностей..............16
1.3. Двумерная аппроксимация сейсморазведочных данных . . 28
ГЛАВА 2. V .
ПРОСТОЙ МЕТОД ТРЕХМЕРНОЙ КИНЕМАТИЧЕСКОЙ МИГРАЦИИ НА ОСНОВЕ СРЕДНЕСКОРОСТНОЙ МОДЕЛИ
СРЕДЫ....................................................32
2 I. Элементарное решение прямой и обратной задач кинематики нормальных лучей для однородной среды с плоской отражающей границей..................................33
2.2. Опробование алгоритма "МЮ-ЗЮ" на теоретических моделях..............................................42
2.3. Дифракция сферической волны на ребре двугранного
угла.............................................51
2.3.1 Кинематические и динамические характеристики
дифрагированных волн........................52
2.3.2. Восстановление пространственного положения
дифрагирующего ребра но алгоритму "МЮ-ЗО".... 63
ГЛАВА 3.
МЕТОД ТРЕХМЕРНОЙ КИНЕМАТИЧЕСКОЙ МИГРАЦИИ НА ОСНОВЕ ПЛАСТОВОЙ СКОРОСТНОЙ
МОДЕЛИ СРЕДЫ.............................................67
3.1. Геофизическая и математическая постановка задачи 67
2
3.2. Решение задачи трехмерной миграции.................70
3.2.1. Определение начальных условий................70
3.2.2. Учет преломления-отражения на границе раздела .. .72
3.2.3. Сглаживания и интерполяции на двумерных произвольных сетках..................................79
ГЛАВА 4
ИССЛЕДОВАНИЕ АЛГОРИТМА "MIG-3D PLAST".......................83
4.1. О выборе парамегров сглаживания....................83
4.2. Контроль точности расчета кинематических
параметров лучей....................................86
4.3. Опробование алгоритма "MIG-3D PLAST" на теоретических моделях....................................89
ГЛАВА 5
ПРИМЕНЕНИЕ РАЗРАБОТАННЫХ АЛГОРИТМОВ К РЕАЛЬНЫМ ДАННЫМ МОГТ............................98
5.1. Геологическая характеристика района работ..........98
5.2. Характеристика исходных экспериментальных материалов..............................................101
5.3. Построение структурных карт основных отражающих горизонтов.........................................104
5.3.1. Построение отражающих границ по алгоритму "MIG-3D"............................................111
5.3.2. Построение отражающих границ по алгоритму "MIG-3D PLAST”......................................118
ЗАКЛЮЧЕНИЕ.................................................120
ЛИТЕРАТУРА
121
Il
координат, в основе которого лежит аппроксимация наблюденного годографа гиперболическим годографом, который затем преобразуется в прямолинейный голограф Усреднение скоростного распределения в среде носит при этом интегральный характер (А1-Chalabi, 1984; Sain, Kaila, 1996).
К усреднению скоростного распределения в среде возможен и другой подход, реализованный в работе (Ампилов, Облогина, 1980), заключающийся в следующем. Вводится понятие локальной эффективной скорости V, определяемой по годографу отраженной волны по формуле V = -Jv'xft ,где х. I имеют смысл координат точки годографа, a v* - кажущаяся скорость в этой же точке. Данная формула позволяет вычислить эффективную (среднюю) скорость фактически для любого луча отраженной водны, приходящего в точку наблюдения, за исключением точки минимума годографа Показано, что в условиях слабонаклонного залегания отражающих границ определение скоростного закона среды с использованием локальной эффективной скорости более предпочтительно, чем по методу квадратичных координат.
Способы, основанные на пересчете кинематических параметров отраженных волн в скоростные параметры среды путем использования формул, связывающих значения пластовых скоростей со значениями предельных эффективных скоростей, предложены в работах (Урупов, Аккуратов,!975; Левин, 1975). Использование этих способов требует предварительного приведения интегральных значений v«-, к их предельным (дифференциальным) значениям. Такое приведение имеет целью учет преломления лучей на промежуточных границах, угла наклона границ и сейсмического сноса (Лозинский и др., 1974; Урупов, Левин, 1985). В работе (Левин, 1982) рассмотрены вопросы,
12
связанные с вычислением предельных значений эффективных скоростей для слоистых сред с криволинейными границами раздела. Показано, что величины предельных значений эффективных скоростей для годографов ОГТ определяются параметрами нормального сейсмического луча.
В работах (Бляс, 1982, 1987) предлагается метод определения пластовых скоростей горизонтально-слоистой среды, не требующий нахождения предельных эффективных параметров. В методе используются параметры гиперболических аппроксимаций годоірафов, которые находятся при анализе спектров скоростей
Способы, основанные на итеративном подборе параметров среды путем минимизации заданного функционала (различия рассчитанных теоретических годографов или их эффективных параметров с исходными) рассмотрены в работах (БаПк^ег, 1965; Гогонснков, Борейко, 1975). К рассматриваемой группе способов относится оценка параметров среды на основе метода оптимизации (Гольдин, 1979). Сущность способа заключается в отыскании такого вектора параметров среды, который минимизирует заданную целевую функцию.
Исходной информацией для способов послойного определения пластовых скоростей являются кривые и ц,, получаемые в результате интерпретации горизонтальных спектров скоростей и временных разрезов. Однако если годограф волны негиперболичен, скоростной спектр резко искажается и соответственно возникают ошибки в определении величины уогт. Эго приводит к ошибкам в вычисленных значениях пластовых скоростей. В работе (Мешбей и др., 1985) предложен алгоритм, в котором решение о соответствующей данному пласту (интервалу) скорости принимается в
- Київ+380960830922