Фильтр скорости - Velocity filter

А фильтр скорости удаляет мешающий сигналов, используя разницу между скоростями перемещения желаемых сейсмический сигнал и нежелательные мешающие сигналы.

Вступление

В геофизический Датчики приложений используются для измерения и регистрации сейсмических сигналов.[1][2] Доступны многие методы фильтрации, в которых одна форма выходного сигнала создается с более высоким соотношение сигнал шум чем отдельные записи датчика. Фильтры скорости предназначены для удаления мешающих сигналов за счет использования разницы между скоростями распространения желаемой сейсмической волны и нежелательных мешающих сигналов.[3] В отличие от одномерного вывода, создаваемого многоканальной фильтрацией, фильтры скорости производят двумерный вывод.

Рассмотрим массив N датчики, которые получают желаемый и M нежелательные широкополосные помехи. Пусть измерение пth датчик можно смоделировать выражением:

Иксп(t) = амлнsм(т-тмлн) + ŋп(t) ___ (1)

где n = 1,2, ..., N; т = 0,1, ..., М; sм(t) - сигналы, проходящие через решетку, а ŋп(t) представляет белый случайный шум с нулевым средним на пth датчик, несвязанный от датчика к датчику. Параметры aмлн и тмлн - усиление по амплитуде и временные задержки сигнала sм(t) при получении в пth датчик.

Без ограничения общности будем считать, что s0(t) - полезный сигнал, а s1(t), с2(t), ..., сM(t) - нежелательные помехи. Дополнительно будем предполагать, что T0n= 0 и a0n = 1. По сути, это означает, что данные были сдвинуты по времени для выравнивания полезного сейсмического сигнала, чтобы он появлялся на всех датчиках одновременно и был сбалансирован, так что полезный сигнал появлялся с равными амплитудами. Мы предполагаем, что сигналы оцифровываются перед записью и что длина K временных последовательностей записанных данных достаточно велика для того, чтобы в записанные данные были включены полные задержанные мешающие сигналы. В дискретной частотной области пth след можно выразить как:

Иксп(k) = S0(k) + амлнSм(k) e−jwkТмлн + Nп(л) ___ (2)

где k = 0,1, ..., K-1; шk= (2π / K) - угловая частота дискретизации.

Используя матричные обозначения, (2) можно выразить в виде:

X (k) = A (k) S (k) + N (k) ___ (3)
Уравнение 4
Уравнение 5

Фильтрация скорости

Многоканальные фильтры в частотной области F1(k), F2(k), ..., FN(k) может применяться к данным для создания одной единственной выходной трассы в форме:[4][5]

Y (k) = FN(k) Xп(л) ___ (6)

В матричной форме приведенное выше выражение можно записать как:

Y (k) = X '(k) F (k) = S' (k) A '(k) F (k) + N' (k) F (k) ___ (7)

где F (k) - вектор размером N x 1, элементами которого являются фильтры отдельных каналов. То есть,

F (k) = [F1(k), F2(k), ..., FN(k)] ___ (8)

Следуя процедуре, описанной в,[6][7] оптимальный вектор фильтра F (k) может быть разработан для ослабления, в смысле наименьших квадратов, нежелательных когерентных помех S1(k), S2(k), ..., SM(k) при сохранении полезного сигнала S0(k) в Y (k). Этот фильтр может быть показан [6], [7] как имеющий вид:

Уравнение 9

куда час - произвольный ненулевой вектор размером N x 1, ты = [1,0,...,0], я - единичная матрица, Bр(k) - подматрица матрицы, полученной отбрасыванием всех линейно зависимых строк, а L (k) - нижнетреугольная матрица, удовлетворяющая:

[L (k) Bр(k)] [L (k) Bр(k)]ЧАС = Я
Уравнение 10

Схема многоканальной обработки, описанная в [6] - [11], дает одномерную выходную трассу. С другой стороны, фильтр скорости - это двумерный фильтр, который производит двумерную выходную запись.

Двумерная запись может быть сгенерирована процедурой, которая включает многократное применение многоканальных оптимальных фильтров к небольшому количеству перекрывающихся подмассивов входных данных.[8][9]

