Method for determining the direction of a source of waterborne sound, a computer program product, a computer or sonar, and a watercraft

Abstract

The invention relates to a method for determining the direction of a source of waterborne sound that emits a waterborne acoustic signal, by means of a hydrophone arrangement which forms a linear antenna or a virtual linear antenna, as well as to a computer program product, a computer, a sonar, and a watercraft.

Claims

1. A method for determining the direction of a source of waterborne sound that emits a waterborne acoustic signal by means of a hydrophone arrangement which forms a linear antenna or a virtual linear antenna, wherein the method comprises the following steps: determining the waterborne acoustic signal incoming at the hydrophone arrangement at a first measurement time t.sub.1, performing an electronic antenna swing analysis which determines a sum signal S of the linear antenna or of the virtual signal antenna at every rotation angle n=0n, where 0 is the starting angle, is the step-width angle, and n=0, 1, 2, 3, 4 . . . , such that a measurement series S() is created, determining a maximum sum signal, S.sub.max(.sub.max), where .sub.max represents the maximum rotation angle of the measurement series S(), determining a fit measurement series S.sub.Fit(), which includes the maximum sum signal S.sub.max(.sub.max) and at least two more sum signals S(.sub.n).sub.m, with the sum number m=2, 3, 4 . . . , determining a mathematical fit function S.sub.Fit(.sub.Fit) for the fit measurement series S.sub.Fit() by means of a mathematical equalization calculation, and determining a fit function maximum f.sub.Fit,Max(.sub.Fit,Max) of the mathematical fit function f.sub.Fit(.sub.Fit), wherein the direction of the waterborne sound source can be determined by means of .sub.Fit,Max.

2. Method according to claim 1, wherein the procedure or sub-steps of the procedure can be carried out at a second measurement time t2, at a third measurement time t3, at a fourth measurement time t4, and/or at further times t, wherein the corresponding measurement series S().sub.t or the determined directions of the waterborne sound source are averaged.

3. Method according to claim 2, wherein the remaining steps of the procedure are carried out according to claim 1.

4. Method according to claim 1, wherein the additional sum signals are sum signals adjacent to the maximum sum signal S.sub.max(.sub.max).

5. Method according to claim 1, wherein the mathematical fit function is a polynomial function, a Gaussian function, a Lorentz function, or a Voigt function.

6. An invention selected from the group consisting of: (a) a computer program product which is set up so that the method according to claim 1 can be implemented with a computer or a sonar; (b) a computer configured for carrying out a method according to claim 1; (c) a sonar system for determining a direction of a source of a waterborne sound, the system comprising: an antenna, wherein the antenna includes a plurality of piezoceramic elements, wherein the piezoceramic elements are configured to: detect the waterborne acoustic sound; and generate, based on the waterborne acoustic sound, a voltage; a control circuit communicatively coupled to the plurality of piezoceramic elements, wherein the control circuit is configured to: receive, from the plurality of piezoceramic elements at multiple points in time, the voltage; perform an electronic swing analysis at a plurality of rotation angles at the multiple points in time; determine, based on the electronic swing analysis, a maximum sum signal; determine a fit measurement series, wherein the fit measurement series includes the maximum sum signal and at least two additional sum signals; determine, based on a mathematical equalization calculation and the fit measurement series, a mathematical fit function; and determine, based on the mathematical fit function, the direction of the source of the waterborne sound; and (d) a sonar system configured for carrying out a method according to claim 1.

7. An invention according to claim 6, wherein the invention is (b) a computer configured for carrying out a method according to claim 1.

8. An invention according to claim 6, wherein the invention is (d) a sonar system configured for carrying out a method according to claim 1.

