JP2020063958A - Position estimation device and method - Google Patents

Abstract

To estimate a position of a short-time, unknown sound source in the ocean environment having depth dependency of sound velocity distribution and distance dependency of sea bottom topography.SOLUTION: A position estimation method obtains each sound ray parameter of direct waves, sea surface reflection waves, and sea bottom reflection waves of natural sound rays of pulse waves passing through two points, namely tentative coordinates of a sound source and coordinates of a receiver for receiving pulse waves generated by the sound source, obtains the amount of perturbation of a difference between a propagation time difference actual measurement value between direct waves and sea surface reflection waves and between direct waves and sea bottom reflection waves and a propagation time difference calculation value between direct waves and sea surface reflection waves and between direct waves and sea bottom reflection waves obtained from sound ray parameters of each wave by the perturbation method, and determines estimated coordinates obtained this time as coordinates of the sound source when an error value between estimated coordinates of the sound source estimated from the amount of perturbation and estimated coordinates obtained last time goes below a determination threshold.SELECTED DRAWING: Figure 1

Description

  The present invention relates to a position estimation device and method for passively receiving a sound emitted from a sound source in the sea and estimating the position of the sound source.

  Generally, sound waves are used to detect an underwater object because electromagnetic waves are greatly attenuated in water. The active method and the passive method are known as methods for detecting an underwater object by sound waves.

  The active method uses echo localization, detects the echo sound reflected from an object by transmitting sound from a sound source, and calculates the distance and direction to the object. You will be exposed. On the other hand, the passive method detects the radiated sound of an underwater object and does not make a sound by itself, so that the underwater object can be secretly detected.

  However, in the passive method, it is easy to find the orientation of the object, but complicated calculations such as target motion analysis and matched field processing are required to find the distance to the object. In the target motion analysis, the distance to the object can be estimated by continuous observation of radiated sound, but the distance to the object cannot be obtained by the transient transient sound. In the case of matched field processing, it is necessary to know the marine environment in advance, but it is difficult to create a database of the complicated characteristics of the marine environment that change in space and time.

  As a method of obtaining the distance to an object by a transient sound, a method of utilizing the propagation time difference of sound waves has been studied. As such a method, there is a method called a brute force method disclosed in Patent Documents 1 and 2 below.

  In the methods of Patent Documents 1 and 2 below, a target position is obtained by using a propagation time difference between a sound wave (direct wave) that directly propagates in water from a sound source and arrives and a sound wave (reflected wave) that reflects and arrives at a boundary surface. The entire ocean space is gridded in Cartesian coordinates and polar coordinates within the assumed range of distances and depths, and the sound wave propagation time difference from the temporary target on the grid points is exhaustively calculated to propagate the sound wave from the actual target. The distance or depth that minimizes the difference from the time difference (actual measurement value) is estimated as the target position. Therefore, the brute force method is a reliable method because it requires an iterative calculation for the number of grid points determined by the estimation range (depth and distance) of the target position and the calculation step size. There is a problem that it takes time.

  Therefore, as another method for solving the above problem, for example, there is a hyperbolic method disclosed in Patent Document 3 below. This method consists of a direct wave that reaches the receiver directly from the sound source, a first reflected wave that reaches the receiver after being reflected once by the sea surface, and a first reflected wave that reaches the receiver after reflecting once at the sea bottom. 2 The reflected wave is received by the wave receiver installed between the sea surface and the sea bottom, and the formula of the first hyperbola whose focus is the sound source and the mirror image of the sound source with respect to the sea surface, and the sound source and the sound source with respect to the sea bottom The common focus coordinates of the first hyperbola and the second hyperbola, which are the coordinates of the sound source, are obtained from the formula of the second hyperbola whose focal point is the mirror image.

JP, 2006-194627, A JP, 2011-149864, A Japanese Patent Laid-Open No. 2015-152376

  However, since the method disclosed in Patent Document 3 assumes a uniform environment, when applied to a marine environment in which the sound velocity distribution has depth dependency or the seabed topography has distance dependency. , There was a problem that the sound source position could not be calculated accurately. Therefore, it is required to establish a new method for estimating the sound source position in a short time in an environment where the sound velocity distribution depends on the depth and the seabed topography depends on the distance.

  Therefore, the present invention has been made in view of the above problems, and a position estimation device and method capable of estimating the position of a sound source in a short time regardless of the depth dependency of the sound velocity distribution and the distance dependency of the seabed topography. The purpose is to provide.

In order to achieve the above-mentioned object, a first aspect of the present invention is a position estimation device that estimates the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed is flat,
Sound ray parameters of a direct wave, a sea surface reflected wave, and a seabed reflected wave that are eigen sound rays of the pulse wave that passes through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source. And the process of
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A process for obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset determination threshold value, a process of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
A position estimation device comprising an arithmetic processing unit for performing the above.

A second aspect of the present invention is a position estimation device that estimates the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed undulates,
A process of obtaining sound ray parameters of a direct wave and a sea surface reflected wave that are eigen sound rays of the pulse wave passing through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source;
Temporary coordinates of the sound source, a process of obtaining sound ray parameters of the seabed reflected wave from the virtual horizontal distance from the receiver to the seabed reflection point,
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A process for obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset determination threshold value, a process of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
A position estimation device comprising an arithmetic processing unit for performing the above.

  A third aspect of the present invention is the position estimation apparatus according to the first or second aspect of the present invention, wherein the sound ray parameter starts from the coordinates of the wave receiver and uses temporary coordinates of the sound source. A position estimation device characterized by being obtained by a numerical solution method of a nonlinear equation.

  A fourth aspect of the present invention is the position estimation apparatus according to any one of the first to third aspects of the present invention, wherein when the sound ray parameter of the direct wave does not have a real number solution, the sound of the direct wave is A position estimation device, characterized in that the line is assumed to pass through a turning point, the sound ray parameter of the direct wave is obtained based on the turning point depth, and the coordinates of the sound source are estimated.

  A fifth aspect of the present invention is, in the position estimation device according to any one of the second to fourth aspects of the present invention, known seabed topographical coordinates and a virtual horizontal distance from the receiver to the seabed reflection point. The position estimating device is characterized in that the seabed reflection point coordinates of the seabed reflected wave and the sound ray parameter of the seabed reflected wave are obtained based on the above.

