# Performance Analysis of Target Depth Classification Algorithm Based on Sea Experiment Data

Authors: Anbang Zhao; Xuejie Bi; Nansong Li; Minghui Zhang; Shengchun Piao
Journal: Journal of Sensors
 (2019)

---

## Abstract

Cross-spectrum signals can be calculated by the pressure signals. The sign distribution of cross-spectrum active component can be effectively used for target depth classification algorithm. The algorithm is applicable for depth classification of targets where frequencies can only excite the first two normal modes. The corresponding research results are mainly based on the theoretical study. There are few researches on the algorithm performance based on experiment results. To overcome this research lack, based on the effective depth model, the effects on various receiving depth, source frequency, and received signal-to-noise ratio on the algorithm performance have been studied in this paper. The influence of sound velocity profile parameters (negative gradient, thermocline intensity, thermocline thickness, and up-boundary depth) on the algorithm performance has also been researched. According to the simulation results, proper adjustment of the receiving depths can effectively improve the algorithm performance. The source frequency primarily affects the position of the ideal receiving depth which can be appropriately adjusted according to the depth classification requirements of the real sea environment. The algorithm performance improves gradually with the increase of signal-to-noise ratio. Moreover, the algorithm can also be applied under the conditions of negative gradient and thermocline. The comprehensive sound velocity profile parameters have a large impact on the depth classification performance of the algorithm. Even in the case of strong negative gradient or strong thermocline, the robustness of the algorithm is still high. The feasibility of our presented method has been verified by sea experiment. The practical application value of the ideal receiving depth has been researched and validated. The factors affecting the algorithm performance including line spectrum continuity and received signal-to-noise ratio have also been analyzed in our simulation and real sea experiments.


---

## Body

## 1. Introduction

The movement of underwater platforms in channel encounters the security threats caused by various targets. To achieve the passive localization of the target without the target prior information is very essential. The target depth classification results directly related to the target category in the passive localization process are particularly important. Effective classification of targets (surface targets or underwater targets) has a profound impact on the security and concealment of underwater platforms.The main research focus of this paper is to study a method for target depth classification, which is actually the method for target category classification. In recent years, this area of research has developed huge interests and a large number of research work are published in the field of target depth classification.The initial purpose of the target category is to classify the target depth, rather than accurately estimating the target depth. There are numerous algorithms for depth classification using array signals. Golden [1] performed maximum-likelihood estimation for localization. Liang et al. [2] and Premus et al. [3] proposed mode filtering approaches to acoustic source depth discrimination. Based on either modal decomposition technique [4] or data-based method [5], source depth estimation was realized by Nicolas and Yang. Target depth classification was achieved using a waveguide invariant adaptive matched filter by Goldhahn [6]. There are various researches about target depth estimation on the basis of matched field processing [7–12]. Premus and Backman [13] proposed a matched subspace [14] approach for depth discrimination. The concept of modified modal scintillation index was introduced for source depth classification by An et al. [15]. By using cepstrum techniques, the determination of source depth was carried out by Wu, Zhang et al., and Mitchell and Bedford [16–18]. Source range and depth are obtained using a ray-based backpropagation algorithm by Felisberto et al. [19]. These depth classification algorithms have high accuracy, but most of these methods are not feasible to be applied in real application due to their high complexity and time consumption. In addition, considering the space limitations and security requirements of underwater platforms, it may be impossible to place or tow an array with a large aperture. However, these research hold high theoretical value but need to be further researched for implying in real application scenario for target depth classification by placing one or two sensors.The cross-spectrum calculation is performed using signals acquired by a single 3D vector sensor [20, 21] or two 2D vector sensors [22, 23], and the target depth classification can be achieved based on the characteristics of the cross-spectrum signals. However, most of the above algorithms are theoretical studies [20–23], and the actual performance of the algorithm is rarely analyzed through the analysis of experiment data processing results [24, 25]. Based on the sea experiment data processing results, this paper analyzes the influence of the ideal receiving depth concept on the performance of the algorithm while verifying the feasibility of the target depth classification algorithm. It shows the real application value of the ideal receiving depth concept. Moreover, the paper makes an in-depth research about the affection of the sound velocity profile parameters (negative gradient, thermocline intensity, thermocline thickness, up-boundary depth, etc.) on the algorithm performance.


## 2. Theoretical Model

