Алгоритм программы AKSP системы СОС-ПС
Процедура для каждой из исходных трасс определяет временной сдвиг относительно адаптирующейся эталонной модельной трассы и показатель качества, т.е. показатель отношения сигнал/помеха, представляющий собой амплитуду ФВК.
На вход процедуры подаются трассы подобранные по ОГТ, с выравненными амплитудами, введенными кинематическими поправками, априорными СТП, после мьютинга и редакции.
На выходе процедуры получаются:
- квазитрассы, в которых записаны определенный по каждой трассе временной сдвиг и амплитуда ФВК, данной трассы с промежуточной эталонной трассой;
- финальные эталонные трассы (модели) для каждой ОГТ.
Алгоритм работы процедуры состоит из следующих шагов:
1. Амплитуды по каждой трассе сейсмограммы ОГТ в пределах используемых окон приводятся к единому уровню путем умножения последних на величину 16000/D, где D – среднеквадратическое значение амплитуды в пределах задаваемого окна; также приводится уровень амплитуд предварительной модели.
2. Вычисляется ФВК каждой трассы данной сейсмограммы ОГТ Fjkh(T) с предварительной эталонной трассой – MPjh(T), где: h – номер заданного в пределах трассы окна, k – номер трассы в сейсмограмме ОГТ, j – номер ОГТ; каждая амплитуда ФВК делится на величину n*D, т.е. на произведение среднеквадратического значения амплитуды (D) на число выборок в заданном окне, чтобы избежать возможного в этом случае переполнения.
3. По ФВК определяется наибольший максимум, далее по трем точкам с помощью параболической интерполяции, т.е. высокоточно, определяется значение амплитуды ФВК и соответствующий ему временной сдвиг – Qjkh и Tjkh. Выбор максимума ФВК осуществляется с использованием весов, величины которых изменяются от 1 при 0-вом сдвиге до 0.5 при крайних сдвигах равным MTS. Это означает, что для выбора больших сдвигов они должны иметь большую амплитуду.
4. Каждая трасса сдвигается на соответствующую ей величину Tjkh.
5. Вычисляется промежуточная эталонная трасса по формуле:
,
где: Rj – это доля запоминания предварительной модели, которая задается
в операторе процедуры и которая может изменяться для разных
ОГТ и по разным окнам;
MPjh(T)- является предварительной моделью;
Wjkh=Qjkh/(1-Qjkh2) – весовой коэффициент, значение которого,
благодаря использованию формулы, существенно более
дифференцировано по сравнению со значением самой величины
Qjkh;
n – общее число трасс в сейсмограмме ОГТ.
6. Вычисляется ФВК каждой трассы данной сейсмограммы ОГТ с найденной промежуточной моделью – MIjh(T) и определяются новые значения амплитуды – Qjkh и соответствующее ей значение сдвига – Tjkh.
7. По новым значениям Tjkh и Qjkh вычисляется средневзвешенная величина сдвига:
8. Определяются остаточные значения сдвигов:
Исходные трассы сдвигаются на величину остаточных сдвигов (T’jkh) и суммируются. Крайние значения T’jkh для каждой ОГТ печатаются в листинге (Tmin, Tmax).
9. Предварительная
модель – MPjh(T) сдвигается
на величину 
и в результате взвешенного
суммирования с суммой сдвинутых трасс находится финальная модель:
10. Предварительная модель для следующей (j+1) ОГТ получается путем весового суммирования финальной и предварительной модели текущей (j) ОГТ:
где:
- сдвинутая предварительная
модель;
- финальная модель;
Wjh=
/(1-2)
где:
– т.е. среднее из Qjkh;
Перечисленные действия повторяются для каждой ОГТ.
Результаты вычислений: Qjkh – показатель качества (отношение сигнал/помеха) и сумма остаточного временного сдвига с «накопленным» сдвигом:
