Поиск

Полнотекстовый поиск:
Где искать:
везде
только в названии
только в тексте
Выводить:
описание
слова в тексте
только заголовок

Рекомендуем ознакомиться

'Лекция'
Проводя определенную политику в области ценообразования, предприятие активно воздействует как на объем продаж на рынке, так и на величину получаемой ...полностью>>
'Автореферат'
Работа выполнена в Государственном образовательном учреждении высше­го профессионального образования «Волгоградский государственный медицинский уни­в...полностью>>
'Документ'
08:30 со Львова, а в 14:00 с Чопа – переезд на север Венгрии – в Эгерский винодельческий регион – родину огненных красных и нежных белых вин. Самое и...полностью>>
'Программа дисциплины'
Дисциплина «Анализ облигаций» изучается студентами факультета «Экономика», обучающимися по магистерским программам «Финансовые рынки», "Страхово...полностью>>

Метод акустического темного поля

Главная > Статья
Сохрани ссылку в одной из сетей:

Метод акустического темного поля.

В. А. Зверев., П. И. Коротин

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

Один из героев Мольера очень удивился, когда узнал, что всю жизнь говорил прозой. Операцию пространственно-временной фильтрации поля подсветки можно выполнять не задумываясь о том, что переходишь в акустическое темное поле, которое подчиняется своим законам, что многие авторы и делают. Но прозаиком можно стать (то есть, по аналогии, пользоваться операцией перехода в темное поле с наибольшим эффектом) только после понимания особенностей жанра. В случае акустического темного поля это связано и с особенностями схем наблюдения, и с нюансами регистрации сигналов, и с требованиями к измерительному оборудованию, выбором излучателя и типа сигнала подсветки и самой процедурой обработки.

Один из самых распространенных примеров пользователя метода темного поля – лягушка. Берем на себя смелость, не вдаваясь глубоко в бионику и, тем более, зоологию, описать принцип действия ее зрения, которым она в состоянии видеть только движущиеся предметы. Сигнал от света, попавшего на ее глазные рецепторы, не передается в мозг полностью. На ее ступени развития или мозг не в состоянии переработать весь поток информации от глаз, или глазной нерв не в состоянии его передать. Передается в мозг только изменение сигнала на рецепторах. Лягушка, чтобы видеть, не должна шевелиться, так как изображение потонет в шумах. Не беремся утверждать за какое время происходит покадровое усреднение, но, по-видимому, оно определяется характерной скоростью полета насекомых на удалении одного её прыжка. Тогда не мешают помехи в виде качающейся травы и прочего. Этот сигнал достаточен для формирования образа цели, ее идентификации и принятия решения на поражение. В радиолокации этот метод обработки сигналов назвали бы зачатком череспериодного вычитания.

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

Переход от оптики к акустике принципиально не меняет сложности проблем приема, передачи и обработки сигналов. В биоакустике более распространен пассивный прием сигналов, то есть принимается и обрабатывается прямой сигнал от источников звука. Эхолокация, то есть подсветка окружающего пространства и прием отражений, используется живыми организмами только импульсная. Обрабатываются ли ими отраженные сигналы когерентно, то есть с учетом фазы несущей частоты лоцирующего сигнала, сказать однозначно трудно. В принципе по результатам деятельности живых эхолокаторов можно предположить что да. Так дельфины используют для целей на больших дистанциях пачки частых коротких импульсов, летучие мыши определяют цели с точностью намного превышающей длины лоцирующих импульсов [1], и то и другое означает, что сверхразрешение в биоакустике достигается. Элементарный механизм реализации такой обработки это использование органа эквивалентного линии задержки, чтобы сравнивать между собой посланный и принятый сигналы или последовательно отраженные импульсы. Как мы покажем далее такое сравнение близко к переходу в темное поле в варианте импульсной подсветки.

В гидроакустике появляется возможность осуществить на низких частотах тональную подсветку акватории с достаточно большими размерами. Это является принципиальным моментом н е только для технической реализации метода в гидролокации. С использованием тональных сигналов подсветки возникает ситуация, когда пассивный режим может отойти на второй план даже при наблюдении за шумящими объектами. Предположим, что уровни шума (первичного поля) и рассеяния (вторичного поля) объекта в какой-то полосе частот одинаковы, что вполне возможно при соответствующей мощности сигнала подсветки. При общем средстве наблюдения – приемной антенне, казалось бы, пороги обнаружения будут одинаковы. Однако это не так. Поле подсветки может быть задано с определенными свойствами, которые при обработке сигнала понизят порог и повысят вероятность обнаружения вторичного поля. Простейший пример – использование для подсветки тонального высокостабильного сигнала с когерентными свойствами намного выше чем у шумовых или дискретных составляющих первичного поля. При прочих равных условиях обнаружить и выделить такой сигнал проще, проще и наблюдать за движущейся целью в темном поле. Особенно упрощается ситуация в условиях сильной реверберации мелкого моря, когда основной помехой является переотраженный сигнал подсветки. Обычные импульсные гидролокаторы, в независимости от их мощности, незащищены от реверберационной помехи. Для тональной подсветки нет необходимости в высоких уровнях потребляемой мощности. Облегчаются требования к измерительным средствам по производительности обработки и скорости передачи данных. Современный уровень развития техники тоже пока (как и у лягушки) не в состоянии обеспечить передачу и обработку данных с многоканальных антенных приемных систем во широкой полосе частот в реальном масштабе времени. Оценить данные преимущества также является целью работы.

Назовём тёмным полем пространственно-временное распределение принятого акустического поля в локационном эксперименте из которого сигнал подсветки удален пространственно-временной фильтрацией без искажений вторичного поля флуктуирующих рассеивателей. Методом темного поля назовем совокупность операций получения темного поля с целью наблюдения за нестационарной составляющей рассеяния. Последующий анализ темного поля может выполняться произвольными алгоритмами обработки сигналов, в частности, различными видами томографических на анализе которых останавливаться здесь мы не будем.

Существует несколько способов обработки сигналов, которые часто используются в гидролокации, эквивалентные в той или иной степени переходу к темному акустическому полю. Основная цель перехода – максимально удалить или ослабить источник прямого сигнала подсветки как помеху с минимальными искажениями рассеяного поля от движущихся или флуктуирующих объектов.

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

Широко в гидролокации используется доплеровская локация, которая за счет режекции несущей частоты лоцирующего тонально-импульсного сигнала или анализа спектра модуляции выделяет движущиеся рассеиватели и поверхностное волнение, которые в обычном режиме выделить не удается. Однако и здесь возможности гидроакустики существенно уступают радиолокации: отношение несущей частоты к полосе сигнала в гидроакустике намного меньше и в импульсном режиме вступают в противоречие требования к параметрам частотного разрешения, длительности импульса и скважности. Возможности когерентного накопления импульсов чаще всего ограничиваются нестационарностью приемо-передатчиков, устанавливаемых на движущихся или качающихся объектах. Колебания приемо-передатчиков приводят к формированию “пьедестала” в спектре принятого сигнала, что, в первую очередь, исключает возможность наблюдения за малыми рассеивателями.

И всё же в гидроакустике можно добиться существенного подавления прямого сигнала, когда и имеет смысл вводить специальное понятие темного акустического поля. В первую очередь метод требует обеспечения высокостабильной подсветки акватории. Это может быть и тональный гармонический сигнал и импульсный сигнал с синхронизацией запуска и тонально-импульсный с привязкой длительности импульса к частоте несущей. Возможности наблюдения за средой от этих подсветок будут различны, что мы покажем далее, но во всех случаях источник и приемник должны быть жестко связаны по опорным частотам и синхронизации, а также установлены стационарно. Данные условия полностью используют возможности временной режекции прямого поля и позволяют выделять с помощью пространственно развитых приемных систем медленные детерминированные процессы, которые без них полностью тонут во флуктуациях.

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

Из тривиального утверждения следуют принципиальные рекомендации для постановки наблюдений и принципиально новые результаты. Проиллюстрируем данное утверждение численным примером в некоторой усредненной ситуации. В мелком море мы хотим создать комплекс для наблюдения за подводной обстановкой. Для подсветки акватории используем тональный источник, к примеру, на частоте 1кГц и линейную кабельную приемную антенну донной постановки. Частота 1кГц выбрана не слишком высокой из соображений малого затухания и освещения достаточно большой области в мелком море и не слишком низкой, чтобы не пропустить мелкие неоднородности. Приемная антенна пусть состоит из 64 эквидистантных элементов через половину длины волны и, соответственно, имеет апертуру 50 метров и угловое разрешение порядка 20.

Если мы используем яркую всенаправленную “лампу” подсветки, то находясь около неё мы можем ослепнуть (не хватит динамического диапазона приемника) или прищурившись (понизив чуствительность приемника), видеть только то, что находится в непосредственной близости от источника. Этот пример поясняет преимущество бистатической схемы наблюдения перед моностатической и, также естественно, приводит к рекомендации о полезности направленного излучателя подсветки в виде “фонарика” с лучом.

Для наблюдения за рассеянным полем нужно выбрать частоту кадров просмотра. Если в кино принято 24 кадра в секунду, то в гидроакустике все изменений происходят медленнее и вполне достаточно 10 отсчетов в секунду, так как процессы с частотой выше 5Гц в море достаточно редки.

Также как в черно-белом кино нам не важно каков цвет у софитов, так и при обработке темного поля значение несущей частоты уже не имеет значения. Отличие лишь в том, что в в акустике мы можем получить квадратуры сигнала полосе +/-5Гц и использовать информацию о фазе рассеянного поля. Таким образом, для нас в полосе от 995Гц до 1005Гц (или квадратурах сигнала с полосой от 0 до 5Гц) находится вся информация о рассеивателях и процессах в среде, которые рассеивают на частоте 1000Гц.

Пусть мы имеем стабильность излучаемого сигнала порядка 10-6. На частоте 1кГц это соответствует ширине линии 10-3Гц. Необходимо иметь реализацию в 1000 секунд, чтобы определить ширину самой узкой спектральной линии в принятом сигнале. Более длинная реализация для когерентного анализа уже не нужна, так как она испорчена источником подсветки.

Отметим, что оптимально использование частоты излучателя как опорной частоты для гетеродинного приема, тогда длина когерентного приема не огранивается стабильностью генератора излучателя, а уровень помех и их характер будет определяться стационаростью его установки. Гетеродинная регистрация квадратур сигнала – оптимальная процедура как по объему данных, так и по затратам вычислительных ресурсов.

