Research on the Range-Frequency Interference Characteristics of Target Scattering Field in a Shallow Water Waveguide

: Based on the target scattering model and the normal mode theory in a shallow water waveguide, a mathematical model of the acoustic intensity under the coupling condition of the target and the environment was deduced, and the interference striations in the monostatic and bistatic conﬁguration were obtained by simulations. Further, a ﬁeld experiment was carried out in a lake, and the data were collected for a spherical target in two frequency bands, i.e., 20 − 40 kHz and 40 − 80 kHz. The experiment results showed good agreement with the simulations. The results of the simulation and the experiment showed that the existence of the target made the interference phenomenon in shallow water waveguides more complex, and the range-frequency interference characteristics were closely related to the conﬁguration of the sonar system, the target scattering function, the frequency range, and the target movement trajectory. These interference phenomena were found and theoretically analyzed in the paper. The research results can be applied to target detection and recognition, signal parameter estimation, target tracking, and other ﬁelds.


Introduction
In the shallow water environment, the images of acoustic intensity versus range and frequency show alternately dark and bright striation structures [1,2].This study started from passive acoustic fields.Bachman [3] and Vianna [4] studied the space−frequency domain distribution of the acoustic intensity of broadband sources at a fixed depth, and Tang [5] and Weng [6] observed the interference striations in the field experiments.All interference striations in these studies are distributed along straight lines.Chuprov [7] introduced the concept of a waveguide invariant, which related the intensity of the wave, the distance from the source to the receiver, and the frequency.Researchers used the waveguide invariant to characterize the interference striations in their studies and applied it to the source depth estimation [8,9], source localization [10,11], and passive ranging [12].
In active acoustic fields, interference striations still exist [13].When there is a target in active acoustic fields, the interference pattern of the sound field will become more complicated due to the coupling of the target and the environment.There have been a number of recent demonstrations of active sonar striations in underwater acoustic fields.Zurk's team [14,15] carried out a series of studies.They studied the coherence characteristics of the echo sound field of the active sonar target and proposed an invariance principle similar to that of passive sonar.In their research, the scattering function of the target is approximated by δ and 1 (the former is used to represent a non-mixed target without coupling between scattering modes; the latter is used to represent the scattering function of a target with coupling between scattering modes).The experiments are carried out with Appl.Sci.2022, 12, 9828 2 of 12 the glass sphere and the flat plate as the target, and the structural characteristics of the interference striations are obtained.Later, Tang [16] studied the interference structure of active sonar reverberation signals and proposed an interference striation enhancement algorithm based on a variable window STFT transformation, which made the reverberation interference structure clearer.Goldhahn [17] carried out a controlled tank experiment.
The interference striations features of the active sound field are of great significance and can be applied to many fields, such as target detection [18,19] and target tracking [20].
In the current research, the δ function and 1 are used to approximate the scattering function of a target.This approximation method does not contain the target parameters (shape, structure, material, etc.) and scattering angle information [21].This leads to the interference striations obtained, which do not reflect the target features.In addition, the impact of a configuration of active sonar and its operating parameter (frequency, wideband, etc.) on the interference striations has not been analyzed theoretically and experimentally in depth.All of the above limit the understanding and application of the active interference striations.Therefore, this paper pays close attention to the following problems: What are the different phenomena of interference striations between the monostatic configuration and the bistatic configuration?What is the effect of the target's scattering function on the interference striations?How do the sonar work parameters and the motion state of the target affect the interference striations?
The theoretical model of the interference striations was derived based on the normal mode theory and the target scattering theory in shallow water waveguides in Section 2. Simulations were performed for two scenarios, i.e., the monostatic configuration and the bistatic configuration, and the target scattering characteristic and the interference patterns for the two scenarios were obtained in Section 3. Four typical phenomena, i.e., the bending phenomena, the density phenomena, the slope phenomena, and the discontinuity phenomena, were obtained and analyzed theoretically in Section 4, and the above problems were answered in this part.Furthermore, the lake experiment was carried out, and the theory in this paper was verified by the measured data in Section 5.