9. A watercraft which includes a sonar system configured for carrying out a method comprising the following steps: determining the waterborne acoustic signal incoming at the hydrophone arrangement at a first measurement time t.sub.1, performing an electronic antenna swing analysis which determines a sum signal S of the linear antenna or of the virtual signal antenna at every rotation angle n=0n, where 0 is the starting angle, is the step-width angle, and n=0, 1, 2, 3, 4 . . . , such that a measurement series S() is created, determining a maximum sum signal S.sub.max(.sub.max), where .sub.max represents the maximum rotation angle of the measurement series S(), determining a fit measurement series S.sub.Fit(), which includes the maximum sum signal S.sub.max(.sub.max) and at least two more sum signals S(.sub.n).sub.m, with the sum number m=2, 3, 4 . . . , determining a mathematical fit function f.sub.Fit(.sub.Fit) for the fit measurement series S.sub.Fit() by means of a mathematical equalization calculation, and determining a fit function maximum f.sub.Fit,Max(.sub.Fit,Max) of the mathematical fit function f.sub.Fit(.sub.Fit), wherein the direction of the waterborne sound source can be determined by means of .sub.Fit,Max.

10. An invention according to claim 6, wherein the invention is (c) a sonar system for determining a direction of a source of a waterborne sound, the system comprising: an antenna, wherein the antenna includes a plurality of piezoceramic elements, wherein the piezoceramic elements are configured to: detect the waterborne acoustic sound; and generate, based on the waterborne acoustic sound, a voltage; a control circuit communicatively coupled to the plurality of piezoceramic elements, wherein the control circuit is configured to: receive, from the plurality of piezoceramic elements at multiple points in time, the voltage; perform an electronic swing analysis at a plurality of rotation angles at the multiple points in time; determine, based on the electronic swing analysis, a maximum sum signal; determine a fit measurement series, wherein the fit measurement series includes the maximum sum signal and at least two additional sum signals; determine, based on a mathematical equalization calculation and the fit measurement series, a mathematical fit function; and determine, based on the mathematical fit function, the direction of the source of the waterborne sound.

11. An invention according to claim 10, wherein in (c) the at least two additional sum signals are adjacent to the maximum sum signal.

12. An invention according to claim 11, wherein in (c) the at least two additional sum signals are symmetric about the maximum sum signal.

13. An invention according to claim 6, wherein in (c) the mathematical fit function is one of a polynomial function, a Gaussian function, a Lorentz function, and a Voigt function.

14. An invention according to claim 6, wherein in (c) the plurality of piezoceramic elements includes between 50 and 200 piezoceramic units.

15. An invention according to claim 6, wherein in (c) the plurality of piezoceramic elements form a side-scan sonar antenna.

16. An invention according to claim 15, wherein the side-scan sonar antenna is one of a linear antenna and a virtual linear antenna.

17. An invention according to claim 6, wherein in (c) there is a common step width between each of the plurality of rotation angles.

18. An invention according to claim 17, wherein in (c) the common step width is approximately two degrees.

19. An invention according to claim 6, wherein in (c) a total rotation for the plurality of rotation angles is one hundred eighty degrees.

20. An invention according to claim 6, wherein the invention is () a computer program product which is set up so that the method according to claim 1 can be implemented with a computer or a sonar.

21. A ship or a submarine having a sonar system for determining a direction of a source of a waterborne sound, wherein the sonar system comprises: an antenna, wherein the antenna includes a plurality of piezoceramic elements, wherein the piezoceramic elements are configured to: (a) detect the waterborne acoustic sound; and (b) generate, based on the waterborne acoustic sound, a voltage; a control circuit communicatively coupled to the plurality of piezoceramic elements, wherein the control circuit is configured to: (i) receive, from the plurality of piezoceramic elements at multiple points in time, the voltage; (ii) perform an electronic swing analysis at a plurality of rotation angles at the multiple points in time; (iii) determine, based on the electronic swing analysis, a maximum sum signal; (iv) determine a fit measurement series, wherein the fit measurement series includes the maximum sum signal and at least two additional sum signals; (v) determine, based on a mathematical equalization calculation and the fit measurement series, a mathematical fit function; and (vi) determine, based on the mathematical fit function, the direction of the source of the waterborne sound.

22. A ship or a submarine according to claim 21, wherein the at least two additional sum signals are symmetric about the maximum sum signal.

