способ определения коэффициента объемного рассеяния звука в океанической среде

Формула изобретения

1. СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ОБЪЕМНОГО РАССЕЯНИЯ ЗВУКА В ОКЕАНИЧЕСКОЙ СРЕДЕ, включающий направленный прием рассеянного звукового поля приемным трактом, установленным на погруженном движущемся аппарате, отличающийся тем, что в качестве звукового поля используют шумы моря при волнении моря два балла и более, осуществляют два вида пространственной фильтрации принятых шумов моря одновременно, посредством тракта с остронаправленной антенной формируют характеристику направленности с раствором главного максимума меньшим угла рефракционного минимума шумов моря в вертикальной плоскости, расположенным в пределах этого угла, и добавочными максимумами, при которых принятая этими добавочными максимумами энергия шумов по крайней мере в десять раз меньше энергии шума моря, принятого в пределах главного максимума, а также посредством приемного тракта со слабонаправленной антенной формируют характеристику направленности, близкую к полусфере с осью, направленной вперед по направлению движения аппарата в горизонтальной плоскости, и измеряют на выходах трактов с остронаправленной и слабонаправленной антеннами средние квадраты электрических напряжений и соответственно, предварительно с помощью тех же трактов производят прием шумов подводного аппарата на скоростях его движения тех же, что и приеме шумов моря, и при волнении моря менее одного балла, и измеряют средние квадраты электрического напряжения на выходе трактов с остронаправленной и слабонаправленной антенной и соответственно, а коэффициент объемного рассеяния определяют по формуле

_об= -1,

где - коэффициент затухания звука в океанической среде, табулированный для данной акватории;

K_пp , K_пс , _пp , _пс- - коэффициент концентрации и чувствительность трактов с остронаправленной и слабонаправленной соответственно.

2. Способ по п. 1, отличающийся тем, что прием шумов моря осуществляют в интервалах времени T, ограниченных неравенством

T<T,

где _сл - требуемая средняя квадратическая погрешность определения шумов моря;

T_ст - минимальный интервал времени, на котором шум моря является стационарным;

F - полоса частот, которое определяют коэффициент объемного рассеяния.

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

4. Способ по пп. 1 и 2, отличающийся тем, что дополнительно производят усреднение средних квадратов напряжений и , ранее усредненных в интервале времени T, по ансамблю n реализации за время T_общ = nТ.

Описание изобретения к патенту

Изобретение относится к акустике океана, а именно к гидроакустическим измерениям рассеяния звука в океанической среде.

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

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

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

W= R², где _c - акустическое сопротивление среды;

R - расстояние, и определение коэффициента обратного объема рассеивания как отношения мощности W, рассеянной статистически однородным объемом V, к произведению интенсивность звука I, падающего на рассеивающий объем V:

m_об= , после чего в предположении, что малые рассеивания, вычисляют коэффициент объемного рассеивания

_об= 4 m_об, имеющий размерность длины в минус первой степени (м^-1). Коэффициент объемного рассеивания, характеризующий мощность, рассеиваемую единичным объемом среды, понимается как среднее по многим неперекрывающимся объемам, выделенным в статистически однородной среде, при размерах этих объемов существенно больше радиуса пространственной корреляции неоднородностей [2] .

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

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

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

_об= -1, (1) где - коэффициент затухания звука в океанической среде, табулированный для данной акватории; К_пр, К_пс, _пр, _пс - коэффициенты концентрации и чувствительности трактов с остронаправленной и слабонаправленной антенной соответственно.

Прием шумов моря производят в интервалах времени Т, ограниченных неравенством

T<T, (2) где _cл - требуемая средняя квадратическая погрешность определения шумов моря; Т_ст - минимальный интервал времени, на котором шум моря является стационарным; F - полоса частот, в которой определяют коэффициент объемного рассеяния.

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

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

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

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

, (3) что позволяет извлечь информацию о коэффициенте объемного рассеивания _об.

Для определения лучевой интенсивности I необходимо отфильтровать шумы в области рефракционного минимума. Для этого необходимо применять направленную антенну с главным максимумом в вертикальной плоскости, не превышающим угол ₀= ₁+₂. Меньшая ширина главного максимума может быть использована, так как в пределах угла ₀ шум практически изотропен. Большая ширина недопустима, так как затрагивает области в вертикальной плоскости, в которых шум моря обусловлен нерассеиванием на объемным неоднородностях. Необходимый эффект достигается только при малых добавочных максимумах характеристики направленности приемной антенны. Уточним понятие малых добавочных максимумов. Условием можно считать

(360-₀)R , где R - средний уровень ореола добавочных максимумов.

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

R= .