Вычислив одно мгновенное распределение квадратур акустического поля по апертуре антенны (всего 128 чисел) мы получаем полностью текущую информацию, которую может дать данная система наблюдения в конкретный момент. Для получения мгновенного распределения квадратур в полосе 5Гц в нашем случае необходимо обработать 64 реализации сигналов длительностью 0,2сек с частотой дискретизации минимум 2,56кГц (это минимум 33280 чисел).

Фильтрация (в широком смысле) подсветки если она остается аддитивной помехой после дискретизации сигналов с антенной системы - операция линейная и не влияет на реализацию последующих операций обработки. Но при обработке дискретных по пространству и времени и ограниченных массивов данных возможны два случая. Первый, когда значение частоты подсветки совпадает с одной из точек отсчета спектра и второй, когда это значение приходится на интервал, находящийся между точками отсчета спектра. Спектр можно понимать и в пространственной [2, 3] области по апертуре антенны и во временной [4] реализации по отдельному каналу. В первом случае стабильная помеха не маскирует собой сигнал при любом уровне помехи и с ней можно не бороться. Во втором случае, в независимости от стабильности помехи, её вклад присутствует во всех точках отсчета спектра. В этом случае помеха при достаточном ее уровне может маскировать любой сигнал.

Второй случай наиболее вероятен. Если помеха действительно занимает всего одну точку отсчета, то второй случай можно свести к первому путем сдвига всего спектра на такую величину, при которой спектр помехи приходился бы целиком на одну из точек отсчета спектра. Однако спектр помехи имеет пусть весьма малую, но конечную ширину. В этом случае в спектре помехи обязательно (при любом сдвиге спектра) появится компонента, не совпадающая ни с одной из точек отсчета. В этом случае сигнал подсветки на 40…60дБ превышающий фон в состоянии маскировать сигнал по всему его спектру.

Как эффективно отфильтровать такой сигнал? С помощью мультипликативного фильтра (например прямоугольного окна в частотной области) в частотной области можно вырезать узкий участок спектра, содержащий помеху. Однако мультипликативный фильтр неэффективен для узкополосной помехи, частота которой не совпадает в точности с точкой отсчета спектра. На этом рисунке видно, что фильтр с прямоугольной частотной характеристикой (рис.1а) не влияет на помеху между отсчетными точками.

Уровень в дБ

Частота в единицах квантования

Рис.1. Частотная характеристика мультипликационного фильтра, полученная только по частотам, взятым в точках отсчета (а) и по частотам, находящимся в промежутках между точками отсчета (b).

Для помехи занимающей одну точку спектра будет хорошо работать адаптивный фильтр, учитывающий амплитуду и фазу помехи и вычитающий её. Этот способ эффективен для узкополосной помехи, занимающей одну точку отсчета спектра уже вне зависимости от того, совпадает эта точка с какой-либо точкой отсчета спектра или находится в промежутке между ними.

В случае, когда спектр помехи шире одной отсчетной точки, предложен метод, названный впервые в [4] методом акустического темного поля. Этот метод использует фильтр, который сочетает достоинства мультипликативного и адаптивного фильтров, устраняя их недостатки. Частотная характеристика фильтра имеет вид синусоиды с нулем, приходящимся на точку максимума помехи. Эта точка определяется адаптивным методом. Такой фильтр может быть сформирован как в частотной области [2, 3], так и в области сигналов [4, 5, 6]. Для формирования такого фильтра в области сигналов достаточно произвести вычитание сигналов, содержащих помеху, сдвинутых по времени на величину , что является неким аналогом череспериодного вычитания в радиолокации. Величина этого временного сдвига определяет период образующейся в спектре синусоиды. Нуль на максимум помехи настраивается путем ввода дополнительного фазового множителя (при вычитании), определяемого с помощью адаптивной процедуры. Результаты показаны на рисунке 2.

Уровень в дБ

Уровень в дБ

Рис.2. Спектр сигнала с помехой (а вверху). Результат применения адаптивного фильтра для подавления помехи (b вверху). Подавление помехи методом темного поля (а и b внизу).

Здесь показана сравнительная эффективность двух методов фильтрации: адаптивного и метода темного поля. Помехой служит узкополосный сигнал, обладающий конечной шириной спектра. Такой сигнал был приготовлен численно так, как это описано в [7]. На рис.2а вверху показан спектр сигнала с помехой. Сигналом служил тоже узкополосный сигнал, спектр которого лежит вблизи спектра помехи. К сигналу добавлен широкополосный шум в качестве дополнительной помехи. На спектре сигнала с помехой сигнал еле-еле виден, так как его закрывает спектр узкополосной помехи, распространяющийся на все точки отсчета спектра. На рис.2б вверху показан результат фильтрации помехи с применением адаптивного метода. Помеха в ее максимуме полностью убирается, однако, спектр вне помехи остался практически прежним и все еще маскирует сигнал.

На нижних кривых рис.2а,б показан результат применения метода темного поля к тому же сигналу. В результате помеха убрана не только в ее максимуме, а и в ближайших окрестностях от него. В результате сигнал отчетливо виден. Стал заметен также и шум, добавленный к сигналу. Этот шум на верхних графиках также маскируется узкополосной помехой.

Амплитуда

Частота в единицах квантования