The acoustic propagation characteristics of shallow water field on the basis of the two-layer fluid-fluid model is shown in Figure1. H is the sea depth; seawater density ρ1 is constant, and the sound velocity of seawater c1 is constant or variable with depth z only. The sediment layer is homogeneous fluid half-space whose sound velocity c2 and density ρ2 are constant. A point source of unit amplitude and angular frequency ω is located at (0,z0) in cylindrical coordinates. The receiver is located at (r,z). The sea surface is an absolutely soft interface whose depth is z=0 [26].Figure 1
Two-layer fluid-fluid model.Ifc1 is constant, this model is known as the Pekeris model. A Pekeris channel of depth H in the shallow water is essentially the same as a waveguide with two pressure-release boundaries whose depth is He [27]. He is called the effective depth.
(1)He=H1+1bk1Hsinαc,where b=ρ1/ρ2; k1=ω/c1; αc=cos−1c1/c2 is the critical grazing angle of the bottom interface. Thus, the pressure field can be expressed as follows:
(2)Pz0,z,r=jHe∑n12πξnrψnz0Ψnze−jξnr−π/4,where n is the serial number of the normal mode; ξn=k12−βn is the nth order eigenvalue; ξn and βn are also called horizontal and vertical wave number, respectively. Ψnz=sinβnz is the mode depth function. In the effective depth model, βn=nγ, γ=π/He.Pz0,z1,r and Pz0,z2,r are signals collected by receiving pressure sensors at different depths z1 and z2. The active component of their cross-spectrum signal is defined as IpA:
(3)IpA=RePz0,z1,rP∗z0,z2,r=1r∑nA1nA2n+∑n,n≠m∑mA1nA2mcosΔξmnr.In the equation,Ainz0,zi=1/He2πξnΨnz0Ψnzi,i=1,2, the superscript ∗ represents the complex conjugate operation; Re⋅ represents the operation of getting real part from the complex number.The cutoff frequency of thenth order normal mode can be calculated as
(4)fn=2n−1c1c24Hc22−c12.If source frequencyf satisfies fn<f<fn+1, f can excite the first nth order normal modes.IfIpA=1/rSrz+Sz, Srz and Sz are two components of IpA:
(5)Srz=∑n,n≠m∑mA1nA2mcosΔξmnr,(6)Sz=∑nA1nA2n.Srz represents the component of IpA whose sign changes in the horizontal and vertical direction. Sz is denoted for the component of IpA whose sign changes in the vertical direction only. According to equation (5), the peak value of Srz is ∑n,n≠m∑mA1nA2m. Srz resembles cosine signal but its peak value varies with source and receiver depth.The frequency can only excite the first two normal modes, and target depth classification can be realized using sign distribution ofIpA. Based on effective depth model, when source frequency can only excite the first two modes, there are ideal receiving depths z1+z2=He resulting in A11A22+A12A21cosΔξ12r≈0, so that Srz is tending to zero. The impact of Srz on the sign distribution of IpA can be ignored. If z1+z2=He, it can guarantee that the sign distribution of IpA will not change in the horizontal direction but change in the vertical direction only.The critical depthhs is defined as the position of boundary of IpA sign when z1+z2=He. All the signs of IpA above hs are identical, but different from the value lesser than hs. The approximate expression of hs is
(7)hs≈1γcos−1ξ216ξ1cos2γz1.The flow chart of target depth classification is presented in Figure2. According to the simulation parameters used in this paper, when the source depth is lesser than hs, the sign of IpA is negative, and sources located at these positions can be identified as surface targets. When the source depth is greater than hs, the sign of IpA is positive, and the source at these positions can be identified as underwater targets. Corresponding to different simulation parameters, the sign of surface targets’ IpA can be positive and the sign of underwater targets’ IpA can be negative.Figure 2
Flow chart of target depth classification.Based on the above theory which have already been published in [22, 23], we can conclude that the value of critical depth changes with different receiving depth, source frequency, sound velocity profile, and seabed. The reason why these factors affect algorithm performance has never been researched. We have to develop a model of these parameters affecting the algorithm effectiveness. The selection method of receiving depths has also been described briefly.Only ifξ2/16ξ1cos2γz1≤1, the existence of critical depth hs can be guaranteed. Thus, the maximum allowable value of the first sensor depth is
(8)z1max=1γcos−1ξ216ξ1.⋅ represents that fractions are rounded down.Ifz1+z2≠He, the sign of Srz changes in the horizontal direction whose variation law depends on cosΔξ12r component. The sign of cosΔξ12r will hop every other π/Δξ12 in the horizontal direction. The sign of Srz changes in the vertical direction whose variation law depends on A11A22+A12A21 component. The sign of A11A22+A12A21 will hop at about He/2. It can be found that whether in the horizontal or vertical direction, the sign hopping of Srz is both not connected with the value of z1,z2, but the amplitude of Srz varies with z1,z2. So, the sign of IpA will vary with z1,z2. At most depths, Srz≠0, so that Srz will exert significant influence on the sign distribution of IpA. On the premise of the existence of the critical depth, the existence of Srz will make the boundary of IpA sign change from straight line to curve that resembles cosine distribution whose peak value changes with z1,z2.Summarizing our presented algorithm, the acoustic field interference structureIpA changes with the variation of receiving depths, source frequency, sound velocity profile, and seabed, which alters the value of hs and the amplitude of Srz accordingly. The change of hs and Srz leads to fluctuate the algorithm performance. The ideal receiving depths corresponding to various conditions are different. So the performance can be improved through proper adjustment of the receiving depths according to the change of these key factors.


## 3. Simulation Results