На практике ₀ 30^о, I/I_экв 0,1, что требует среднего уровня добавочных максимумов меньше одного процента, что технически достижимо.

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

I= (4) где _c - среднее акустическое сопротивление морской среды.

Для определения I_экв необходимо осуществить всенаправленный прием шумов моря. Для этого используется примерно полусферическая характеристика направленности. Проведенное математическое моделирование показало, что для малонаправленного приема с коэффициентом концентрации от двух до трех и при главном максимуме характеристики направленности, ориентированном в горизонтальном направлении, коэффициент концентрации К_nс(t) совпадает с коэффициентом помехоустойчивости в анизотропном поле шумов моря с точностью до 5-7% . Это позволяет использовать коэффициент концентрации для определения I_экв. После измерения напряжения на выходе этого приемного тракта и компенсации помехи от шумов объектов появляется возможность определить эквивалентную лучевую интенсивность шумов моря как

I_экв= , (5) где U_оок, К_nс, _nc - соответственно уровень корабельных помех на выходе тракта, коэффициент концентрации и чувствительность. Коэффициент объемного рассеяния _об получают подстановкой в формулу (2) выражений (4) и (5).

Прием шумов моря целесообразно производить при

_сл= , причем Т не может быть меньше (при заданном _cл.) величины

T= , где _cл - требуемая средняя квадратическая погрешность определения шумов моря.

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

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

₀= ₁+₂, где

₁= arccos;

₂= arccos. На фиг. 1 показано распределение скоростей; на фиг. 2 - расчетное пространственное распределение интенсивности динамического шума моря в вертикальной плоскости (без учета объемного рассеивания в области рефракционного минимума): угол отсчитывается от горизонта, f_о= 1 кГц; на фиг. 3 - устройство для реализации способа, где 1 - электроакустический преобразователь, 2 - гидроакустическая антенна из n электроакустических преобразователей-приемников звука, 3 - n предварительных усилителей и фильтров, ограничивающих полосу F, 4 - формирователь характеристики направленности, 5 - первый аналого-цифровой преобразователь (АЦП), 6 - вычислитель - среднего квадрата напряжения на выходе формирователя характеристики направленности, 7 - вычислитель лучевой интенсивности I, 8 - вычислитель коэффициента объемного рассеивания _об, 9 - второй АЦП, 10 - вычислитель - среднего квадрата напряжения на выходе приемного тракта с полусферической характеристикой направленности, 11 - вычислитель I_экв - эквивалентной лучевой интенсивности, 12 - блок памяти помех аппарата-носителя гидролокатора - на выходе приемного остронаправленного тракта, 13 - блок памяти параметра антенны, отнесенного к акустическому сопротивлению среды, К_пр/_пр² 4 _c , 14 - блок памяти среднего квадрата напряжения помех аппарата-носителя гидролокатора на выходе приемного тракта с полусферической характеристикой направленности, 15 - блок памяти параметра приемника с полусферической характеристикой направленности, отнесенного к акустическому сопротивлению среды , 16 - блок памяти затухания звука в морской среде, 17 - блок управления устройством определения _об (на фиг. 3 для упрощения графика не показаны синхросвязи блока 17 управления с блоками устройства, в качестве слабонаправленной антенны используются один (на схеме четвертый) электроакустический преобразователь антенны); на фиг. 4 приведена блок-схема 17 управления, где 18 - генератор тактовых сигналов, 9, 20, 21, 22, 23 - соответственно первая, вторая, третья, четвертая, пятая линии задержки (цифры у выходов блока 17 следует считать связями).

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

I_об= e, (6) где Р_а - акустическая мощность излучателя;

с - скорость звука;

- продолжительность импульса;

_o- коэффициент, учитывающий влияние направленных свойств излучающей и приемной антенн на интенсивность реверберации.

В формуле (6) все параметры: акустическая мощность излучателя Р_а, коэффициент направленности _o, километрическое затухание звука , скорость звука с, длительность импульса , расстояние R заведомо известны. Необходимо только определить коэффициент объемного рассеяния _об.

В соответствии с предлагаемым способом измерено вертикальное распределение скорости звука с помощью скоростомера (фиг. 1). Было определено, что при глубине погружения аппарата-носителя гидролокатора 200 м, значение угла рефракионного минимума составляет _o = 30^о (от 75 до 105^о относительно направления на поверхность). Расчетные пространственные спектры без учета объемного рассеивания приведены на фиг. 2. Производился прием шума моря в полосе частот F= 100 Гц в течение интервалов времени Т= 60с. Выбор интервала Т сделан, исходя из неравенства (3) из заданного значения _сл= 1% . Была сформирована характеристика направленности приемной антенны с раствором 5^о на средней частоте работы гидролокатора f_о= 1200 Гц при расположении главного максимума характеристики направленности по горизонту. На выходе этого тракта измерялся средний квадрат напряжения, по формуле (4) была определена I. Предварительно были измерены помехи и , обусловленные движением и жизнеобитанием подводного аппарата. В соответствии с руководствами по нормированию и контролю этих помех измерения осуществлялись при волнении моря менее 1 балла, при удалении от берега не менее 30 миль в акватории с илистым грунтом глубиной не менее 200 м. Эти помехи определяются для разных скоростей хода аппарата. Одновременно формировалась характеристика направленности, близкая к полусфере (использовался один приемный элемент антенны), и была определена I_экв по формуле (5). По формуле (1) было найдено, что коэффициент объемной реверберации