Theoretical Model
A complete active sonar system consists of a source, a receiver, and a target.For the active sonar system in shallow water waveguides, the source transmits an acoustic signal, and the acoustic signal propagates forward.When it encounters the target, it scatters and the scattered field is received by the receiver.A front view of the relative positional relationship between the source, the target, and the receiver is shown in Figure 1, which includes two cases: Figure 1a is the monostatic configuration, and Figure 1b is the bistatic configuration.In Figure 1, r 1 is the distance from the source to the target, r 2 is the distance from the target to the receiver, the ocean depth is denoted as H, z s is the depth of the source, z t is the depth of the target, and z r is the depth of the receiver.
get is approximated by δ and 1 (the former is used to represent a non-mixed target without coupling between scattering modes; the latter is used to represent the scattering function of a target with coupling between scattering modes).The experiments are carried out with the glass sphere and the flat plate as the target, and the structural characteristics of the interference striations are obtained.Later, Tang [16] studied the interference structure of active sonar reverberation signals and proposed an interference striation enhancement algorithm based on a variable window STFT transformation, which made the reverberation interference structure clearer.Goldhahn [17] carried out a controlled tank experiment.
The interference striations features of the active sound field are of great significance and can be applied to many fields, such as target detection [18,19] and target tracking [20].In the current research, the δ function and 1 are used to approximate the scatter- ing function of a target.This approximation method does not contain the target parameters (shape, structure, material, etc.) and scattering angle information [21].This leads to the interference striations obtained, which do not reflect the target features.In addition, the impact of a configuration of active sonar and its operating parameter (frequency, wideband, etc.) on the interference striations has not been analyzed theoretically and experimentally in depth.All of the above limit the understanding and application of the active interference striations.Therefore, this paper pays close attention to the following problems: What are the different phenomena of interference striations between the monostatic configuration and the bistatic configuration?What is the effect of the target's scattering function on the interference striations?How do the sonar work parameters and the motion state of the target affect the interference striations?
The theoretical model of the interference striations was derived based on the normal mode theory and the target scattering theory in shallow water waveguides in Section 2. Simulations were performed for two scenarios, i.e., the monostatic configuration and the bistatic configuration, and the target scattering characteristic and the interference patterns for the two scenarios were obtained in Section 3. Four typical phenomena, i.e., the bending phenomena, the density phenomena, the slope phenomena, and the discontinuity phenomena, were obtained and analyzed theoretically in Section 4, and the above problems were answered in this part.Furthermore, the lake experiment was carried out, and the theory in this paper was verified by the measured data in Section 5.