To verify the feasibility of our presented technique, we conducted a simulation experiment. Sea depth is selected asH=22m; densities of seawater and sediment are ρ1=1.0g/cm3,ρ2=1.75g/cm3; sound velocity in the sediment is c2=1750m/s; source frequency is f=120Hz, and the range of source depth is 1~22 m. The range of horizontal distance is 2~8 km. f2=100Hz,f3=165Hz; thus, the source frequency can only excite the first two modes. The impact of sound velocity profile on the algorithm performance has been analyzed considering three common types of sound velocity profile (isovelocity, negative gradient, and thermocline) in shallow water. Corresponding to isovelocity profile, sound velocity in the seawater is c1=1500m/s.With respect to different negative gradient profiles, sound velocity in the seawater isc1=1540m/s at z=0m and c1=1520,1500,1480m/s at z=22m; the gradient is negative and constant, namely, the gradient of different conditions is Gc=0.9,1.8,2.7s−1.Thermocline profiles are effected mainly due to thermocline intensity, thermocline thickness, and up-boundary depth. The sound velocity in the seawater isc1=1540m/s at z=0m, and c1=1530,1520,1510,1500m/s at z=22m corresponding to different conditions. Thus, the thermocline intensity can be Gc=2.5,5s−1, the thermocline thickness may be considered as 4m,8m,12m, and the up-boundary depth may be 5m,10m.


### 3.1. Impact of Receiving Depth on Algorithm Performance

Under the condition of isovelocity profile, according to equation (1), He≈29m. According to equation (8), only if z1≤12m, target depth classification can be realized. The relationship between Srz,hs and z1,z2 is shown in Figure 3. It can be witnessed from Figure 3 that with the increase of z1, the peak value of Srz gradually increases. The peak value of Srz will also vary with the distance between z1+z2 and He, so as to affect the sign distribution of IpA; hs will reduce with the increase of z1. Consequently, in order to get a better classification of target depth, the value of z1,z2 should be considered comprehensively. Firstly, the ideal range of z1,z2 should be obtained. Then, conferring to the requirement of hs, the ideal range of z1 should be confirmed. Finally, in conditions permitting, the value of z1 should be reduced appropriately to improve the performance of the depth classification algorithm.Figure 3
Relationship betweenSrz,hs and z1,z2.Figure4 illustrates the impact of z1 on the sign distribution of Sz and IpA. In this article, all sign distributions are indicated by white and grey as positive sign and black indicates negative sign. It can be observed from Figure 4 that the critical depth hs is inversely proportional to z1. When z1≤12m, the target depth classification can be achieved by using the sign distribution of IpA. If z1>12m, the sign of IpA does not change in the horizontal direction, but the signs are all positive. The sign boundary does not exist and the target depth classification cannot be performed.Figure 4
Impact ofz1 on the sign distribution of Sz and IpA when z1+z2=28m. (a) Relationship between Sz and z1. (b) Relationship between IpA and z1.
(a)
(b)Figure5 expresses the impact of z1+z2 on the sign distribution of Sz and IpA when z1=11m. As it can be witnessed from Figure 5, when z1+z2=26,27m, the sign of IpA changes with horizontal distance r, which coincides with the change rule shown in Figure 3. With the change of the distance between z1+z2 and He, Srz does not tend to 0. Srz affects the sign distribution of IpA far more than Sz resulting in the sign distributions as shown in Figure 5 when z1+z2=26,27m. If z1+z2=28,29m, Srz tends to 0, the sign of IpA hardly changes with r. Figure 6 is the sign distribution of Srz under the condition of isovelocity profile. The position of the sign hopping in the vertical direction is about 14 m, roughly equal to He/2.Figure 5
Impact ofz1+z2 on the sign distribution of Sz and IpA when z1=11m. (a) Relationship between Sz and z1+z2. (b) Relationship between IpA and z1+z2.
(a)
(b)Figure 6
The sign distribution ofSrz.


### 3.2. Impact of Source Frequency on Algorithm Performance

When the source frequency can only excite the first two modes, the change of the source frequency will affect the target depth classification effect. It can be known from equation (1) that, in case of isovelocity profile, the effective depth He is inversely proportional to the frequency. The maximum depth of the first receiving sensor (z1max) is related to He,ξ1,ξ2, and the change of the frequency causes the change of the above three variables. Therefore, there is no obvious changing law of z1max along with frequency. Similarly, the critical depth hs is also related to He,ξ1,ξ2, so the relationship between hs and frequency does not follow any obvious regular pattern. When f∈100,165, He∈26.9,30.2, the change scope of He is small. hs and z1max are both directly related to ξ2/ξ1. When f∈100,165, ξ2/ξ1∈0.8875,0.9520, the change scope of ξ2/ξ1 is also not large. Therefore, as the frequency changes, hs and z1max will vary, but the impact of frequency on the performance of the algorithm is insignificant. When f=110,150Hz, He=28m, the sign distributions of IpA are graphed in Figure 7. Through the analysis of Figure 7, the influence of source frequency on the algorithm performance has been analyzed.Figure 7
The impact of frequency on the sign distribution ofIpA under condition of isovelocity profile. (a) f=110Hz. (b) f=150Hz.
(a)
(b)It can be witnessed from Figure7 that the change of the frequency affects the value of hs and z1max and the sign distribution of IpA changes too. However, by adjusting the receiving depth appropriately, based on the sign distribution of IpA, the depth classification algorithm can still obtain ideal classification results under different frequency conditions. Therefore, the change in frequency has minor impact on the performance of the algorithm.


### 3.3. Impact of Signal-Noise Ratio on Algorithm Performance