23. A ship or submarine according to claim 21, wherein the plurality of piezoceramic elements form a side-scan sonar antenna.

Description

BRIEF DESCRIPTIONS OF THE DRAWINGS

(1) FIG. 1 is a very schematic representation of a side-scan sonar (flank array) with 151 individual hydrophones, which is arranged on a submarine.

(2) FIG. 2 is a very schematic representation of a computational rotation of the side-scan sonar from FIG. 1 and the underwater acoustic signals emitted by a waterborne sound source.

(3) FIG. 3 shows a functional relationship between a cumulated signal strength S and a rotation angle , as well as three associated measured values, and the determined mathematical fit function, as well as a global maximum associated with the mathematical fit function and the rotation angle which can be derived from it.

DETAILED DESCRIPTION

(4) It will be useful to clarify certain terms:

(5) By determination of direction is meant determining the direction of the underwater sound source with respect to a reference point for example, the center and/or the pivot point of the linear antenna. In this context, the term bearing is also used instead of the term direction. Direction and bearing are to be understood here as synonymous.

(6) A waterborne sound source is, in particular, a natural or artificial object which emits sound waves from a localized region. This could be, for example, a ship's engine and/or a point source. In addition, a waterborne sound source is also an artificial or natural object which reflects a previously-emitted, underwater acoustic signal. For example, an underwater vehicle lying on the lookout, which is pinged by waterborne sound from an active sonar of another watercraft and thus reflects the waterborne acoustic signal, is here understood as a waterborne sound source.

(7) A waterborne acoustic signal is, in particular, an acoustic transient signal in the water with a time duration of 10 ms to 200 ms, which can extend, in particular, to a range from 2 Hz to several hundred kHz, with a signal shape which is frequency-modulated and/or has a constant wavelength.

(8) A hydrophone arrangement includes, in particular, at least two piezoceramic elements in which the waterborne acoustic signal creates a voltage change which is processed by signal technology. The more hydrophones are fitted, the better can frequencies and/or directions be resolved. In the present case, such a hydrophone arrangement comprises, in particular, between 50 and 300 hydrophones.

(9) A linear antenna is, in particular, present when the individual hydrophones of the hydrophone arrangement are arranged along a straight line. This is, for example, the case with a side-scan sonar. However, even curved hydrophone arrangements or, for example, hydrophone arrangements with a cylindrical basis can form a linear antenna, which in the present case is referred to as a virtual linear antenna, since such arrangements are converted by signal technology or by computer technology into a linear antenna. In this regard, see Manfred Siegel, Einfhrung in die Physik und Technik der Unterwasserschallsysteme [Introduction to the Physics and Technology of Underwater Sound Systems] (ISBN 3-936799-28-8), pages 5-52 to 5-60, the relevant content of which forms a constituent part of the present application.

(10) A first measurement time t.sub.1 covers one time period or one time interval, so that the antenna swing analysis can be carried out. During this first measuring time, the voltage, in particular, present at each hydrophone due to the waterborne acoustic signal is determined by measurement technology. Consequently, a measurement signal which varies over the measurement time is present for the measurement time.

(11) An electronic antenna swing analysis is a delay which differs computationally or electronically for each individual hydrophone, so that the linear antenna is rotated or curved electronically or computationally. In this regard, see Manfred Siegel, Einfhrung in die Physik und Technik der Unterwasserschallsysteme [Introduction to the Physics and Technology of Underwater Sound Systems] (ISBN 3-936799-28-8), pages 5-41 to 5-60, the relevant content of which forms a constituent part of the present application. This rotation or curvature is effected, in particular, by the individual hydrophones being read with a different time delay, so that time-delayed hydrophones are obtained.

(12) The sum signal S is, in particular, the cumulative signal of all (time-delayed) hydrophones. The sum signal can be determined, not only for a single frequency, but also for a single frequency interval. Furthermore, the signal levels at the hydrophones can simply be added together. The sum signal is thus a measure of the sound energy present at the linear antenna or the virtual linear antenna. The sum signal is determined at each individual rotation angle .sub.n=.sub.0n so that the sum signal S is present at each rotation angle .sub.n.