A sixth aspect of the present invention is a position estimating device for estimating the coordinates of a sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver for receiving the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the speed of sound, the depth coordinate is known. It is a position estimation device characterized by estimating coordinates of the sound source by determining a distance to the sound source.

A seventh aspect of the present invention is a position estimation device that estimates the coordinates of a sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver that receives the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the distance coordinates are known. It is a position estimation device characterized by estimating the coordinates of the sound source by determining the depth of the sound source.

An eighth aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed is flat,
Sound ray parameters of a direct wave, a sea surface reflected wave, and a seabed reflected wave that are eigen sound rays of the pulse wave that passes through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source. And the step of finding
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
It is a position estimation method characterized by including.

A ninth aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed undulates,
A step of obtaining sound ray parameters of a direct wave and a sea surface reflected wave which are characteristic sound rays of the pulse wave passing through two points of temporary coordinates of the sound source and coordinates of a receiver for receiving the pulse wave emitted by the sound source;
Temporary coordinates of the sound source, and a step of obtaining a sound ray parameter of the seabed reflected wave from the virtual horizontal distance from the receiver to the seabed reflection point,
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
It is a position estimation method characterized by including.

  A tenth aspect of the present invention is the position estimation method according to the eighth or ninth aspect of the present invention, wherein the sound ray parameter starts from the coordinates of the receiver and uses temporary coordinates of the sound source. It is a position estimation method characterized by being obtained by a numerical solution of a nonlinear equation.

An eleventh aspect of the present invention is the position estimation method according to any one of the eighth to tenth aspects of the present invention, wherein when the sound ray parameter of the direct wave does not have a real number solution, the sound of the direct wave is It is a position estimation method characterized in that the sound ray parameter of the direct wave is obtained based on the turning point depth on the assumption that the line passes through the turning point.

A twelfth aspect of the present invention is, in the position estimation method according to any one of the eighth to tenth aspects of the present invention, a known seabed topographic coordinate and a virtual horizontal distance from the wave receiver to the seabed reflection point. The position estimating method is characterized in that the sea bottom reflection point coordinates of the sea bottom reflected wave and the sound ray parameter of the sea bottom reflected wave are obtained based on the above.

A thirteenth aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver for receiving the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the speed of sound, the depth coordinate is known. A position estimation method characterized in that the coordinates of the sound source are estimated by determining the distance to the sound source.

A fourteenth aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver that receives the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the distance coordinates are known. A position estimation method characterized in that the coordinates of the sound source are estimated by determining the depth to the sound source.

  According to the present invention, in an environment in which the sound velocity distribution has depth dependency and the seabed topography has distance dependency, an echo that grasps a behavior while observing a pulse wave emitted by a ship or a large marine organism underwater vehicle and estimating a position If it is applied to a location or the like, it is possible to estimate the sound source position with high accuracy in a short time without requiring a huge amount of arithmetic processing like the conventional brute force method.

It is a schematic diagram showing an example of an eigen sound line used by a position estimating method in a position estimating device of the present invention. It is a schematic diagram which shows the relationship between the depth distribution of sound velocity, and the sound ray of a propagation sound. It is a schematic functional block diagram which shows the electric constitution of the position estimation apparatus of this invention. It is a flow chart which shows a flow of position presumption processing of a position presumption device of a 1st embodiment. It is a figure showing an example of a turning point. It is a schematic diagram which showed the incident elevation angle of a sound wave to the undulating sea bottom, and the relationship of the reflection depression angle. FIG. 3 is a schematic diagram showing an example of sound rays of a pulse wave emitted from a sound source S and shadow zones of the sound rays.

  Hereinafter, embodiments for carrying out the present invention will be described in detail with reference to the accompanying drawings. Further, the present invention is not limited to the embodiments, and other forms, examples, operational techniques and the like that can be implemented by those skilled in the art based on the embodiments are all included in the scope of the present invention. I shall.

[Outline of the present invention]
FIG. 1 is a schematic diagram showing an example of the characteristic sound ray used in the position estimation device 1 of the present invention. In FIG. 1, the height direction (depth direction) is represented by the z-axis, and the horizontal distance is represented by the r-axis. The following is an example of obtaining the coordinates (r, z) of the sound source S when the wave receiver R is at the coordinates (r ₀ , z ₀ ) as shown in FIG.

  Normally, when a sound wave (pulse wave) is radiated from the sound source S, the direct wave dw traveling on a straight line connecting the wave receiver R and the sound source S reaches the receiver R, and then the sea surface. The sea surface reflected wave sw reflected once by W1 and the sea bottom reflected wave bw reflected once by the sea bottom W2 reach the receiver R, and then the multiple reflected waves at the sea water surface W1 and the sea bottom W2 are received. Reach vessel R.

The present invention, in an marine environment in which the sound velocity distribution has depth dependency and the seabed topography have distance dependency, determines whether or not there is depth dependence of the sound velocity distribution from the underwater sound source S and distance dependence of the seabed topography by gridding only the depth direction. However, in order to estimate the position of the sound source S in a short time, the direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw, which are pulse waves from the unknown sound source S in the sea, are received by the receiver R, and The position (coordinates) of the sound source S is estimated using the sound ray parameter of the wave.

[First Embodiment]
Next, a first embodiment of the present invention will be described.
In the first embodiment, the speed of sound changes in the depth direction, and the sea surface W1 and the sea bottom W2 are planes having no unevenness and are parallel to each other as an ocean environment. The space between the sea surface and the sea bottom is filled with a medium (sea water) that transmits sound waves.

  FIG. 2 is a diagram showing the relationship between the depth distribution of sound velocity and the sound ray of propagating sound. When the sea surface position is the origin O, the height direction (depth direction) is represented by z, and the horizontal direction is left. FIG. 2A shows the sound velocity c (z), and FIG. 2B shows the distance r.

For sound propagating along one sound line from a depth z ₀ at a vertical elevation angle θ ₀ to a depth z that is a depression angle θ (z), the following formula 1 is established according to Snell's law. , P have values that do not change on one sound ray path, so this value p is used as a sound ray parameter, and the display using this p is used as a parameter display.

The horizontal component of the sound ray path length at this time is the distance r, which can be expressed by the following equation 2 in the parameter display. The functional on the right side represented by the parametric representation is represented as R (z ₀ , z, p).

