This application is based upon and claims the benefit of priority from Japanese Patent Application No. 2016-123200, filed on Jun. 22, 2016, the disclosure of which is incorporated herein in its entirety by reference.
The disclosed subject matter relates to an active sonar and a control method for the active sonar.
In underwater monitoring, an active sonar employing a cross fan beam method is known as a method for obtaining vertical-direction information of a target.
Such an active sonar employing the cross fan beam method which transmits transmitted fan beams, each transmitted fan beam is narrow in a vertical direction (i.e., a water depth direction) and wide in a horizontal direction. And the active sonar transmits a plurality of transmitted fan beams whose elevation angles are mutually different. And receives received fan beams, each being, reversely to the case of the transmitted fan beams, wide in the vertical direction and narrow in the horizontal direction. Such a configuration enables the active sonar to receive an echo sound having significantly high directivity from each of regions that is surrounded by vertical width of the transmitted fan beam and horizontal width of received fan beam. Further, synthesizing of such echo sound signals obtained from the respective regions enables to obtain two-dimensional information. Moreover, the visualization of this two-dimensional information enables the two-dimensional information to be utilized as a two-dimensional image, and also enables to form a three-dimensional shape of the target by using a plurality of such two-dimensional images. It is also known a method that modifies the frequency of the transmitted fan beams together with the elevation angle.
When the above-described cross fan beam method is simply used, there are some cases that accurate information cannot be obtained because of the modulations of the echo sounds or the interfusion of noises. In such background, efforts to eliminate the influence of these unfavorable factors to obtain an accurate information have been made.
For example, in patent document 1 (Japanese Unexamined Patent Application Publication No. 2015-22226), discloses a technique that corrects a frequency changes through the use of a Doppler effect. This technique enables to obtain the position and the velocity of a target by transmitting two kinds of transmitted fan beams, the one whose frequency is increasing along with the changes of its elevation angle, and the other one whose frequency is decreasing along with the changes of its elevation angle. Then the Doppler shift caused by the movement of the target is calculated, and the position and the velocity of the target can be calculated by using the Doppler shift.
Further, in patent document 2 (Japanese Unexamined Patent Application Publication No. 2012-189499) discloses a technique that can eliminate side lobes and multipath signals from a received signal through the use of the phase variance of the received signal. This technique eliminates side lobes and multipath signals by phasing a received signal, and calculating a variance of a received signal, and calculating the correlation between the phase variance and the amplitude of the received signal. Through these processes, the side lobes and the multipath echo sounds are removed. Here, the side lobes means signals that appear, as a result of beam foaming, in a direction different from the direction of a main lobe (a sound axis). And, the multipath signals means signals having arrived at a receiver after passing through a plurality of paths due to reflections or refractions.
The use of the technique disclosed in patent document 1 or patent document 2 enables the acquisition of highly accurate information in relation to a short-distance target. In general cross fan beam methods, however, there has been a problem in that it is difficult to acquire information in relation to a long-distance target. In the case of the long-distance target, a plurality of transmitted fan beams transmitted at mutually different elevation angles propagate through mutually different propagation paths. This is because each transmitted fan beam is reflected or refracted in a deferent way depend on their elevation angle. As a result, for each transmitted fan beam, a difference occurs in the time to arrive at the same target. Further, there is a case where the correspondence between the elevation angles of the transmitted fan beams and the vertical relationship of the irradiation positions at the target arrival becomes unclear. For these reasons, the cross fan beam method has not been used in the acquisition of images of a long-distance target, and the application of the cross fan beam method has been limited to the acquisition of images of a short-distance target. Such a problem of the distance limitation has not been able to be eliminated even through the use of the techniques disclosed in patent documents 1 and 2.
The disclosed subject matter has been made in view of the above problem, and is intended to provide an active sonar capable of acquiring echo sound information in relation to a long-distance target.
In order to solve the above problem, an active sonar according to an aspect of the disclosed subject matter includes a fan beam transmitter, a fan beam receiver, a propagation path calculator, a path time calculator, and a horizontal distance calculator. The fan beam transmitter transmits a plurality of transmitted fan beams whose elevation angles are mutually deferent. Each transmitted fan beam spreading widely in a horizontal direction and spreading narrowly in a vertical direction. The fan beam receiver receives a plurality of received fan beams. Each received fan beam is spreading widely in the vertical direction and spreading narrowly in the horizontal direction. The propagation path calculator calculates a propagation path of each of the transmitted fan beams based on an acoustic velocity profile for a medium and the elevation angle of each of the transmitted fan beams. The path time calculator calculates a path time based on the time point of transmission of a transmitted fan beam and the receiving time point of the echo of the transmitted fan beam. The horizontal distance calculator calculates a horizontal distance from the active sonar to a generation source point of each of echo sounds based on the propagation path and the path time.
As described above, according to an aspect of the disclosed subject matter, an active sonar capable of acquiring echo sound information in relation to a long-distance target is provided.
Exemplary features and advantages of the disclosed subject matter will become apparent from the following detailed description when taken with the accompanying drawings in which:
Hereinafter, example embodiments of the disclosed subject matter will be described in detail referring to the drawings. It should be noted, here, that example embodiments described below contain limitations that are technically preferable in the practice of the disclosed subject matter, but the scope of the disclosed subject matter is not limited to such example embodiments described below. Additionally, note that, in the description below, like constituent elements across figures are denoted by the same reference number, and thereby, the description of one or more of the like constituent elements may be omitted.
The active sonar includes a fan beam transmitter 1, a fan beam receiver 2, a propagation path calculator 3, a path time calculator 4, and a horizontal distance calculator 5.
The fan beam transmitter 1 transmits each of transmitted fan beams spreading widely in a horizontal direction and spreading narrowly in a vertical direction, with the elevation angle of each of the of the transmitted fan beams set to a corresponding one of mutually different elevation angles.
The fan beam receiver 2 receives each of received fan beams spreading widely in the vertical direction and spreading narrowly in the horizontal direction.
The propagation path calculator 3 calculates a propagation path of each of the transmitted fan beams based on the elevation angle of each of the transmitted fan beams and the profile of a medium. A typical medium is water (fresh water or seawater), and the profile in this case includes various kinds of profile elements, such as a water pressure, a water temperature, and a salinity concentration. Further, an acoustic velocity profile regarding an acoustic velocity in a medium can be calculated from these kinds of profile elements. Further, acoustic velocities at each position in the medium is obtained from this acoustic velocity profile, and through the use of the obtained acoustic velocities, the propagation path of each of the transmitted fan beams can be calculated. The calculation of the propagation path can be performed in accordance with, for example, an acoustic ray theory employing the well-known Snell's law.
The path time calculator 4 calculates a path time that is a period of time from the transmission of each of the transmitted fan beams until the reception of each of echo sounds generated due to each of the transmitted fan beams. The path time corresponds to a period of time during which each of the transmitted fan beams and each of the echo sounds, which is included in one of the received fan beams that is associated with each of the transmission fan beams, propagate through the propagation path.
The horizontal distance calculator 5 calculates a horizontal distance from the active sonar to a generation point of each of the echo sounds (i.e., a target) based on the propagation path and the path time, which have been obtained through the two kinds of calculations.
As described above, the configuration according to the present example embodiment enables the acquisition of a horizontal distance from the active sonar to a generation point of each of the echo sounds. The acquisition can be done even though each of the transmitted fan beams and the received fan beam, which is associated with each of the transmitted fan beams, are subjected to reflections and/or refractions.
Thus, the active sonar according to the present example embodiment is capable of acquiring a position of a long-distance target.
Before the description of an active sonar according to this second example embodiment, the outline of the cross fan beam method will be described below.
Moreover, similar signals are acquired with respect to each of the transmitted fan beams 10, which have been transmitted toward the plurality of mutually different orientations, and reception time differences are corrected using transmission time differences to obtain received signals 32. These received signals 32 correspond to one frame (
A period of time from the transmission of a transmitted fan beam until the reception of an echo corresponding to the transmitted fan beam can be converted into a distance from the active sonar to a generation point of each of the echo sounds, through the use of an acoustic velocity in a medium.
Thus, three-dimensional information in relation to a target can be obtained through the use of a plurality of frames having been received at mutually different time points.
In the above description of the cross fan beam method, it is assumed that the upper/lower relations among the plurality of fan beams 10 are not reversed until the plurality of fan beams 10 arrive at the target (the echo sound generation point). This assumption corresponds to short-distance searching in which any reflection does not occur at the surface and the bottom of water. On the other hand, in the case where long-distance searching is performed, the received fan beams may penetrate through paths that are reflected at the surface and/or the bottom of water. Consequently, it is difficult to acquire accurate echo sound information in relation to a target merely by employing the same method as a general cross fan beam method.
The above situation will be described using a specific example.
As shown in
Thus, accurate echo sound information in relation to the short-distance target 300 can be obtained merely by correcting the reception time differences among received signals using transmission time differences.
On the other hand, in the case of
Next, an active sonar according to the present example embodiment will be described.
The fan beam transmitter 210 transmits a transmitted fan beam having a wide spread in the horizontal direction and a narrow spread in the vertical direction. At this time, the fan beam transmitter 210 transmits a plurality of the transmitted fan beams with mutually different elevation angles.
The fan beam receiver 220 receives received fan beams each having been subjected to beam forming so as to have a wide spread in the vertical direction and a narrow spread in the horizontal direction.
The signal processor 230 acquires echo sound information in relation to echo sounds from a target on the basis of received signals included in the received fan beams. In order to achieve this acquisition, the signal processor 230 includes a profile data acquisition unit 231, a propagation path calculation unit 232, a path time calculation unit 233, a horizontal distance calculation unit 234, and a synthesis processing unit 235.
The profile data acquisition unit 231 acquires various kinds of profile data about a medium for fan beams, such as seawater. The profile data includes, for example, water temperature profile data, salinity profile data, and seafloor topography data. These kinds of profile data can be acquired from measurement tools or databases.
The propagation path calculation unit 232 calculates the propagation path of each of the transmitted fan beams and the received fan beams based on the various kinds of profile data having been acquired by the profile data acquisition unit 231 and the elevation angles of the respective transmitted fan beams. The calculation at this time can be performed by, for example, obtaining acoustic velocities at individual positions based on the profile data and applying the acoustic velocities to an acoustic ray theory employing the Snell's law. Here, it is assumed that each of the transmitted fan beams and a received fan beam corresponding to each of the transmitted fan beams propagate through the same propagation path.
The path time calculation unit 233 calculates a path time from a time point at which each of the transmitted fan beams is transmitted until a time point at which each of echo sounds corresponding to each of the transmitted fan beams is received by the receiver.
The horizontal distance calculation unit 234 calculates a horizontal distance from the active sonar 200 to an echo sound generation point of each of the echo sounds based on the propagation path and the path time. A method for the calculation of the horizontal distance will be described later.
The synthesis processing unit 235 synthesizes signals from same horizontal distance and orientation. The horizontal distances are calculated by the horizontal distance calculation unit 234. At this time, the grouping of signals may be performed based on only the horizontal orientations of the signals without taking the vertical orientations of the signals into consideration. This is because it is difficult to acquire detailed information in relation to the vertical orientations of the signals. Since, as having been mentioned in the description of
Next, the method for the calculation of the horizontal distance will be described.
Here, the acoustic velocity c is expressed as a function of a position and is calculated, for each point, from the profile data. The elevation angle θ is obtained from acoustic velocities at both sides of a boundary based on the Snell's law. Further, a distance along an acoustic ray is denoted by s; a horizontal distance is denoted by r; and a period of time is denoted by t.
When a distance (the length of an acoustic ray) by which an acoustic wave proceeds at the i point during a minute time Δti is denoted by Δsi, and a horizontal distance by which the acoustic wave proceeds at this time is denoted by Δri, the following two formulas are satisfied.
Δsi=ci·Δti (1)
Δri=ci·Δti·cos θi (2)
Similarly, for the sound ray at the (i+1) point, the following two formulas are satisfied.
Δsi+1=ci+1·Δti+1 (3)
Δri+1=ci+1·Δti+1·cos θi+1 (4)
Similarly, a horizontal distance corresponding to a distance by which the acoustic wave proceeds at each of points on the propagation path during a minute time corresponding to each of the points can be calculated. Here, a path time, that is, a period of time from the transmission of a certain transmitted fan beam until the reception of a corresponding echo sound, is denoted by T, and a horizontal distance to a target that is an echo sound generation point is denoted by R.
In this case, since the path time is a period of time during which an acoustic wave reciprocates, T is twice the sum of Δt having been described above. Thus, Δti should be added so as to allow the following formula to be satisfied.
A number n, that is, the total number of pieces of data to be added, is obtained from the formula (5). Thus, the horizontal distance R to the target is obtained by summing n pieces of the minute horizontal distance Δri. That is, the horizontal distance R can be represented by the following formula.
Through the above calculations, a horizontal distance to a generation source point of each received signal can be obtained.
Further, an echo sound signal from the position of a specific horizontal orientation and a specific horizontal distance can be obtained by synthesizing the signals from the same orientation and the same distance among the different received fan beams.
Next, the operation of the active sonar will be described.
First, the fan beam transmitter transmits a plurality of the transmitted fan beams with mutually different elevation angles (S1). Next, the fan beam receiver receives echo sound signals through a received fan beam corresponding to each of the transmitted fan beams (S2). Next, the signal processor calculates horizontal distance to the generation source point of each of the echo sound signals. This calculation is repeatedly performed until the completion of the calculations on desired received signals (loop 1, S3 to S5). Further, the signal processor synthesize the echo sound data whose generation source points are located in the same orientations and same horizontal distances (S6).
As described above, the configuration according to the present example embodiment enables the acquisition of a target located at a long-distance position through the use of the cross fan beam method.
In this example embodiment, a specific and detailed configuration example of the active sonar will be described.
The transmission unit 610 includes a wave transmitter 611, a transmitted wave processing unit 612, a DAC 613, and a signal generator 614. Here, the DAC means a digital-to-analog converter.
The wave transmitter 611 converts an electrical signal into an acoustic wave signal, and transmits the acoustic wave signal into the water. The transmitted wave processing unit 612 cuts high frequency noise of a digital signal, and amplifies the electrical signal. The DAC 613 converts the digital signal output from the signal generator 614 into an analog signal. The signal generator 614 generates a digital waveform, such as a PCW (Pulsed Continuous Wave) or an LFM (Linear Frequency Modulation). The wave transmitter 611 transmits a plurality of transmitted fan beams having mutually different elevation angles.
The reception unit 620 includes a wave receiver 621, a received wave processing unit 622, an ADC 623, and a received signal processing unit 624. Here, the ADC means an analog-to-digital converter.
The wave receiver 621 receives an acoustic wave, and converts the received acoustic wave into an electrical signal. The received wave processing unit 622 amplify the electrical signals output from the wave receiver 621 and filtering the signals to eliminate the bands other than a needed band. The ADC 623 converts an analog signal output from the received wave processing unit 622 into a digital signal. The received signal processing unit 624 performs signal processing, such as a baseband conversion process, a phasing process, an AGC (Automatic Gain Control) process, a pulse compression process, and an automatic detection process, on the output of the ADC 623.
The signal processing unit 630 includes an environmental data retention unit 631, an acoustic velocity profile conversion unit 632, a propagation path calculation unit 633, a path time calculation unit 634, a synthesis processing unit 635, a tracking processing unit 636, a horizontal monitoring display 637, and a vertical signal display 638. The environmental data retention unit 631 retains environment data, such as water temperature profile data observed by a water temperature measurement tool or the like through the use of the relations between water depths and water temperatures, salinity profile data observed by a salinity measurement tool or the like, and seafloor topography data. The acoustic velocity profile conversion unit 632 converts a water temperature profile and a salinity profile into an acoustic velocity profile represented by the relations between water depths and acoustic velocities. The propagation path calculator 633 calculates a propagation path based on the acoustic velocity profile output from the acoustic velocity profile conversion unit 632 and the seafloor topography data input from the environmental data. The path time calculation unit 634 calculates a path time of a fan beam based on the elevation angle and the respective propagation path which have been calculated by the propagation path calculation unit 633. The synthesis processing unit 635 calculates a horizontal distance from the active sonar to each of the echo sound generation points based on the calculated propagation path and the path time, and performs synthesis processing on detection results output from the received signal processing unit 624 based on the calculated horizontal distance.
The tracking processing unit 636 associates pieces of data with one another based on the synthesized detection results of the synthesis processing unit 635 and predicts next movement of the target. In the tracking processing, time-series of the predicted positions of a moving target are calculated by means of a state estimation method, such as a Kalman filter or a particle filter, and further, refine the estimated position by taking correlations among a plurality of predicted positions.
The horizontal monitoring display 637 displays the tracking results of the tracking processing unit 636 as a horizontal figure. The vertical signal display 638 vertically displays acoustic data output from the synthesis processing unit 635 by means of a cross fan beam method.
Next, the operation of the active sonar will be described. First, the signal generator 614 forms a digital waveform which is effective for a sonar. A waveform image at this time is assumed to be, for example, a waveform shown in
When a transmission frequency is denoted by Fn; the maximum velocity of an target is denoted by V; and an underwater acoustic velocity is denoted by C, an obtained Doppler shift amount ΔF of the target is represented by the following formula.
In the transmission of the fan beams, a transmission time point of each transmission elements of wave transmission unit 611 (for example, L elements) should be delayed in accordance with the value of each of transmission elevation angles. Here, it is assumed that the number of the transmission elevation angles is two or more.
Further, digital waveforms having been generated by the signal generator 614 are converted by the DAC 613 into L analog signals each associated with a corresponding one of the L elements included in the wave transmitter 611. The L transmitted signals having been converted into the analog signals are each subjected to analog filtering by the transmitted wave processing unit 612 so as to allow only its transmission frequency band (a band resulting from adding a modulation bandwidth to each of the frequencies F1 to Fn) to be passed. Further, the amplitudes of the transmitted signals are amplified by power amplifiers, and then, the transmitted signals are input into the wave transmitter 611. The wave transmitter 611 is constituted by a plurality of wave transmission elements arranged in the vertical direction. Each of the wave transmission elements converts a corresponding one of the electrical signals output from the transmitted wave processing unit 612 into a corresponding one of acoustic waves. Thereafter, each transmission elements transmits a transmitted fan beam with the respective elevation set by vertical beam shifting into the water. Each of the transmitted fan beams having been transmitted into the water propagates toward a target 700.
Meanwhile, reflection waves having been reflected at the target 700 are received and converted into electrical signals by M wave reception elements included in the wave receiver 621. The electrical signals output from the wave receiver 621 are amplified in the received wave processing unit 622. Further, the received wave processing unit 622 performs a filtering process to eliminate bands other than a reception band for the electrical signals. Here, the reception bandwidth is generally wider than a transmission bandwidth. The filtered electrical signals are input into the ADC 623. The ADC 623 converts each of the M analog signals into a corresponding one of digital signals by performing level quantization and temporal sampling on the M analog signals. The digital signals are converted into baseband signals by the received signal processing unit 624 to lower the sampling rate, and then are separated into n bands that are divided in accordance with the elevation angles by filters. The received signal processing unit 624 performs phasing process for the separated M pieces of received data so as to the received data being converted into horizontal beams. Hereinafter, received signals corresponding to n kinds of elevation angles will be referred to as a layer 1, a layer 2, . . . , and a layer n.
When the waveform is any other modulation waveform, such as the LFM, pulse compression processing is performed. The pulse compression process have a procedure of taking the correlation of time point of each element of received signals with the corresponding transmitted signal. Moreover, the absolute level is converted into a corresponding relative level through AGC (Automatic Gain Control) processing, since the absolute levels of each of the received signals greatly varies in accordance with the distance of its propagation path. Through the processes having been described so far, the received signals having been obtained by the wave receiver 621, constituted by the M elements, are converted, for each of n layers of bands, into pieces of position data constituted by an array consisting of rows each associated with a corresponding one of a plurality of horizontal beams (whose number is denoted by N) and columns each associated with a corresponding one of P distances in a distance direction. Here, each of the pieces of position data is denoted by (an orientation number N, a distance number P). That is, n represents the number of vertical planes whose angles toward the vertical direction are made different from one another, and for each of the vertical planes, pieces of received data constituted by grid squares forming a P×N array are obtained. At this time, the total number of all of the pieces of received data is equal to a number obtained by performing a multiplication operation: n (kinds)*P (distances)*N (orientations).
Next, the signal processing unit 630 will be described. The environmental data retention unit 631 retains environmental data that is observed in advance, such as a water temperature profile and a salinity profile, and inputs these kinds of data into the acoustic velocity profile conversion unit 632. A formula for the conversion of an acoustic velocity is associated with a water temperature and a salinity level, and some conversion formulas (for example, an acoustic velocity formula by SOOP) are known. The disclosed subject matter, however, does not depend on the kind of such an acoustic velocity conversion formula. The acoustic velocity profile is input to the propagation path calculator 633. The calculation of the propagation path is performed in accordance with, for example, an acoustic ray theory employing the well-known Snell's law. A propagation path is calculated for each of the layers n. The path time calculator 634 calculates path time sets T1, T2, . . . , and Tn corresponding to the respective n layers based on a transmission time point of each of the transmitted fan beams and a reception time point of each of echo sound signals corresponding to each of the transmitted fan beams, and inputs the results into the synthesis processing unit 635. The synthesis processing unit 635 calculates a horizontal distance to the generation source point of each of the received signals based on the calculation results of the propagation paths and the path times. Further, the synthesis processing unit 635 synthesizes pieces of data included in a corresponding one of the n layers which have the same horizontal direction and the same horizontal distance. Moreover, pieces of data each resulting from this synthesis process are input to the tracking processing unit 636.
As described above, the configuration according to the present example embodiment enables not only the acquisition of echo sound information in relation to a long-distance target through the use of the cross fan beam method, but also the utilization in tracking processing and the like.
In the present example embodiment, processes from the received signal processing unit 624 to the tracking processing unit 636 shown in
Next, the received signal processing unit 624 outputs signals having been subjected to the above-described processing to the synthesis processing unit 635. The synthesis processing unit 635 performs processing for synthesizing groups of pieces of data each of which is included in a corresponding one of the mutually different layers and which have the same horizontal orientation and the same horizontal distance. A specific examples are illustrated in
As described above, this configuration enables the acquisition of echo sound data with a higher certainty, in the case where signals from the same target can be acquired through the deferent fan beams whose elevation angles and propagation paths are mutually deferent. That is, this configuration enables the enhancement of the probability of detecting of a desired target.
As described above, the configuration according to the present example embodiment enables monitoring throughout a wide range spreading from short distance to long distance and spreading from a sea surface to a sea floor, and eliminates or minimizes omissions of detection of a moving target.
A program that allows a computer to execute the above-described processes according to the first to fourth example embodiments and a recording medium that stores in itself the program are also included in the scope of the present invention. For the recording medium, for example, a magnetic disc, a magnetic tape, an optical disc, a magneto-optical disk, or a semiconductor memory device may be used.
The previous description of embodiments is provided to enable a person skilled in the art to make and use the present invention. Moreover, various modifications to these example embodiments will be readily apparent to those skilled in the art, and the generic principles and specific examples defined herein may be applied to other embodiments without the use of inventive faculty. Therefore, the present invention is not intended to be limited to the example embodiments described herein but is to be accorded the widest scope as defined by the limitations of the claims and equivalents.
Further, it is noted that the inventor's intent is to retain all equivalents of the claimed invention even if the claims are amended during prosecution.
Number | Date | Country | Kind |
---|---|---|---|
2016-123200 | Jun 2016 | JP | national |
Number | Name | Date | Kind |
---|---|---|---|
5359334 | Gutman | Oct 1994 | A |
5420827 | Feintuch | May 1995 | A |
6198692 | Sekine | Mar 2001 | B1 |
6445646 | Handa | Sep 2002 | B1 |
6772100 | Lipp | Aug 2004 | B2 |
7362655 | Hester | Apr 2008 | B1 |
20160131760 | Kozuki | May 2016 | A1 |
20170031022 | Ivanov | Feb 2017 | A1 |
20170031023 | Ivanov | Feb 2017 | A1 |
20170350976 | Hauschildt | Dec 2017 | A1 |
Number | Date | Country |
---|---|---|
2012-189499 | Oct 2012 | JP |
2015-222226 | Dec 2015 | JP |
Number | Date | Country | |
---|---|---|---|
20170371037 A1 | Dec 2017 | US |