Рис.3. Частотные характеристики фильтров темного поля, полученные путем одно и двукратного вычитания задержанных сигналов (а) и вычитанию текущего среднего значения (по [6]) (b). Характеристики определены по частотам, находящимся как в точках отсчета, так и в промежутках между ними.

На рис.3 показаны варианты частотных характеристик фильтров, осуществляющих фильтрацию по методу темного поля. Фильтры, сформированные вычитанием смещенных реализаций, особенно удобны, когда роль узкополосной помехи выполняет сигнал, имеющий спектр, близкий к нулю частот. Эти случаи весьма распространены, например, когда выполнено гетеродинирование сигнала, но опорная частота гетеродина несвязана с несущей частотой генератора излучателя подсветки.

На рис.3а показана частотная характеристика фильтра, полученная вычитанием сигналов через некоторый временной интервал. Этот временной интервал является параметром фильтра, который может быть оптимизирован с целью выделения нужного сигнала или группы сигналов. На рис.3б показана частотная характеристика фильтра, полученного путем вычитания из сигнала его скользящего среднего. Время интегрирования при получении скользящего среднего тоже является параметром фильтра, который может быть оптимизирован. Такой фильтр предложен в [5,6].

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

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

Наиболее легко интерпретируемую информацию дает горизонтальная линейная антенна. Следить за подводной обстановкой в мелком море для тональной подсветки с помощью горизонтальной антенны можно либо в плоскости “время - пеленг”, либо “частота-пеленг”. Для вертикальной антенны остаются, например, возможности слежения за неоднородностями в просветном режиме [6] и выделения их по вкладу в рассеяние в отдельные моды [8], что однако требует специальных фазированных маломодовых излучающих систем и накладывает жесткие ограничения на особенности распространения. Исключение возникновения зон тени может потребовать комбинации горизонтальной антенны с другими или выбирать место её установки по глубине.

Для нашего примера при гетеродинном приеме поток данных невелик и можно копить реализацию непрерывно до нескольких суток (это десятки ГБт, то есть на уровне возможностей современных носителей) и выделять тем самым частоты флуктуаций связанные даже с приливно-отливными течениями и внутренними волнами. Процессы с флуктуациями 10-3…10-4Гц при наблюдениях и обработке одно-двух часовых реализаций выглядят в спектре как эффективное уширение дискреты подсветки. Принципиальных ограничений для их изучения нет, но практический интерес представляют более быстрые явления. Например движение сосредоточенной неоднородности в диапазоне скоростей от 1м/с до 10м/с. Это соответствует доплеровским частотам порядка герца и на размерах акватории порядка нескольких десятков километров, в которой удается создать достаточный уровень подсветки, изменениям пеленга порядка градуса в минуту, то есть эти процессы не пересекаются по частотному диапазону и не мешают наблюдению друг за другом.

Приведенный ниже результат эксперимента показывает проявление в темном поле небольшого движущегося рассеивателя - трубы (5 метров длиной и 40 см в диаметре) с балластом.

а)

б)

в)

г)

д )

Рис.4. Результаты наблюдения за движущейся неоднородностью в темном поле: а - амплитуда сигналов с антенны в полосе 5Гц, б – пространственно-временной спектр, в – маска для пространственно-временного спектра, г – траектория движения рассеивателя в плоскости «время-пеленг» после действия маски, д – восстановленное из «4г» распределение вторичного поля рассеивателя на апертуре приемной антенны.

Длина рассеивателя в двести раз и диаметр в полторы тысячи раз меньше базы разнесения (около 500 метров) приемника и передатчика в бистатической схеме наблюдения. Эксперимент соответствует рассматриваемому выше примеру (с масштабом 1/10) и подробно описан в [9]. Рассеиватель пересекал трассу распространения звука примерно перпендикулярно со скоростью около 0,5м/сек. Пеленг горизонтальной приемной антенны на излучатель - 00. Приемная 64-х элементная антенна установлена в полводы на дне, излучатель подсветки установлен на дне (глубина ~15м). Излучатель ненаправленный с частотами в несколько килогерц. В полосе 1Гц он обеспечивает на приеме отношение сигнал/шум порядка 60дБ. Запись сигналов - гетеродинная со сведением частот излучателя и приемника, т.е выполнены условия метода темного поля.

На рис.4а временная реализация длиной ~850 секунд, в ней 2048 точек, то есть реализация прорежена до 3-х отсчетов в секунду для того чтобы можно было в одну запись обработать один галс движения рассеивателя. Маска (рис.4в) выбрана так, чтобы в первую очередь подавить прямой сигнал и сохранить в пространственно-временном спектре вторичное поле движущегося рассеивателя. По рис.4б видно, что форма маски не является критичной для выделения отметки рассеивателя. Маска в пространственно-временном спектре необходима в случае сильной поверхностной реверберации. Она позволяет выделить движущуюся дискретную неоднородность даже если ее сила цели меньше силы цели поверхностной реверберации. На рис.5 представлен случай сильной поверхностной реверберации.

Рис.5 Вверху пространственно-временной спектр поля на апертуре горизонтальной антенны. Внизу временной спектр для сигнала с пеленга -60. На графиках вырезана полоса частот +/-0,0125Гц.