Here, in estimating the coordinates of the sound source S, the water depth h is fixed, the coordinates (r ₀ , z ₀ ) of the receiver R (receiving point) are known, and the coordinates (r, z) of the sound source S are unknown. Assuming that a sound ray passing through the sound source S and the receiving point is calculated, the sound ray path from the wave receiver R to the sound source S and the sound ray path from the sound source S to the wave receiver R are calculated from the reciprocity theorem. Will be the same. Therefore, in the first embodiment, a process of estimating the position coordinates of the sound source S is performed by using the wave receiver R, which is the wave receiving point, as a starting point.

Next, the system configuration of the position estimation device 1 according to the first embodiment will be described.
As shown in FIG. 3, the position estimation device 1 includes a wave receiver R that receives a pulse wave (direct wave dw, a sea surface reflected wave sw, and a seabed reflected wave bw) radiated from a sound source S, and a wave receiver R. Based on the data input unit 10 for inputting various information necessary for the position estimation processing by the received pulse wave, the information obtained from each pulse wave received by the receiver R and the various information input by the data input unit 10 The calculation processing unit 20 is configured to calculate the estimated coordinates of the sound source S.

<Data input section>
The data input unit 10 inputs various data such as an initial value and an actual measurement value necessary for estimating the position of the sound source S. This data includes, for example, the sound velocity distribution c (z) indicating the depth dependence of the sound velocity, the coordinates (r ₀ , z ₀ ) of the receiver R, the direct wave dw received by the receiver R, and the sea surface reflected wave sw. Of the measured propagation time difference (τ ₁₂ ^(ob) ) and the measured value of the propagation time difference (τ ₁₃ ^(ob) ) between the direct wave dw and the ocean bottom reflected wave bw. The ocean environment data including the ocean topographical coordinates in which the latitude / longitude of each mesh position and the water depth data of each mesh position are associated with each other by being divided into meshes, and temporary coordinates of the sound source S are included.
The temporary coordinates are the coordinates set when the estimated coordinates of the sound source S are calculated, and when the estimation process is the first time, the initial value coordinates (r ⁽⁰⁾ and z ⁽⁰⁾ ) are set.

<Arithmetic processing unit>
The arithmetic processing unit 20 includes a sound ray parameter calculation unit 21, a propagation time difference calculation unit 22, a perturbation calculation unit 23, a coordinate calculation unit 24, and a determination output unit 25. The direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw that reach the receiver R from the sound source S identify the type of pulse wave based on the arrival order of each pulse wave and the reception angle.

The sound ray parameter calculation unit 21 calculates sound ray parameters p (p ₁ , p ₂ , p ₃ ) of the pulse waves (direct wave dw, sea surface reflected wave sw, sea bottom reflected wave bw) emitted from the sound source S, respectively. . Here, "p _1" is sound ray parameters of the direct wave dw, "p _2" is the sound ray parameter of sea reflected sw, "p _3" represents the sound ray parameters of the seabed reflected wave bw.

The sound ray parameter calculation unit 21 uses initial values coordinates (r ⁽⁰⁾ , z ^{(0 )} of the sound source S as temporary coordinates when the direct wave dw with the receiver R as a starting point is within the section [z ₀ , z]. ⁾ ) And the above equation 2, the r coordinate of the direct wave dw is given by the following equation 3.

Although this functional is non-linear, it is a monotonically increasing function of the sound ray parameter p ₁ of the direct wave dw. Therefore, by giving the initial value coordinates of the sound source S as the temporary coordinates to the values of r and z, the above equation 3 becomes a function of p ₁ , and the sound ray parameter of the direct wave dw is obtained by the numerical solution of the nonlinear equation represented by the Newton method. p ₁ is obtained.

  Further, similarly to the direct wave dw, the r coordinate of the sea surface reflected wave sw from the sound source S is given by the following formula 4, and the r coordinate of the seabed reflected wave bw is given by the following formula 5.

Although these functionals are also non-linear, they are monotonically increasing functions of the sound ray parameters p ₂ and p ₃ , and by giving the initial value coordinates of the sound source S as temporary coordinates to the values of r and z, the above equations 4 and 5 are given. Becomes a function of p ₂ and p ₃ , and sound ray parameters p ₂ and p ₃ of the sea surface reflected wave sw and the seabed reflected wave bw are obtained by a numerical solution of a nonlinear equation represented by the Newton method.

Of the sound ray parameters calculated by the sound ray parameter calculation unit 21, immediately after the position estimation processing is started, the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) preset as temporary coordinates and the above number are calculated. 3 to 5 are appropriately used to calculate the initial values (p ₁ ⁽⁰⁾ , p ₂ ⁽⁰⁾ , p ₃ ⁽⁰⁾ ) of the sound ray parameters of each wave, and the coordinate calculation unit 24 calculates the estimated coordinates of the sound source S. Once calculated, the calculated estimated coordinates of the sound source S are updated and used as temporary coordinates. That is, when the estimation process is the first time (first time), the preset initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) are used, and from the second time onward, the sound source S calculated by the coordinate calculation unit 24 is calculated. The estimated coordinates (r ^(n-1) , z ^(n-1) ) are used. Here, n is the number of estimation processes. When the sound ray parameter p ₂ of the sea surface reflected wave sw does not have a real number solution, it is assumed that the sea surface reflected wave sw does not exist, and the coordinate estimation processing of the sound source S, which is the processing of the present device, is terminated.

Transit time calculation unit 22, a direct wave dw, with calculating the propagation time of the sea reflected sw and seabed reflected wave bw, the propagation time difference calculated value tau ₁₂ between the direct wave dw and sea reflected sw, and direct wave dw The calculated propagation time difference τ ₁₃ of the sea bottom reflected wave bw is calculated.

The propagation time from the depth z ₀ to z in the propagation path s along the sound ray is represented by t, and this is expressed as the following equation 6, and the functional on the right side is T (z ₀ , z, p). Represent

Then, the propagation time T ₁ of the direct wave dw, the propagation time T ₂ of the sea reflected sw, the propagation time T ₃ of the submarine reflected wave bw is respectively obtained as the following Equation 7-9, the direct wave dw as the difference between these propagation times And a propagation time difference τ ₁₂ (= T ₂ −T ₁ ) of the sea surface reflected wave sw, and a propagation time difference τ ₁₃ (= T ₃ −T ₁ ) of the direct wave dw and the seabed reflected wave bw are obtained as propagation time difference calculated values.