Рис. 1. Скользящий подмассив многоканальных фильтров.

Более конкретно, рассмотрим подмассив из W каналов, где W << N, который скользит по входным данным, как показано на рисунке 1. Для каждой позиции подмассива может быть спроектирован оптимальный многоканальный фильтр на основе (9) так, чтобы нежелательные помехи подавляются из соответствующей выходной трассы. При разработке этого фильтра мы используем W вместо N в выражении (9). таким образом, трассы 1,2, ..., W входной записи создают первую трассу выходной записи, трассы K, K + 1, ... K + W-1 входной записи создают Kth трассировка выходной записи, а трассы N-W + 1, N-W + 2, ..., N входной записи создают (N-W + 1)ул трассировка, которая является последней трассой выходной записи. Для больших N и малых W, как это обычно бывает в геофизических данных, выходная запись может рассматриваться как сопоставимая по размерам с входной записью. Очевидно, чтобы такая схема работала эффективно, W должно быть как можно меньше; в то же время он должен быть достаточно большим, чтобы обеспечить необходимое ослабление нежелательных сигналов. Обратите внимание, что с помощью такой схемы можно полностью подавить максимум нежелательных помех W-1.[10][11]

Рекомендации

  1. ^ 1. Дж. Х. Джастис, "Обработка массивов в разведочной сейсмологии", в Array Signal Processing, С.Хайкин. Эд. Энглвуд Клиффс, Нью-Джерси: Прентис-Холл, 1985, гл. 2. С. 6-114.
  2. ^ Робинсон Э.А., Робинсон Т. Дуррани, "Обработка геофизических сигналов". Энглвуд Клиффс, Нью-Джерси: Прентис-Холл, 1986.
  3. ^ Р.Л. Сенгбуш, М.Р. Фостер, "Оптимальные многоканальные скоростные фильтры", Геофизика, т. 33, стр. 11-35, февраль 1968 г.
  4. ^ M.T.Hanna и M. Simaan, "Абсолютно оптимальные фильтры для матриц датчиков", IEEE Trans. Акуст., Обработка речевых сигналов, т. ASSP-33, стр. 1380-1386, декабрь 1985 г.
  5. ^ M.T. Ханна и М. Симан. «Матричные фильтры для устранения боковых лепестков», IEEE J. Oceanic Eng., Vol. OE-10, стр. 248-254, июль 1985 г.
  6. ^ СМ. Чен и М. Симаан, "Фильтры частотной области для подавления множественных помех в массиве данных", в Proc. 1990 IEEE COnf. Acoust., Обработка речевых сигналов (Альбукерке, Нью-Мексико), 3-6 апреля 1990 г., стр. 1937-1940.
  7. ^ СМ. Чен, «Оптимальные многоканальные фильтры для множественных нежелательных помех на массивах датчиков», доктор философии. дисс., Сигнальный процесс. Интерпретировать. Лаб., Унив. Питтсбург. Респ. СПИЛ № 91-01, 1991.
  8. ^ М. Симаан, П.Л. Лав, «Оптимальное подавление когерентных сигналов с линейным сдвигом в сейсмических данных», Геофизика, т. 49, стр. 215-226, март 1984 г.
  9. ^ M.T. Ханна и М. Симаан, «Разработка и реализация скоростных фильтров с использованием методов обработки многоканальных массивов», IEEE Trans. Акуст., Обработка речи и сигналов, т. ASSP-35, стр. 864-877, июнь 1987 г.
  10. ^ Чи-Мин Чен и Марван Симаан, «Скоростные фильтры для множественных интерфейсов в двумерных геофизических данных», IEEE Trans. по геонаукам и дистанционному зондированию, Том 29, №4, с. 563-570, июль 1991 г.
  11. ^ Магди Т., Ханна, "Скоростные фильтры для подавления множественных помех в данных геофизических массивов", IEEE Trans. по геонаукам и дистанционному зондированию, Том 26, вып. 6, с. 741-748, ноябрь 1998 г.