[go: up one dir, main page]

RU2740334C1 - Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation - Google Patents

Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation Download PDF

Info

Publication number
RU2740334C1
RU2740334C1 RU2020123339A RU2020123339A RU2740334C1 RU 2740334 C1 RU2740334 C1 RU 2740334C1 RU 2020123339 A RU2020123339 A RU 2020123339A RU 2020123339 A RU2020123339 A RU 2020123339A RU 2740334 C1 RU2740334 C1 RU 2740334C1
Authority
RU
Russia
Prior art keywords
receiver
components
hydroacoustic
component
wave
Prior art date
Application number
RU2020123339A
Other languages
Russian (ru)
Inventor
Владимир Ильич Коренбаум
Сергей Владимирович Горовой
Вениамин Мефодьевич Дорожко
Алексей Евгеньевич Бородин
Original Assignee
Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) filed Critical Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН)
Priority to RU2020123339A priority Critical patent/RU2740334C1/en
Application granted granted Critical
Publication of RU2740334C1 publication Critical patent/RU2740334C1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/162Details
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/181Geophones
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/186Hydrophones
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/189Combinations of different types of receiving elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808Seismic data acquisition, e.g. survey design

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

FIELD: physics.SUBSTANCE: invention relates to hydroacoustics and geophysics and can be used to detect sound sources or obtain information on the structure and characteristics of the sea bottom (soil) by analyzing acoustic waves propagating in the ocean (geoacoustic inversion). Method involves directed reception of aggregate wave coaxially located inside a single rigid housing by 3-component vibration receiver and located on its surface of 3-component pressure gradient receiver of power type, equipped with acoustic pressure receiver, with subsequent joint processing of responses of coaxial components of 3-component receivers of pressure gradient and receiver of vibrations with extraction of hydroacoustic component.EFFECT: extraction of a hydroacoustic component from a set of seismic-acoustic and hydroacoustic waves affecting a receiving device, which further improves the quality of detecting a radiation source or estimating parameters of the bottom.5 cl, 1 dwg

Description

Изобретение относится к гидроакустике и геофизике и может быть использовано для обнаружения источника звука или получения информации о структуре и характеристиках морского дна (грунта) путем анализа акустических волн, распространяющихся в океане (геоакустическая инверсия).The invention relates to hydroacoustics and geophysics and can be used to detect a sound source or obtain information on the structure and characteristics of the seabed (soil) by analyzing acoustic waves propagating in the ocean (geoacoustic inversion).