График рис.5 получен в той же измерительной схеме, что и на рис.4б, но получен при сильном поверхностном волнении, которое вызвало появление модуляции с частотой в районе 1 – 1,5Гц. Отметка цели в виде наклонной линии имеет несколько меньшую силу (на 3-6дБ) чем поверхостная реверберация. Реверберация весьма неравномерно распределена по углам, что говорит о анизотропии поверхностного волнения.

Выбирая соответствующую маску в рамках данного эксперимента можно выделить и наблюдать отдельно движущийся рассеиватель, либо исследовать свойства реверберации и поверхностного волнения. Пространственно-временная структура вторичного поля практически не искажена колебаниями приемо-передающих элементов и мы в состоянии измерять параметры полей на 70-80дБ меньших поля подсветки. К примеру в [10] доплеровскую томографию взволнованной поверхности пытаются проводить с движущегося судна при максимум 40дБ диапазоне и ширине линии несущей частоты около 0,05Гц.

Увидеть в суммарном (рис.4а) поле вторичное поле дискретной неоднородности (рис.4д) практически невозможно. Без операции перехода в темное поле расеиватель не виден и в зависимости пространственного спектра поля от времени. На рис.6 дана зависимость пространственного спектра поля на апертуре горизонтальной антенны (по вертикальной оси синус пеленга) от времени (единица по горизонтальной оси соответствует примерно минуте).

а)

б)

Рис.6.

Рис.6 получен в том же эксперименте и при слабом волнении, но для четырех галсов движения рассеивателя поперек трассы распространения звука. Для получения графика рис.6б использован фильтр с двойным вычитанием, характеристика которого показана на рис.3а, без маски в области пространственных частот.

Отметим еще раз, что отношение амплитуд прямого и вторичного поля на рис.4а и 4д составляет от 60 до 70дБ. Наблюдения за столь малыми рассеивателями возможны только при высоких требованиях к шумам измерительной аппаратуры, к разрядности АЦП ( в данном случае 16) и шумам цифровых гетеродинных фильтров (64-х разрядные операции).

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

Частотная характеристика фильтра в методе темного поля имеет нули, определяемые временной задержкой вычитаемой реализации. При переходе к вычитанию через период следования задержка должна равняться периоду следования импульсов подсветки. Для дискретных сигналов в этом случае становится принципиальным требование кратности периода следования импульсов частоте дискретизации, что обеспечит возможность когерентного накопления сигнала с данной дистанции.

В гидроакустике имеют практическое значение только тонально-импульсные сигналы и применение гетеродинного ввода для уменьшения потока данных с антенны на регистрацию и обработку остается возможным, но заметно снизить поток данных он уже не в состоянии вследствии малого отношения ширины полосы сигнала и несущей частоты. В нашем примере для несущей частоты 1кГц и длительности импульса 2-3 мсек минимальный объем данных с которым можно работать с когерентным накоплением по 64 импульсам при наблюдении за пространством размерами порядка 1км (период следования 2сек) составит около 2,5МБт.

Основное преимущество которое дает импульсная подсветка заключается в возможности представления темного поля в координатах “пеленг-дальность” и “время-дальность” для локализации и выделения рассеивателей. Однако возможности когерентного накопления для импульсного режима уже ограничиваются техническими проблемами. Напомним, что при обработке темного поля на рис.4а мы работали с массивами квадратур сигнала размерами 2048 на 64, что невызывает сложностей на стандартном персональном компьютере. В импульсном режиме с реализациями такой длины работать уже затруднительно.

Сформулируем алгоритм обработки импульсных сигналов эквивалентный переходу в темное поле. Переходим в принятом сигнале к квадратурным составляющим с опорной частотой излучателя для каждого приемника антенны. Для одиночного приемника формируем массив в координатах «номер импульса (время) – дальность (задержка от начала посылки)». При выполнении условия кратности периода следования импульсов частоте дискретизации получаем возможность работать с реализацией сигнала рассеянного на определенной дальности. К этой реализации применимы операции фильтрации рассмотренные ранее по устранению первичного поля и для выделения определенной области частот. После выделения нужного сигнала по всем элементам антенны в координатах «время-дальность» можно выполнить ее фазирование и перейти в координаты «пеленг-дальность» для последовательных моментов времени.

Другой вариант обработки заключается в первоначальном фазировании антенны на выбранный пеленг и формировании массива данных в координатах «номер импульса (время) – дальность (задержка)» для одного пеленга и получению реализации от локализованного объема как по дальности, так и по пеленгу. С ней также выполняются операции фильтрации прямого поля. Выполнив последовательно операции для всех приемных пеленгов можно перейти к слежению за подводной обстановкой в координатах «пеленг- дальность» в виде просмотра отдельных кадров по времени.

Представим результат наблюдения за движением рассеивателя в гидролокационном эксперименте в соответствии с последними алгоритмами.

Рис.7 Результаты наблюдения за движущейся неоднородностью в темном поле в импульсном режиме: амплитуда отраженного сигнала в градациях серого на плоскости «дальность (задержка) – время (номер импульса)» для одного пеленга антенны (–500). Ширина полосы сигнала 1кГц, длительность импульса – 2мс, несущая частота – 2кГц, период следования 200мс. Апертура антенны – 6 метров. Подсветка ненаправленная.