(13) The rotation angle is created by the rotation under electronic or computer control of the linear antenna or the virtual linear antenna about a pivot point. In particular, the antenna is rotated in each case with a step-width angle of =180/hydrophones_number/k (in degrees; full circle 360) (k=3, 5, 7 . . . is the fine bearing factor), (for example, for an antenna with 18 hydrophones and k=5, will be 2), so that, with an angular segment of 180, a total of 90 sum signal angles S will be determined. These 90 sum signal values S, together with the associated rotation angles, will then, in particular, be saved under measurement series S().

(14) In order to determine the maximum sum signal S.sub.max(.sub.max), the quantitatively largest sum signal S from the measurement series S() is determined. This can be done with, for example, a so-called quicksort or bubble sort algorithm. The maximum rotation angle .sub.max is then the angle of the maximum signal S.sub.max thus determined.

(15) To determine the fit measurement series S.sub.Fit(), at least two more sum signals S(.sub.n).sub.m are determined in addition to the maximum sum signal S.sub.max(.sub.max). The m subscript here denotes the number of additional sum signals present in the fit series S.sub.Fit() alongside the maximum sum signal S.sub.max(.sub.max).

(16) To achieve a particularly good result, the sum signals directly adjacent to the maximum sum signal S.sub.max(.sub.max), in particular, are also included in the fit series. In particular, the sum signals symmetrically flanking the maximum sum signal S.sub.max(.sub.max) are used. The fit measurement series S.sub.Fit() thus includes at least three measured values. In case of the symmetrical use of additional sum signals, the fit measurement series 5, 7, 9, . . . includes sum signals.

(17) A mathematical fit function S.sub.Fit() is determined for the measured values saved in the fit measurement series S.sub.Fit(), including the maximum sum signal S.sub.max(.sub.max). The determination of a mathematical fit function of this kind is also referred to as regression. For example, the assumption is made that three measured values saved in the fit measurement series can be described by a second-degree polynomial (.sub.0+.sub.1x+.sub.2x.sup.2). The parameters of the second-degree polynomial are thus defined such that, ideally, all three points lie on the determined mathematical fit function.

(18) In particular, it may be the case here that the metrologically-determined, maximum rotation angle does not agree with a maximum rotation angle of a local minimum of the mathematical fit function, but it also, surprisingly, turns out that the local or global maximum often results in a much more correct rotation angle for determining the direction of an underwater sound source. In this way, a fit measurement function maximum f.sub.Fit,Max(.sub.Fit,Max) of the mathematical fit function f.sub.Fit(.sub.Fit) is determined, in particular, numerically or analytically. In particular, the direction of the waterborne sound source is orthogonal to a linear antenna or virtual linear antenna with the rotation angle .sub.Fit,Max determined by the fit function, wherein the orthogonal here intersects the pivot point, in particular, of the linear antenna or virtual linear antenna.

(19) In order to improve the quality and, thus, the reliability of the direction of the waterborne sound source, the procedure or sub-steps of the procedure can be carried out at a second measurement time t.sub.2, at a third measurement time t.sub.3, at a fourth measurement time t.sub.4, and/or at further times t, wherein the corresponding measurement series S().sub.t or the determined directions of the waterborne sound source are averaged.

(20) This approach delivers excellent results, particularly in the case where the waterborne sound source is not moving. The corresponding measurement times emerge from the corresponding repetition rate, in particular, at which the hydrophone arrangement is interrogated electronically or computationally, or from an overlaid window function which is superimposed as a filter on the corresponding measurement curves of the hydrophone.

(21) To also allow determination of the direction of the waterborne sound source at the end of the procedure, in this embodiment, the remaining steps of the method according to the invention are also carried out.

(22) In another embodiment, the additional sum signals are sum signals adjacent to the maximum sum signal S.sub.max(.sub.max).

(23) This makes it possible to obtain a high measured value density around the maximum sum signal, so that the area around the maximum signal is particularly well-reproduced by the mathematical fit function.