In this section, we will discuss the effect of SNR on the performance of our presented algorithm. The research object changes from point source to radiation noise of surface or underwater targets. The parameters considered are as follows: heading angle isψ=50∘, target velocity is vt=11m/s, closest distance is rmin=1200m, initial distance is r0=2000m, platform velocity is vs=2m/s, and range of SNR (signal-noise ratio) is SNR=−30~0dB. Assuming that the target radiated noise is a single-frequency line spectrum superimposed on the continuum, two pressure sensors are used to receive signals at depths of 45 m and 90 m, respectively, and the total sailing time is 400 s. The target is kept in motion away from the receiving platform. With respect to the abovementioned simulation parameters, the influence of SNR on the performance of the algorithm is analyzed. Figure 8 expresses the relationship between the performance of the target depth classification algorithm and SNR under isovelocity condition. In the process of target depth classification, the index that implies the performance of the algorithm is the probability of identified as underwater target in percentage. The probability of identified as underwater target is defined as the proportion of the cases in which the sign of the pressure cross-spectrum active component is positive. It is equal to ∑i=1N1Li/M1N1×100%. N1 is the total number of times of Monte Carlo simulation; M1 is the total sampling numbers in each Monte Carlo simulation; Li (i=1,2,⋯,N1) is the number of times of the cases in which the sign is positive when discriminating the sign of the pressure cross-spectrum active component in each simulation, 0≤Li≤M1. It can be found from Figure 8 that when SNR≥−20dB, the target depth classification effect is good, and there is an obvious sign boundary at about 7 m. When SNR=−30dB, the algorithm can achieve the target depth classification, but the classification effect is poor.Figure 8
The impact of SNR on the algorithm performance under condition of isovelocity profile.


### 3.4. Impact of Sound Velocity Profile on Algorithm Performance

#### 3.4.1. Negative Gradient Sound Velocity Profile

Receiving depths are considered asz1=11m,z2=17m. The negative gradient sound velocity profiles corresponding to different negative gradient Gc=0.9,1.8,2.7s−1 are shown in Figure 9(a). Under different conditions of Gc=0.9,1.8,2.7s−1, the sign distributions of IpA are shown in Figure 9(b). It can be witnessed from Figure 9 that the sound velocity gradient is inversely proportional to the depth classification effect and is directly proportional to the critical depth. Even in the case of strong negative gradient, the robustness of the algorithm is high.Figure 9
The impact of gradient on the algorithm performance under condition of negative gradient sound velocity profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)


#### 3.4.2. Thermocline Existing in Sound Velocity Profile

When the thermocline thickness is fixed to8m and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA under different conditions of Gc=2.5,5s−1 are displayed in Figure 10. When the thermocline intensity is fixed to Gc=2.5s−1 and the thermocline thickness is fixed to 8m, the sound velocity profiles and the sign distributions of IpA when the up-boundary depths are 5m,10m are shown in Figure 10.Figure 10
The impact of thermocline intensity and up-boundary depth on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be witnessed from Figure10 that the algorithm in this paper is also applicable when thermocline exists in the sound velocity profile. The depth classification effect of the algorithm deteriorates as thermocline intensity increases, and the critical depth value is nearly constant. While the depth classification effect of the algorithm deteriorates as the depth of the up-boundary increases, the critical depth value almost keeps the same. In the presence of thermocline, both the thermocline intensity and up-boundary depth have a large impact on the depth classification performance of the algorithm, but the critical depth value remains almost the same. Even with the presence of a strong thermocline, the depth classification performance of the algorithm is good.When the thermocline intensity is fixed toGc=2.5s−1 and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA when the thermocline thicknesses are 4m,8m,12m are displayed in Figure 11.Figure 11
The impact of thermocline thickness on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be seen from Figure11 that the depth classification effect of the algorithm deteriorates with the increase of the thermocline thickness and the critical depth value increases with the increase of the thermocline thickness. Therefore, the thermocline thickness has a great impact on the depth classification performance of the algorithm.


## 3.1. Impact of Receiving Depth on Algorithm Performance

Under the condition of isovelocity profile, according to equation (1), He≈29m. According to equation (8), only if z1≤12m, target depth classification can be realized. The relationship between Srz,hs and z1,z2 is shown in Figure 3. It can be witnessed from Figure 3 that with the increase of z1, the peak value of Srz gradually increases. The peak value of Srz will also vary with the distance between z1+z2 and He, so as to affect the sign distribution of IpA; hs will reduce with the increase of z1. Consequently, in order to get a better classification of target depth, the value of z1,z2 should be considered comprehensively. Firstly, the ideal range of z1,z2 should be obtained. Then, conferring to the requirement of hs, the ideal range of z1 should be confirmed. Finally, in conditions permitting, the value of z1 should be reduced appropriately to improve the performance of the depth classification algorithm.Figure 3
Relationship betweenSrz,hs and z1,z2.Figure4 illustrates the impact of z1 on the sign distribution of Sz and IpA. In this article, all sign distributions are indicated by white and grey as positive sign and black indicates negative sign. It can be observed from Figure 4 that the critical depth hs is inversely proportional to z1. When z1≤12m, the target depth classification can be achieved by using the sign distribution of IpA. If z1>12m, the sign of IpA does not change in the horizontal direction, but the signs are all positive. The sign boundary does not exist and the target depth classification cannot be performed.Figure 4
Impact ofz1 on the sign distribution of Sz and IpA when z1+z2=28m. (a) Relationship between Sz and z1. (b) Relationship between IpA and z1.
(a)
(b)Figure5 expresses the impact of z1+z2 on the sign distribution of Sz and IpA when z1=11m. As it can be witnessed from Figure 5, when z1+z2=26,27m, the sign of IpA changes with horizontal distance r, which coincides with the change rule shown in Figure 3. With the change of the distance between z1+z2 and He, Srz does not tend to 0. Srz affects the sign distribution of IpA far more than Sz resulting in the sign distributions as shown in Figure 5 when z1+z2=26,27m. If z1+z2=28,29m, Srz tends to 0, the sign of IpA hardly changes with r. Figure 6 is the sign distribution of Srz under the condition of isovelocity profile. The position of the sign hopping in the vertical direction is about 14 m, roughly equal to He/2.Figure 5
Impact ofz1+z2 on the sign distribution of Sz and IpA when z1=11m. (a) Relationship between Sz and z1+z2. (b) Relationship between IpA and z1+z2.
(a)
(b)Figure 6
The sign distribution ofSrz.


