Процедура FLTRGM
Автоматическая рекурсивная остронаправленная режекторная фильтрация.
Очень часто промышленные электрические наводки не равны строго 50 или 60 Гц и их частота может меняться от сейсмограммы к сейсмограмме. Особенно это актуально, если наводки осуществляются находящейся поблизости дизельной электростанцией. Кроме того, помеха может присутствовать не на всех трассах.
Для этого случая и разработана процедура FILTRGM.
Параметры:
FREQ - частота режекции в герцах в первом приближении (целое или действительное число).
Частота задается именно в герцах (независимо от значения делителя времен).
DF - частотная полоса, в пределах которой программа будет уточнять значение частоты режекции.
Значение также задается именно в герцах (независимо от значения делителя времен).
Программа уточняет значение частоты режекции в диапазоне частот (FREQ-DF - FREQ+DF)
Coef - Пороговый коэффициент.
Задается как целое или вещественное положительное число.
Решение на применение режекторной фильтрации принимается только в том случае, если найденное максимальное значение амплитудного спектра в диапазоне частот (FREQ-DF - FREQ+DF) превышает средний по модулю уровень амплитудного спектра в этом же диапазоне частот в Coef раз. В противном случае сейсмотрасса остается неизменной.
Start Time - Начало окна для определения амплитудного спектра.
Stop Time - Конец окна для определения амплитудного спектра.
Окно определения амплитудного спектра рекомендуется задавать в конце записи, где уровень наводок превышает уровень отраженных волн. В этом случае достаточно значения Coef=3. Если окно определения амплитудного спектра захватывает интервал полезных волн, то Параметр Coef следует снижать до 1 или меньше.
Значения Start Time и Stop Time должны быть умножены на делитель времен, указанный в Geom2D.
Nmax - Число отсчетов,
до которого будет удлинено окно анализа нулями с целью увеличения разрешенности спектра. Чем больше это число, тем точнее будет определена частота режекции.
Собственно режекторная фильтрация осуществляется по алгоритму FILRG.
Программу следует применять по исходным сейсмограммам до каких-либо регулировок амплитуд.
Можно предварительно применить процедуру Gain для приведения общего уровня амплитуд к уровню 1000 (весьма приближенно).