Техника и методика геофизического эксперимента : сборник научных трудов / Рос. акад. наук, Кол. науч. центр, Поляр. геофиз. ин-т. – Апатиты : [б. и.], 2003. – 194 с.

Для каждого обнаруженного описанным выше способом вступления сейсмической волны из буфера извлекается 1 с данных и вновь фильтруется тем из полосовых фильтров, для которого отношение STA/LTA оказалось максимальным. Все три канала фильтруются по отдельности. Затем рассчитываются две оценочные функции для определения поляризации, зависящие от предполагаемого направления на событие а . Первая из них R(a) рассчитывается по горизонтальным каналам ( a -угол от 0 до 360°): R(a) изменяется от - 1 до 1 и максимально для Р и Rg-волн в том случае, когда a или а+180° - азимут на событие, для S и Lq-волн - когда а перпендикулярно азимуту на событие. Однако функция R не позволяет различить тип волны и определить знак азимута. Поэтому вводится вторая функция - корреляция между вертикальной компонентой и проекцией горизонтальной на определенное направление: CZ(a) = Corr( E-cosa+N-sina, Z), где E, N, Z - отсчеты восточного, северного и вертикального каналов. Если a - азимут на событие, то для P-волны CZ(a) максимально, для Rg CZ(a)=0, для S - CZ(a+90°)=0. Указанные функции рассчитываются для углов от 0 до 360 с шагом в 3°, значения запоминаются (рис.2). S(a) - S(a+90°) R(a) = S(a) + S(a+90°) где (Е-отсчеты восточного, N-северного каналов) S(a) = £ | EjCosa + Njsina | . I Рис.2. Типичный вид функций R u C Z для P и S - волн 64