Статья: Разрешение и идентификация нескольких шумовых источников

Разрешение и идентификация нескольких шумовых источников

В.М. Кузькин¹, С.А. Пересёлков^1,2, Г.Н. Кузнецов¹, И. В. Казначеев²

¹ Институт общей физики им. А.М. Прохорова РАН (ИОФ РАН), Москва, Россия.

² Воронежский государственный университет (ВГУ), Воронеж, Россия.

Аннотация

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

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

RESOLUTION AND IDENTIFICATION OF FEW NOISE SOURCES

V.M. Kuz'kin¹, S.A. Pereselkov^1,2, G.N. Kuznetsov¹, I.V. Kaznacheev²

¹ General Physics Institute of RAS (GPI RAS), Moscow, Russia

² Voronezh State University (VSU), Voronezh, Russia

Abstract. The results of numerical experiment of resolution and identification of three noise sources with different intensity at presence of isotropic noise in the ocean waveguide are presented. The numerical simulation is based on the interferometric method of source localization. The comparative analysis of estimation accuracy of coordinates sources: direction, radial speed, range, depth is carried out. The numerical simulation is implemented by using vector-scalar receivers.

Keywords: hydroacoustics, sound field, moving broadband source, noise source, interference structure, hydrolocation.

Введение

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

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

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

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

Интерферометрический метод локализации источника

шумовой источник изотропный помеха

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

Применительно к скалярной компоненте поля (давление P) спектрограмма имеет вид

Здесь и - циклическая частота и время спектрограммы; радиальная скорость источник; время наблюдения; полоса сигнала, и t циклическая частота и время интерферограммы. По отношению к другим компонентам поля ВСП формирование спектрограмм аналогично. Интерферограммы разных компонент поля ВСП и их комбинаций подобны по форме и различаются только значениями интенсивности.

Спектральная плотность (1) локализована в двух полосах плоскости , зеркально перевернутых относительно начала координат, в виде отдельных фокальных пятен. Их изображения расположены в первом и третьем квадрантах, если радиальная скорость , т.е. источник приближается к приемнику и во втором и четвертом квадрантах(), когда источник удаляется от приемника. Область локализации содержит главных максимумов с координатами расположенными на прямой . Здесь M число мод, формирующих поле, номер фокального пятна. Ближайший к началу координат пик, обусловленный интерференцией соседних мод, расположен в точке . Координаты пика, вызванного интерференцией мод номеров , и т.д. И, наконец, координаты самого удаленного пика, порожденного интерференцией первой и последней моды . Фокальные пятна сосредоточены в полосе, ограниченной прямыми

где полуширина фокальных пятен в направлении оси ; в направлении оси полуширина . Вне этой полосы спектральная плотность практически подавлена.

Положения максимумов фокальных пятен пропорциональны радиальной скорости и удалению источника от приемника

Где

параметры, определяющие характерные пространственные и частотные масштабы изменчивости передаточной функции волновода. Здесь qвходное отношение сигнал/помеха (с/п); , постоянная распространения m-й моды; , групповая скорость m-й моды; черта сверху означает усреднение по номерам мод, точка сверху обозначает оценку координат источника. Соотношение (4) справедливо при выполнении условия .

В качестве критерия обнаружения источника принимается наличие выраженного пика сигнала, преобладающего над помеховыми пиками, функции обнаружения

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

Минимальное (предельное) отношение с/п, при котором реализовывается когерентное накопление спектральных максимумов вдоль интерференционных полос, т.е. источник обнаруживается, а оценки его координат близки истинным значениям, для шумового источника на фоне изотропной помехи оценивается как (по мощности) [2, 5]. Здесь Jчисло временных отсчетов при преобразовании Фурье (1).

Итерационная процедура оценки пеленга осуществляется с использованием интерферограмм x-й составляющей колебательной скорости и комбинационной компоненты.На первом шаге определяется отношение максимумов функций обнаружения

Алгоритму (8) свойственна неоднозначность определения направления на источник: значения пеленгов , равнозначны.

На втором шаге область спектрограмм за пределами зоны фокусировки очищается от помехи и выполняется двукратное обратное преобразование Фурье по восстановлению интерферограмм. Анализируется корреляция между интерферограммами и в фиксированный момент времени t. Если корреляция положительная, то пеленг расположен в первом или третьем квадрантах горизонтальной плоскости ВСП, если корреляция отрицательная, то во втором или четвертом квадрантах. Данное условие можно формализовать в форме соотношения

Условие соответствует положительной корреляции, а условие отрицательной корреляции.

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

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

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

обращается в нуль, . Здесь

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

Разрешение нескольких источников

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

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

Численный эксперимент

Волновод полагается горизонтально-однородным. Распределение скорости звука по глубине приведено на рис. 1.

Рис. 1. Профиль скорости звука.

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

ВСП размещены на глубине м, расстояние между ними м. Положительные направления осей xВСП и системы координат совмещены. Шумовые точечные источники на глубине движутся по прямолинейной траектории с постоянной скоростью , . Источник приближается к ВСП, источники -удаляются. Начальные координаты источника : удаление км, скорость м/с, глубина м, пеленг , угол траектории , радиальная скорость м/с, входное отношение с/п (-22 дБ). Начальные координаты источника:удаление км, скорость м/с, глубина м, пеленг , угол траектории , радиальная скорость м/с, входное отношение с/п (0 дБ). Начальные координаты источника:удаление км, скорость м/с, глубина м, пеленг , угол траектории , радиальная скорость м/с, входное отношение с/п (-10 дБ). Геометрия задачи показана на рис. 2.

Рис. 2. Схема расположения источников по отношению к ВСП.

Сигналы источников и распределенная помеха являются белым шумом. Время наблюдения мин. Ширина полосы Гц. Шумовое поле регистрируется в течение с, временной интервал между регистрациями с. Число временных точек отсчета , так что предельное входное отношением с/п (-23 дБ).

При моделировании использовалось семь мод, их постоянные распространения и групповые скорости приведены в табл. 1.

Таблица 1.Постоянные распространения и групповые скорости мод на частоте Гц.

Для иллюстрации линейной суперпозиции формирования изображений источников на рис. 3 представлены нормированные интерферограмма, спектрограмма и функция обнаружения суммарного поля трех шумовых источников одинаковой интенсивности в отсутствие помехи. Нормированные величины обозначаются сверху знаком «крышка».

(а) (б) (в)

Рис. 3.Нормированнаяинтерферограмма (а), спектрограмма (б) и функция обнаружения (в) трех шумовых источников одинаковой интенсивности в отсутствие помехи.