(24) In order to compensate for antenna-specific properties, the mathematical fit function can be a polynomial function, a Gaussian function, a Lorentz function, or a Voigt function.

(25) Both the Gaussian function and the Lorentz function describe bell curves, while, for example, the Voigt function is, in particular, a convolution of the Gaussian function and the Lorentz function in which, for example, broadening effects on the underwater sound due to the water and its properties can be compensated for.

(26) In another aspect, the aim is achieved by a computer program product which is set up so that the previously described method can be implemented with a computer or a sonar.

(27) A software application can thus be provided on a data medium, said software sending instructions to a computer or a sonar such that the above-described method according to the invention can be implemented. In this way, existing computer systems and/or sonars can therefore be adapted to enable them to carry out the procedure.

(28) In another aspect, the aim is achieved by a computer which is set up in such a way that a previously described method can be implemented.

(29) In addition, in another aspect, the aim is achieved by a sonar, which is set up in such a way that a previously described method can be implemented, or which includes a previously described computer.

(30) Finally, the aim is achieved by a watercraft equipped with a sonar as previously described. Watercraft of this kind may, in particular, be ships or submarines.

(31) The invention is explained in more detail below with reference to an exemplary embodiment. The figures show the following:

(32) FIG. 1 a very schematic representation of a side-scan sonar (flank array) with 151 individual hydrophones, which is arranged on a submarine,

(33) FIG. 2 a very schematic representation of a computational rotation of the side-scan sonar from FIG. 1 and the underwater acoustic signals emitted by a waterborne sound source, and

(34) FIG. 3 a functional relationship between a cumulated signal strength S and a rotation angle , as well as three associated measured values, and the determined mathematical fit function, as well as a global maximum associated with the mathematical fit function and the rotation angle which can be derived from it.

(35) A side-scan sonar antenna 101 is arranged on a submarine 111. The side-scan sonar antenna 101 includes individual hydrophones 103, which have a piezoceramic element. The pressure differences impressed on the piezoceramic elements by the waterborne acoustic signal generate a voltage, which is metrologically processed and evaluated. This processing and evaluation takes place almost continuously. This results over time in an almost continuous signal voltage curve at a hydrophone. These voltage curves are convoluted with a temporal window filter, so that a measurement interval results.

(36) The individual hydrophone signals are evaluated in each case at different times within the temporal filter window, so that what is effectively a computational rotation of the side-scan sonar antenna 101 results.

(37) In the present case, the step width for rotation is =2. In addition, a cumulated sum signal covering all hydrophones 103 is formed for each angle formed thereby. The rotation is in the direction of rotation 225. Here, the rotation is about the pivot point 227, which in the present case is the center of the hydrophonesin this case, the central hydrophone with the number 76. As soon as the (rotated) side-scan sonar antenna 223 is rotated (.sub.max), such that it essentially emits in parallel to the propagation waves 231 of the propeller signal 229, the maximum sum signal S.sub.max is obtained.

(38) The direct sum signals around the maximum sum signal S.sub.max at the rotation angle .sub.max are extracted and fitted mathematically by means of a second-degree polynomial.

(39) FIG. 3 shows not only the extracted measured values 333 with the corresponding rotation angles , but also a graphical representation of the determined fit function 341.

(40) Here, the rotation angle is shown in the graph on the independent axis, and the sum signal S on the dependent axis.

(41) The maximum .sub.fit,max is identified for the determined mathematical fit function 341. This differs somewhat from the value (.sub.max) determined by the rotation of the side-scan sonar. The maximum value .sub.fit,max thus obtained is a marked improvement upon the angular value .sub.max and gives the direction from the pivot point 227 to the point location of the sound source 229, thereby providing a bearing.

LIST OF REFERENCE SYMBOLS

(42) 101 Side-scan sonar antenna 103 Hydrophones 111 Submarine 221/223 Antenna pivoted electronically to antenna 101 225 Direction of rotation 227 Pivot point 229 Underwater sound sources 231 Underwater sound waves 333 Measured values 341 Graph of a fit function