When the estimation process is performed for the first time, the initial value (p ₁ ⁽⁰⁾⁾ of the sound ray parameter of each wave obtained from the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) of the sound source S which is the temporary coordinate. , P ₂ ⁽⁰⁾ , p ₃ ⁽⁰⁾ ), the initial values (τ ₁₂ ⁽⁰⁾ , τ ₁₃ ⁽⁰⁾ ) of the calculated values of the propagation time difference of each wave are calculated. In addition, after the second time, the calculated estimated propagation coordinates (r ^(n-1) , z ^(n-1) ) are used as temporary coordinates, and the calculated propagation time difference of each wave (τ ₁₂ ^(n-1) , τ ₁₃ ^{( n-1)} ) is calculated.

The perturbation calculator 23 calculates the perturbation amount of the estimated coordinates of the sound source S by the perturbation method. Since there are 5 equations for the above 5 variables (r, z, p ₁ , p ₂ , p ₃ ), it is possible to solve the equations in principle, but the series of equations is non-linear. Therefore, it is difficult to obtain a solution analytically. In addition, there are probabilistic solution methods such as the Simulated annealing method and the Genetic Algorithm method in order to find the solution of a complex nonlinear equation system, but these methods generally require calculation time and have difficulty in determining the convergence of the solution. .
Therefore, in the present invention, the solution of the nonlinear equation is obtained by the perturbation method.

Let the variables r, z, p ₁ , p ₂ , p ₃ be slightly changed and be δr, δz, δp ₁ , δp ₂ , δp ₃ , respectively, and consider the perturbation up to the first term of the Taylor expansion. It is applied to the above expressions 3 to 5 and expressions 7 to 9. There are two independent variables among the variation quantities δr, δz, δp ₁ , δp ₂ , and δp ₃ , and δp ₁ , δp ₂ , and δp ₃ can be represented by δr and δz. Is obtained.

Here, “δr” in the equation 10 is a column vector whose components are the variation amounts δr and δz as shown in the following equation 11, and G is a matrix of 2 rows and 2 columns and z, p _{i of the} functionals R and T. It is composed of partial differential coefficients with respect to (i = 1, 2, 3). In addition, “δτ” in Expression 10 is the difference between the actual measurement value of the propagation time difference, which is the observed value, and the calculated value of the propagation time difference, which is the calculated value calculated by the propagation time difference calculation unit 22. For example, when the estimation process is the first time, the measured values of the propagation time difference (τ ₁₂ ^(ob) , τ ₁₃ ^(ob) ) that are observed values and the initial value (τ _It is the difference between ₁₂ ⁽⁰⁾ and τ ₁₃ ⁽⁰⁾ ).

  Further, since the column vector is first-order independent, an inverse matrix of G exists, and the following Expression 13 is obtained.

That is, the perturbation calculation unit 23 calculates the propagation time difference calculated by the propagation time difference calculation unit 22 and the measured propagation time difference between the direct wave dw and the sea surface reflected wave sw input by the data input unit 10 (τ ₁₂ ^(ob) ) And the difference between the measured propagation time difference (τ ₁₃ ^(ob) ) between the direct wave dw and the reflected sea wave bw (τ ₁₃ ^(ob) ), and calculates the matrix of Equation 13 from the partial differential coefficient of the functional R, T to obtain the inverse matrix. And the perturbation amount of the estimated coordinates of the sound source S is calculated from the above equation 13.

The coordinate calculation unit 24 obtains the estimated coordinates (r ⁽ⁿ⁾ , z ⁽ⁿ⁾ ) of the sound source S using the perturbation amount calculated by the perturbation calculation unit 23 and the following Expression 14. The coordinates obtained by the coordinate calculation unit 24 are estimated coordinates estimated as the current position of the sound source S.

The determination output unit 25 determines whether or not the estimated coordinates of the sound source S calculated by the coordinate calculation unit 24 are accurate coordinates. The contents of the convergence determination processing are as follows: The estimated coordinates for the nth time are (r ⁽ⁿ⁾ , z ⁽ⁿ⁾ ), the error from the value for the n-1th time is calculated, and this error value is obtained in advance by experiments or the like. It is determined whether or not it falls below the determination threshold. For example, as shown in the following Expression 15, the absolute value of the difference (L2 norm) is taken to normalize the distance to calculate an error value, and it is determined whether or not this error value is below the determination threshold value.

When the determination output unit 25 determines that the calculated error value exceeds the determination threshold value, it outputs an instruction to calculate the estimated coordinates of the sound source S again, and determines that the error value falls below the determination threshold value. Outputs the estimated coordinates of the sound source S calculated by the coordinate calculation unit 24 to an external output destination. Here, the instruction given by the determination output unit 25 specifically means to the sound ray parameter calculation unit 21 that the estimated coordinates (r ⁽ⁿ⁾ , z ⁽ⁿ⁾ ) obtained by the coordinate calculation unit 24 are tentative coordinates. It is to make the calculation again.

When estimating the coordinates of the unknown sound source S, the estimation obtained by performing the calculation by the sound ray parameter calculation unit 21, the propagation time difference calculation unit 22, the perturbation calculation unit 23, and the coordinate calculation unit 24 with the provisional coordinates given. Since the reliability of the coordinate value is low, the estimation value approaches the true value and the reliability is increased by repeating the calculation of the above Expression 15 a plurality of times (for example, about 10 times) to converge the error value.

Although the judgment threshold value is used as the judgment criterion by the judgment output unit 25, it is determined whether or not the number of calculations (for example, about 10 times) at which reliability is obtained as an estimated value is satisfied, or whether or not this calculation is performed. Whether or not to output the estimated coordinate may be determined based on the coordinate that gives the minimum value of Expression 15. The number of calculations is not particularly limited as long as the reliability of the obtained estimated coordinate is high. The number of times of this calculation is also obtained by an experiment or the like like the judgment threshold value.

Next, the processing operation (position estimation method of the sound source S) of the position estimation device 1 of the above-described first embodiment will be described with reference to FIG.
Note that the following processing is the processing content when the position estimation device 1 observes the first pulse wave from the sound source S that is the position estimation target.

