Акустический журнал, 2020, T. 66, № 4, стр. 437-445
Интерференционный метод оценки координат движущегося шумового источника в мелком море с использованием высокочастотных сигналов
С. А. Пересёлков a, b, *, В. М. Кузькин a, **, Г. Н. Кузнецов a, Д. Ю. Просовецкий b, С. А. Ткаченко b
a Институт общей физики им. А.М. Прохорова Российской академии наук
119991 Москва, ул. Вавилова 38, Россия
b Воронежский государственный университет
394006 Воронеж, Университетская пл. 1, Россия
* E-mail: pereselkov@yandex.ru
** E-mail: kumiov@yandex.ru
Поступила в редакцию 07.08.2019
После доработки 21.02.2020
Принята к публикации 25.02.2020
Аннотация
Приведены результаты высокочастотного эксперимента по локализации движущегося шумового источника с использованием цилиндрической малогабаритной векторно-скалярной антенны. Использовалась частотно-временная обработка, согласованная с интерференционной картиной, формируемой источником. Восстановлены временные зависимости пеленга, скорости, удаленности и глубины источника. Дано качественное и количественное объяснение экспериментальным данным на основе модельной двухлучевой интерферограммы, образованной прямым лучом и лучом, отраженным от свободной поверхности.
ВВЕДЕНИЕ
Обнаружение и оценка координат источников по их шумовому полю является одной из актуальных проблем пассивной гидроакустики [1]. В настоящее время для ее решения в низкочастотном диапазоне (десятки и сотни герц), когда применимо модовое описание структуры поля, разработан и успешно апробирован интерферометрический метод локализации шумового источника [2–8]. В его основе лежат механизмы межмодовой интерференции и волноводной дисперсии, обусловливающие формирование частотно-временной (пространственной) интерференционной картины (интерферограммы) [9–11]. В отличие от традиционной согласованной со средой обработки [12, 13] интерферометрическая обработка позволяет реализовать одновременно обнаружение источника и определение его параметров (пеленг, расстояние, глубина, скорость). При этом она может быть реализована в отсутствие знаний о передаточной функции волновода [7, 8], что значительно расширяет ее область применения. В настоящее время, на примере интенсивных внутренних волн, экспериментально показана возможность применения интерферометрии для восстановления передаточной функции невозмущенного волновода и диагностики временной изменчивости океанической среды [14, 15].
В области высокочастотного (килогерцового) диапазона интерферометрическая обработка применяется впервые и нуждается в экспериментальной проверке, так как изменяется механизм формирования интерферограммы. Если в области низких частот физико-математические принципы, лежащие в основе локализации малошумных источников, ясны [2–4], то остается открытым вопрос: насколько успешно этот метод может применяться в высокочастотном диапазоне?
В настоящей работе представлены результаты обработки высокочастотного эксперимента по локализации движущегося шумового источника в окрестности траверса, когда одновременно изменяются расстояние и пеленг. Прием осуществлялся на цилиндрическую вертикальную векторно-скалярную антенну (далее антенна). Изложен алгоритм интерферометрической обработки шумового поля и на его основе получены временные оценки параметров источника. Теоретически проанализирована динамика модельной двухлучевой интерферограммы, адекватно передающей основные свойства экспериментальной интерферограммы.
РЕЗУЛЬТАТЫ НАТУРНОГО ЭКСПЕРИМЕНТА
Цель эксперимента заключалась в проверке работоспособности интерферометрического метода локализации шумового источника в высокочастотной области. Эксперимент проводился в мелководной акватории с глубиной волновода $H\sim 87$ м и скоростью звука, практически не зависящей от глубины, $c\sim 1400$ м/с. Обрабатывался сигнал с одиночного приемника антенны, расположенного на глубине ${{z}_{q}} = 30$ м. Интерферограмма источника регистрировалась в полосе $\Delta f = 1{\kern 1pt} --{\kern 1pt} 15$ кГц. Длительность принимаемой шумовой реализации $\delta t = 0.2$ с. Входное отношение сигнал/помеха (по мощности) $q\sim 8$. Априорная информация о гидрофизических характеристиках акватории отсутствовала.
Схема движения источника относительно антенны в окрестности траверса показана на рис. 1. Оси x и y векторно-скалярного приемника (далее приемника) Q совпадают с положительными направлениями системы координат. Из начальной точки A источник приближался к приемнику, проходил точку B траверса и далее от нее удалялся по направлению к точке C. После прохождения точки C источник совершал петлеобразные траектории (на рис. 1 они не отображены). Горизонтальное расстояние от источника до антенны обозначено r, пеленг – θ.
На рис. 2 приведена экспериментальная интерферограмма ${{I}_{e}}\left( {f,~t} \right) = {{\left| {P\left( {f,~t} \right)} \right|}^{2}}$ скалярной компоненты (давление P) поля источника. С целью повышения контрастности и информативности на ней вырезано среднее значение. В момент времени $~t = 10$ с, отвечающий прохождению точки траверса, частотный масштаб изменчивости интерферограммы оценивается как $\Lambda = 404$ Гц. При удалении от точки траверса частотный масштаб возрастает, что, естественно, приводит к увеличению ширины интерференционной полосы. При этом он практически от частоты не зависит. Временной масштаб D изменчивости интерферограммы с возрастанием частоты уменьшается и увеличивается с возрастанием времени наблюдения (расстояния). В точке траверса, например, на частоте $f = 3.5$ кГц он равен $D = 2.9$ c. Крутизна интерференционных полос в области траверса возрастает со временем и частотой, и по мере удаления от нее уменьшается. Данные закономерности характерны до момента времени $t{\kern 1pt} '\sim 70$ c, т.е. по истечении 60 с после удаления от точки траверса. Для значений $t > t{\kern 1pt} '$ интерферограмма носит осциллирующий характер, вызванный многократным изменением направления движения источника (см. рис. 8).
Спектральные плотности, получаемые двукратным преобразованием Фурье интерферограммы $I\left( {\omega ,~t} \right)$ скалярной компоненты поля,
(1)
${{F}_{i}}\left( {\tilde {\nu },~\tau } \right) = \int\limits_{{{t}_{i}}--\frac{{\Delta t}}{2}}^{{{t}_{i}} + \frac{{\Delta t}}{2}} {\int\limits_{{{\omega }_{0}}--\frac{{\Delta \omega }}{2}}^{{{\omega }_{0}} + \frac{{\Delta \omega }}{2}} {I\left( {\omega ,~t} \right)} } \exp \left[ {i\left( {\tilde {\nu }t--\omega \tau } \right)} \right]dtd\omega ,$Для оценки направления на источник применялся интерферометрический метод пеленгования, основанный на отношении очищенных от помех интерферограмм
(2)
$\frac{{G\left( {f,~{{t}_{i}}} \right)}}{{R\left( {f,~{{t}_{i}}} \right)}} = \sin \theta \left( {{{t}_{i}}} \right),\,\,\,\,~\frac{{U\left( {f,~{{t}_{i}}} \right)}}{{R\left( {f,~{{t}_{i}}} \right)}} = \cos \theta \left( {{{t}_{i}}} \right)$Определение параметров источника осуществлялось с использованием согласованной обработки по отношению к интерферограмме, формируемой движущимся источником. Горизонтальное расстояние до приемника r, глубина z и скорость w оценивались как координаты основного максимума трехмерной взаимокорреляционной функции экспериментальной ${{I}_{e}}$ и модельной ${{I}_{m}}$ интерферограмм
(3)
$E\left( {r*,~z*,~w{\text{*}}} \right) = \sum\limits_k {\sum\limits_i {{{I}_{e}}} } \left( {{{f}_{k}},~{{t}_{i}}} \right)I_{m}^{*}\left( {{{f}_{k}},~{{t}_{i}}} \right),$Отметим, что ранее в килогерцовом диапазоне с применением двух- и многолучевой моделей поля, используя многоэлементные приемные вертикальную или горизонтальную антенны, основываясь на алгоритме традиционной согласованной обработки, экспериментально показана возможность оценки глубины и расстояния излучателя от антенны [16–18], однако, источник был неподвижным.
Для отслеживания временного восстановления параметров источника использовалось каскадное вычисление взаимокорреляционной функции. В начальный момент времени ${{t}_{0}} = 0$ осуществлялась обработка (3) во временном интервале $\Delta t = 20$ с. В момент времени ${{t}_{1}} = {{t}_{0}} + \,\,~\delta t{\kern 1pt} '$ проводилась следующая обработка в течение времени $\Delta t$, $\delta t{\kern 1pt} ' = 5$ с, и т.д. В результате за время наблюдения $T = 275$ с получено $N = \left[ {{{\left( {T - \Delta t} \right)} \mathord{\left/ {\vphantom {{\left( {T - \Delta t} \right)} {\delta t{\kern 1pt} '}}} \right. \kern-0em} {\delta t{\kern 1pt} '}}} \right] + 1 = 52$ оценок параметров источника.
Модельная интерферограмма ${{I}_{m}}\left( {f,~t} \right)$ (рис. 5), адекватно передающая основные свойства экспериментальной интерферограммы, строилась на основе интерференции полей двух лучей с использованием лучевой программы. При этом предполагалось, что в зоне траверса источник двигался по прямолинейной траектории.
Для иллюстрации на рис. 6 изображены двумерные разрезы нормированной трехмерной взаимокорреляционной функции (3) для двух моментов времени. Нормированные восстановленные параметры источника обозначены сверху прямой линией. Скорость нормирована на значение 30 м/с; расстояние и глубина нормированы на значение 87 м. По результатам моделирования параметры источника оцениваются как $r = 189$ м, $z = 10.9$ м, $w = 27.6$ м/с для ${{t}_{i}} = 10$ с; $r = 1049$ м, $z = 12.8$ м, $w = 26.9$ м/с для ${{t}_{i}} = 50$ с. За время наблюдения $T = 0{\kern 1pt} - {\kern 1pt} 275$ c были обработаны интерферограммы для всех временных интервалов и с использованием алгоритма (3) восстановлены временные зависимости параметров источника $~z\left( t \right)$, $r\left( t \right)$, $w\left( t \right)$.
На рис. 7 показана нормированная зависимость глубины $z\left( t \right)$. Как видно, во время эксперимента глубина источника изменялась в пределах ~1.5 м (в основном, при поворотах) и в среднем оценивается как $z = 11.4$ м.
С использованием зависимостей $\theta \left( t \right)$ и $r\left( t \right)$ вычислены координаты источника
(4)
$X\left( t \right) = r\left( t \right)\sin \theta \left( t \right),~\,\,\,\,Y\left( t \right) = r\left( t \right)\cos \theta \left( t \right)$Скорость источника w на временном интервале $\Delta t$ определялась как
(5)
$w\left( {\Delta t} \right) = \frac{{\sqrt {{{X}^{2}}\left( {\Delta t} \right) + {{Y}^{2}}\left( {\Delta t} \right)} }}{{\Delta t}}.~$Нормированная временная зависимость $w\left( t \right)$ отображена на рис. 9. Видно, что на поворотах не только изменяется глубина, но и уменьшается скорость. Как следствие, зависимость скорости от времени имеет сложный изрезанный характер. Восстановленные координаты источника соответствуют реальным значениям, погрешность не превышает 10%.
ТЕОРЕТИЧЕСКАЯ МОДЕЛЬ ИНТЕРФЕРОГРАММЫ
Задача теоретического рассмотрения состояла в том, чтобы разобраться, как характеристики интерферограммы (или положения максимумов фокальных пятен голограммы) связаны с параметрами движущегося источника. Обсудим это на простом примере удаления источника от точки траверса B по прямолинейной траектории параллельно оси x с постоянной скоростью w и на фиксированной глубине z (рис. 1, пунктир). Допустим, что коэффициент отражения от свободной поверхности не зависит от частоты, угла падения и равен –1.
Рассмотрим для примера зону, в которой источник начинал движение из точки траверса в сторону удаления от приемника. Тогда двухлучевая интерферограмма принимает вид [19]
(6)
${{I}_{t}}\left( t \right) = \frac{{2{{W}^{2}}}}{{{{{\bar {R}}}^{2}}\left( t \right)}}{\text{cos}}\Phi \left( {{\omega },~t} \right),~$Здесь $\eta = {{2{{z}_{q}}z} \mathord{\left/ {\vphantom {{2{{z}_{q}}z} c}} \right. \kern-0em} c}$, c – скорость звука, ${{z}_{q}}$ – глубина приемника; $~r\left( t \right) = \sqrt {r_{0}^{2} + {{{\left( {wt} \right)}}^{2}}} $ – горизонтальное расстояние от приемника до источника в момент времени t, $~{{r}_{0}}$ – расстояние в начальный момент времени $t = 0$ от точки траверса до приемника; $\bar {R}\left( t \right) = \sqrt {{{r}^{2}}\left( t \right) + z_{q}^{2}} $; ${{W}^{2}}$ – коэффициент, характеризующий мощность излучения. Принимая в среднем глубину источника $z = 11$ м, имеем $\eta = 0.47$ м с.
Кривая интерференционных максимумов определяется условием
Используя (7), из (8) для производной кривой постоянной фазы получаем
Отсюда следует, что крутизна интерференционных полос в области траверса возрастает при увеличении времени, частоты, скорости источника и по мере удаления от точки траверса уменьшается. Подобный результат зарегистрирован и в эксперименте. Решая дифференциальное уравнение (9) при начальном условии $f\left( 0 \right) = {{f}_{0}}$, находим уравнение кривой интерференционных максимумов поля источника
Частотный масштаб изменчивости интерферограммы в момент времени t, как следует из (7), равен
так что он возрастает с возрастанием скорости источника и удалением от приемника, уменьшается с увеличением заглубления и не зависит от частоты. Подобные закономерности наблюдались и в эксперименте. Например, в точке траверса (${{r}_{0}} = 189$ м) получаем $\Lambda \left( 0 \right) = 402$ Гц, что согласуется с экспериментальным значением $\Lambda = 404$ Гц. Если ширину интерференционной полосы $\Delta f\left( t \right)$ определить как ширину, в пределах которой фаза (7) меняется на ${{\pi } \mathord{\left/ {\vphantom {{\pi } 2}} \right. \kern-0em} 2}$, то $\Delta f\left( t \right) = {{\Lambda \left( t \right)} \mathord{\left/ {\vphantom {{\Lambda \left( t \right)} 4}} \right. \kern-0em} 4}$. Величина $\Delta {\tau }$, обратная частотному масштабу, $\Delta {\tau } = {1 \mathord{\left/ {\vphantom {1 \Lambda }} \right. \kern-0em} \Lambda }$, представляет собой разность времен прихода прямого и отраженного от свободной поверхности лучей.Временной масштаб D изменчивости интерферограммы, как следует из (7), определяется условием
Аналитического решения эта задача не имеет, поэтому обратимся к тем результатам, которые можно извлечь из (12) в некоторых предельных случаях.
1. Пусть $r_{0}^{2} \gg {{\left( {wt} \right)}^{2}}$, т.е. рассматривается окрестность точки траверса, когда радиальная скорость равна нулю. Тогда, как следует из (12),
где $a = {{2r_{0}^{3}} \mathord{\left/ {\vphantom {{2r_{0}^{3}} {f{{w}^{2}}}}} \right. \kern-0em} {f{{w}^{2}}}}$. На практике можно ограничиться более слабым неравенством $r_{0}^{2} \geqslant 3{{\left( {wt} \right)}^{2}}$. Для условий эксперимента данное приближение реализуется уже при $t \leqslant 4$ с. В точке траверса $t = 0$ имеемДля полученных оценок (${{r}_{0}} = 189$ м, $\eta = 0.47$ м с), задаваясь частотой $f = 3.5$ кГц, имеем $D\left( 0 \right) = 3.3$ с, что близко к экспериментальному значению $D = 2.9$ с.
2. Пусть $r_{0}^{2} \ll {{\left( {wt} \right)}^{2}}$, т.е. расстояние от источника до приемника больше траверсного расстояния. В этом случае, согласно (12),
где $b = {w \mathord{\left/ {\vphantom {w {\eta f}}} \right. \kern-0em} {\eta f}}$. С ростом времени слагаемое bt растет и при значении ${{t}_{{{\text{cr}}}}} = {1 \mathord{\left/ {\vphantom {1 b}} \right. \kern-0em} b}$ временной масштаб становится бесконечным, так что локализация интерференционных полос производится вдоль горизонтальных линий. Предельному значению ${{t}_{{{\text{cr}}}}}$ отвечает удаление источника от антенны ${{r}_{{{\text{cr}}}}} = f$. Задаваясь значениями $\eta = 0.47$ м с, $f = 3.5$ кГц, получаем: ${{t}_{{{\text{cr}}}}} = 60$ с, ${{r}_{{{\text{cr}}}}} = 1.64$ км. Данный эффект наблюдается на экспериментальной (см. рис. 2а) и модельной (см. рис. 5) интерферограммах и полностью согласуется с известными характеристиками поля для модели, когда справедлив эффект Ллойда [20].Очевидным условием наблюдаемости интерферограммы является требование, чтобы ширина спектра $\Delta f$ хотя бы в несколько раз превышала частотный масштаб $\Lambda \left( t \right)$ (11). В качестве критерия наблюдаемости интерферограммы примем неравенство
которое эквивалентно условию, что наблюдается одна и более интерференционных полос. Ширина спектра ограничивает удаление источника от приемника, при этом максимальное удаление, соответствующее условию реализации одной полосы, оценивается как так что с увеличением полосы анализа и заглубления источника максимальное удаление возрастает. Для условий эксперимента ${{r}_{{{\text{max}}}}} = 3.29$ км.Можно показать также, что координаты максимума фокального пятна на голограмме (1), отвечающие действительному изображению источника, в точке траверса равны ${\tau } = {\eta \mathord{\left/ {\vphantom {\eta {{{r}_{0}}}}} \right. \kern-0em} {{{r}_{0}}}}$, $\nu = 0$. Для оценок $\eta = 0.47$ м с, ${{r}_{0}} = 189$ м получаем ${\tau } = 2.5$ мс, что близко к экспериментальному значению ${\tau } = 2.4$ мс (см. рис. 3).
Таким образом, модельная интерферограмма адекватно передает основные свойства экспериментальной интерферограммы до момента времени, когда источник начал совершать петлеобразные траектории.
ЗАКЛЮЧЕНИЕ
Изложен высокочастотный интерферометрический метод локализации шумового источника, использующий сигнал с одного из одиночных приемников антенны. В основе метода лежит частотно-временная обработка шумовых сигналов, согласованная с интерферограммой, формируемой движущимся источником. Модельная интерферограмма строилась на основе интерференции полей двух лучей: прямого и отраженного от свободной поверхности лучей. Метод апробирован на данных натурного эксперимента. По результатам обработки восстановлены временные зависимости пеленга, расстояния до антенны, скорости и глубины источника.
Двухлучевая модельная интерферограмма не требует знаний о параметрах дна и устойчива к поверхностному волнению, поэтому предложенная согласованная обработка оказалась работоспособной, позволив получить правильные временные зависимости параметров шумового источника. Совокупная картина фокальных пятен голограммы позволяет определять число лучей, формирующих поле, в различных акваториях. С удалением источника от приемника число лучей может увеличиваться, что будет приводить к усложнению интерференционной картины волнового поля. Если фокальные пятна не перекрываются, то на голограмме можно будет вырезать отдельные фокальные пятна и по отношению к ним реализовать обратное двукратное преобразование Фурье. Это даст возможность работать с двухлучевой интерферограммой, устойчивой по отношению к вариациям параметров дна.
Работа выполнена при финансовой поддержке грантов РФФИ (№ 19-08-00941, № 19-29-06075). Научно-исследовательская работа Д.Ю. Просовецкого поддержана грантом Президента РФ (№ МК-933.2019.8).
Список литературы
Корякин Ю.А., Смирнов С.А., Яковлев Г.В. Корабельная гидроакустическая техника: состояние и актуальные проблемы. СПб.: Наука, 2004. 410 с.
Кузнецов Г.Н., Кузькин В.М., Пересёлков С.А. Спектрограмма и локализация источника звука в мелком море // Акуст. журн. 2017. Т. 63. № 4. С. 406−418.
Казначеев И.В., Кузнецов Г.Н., Кузькин В.М., Пересёлков С.А. Интерферометрический метод обнаружения движущегося источника звука векторно-скалярным приемником // Акуст. журн. 2018. Т. 64. № 1. С. 33−45.
Kuz’kin V.M., Kuznetsov G.N., Pereselkov S.A., Grigor’ev V.A. Resolving power of the interferometric method of source localization // Phys. Wave Phenom. 2018. V. 26. № 2. P. 150−159.
Kuz’kin V.M., Pereselkov S.A., Kuznetsov G.N., Kaznacheev I.A. Interferometric direction finding by a vector-scalar receiver // Phys. Wave Phenom. 2018. V. 26. № 1. P. 63−73.
Kaznacheeva E.S., Kuznetsov G.N., Kuz’kin V.M., Lyakhov G.A., Pereselkov S.A. Measurement capability of the interferometric method of sound source localization in the absence of data on the waveguide transfer function // Phys. Wave Phenom. 2019. V. 27. № 1. P. 73–78.
Kuznetsov G.A., Kuz’kin V.M., Lyakhov G.A., Pereselkov S.A., Prosovetskiy D.Yu. Direction finding of a noise sound source // Phys. Wave Phenom. 2019. V. 27. № 3. P. 237–241.
Кузькин В.М., Пересёлков С.А., Рыбянец П.В., Ткаченко С.А. Реализация адаптивного интерферометрического метода локализации источника звука. Ч. I // Вестн. ВГУ. Серия: Физика. Математика. 2019. № 4. С. 19−29.
Чупров С.Д. Интерференционная структура звукового поля в слоистом океане / Акустика океана. Современное состояние. М.: Наука, 1982. С. 71–82.
Орлов Е.Ф. Интерференционная структура широкополосного звука в океане / Проблемы акустики океана. М.: Наука, 1984. С. 85–93.
Орлов Е.Ф., Шаронов Г.А. Интерференция звуковых волн в океане. Владивосток: Дальнаука, 1998. 195 с.
Baggeroer A.B., Kuperman W.A., Mikhalevsky P.N. An overview of matched field methods in ocean acoustics // IEEE. J. Oceanic Eng. 1993. V. 18. № 4. P. 401–423.
Сазонтов А.Г., Малеханов А.И. Согласованная пространственная обработка сигналов в подводных звуковых каналах (обзор) // Акуст. журн. 2015. Т. 61. № 2. С. 233–253.
Kuz’kin V.M., Pereselkov S.A., Zvyagin V.G., Malykhin A.Yu., Prosovetskiy D.Yu. Intense internal waves and their manifestation in interference patterns of received signals on oceanic shelf // Phys. Wave Phenom. 2018. V. 26. № 2. P. 160−167.
Badiey M., Kuz’kin V.M., Lyakhov G.A., Pereselkov S.A., Prosovetskiy D.Yu., Tkachenko S.A. Intense internal waves and their manifestation in interference patterns of received signals on oceanic shelf. Part II // Phys. Wave Phenom. 2019. V. 27. № 4. P. 313–319.
Вировлянский А.Л., Казарова А.Ю., Кенигсбергер Г.В., Колодиев О.В., Коротин П.И., Любавин Л.Я., Моисеенков В.И., Орлов Д.А., Потапов О.А., Турчин В.И. Эксперимент по оценке координат источника звука на шельфе черного моря // Акуст. журн. 2015. Т. 61. № 1. С. 1−9.
Аверьянов А.В., Глебова Г.М., Кузнецов Г.Н., Смирнов Н.М. Экспериментальная оценка пространственных координат источника шумового сигнала // Гидроакустика. 2013. Вып. 17 (1). С. 54–60.
Полканов К.И., Кузнецов Г.Н., Михнюк А.Н., Смирнов Н.М. Использование буксируемого векторно-скалярного модуля и согласованной фильтрации для однозначной оценки координат широкополосного источника в пассивном режиме // Гидроакустика. 2015. Вып. 24 (4). С. 36–51.
Бреховских Л.М., Лысанов Ю.П. Теоретические основы акустики океана. М.: Наука, 2007. 370 с.
Урик Р.Дж. Основы гидроакустики. Пер с англ. Л.: Судостроение, 1978. 448 p.
Дополнительные материалы отсутствуют.
Инструменты
Акустический журнал