На рис.7 показана совокупность рассеянных сигналов в координатах «время - дальность» для одного приемного пеленга антенны за часть галса движения рассеивателя. Вверху сигнал без обработки, внизу сигнал после фильтрации низких частот в реализациях равной дальности (горизонтальных секущих графика). Градации серого на обоих графиках отнормированы к максимуму. Эксперимент описан подробно в [11].

Масштабы рассеивателя, антенны и акватории примерно теже, что в эксперименте работы [9] результаты которого на рис.4 и 5. По горизонтали рис.6 отложено время (примерно 200 секунд или 1024 посылки). По вертикали отстробированное время задержки (в отсчетах через 0,8мсек от 126 до 242) от прихода прямого импульса посылки на приемную антенну до задержки соответствующей дальности примерно 150 метров. Таким образом каждая вертикальная секущая рисунка формируется эхом одной посылки. Две темные параллельные полосы вверху рис.6а соответствуют прямому сигналу и донной реверберации, серые горизонтальные полосы на больших задержках – реверберации и неподвижным отражателям. На рис.7б тот же сигнал после фильтрации низких частот в сигналах равной дальности (по горизонтальным секущим рис.7а) для устранения вклада неподвижных рассеивателей, донной и части поверхностной реверберации. Если на рис.7а рассеиватель проявляется очень слабо, то на рис.7б полезный сигнал уверенно выделяется в виде наклонной линии - смещения рассеивателя вдоль одного пеленга на фоне поверхностной реверберации.

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

В эксперименте использовался электродинамический излучатель генерирующий импульс длительностью 1 мсек с частотой заполнения 5,3 кГц. Приемный микрофон находился на удалении около пяти метров от излучателя. На рис.8а – исходный сигнал, 8б – фон в темном поле, 8в – темное поле при одновременном колебании двух рассеивателей. По вертикальной оси отложен номер импульса (от 1 до 64; всего примерно 3,5 секунды) по горизонтальной – время запаздывания эха (16 мсек или 1024 отсчета; отрезки до момента прихода прямого сигнала на микрофон и в конце - отброшены).

Два рассеивателя (жесткие сферы радиусом около 5см) раскачивались между излучателем и приемником примерно по середине, но на разной высоте и разными периодами колебаний.

Период следования импульсов локации - 64 мсек (4096 отсчетов на каждый период следования, 64 отсчета на каждый импульс посылки).

а)

б)

в)

Рис.8 Результаты наблюдения за двумя движущимися неоднородностями в импульсном режиме с помощью одиночного приемника. Показана амплитуда отраженного сигнала в градациях серого на плоскости « время (номер импульса) - дальность (задержка)».

Условия метода темного поля выполняются в эксперименте засчет общего опорного генератора для ЦАПа излучателя и АЦП приемника. На рис.8а нет никаких признаков движущихся рассеивателей (аналогично рис.4а), рис.8б приведен как фоновый уровень помех, когда рассеиватели неподвижны, и на рис.8в легко видеть сигнал от движущихся рассеивателей в координатах «дальность-время».

Стационарная реверберация, как видим из рис.7 и 8, устраняется легко. Для флуктуирующей помехи частота следования импульсов для когерентной обработки определяется частотой флуктуаций fф, то есть, если надо отфильтровать поверхностную реверберацию с частотой 0,1Гц, надо за 10 секунд послать минимум три импульса. Это требование ограничивает размер зоны наблюдения D (так как она тоже определяется периодом следования импульсов) - D<c/fф. В мелком море максимум флуктуаций около 0,2Гц и область наблюдения D, соответственно, порядка 7км. При несоблюдении этого условия реверберационная помеха в накопленном с определенного удаления сигнале засчет редкого стробирования становится широкополосной помехой. Это таже ситуация, когда в оптике при стробоскопическом освещении рассеиватель колеблется быстрее, чем частота стробоскопа. В такой ситуации выделить медленые процессы, например движущийся дискретный рассеиватель, можно, но поверхностная реверберация будет уже искажена.

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

Дадим грубую оценку мощности источника подсветки. Вернемся к численному примеру рассмотреному выше. Пусть источник подсветки развивает на расстоянии 1 метр акустическое давление 1000Па (180дБ отн.мкПа). Это совсем немного – около 5Вт акустической мощности. Возьмем за грубую оценку потерь на распространение на 10км величину 90дБ (для глубины акватории 100 метров взят сферический закон спадания для 100м – 40дБ, затем цилиндрическое спадание от 100 метров до 10км ещё 40дБ и потери 1дБ/км). На приемной антенне уровень поля составит 90дБ. Пусть уровень шума моря на частоте 1кГц составляет 70дБ в полосе 1Гц, что является заведомо высоким уровнем помех. Считаем шум изотропным и широкополосным, тогда в полосе 10-3Гц, соответствующей стабильности подсветки, с выхода одного приемного пеленга антенны уровень шума составит 22дБ. При неизотропных шумах выигрыш антенны для пеленгов без сосредоточенных источников будет только больше (чем корень из 64 или 18дБ), а при наличии дискретной помехи около несущей частоты всегда можно перестроить частоту подсветки. То есть 22дБ это оценка для помехи сверху и мы получим на приеме в узкой полосе сигнал/шум порядка 60дБ, что примерно соответствует условиям экспериментов [9, 11] и [16, 17].