## 3.2. Impact of Source Frequency on Algorithm Performance

When the source frequency can only excite the first two modes, the change of the source frequency will affect the target depth classification effect. It can be known from equation (1) that, in case of isovelocity profile, the effective depth He is inversely proportional to the frequency. The maximum depth of the first receiving sensor (z1max) is related to He,ξ1,ξ2, and the change of the frequency causes the change of the above three variables. Therefore, there is no obvious changing law of z1max along with frequency. Similarly, the critical depth hs is also related to He,ξ1,ξ2, so the relationship between hs and frequency does not follow any obvious regular pattern. When f∈100,165, He∈26.9,30.2, the change scope of He is small. hs and z1max are both directly related to ξ2/ξ1. When f∈100,165, ξ2/ξ1∈0.8875,0.9520, the change scope of ξ2/ξ1 is also not large. Therefore, as the frequency changes, hs and z1max will vary, but the impact of frequency on the performance of the algorithm is insignificant. When f=110,150Hz, He=28m, the sign distributions of IpA are graphed in Figure 7. Through the analysis of Figure 7, the influence of source frequency on the algorithm performance has been analyzed.Figure 7
The impact of frequency on the sign distribution ofIpA under condition of isovelocity profile. (a) f=110Hz. (b) f=150Hz.
(a)
(b)It can be witnessed from Figure7 that the change of the frequency affects the value of hs and z1max and the sign distribution of IpA changes too. However, by adjusting the receiving depth appropriately, based on the sign distribution of IpA, the depth classification algorithm can still obtain ideal classification results under different frequency conditions. Therefore, the change in frequency has minor impact on the performance of the algorithm.


## 3.3. Impact of Signal-Noise Ratio on Algorithm Performance

In this section, we will discuss the effect of SNR on the performance of our presented algorithm. The research object changes from point source to radiation noise of surface or underwater targets. The parameters considered are as follows: heading angle isψ=50∘, target velocity is vt=11m/s, closest distance is rmin=1200m, initial distance is r0=2000m, platform velocity is vs=2m/s, and range of SNR (signal-noise ratio) is SNR=−30~0dB. Assuming that the target radiated noise is a single-frequency line spectrum superimposed on the continuum, two pressure sensors are used to receive signals at depths of 45 m and 90 m, respectively, and the total sailing time is 400 s. The target is kept in motion away from the receiving platform. With respect to the abovementioned simulation parameters, the influence of SNR on the performance of the algorithm is analyzed. Figure 8 expresses the relationship between the performance of the target depth classification algorithm and SNR under isovelocity condition. In the process of target depth classification, the index that implies the performance of the algorithm is the probability of identified as underwater target in percentage. The probability of identified as underwater target is defined as the proportion of the cases in which the sign of the pressure cross-spectrum active component is positive. It is equal to ∑i=1N1Li/M1N1×100%. N1 is the total number of times of Monte Carlo simulation; M1 is the total sampling numbers in each Monte Carlo simulation; Li (i=1,2,⋯,N1) is the number of times of the cases in which the sign is positive when discriminating the sign of the pressure cross-spectrum active component in each simulation, 0≤Li≤M1. It can be found from Figure 8 that when SNR≥−20dB, the target depth classification effect is good, and there is an obvious sign boundary at about 7 m. When SNR=−30dB, the algorithm can achieve the target depth classification, but the classification effect is poor.Figure 8
The impact of SNR on the algorithm performance under condition of isovelocity profile.


## 3.4. Impact of Sound Velocity Profile on Algorithm Performance

### 3.4.1. Negative Gradient Sound Velocity Profile

Receiving depths are considered asz1=11m,z2=17m. The negative gradient sound velocity profiles corresponding to different negative gradient Gc=0.9,1.8,2.7s−1 are shown in Figure 9(a). Under different conditions of Gc=0.9,1.8,2.7s−1, the sign distributions of IpA are shown in Figure 9(b). It can be witnessed from Figure 9 that the sound velocity gradient is inversely proportional to the depth classification effect and is directly proportional to the critical depth. Even in the case of strong negative gradient, the robustness of the algorithm is high.Figure 9
The impact of gradient on the algorithm performance under condition of negative gradient sound velocity profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)


### 3.4.2. Thermocline Existing in Sound Velocity Profile

