Method of using a parabolic equation model for range-dependent seismo-acoustic problems

Patent metadata

Field	Value
Publication number	US-9945971-B2
Application number	US-201615015425-A
Country	US
Kind code	B2
Filing date	Feb 4, 2016
Priority date	Feb 10, 2015
Publication date	Apr 17, 2018
Grant date	Apr 17, 2018

Abstract

A method of modeling an aquatic environment or locating an acoustic source in the aquatic environment. A range-dependent medium is approximated in terms of a series of range-independent regions and obtaining single-scattering solutions across the vertical interfaces between regions. One or more acoustic waves are propagated from a known acoustic source through the range-dependent medium to one or more known seismoacoustic receivers to model iteratively the various solid and liquid layers of the range-dependent medium. Alternatively, one or more acoustic waves are reverse-propagated from one or more known seismoacoustic receivers through the range-dependent medium to determine whether an acoustic source is present within a user-defined range.

First claim

What is claimed as new and desired to be protected by letters Patent of the United States is:

1. A method of locating an acoustic source within a user-defined maximum range, the method comprising:
receiving an acoustic signal from the acoustic source through an aquatic environment using at least one seismoacoustic sensor, the acoustic signal comprising an acoustic pressure and an acoustic signal frequency, the acoustic signal comprising compressional waves and shear waves, the environment comprising a plurality of layers, the plurality of layers comprising at least one of solid layers and liquid layers, the solid layers propagating the compressional waves and the shear waves, the liquid layers propagating the compressional waves;
determining from the received acoustic signal a plurality of initial conditions, the plurality of initial conditions being based on the acoustic pressure and the acoustic signal frequency;
digitizing the aquatic environment into a plurality of range-independent regions and a plurality of vertical interfaces, each vertical interface of the plurality of vertical interfaces being located between adjacent range-independent regions of the plurality of range-independent regions;
determining an initial transmitted field for a range-independent region of the plurality of range-independent regions based on the plurality of initial conditions and determining a subsequent incident field;
determining a subsequent transmitted field for an adjacent vertical interlace of the plurality of vertical interfaces based on the subsequent incident field and determining a next incident field on an adjacent range-independent region of the plurality of range-independent regions;
determining the subsequent transmitted field for another range-independent region of the plurality of range-independent regions based on the next incident field and determining the subsequent incident field on another vertical interface of the at least one vertical interface;
repeating for the user-defined maximum range determining the subsequent transmitted field for another range-independent region of the plurality of range-independent regions based on the next incident field and determining the subsequent incident field on another vertical interface of the at least one vertical interface, thereby approximating a propagation of the received acoustic signal propagating from the acoustic source via the environment to the at least one seismoacoustic sensor;
determining an acoustic source location based on the approximation.

2. The method according to claim 1 , wherein said determining the subsequent transmitted field for an adjacent vertical interface of the plurality of vertical interfaces based on the subsequent incident field and determining a next incident field on an adjacent range-independent region of the plurality of range-independent regions comprises using single scattering with a parabolic wave equation.

3. The method according to claim 2 , wherein said determining the subsequent transmitted field for an adjacent vertical interlace of the plurality of vertical interfaces based on the subsequent incident field and determining a next incident field on an adjacent range-independent region of the plurality of range-independent regions comprises conserving horizontal displacement and tangential stress, and vertical displacement and normal stress.

4. The method according to claim 3 , wherein said transmitted field comprises an average of the conserved horizontal displacement and the tangential stress, and the conserved vertical displacement and the normal stress.

5. The method according to claim 1 , wherein the at least one seismoacoustic sensor comprises at least one of a one-dimensional vertical hydrophone array, a two-dimensional vertical hydrophone array, a three-dimensional vertical hydrophone array, a one-dimensional horizontal hydrophone array, a two-dimensional horizontal hydrophone array, a three-dimensional horizontal hydrophone array, a vertical water-body bottom-mounted hydrophone array, a horizontal water-body bottom-mounted hydrophone array, and a seismometer array.

6. The method according to claim 1 , wherein the acoustic signal further comprises a solid layer motion, the plurality of initial conditions being further based on the solid layer motion.

7. The method according to claim 1 , wherein the at least one seismoacoustic sensor is located one of on land, on a body of water, in the body of water, and under the body of water.

8. The method according to claim 1 , wherein the acoustic source is located one of on land, in a body of water, and under the body of water.

9. A method of locating at least one object of interest in an aquatic environment within a user-defined maximum range, the method comprising:
receiving an actual acoustic signal from an acoustic source through the aquatic embroilment using at least one seismoacoustic sensor, the actual acoustic signal comprising an acoustic pressure, an acoustic signal frequency, and an acoustic source location, the actual acoustic signal comprising compressional waves and shear waves, the aquatic environment comprising a plurality of actual layers, the plurality of actual layers comprising at least one of solid actual layers and liquid actual layers, the solid actual layers propagating the compressional waves and the shear waves, the liquid actual layers propagating the compressional waves;
determining from the received actual acoustic signal a plurality of initial conditions, the plurality of initial conditions being based on the acoustic pressure, the acoustic signal frequency, and the acoustic source location;
digitizing an aquatic environment, model into a plurality of range-independent regions and a plurality of vertical interfaces, each vertical interface of the plurality of vertical interfaces being located between adjacent range-independent regions of the plurality of range-independent regions, the initial aquatic environment model comprising a plurality of model layers, the plurality of model layers comprising at least one of solid model layers and liquid, model layers, the plurality of model layers comprising a plurality of respective thickness geometries and a plurality of respective material densities;
determining an initial transmitted field for a range-independent region of the plurality of range-independent regions based on the plurality of initial conditions and determining a subsequent incident field;
determining a subsequent transmitted field for an adjacent vertical interface of the plurality of vertical interlaces based on the subsequent incident field and determining a next incident field on an adjacent range-independent region of the plurality of range-independent regions;
determining the subsequent transmitted field for another range-independent region of the plurality of range-independent regions based on the next incident field and determining the subsequent incident field on another vertical interface of the at least one vertical interface;
repeating for the user-defined maximum range determining the subsequent transmitted field for another range-independent region of the plurality of range-independent regions based on the next incident field and determining the subsequent incident field on another vertical interface of the at least one vertical interface, thereby approximating a propagation of the acoustic signal propagating from the acoustic source via the aquatic environment to the at least one seismoacoustic sensor;
determining a model acoustic signal based on the approximation;
determining whether the model acoustic signal is converging toward the received actual acoustic signal;
adjusting the aquatic environment model, if the model acoustic signal is not converging

Assignees

• Us Navy

Inventors

• Collins Michael D

Classifications

• G01V1/38 Primary

  specially adapted for water-covered areas (G01V1/28 takes precedence) · CPC title

• G01V2210/123

  Passive source, e.g. microseismics · CPC title

• G01S5/18

  using ultrasonic, sonic or infrasonic waves · CPC title

• G01V1/186

  Hydrophones · CPC title

• G01V2210/1232

  Earthquakes · CPC title