As shown in FIG. 4, first, various data necessary for the estimation process is input from the data input unit 10 (ST1). Here, the sound velocity distribution c (z), the coordinates (r ₀ , z ₀ ) of the receiver R, and the measured propagation time difference between the direct wave dw received by the receiver R and the sea surface reflected wave sw (τ ₁₂ ^{( ob)} ) and the actual measurement value of the propagation time difference between the direct wave dw and the sea bottom reflected wave bw (τ ₁₃ ^(ob) ), the marine environment data, and the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ of the sound source S which are temporary coordinates ^). ) Is entered.

Then, sound ray parameters using any of the formula in Formula 2-5, i.e. sound ray parameter p ₁ of the direct wave dw, sound ray parameter of the sound ray parameter p ₂ and seabed reflected wave bw of sea reflected sw Calculate p ₃ (ST2).
Since the estimation process is the first time, the position of the sound source S is unknown, so the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) input by the data input unit 10 and the above-mentioned ^equations 3 to 5 are used as appropriate for each wave. The initial values (p ₁ ⁽⁰⁾ , p ₂ ⁽⁰⁾ , p ₃ ⁽⁰⁾ ) of the sound ray parameters of are calculated.

Next, the propagation time is calculated using any one of the equations 6 to 9 (ST3).
Here, the propagation time difference calculation unit 22, the propagation time T ₁ of the direct wave dw, the propagation time T ₂ and the seabed reflection of sea reflected sw calculates the propagation time T ₃ of bw.

Next, the propagation time difference calculation value τ ₁₂ which is the propagation time difference between the direct wave dw and the sea surface reflected wave sw from the propagation time (T ₁ , T ₂ , T ₃ ) of each wave obtained in ST3, and the direct wave dw and the seabed A propagation time difference calculation value τ ₁₃ which is a propagation time difference from the reflected wave bw is calculated (ST4).
Since the estimation process is the first time, the initial values τ ₁₂ ⁽⁰⁾ and τ of the propagation time difference calculation values are calculated based on the initial values (p ₁ ⁽⁰⁾ , p ₂ ⁽⁰⁾ and p ₃ ⁽⁰⁾ ) of the sound ray parameters. Calculate ₁₃ ⁽⁰⁾ .

Next, the propagation time difference actual measurement value (τ ₁₂ ^(ob) ) and the propagation time difference actual measurement value (τ ₁₃ ^(ob) ) input at the data input unit 10 and the propagation time difference calculation value τ ₁₂ ⁽⁰⁾ calculated at ST4, The perturbation amount of the estimated coordinates of the sound source S is calculated using the above ^equations 10 and 13 in which the difference from τ ₁₃ ⁽⁰⁾ is replaced by δτ (ST5).

  Next, the estimated coordinates of the sound source S are calculated using the above equation 14 based on the perturbation amount calculated in ST5 (ST6).

Next, the estimated coordinates at the n-th time are set to (r ⁽ⁿ⁾ , z ⁽ⁿ⁾ ), and the error with the value at the ⁽ n-1) -th time is calculated as shown in the above equation 15, and this error value is calculated in advance. It is determined whether or not the obtained value is below the determination threshold value (ST7).
Since the estimation procedure is the first time, the estimated coordinates (r ⁽¹⁾ , z ⁽¹⁾ ) calculated in ST6 and the initial value coordinates (r ⁽⁰⁾ , z ^{(0) of} the sound source S input by the data input unit 10 are used. ) To calculate the error value.

In ST7, if the error value is below the threshold for determination (ST7-Yes), the calculated estimated coordinates (r ⁽¹⁾ , z ⁽¹⁾ ) are output to the outside as the coordinates of the final sound source S and the processing is performed. finish.
On the other hand, when the error value exceeds the determination threshold value (ST7-No), the process returns to ST2 again to perform the second estimation process. At this time, the estimated coordinates (r ⁽¹⁾ , z ⁽¹⁾ ) calculated in the first time are updated as temporary coordinates of the sound source S, and a series of calculation processing is performed. Further, the processing loop of ST2 to ST7 is repeated a predetermined number of times while updating the provisional coordinates until the obtained error value falls below the determination threshold value.

In ST2, if the sound ray parameter p ₁ of the direct wave dw or the sound ray parameter p ₂ of the sea surface reflected wave sw does not have a real number solution, this processing is terminated without performing the subsequent processing.

[Second Embodiment]
Next, a second embodiment of the present invention will be described.
In the following description, the same components as those in the above-described first embodiment will be denoted by the same reference numerals, description of the components and processing contents will be omitted, and different components and processing functions will be mainly described.

In the position estimation device 1 of the second embodiment, in the calculation processing by the sound ray parameter calculation unit 21, when the sound ray parameter p ₁ of the direct wave dw does not have a real number solution, the rotation outside the section [z ₀ , z] is performed. This configuration takes into consideration the possibility that there are sound ray paths that reach via points. The calculation processing of the sound ray parameter p ₁ of the direct wave dw in the sound ray parameter calculation unit 21 of the second embodiment and the propagation time calculation processing of the direct wave dw in the propagation time difference calculation unit 22 are different from the processing of the first embodiment.

FIG. 5 is a schematic diagram showing an example of a turning point. As shown in the figure, at the turning point, the maximum / minimum point of the depth of the sound ray in which the direction of the sound ray changes from bottom to top or from top to bottom is represented by z _tp , and pc (z) = 1 is satisfied at the turning point. . In FIG. 5, (a) shows the minimum point of the turning point, and (b) shows the maximum point. Then, assuming that the turning point depth is z _tp as illustrated, the r coordinate of the sound source S is given by the following Expression 16.

The sound ray parameter calculation unit 21 of the second embodiment assumes that the direct wave dw has passed through the turning point when the value of the sound ray parameter p ₁ does not have a real number solution, and instead of the above Equation 3, the above equation 16 is used to find the sound ray parameter p _1tp of the direct wave dw. In the calculation process, if the estimation process is the first time, the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) are used as the temporary coordinates, and if the second time or later, the estimated coordinates (r ^{( n-1)} , z ^(n-1) ) are used as temporary coordinates for calculation.