Можно считать экспериментальным фактом, что для данного отношения сигнал/шум в полосе до +/-2Гц от несущей частоты фон в большинстве ситуаций будет определяться модуляцией сигнала поверхностным волнением. Глубина модуляции тонального сигнала поверхностным волнениемна стацтрассе является функцией большого числа переменных: гидрологии, скорости ветра, частоты сигнала, длины трассы, профиля дна, угла между трассой и направлением ветра и т.д. Теоретические оценки следует делать в привязке к конкретному району размещения системы наблюдения, так как шум моря формируют и факторы не связанные с поверхностным волнением: судоходство, прибой, биологические шумы. Однако укажем, что у поверхностной реверберации и у предельно низкого шума моря – причина общая – ветер, т.е. при ослаблении ветра падает и уровень шума и поверхностная реверберация. Выделив на стацтрассе при каком-то ветре модуляцию тональной подсветки поверхностным волнением мы с большой долей вероятности можем утверждать, что будем наблюдать ее при любой скорости ветра при неизменных остальных условиях.

По большому количеству экспериментальных работ в море [13, 14, 15, 16] с различными измерительными средствами, условиями распространения, дистанциями и схемами измерений глубина модуляции тонального сигнала изменяется в пределах от –15дБ до –50дБ. Например из [14], на стацтрассе “Сахалин – Итуруп” длиной 300км, получено –15дБ зимой и -40дБ там же, но летом и со сформировавшимся ПЗК. Примеры более коротких трасс в мелких морях: -30дБ на 10км в Балтийском море [16], -45дБ на 13км в Баренцевом море [17]. Это подтверждает то, что для нашего численного примера, когда сигнал превышает шум на приеме не менее чем на 60дБ, фон темного поля определяется рассеянием на поверхностном волнении. Данная ситуация будет сохраняться до дистанций заведомо не менее 30-ти километрового разнесения излучателя и приемной антенны. Таким образом до данной дистанции помехой выделению движущейся сосредоточенной неоднородности в темном поле будет являться только поверхностная реверберация. При этом уровень помехи (соответственно и достижимая дальность наблюдения) может меняться в зависимости от погодных условий на величину до 20…25дБ, настолько, насколько меняется сила цели поверхностной реверберации [9] в совокупности с изменением потерь на распространение [14].

Оценку предела обнаружения неоднородности можно получить из теоретической оценки силы цели поверхностной реверберации. Всевозможные операции обработки (типа применения маски - рис.4в) или привлечения статистических характеристик эхо-сигналов [11], естественно, понижают данный порог.

Теоретические оценки силы цели поверхностного волнения дают достаточно хорошие результаты, например [14], но, однако, очень ценны прямые измерения, дающие непосредственно количественный результат. Так в эксперименте, описанном в [9], специально измерялась сила цели тестового рассеивателя, по которой можно сопоставить с помехой (рис.5) и оценить уровень реверберационной помехи количественно. Для нашего численного примера, который примерно соответствует данному эксперименту в масштабе 1/10, на расстоянии 10км порядок величины силы цели поверхностного волнения составит по данным [9] при полном штиле в полосе 0,1Гц с выхода антенны по одному пеленгу не более 100м2, при скорости ветра 5м/с не более 500м2. Эту оценку можно использовать для анализа результатов [16, 17].

Максимумы поверхностного волнения для мелкого моря имеют частоты от 0,15 до 0,4Гц и волнение ближе к изотропному (особенно в закрытых акваториях), в глубоком море это частоты от 0,1Гц до 0,2Гц и волнение чаще имеет выделенное направление, что делает помеху неизотропной. Это создает диапазоны параметров по скорости и направлению движения сосредоточенных неоднородностей неудобные для наблюдения за ними. Для разработки схем наблюдения для конкретных районов и выбора частот и интенсивности подсветки существует необходимость предварительных измерений. Под схемой понимается также место установки приемной антенны и излучателей подсветки и их тип (апертура приемника, направленность подсветки и т.д.). Мультистатическая схема наблюдения позволит устранить зоны тени, а для упрощения технологии удобно размещать несколько подсвечивающих излучателей и одну приемную антенну (или две горизонтальные в одном месте, повернутые на 900).

Оценка значения порога обнаружения по теоретическому значению силы цели рассеивателя дает завышенный результат. Дело в том, что реверберационное поле само подсвечивает рассеиватель и способно увеличивать его эффективную силу цели. Экспериментально обнаружено увеличение отношения сигнал/шум при увеличении скорости ветра от 1-2м/с до 6-8м/с (что соответствует росту боковиков, вызванных поверхностной реверберацией, на 12-15дБ) на величину до 6дБ [18]. В бистатической схеме сила цели максимальна в просветной ситуации, когда рассеиватель пересекает трассу распространения от излучателя к приемнику. На удалении от просветной точки сила цели значительно падает (на 20-30дБ для волновых размеров рассеивателя больше единицы), однако рассеиватель всегда находится в реверберационном поле, пусть и меньшем прямого, но приходящем на него со всех сторон, в частности и в просветном для приемника ракурсе. Это приводит к ситуации когда в бистатической схеме наблюдений при увеличении волнения рассеиватель на больших углах от направления на излучатель начинает лучше выделяться в пространственно-временном спектре темного поля (сравните рис.4а и рис.5), несмотря на рост помех. В оптике такая ситуация естественна – темный объект хорошо виден на светлом фоне, а при отсутствии фона только при пересечении им прямого луча подсветки.

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