When the thermocline thickness is fixed to8m and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA under different conditions of Gc=2.5,5s−1 are displayed in Figure 10. When the thermocline intensity is fixed to Gc=2.5s−1 and the thermocline thickness is fixed to 8m, the sound velocity profiles and the sign distributions of IpA when the up-boundary depths are 5m,10m are shown in Figure 10.Figure 10
The impact of thermocline intensity and up-boundary depth on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be witnessed from Figure10 that the algorithm in this paper is also applicable when thermocline exists in the sound velocity profile. The depth classification effect of the algorithm deteriorates as thermocline intensity increases, and the critical depth value is nearly constant. While the depth classification effect of the algorithm deteriorates as the depth of the up-boundary increases, the critical depth value almost keeps the same. In the presence of thermocline, both the thermocline intensity and up-boundary depth have a large impact on the depth classification performance of the algorithm, but the critical depth value remains almost the same. Even with the presence of a strong thermocline, the depth classification performance of the algorithm is good.When the thermocline intensity is fixed toGc=2.5s−1 and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA when the thermocline thicknesses are 4m,8m,12m are displayed in Figure 11.Figure 11
The impact of thermocline thickness on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be seen from Figure11 that the depth classification effect of the algorithm deteriorates with the increase of the thermocline thickness and the critical depth value increases with the increase of the thermocline thickness. Therefore, the thermocline thickness has a great impact on the depth classification performance of the algorithm.


## 3.4.1. Negative Gradient Sound Velocity Profile

Receiving depths are considered asz1=11m,z2=17m. The negative gradient sound velocity profiles corresponding to different negative gradient Gc=0.9,1.8,2.7s−1 are shown in Figure 9(a). Under different conditions of Gc=0.9,1.8,2.7s−1, the sign distributions of IpA are shown in Figure 9(b). It can be witnessed from Figure 9 that the sound velocity gradient is inversely proportional to the depth classification effect and is directly proportional to the critical depth. Even in the case of strong negative gradient, the robustness of the algorithm is high.Figure 9
The impact of gradient on the algorithm performance under condition of negative gradient sound velocity profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)


## 3.4.2. Thermocline Existing in Sound Velocity Profile

When the thermocline thickness is fixed to8m and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA under different conditions of Gc=2.5,5s−1 are displayed in Figure 10. When the thermocline intensity is fixed to Gc=2.5s−1 and the thermocline thickness is fixed to 8m, the sound velocity profiles and the sign distributions of IpA when the up-boundary depths are 5m,10m are shown in Figure 10.Figure 10
The impact of thermocline intensity and up-boundary depth on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be witnessed from Figure10 that the algorithm in this paper is also applicable when thermocline exists in the sound velocity profile. The depth classification effect of the algorithm deteriorates as thermocline intensity increases, and the critical depth value is nearly constant. While the depth classification effect of the algorithm deteriorates as the depth of the up-boundary increases, the critical depth value almost keeps the same. In the presence of thermocline, both the thermocline intensity and up-boundary depth have a large impact on the depth classification performance of the algorithm, but the critical depth value remains almost the same. Even with the presence of a strong thermocline, the depth classification performance of the algorithm is good.When the thermocline intensity is fixed toGc=2.5s−1 and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of IpA when the thermocline thicknesses are 4m,8m,12m are displayed in Figure 11.Figure 11
The impact of thermocline thickness on the algorithm performance under condition of thermocline existing in profile. (a) Sound velocity profile. (b) The sign distribution ofIpA.
(a)
(b)It can be seen from Figure11 that the depth classification effect of the algorithm deteriorates with the increase of the thermocline thickness and the critical depth value increases with the increase of the thermocline thickness. Therefore, the thermocline thickness has a great impact on the depth classification performance of the algorithm.


## 4. Experiment Results