Further, the propagation time difference calculation unit 22 of the second embodiment assumes that the direct wave dw has passed through the turning point when the value of the sound ray parameter p ₁ does not have a real number solution, and instead of the above Equation 7, The propagation time T ₁ of the direct wave dw is obtained by using Expression 17.

  Note that the sound ray parameter calculation unit 21 in the second embodiment performs the above-described processing about a few times while determining whether the direct wave dw is a real number solution, and the sound ray parameter of the direct wave dw is a real number. If there is no solution, it is assumed that the direct wave dw does not exist, and the coordinate estimation process of the sound source S, which is the process of the present device, ends.

As the processing operation of the position estimation device 1 of the second embodiment, it is assumed that the direct wave dw passes through the turning point when the direct wave dw does not have a real solution in the processing of the sound ray parameter calculation unit 21 as described above. The processing corresponding to it is performed.
That is, the processing operation of the second embodiment is basically the same as the processing operation of the first embodiment described above, although the mathematical expressions used are slightly different, and the sound ray parameter processing of the direct wave dw in ST2 is different.

In the sound ray parameter calculation process (ST2) in the second embodiment, when the sound ray parameter p ₁ of the direct wave dw does not have a real number solution, the sound ray parameter p _1tp is calculated using the above equation 16. At this time, if the estimation process is the first time, the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) are used as the temporary coordinates, and if it is the second time or later, the estimated coordinates (r ^{(n-1) )} , Z ^(n-1) ) as temporary coordinates.

Then, if the sound ray parameter p _1tp of the direct wave dw assumed to _pass through the turning point becomes a real number, the process proceeds to the next step (ST4). When the calculated result does not have a real number solution, the sound ray parameter p _1tp of the direct wave dw is calculated again. As a result of repeating this process a few times, when the sound ray parameter p _1tp of the direct wave dw does not have a real number solution, it is determined that the direct wave dw does not exist, and the subsequent processes are ended.

[Third Embodiment]
Next, a third embodiment of the present invention will be described.
In the following description, the same configurations as those of the above-described first and second embodiments will be denoted by the same reference numerals, description of the configurations and processing contents will be omitted, and different configurations and processing functions will be mainly described.

When the seabed is uneven, the sound ray path of the seabed reflected wave bw is different from when the seabed is flat. Therefore, in the third embodiment, when the seabed has undulations and the seabed topography is known and can be obtained by the function h (r), the process of calculating the sound ray parameter p ₃ of the seabed reflected wave bw and the propagation time T ₃ is performed. Different from the first and second embodiments.

FIG. 6 is a schematic diagram showing the relationship between the incident elevation angle and the reflected depression angle of a sound wave on the undulating sea bottom. The seabed topography based on the receiver R is known and is given by the function h (r). At this time, the reflected depression / elevation angle θ _out of the sound wave incident on the sea bottom h (r) having the inclination angle α = dh / dr with respect to the horizontal plane at the depression angle θ _in is given by the following formula 18.

  Therefore, since the sound ray parameter changes before and after the reflection, the sound ray parameter after the reflection is given by the above Expression 1 and Expression 18 from the following Expression 19. Here, the r coordinate from the receiving point to the seabed reflection point is expressed by the following formula 20, and the sound velocity at the seabed reflection point is expressed by the following formula 21.

Further, regarding the characteristic sound line of the seabed reflected wave bw that passes through the receiving point, the seabed reflection point, and the sound source S, the distance to the sound source S is given by the following formula 22, so the propagation time T ₃ of the seabed reflected wave bw is Required from 23.

As described above, in the third embodiment, since the sound ray parameter p ₃ when the seabed has undulations is determined, the temporary coordinates of the sound source S (for example, initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ )). Also, the initial value of the sound ray parameter p ₃ is obtained by calculating the above equation 20 based on the virtual horizontal distance from the wave receiving point to the seabed reflection point (virtual distance obtained in advance by experiments or the like). In addition, the value of the horizontal distance from the receiving point to the seafloor reflection point is updated by solving the above equation 22, and the calculation of the above equations 20 and 22 is repeated, and the error is a threshold value obtained in advance by experiments or the like. The horizontal distance is determined when it falls within the range or reaches a predetermined number of calculations (about 10 times). Therefore, not only the sound ray parameter p ₃ of the seabed reflected wave bw but also the coordinates of the seabed reflection point shown in the following Expression 24 can be obtained.

The sea bottom reflection point can also be obtained by calculating the intersection of the sound ray represented by the sound ray parameter p ₃ and the sea bottom.

Therefore, when the sound ray parameter calculation unit 21 of the third embodiment determines from the ocean environment data from the data input unit 10 that the seafloor topography is not flat and is changing, it uses Equation 22 instead of Equation 5. If the estimation process is the first time, the initial value coordinates (r ⁽⁰⁾ , z ⁽⁰⁾ ) are used as temporary coordinates, and if the second and subsequent times, the estimated coordinates (r ^(n-1) , obtained by the coordinate calculation unit 24 are calculated. The sound ray parameter p ₃ is calculated using z ^(n-1) ) as temporary coordinates.

Further, when the propagation time difference calculation unit 22 of the third embodiment determines that the seafloor topography is not flat and is changing, the propagation time T of the seafloor reflected wave bw is calculated by using the above Formula 23 instead of the above Formula 9. Ask for ₃ .

Note that the sound ray parameter calculation unit 21 of the third embodiment also executes the processing of the above-described second embodiment when the value of the sound ray parameter p ₁ does not have a real number solution. That is, when the value of the sound ray parameter p ₁ does not have a real number solution, it is estimated that the direct wave dw has passed through the turning point, and instead of the above Equation 3, the above Equation 16 is used to calculate the sound ray of the direct wave dw. Find the parameter p ₁ .

In addition, the propagation time difference calculation unit 22 of the third embodiment also executes the processing of the above-described second embodiment when the value of the sound ray parameter p ₁ does not have a real number solution. That is, the propagation time T ₁ of the direct wave dw is obtained by using the sound ray parameter p ₁ of the direct wave dw passing through the turning point and using the above expression 17 instead of the above expression 7.

As the processing operation of the position estimation device 1 of the third embodiment, the sound ray parameter p ₃ of the seabed reflected wave bw that reaches the receiver R from the sound source S via the seabed reflection point when the seabed has undulations and propagation. Corresponding processing for obtaining the time T ₃ is performed.

