Доклады БГУИР
2011 №4 (58)
УДК 621.391
ОСОБЕННОСТИ КОРРЕЛЯЦИОННОЙ ОБРАБОТКИ СЕЙСМИЧЕСКИХ
СИГНАЛОВ
Н.Г. КУХАЛЬСКИЙ, ПИ. КОЗЕКА*
Институт пограничной службы Республики Беларусь Калиновского, 4, Минск, 220103, Беларусь
*Войсковая часть 01549 Бобруйская, 3, Минск, 220600, Беларусь
Поступила в редакцию 31 января 2011
Рассмотрены варианты корреляционной обработки принятого сейсмического сигнала движущихся наземных объектов и их влияние на величину отношения сигнал/шум на выходе обнаружителей, применяемых в схемах технических средств охраны различного назначения.
Ключевые слова: сейсмология, корреляция, спектр.
Введение
Физический смысл корреляционной обработки заключается в том, что принятые колебания перемножаются с опорным сигналом, соответствующим по форме и параметрам точно ожидаемому сигналу, а результаты произведения интегрируются [1]. Корреляционная обработка выявляет сходство (корреляцию) принимаемых колебаний с ожидаемыми.
Сейсмический сигнал движущегося наземного объекта, как правило, имеет неизвестные параметры по времени появления, длительности и форме огибающей, его преобладающей частоте. Такая неопределенность значительно усложняет задачу обнаружения по сравнению с вариантами выделения радиолокационного сигнала с известными параметрами, когда сигнал сначала формируется, излучается, а потом принимается и проходит корреляционную обработку относительно известной модели сформированного сигнала. Среда распространения сейсмических волн является частотнозависимой и значение преобладающей частоты сейсмического полезного сигнала может изменяться в зависимости от среды и дальности до цели от единиц до десятков Герц [2]. Поэтому примеры корреляционных обнаружителей, которые описаны во многих работах [3-5], посвященных методам обработки принятых радиолокационных сигналов, могут быть использованы для проектирования устройств обнаружения сейсмических сигналов после учета их спектрально-временных особенностей.
Особенности корреляционной обработки сейсмических сигналов
Рассмотрим пример корреляционной обработки сейсмического сигнала движущегося человека, когда в качестве опорного сигнала коррелятора используется сигнал, по положению и форме совпадающий с огибающей ожидаемого сигнала шага человека с внутренним заполнением колебаниями на преобладающей частоте принимаемого сигнала (рис. 1,а) и без (рис. 1,б).
а б
Рис. 1. Опорный сигнал: а - с заполнением; б - без заполнения
На рис. 2,а показан сейсмический сигнал шага человека, принятый с дальности 15 м, а на рис. 2, б его энергетический спектр.
а б
Рис. 2. Принятый сейсмический сигнал одного шага человека (а - сейсмограмма полезного сигнала с шумом; б - энергетический спектр полезного сигнала и шума)
Преобладающая частота полезного сигнала, изображенного на рис. 2,а составляет около
60 Гц.
Рассмотрим обработку принятого сейсмического сигнала корреляционным обнаружителем с одним квадратурным каналом (рис. 3).
/.»г.+'(>,) X Накопитель ДПФ 1 |2 *
т
Рис. 3. Коррелятор с одним квадратурным каналом В этом случае при совпадении частоты fop и фазы фор опорного сигнала и преобладающей частоты fv и фазы фу принятого (рис. 4,а), после перемножения происходит перенос сигнала по частоте (рис. 4,6) на величину разности и суммы частот перемножаемых сигналов согласно выражению:
Uf 8 g sin (2nfvnT )• UoPn 8 g„ sin (2fpnT )= = X UfUopn 8 L (cos2™T f + fop )-С082тагд f - f0p )),
где Uf - принятый сейсмический сигнал; 8g - огибающая сигнала на выходе сейсмоприем-
ника; n - номер дискретного отсчета; Тд - частота дискретизации; Uop - опорный сигнал.
Отношение сигнал/фон в данном случае увеличивается на 6-7 дБ, в том числе за счет совпадения огибающих принятого и опорного сигналов по положению.
На рис. 8 (график 1) показана зависимость изменения отношения сигнал/фон рsf после перемножения от величины расстройки частоты опорного сигнала. Анализ графика показывает,
что при отклонении частоты опорного сигнала от частоты принятых колебаний на Д[в = ±3 Гц, отношение сигнал/фон увеличивается на 5-6 дБ; при отклонении на А/в = ±5 Гц отношение сигнал/фон увеличивается на 3-4 дБ, при отклонении частоты на А/в = ±7 Гц - на 0,8-1 дБ, при
отклонении фазы опорного сигнала от фазы принятых колебаний на Дфор = ± ^ радиан отношение сигнал/фон увеличивается на 2,9-3,1 дБ, а при отклонении на Дфор =±^2 - на 0,30,5 дБ соответственно (рис. 9, график 1).
а б
Рис. 4. Сейсмический сигнал после перемножения (а - принятый сигнал и сигнал опорного генератора; б - энергетический спектр сигнала после перемножения)
Для исключения потери полезного сигнала из-за априорной неопределенности его начальной фазы в схеме обнаружителя могут использоваться два квадратурных канала (рис. 5).
Рис. 5. Коррелятор с двумя квадратурными каналами Принятый сигнал перемножается с опорным сигналом с ожидаемой частотой и фазой в первом канале и с фазой, сдвинутой на ^, во втором канале. На вход устройства ДПФ (дискретного преобразования Фурье) поступают две квадратурные составляющие g и g1. Над ними проводится операция ДПФ согласно выражениям:
^ БТ-1
^ = Re G„ = £
и г cos
41
Г 2пп1Л N
\1УБТ У
+ иг sm
41 ±
N
V i убт У у
"БТ-1
^ != 1т Gn = £
и, cos
41 ±
N
\ 1 У БТ
- и , sm
41
N
у1 УБТ У У
где Gn - значения спектральных составляющих обрабатываемого сигнала; NFT - размер спектральной выборки.
Совокупность значений двух ортогональные составляющих gn и gn ± определяет значения спектральных составляющих обрабатываемого сигнала Gn.
Как и в предыдущем случае при совпадении частоты и фазы опорного сигнала и преобладающей частоты и фазы принятого, после перемножения происходит перенос сигнала по частоте на величину разности и суммы частот перемножаемых сигналов, отношение сигнал/фон в данном случае увеличивается на 6-7 дБ (рис. 6).
На рис. 8 (график 2) показана зависимость изменения отношения сигнал/фон после перемножения от величины расстройки частоты опорного сигнала для схемы обнаружителя с двумя квадратурными каналами. Анализ графика показывает, что при отклонении частоты опорного сигнала от частоты принятых колебаний на Л/в = ±3 Гц, отношение сигнал/фон увеличивается на 5,8-6 дБ, при отклонении на Л/в = ±5 Гц отношение сигнал/фон увеличивается на 5,2-5,5 дБ, при отклонении частоты на Л/в = ±7 Гц - на 4,5-5 дБ, при отклонении фазы
опорного сигнала от фазы принятых колебаний на Лфор = ± радиан отношение сигнал/фон
увеличивается на 5,7-6 дБ, при отклонении на Лфор = ±/2 - на 5,4-5,6 дБ соответственно (рис. 9, график 2).
а б
Рис. 6. Сейсмический сигнал после перемножения (а - принятый сигнал и опорные сигналы; б -энергетический спектр сигнала после перемножения)
Для устранения априорной неопределенности относительно значений преобладающей частоты и фазы принимаемого сейсмического сигнала после перемножения в схеме устройства корреляционной обработки в первом и во втором случае необходимо предусмотреть многоканальную обработку.
Для снижения аппаратных затрат, позволяющих устранить априорную неопределенность относительно преобладающей частоты принимаемого сейсмического сигнала, а так-же его начальной фазы, предлагается использовать в качестве опорного сигнала коррелятора сигнал, по форме и положению совпадающий с огибающей ожидаемого сигнала шага человека без внутреннего заполнения (рис. 7,а).
а б
Рис. 7. Сейсмический сигнал после перемножения (а - принятый сигнал и сигнал опорного генератора без заполнения; б - энергетический спектр сигнала после перемножения)
В этом случае значение отношения сигнал/фон после перемножения инвариантно преобладающей частоте (рис. 8, график 3) и фазе (рис. 9, график 3)принимаемого сейсмического сигнала. Отношение сигнал/фон в данном случае увеличивается на 3-4 дБ за счет совпадения огибающих принятого и опорного сигналов по положению (рис. 7, б).
4 2
°-7 "3.5 0 3.5 7 Щ, ' Ц
Рис. 8. Зависимость отношения сигнал/фон после перемножения от расстройки по частоте опорного сигнала: 1 - график, соответствующий опорному сигналу с заполнением в схеме с одним квадратурным каналом; 2 - график, соответствующий опорному сигналу с заполнением в схеме с двумя квадратурными каналами; 3 - график, соответствующий опорному сигналу без заполнения
Анализ графика показывает, что при отклонении частоты опорного сигнала от частоты принятых колебаний на А/в = ±3 Гц отношение сигнал/фон увеличивается на 2,8-3 дБ, при отклонении на А/в = ±5 Гц отношение сигнал/фон увеличивается на 2,2-2,5 дБ.
Р„,ДБ 8
б
4
2
-1 57 Афор, радиан
Рис. 9. Зависимость отношения сигнал/фон после перемножения от расстройки по фазе опорного сигнала: 1 - график, соответствующий опорному сигналу с заполнением в схеме с одним квадратурным каналом; 2 - график, соответствующий опорному сигналу с заполнением в схеме с двумя квадратурными каналами; 3 - график, соответствующий опорному сигналу без заполнения
Заключение
Проведенные исследования показали, что максимально возможные потери в отношении сигнал/фон при использовании в качестве опорного сигнала без внутреннего заполнения, позволяющего снизить аппаратные затраты при построении коррелятора, по сравнению с опорным сигналом, совпадающим по форме и положению и имеющим внутреннее заполнение, составляют около 5 дБ, а в среднем, из практики, 1-3 дБ и менее.
Для осуществления обработки сигнала с неизвестным временем запаздывания волнового пакета необходимо использовать многоканальные схемы, устраняющие априорную неопределенность относительно этого параметра.
FEATURES OF SEISMIC SIGNALS CERELATION PROCESSING
N.G. KUKHALSKI, P.I. KOZEKA
Abstract
The variants of the seismic signal correlation processing of moving ground targets and their impact on the ratio of signal to noise at the output of detectors used in the schemes of technical means of protection for different purposes, are studied.
Литература
1. Ширман Я.Д., Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех. М., 1981.
2. Кухальский Н.Г. // Докл. БГУИР. 2008ю № 3. С. 116-118.
3. Вайнштейн Л.А., Зубаков В.Д. Выделение сигналов на фоне случайных помех. М., 1960.
4. Охрименко А.Е. Основы радиолокации. М., 1983.
5. Ширман Я.Д., Лосев Ю.И., Минервин Н.Н.и др. Радиоэлектронные системы: основы построения и теория. М., 1998.