_об= 10^-8 1/м.

Измерительные тракты предварительно были отградуированы по чувствительности методом сравнения с рабочим гидрофоном. Кроме того, были измерены характеристики направленности и по ним графическим методом определены значения коэффициентов концентрации.

Анализ погрешности показал, что случайная составляющая погрешности _сл при выбранном усреднении Т= 60 с пренебрежимо мала:

_сл= 1,5% , методические погрешности обусловлены принятыми предположениями о неизменности уровня помех от шумов аппарата, предположением о сферической индикатрисе рассеивания. Они характерны и для способа-прототипа и не превышают _м= 8% .

Кроме того, оценивались и погрешности _уи из-за изменения условий измерений, в частности при внезапном набеге шквала суммарная погрешность

= 16% .

Сопоставление определения _об по предлагаемому способу и по способу-прототипу показали удовлетворительное совпадение (20% ).

Таким образом, предлагаемый способ практически реализуем.

Реализация способа осуществлялась с помощью устройства, блок-схема которого приведена на фиг. 3. Устройство содержит параллельно-последовательно соединенные гидроакустическую антенну 2, состоящую из n электроакустических преобразователей 1, блок 3 предварительных усилителей и фильтров, блок 4 формирования характеристики направленности, первый АЦП 5, вычислитель 6 , вычислитель 7 I, вычислитель 8, а также содержит последовательно соединенные второй АЦП 9, вход которого соединен с выходом предварительного усилителя центрального в антенне электроакустического преобразователя, вычислитель 10 , вычислитель 11 I_экв, выход которого соединен с вторым входом вычислителя 8 _об. Кроме того, устройство содержит блок 12 памяти и блок 13 памяти , выходы которых соединены соответственно с вторым и третьим входами вычислителя 7, блок 14 памяти и блок 15 памяти , выходы которых соединены соответственно с вторым и третьим входами вычислителя 11 Iэкв, блок 16 памяти , выход которого соединен с третьим входом вычислителя _об, а также блок 17 управления, синхровходы которого соединены соответственно с синхровходами АЦП 5 и 9, вычислителей 6 и 10, блоков 12 и 13 и 14 и 15, вычислителей 7 и 11, блока 16, вычислителя 8.

Блок 17 управления содержит последовательно соединенные генератор 18 пеленговых сигналов, синхровыход которого соединен с синхровходами АЦП 5 и 9, первую линию 19 задержки, синхровыход которой соединен с синхровходами вычислителей 6 и 10, вторую линию 20 задержки, синхровыход которой соединен с синхровходами блоков 12, 13, 14 и 15, третью линию 21 задержки, синхровыход которой соединен с синхровходами вычислителей 7 и 11, четвертую линию 22 задержки, синхровыход которой соединен с синхровходом блока 16, пятую линию 23 задержек, синхровыход которой соединен с синхровходом вычислителя 8.

Работает устройство следующим образом.

По результатам измерений помех аппарата-носителя гидролокатора по правилам руководства по контролю и нормированию шумности в блоки 12 и 14 памяти заносятся данные о помехах-шумах аппарата для разных скоростей хода. В блоки 13 и 15 памяти заносятся паспортные параметры градуировки приемных трактов, а также средние значения _c для данной акватории. В блок 16 памяти заносятся табличные данные о затухании звука в море . В вычислителе 7 формируется характеристика направленности с малым уровнем добавочных максимумов обычно за счет использования Чебышевского распределения. Раствор главного максимума располагается по горизонту. Он выбирается заведомо уже, чем угол рефракционного минимума _o. В вычислителях 6 и 10 дискретизованные значения преобразуются в и .

В вычислителях 7 и 11 с учетом данных блоков 12, 13, 14, 15 памяти вычисляются I и I_экв. В вычислителе вычисляется _об с учетом данных о I и I_экв и значения для средней частоты гидролокации. Блок 17 управления обеспечивает синхронизацию работы всего устройства. Данные об _об подаются на блок прогнозирования дальности действия гидролокатора.