In order to examine the influence of the ideal receiving depth concept on the performance of target depth classification algorithm, an experiment was carried out in Yellow Sea in Qingdao, China. The layout of the sea trial is displayed in Figure12. The source is suspended into the water by the sending ship at about 8 m. The sound velocity profile near the source position is calculated by CTD (Conductivity-Temperature-Depth System). The CTD device is shown in Figure 13(a), and the collected sound velocity profile is graphed in Figure 14. It can be witnessed from Figure 14 that the type of the sound velocity profile is negative gradient sound velocity profile.Figure 12
Sea experiment layout diagram.Figure 13
Devices used in the sea experiment. (a) CTD. (b) Receiving sensor and TD.
(a)
(b)Figure 14
Sound velocity profile.Four receiving pressure sensors are equally spaced into the water by the receiving ship, and one TD (depth sensor) is fixed at each sensor. The signal collection device is fixed between the second and third sensors, and all sensors are connected to the signal collection device. The effect of the ocean current on the experiment results is reduced by fixing the weight on the rope. The receiving device is shown in Figure13(b).According to the CTD, the sea depth is approximately 22 m. It can be obtained from equation (1) that He∈29,30m. According to equation (8), z1max<13m. The length of the pressure sensor cable is limited; in order to ensure the stability of the sensor power supply, the overall length of pressure array length cannot exceed 6 m. The four array elements are laid at equal intervals to keep the interval between two neighboring elements of 2 m. The designed depths of the four array elements are 11, 13, 15, and 17 m. According to the TD record, the actual depths of the four elements are z1,z2,z3,z4≈10.8,12.5,14.6,16.5m. Due to the existence of the sea current in the sea, the four elements are deviated slightly from their designed depth. When the sending ship is about 3, 5, and 7 km away from the receiving ship, the single-frequency continuous sine transmitted by the source is repeated for two times where its frequency is 120 Hz and its time length is 1 min. The serial number of the corresponding signals is signal 1 to signal 6. As there are a total of four array elements, during the processing of the transmitted signal, signals collected by two elements are selected to perform the target depth classification using the algorithm described in this paper. There are six different combinations of elements. Different combinations of elements can result in different target depth classification results. A sort is defined for the six element combinations, and the serial number of element combinations is 1 to 6, corresponding to z1,z2, z1,z3, z1,z4, z2,z3, z2,z4, and z3,z4, respectively.If the parameters of the negative gradient sound velocity profile in the experiment are approximated asc1=1540m/s at z=0m and c1=1520m/s at z=22m, the simulation analyzes the effect of receiving depths on depth classification performance under negative gradient conditions, as shown in Figure 15.Figure 15
Relationship between the effect of depth classification and receiving depths under the condition of negative gradient profile. (a) First element combination. (b) Second element combination. (c) Third element combination. (d) Fourth element combination. (e) Fifth element combination. (f) Sixth element combination.
(a)
(b)
(c)
(d)
(e)
(f)The algorithm presented in this paper is applicable to the isovelocity condition as well as to the negative gradient condition. Based on the signals collected by the third and fifth element combinations, the sign distributions of the calculatedIpA can be used for depth classification of both surface targets and underwater targets. The other four combinations can only accurately classify the depth of the underwater targets but cannot accurately classify the depth of the surface targets. In addition, although the critical depth of the fifth combination is smaller than that of the third combination, the sign boundary of the third combination has a smaller fluctuation than the fifth combination. If the above two combinations can both meet the actual depth classification requirements, it can be considered that the depth classification effect of the third combination is better than that of the fifth combination. The change of the sound velocity profile type will have a certain impact on the depth classification performance, but the impact of the receiving depths on the depth classification performance is similar to that of the isovelocity condition. An improved depth classification effect can be obtained by appropriately adjusting the depths of the two receivers. In addition, according to the simulation results, the source used in this experiment should theoretically be identified as an underwater target.The PSD has been achieved through smoothed FFT spectrum calculation. The Hamming window is the corresponding weighted window used for smoothing. Its window length and step are 4 s and 2 s, respectively.The spectra of the signal 1 to signal 6 received by the four elements in the actual sea experiment are graphed in Figure16. It can be witnessed from Figure 16 that when the horizontal distance between the source and the receiver is r≈3km, each spectrum of signal 1 received by four elements has a clear line spectrum at f=120Hz, but only the spectra of signal 2 received by the first three elements have a clear line spectrum at f=120Hz. There is not a clear line spectrum at f=120Hz in the signal 2 bandwidth received by the fourth element. When r≈5km, only the spectra of the signals 3 and 4 received by the first element have a clear line spectrum at f=120Hz. The spectra of the signals 3 and 4 received by the other three elements do not have a clear line spectrum at f=120Hz. When r≈7km, the spectra of signals 5 and 6 received by all the four elements do not have a clear line spectrum at f=120Hz.Figure 16
The frequency spectra of the signals collected by four elements. (a) Signal 1 (r≈3km). (b) Signal 2 (r≈3km). (c) Signal 3 (r≈5km). (d) Signal 4 (r≈5km). (e) Signal 5 (r≈7km). (f) Signal 6 (r≈7km).
(a)
(b)
(c)
(d)
(e)
(f)Line spectrum intensity and received SNR have a large impact on the depth classification effect. The target depth classification is performed using the signals 1 to 6, and the probability that the source is accurately discriminated as an underwater target is shown in Table1. The serial number in the table is the serial number of element combinations.Table 1
Depth classification accuracy of underwater target.
Combination number
1
2
3
4
5
6
Signal 1
99.1%
96.5%
93.8%
89.4%
88.5%
90.3%
Signal 2
86.7%
74.3%
77.9%
69.9%
73.5%
75.2%
Signal 3
57.5%
66.4%
33.6%
57.5%
58.4%
48.7%
Signal 4
60.2%
61.1%
58.4%
61.1%
46.0%
53.1%
Signal 5
77.0%
56.6%
56.6%
41.6%
68.1%
43.4%
Signal 6
73.5%
69.9%
45.1%
62.0%
54.9%
49.6%According to the data in Table1, when r≈3km, the source can be accurately identified as an underwater target using signals collected by each of the six array element combinations. This outcome is identical to the simulation results. When r≈5,7km, signals collected by only some of the element combinations are useful for target depth classification, and the corresponding classification effects are poor. This is mainly because the source level of the low-frequency source is limited, and the amplification factor of the power amplifier is also limited. So, when the transmitted signal propagates to a distant place the signal becomes weak. The low SNR affects the performance of the target depth classification algorithm. In addition, when r≈3km, the obtained depth classification results cannot prove the feasibility of the method. This method accurately classifies underwater targets and perfectly classifies surface targets at the same time. In this experiment, only one source depth was tested.In the actual experiment data processing, there are surface targets such as fishing boats and patrol boats in the sea area. On the contrary, the source carried by the sending ship is the only underwater target which can radiate strong line spectrum signal. Therefore, by performing spectrum analysis on the received signal, if there are clear continuous line spectrum outside the time range of the underwater target radiated signals, they are the spectrum of the signals radiated by surface targets. After time-frequency analysis of the signal in the whole observation time range, it is found that there are evident continuous line spectrum corresponding to the four different segments of received signal as it can be witnessed in Figure17.Figure 17
Time-frequency distributions of surface target radiated signals. (a) First segment of signal. (b) Second segment of signal. (c) Third segment of signal. (d) Fourth segment of signal.
(a)
(b)
(c)
(d)The first segment of signal contains two line spectra whose frequencies aref1,f2=106,121.5Hz. The line spectrum frequency of the second signal is f3=142.5Hz. There are two line spectra in the third signal, whose frequencies are f4,f5=119.8,159.5Hz. f6,f7,f8,f9=123,138.5,153.5,169Hz are line spectrum frequencies of the fourth signal. Each segment of signal containing multiple clear continuous line spectrum may be harmonic signals radiated by the same surface target, but it may also be single or harmonic signals radiated by different surface targets. Therefore, the order of these line spectra is represented by the line spectrum number f1~f9. The corresponding depth classification results are shown in Table 2. The data in the table is the probability that the target is accurately classified as the surface target. The serial number in Table 2 is the serial number of the element combinations.Table 2
Depth classification accuracy of surface target.
Combination number
1
2
3
4
5
6
f1
13.8%
23.5%
77.8%
14.8%
56.8%
25.6%
f2
12.3%
21.8%
84.3%
14.7%
68.9%
25.2%
f3
15.6%
49.0%
46.8%
41.4%
38.7%
40.3%
f4
2.6%
29.8%
87.5%
15.3%
78.2%
40.5%
f5
5.1%
89.8%
97.5%
20.1%
86.1%
7.7%
f6
7.1%
15.9%
81.4%
7.1%
67.3%
30.1%
f7
16.8%
42.5%
77.0%
21.2%
69.9%
38.9%
f8
11.5%
86.7%
94.7%
11.5%
76.1%
14.2%
f9
8.0%
99.1%
97.4%
34.5%
95.6%
0.9%According to the data in Table2, except for the third line spectrum, only using signals collected by the third and fifth elements can accurately identify the target as a surface target. It is only based on signals collected by these two combinations, whether it is a surface target or an underwater target, its depth can be accurately classified. Among the six combinations, the depth summation of these two combinations is closest to 29 m exactly, so that the correctness and practical application value of the ideal receiving depth concept have been verified. In addition, the depth classification effect based on the signals collected by the third element combination is better than the fifth array combination. It is mainly because the depth classification effect is related to the depth summation of the element combination and related to the depth of the first element. Extremely large depth of the first element results in poor classification effect and even it is not able to classify the target depth. Considering the influence of the depth summation of the element combination and the depth of its first element on the performance of the algorithm, the depth classification effect using the third element combination is the best. The above results are consistent with the outcome obtained from the simulation results.However, using the third line spectrum signal for target depth classification, a good classification effect cannot be obtained because this segment of signal is less continuous than the other received signals, and the received SNR is lower resulting in poor performance of the source depth classification algorithm. So, accurate depth classification results cannot be obtained.Through the analysis of the data in Table2 and Figure 17, it can be found that the target depth classification effect is directly related to the line spectrum continuity and the received SNR. The better line spectrum continuity or the larger received signal-to-noise ratio results in a better target depth classification effect.