Theoretical Model
A complete active sonar system consists of a source, a receiver, and a target.For the active sonar system in shallow water waveguides, the source transmits an acoustic signal, and the acoustic signal propagates forward.When it encounters the target, it scatters and the scattered field is received by the receiver.A front view of the relative positional relationship between the source, the target, and the receiver is shown in Figure 1, which includes two cases: Figure 1a is the monostatic configuration, and Figure 1b is the bistatic configuration.In Figure 1, r1 is the distance from the source to the target, r2 is the distance from the target to the receiver, the ocean depth is denoted as H, zs is the depth of the source, zt is the depth of the target, and zr is the depth of the receiver.The acoustic path in the sonar configuration has two parts: the source to the target path (r 1 ) and the target to the receiver path (r 2 ).For a target in the far field, the field point acoustic pressure for small grazing angles can be expressed as: p(r 1 , z t ; r 2 , z r ; z s The interpretation of Equation ( 1) is straightforward.The acoustic pressure amplitude of the scattered field is the product of the far-field scattering function and the mode function of the incident and scattered fields.The azimuth of the scattering function is obtained jointly from the source and the field point.Where C = ie iπ/4 /ρ √ 8π is the normalization constant, k m and k n are the horizontal wavenumbers of the mth and nth normal modes, respectively.
, it is the effect of the mth propagating mode on the sound field, when the acoustic signal propagates between the source and the target.
is the effect of the nth scattering mode, when the acoustic signal is coupled with the target and propagates to the receiver.The scattering function [21] S(α m , Φ 0 ; α n , Φ), where α m and Φ 0 are the grazing and azimuth angles of the incident mode m, respectively, and α n and Φ are the grazing and azimuth angles, respectively, corresponding to the scattered mode, is an important property of the target that can reflect the shape, structure, and material of the target, as well as the scattering angle information.
The acoustic intensity of the field scattered from a far-field target can be written as I =|p| 2 .In order to simplify the expression, and S mn = S(α m , Φ 0 ; α n , Φ) are defined.The acoustic pressure can be expressed as: p(r 1 , z t ; r 2 , z r ; z s ; ω) = ∑ m ∑ n B mn S mn e i(k m r 1 +k n r 2 ) . ( Using this equation, the sound field intensity I =|p| 2 at the receiver can be expressed as: B mn B pq S mn S pq e i(k m r 1 +k n r 2 −k p r 1 −k q r 2 ) . ( In this way, the intensity of the sound field can be quickly calculated with the help of the above equation to produce the range-frequency image of the acoustic intensity.In order to more intuitively explain how the interference pattern is formed, the acoustic pressure representation is converted into the following form: The cosine term in Equation ( 4) represents the interference between any pair of modes in the propagation paths (r 1 and r 2 ) of the two acoustic waves.It includes the mutual interference of any two incident modes, the mutual interference of any two scattered modes, and the mutual interference between any incident and scattered modes, which cause the interference phenomenon.
In the interference striations, the acoustic intensity satisfies dI = ∂I ∂w dw + ∂I ∂r dr.Referring to the definition of the sonar waveguide invariant, the active sonar waveguide invariant γ is obtained as: where is a function of the phase velocity, the target scattering properties, and the mode function.

Simulation
The sound field interference pattern in the ocean waveguide with a pressure-release surface and an ideal rigid bottom can be calculated from Equations ( 1) and ( 3).The ocean depth was 50 m.The ocean environment was modeled as an isovelocity channel with a sound speed c of 1500 m/s.The source transmitted the linear frequency-modulated signal (LFM) in three bands: 400−800 Hz, 4000−8000 Hz, and 40−80 kHz.Based on the above model, simulations were performed for two scenarios, i.e., the field with a target observed by an active sonar in monostatic configuration and the field with a target observed by an active sonar in bistatic configuration.The interference patterns for the two scenarios were obtained.

The Target Scattering Simulation
The target was a sphere with a diameter of 0.5 m in this paper.Figure 2a-c where  is a function of the phase velocity, the target scattering properties, and the mode function.

Simulation
The sound field interference pattern in the ocean waveguide with a pressure-release surface and an ideal rigid bottom can be calculated from Equations ( 1) and ( 3).The ocean depth was 50 m.The ocean environment was modeled as an isovelocity channel with a sound speed c of 1500 m/s.The source transmitted the linear frequency-modulated signal (LFM) in three bands: 400−800 Hz, 4000−8000 Hz, and 40−80 kHz.Based on the above model, simulations were performed for two scenarios, i.e., the field with a target observed by an active sonar in monostatic configuration and the field with a target observed by an active sonar in bistatic configuration.The interference patterns for the two scenarios were obtained.

The Target Scattering Simulation
The target was a sphere with a diameter of 0.5 m in this paper.Figure 2a-c  In Figure 2, it can be seen that the magnitude of the scattering function was a function of the frequency and the angle.At low frequencies, the scattering function amplitude did not change much between different angles, and the scattering function amplitude did not change much during the movement of the target.The greater the frequency, the greater the difference in amplitude between different angles, and small changes in angle led to greater changes in amplitude, as shown in Figure 1c.

Monostatic Mode Simulation
The monostatic mode simulation is shown in Figure 1a.The target, the source, and the receiver were all at the same depth of 25 m.The source and the receiver were fixed on the same platform, and the distance between the two was small (0.5 m).The target moved in a direction away from the source.In Figure 2, it can be seen that the magnitude of the scattering function was a function of the frequency and the angle.At low frequencies, the scattering function amplitude did not change much between different angles, and the scattering function amplitude did not change much during the movement of the target.The greater the frequency, the greater the difference in amplitude between different angles, and small changes in angle led to greater changes in amplitude, as shown in Figure 1c.

Monostatic Mode Simulation
The monostatic mode simulation is shown in Figure 1a.The target, the source, and the receiver were all at the same depth of 25 m.The source and the receiver were fixed on the same platform, and the distance between the two was small (0.5 m).The target moved in a direction away from the source.
Figure 3 shows an interference pattern observed by the above configuration for frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.When the target moved a range of 90 m to a range of 60 m, the shape of the striations was similar to the trajectory of the target.The interference pattern of the high-frequency part of the band was wider than that of the low-frequency part, and the striations in the high-frequency band were close.The interference pattern of the active sonar monostatic configuration was basically the same as that of the non-target experiment field.
frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.When the target moved a range of 90 m to a range of 60 m, the shape of the striations was similar to the trajectory of the target.The interference pattern of the high-frequency part of the band was wider than that of the low-frequency part, and the striations in the high-frequency band were close.The interference pattern of the active sonar monostatic configuration was basically the same as that of the non-target experiment field.

Bistatic Mode Simulation
For the bistatic configuration in Figure 1b, the target, the source and the receiver were all at the same depth of 25 m.The distance between the fixed source and the fixed receiver was 60 m.The spherical target with a diameter of 0.5 m moved along the line which was 60 m away from the line where the source and the receiver were located.
Figure 4 shows the interference pattern observed by the above configuration for frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.As the target moved from 78 m to 60 m, the shape of the striations followed the shape of the target's trajectory.The striations bent more in this configuration, compared to that in the monostatic configuration.The striations were discontinuous, and the sound field intensity of the high-frequency part was more pronounced than that of the low-frequency part.The interference striations were denser at high-frequency bands.In addition, the striations were wide in the high-frequency part and narrow in the low-frequency part at one frequency band.

Bistatic Mode Simulation
For the bistatic configuration in Figure 1b, the target, the source and the receiver were all at the same depth of 25 m.The distance between the fixed source and the fixed receiver was 60 m.The spherical target with a diameter of 0.5 m moved along the line which was 60 m away from the line where the source and the receiver were located.
Figure 4 shows the interference pattern observed by the above configuration for frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.As the target moved from 78 m to 60 m, the shape of the striations followed the shape of the target's trajectory.The striations bent more in this configuration, compared to that in the monostatic configuration.The striations were discontinuous, and the sound field intensity of the high-frequency part was more pronounced than that of the low-frequency part.The interference striations were denser at high-frequency bands.In addition, the striations were wide in the high-frequency part and narrow in the low-frequency part at one frequency band.
frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.When the target moved a range of 90 m to a range of 60 m, the shape of the striations was similar to the trajectory of the target.The interference pattern of the high-frequency part of the band was wider than that of the low-frequency part, and the striations in the high-frequency band were close.The interference pattern of the active sonar monostatic configuration was basically the same as that of the non-target experiment field.

Bistatic Mode Simulation
For the bistatic configuration in Figure 1b, the target, the source and the receiver were all at the same depth of 25 m.The distance between the fixed source and the fixed receiver was 60 m.The spherical target with a diameter of 0.5 m moved along the line which was 60 m away from the line where the source and the receiver were located.
Figure 4 shows the interference pattern observed by the above configuration for frequencies: 400−800 Hz, 4−8 kHz, and 40−80 kHz.As the target moved from 78 m to 60 m, the shape of the striations followed the shape of the target's trajectory.The striations bent more in this configuration, compared to that in the monostatic configuration.The striations were discontinuous, and the sound field intensity of the high-frequency part was more pronounced than that of the low-frequency part.The interference striations were denser at high-frequency bands.In addition, the striations were wide in the high-frequency part and narrow in the low-frequency part at one frequency band.

Experimental Verification 4.1. Overview
In November 2020, a test was carried out at the Moganshan Lake (Figure 5d) in Huzhou City, Zhejiang Province, using a spherical target in a monostatic configuration.The device and the test scene are shown in Figure 5.The lake depth was 14 m.The environment was modeled as an isovelocity channel with a sound speed c of 1500 m/s.The target (Figure 5c), the fixed source (Figure 5a), and the fixed receiver (Figure 5b) were at the same depth of 7 m.The distance between the source and the receiver was 4 m.The spherical target with a diameter of 1.2 m moved in the direction indicated by the arrow.This test was performed at frequency bands of 20−40 kHz and 40−80 kHz.The test configuration is shown in Figure 6.

Overview
In November 2020, a test was carried out at the Moganshan Lake (Figure 5d) in Huzhou City, Zhejiang Province, using a spherical target in a monostatic configuration.The device and the test scene are shown in Figure 5.The lake depth was 14 m.The environment was modeled as an isovelocity channel with a sound speed c of 1500 m/s.The target (Figure 5c), the fixed source (Figure 5a), and the fixed receiver (Figure 5b) were at the same depth of 7 m.The distance between the source and the receiver was 4 m.The spherical target with a diameter of 1.2 m moved in the direction indicated by the arrow.This test was performed at frequency bands of 20−40 kHz and 40−80 kHz.The test configuration is shown in Figure 6.The target echo signal received by the receiver was processed using the following steps: the signal was preprocessed to remove the interference of the noise and the clutter, the track of the target was obtained, a spectrum analysis was performed to obtain the range-frequency image of the echo signal, and the striations characteristics were analyzed.In November 2020, a test was carried out at the Moganshan Lake (Figure 5d) in Huzhou City, Zhejiang Province, using a spherical target in a monostatic configuration.The device and the test scene are shown in Figure 5.The lake depth was 14 m.The environment was modeled as an isovelocity channel with a sound speed c of 1500 m/s.The target (Figure 5c), the fixed source (Figure 5a), and the fixed receiver (Figure 5b) were at the same depth of 7 m.The distance between the source and the receiver was 4 m.The spherical target with a diameter of 1.2 m moved in the direction indicated by the arrow.This test was performed at frequency bands of 20−40 kHz and 40−80 kHz.The test configuration is shown in Figure 6.The target echo signal received by the receiver was processed using the following steps: the signal was preprocessed to remove the interference of the noise and the clutter, the track of the target was obtained, a spectrum analysis was performed to obtain the range-frequency image of the echo signal, and the striations characteristics were analyzed.The target echo signal received by the receiver was processed using the following steps: the signal was preprocessed to remove the interference of the noise and the clutter, the track of the target was obtained, a spectrum analysis was performed to obtain the range-frequency image of the echo signal, and the striations characteristics were analyzed.Figures 7 and 8 show that the striations shape of the simulation results was consistent with that of the experimental results.Figures 7b and 8b show the intensity of the experimental striations varied with frequency.The target moved from near to far, and due to the attenuation characteristics of the signal in water, the intensity of the signal trace was weak at a longer range.The shape of the striations followed the shape of the target's trajectory.The interference striations at the high frequency were denser than those at the low frequency.

The Striations Interference Phenomenon and Its Formation Mechanism
From the above simulation and experimental results, four typical phenomena were observed, namely the bending, the density, the slope, and the discontinuities.These mechanisms were analyzed theoretically as follows.

The Phenomenon of Striations Bending and Its Generation Mechanism
In the simulation results for the monostatic configuration (Figure 3), the striations were straight in all frequency bands.In Section 4, the striations obtained from the monostatic configuration experiments were also straight.However, in the bistatic configuration (Figure 4), the simulation striations became bent at every frequency band.In order to explain the phenomenon, we started with the definition of the waveguide invariant.From Equation (5), we obtained the relationship between ω and r, that is: Figures 7 and 8 show that the striations shape of the simulation results was consistent with that of the experimental results.Figures 7b and 8b show the intensity of the experimental striations varied with frequency.The target moved from near to far, and due to the attenuation characteristics of the signal in water, the intensity of the signal trace was weak at a longer range.The shape of the striations followed the shape of the target's trajectory.The interference striations at the high frequency were denser than those at the low frequency.

The Striations Interference Phenomenon and Its Formation Mechanism
From the above simulation and experimental results, four typical phenomena were observed, namely the bending, the density, the slope, and the discontinuities.These mechanisms were analyzed theoretically as follows.

The Phenomenon of Striations Bending and Its Generation Mechanism
In the simulation results for the monostatic configuration (Figure 3), the striations were straight in all frequency bands.In Section 4, the striations obtained from the monostatic configuration experiments were also straight.However, in the bistatic configuration (Figure 4), the simulation striations became bent at every frequency band.In order to explain the phenomenon, we started with the definition of the waveguide invariant.From Equation (5), we obtained the relationship between ω and r, that is: The invariants γ and ∆ω/∆r b were constants, and the relationship between ω and the distance r was linear, indicating that the striations shape in the range-frequency spectrogram was approximately consistent with the target motion trajectory.This was verified by the simulation for two scenarios (Figures 3 and 4).The specific factor was determined by the value of the invariant and ∆ω/∆r b .The white dotted line in Figures 4a, 7b and 8b were the striations predicted by Equation (6).The prediction based on the invariant showed good agreement with the interference pattern obtained by the mathematical model of the sound field interference.
The acoustic path in active sonar had two parts: the source to the target path and the target to the receiver path.For bistatic configuration, the source and the receiver were at two separate locations, and there was a big difference between the distance from the source to the target and the distance from the target to the receiver.The variation of the total distance was nonlinear.When the two paths were combined, the motion trajectory appeared curved.The shape of the striation was consistent with the target trajectory, so the striations appeared the bending phenomenon.However, the devices were in the same location in the monostatic configuration, so the variations of the distance from the target to the source and from the target to the receiver were the same.The variation of the total distance was linear, so the striation patterns were straight.Figures 7 and 8 show the experimental striations had a slight bending compared with the simulation results, because during the lake test, the target sphere could not move completely along the specified straight line due to the actual conditions.The sound source, the receiver, and the target were no longer strictly distributed according to the monostatic configuration.The nonlinear variation of the distance made the striations slightly curved.

The Phenomenon of the Stripe Density and Its Generation Mechanism
Figures 3 and 4 show the simulation striations for the monostatic configuration and the bistatic configuration.The striations became denser from low frequency to high frequency.The striations were sparse at 400−800 Hz and denser at the 4−8 kHz frequency range.When the frequency was in the range of 40−80 kHz, the striations were denser than those at 4−8 kHz.The density of the striations increased with the frequency.Figures 7 and 8 show the simulation striations and the experimental striations for the monostatic configuration.It can also be observed that the striations in the 20−40 kHz were denser than those in the 40−80 kHz frequency range.
In order to explain this, we took the sound propagation process from the target to the receiver.The number of modes N excited by the signal was related to the frequency f, the depth of the environment H, and the underwater sound speed c.The relationship among them can be expressed as: The number of the modes is proportional to the frequency of the signal.In the same environment, the number of the modes of the high-frequency part is more than that of the low-frequency part, so the interference pattern of the high-frequency part is more complicated.For the simulation scenarios in Section 3, the maximum numbers of modes excited by the source at 400−800 Hz, 4−8 kHz, and 40−80 kHz frequency bands were about 53, 533, and 5333, respectively.Similar to simulation experiments, the maximum numbers of modes excited by the source at 20−40 kHz and 40−80 kHz frequency bands were about 746 and 1493, respectively, for the monostatic experimental configuration in Section 4. In comparison, the number of modes excited by signals at high-frequency bands was larger than that of signals at low-frequency bands.The mutual interference between the modes was more obvious, and the interference striations were denser.
In both Sections 3 and 4, the striations were wide in the high part and narrow in the low part at one frequency band.For example, in Figure 4a, the simulation striations interference length in the range of 400−500 Hz was shorter than that in the range of 600−700 Hz.In Figure 7b, the experimental striations interference length in the range of 20−27 kHz was shorter than that in the range of 27−40 kHz.According to Equation (4), the oscillation period of a sound field with a second-order mode is written as: This oscillation period is the interference length, and the interference length between modes is inversely proportional to the wavenumber difference.According to the definition of the horizontal wavenumbers, the wavenumber difference between two adjacent modes at the high-frequency part is smaller than the low-frequency part.The corresponding striations interference length is also longer.

The Phenomenon of Striations Discontinuity and Its Generation Mechanism
In the simulation experiment of the bistatic mode shown in Figure 4, there were discontinued in the interference pattern and the amplitude changed more drastically with the increase in frequency.We did not see these changes in the monostatic mode given in Figure 9 is the interference striation diagram obtained in the monostatic and bistatic configuration simulation experiments and the variation relationship diagram of the scattering function amplitude during the movement of the target.The variation of the sound field intensity showed agreement with the change of the scattering function amplitude.In Figure 9b, the amplitude of the scattering function was almost zero around 4000-5000 Hz, so the striations appeared discontinuous near the corresponding frequency.
responding striations interference length is also longer.

The Phenomenon of Striations Discontinuity and Its Generation Mechanism
In the simulation experiment of the bistatic mode shown in Figure 4, there were discontinued in the interference pattern and the amplitude changed more drastically with the increase in frequency.We did not see these changes in the monostatic mode given in Figure 3.The discontinuity of the striations was related to the magnitude of the scattering function shown in Figure 2.This can explain the reason for the difference in the discontinuity of the striations between the monostatic and bistatic configurations.For the monostatic configuration in Figure 9a, the source and the receiver were located on the same platform, the range between them was extremely small, and the scattering angle was basically kept around 180 degrees during the movement of the target, so the amplitude of the scattering function changed almost smoothly with the change of frequency.For the bistatic configuration shown in Figure 9b,c, the source and the receiver were located in two parts.When the target moved, the scattering angle changed, so the amplitude of the scattering function also changed.The sound field interference striations appeared discontinuous.Compar- This can explain the reason for the difference in the discontinuity of the striations between the monostatic and bistatic configurations.For the monostatic configuration in Figure 9a, the source and the receiver were located on the same platform, the range between them was extremely small, and the scattering angle was basically kept around 180 degrees during the movement of the target, so the amplitude of the scattering function changed almost smoothly with the change of frequency.For the bistatic configuration shown in Figure 9b,c, the source and the receiver were located in two parts.When the target moved, the scattering angle changed, so the amplitude of the scattering function also changed.The sound field interference striations appeared discontinuous.Comparing Figure 9b,c, the degrees of discontinuities striations with different frequencies were also different.For the bistatic configuration, the scattering angle of the object constantly changed during the target movement, the scattering function amplitude changed drastically at high frequencies.In addition, the degree of discontinuity striations became more obvious.Therefore, the striations discontinuity phenomenon under the condition of 40-80 kHz was more obvious than that at 4-8 kHz.The experimental striations shown in Figures 7b and 8b were the results of the monostatic configuration test, but they showed the phenomenon of striations discontinuity.On the one hand, the signal-to-noise ratio of the signal was large at high frequencies, so the structure of the striations was clearer than at the low frequency.On the other hand, during the lake test, the target sphere did not move along the specified straight line due to the actual conditions.The scattering angle changed while the target sphere moved, so the experimental striations showed the phenomenon of striations discontinuity.

The Relationship between Striations Slope and Its Generation Mechanism
In either the bistatic or monostatic configuration, the distribution of the interference striations was very regular, and there seemed to be a constrained relationship between the slope of the striations and the range and frequency of the source.Assuming the target was far away from the source, Figure 10a shows the relationship between the intensity of the sound field and the range when the frequencies were different.Figure 10b shows the relationship between the intensity of the sound field and the frequency of the source signal when the range changed.
In either the bistatic or monostatic configuration, the distribution of the interference striations was very regular, and there seemed to be a constrained relationship between the slope of the striations and the range and frequency of the source.Assuming the target was far away from the source, Figure 10a shows the relationship between the intensity of the sound field and the range when the frequencies were different.Figure 10b shows the relationship between the intensity of the sound field and the frequency of the source signal when the range changed.With the increase of the distance, the frequency increased.The value of  in Equation ( 5) was zero for the simulation results of the monostatic mode.The mathematical relationship for the waveguide invariant under the monostatic configuration is the same as that obtained for the passive sonar.For the monostatic simulation results in Figure 3, we obtained the values of the waveguide invariant under different distances and different frequencies as shown in Table 1.With the increase of the distance, the frequency increased.The value of in Equation ( 5) was zero for the simulation results of the monostatic mode.The mathematical relationship for the waveguide invariant under the monostatic configuration is the same as that obtained for the passive sonar.For the monostatic simulation results in Figure 3, we obtained the values of the waveguide invariant under different distances and different frequencies as shown in Table 1.For the bistatic simulation results in Figure 4, the calculation results of the waveguide invariants are shown in Table 2.In the bistatic configuration, the slope value ∆ω/∆r is replaced by the slope of the tangent line for curved striations.Similar to the simulation, the value of the lake test striations waveguide invariant was calculated.Taking Figure 7b as an example, for the striations A (f was 40 kHz, r was 51 m, and ∆ω/∆r was 817 Hz/m), the calculated value of the waveguide invariant was 1.042.In the same way, the waveguide invariant value of the striations B was calculated to be 1.057.According to the above idea, in Figure 7, the value of the waveguide invariant for the lake test was about 1.046 at 20−40 kHz, and in Figure 8, the value at 40−80 kHz was about 1.104.
Comparing the data in the above table, it can be seen that the value of the waveguide invariant at long distances tended to be 1.When the frequency was the same, the slope ∆ω/∆r was small when the distance was large.When the distance was the same, the slope ∆ω/∆r was large when the frequency was large.That mean if any element between the distance and the frequency changed, the other element compensated for it, so that there was a restrictive relationship between the changes.

Conclusions
In this paper, based on the target scattering characteristics and the normal mode theory in a shallow water waveguide, a mathematical model of the acoustic intensity under the coupling condition of the target and the environment was deduced.Through simulation and lake tests, the range-frequency interference patterns of the target scattering field in the shallow water waveguide were obtained, and four typical phenomena named the bending, the density, the slope, and the discontinuities were found and explained, as follows:

•
The relationship between the range and the frequency changed linearly, so the shape of the striations followed the shape of the target's trajectory.The propagation path of the sound wave in the system and the position distributions of the source and receiver made the interference striations of the bistatic configuration more complicated than that of the monostatic configuration;

•
The number of the modes of the high-frequency part was larger than that of the lowfrequency part, so the interference striations of the high-frequency part were denser than those of the low-frequency part; • There was a constraint relationship between the slope of the interference striations and the range and frequency.If any element between the range and the frequency changed, other elements compensated for it;

•
The interference striations were affected by the amplitude of the scattering function.The higher the frequency, the more severe the changes in the amplitude of the scattering function with the incident and reflected angles.The interference striations in the acoustic intensity appeared discontinuous due to the change in the amplitude of the scattering function.
In conclusion, the interference striations of the target scattering field in the shallow water waveguide were closely related to the configuration of the sonar system, the target scattering function, the frequency range, the target motion trajectory, etc.An in-depth analysis of the mechanisms causing these interference phenomena is of great value to target detection and recognition, signal parameter estimation, and target tracking.

Figure 3 .
Figure 3.The simulation results in the monostatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 4 .
Figure 4.The simulation results in the bistatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 3 .
Figure 3.The simulation results in the monostatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 3 .
Figure 3.The simulation results in the monostatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 4 .
Figure 4.The simulation results in the bistatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 4 .
Figure 4.The simulation results in the bistatic configuration: (a) the interference pattern for 400−800 Hz; (b) the interference pattern for 4−8 kHz; (c) the interference pattern for 40−80 kHz; (d) the motion trajectory.

Figure 5 .
Figure 5. Test device and test scene: (a) the source; (b) the receiver; (c) the target sphere; (d) the test scene.

Figure 6 .
Figure 6.Schematic diagram of the experiment.

Figures 7 Figure 5 .
Figures 7 and 8 show a comparison between the experiment results and the simulation results under the condition shown in Section 4.1.

Figure 5 .
Figure 5. Test device and test scene: (a) the source; (b) the receiver; (c) the target sphere; (d) the test scene.

Figure 6 .
Figure 6.Schematic diagram of the experiment.

Figures 7 Figure 6 .
Figures 7 and 8 show a comparison between the experiment results and the simulation results under the condition shown in Section 4.1.

Figures 7
Figures 7 and 8 show a comparison between the experiment results and the simulation results under the condition shown in Section 4.1.Appl.Sci.2022, 12, x FOR PEER REVIEW 7 of 13

Figure 7 .
Figure 7.Comparison between the experiment result and the simulation result at 20−40 kHz: (a) simulation result; (b) experiment result; (c) the motion trajectory.

Figure 7 .
Figure 7.Comparison between the experiment result and the simulation result at 20−40 kHz: (a) simulation result; (b) experiment result; (c) the motion trajectory.

Figure 7 .Figure 8 .
Figure 7.Comparison between the experiment result and the simulation result at 20−40 kHz: (a) simulation result; (b) experiment result; (c) the motion trajectory.

Figure 8 .
Figure 8.Comparison between the experiment result and the simulation result at 40−80 kHz: (a) simulation result; (b) experiment result; (c) the motion trajectory.

Figure 3 .
Figure 3.The discontinuity of the striations was related to the magnitude of the scattering function shown in Figure 2.Figure9is the interference striation diagram obtained in the monostatic and bistatic configuration simulation experiments and the variation relationship diagram of the scattering function amplitude during the movement of the target.The variation of the sound field intensity showed agreement with the change of the scattering function amplitude.In Figure9b, the amplitude of the scattering function was almost zero around 4000-5000 Hz, so the striations appeared discontinuous near the corresponding frequency.

Figure 9 Figure 9 .
Figure 9.Comparison of interference striations and scattering function amplitudes: (a) the result of the monostatic configuration for 4−8 kHz; (b) the result of the bistatic configuration for 4−8 kHz; (c) the result of the bistatic configuration for 40−80 kHz.

Figure 9 .
Figure 9.Comparison of interference striations and scattering function amplitudes: (a) the result of the monostatic configuration for 4−8 kHz; (b) the result of the bistatic configuration for 4−8 kHz; (c) the result of the bistatic configuration for 40−80 kHz.

Figure 10 .
Figure 10.The sound field intensity variation: (a) the variations at different distances; (b) the variations at different frequencies.

Figure 10 .
Figure 10.The sound field intensity variation: (a) the variations at different distances; (b) the variations at different frequencies.

Table 1 .
The slopes and waveguide invariants for the monostatic configuration.

Table 1 .
The slopes and waveguide invariants for the monostatic configuration.

Table 2 .
The slopes and invariants for the bistatic configuration.