Therefore, the processing operation of the third embodiment is basically the same as the processing operation of the first embodiment described above, although the mathematical expressions used are slightly different, and the sound ray parameter of the seabed reflected wave bw is acquired in the processing of ST2. In addition to the calculation processing of the sound ray parameter p ₃ in consideration of the seabed topography, the processing of calculating the coordinates of the seabed reflection point is added.

[Fourth Embodiment]
Next, a fourth embodiment of the present invention will be described.
In the following description, the same configurations as those of the above-described first to third embodiments will be denoted by the same reference numerals, description of the configurations and processing contents will be omitted, and different configurations and processing functions will be mainly described.

Depending on the depth of the sound source S, there is a region where the direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw do not reach the receiver R. FIG. 7 is a schematic diagram showing an example of a sound ray of a pulse wave emitted from the sound source S and a shadow zone of the sound ray.
As shown in the figure, the pulse wave emitted from the sound source S propagates while refracting up and down in the sea, but the pulse wave does not propagate while refracting in the entire undersea area, and in some areas (shadow zones) where sound waves do not reach. Exists. In other words, if the wave receiver R is in the shadow zone, the pulse wave required for implementing the first to third embodiments cannot be observed, and thus the coordinates of the sound source S due to the propagation time difference cannot be estimated.

  Therefore, in the fourth embodiment, when the depth of the sound source S or the distance from the wave receiver R to the sound source S is approximately known, any one of the direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw is propagated. The r coordinate or z coordinate of the sound source S is estimated by observing the angle θ. In the first to third embodiments described above, the coordinate estimation processing is performed when the coordinates (distance, depth) of the sound source S are both unknown, but in the fourth embodiment, the depth of the sound source S or the sound source S from the wave receiver R is calculated. This is coordinate estimation processing when the distance to is known.

Letting θ be the observed value of the propagation angle of elevation of a certain pulse wave, the sound ray parameter p is given by the following equation 25 from the above equation 1. Here, the observation value θ of the propagation depression angle is obtained by arranging a plurality of receivers R in parallel in the sea and correcting the reception time (phase shift) of the pulse wave received by each receiver R.

  Further, in the fourth embodiment, in order to identify which of the direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw is the pulse wave from the sound source S, the pulse received from above the receiver R is identified. The wave is identified as the direct wave dw or the sea surface reflected wave sw, and the pulse wave received from below the receiver R is identified as the seabed reflected wave bw. Therefore, in the fourth embodiment, it is necessary to identify what kind of wave the pulse wave from the sound source S is. Therefore, in the present embodiment, the processing can be performed when it can be identified as the seabed reflected wave bw.

When the depth z of the sound source S is approximately known, the sound ray parameter calculation unit 21 of the fourth embodiment specifies the type of sound ray by the above processing, and when it can be identified as the seabed reflected wave bw, the seabed reflected wave. By substituting the above equation 25 into the above equation 5 for obtaining the r coordinate of bw, the r coordinate of the sound source S is obtained. If the r-coordinate is obtained, the coordinates of the sound source S can be naturally obtained. That is, the sound ray parameter p ₃ of the seabed reflected wave bw is obtained based on the actual measurement value of the propagation angle of elevation θ of the seabed reflected wave bw at the coordinate of the receiver R and the value of the sound velocity, and the sound source S whose depth coordinate is known is obtained. The coordinates of the sound source S can be estimated by determining the distance of.

When the distance r to the sound source S is roughly known, the sound ray parameter calculation unit 21 of the fourth embodiment specifies the type of sound ray by the above processing, and when it can be identified as the seabed reflected wave bw, the seabed reflection The above-mentioned equation 25 is substituted into the above-mentioned equation 5 which gives the r-coordinate of the wave bw, and the z-coordinate of the sound source S is obtained by numerically solving the z-coordinate. If the z coordinate is obtained, the coordinate of the sound source S can be obtained naturally. That is, the sound ray parameter p ₃ of the seabed reflected wave bw is obtained based on the actual measurement value of the propagation angle of elevation θ of the seabed reflected wave bw at the coordinate of the receiver R and the value of the sound speed, and the distance coordinate of the sound source S whose distance coordinate is known is obtained. The coordinates of the sound source S can be estimated by determining the depth.

  As described above, the position estimation device 1 according to the present invention receives the direct wave dw, the sea surface reflected wave sw, and the seabed reflected wave bw of the pulse wave generated by the sound source S with the receiver R as a starting point, and The propagation time difference is calculated from the sound ray parameters, the perturbation amount of the difference between the propagation time difference calculation and the actual measurement value of the propagation time difference in each received wave is obtained, and the error value of the estimated coordinate obtained based on this perturbation amount is converged. The coordinates of the sound source S are estimated.

  For this reason, unlike the conventional brute force method, it is not necessary to grid the entire ocean space in Cartesian coordinates or polar coordinates and calculate it using a numerical solution method using a differential equation. Can be estimated.

DESCRIPTION OF SYMBOLS 1 ... Position estimation device 10 ... Data input part 20 ... Arithmetic processing part 21 ... Sound ray parameter calculation part 22 ... Propagation time difference calculation part 23 ... Perturbation calculation part 24 ... Coordinate calculation part 25 ... Judgment output part R ... Wave receiver S ... Sound source dw ... direct wave sw ... sea surface reflected wave bw ... sea bottom reflected wave
W1 ... Sea surface
W2 ... sea bottom

A sixth aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed is flat,
Sound ray parameters of a direct wave, a sea surface reflected wave, and a seabed reflected wave that are eigen sound rays of the pulse wave that passes through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source. And the step of finding
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
It is a position estimation method characterized by including.

A seventh aspect of the present invention is a position estimation method for estimating the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed undulates,
A step of obtaining sound ray parameters of a direct wave and a sea surface reflected wave which are characteristic sound rays of the pulse wave passing through two points of temporary coordinates of the sound source and coordinates of a receiver for receiving the pulse wave emitted by the sound source;
Temporary coordinates of the sound source, and a step of obtaining a sound ray parameter of the seabed reflected wave from the virtual horizontal distance from the receiver to the seabed reflection point,
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
It is a position estimation method characterized by including.