В заключении выскажем благодарность принимающим участие в обсуждении проблемы, экспериментах и разработке измерительного оборудования для них Матвееву А.Л., Салину Б.М., Турчину В.И., Цибереву А.В, Чащину А.С.

    Литература

1.Морозов В.П.. Занимательная биоакустика. М.: Знание, 1987, 206 с.

    2.Зверев В.А. Антенна темного поля // Акуст.ж.1994 т.40, №3. С. 401-404.

    3.Зверев В.А., МатвеевА.Л., Славинский М.М., Стромков А.А. Фокусируемая антенна темного поля // Акуст.ж. 1997 т.43,№4

    4.Зверев В. А. Акустическое темное поле // Акуст.ж. 2000 т.46, №.1. С.103–111.

5.Горский С.М., Зверев В.А, Матвеев А.Л.,.Митюгов В.В. Некогерентное накопление сигналов акустической дифракции. Акустический журнал, т.41, N 2, 1995, с.223-231.

6.Зверев В.А., Матвеев А.Л., Митюгов В.В. Согласованная фильтрация откликов акустической дифракции при некогерентом накоплении на вертикальной антенне. Акустический журнал, т.41, N 4, 1995 г.,c. 591-595

7.Зверев В.А., Стромков А.А. Выделение сигналов из помех численными методами. – Нижний Новгород: ИПФ РАН, 2001, 188с.

8.Нечаев А.Г., Хилько А.И. Дифференциальная акустическая диагностика случайных неоднородностей океана // Акуст.ж., 1988, т.34, №2, с.285-289.

9.Зверев В.А., Коротин П.И., Матвеев А.Л.,Митюгов В.В.,Орлов Д.А., Салин Б.М., Турчин В.И. Экспериментальные исследования дифракции звука на движущихся неоднородностях в мелководных условиях // Акуст.ж, 2001, т.47, №2, с.227-237.

10. Бурлакова И.Б., Дубовой Ю.А., Зейгман А.Л. и др. О возможности акустической томографии взволнованной поверхности океана //Акуст.ж., 1988, т.34, №3, с.423-430.

11. Зверев В.А., Коротин П.И. Качественное отличие полезного сигнала от поверхностной реверберации при когерентном накоплении импульсов активной локации //Акуст.ж., 2002, т.48, №2, с. …..

12. Зверев В.А., Коротин П.И. Интерференционная локация высокого разрешения // Акуст.ж., 2001, т.47, №4, с.489-493.

13. Бреховских Л.М., Лысанов Ю.П. Теоретические основы акустики океана. Л.: Гидрометеоиздат, 1982, 264с.

14. Артельный В.В., Вировлянский А.Л. и др. Статистическая модель низкочастотных акустических полей в приповерхностных волноводах в глубоком море// Н.Новгород, в Сборнике трудов ИПФ РАН «Формирование акустических полей в океанических волноводах», 1991, с.149-173

15. Гончаров В.В, Зайцев В.Ю., Куртепов В.М. и др. Акустическая томография океана. Н.Новгород: ИПФ РАН, 1997, 254с.

16. Gershman A.V., Turchin V.I., Zverev V.I. Experimental results of localization of moving underwater signal by adaptive beamforming //IEEE Transaction on Signal Processing, 1995, v.43(10), p.2249-2257

17. Матвеев А.Л., Митюгов В.В. Отчет по 7-му рейсу НИС “Ак. Сергей Вавилов”, Москва, ИОФ АН СССР, 1990, гл.7

18. Zverev V.A., Korotin P.I. Effect of Surface Reverberation on the Selection of a Diffracted Signal // Acoustical Physics, 2002, v.48, №3, p.340-345

19. Зверев В.А., Коротин П.И., Матвеев А.Л., Салин Б.М., Турчин В.И. Обращенный апертурный синтез в акустическом темном поле // Акуст.ж., 2000, т.46, №5, с.650-657.



Скачать документ

Похожие документы:

  1. Определение пола птиц по ДНК

    Документ
      Многие виды птиц не имеют ярко выраженного полового диморфизма или же приобретают его в только в половозрелом возрасте. Морфологические промеры, определение по поведению не всегда надежны.
  2. «Методы микроскопии»

    Учебная программа курса
    2. Разработка сетевых программ обучения и подготовки слушателей и студентов для сети Образовательных учреждений на базе межрегионального отраслевого ресурсного центра
  3. І. П. Основи дефектоскопії-К.: «Азимут-Україна», 2004. 496 с. Ермолов И. Н., Останин Ю. А. Методы и средства неразрушающего контроля качества. М.: Высшая школа, 1988. Неразрушающий контроль металлов и изделий. Справочник

    Справочник
    Основные принципы и положения диагностики, дефектоскопии и неразрушающего контроля качества металлических изделий. Приобрести практические навыки относительно их выбора и применения.
  4. Возможность классификации отражателей акустических волн по частотным характеристикам эхо-импульсов

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

    Документ
    Психологическая наука исключительна не только в силу уникальности своего предмета – самого тонкого, сложного, возвышенного и противоречивого, чем является душа, мысль, сознание человека, но и в силу высочайшей степени свободы, которую

Другие похожие документы..