Расчет синтетических сейсмограмм 2-Д
Решается прямая задача сейсморазведки в неоднородной среде.
Программа рассчитывает волновое поле в среде на основе решения уравнения эйконала (метод полей времен).
Должна быть описана геометрия профиля, схема наблюдений и создан паспорт.
Внимание:
Границы планшета модели среды определяются пикетами профилей взрыва, приема и ОГТ, описанными в процедурах Design и Binning.
Левая граница планшета определяется минимальным пикетом взрыва, приема и ОГТ.
Правая граница планшета определяется максимальным пикетом взрыва, приема и ОГТ.
Поэтому, невозможно нарастить (удлинить) модель, не удлинив предварительно описание профилей.
Границы планшета должны быть кратными шагу сетки по X.
В основу алгоритма положено
Cohen, J. K. and Stockwell, Jr. J. W., (1998), CWP/SU: Seismic Unix Release 31:
A free package for seismic
research and processing,
Center for Wave Phenomena,
Colorado School of Mines.
Алгоритм основан на формуле (58) в статье N. Bleistein.
Geo. Pros. 34, 686-703.
Для расчета сейсмограмм необходимо создать модель среды, нажав кнопку Edit Model.
Обязательно проверяйте при создании и редактировании модели скорости.
Внимание: процедура игнорирует (пока) делители времен, координат и глубин.
Поэтому при моделировании устанавливайте одинаковыми делители времен, координат и глубин.
Параметры:
Save Model in File - файл, содержащий описание модели среды.
Output Seismic Data Set - файл для записи результата (сейсмограмм).
Start Shot Point - начальный номер магнитограммы для расчета.
End Shot Point - конечный номер магнитограммы.
Можно рассчитывать сейсмотрассы порциями.
Width of Riccer Impulse (ms) - Видимый период волнового поля в мсек
(Длительность импульса Риккера).
Number of Diffractors per Fresnel Zone - Значимое количество колец в зоне Френеля.
Horizontal Velocity range - Расстояние по горизонтали в метрах,
на котором, может происходить значимое изменение скоростей.
Программа будет сглаживать скоростную модель (не границы) на этой базе.
Vertical Velocity range - Расстояние по глубине в метрах,
на котором, может происходить значимое изменение скоростей.
Программа будет сглаживать скоростную модель (не границы) на этой базе.
Параметры Horizontal Velocity range и Vertical Velocity range следует увеличивать, если в информационном окне прохождения процедуры моделирования будет много сообщений об ошибках решения уравнения эйконала. Особенно критичны ошибки для источников, поскольку при этом исключается из расчета целая магнитограмма.
Настоятельно рекомендуется испытать достаточную гладкость скоростей при редактировании модели в режиме Draw->Eikonal Times.
Apply statics - если включить,
то в сейсмотрассы будут введены статические поправки из 1-го раздела паспорта с обратным знаком, (положительная статика увеличит времена). Иначе статика не вводится.
Ввод статики полезен для дальнейшего тестирования программ коррекции статики. Напоминаем, что разделы статики SSP1 и SOP1 заносятся в паспорт профиля программой Geom2D при создании паспорта, если в данном проекте:
- сформирован рельеф и скорости в ВЧР;
- заполнены поля глубин заряда и вертикальных времен в описании магнитограмм.
Smooth Reflectors - если включить,
то отражающие границы будут интерполироваться и сглаживаться кубическим сплайном.
При выключенном флажке границы воспринимаются как прямолинейные отрезки между узловыми точками.
Line Source - если включить,
то предполагается линейный поверхностный источник.
Иначе - точечный источник упругих колебаний.
Для линейного источника потери энергии за геометрическое расхождение будут меньшими.
Output Trace Format - один из форматов записи: СЦС-3 или SEG-Y.
Выходной формат сейсмотрасс (I2, R2, I4, R4) указывается ранее при создании паспорта.
Для расчетных магнитограмм потери за геометрическое расхождение могут быть очень большими (особенно для точечных источников). Поэтому попытку записи в формате не R4 опробуйте предварительно на нескольких сейсмограммах.
Шаг дискретизации трассы и Длина выходной трассы также определяется значениями в паспорте профиля.
Для формата SEG-Y формат R2 тождественен формату R4.
Нажатие кнопки Save JOB позволяет сохранить задание, не выходя из диалогового окна.
Нажатие кнопки Edit Model позволяет создавать и редактировать модель среды.
Нажатие кнопки Export Velocities позволяет экспортировать сетку мгновенных скоростей в формате GS Binary. Сетка вычисляется по данным в описании горизонтов.
Нажатие кнопки Import Velocities позволяет импортировать сетку мгновенных скоростей в формате GS Binary. При этом скорости, указанные в описании горизонтов используются только для расчета коэффициентов отражения.
При включенном флажке From Elevalion источники и приемники будут располагаться на рельефе дневной поверхности. При этом линия приведения должна быть выше самой высокой точки рельефа, а скоростная модель должна быть определена от линии приведения. Лучше всего сделать рельеф отрицательным, а линию приведения – нулевой.
Нажатие кнопки OK запускает выполнение процедуры расчета синтетических сейсмограмм.
Если модель среды не была создана, то расчет синтетических сейсмограмм не выполняется.
Обязательно проверяете при создании и редактировании модели скорости.
Примечания:
1. Допускается исключительно продольный линейный профиль (не Slalom Line).
2. Если в информационном окне прохождения процедуры моделирования будет слишком много сообщений об ошибках решения уравнения эйконала, следует сгладить скорости увеличением параметров Horizontal Velocity range и Vertical Velocity range. В большинстве случаев даже небольшое сглаживание (3-5 шагов) сетки устраняет проблему.
Если увеличение Horizontal Velocity range и Vertical Velocity range не спасает или нецелесообразно, то следует упростить скоростную модель среды.
3.Для расчетных магнитограмм потери за геометрическое расхождение могут быть очень большими (особенно для точечных источников). Поэтому попытку записи в формате не R4 опробуйте предварительно на нескольких сейсмограммах.
4. Границы планшета модели среды определяются пикетами профилей взрыва, приема и ОГТ, описанными в процедурах Design и Binning.
Левая граница планшета определяется минимальным пикетом взрыва, приема и ОГТ.
Правая граница планшета определяется максимальным пикетом взрыва, приема и ОГТ.
Поэтому, невозможно нарастить (удлинить) модель, не удлинив предварительно описание профилей.