An eighth aspect of the present invention is the position estimation method according to the sixth or seventh aspect of the present invention, wherein the sound ray parameter uses the coordinates of the wave receiver as a starting point and uses temporary coordinates of the sound source. It is a position estimation method characterized by being obtained by a numerical solution of a nonlinear equation.

A ninth aspect of the present invention is the position estimation method according to any one of the sixth to eighth aspects of the present invention, wherein when the sound ray parameter of the direct wave does not have a real number solution, the sound of the direct wave is The position estimating method is characterized in that the sound ray parameter of the direct wave is obtained based on the turning point depth on the assumption that the line passes through the turning point.

A tenth aspect of the present invention is, in the position estimation method according to any one of the sixth to eighth aspects of the present invention, a known seabed topographical coordinate and a virtual horizontal distance from the wave receiver to the seabed reflection point. The position estimating method is characterized in that the sea bottom reflection point coordinates of the sea bottom reflected wave and the sound ray parameter of the sea bottom reflected wave are obtained based on the above.

Claims (14)

In the position estimation device that estimates the coordinates of the sound source existing in the known ocean where the speed of sound changes in the depth direction and the seabed is flat,
Sound ray parameters of a direct wave, a sea surface reflected wave, and a seabed reflected wave that are eigen sound rays of the pulse wave that passes through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source. And the process of
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A process for obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
A process of obtaining estimated coordinates of the sound source estimated from the perturbation amount;
When the error value between the estimated coordinate of the sound source obtained this time and the estimated coordinate obtained last time is less than a preset threshold for determination, a process of determining the estimated coordinate obtained this time as the coordinate of the sound source,
A position estimation device comprising an arithmetic processing unit for performing the following. In the position estimation device that estimates the coordinates of the sound source existing in the known ocean where the speed of sound changes in the depth direction and the seabed undulates,
A process of obtaining sound ray parameters of a direct wave and a sea surface reflected wave that are eigen sound rays of the pulse wave passing through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source;
Temporary coordinates of the sound source, a process of obtaining sound ray parameters of the seabed reflected wave from the virtual horizontal distance from the receiver to the seabed reflection point,
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A process for obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a process of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
A position estimation device comprising an arithmetic processing unit for performing the following.   3. The position estimating apparatus according to claim 1, wherein the sound ray parameter is obtained by a numerical solution of a non-linear equation using a coordinate of the wave receiver as a starting point and a temporary coordinate of the sound source. When the sound ray parameter of the direct wave does not have a real number solution, it is assumed that the sound ray of the direct wave passes through a turning point, and the sound ray parameter of the direct wave is obtained based on the turning point depth. The position estimation device according to claim 1, wherein the coordinates of the sound source are estimated.   The seabed reflection point coordinates of the seabed reflected wave and the sound ray parameter of the seabed reflected wave are obtained based on known seabed topographic coordinates and a virtual horizontal distance from the receiver to the seabed reflected point. The position estimation device according to any one of 2 to 4. In the position estimation device that estimates the coordinates of the sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver for receiving the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the depth coordinate is known. A position estimating device for estimating coordinates of a sound source by determining a distance to the sound source. In the position estimation device that estimates the coordinates of the sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver that receives the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the distance coordinates are known. A position estimating device for estimating coordinates of a sound source by determining a depth of the sound source. In the position estimation method for estimating the coordinates of a sound source existing in a known ocean where the speed of sound changes in the depth direction and the seabed is flat,
Sound ray parameters of a direct wave, a sea surface reflected wave, and a seabed reflected wave that are eigen sound rays of the pulse wave that passes through two points of the temporary coordinates of the sound source and the coordinates of a receiver that receives the pulse wave emitted by the sound source. And the step of finding
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
Determining estimated coordinates of the sound source estimated from the perturbation amount,
When the error value between the estimated coordinate of the sound source obtained this time and the estimated coordinate obtained last time is below a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source,
A position estimation method comprising: In the position estimation method for estimating the coordinates of the sound source existing in the known ocean where the speed of sound changes in the depth direction and the seabed undulates,
A step of obtaining sound ray parameters of a direct wave and a sea surface reflected wave which are characteristic sound rays of the pulse wave passing through two points of temporary coordinates of the sound source and coordinates of a receiver for receiving the pulse wave emitted by the sound source;
Temporary coordinates of the sound source, and a step of obtaining a sound ray parameter of the seabed reflected wave from the virtual horizontal distance from the receiver to the seabed reflection point,
By the perturbation method, the direct time and the sea surface reflected wave and the propagation time difference measured value of the direct wave and the seabed reflected wave, the direct wave and the sea surface reflected wave and the direct wave obtained from the sound ray parameter of each wave A step of obtaining a perturbation amount of the difference between the propagation time difference calculation value of the seabed reflected wave,
When the error value between the estimated coordinate of the sound source estimated from the perturbation amount and the estimated coordinate obtained last time is less than a preset threshold for determination, a step of determining the estimated coordinate obtained this time as the coordinate of the sound source, ,
A position estimation method comprising:   10. The position estimating method according to claim 8, wherein the sound ray parameter is obtained by a numerical solution of a nonlinear equation using the coordinates of the wave receiver as a starting point and the temporary coordinates of the sound source. When the sound ray parameter of the direct wave does not have a real number solution, it is assumed that the sound ray of the direct wave passes through a turning point, and the sound ray parameter of the direct wave is obtained based on the turning point depth. The position estimation method according to any one of claims 8 to 10, which is characterized.   The seabed reflection point coordinates of the seabed reflected wave and the sound ray parameter of the seabed reflected wave are obtained based on known seabed topographic coordinates and a virtual horizontal distance from the receiver to the seabed reflected point. The position estimation method according to any one of 8 to 10. In the position estimation method that estimates the coordinates of the sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver for receiving the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the depth coordinate is known. A position estimating method comprising estimating coordinates of the sound source by determining a distance to the sound source. In the position estimation method that estimates the coordinates of the sound source existing in the ocean,
The actual measurement value of the propagation angle of elevation and the propagation angle of the seabed reflected wave at the coordinates of the receiver that receives the pulse wave emitted by the sound source, and the sound ray parameter of the seabed reflected wave based on the value of the sound velocity, the distance coordinates are known. A position estimating method comprising estimating coordinates of the sound source by determining a depth to the sound source.

Images