Obrabotka_i_interpretatsia_dannykh_seysmorazvedki
.pdfДля определения статических поправок используются годографы первых вступлений полевых сейсмограммам ОГТ или проводятся специальные работы по изучению верхней части разреза. При расчёте принимается, что волны проходят область выше уровня приведения по вертикали, что является хорошим приближением при небольшой скорости в ЗМС и достаточной глубине отражающих границ.
2. Модель ВЧР
Рисунок 12 Модель ВЧР для определения статических поправок
z, м
|
tВ |
hВ |
|
|
|
hР hЗ |
VЗ |
|
|
|
VП |
|
|
Линия приведения |
|
|
|
|
х, м |
Общие параметры |
|
||
hЗ |
– |
глубина подошвы ЗМС; |
|
hР |
– |
глубина линии приведения; |
|
VЗ |
– |
скорость в ЗМС; |
|
VП |
– |
скорость в слое между подошвой ЗМС и линией |
|
приведения. |
|
|
|
Для взрывной сейсморазведки |
|
||
hВ |
– |
глубина заложения заряда ВВ; |
|
tВ |
– |
время пробега |
волны от поверхности до уровня |
заложения ВВ по вертикали.
3. Формулы расчёта статических поправок
Статическая поправка сейсмической трассы складывается из статических поправок за пункт приёма и пункт возбуждения трассы:
СтП СтППП СтППВ
11
Суммарная статическая поправка вычитается из времени вступления всех волн соответствующей трассы:
|
t' t |
СтП |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Таблица 3 Формулы расчёта статических поправок |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Вид поправки |
|
|
|
|
Формула |
|
|
|
|
|
|
||||||||||
|
За пункт приема |
СтП |
|
|
|
|
hЗ |
|
|
hР hЗ |
|
|
||||||||||
|
|
ПП |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
VЗ |
Vп |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
За пункт возбуждения |
СтП |
|
|
|
hЗ |
|
hР hЗ |
|
|
|
|
||||||||||
|
hВ = 0, источник на |
|
|
|
|
|
||||||||||||||||
|
ПВ |
VЗ |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
VП |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
поверхности земли |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
За пункт возбуждения |
СтП |
t |
|
hЗ hВ |
|
hР hЗ |
|
||||||||||||||
|
hВ < hЗ, забой скважины |
|
||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||
|
ПВ |
|
|
|
В |
|
|
|
|
VЗ |
|
VП |
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
не достигает подошвы |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ЗМС |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
За пункт возбуждения |
СтП |
|
|
t |
|
|
hР hВ |
|
|
|
|||||||||||
|
hВ ≥ hЗ, забой скважины |
|
|
|
|
|||||||||||||||||
|
ПВ |
В |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
VП |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
пересекает подошву |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ЗМС |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дано
По результатам проведения микросейсмокаротажа были определены параметры ВЧР с шагом 500 метров на профиле длиной 10000 м (таб 5). Верхняя часть разреза характеризуется превышением рельефа hP, мощностью ЗМС hЗ, глубиной заложения взрывных скважин hВ и вертикальным временем tВ (для источников). Зона малых скоростей и подстилающие породы имеют скорости VЗ и VП соответственно.
На указанном профиле проводились работы методов общей глубинной точки с кратностью 12 по системе наблюдений с параметрами, указанными в таблице 6.
Таблица 4 Результаты определения параметров ВЧР по данным микросейсмокаротажа
Пикет |
hР, м |
hЗ, м |
VЗ, м/с |
VП, м/с |
hВ, м |
tВ, мс |
000 |
50.0 |
20.00 |
700 |
2000 |
15 |
21 |
500 |
47.0 |
19.00 |
800 |
2100 |
13 |
16 |
1000 |
42.0 |
16.00 |
600 |
2200 |
10 |
17 |
1500 |
40.0 |
15.50 |
400 |
1900 |
12 |
30 |
2000 |
37.0 |
13.00 |
450 |
1800 |
15 |
30 |
2500 |
35.0 |
11.50 |
600 |
1700 |
17 |
22 |
3000 |
33.0 |
10.00 |
650 |
1600 |
16 |
19 |
3500 |
32.0 |
9.50 |
700 |
2100 |
12 |
15 |
|
|
|
12 |
|
|
|
4000 |
30.0 |
8.00 |
850 |
2200 |
10 |
10 |
4500 |
35.0 |
14.00 |
600 |
1900 |
16 |
24 |
5000 |
40.0 |
19.00 |
650 |
1800 |
16 |
25 |
|
|
|
|
|
|
|
5500 |
42.0 |
21.50 |
600 |
2100 |
17 |
28 |
|
|
|
|
|
|
|
6000 |
44.0 |
23.00 |
450 |
2300 |
18 |
40 |
6500 |
44.0 |
22.50 |
400 |
2200 |
19 |
48 |
|
|
|
|
|
|
|
7000 |
42.0 |
20.00 |
700 |
2100 |
18 |
26 |
|
|
|
|
|
|
|
7500 |
40.0 |
16.00 |
750 |
1950 |
15 |
20 |
8000 |
35.0 |
9.00 |
850 |
2000 |
17 |
15 |
8500 |
35.0 |
8.00 |
900 |
2100 |
19 |
14 |
9000 |
33.0 |
5.00 |
800 |
2400 |
10 |
8 |
9500 |
37.0 |
8.00 |
600 |
2600 |
14 |
16 |
10000 |
40.0 |
10.00 |
450 |
2700 |
16 |
24 |
Задание
1.Произвести выборку 12-ти ПВ и ПП для указанной в задании ОГТ.
2.Построить модель ВЧР в пределах выбранных ПВ и ПП
3.Рассчитать и построить на графике суммарные статические поправки
Таблица 5 Варианты задания
Вар |
1 |
|
2 |
3 |
|
4 |
5 |
|
6 |
7 |
|
8 |
9 |
|
10 |
11 |
|
12 |
ОГТ |
1400 |
|
2000 |
2600 |
|
3200 |
3800 |
|
4400 |
5000 |
|
5600 |
6200 |
|
6800 |
7400 |
|
8000 |
Система |
Фланговая прямая |
|
Центральная |
|
|
Центральная |
|
|||||||||||
наблюд |
|
dXПП = 50 м |
|
симметричная |
асимметричная |
|||||||||||||
ений |
|
|
||||||||||||||||
s = 24; d = 1; dR = 1; |
|
dXПП = 100 м |
|
|
dXПП = 100 м |
|
||||||||||||
|
|
|
|
|
||||||||||||||
|
|
|
SПВ =1 |
|
s = 24; d = 1; dR = 0; |
s = 24; d = 1; dR = 0; |
||||||||||||
|
|
|
|
|
|
|
|
|
SПВ =12 |
|
|
|
SПВ =6 |
|
Общая длина профиля – 10000 м
Ход решения
1.Введите параметры системы наблюдений в программу Geometry и выберите из файла SHOT_CDP в листе CDP 12 пар источник/приемник: X_shot1 и X_rec1; X_shot2 и X_rec2 и т.д. для ОГТ(CDP) вашего варианта. Данные занесите в таблицу 6:
Таблица 6 Пары источник/приемник для заданной ОГТ
Номер |
ХПВ |
ХПП |
1 |
|
|
2 |
|
|
3 |
|
|
… |
|
|
N |
|
|
|
13 |
|
2.Для части профиля, на которой располагаются пункты возбуждения и пункты приема составленной выборки ОГТ установите параметры ВЧР (см таб 5).
3.Постройте по данным таблицы 5:
Глубинный разрез ВЧР – зависимости hP, hЗ, hВ от координат пикетов
Скоростную модель – зависимости VП, VЗ от координат пикетов Изменение tВ по профилю – зависимость tВ от координаты
пикетов (рис 13)
4.Определите по построенным графикам при помощи графической интерполяции параметры ВЧР для всех пунктов возбуждения и пунктов приёма.
5. Рассчитайте суммарные статические поправки |
для пар |
источник/приемник по таблице 6 |
|
СтП СтППП СтППВ
6.Постройте график зависимости суммарных статических поправок от удаления
СтП f ( Х ПП Х ПВ )
Контрольные вопросы
1.Что такое статические поправки?
2.Статические поправки за пункты возбуждения, приема и суммарные.
3.Верхняя часть разреза и ее параметры.
4.Способы определения статических поправок.
14
Рисунок 13 Графики для определения параметров ВЧР
Глубинный разрез
z, м
hВ
hР
hР - hЗ
Линия приведения
V, м Скоростная модель
VП
VЗ
t, мc Изменение tВ по профилю
tВ
х, м
15
3. Скоростной анализ и кинематические поправки
Введение
1.Скоростной закон
Вмногослойной среде годографы отраженных волн - сложные
зависимости, приближенно соответствующие гиперболами с
эффективными параметрами - t0 ,VОГТ.
t |
x2 |
t 2 |
|
2 |
|||
|
0 |
||
|
VОГТ |
|
Это выражение выводится для горизонтальной отражающей границы и однородной покрывающей толщи достаточно точно описывает реальные отражения. Определяемая по годографу ОГТ скорость является фиктивной - зависящей от эффективной скорости покрывающей толщи V и угла наклона границы φ.
VОГТ V cos
Зависимость скорости VОГТ полезных отражений от t0 называется вертикальным спектром или скоростным законом. Определение скоростного закона - скоростной анализ - происходит путем подбора расчетных годографов. На исходных сейсмограммах ОГТ производится перебор по гиперболам, кривизна которых зависит от последовательно увеличиваемых параметров VОГТ и t0. Эффективные свойства отраженных волн находятся при достаточном сходстве одной из расчетных гипербол и наблюдаемой волны.
Рисунок 14 Определение параметров отражения по годографу ОГТ
2. Модель пластовых скоростей |
|
|
|
|
Для |
горизонтально-слоистой |
среды |
VОГТ |
является |
удовлетворительной оценкой предельной эффективной скорости VЕ. |
||||
Используя |
VОГТ вместо предельной |
эффективной |
скорости, можно |
рассчитать модель пластовых скоростей по формуле Урупова–Дикса:
16
|
|
V |
ПЛ |
m |
V ОГТ m |
2t |
0m |
V ОГТ m 12t |
0m 1 |
|||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
t0m t0m 1 |
|
|||
ПЛ |
|
– |
|
пластовая скорость m-го слоя |
|
|||||||||
V |
m |
|
|
|||||||||||
VОГТm |
– |
|
фиктивная скорость m–1 го слоя |
|||||||||||
t0 m |
|
– |
|
время t0 |
для m–го слоя |
|
|
|||||||
VОГТ m-1 |
– |
|
фиктивная скорость m–1 го слоя |
|||||||||||
t0 m-1 |
– |
|
время t0 |
для m-1–го слоя |
|
|||||||||
Для первого слоя t0 m-1 = 0; VПЛ |
1 = VОГТ1 |
|
|
|||||||||||
Мощность слоев hm находится по формуле |
|
|||||||||||||
|
|
|
|
|
h m |
|
t0m 1 |
t0m |
VПЛm |
|
||||
|
|
|
|
|
|
|
|
2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3. Выделение кратных волн
Отраженные волны разделяются на полезные и кратные. Полезные волны образуются в результате однократного отражения от целевых границ раздела. Кратные волны несколько раз отражаются от границ изучаемой толщи. Полезные волны содержат информацию о строении геологической среды, а кратные таких данных не несут и могут представлять серьезную помеху для полезных волн и выделяться как ложные границы на временном разрезе.
Интенсивные кратные волны образуются как правило в верхней части разреза, где контрастность акустических свойств наиболее велика. Фиктивная скорость таких волн невелика, а кривизна годографа достигает больших значений. Полезные волны формируются на более глубоких границах раздела и потому имеют большую скорость. Скоростные параметры полезных волн в условиях осадочного чехла постепенно растут с глубиной, так как отражают объективные свойства изучаемого геологического разреза.
Если при примерно одинаковом t0 наблюдаются две волны, то кратные волны выделяются по большей крутизне годографа. При скоростном анализе это различие выглядит как локальное уменьшение скорости на отдельных участках сейсмограммы, когда как в целом фиктивная скорость с увеличением времени возрастает.
4. Кинематические поправки
Кинематические поправки исправляют кривизну осей синфазности полезных отраженных волн для последующего суммирования и приведения к форме отражающих границ.
17
В результате введения кинематических поправок оси синфазности полезных отраженных волн распрямляются, а кратных – за счёт их меньшей скорости и глубины кратнообразующих границ – остаются недоспрямленными. Исправленные сейсмограммы суммируются по способу ОГТ, что приводит к выделению полезных отраженных волн на фоне ослабленных кратных отражений.
Кинематические поправки могут быть введены только после учета влияния верхней части разреза и проведения скоростного анализа. Вместе со статикой кинематика при обработке данных вводится несколько раз, с уточнением значений поправок на каждом шаге.
Дано
Сейсмограмма ОГТ, рассчитанная по случайной модели, состоящей из 6 слоев и 5-и горизонтальных отражающих границ. Сейсмограмма составлена из N = 12 трасс, расстояние между ОГТ dXОГТ = 100 м, вынос пункта возбуждения не применяется R = 0.
Задание
Опередите закон изменения скорости VОГТ и по полученным данным рассчитайте модель пластовых скоростей. Для 3-ей отражающей границы определите таблицу кинематические поправок.
Ход решения
1. Проведите моделирование
Создайте случайную модель в программе Velocity и выведите на экран сейсмограмму ОГТ.
2. Определите скоростной закон
Выделите на сейсмограмме полезные отражения по признаку монотонного возрастания фиктивной скорости и определите данные скоростного закона. Полученные результаты занесите в таблицу 7 и постройте на графике (рис 15).
Таблица 7 Таблица скоростного закона
Номер границы |
t0 |
VОГТ |
1 |
|
|
….. |
….. |
….. |
18
Рисунок 15 Скоростной закон
1 VОГТ, м/с
2
3
4
t0, мс 5
3. Определите модель пластовых скоростей
Используя данные скоростного закона и формулы 1 и 2, заполните таблицу модели пластовых скоростей. Средняя скорость определяется через пластовые для подошвы каждого слоя. Диаграмму пластовых скоростей постройте на графике.
Таблица 8 Модель пластовых скоростей
Номер слоя |
h, м |
VОГТ, м/с |
VПЛ, м/с |
VCР, м/с |
m |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
… |
… |
… |
… |
… |
|
|
|
|
|
Рисунок 16 Модель пластовых скоростей
4. Рассчитайте кинематические поправки
Для третьей отражающей границы вычислите значения кинематических поправок N трасс по формулам:
19
|
|
|
l |
|
|
2 |
|
2 |
|
|||
КнП |
|
|
|
|
|
|
|
t0 |
t0 (3) |
|||
|
|
ОГТ |
|
|
||||||||
|
V |
|
|
3 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
КнП |
|
|
|
l |
|
|
2 |
t |
2 |
t0 (4) |
||
|
|
|
|
|
|
|
|
0 |
||||
|
|
СР |
|
|
||||||||
|
|
V |
|
3 |
|
|
|
|
|
В результате получится таблица с двумя наборами кинематических поправок. Результаты занесите в таблицу и сравните результаты :
Таблица 9 Рассчитанные кинематические поправки
Трасса |
l, м |
КнП (по VОГТ), |
КнП (по VСР), |
|
|
мс |
мс |
|
|
|
|
1 |
0 |
|
|
… |
… |
…. |
…. |
12 |
1100 |
|
|
Постройте на одном графике зависимости КнП (по VОГТ) и КнП (по VСР) от удаления l.
Контрольные вопросы
1.Фиктивная скорость VОГТ
2.Скоростной закон и его определение
3.Выделение кратных волн на скоростном законе.
4.Модель пластовых скоростей.
5.Цель ввода кинематических поправок.
20