## 5. Conclusions

This paper analyzes the impact of receiving depth on the performance of the target depth classification algorithm through simulation and sea experiment. After a sufficient theoretical study on the effective depth model, it is verified that when the sensors are at the ideal receiving depth, the performance of the target depth classification algorithm based on the first two modes is significantly better than that of the sensors placed at the nonideal receiving depth. The depth summation of the sensors and the depth value of the first sensor will affect the depth classification effect and critical depth value. The impact of the source frequency on the performance of the algorithm has also been researched. The influence of SNR and sound velocity profile parameters (negative gradient, thermocline intensity, thermocline thickness, and up-boundary depth) on the presented method has been deeply investigated. Considering the above factors comprehensively, the performance of the depth classification algorithm can be guaranteed to meet the actual requirements by suitably adjusting the receiving depths. The sea experiment was conducted for data processing. The improved performance of the target depth classification algorithm has been verified. The factors affecting the performance of the algorithm in real applications include line spectrum continuity and SNR at the receiver.There are still some drawbacks and limitations for our presented algorithm in real application scenario. Firstly, the condition of using this algorithm is that source frequency must only excite the first two modes. If the source frequency can excite three or more modes, the algorithm described in this paper is not applicable. This limitation can be overcome by using matched field processing theory in our other study. Secondly, the algorithm has certain requirements for depths of the receivers. The application value and performance of the algorithm will be limited by the actual underwater platform installation requirements. Therefore, according to the installation requirements of the actual underwater platform, it may be necessary to further improve the algorithm through more in-depth research.



---

*Source: 5928401-2019-12-04.xml*