Расчет спектров пластовых скоростей по отраженным волнам.

 

Этот вид работы возможен только в случаях:

- горизонты были пропикированы по немигрированному временному разрезу.

- горизонты не содержат «петель», т.е. зависимость t0(x) является однозначной функцией.

 

В случаях с петлями следует сперва выполнить процедуры временной миграции, затем пропикировать горизонты по мигрированному разрезу (без петель), затем получить априорную слоистую модель через пункт меню Depth Migr->Ray Tracing from Migrated Horizon и приступить к уточнению модели пластовых скоростей по дифрагированным волнам.

 

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

 

В диалоговом окне Calc Horizon Velocity Spectrum следует указать:

 

Horizon Number – текущий номер горизонта.

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

Если будет выбран горизонт из ранее обработанных, то, при успешном завершении работы, удаляются из проекта все результаты обработки этого горизонта.

Start Velocity – начальное значение пластовой скорости (м/сек).

Stop Velocity – конечное значение пластовой скорости (м/сек).

Velocity Step – приращение значений пластовой скорости (м/сек).

Time Window – величина временного окна в мсек.

CMP Base – величина базы анализа в точках ОГТ.

Processing From Alop – при включенном флажке предполагается,

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

Start CMP picket – Начальный пикет ОГТ в метрах для обработки и построения модели.

Целое число.

Stop CMP picket – Конечный пикет ОГТ в метрах для обработки и построения модели.

Целое число.

Distance step of Model – Шаг построения модели в метрах (число с плавающей точкой).

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

Modify Current Delta T0 – это очень важная кнопка.

Без нажатия этой кнопки Вы не сможете продолжить работу.

Эта кнопка позволяет подправить (сгладить) производные T0 по текущему горизонту.

 

Результаты будут сохранены в файлах:

Output Spectrum File Name  - Файл для записи скоростных спектров.

Имя файла определяется по имени файла задания. Имя формируется добавление к этому имени строки _HORxxx, где xxx - порядковый номер горизонта (000, 001, ...) Расширение файла SVV. Имя и расширение файла не может быть изменено.

Output Horizon File Name - Файл для записи всех определяемых данных по текущему горизонту.

Такими данными являются:

- скорости ОГТ (по первому горизонту).

- значения T0;

- значения производных линии T0;

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

- углы наклона границы;

- X,Z координаты точек границы.

Имя файла определяется по имени файла задания. Имя формируется добавление к этому имени строки _HORxxx, где xxx - порядковый номер горизонта (000, 001, ...) Расширение файла HVV. Имя и расширение файла не может быть изменено.

 

Описание алгоритма:

Скорости определяются как ординаты максимумов энергетических спектров.

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

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

По первому горизонту суммирование значений отсчетов сейсмотрасс осуществляется по набору гипербол Vогт.

где x - пикет ОГТ, l - удалении ПП от ПВ, V - значение скорости из набора (Start Velocity, Stop Velocity, Velocity Step).

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

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

- координаты точки прихода центрового луча на последнюю из ранее определенных границ.

- угол наклона последней границы.

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

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

 - угол падения центрового луча на последнюю границу.

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

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

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

 - пластовая скорость над последней границей.

 - пластовая скорость в текущем слое, определяемая из набора скоростей (Start Velocity, Stop Velocity, Velocity Step).

  ;   ;      ;

Длина луча в текущем слое для заданного V определяется как

Координаты точки прихода центрового луча на текущую границу определяются как.

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

(j2 - j)

не принимает разные знаки на концах отрезка. Тут  j2 – угол выхода с вертикалью в новом слое с Vconst, j - наклон новой границы под этим слоем. Для центрового луча эти углы должны быть равными. Затем ищется ноль методом Ньютона. С погрешностью 0.00001 радиан  метод сходится за 5-6 итераций. Наклон новой границы определяется как

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

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

Границы могут быть сглажены в скользящем окне.

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

 

Внимание:

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

2. При успешном завершении скоростного анализа программа автоматически вызывает пункт меню Processing->Velocity.

 

К началу раздела.

 

К началу инструкции