Процедура SPMSTAT
Автоматическая коррекция остаточных статических поправок путем максимизации энергии суммирования.
Алгоритм процедуры описан в статье
Ronen, J. and Claerbout, J., 1985,
Surface-consistent residual statics estimation by stack-power maximization:
Geophysics, vol. 50,No. 12 pp 2759-2767.
1. Все трассы в заданном временном окне начитываются в оперативную память в массив data(I,j,t). Поэтому окно корреляции не должно быть очень большим.
2. Начинается итерационный процесс.
2.1. Строим суммотрассы ОГТ. (Суммируем все трассы в пределах каждой сейсмограммы ОГТ)
2.2. Цикл по всем ПВ
2.1.1. Для каждого ПВ получаем т.н. «супертрассу ПВ» - последовательный массив временных интервалов исходных данных всех трасс данного ПВ.
Для этого же ПВ получаем «супертрассу ОГТ» - последовательный массив временных интервалов суммотрасс ОГТ, к которым принадлежат исходные данные в «супертрассе ПВ».
Длина супертрассы определяется как nt_super=(nt+msxh*2)*NTR, где nt – длина временного интервала msxh – максимально допустимый сдвиг, NTR количество трасс в данной сейсмограмме ОПВ.
2.1.2. Вычисляем Функцию взаимной корреляции между супертрассой ПВ и супертрассой ОГТ.
2.1.3. По ординате максимума ФВК определяем задержку для данного ПВ, которую прибавляем к статической поправке за данный ПВ.
2.3. Вводим массив задержек ПВ во все исходные данные.
2.4. Строим новые суммотрассы ОГТ.
2.5. Повторяем цикл 2.2. для всех ПП.
2.6. Вводим массив задержек ПП во все исходные данные.
3. Сохраняем полученные статические поправки в указанных разделах паспорта профиля.
Очевидно, что здесь используется двухфакторная модель статических поправок. Погрешности кинематики и времен t0 процедура пытается загнать в статические поправки за ПВ и ПП.
Параметры процедуры:
Maximum shift (ms) - максимально допустимый сдвиг в статических поправках (мсек).
С увеличением этого параметра существенно увеличивается время работы процедуры.
Number of iterations – Число итераций
С увеличением этого параметра существенно увеличивается время работы процедуры.
Sub - Если 1 то супертрасса ПВ (ПП) вычитается из супертрассы ОГТ перед корреляцией.
Mode – Если =0, то ищется глобальный максимум ФВК, если =1, ближайший к нулю локальный максимум, превышающий заданный процент от значения глобального максимума.
Percent - Процент от значения глобального максимума (Используется только при Mode=1)
Start Time – Начало временного окна в мсек.
Stop Time – Конец временного окна в мсек.
NSP – Номер раздела в паспорте для статических поправок за ПВ.
NOP – Номер раздела в паспорте для статических поправок за ПП.
Horizont Number – Если этот параметр больше нуля, то он воспринимается как номер горизонта в файле горизонтов, вдоль которого будут выбираться участки трасс для обработки.
Horizon Ext - Требуется только если Horizon Number не равен 0. Содержит прослеженные в программе FINES горизонты. Обработка с указанием горизонта возможна только при обработке профилей 2-Д. В 3-Д 6-й параметр воспринимается как ошибка.
Внимание: Горизонт может вообще не придерживаться фазы. Он просто задает переменное по профилю окно выборки фрагментов трасс. При этом параметры Start Time и Stop Time воспринимаются как часть окна над горизонтов (Time Above) и под горизонтом (Time Below). Величина окна при отсутствии горизонта равна разности, а при наличии –сумме
6-го и 7-го параметров.
Процедура еще имеет незавершенный характер и работает долго.
При обработке данных 3-Д создаются громадные временные файлы (1-2 гБайт на 10 кв.км площади) в папке TEMP текущего полтьзователя (переменные среды TEMP и TMP для текущего пользователя). Следите, чтобы в этой папке было достаточно свободного пространства. При нормальном завершении эти файлы автоматически удаляются. При аварийном завершении программы CDPS эти временные файлы сохраняются и должны быть удалены вручную.
При обработке данных 3-Д пока не допускается задание горизонта. Параметр Horizont Number должен быть нулевым.