При обнаружении источников низкочастотных звуков в море и других акваториях часть низкочастотного излучения источника передается гидроакустической плоской волной, бегущей в водной среде, а часть может передаваться за счет различных механизмов в осадочный слой дна и возбуждать в нем донные сейсмоакустических волны (в виде нормальных волн, бегущих в осадочном слое, и/или вытекающих в воду, например, волн Стоунли) [Bevans D.A., Buckingham M.J. Estimating the sound speed of a shallow-water marine sediment from the head wave excited by a low-flying helicopter // J. Acoust. Soc. Am. 2017. 142 (4), P. 2273-2287. doi: 10.1121/1.5007953]. Обе составляющие волны - гидроакустическая и сейсмоакустическая - могут распространяться достаточно далеко и обнаруживаться различными приемными устройствами. Однако накладываясь друг на друга они искажают излученный сигнал, что препятствует обнаружению источника или получению информации о структуре и характеристиках морского дна.When low-frequency sound sources are detected in the sea and other water areas, part of the low-frequency radiation of the source is transmitted by a hydroacoustic plane wave traveling in the aquatic environment, and part can be transmitted through various mechanisms into the sedimentary layer of the bottom and excite bottom seismoacoustic waves in it (in the form of normal waves traveling in the sedimentary layer, and / or flowing into the water, for example, Stoneley waves) [Bevans DA, Buckingham MJ Estimating the sound speed of a shallow-water marine sediment from the head wave excited by a low-flying helicopter // J. Acoust. Soc. Am. 2017.142 (4), P. 2273-2287. doi: 10.1121 / 1.5007953]. Both components of the wave - hydroacoustic and seismoacoustic - can propagate far enough and be detected by various receivers. However, when superimposed on each other, they distort the emitted signal, which prevents the source from being detected or from obtaining information about the structure and characteristics of the seabed.

Актуальной задачей является разделение сейсмоакустической и гидроакустической составляющих волны, позволяющее повысить качество обнаружения источника или оценки параметров дна водоема.An urgent task is to separate the seismoacoustic and hydroacoustic components of the wave, which makes it possible to improve the quality of source detection or the estimation of the parameters of the bottom of the reservoir.

Прием совокупности сесмоакустической и гидроакустической волн обычно производится как углубленными в осадочный слой приемниками вибраций (приемниками колебательного смещения, колебательной скорости или колебательного ускорения) [Shi J., Dosso S.E., Sun D., Qingyu L. Geoacoustic inversion of the acoustic-pressure vertical phase gradient from a single vector sensor // J. Acoust. Soc. Am. 2019. 146 (5), P. 3159-3173. doi: 10.1121/1.5131235], так и расположенными в воде у дна приемниками звукового давления - гидрофонами [Gervaise С., Kinda B.G., Bonnel J., Stephan Y., Vallez S. Passive geoacoustic inversion with a single hydrophone using broadband ship noise // J. Acoust. Soc. Am. 2012. 131 P. 1999-2010]).Reception of a combination of sismoacoustic and hydroacoustic waves is usually performed as receivers of vibrations deepened into the sedimentary layer (receivers of vibrational displacement, vibrational velocity or vibrational acceleration) [Shi J., Dosso SE, Sun D., Qingyu L. Geoacoustic inversion of the acoustic-pressure vertical phase gradient from a single vector sensor // J. Acoust. Soc. Am. 2019.146 (5), P. 3159-3173. doi: 10.1121 / 1.5131235], and located in the water at the bottom of the sound pressure receivers - hydrophones [Gervaise S., Kinda BG, Bonnel J., Stephan Y., Vallez S. Passive geoacoustic inversion with a single hydrophone using broadband ship noise / / J. Acoust. Soc. Am. 2012. 131 P. 1999-2010]).

Известны способ и устройство, в которых волны принимают приемниками вибраций (колебательной скорости) и звукового давления, размещенными над и под поверхностью дна водоема и рассчитывают передаточные функции между датчиками, размещенными над и под поверхностью дна водоема (Crocker S.E., Miller J.H., Potty G.R., Osier J.C., Hines P.С.Nonlinear inversion of acoustic scalar and vector field transfer function // IEEE J. Oceanic Eng. 2012. 37. P. 589-607).There is a known method and device in which the waves are received by receivers of vibrations (oscillatory speed) and sound pressure placed above and below the surface of the bottom of the reservoir and calculate the transfer functions between the sensors located above and below the surface of the bottom of the reservoir (Crocker SE, Miller JH, Potty GR, Osier JC, Hines P.S. Nonlinear inversion of acoustic scalar and vector field transfer function // IEEE J. Oceanic Eng. 2012.37.P. 589-607).

Однако данные описанные способы и реализующие их устройства не позволяют разделить сейсмоакустическую и гидроакустическую составляющие волны.However, these described methods and the devices implementing them do not allow separating the seismoacoustic and hydroacoustic wave components.

Наиболее близким к заявляемому относится устройство для приема совокупности сейсмоакустической и гидроакустической волны, представляющий собой смешанный приемник с единым жестким корпусом, включающим в себя одновременно приемник вибраций (геофон) и приемник звукового давления (гидрофон), при этом в качестве приемника вибраций используется датчик ускорения с инерционной массой, вокруг которого коаксиально расположен приемник звукового давления. С приемником вибраций жестко связан приемный сейсмический штырь. В зависимости от условий эксплуатации прием совокупности сейсмоакустической и гидроакустической составляющих волны выполняют или приемником вибраций (геофоном) или ненаправленным приемником звукового давления (гидрофоном) в зависимости от условий эксплуатации, однако отклики датчиков совместно не обрабатывают с целью выделения сейсмоакустической и гидроакустической составляющих волны (п. FR №2543692).The closest to the claimed device is a device for receiving a set of seismoacoustic and hydroacoustic waves, which is a mixed receiver with a single rigid body, including both a vibration receiver (geophone) and a sound pressure receiver (hydrophone), while an acceleration sensor with the inertial mass, around which the sound pressure receiver is coaxially located. A receiving seismic pin is rigidly connected to the vibration receiver. Depending on the operating conditions, the reception of a set of seismoacoustic and hydroacoustic wave components is performed either by a vibration receiver (geophone) or a non-directional sound pressure receiver (hydrophone), depending on the operating conditions, however, the sensor responses are not jointly processed in order to isolate the seismoacoustic and hydroacoustic wave components (p. FR No. 2543692).

Таким образом, с помощью известного геогидрофона возможно работать с одним и тем же прибором, как в качестве приемника звукового давления, так и в качестве приемника вибраций в зависимости от условий его эксплуатации, то есть так же, как и в приведенных выше аналогах, без какого-либо совместного включения выходов датчиков, то есть без возможности разделения сейсмоакустической и гидроакустической составляющих волны.Thus, with the help of the known geohydrophone, it is possible to work with the same device, both as a sound pressure receiver and as a vibration receiver, depending on the conditions of its operation, that is, the same as in the above analogs, without any - or joint switching on of the sensor outputs, that is, without the possibility of separating the seismoacoustic and hydroacoustic wave components.

Проблема требующая решения - разделение сейсмоакустической и гидроакустической составляющих совокупной волны.The problem to be solved is the separation of the seismoacoustic and hydroacoustic components of the aggregate wave.

Для ее решения предлагается способ приема совокупности сейсмоакустической и гидроакустической составляющих волны у сложенного осадочными породами грунта дна водоема, включающий направленный прием совокупной волны одновременно соосно расположенным внутри единого жесткого корпуса 3-х компонентным приемником вибраций и расположенным на его поверхности 3-х компонентным приемником градиента давления силового типа,, снабженным приемником звукового давления, совместную обработку откликов соосных компонент 3-х компонентных приемников градиента давления и приемника вибраций с выделением гидроакустической составляющей путем вычисления разности откликов компонент приемника градиента давления и откликов соосных компонент приемника вибраций с предварительной настройкой коэффициента усиления каналов приемника вибраций на максимум подавления модельных вибраций, при этом приемник вибраций размещен в центре масс корпуса, расположенного в грунте, а приемник градиента давления с приемником звукового давления - на поверхности корпуса в воде над грунтом.To solve it, a method is proposed for receiving a combination of seismoacoustic and hydroacoustic components of a wave at the bottom of a reservoir folded by sedimentary rocks of the soil, including directional reception of the aggregate wave simultaneously coaxially located inside a single rigid body with a 3-component vibration receiver and a 3-component receiver of pressure gradient located on its surface power type, equipped with a sound pressure receiver, joint processing of the responses of the coaxial components of the 3-component pressure gradient receivers and the vibration receiver with the extraction of the hydroacoustic component by calculating the difference in the responses of the pressure gradient receiver components and the responses of the coaxial components of the vibration receiver with presetting the amplification factor of the vibration receiver channels to the maximum suppression of model vibrations, while the vibration receiver is located in the center of mass of the housing located in the ground, and the pressure gradient receiver with the sound pressure receiver enia - on the surface of the hull in water above the ground.

При необходимости в дальнейшем можно разделить выделенную гидроакустическую волну на, преимущественно, плоскую и неоднородную составляющие, применив интенсиметрическую мультипликативную обработку откликов. Для этого осуществляют операцию усиления разности откликов соосных компонент приемников градиента давления и вибраций с предварительно определенным коэффициентом усиления и вычисляют взаимный спектр полученной разности откликов по каждой компоненте с откликом приемника звукового давления, а из вычисленного взаимного спектра выделяют вещественную часть, характеризующую плоскую гидроакустическую волну от удаленного источника, и мнимую часть, характеризующую неоднородную гидроакустическую волну.If necessary, in the future, it is possible to divide the selected hydroacoustic wave into mainly flat and inhomogeneous components, using intensimetric multiplicative processing of responses. To do this, the operation of amplifying the difference in the responses of the coaxial components of the pressure gradient and vibration receivers with a predetermined gain is carried out and the mutual spectrum of the obtained difference in responses is calculated for each component with the response of the sound pressure receiver, and from the calculated mutual spectrum, the real part characterizing a plane hydroacoustic wave from the remote source, and an imaginary part characterizing an inhomogeneous hydroacoustic wave.

Проблема также решается устройством приема совокупности сейсмоакустической и гидроакустической волн у дна водоема, сложенного осадочными породами грунта, состоящим из единого жесткого корпуса, со средней плотностью близкой к плотности грунта, в центре масс которого размещен 3-х компонентный приемник вибраций, а на поверхности соосно компонентам приемника вибраций установлен 3-х компонентный приемник градиента давления силового типа, дополнительно снабженный приемником звукового давления, установленным с обеспечением единого фазового центра с компонентами приемника градиента давления, при этом электрические выходы соосных компонент приемника градиента давления и приемника вибраций соединены через усилители с блоками вычитания, а выходы блоков вычитания через усилители совместно с выходами приемника звукового давления присоединены к блоку обработки данных.The problem is also solved by a device for receiving a set of seismoacoustic and hydroacoustic waves at the bottom of a reservoir composed of sedimentary rocks of the soil, consisting of a single rigid body, with an average density close to the density of the soil, in the center of mass of which there is a 3-component vibration receiver, and on the surface it is coaxial to the components of the vibration receiver, a 3-component force-type pressure gradient receiver is installed, additionally equipped with a sound pressure receiver installed to ensure a single phase center with the components of the pressure gradient receiver, while the electrical outputs of the coaxial components of the pressure gradient receiver and vibration receiver are connected through amplifiers with subtraction units, and the outputs of the subtraction units through amplifiers, together with the outputs of the sound pressure receiver, are connected to the data processing unit.

Технический результат совместного использования заявляемых решений -выделение гидроакустической составляющей из совокупности сейсмоакустической и гидроакустической волны, воздействующей на приемное устройство, что позволяет в дальнейшем при необходимости использовать полученный результат для последующего выделения из гидроакустической волны преимущественно плоской (бегущей от удаленного источника звука) и неоднородной составляющей, что в дальнейшем совместно с сейсмоакустической составляющей принимаемого совокупного сигнала может быть использовано, в том числе, и для обнаружения/пеленгования удаленного источника звука или получения информации о структуре и характеристиках дна водоема (геоакустическая инверсия) известными методамиThe technical result of the joint use of the proposed solutions is the separation of the hydroacoustic component from the totality of the seismoacoustic and hydroacoustic waves acting on the receiving device, which allows, if necessary, to use the obtained result for the subsequent separation from the hydroacoustic wave of a predominantly flat (traveling from a distant sound source) and inhomogeneous component, that in the future, together with the seismoacoustic component of the received aggregate signal, it can be used, inter alia, for detection / direction finding of a remote sound source or obtaining information about the structure and characteristics of the bottom of a reservoir (geoacoustic inversion) by known methods

На Фиг. представлена конфигурация заявляемого устройства и его размещение на дне акватории, где 1 - донные осадочные породы грунта дна акватории, 2 - жесткий корпус, 3-3-х компонентный приемник вибраций, 4-3-х компонентный приемник градиента давления (ПГД) силового типа с встроенным каналом приемника звуквого давления; pwa, vwa - звуковое давление и колебательная скорость в плоской гидроакустической волне, распространяющеся в водной среде; asa - три компоненты ускорения, вызываемые сейсмоакустической волной, распространяющейся в грунте; pwS, vwS - звуковое давление и колебательная скорость неоднородной гидроакустической волны, переизлучаемой сейсмоакустической волной из грунта в воду, в частности, например, волна Стоунли.FIG. shows the configuration of the claimed device and its placement at the bottom of the water area, where 1 - bottom sedimentary rocks of the bottom of the water area, 2 - rigid body, 3-3-component vibration receiver, 4-3-component pressure gradient receiver (PGD) of power type with built-in sound pressure receiver channel; p wa , v wa - sound pressure and vibrational velocity in a plane hydroacoustic wave propagating in an aqueous medium; a sa - three components of acceleration caused by a seismoacoustic wave propagating in the ground; p wS , v wS - sound pressure and vibrational velocity of an inhomogeneous hydroacoustic wave re-emitted by a seismoacoustic wave from the ground into water, in particular, for example, a Stoneley wave.

Работоспособность способа и устройства подтверждается теоретической моделью его функционирования.The efficiency of the method and device is confirmed by the theoretical model of its functioning.

В предположении, что корпус 2 существенно меньше длины волны и практически полностью находится в грунте, на него действуют лишь сейсмоакустические волны (как продольные, так и поперечные), но не действует плоская гидроакустическая волна, бегущая от удаленного источника. В данной системе (Фиг) мы имеем следующие отклики датчиков.On the assumption that body 2 is significantly less than the wavelength and is almost completely in the ground, only seismic acoustic waves (both longitudinal and transverse) act on it, but a plane hydroacoustic wave traveling from a distant source does not act. In this system (Fig), we have the following sensor responses.

Электрический отклик 3-х компонентного приемника вибраций 3 в виде акселерометра по ускорению может быть записан в виде:The electrical response of a 3-component vibration receiver 3 in the form of an acceleration accelerometer can be written in the form:

Uacsaacsaasa, где νacsa - чувствительность акселерометра по ускорению, иU acsa = ν acsa a sa , where ν acsa is the accelerometer acceleration sensitivity, and

asa - ускорение, вектора с компонентами (ах, ау, az), x,y,z - Декартовы координаты.a sa - acceleration, vectors with components ( ax , y , a z ), x, y, z - Cartesian coordinates.

В случае гармонических сигналов этот электрический отклик акселерометра может быть преобразован к формеIn the case of harmonic signals, this electrical response of the accelerometer can be converted to the form

Figure 00000001
Figure 00000001

где i - мнимая единица, ω - круговая частота сигнала.where i is the imaginary unit, ω is the circular frequency of the signal.

При распространении в осадочных породах продольной сейсмоакустической составляющей волны по значению отношений полученных компонент (νacsaxax, νacsayay, νacsazaz) может быть известным образом получена информация о направлении ее прихода, т.е. направлении на источник сигнала. Однако, при распространеии в осадочных породах поперечной сейсмоакустической волны вся информация о направлении ее прихода теряется.When the longitudinal seismoacoustic component of the wave propagates in sedimentary rocks, information on the direction of its arrival can be obtained in a known manner by the value of the ratios of the obtained components (ν acsax a x , ν acsay a y , ν acsaz a z ), i.e. direction to the signal source. However, when a transverse seismoacoustic wave propagates in sedimentary rocks, all information about the direction of its arrival is lost.

С другой стороны, электрический отклик 3-х компонентного ПГД силового типа 4 складывается из нескольких составляющихOn the other hand, the electrical response of a 3-component PGD of power type 4 consists of several components

Figure 00000002
Figure 00000002

где νpgsa - чувствительность ПГД по ускорению, νpgsv - чувствительность ПГД по колебательной скорости,where ν pgsa is the acceleration sensitivity of PGD, ν pgsv is the oscillatory velocity sensitivity of PGD,

asa/iω вектор, который соответствует колебательной скорости в сейсмоакустической волне,a sa / iω vector that corresponds to the vibrational velocity in the seismoacoustic wave,

vwa вектор, который соответствует колебательной скорости в плоской гидроакустической волне,v wa vector, which corresponds to the vibrational velocity in a plane hydroacoustic wave,

vwS - вектор, который соответствует колебательной скорости неоднородной гидроакустической волны, вытекающей из грунта в воду.v wS is a vector that corresponds to the vibrational velocity of an inhomogeneous hydroacoustic wave flowing from the ground into the water.

Конструируя разность откликов Upgsv - Uacsv, получаем: By constructing the difference of responses U pgsv - U acsv , we obtain:

Figure 00000003
Figure 00000003

В случае выполнения равенства νpgsaacsa (по каждой из 3-х компонент обоих приемников, что достигается предварительной настройкой коэффициента усиления акселерометра νacsa на максимум подавления модельных вибраций при подвешивании устройства на вибростоле над водой в каждом из направлений компонент приемников градиента давления и вибрации, например, как предложено в п. РФ №2624791, член в круглых скобках (3) обнуляется, и само выражение (3) редуцируется к видуIf the equality ν pgsa = ν acsa (for each of the 3 components of both receivers is fulfilled , which is achieved by presetting the accelerometer gain ν acsa to the maximum suppression of model vibrations when the device is suspended on a vibration table above water in each direction of the components of the pressure gradient receivers and vibration, for example, as suggested in p. RF No. 2624791, the term in parentheses (3) is reset to zero, and expression (3) itself is reduced to the form

Figure 00000004
Figure 00000004

Таким образом, полученная разность (4) описывает отклик только на гидроакустическую составляющую совокупной волны, падающей на устройство.Thus, the obtained difference (4) describes the response only to the hydroacoustic component of the total wave incident on the device.

В то же время, известно, что в гидроакустической составляющей имеются две волны - плоская (первое слагаемое) и неоднородная, вытекающая из дна (второе слагаемое). При этом канал звукового давления, совмещенный с ПГД, будучи нечувствителен к вибрациям, регистрирует отклик скалярных величин:At the same time, it is known that there are two waves in the hydroacoustic component - plane (first term) and inhomogeneous, flowing from the bottom (second term). In this case, the sound pressure channel, combined with the PGD, being insensitive to vibrations, registers the response of scalar quantities:

Figure 00000005
Figure 00000005

Для разделения плоской гидроакустической волны, описываемой параметрами pwa, vwa и неоднородную вытекающую гидроакустическую волну, описываемую параметрами pwS, vwS применим интенсиметрическую мультипликативную обработку откликов. Известно, что выделение вещественной части взаимного спектра (потока мощности) через быстрое преобразование Фурье (FFT) позволяет подавить неоднородные гидроакустические волны (Korenbaum VI, Tagiltsev АА. Flow noise of an underwater vector sensor embedded in a flexible towed array // J. Acoust. Soc. Am. 2012. 131(5):3755-3762. DOI: 10.1121/1.3693647) и вычислить спектры пространственных компонент вектора активной части потока мощности (интенсивности) плоской гидроакустической волны. Таким образом, получимTo separate the plane hydroacoustic wave, described by the parameters p wa , v wa, and the inhomogeneous flowing hydroacoustic wave, described by the parameters p wS , v wS, we use the intensimetric multiplicative processing of the responses. It is known that the separation of the real part of the mutual spectrum (power flow) through the fast Fourier transform (FFT) allows suppressing inhomogeneous hydroacoustic waves (Korenbaum VI, Tagiltsev AA. Flow noise of an underwater vector sensor embedded in a flexible towed array // J. Acoust. Soc. Am. 2012. 131 (5): 3755-3762. DOI: 10.1121 / 1.3693647) and calculate the spectra of the spatial components of the vector of the active part of the power flux (intensity) of a plane hydroacoustic wave. Thus, we get

Figure 00000006
Figure 00000006

Figure 00000007
Figure 00000007

Соответственно, наоборот, вычисление мнимой части взаимного спектра позволяет выделить спектры пространственных компонент реактивной части потока мощности:Accordingly, on the contrary, the calculation of the imaginary part of the mutual spectrum makes it possible to single out the spectra of the spatial components of the reactive part of the power flow:

Figure 00000008
Figure 00000008

Таким образом, с помощью выражения (1) может быть зарегистрирована преимущественно сейсмоакустическая составляющая совокупной волны, с помощью выражениий (3, 4) - преимущественно гидроакустическая составляющая совокупной волны, с помощью выражения (6) - преимущественно плоская гидроакустическая волна, а с помощью выражения (7) - преимущественно неоднородная гидроакустическая волна.Thus, using expression (1), the predominantly seismoacoustic component of the aggregate wave can be recorded, using expressions (3, 4), the predominantly hydroacoustic component of the aggregate wave, using expression (6), a predominantly plane hydroacoustic wave, and using the expression ( 7) - predominantly inhomogeneous hydroacoustic wave.

Таким образом, могут быть разделены все гидроакустические и сейсмоакустические составляющие принимаемого совокупного сигнала. А это может быть использовано для обнаружения/пеленгования удаленного источника звука или получения информации о структуре и характеристиках дна водоема (геоакустическая инверсия) известными методами.Thus, all hydroacoustic and seismoacoustic components of the received aggregate signal can be separated. And this can be used to detect / find a distant sound source or obtain information about the structure and characteristics of the bottom of a reservoir (geoacoustic inversion) using known methods.

По уже выделенной из гидроакустической составляющей - плоской звуковой волне (6) - в каждом из направлений компонент приемника градиента давления по соотношению амплитуд откликов в этих компонентах, используя дополнительно канал приемника звукового давления для устранения неоднозначности пеленгования, можно определить направление на источник сигнала известным векторно-фазовым методом. Таким образом, осуществив операцию вычисления вещественной части взаимного спектра по направлениям каждой из соосных компонент приемников градиента давления и вибраций определяют пеленг на источник гидроакустической волны по отношению амплитуд компонент вещественной части взаимного спектра.Already isolated from the hydroacoustic component - a plane sound wave (6) - in each direction of the components of the pressure gradient receiver according to the ratio of the amplitudes of the responses in these components, using the additional channel of the sound pressure receiver to eliminate the ambiguity of direction finding, it is possible to determine the direction to the signal source by a known vector phase method. Thus, having performed the operation of calculating the real part of the mutual spectrum in the directions of each of the coaxial components of the receivers of the pressure gradient and vibrations, the bearing to the source of the hydroacoustic wave is determined in relation to the amplitudes of the components of the real part of the mutual spectrum.

Аналогично, для повышения помехоустойчивости приема гидроакустической части сигнала от удаленного источника или сканирования по пространству в поисках источника сигнала, по уже выделанной разности откликов, характеризующей составляющие гидроакустического сигнала, могут формироваться однонаправленные диаграммы направленности типа кардиоидной или произведения косинусоиды на кардиоиду, например, как в п. РФ 2687301. Для этого осуществляют операцию усиления разности откликов с предварительно определенным коэффициентом усиления разности, и совместно с откликом приемника звукового давления формируют однонаправленные диаграммы направленности типа кардиоидной или произведения косинусоиды на кардиоиду в направлении каждой из компонент приемников градиента давления и вибраций.Similarly, to increase the noise immunity of receiving the hydroacoustic part of the signal from a remote source or scanning through space in search of a signal source, using the already identified difference in responses characterizing the components of the hydroacoustic signal, unidirectional radiation patterns such as a cardioid or a product of a cosine and a cardioid can be formed, for example, as in p. RF 2687301. For this, the operation of amplifying the difference in responses with a predetermined gain of the difference is carried out, and together with the response of the sound pressure receiver, unidirectional radiation patterns of the cardioid type or the product of the cosine wave and the cardioid are formed in the direction of each of the components of the pressure gradient and vibration receivers.

Можно также для повышения эффективности обнаружения удаленного источника сигнала попытаться суммировать сфазированно все выделенные составляющие волны с регулируемой задержкой во времени. При этом задержки во времени можно определять по одному из известных критериев оптимального обнаружения сигнала на фоне помех.You can also try to sum up in phased all the selected wave components with an adjustable time delay to improve the efficiency of detecting a remote signal source. In this case, time delays can be determined by one of the known criteria for optimal signal detection against the background of interference.

Для задач геоакустической инверсии может оказаться полезной информация о неоднородной вытекающей гидроакустической волне, получаемой в соответствии с выражением (7). По соотношению амплитуд сигналов в каналах x,y,z компонент можно получить информацию о вертикальном угле вытекания этой волны из дна, который определяет сейсмоакустические свойства осадочного слоя грунта.For the tasks of geoacoustic inversion, information on the inhomogeneous flowing hydroacoustic wave obtained in accordance with expression (7) may be useful. By the ratio of signal amplitudes in the channels of the x, y, z components, one can obtain information about the vertical angle of the outflow of this wave from the bottom, which determines the seismoacoustic properties of the sedimentary layer of the soil.

При сейсмоакустическом зондировании дна с наклонных направлений близких к вертикальному, как правило, имеет место три акустических волны: гидроакустические падающая, отраженная и сейсмоакустическая прошедшая. Прошедшая сейсмоакустическая волна вызывает колебания корпуса устройства, поэтому для выделения падающей и отраженной волн неоходимо использовать разности откликов компонент приемников градиента давления и приемника вибраций. Только для выделения падающей волны следует из указанной разности откликов вертикальных компонент ПГД и приемника вибраций, а также отклика канала звукового давления, сформировать кардиоидную характеристику направленности ориентированную максимумом вверх, а нулем - вниз. Таким образом, отраженная гидроакустическая волна в значительной степени подавляется.During seismoacoustic sounding of the bottom from inclined directions close to vertical, as a rule, there are three acoustic waves: hydroacoustic incident, reflected and seismoacoustic transmitted. The transmitted seismoacoustic wave causes vibrations of the device body; therefore, to separate the incident and reflected waves, it is necessary to use the differences in the responses of the components of the pressure gradient receivers and the vibration receiver. Only to isolate the incident wave, it follows from the indicated difference in the responses of the vertical components of the PGD and the vibration receiver, as well as the response of the sound pressure channel, to form a cardioid directivity characteristic oriented with its maximum upward, and zero - downward. Thus, the reflected hydroacoustic wave is largely suppressed.

Тогда как для выделения отраженной от дна гидроакустической волны необходимо сформировать противоположным образом ориентированную кардиоидную характеристику направленности из тех же откликов. В этом случае, напротив, падающая гидроакустическая волна подавляется.Whereas, in order to isolate the hydroacoustic wave reflected from the bottom, it is necessary to form an oppositely oriented cardioid directivity characteristic from the same responses. In this case, on the contrary, the incident hydroacoustic wave is suppressed.

Прошедшая же в дно сейсмоакустическая волна, если она продольная, определяется по отклику вертикальной компоненты приемника вибраций, а, если поперечная, то по откликам горизонтальных компонент приемника вибраций. При обеспечении соответствующих коэффициентов усиления, выравнивающих чувствительности датчиков, теперь можно точно определить коэффициенты отражения и прохождения зондирующей волны. Причем в данном случае зондирование может производится не только импульсным, но и непрерывным, в том числе шумовым, сигналом.The seismoacoustic wave that has passed into the bottom, if it is longitudinal, is determined from the response of the vertical component of the vibration receiver, and if it is transverse, then from the responses of the horizontal components of the vibration receiver. By providing appropriate gains that equalize the sensitivities of the sensors, you can now accurately determine the reflection and transmission coefficients of the probe wave. Moreover, in this case, sounding can be performed not only by impulse, but also by continuous, including noise, signal.

Заявляемое изобретение проиллюстрируем следующим примером.The claimed invention will be illustrated by the following example.

В корпусе 2, (выполненном, например, из морского бетона), имеющем среднюю плотность, соответствующую плотности грунта 1 донных осадочных пород, устанавливают в центре масс в качестве приемника вибраций трехкомпонентный акселерометр 3, а сверху на оси симметрии соосный компонентам акселерометра трехкомпонентный ПГД силового типа с каналом звукового давления, например, в виде крестообразного устройства 4, описанного в патенте РФ №2679931. При этом корпус 2 устанавливают днищем в осадочные породы грунта 1 с углублением, но так, чтобы ПГД оставался в водной среде (Фиг).In case 2, (made, for example, of sea concrete), having an average density corresponding to the density of the soil 1 of bottom sedimentary rocks, a three-component accelerometer 3 is installed in the center of mass as a vibration receiver, and on top of the symmetry axis, a three-component power-type PGD is coaxial to the accelerometer components with a sound pressure channel, for example, in the form of a cruciform device 4 described in RF patent No. 2679931. In this case, the body 2 is installed with its bottom in the sedimentary rocks of the soil 1 with a depression, but so that the PGD remains in the aquatic environment (Fig).

Для приема сейсмоакустических волн, распространяющихся в осадочных породах дна используют компоненты акселерометра 3, тогда как для приема гидроакустических волн, распространяющихся в водной среде используют разность откликов компонент ПГД 4 и акселерометра 3, формируемую в блоке вычитания (выполненном, например, на операционном усилителе, расположенном в блоке обработки, согласно выражениям (3, 4), предварительно отрегулировав равенство их чувствительностей по виброускорению νpgsaacsa по каждой из 3-х компонент обоих датчиков. Это достигается, например, предварительной настройкой коэффициента усиления акселерометра νacsa на максимум подавления модельных вибраций при подвешивании устройства на вибрирующем вибростоле над водой с ориентацией в каждом из направлений компонент приемников градиента давления и вибрации (например, п. РФ №2624791. Данное включение обеспечивает компенсацию откликов компонент ПГД на колебания единого жесткого корпуса, вызванные сейсмоакустической волной.Accelerometer 3 components are used to receive seismoacoustic waves propagating in sedimentary rocks of the bottom, while for receiving hydroacoustic waves propagating in an aquatic environment, the difference in responses of the components of PGD 4 and accelerometer 3 is used, formed in the subtraction unit (made, for example, on an operational amplifier located in the processing unit, according to expressions (3, 4), having previously adjusted the equality of their vibration acceleration sensitivities ν pgsa = ν acsa for each of the 3 components of both sensors.This is achieved, for example, by presetting the accelerometer gain ν acsa to the maximum suppression of model vibrations when the device is suspended on a vibrating vibration table above the water with orientation in each direction of the components of the pressure gradient and vibration receivers (for example, p. RF No. 2624791. This inclusion provides compensation for the responses of the PGD components to the vibrations of a single rigid body caused by a seismoacoustic wave.

При этом для отделения плоской гидроакустической волны от вытекающей части неоднородной волны, например, волны Стоунли, вычисляют вещественную часть взаимного спектра (поток мощности) по разности откликов компонент ПГД и компонент акселерометра (после интегрирования) и канала звукового давления согласно выражению (6).In this case, to separate a plane hydroacoustic wave from the outgoing part of an inhomogeneous wave, for example, a Stoneley wave, the real part of the mutual spectrum (power flux) is calculated from the difference between the responses of the PGD components and the components of the accelerometer (after integration) and the sound pressure channel according to expression (6).

Напротив, для выделения неоднородной части гидроакустической волны, в частности вытекающей, например, волны Стоунли, вычисляют мнимую часть взаимного спектра по разности откликов компонент ПГД и компонент акселерометра (после интегрирования) и канала звукового давления, согласно выражению (7).On the contrary, to isolate the inhomogeneous part of the hydroacoustic wave, in particular the flowing out, for example, Stoneley wave, the imaginary part of the mutual spectrum is calculated from the difference in responses of the PGD components and the components of the accelerometer (after integration) and the sound pressure channel, according to expression (7).

Вычисления взаимного спектра осуществляются в блоке обработки с помощью программы, например, SpectraPLUS (Pioneer Hill Software), находящемся либо рядом с устройством, либо на поверхности воды в специальном буе, либо на берегу. Вещественная и мнимая часть взаимного спектра вычисляется, например, с помощью переноса рассчитанных в пакете SpectraPLUS (Pioneer Hill Software) модуля и фазы в виде столбцов в любую электронную таблицу, например, Excel (Microsoft), путем умножения модуля на косинус или синус фазы соответственно.Cross-spectrum calculations are performed in a processing unit using a program, for example, SpectraPLUS (Pioneer Hill Software), located either next to the device, or on the water surface in a special buoy, or on the shore. The real and imaginary parts of the mutual spectrum are calculated, for example, by transferring the module and phase calculated in the SpectraPLUS (Pioneer Hill Software) package in the form of columns to any spreadsheet, for example, Excel (Microsoft), by multiplying the modulus by the cosine or sine of the phase, respectively.

Таким образом, предложенное решение позволяет решить поставленную проблему разделения гидроакустической и сейсмоакустической составляющих из совокупной волны, воздействующей на приемное устройство, и значительно расширить функциональные возможности заявляемого решения, позволив затем решать множество задач, стоящих перед гидроакустикой и геофизикой, в том числе может быть использовано для обнаружения удаленного источника звука или получения информации о структуре и характеристиках морского днаThus, the proposed solution allows us to solve the posed problem of separating the hydroacoustic and seismoacoustic components from the aggregate wave acting on the receiving device, and significantly expand the functionality of the proposed solution, allowing then to solve many problems facing hydroacoustics and geophysics, including can be used for detecting a distant sound source or obtaining information about the structure and characteristics of the seabed

Claims (5)

1. Способ приема сейсмоакустической и гидроакустической составляющих волны у сложенного осадочными породами дна водоема, включающий направленный прием совокупной волны одновременно соосно расположенными в едином жестком корпусе 3-компонентным приемником градиента давления силового типа, снабженным приемником звукового давления, и 3-компонентным приемником вибраций, совместную обработку отклика соосных компонент 3-компонентных приемников и приемника вибраций с последующим выделением гидроакустической составляющей путем вычисления разности откликов компонент приемника градиента давления и откликов соосных компонент приемника вибраций с предварительной настройкой коэффициента усиления каналов приемника вибраций на максимум подавления модельных вибраций, при этом приемник вибраций размещен в центре масс корпуса, расположенного в грунте, а приемник градиента давления с приемником звукового давления - на поверхности корпуса, находящейся в воде.1. A method of receiving seismoacoustic and hydroacoustic components of a wave at the bottom of a reservoir folded by sedimentary rocks, including a directional reception of the aggregate wave simultaneously coaxially located in a single rigid body with a 3-component force-type pressure gradient receiver equipped with a sound pressure receiver and a 3-component vibration receiver, joint processing the response of the coaxial components of the 3-component receivers and the vibration receiver, followed by the extraction of the hydroacoustic component by calculating the difference in the responses of the pressure gradient receiver components and the responses of the coaxial components of the vibration receiver with preliminary adjustment of the amplification factor of the vibration receiver channels to the maximum suppression of model vibrations, while the vibration receiver is placed in the center of mass of the housing located in the ground, and the pressure gradient receiver with the sound pressure receiver is on the surface of the housing in the water. 2. Способ по п. 1, отличающийся тем, что осуществляют операцию усиления разности откликов компонент приемника градиента давления и соосных компонент приемника вибраций с предварительно определенным коэффициентом усиления, а результат суммируют и/или вычитают с откликом приемника звукового давления, формируя однонаправленные диаграммы направленности типа кардиоиды в направлении каждой из компонент.2. The method according to claim 1, characterized in that the operation of amplifying the difference between the responses of the components of the pressure gradient receiver and the coaxial components of the vibration receiver with a predetermined gain is performed, and the result is summed and / or subtracted with the response of the sound pressure receiver, forming unidirectional radiation patterns of the type cardioids in the direction of each of the components. 3. Способ по п. 1, отличающийся тем, что осуществляют операцию усиления разности откликов соосных компонент приемников градиента давления и вибраций с предварительно определенным коэффициентом усиления и вычисляют взаимный спектр полученной разности откликов по каждой компоненте с откликом приемника звукового давления, из которого выделяют вещественную часть, характеризующую плоскую гидроакустическую волну от удаленного источника, и мнимую часть, характеризующую неоднородную гидроакустическую волну.3. The method according to claim 1, characterized in that the operation of amplifying the difference in responses of the coaxial components of the pressure gradient and vibration receivers with a predetermined gain is carried out and the mutual spectrum of the obtained difference in responses is calculated for each component with the response of the sound pressure receiver, from which the real part is extracted , characterizing a plane hydroacoustic wave from a distant source, and an imaginary part characterizing an inhomogeneous hydroacoustic wave. 4. Способ по п. 3, отличающийся тем, что вычисляют отношения амплитуд компонент вещественной части взаимного спектра и определяют пеленг на источник гидроакустической волны.4. The method according to claim 3, characterized in that the ratio of the amplitudes of the components of the real part of the mutual spectrum is calculated and the bearing to the source of the hydroacoustic wave is determined. 5. Устройство приема сейсмоакустической и гидроакустической волн у дна водоема, сложенного осадочными породами, состоящее из единого жесткого корпуса со средней плотностью, близкой к плотности грунта, в центре масс которого размещен 3-компонентный приемник вибраций, а на поверхности соосно компонентам приемника вибраций установлен 3-компонентный приемник градиента давления силового типа, дополнительно снабженный приемником звукового давления, установленным с обеспечением единого фазового центра с компонентами приемника градиента давления, при этом электрические выходы соосных компонент приемника градиента давления и приемника вибраций соединены через усилители с блоками вычитания, а выходы блоков вычитания через усилители совместно с выходами приемника звукового давления присоединены к блоку обработки данных.5. A device for receiving seismoacoustic and hydroacoustic waves at the bottom of a reservoir composed of sedimentary rocks, consisting of a single rigid body with an average density close to the density of the soil, in the center of mass of which there is a 3-component vibration receiver, and on the surface, coaxially with the components of the vibration receiver, 3 - a power-type component pressure gradient receiver, additionally equipped with a sound pressure receiver installed to ensure a single phase center with the components of the pressure gradient receiver, while the electrical outputs of the coaxial components of the pressure gradient receiver and vibration receiver are connected through amplifiers with subtraction units, and the outputs of the subtraction units through the amplifiers together with the outputs of the sound pressure receiver are connected to the data processing unit.
RU2020123339A 2020-07-07 2020-07-07 Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation RU2740334C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2020123339A RU2740334C1 (en) 2020-07-07 2020-07-07 Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2020123339A RU2740334C1 (en) 2020-07-07 2020-07-07 Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation

Publications (1)

Publication Number Publication Date
RU2740334C1 true RU2740334C1 (en) 2021-01-13

Family

ID=74183876

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2020123339A RU2740334C1 (en) 2020-07-07 2020-07-07 Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation

Country Status (1)

Country Link
RU (1) RU2740334C1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4134097A (en) * 1977-06-13 1979-01-09 Shell Oil Company Combination geophone-hydrophone
US4599713A (en) * 1983-03-30 1986-07-08 Compagnie General De Geophysique Seismic detector
EP0762146A1 (en) * 1995-09-12 1997-03-12 PRAKLA-SEISMOS GmbH Seismic sensor
RU2677097C1 (en) * 2018-04-05 2019-01-15 Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) Three-component vector-scalar receiver

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4134097A (en) * 1977-06-13 1979-01-09 Shell Oil Company Combination geophone-hydrophone
US4599713A (en) * 1983-03-30 1986-07-08 Compagnie General De Geophysique Seismic detector
EP0762146A1 (en) * 1995-09-12 1997-03-12 PRAKLA-SEISMOS GmbH Seismic sensor
RU2677097C1 (en) * 2018-04-05 2019-01-15 Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) Three-component vector-scalar receiver

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Shi J., Dosso S.E., Sun D., Qingyu L., "Geoacoustic inversion of the acoustic-pressure vertical phase gradient from a single vector sensor", J. Acoust. Soc. Am., 2019, 146 (5), P. 3159-3173. *

Similar Documents

Publication Publication Date Title
Schock et al. Buried object scanning sonar
NO303033B1 (en) System for attenuation of water column reverberations
Najeem et al. Open lake experiment for direction of arrival estimation using acoustic vector sensor array
CN109239712B (en) Noise detection method based on underwater sound field and sound energy flow
JPS625301B2 (en)
NO328966B1 (en) Procedure for collecting and processing data from seismic survey using subsea cables and streamers
RU2569414C2 (en) Noise suppression system with dual transducer for underwater hydrophone streamer
US5930201A (en) Acoustic vector sensing sonar system
Châtillon et al. SAMI: A low-frequency prototype for mapping and imaging of the seabed by means of synthetic aperture
Li et al. Passive acoustic localisation of undersea gas seeps using beamforming
US9453900B2 (en) Method and apparatus for three dimensional wavenumber-frequency analysis
US6408978B1 (en) Non-acoustic self noise canceller for sensor arrays
CN114563075B (en) A Separation Method for Multi-path Arrival Structure of Deep Sea Sound Field Based on Single Vector Hydrophone
CN109632258B (en) Marine internal wave acoustic detection method based on transmitting-receiving separation of vector sensor
Zou et al. Vector hydrophone array development and its associated DOA estimation algorithms
RU2740334C1 (en) Method of receiving seismic-acoustic and hydroacoustic waves at the bottom of a water reservoir and device for its implementation
CN114488010A (en) Virtual receiving array amplification method for linear array underwater sound source spacing
Abraham Low‐cost dipole hydrophone for use in towed arrays
Liu et al. Research on DOA estimation method of single MEMS vector hydrophone based on pulse signal
JP5448153B2 (en) Artifact detection system, and artifact detection method used in the system
CN115220026A (en) A method for cross-ice localization in matching field based on horizontal array of ice surface
Huerta-Lopez et al. In situ evaluation of shear-wave velocities in seafloor sediments with a broadband ocean-bottom seismograph
JP2558637B2 (en) Sound source position measurement method
Jesus A sensitivity study for full-field inversion of geo-acoustic data with a towed array in shallow water
Sanders et al. Direction of Arrival (DOA) Estimation using Vector Acoustic Data