A system and method for estimating the motion parameters of an elastic seabed sound source based on the interference phenomenon of sound field

Abstract

This invention proposes a system and method for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena. First, the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide are analyzed. Then, based on the analyzed interference structure characteristics, a two-dimensional Fourier transform of the time-frequency interference structure is performed to extract the spectral density peaks. Finally, based on the analyzed normal mode characteristics of the elastic seabed waveguide and the calculated coordinates of the spectral density peak points, the motion parameters of the acoustic source are estimated. This invention analyzes the surface wave characteristics of this ubiquitous marine environment based on the normal mode characteristics of the elastic seabed waveguide, identifies the normal mode that plays a major role in the acoustic field interference structure of the elastic seabed waveguide, optimizes the original two-dimensional Fourier transform algorithm for estimating the motion parameters of the acoustic source, and accurately achieves the estimation of the motion parameters of the acoustic source in this marine environment waveguide.

Description

Technical Field

[0001] This invention belongs to the field of artificial intelligence technology, specifically, it relates to a system and method for estimating the motion parameters of an elastic seabed sound source based on the phenomenon of acoustic field interference. Background Technology

[0002] Compared to liquid seabed models, the ocean covering the solid crust is more consistent with ocean waveguide models of elastic seabeds. However, due to its many unique properties, it remains a subject of ongoing research. According to seismological theory, a special wave mode exists in elastic seabed waveguides: surface waves. As early as the 19th century, Lord Rayleigh theoretically elucidated the possibility of surface acoustic waves propagating along the free surface of a solid. In 1924, Stoneley studied surface waves existing between two layers of elastic solid half-spaces and provided the corresponding characteristic equations. Waves that generally appear at the fluid-elastic solid boundary are called Scholte waves, while those appearing at the interface between elastic solids are called Stoneley waves.

[0003] In 2017, Russian scientists GNKuznetsov et al. proposed a velocimetry and ranging algorithm using two-dimensional Fourier transform on time-frequency domain interferometric images. This algorithm requires only a single-vector hydrophone and achieved good results under low signal-to-noise ratio conditions. Subsequently, SAPerselkov et al. verified the applicability of this algorithm at higher frequencies, broadening its applicability to the estimation of sound source motion parameters.

[0004] However, the marine environment studied by Kuznetsov, Pereselkov, and others was a single-layer Pekeris waveguide located on a semi-infinite liquid seabed medium. They did not conduct in-depth analysis of the normal mode characteristics, nor did they study the applicability of this algorithm in marine conditions where sedimentary layers and rocky substrates with elastic medium characteristics are widely present. Summary of the Invention

[0005] This invention addresses the problem that the modal characteristics of normal waves under elastic seabed conditions differ significantly from those under liquid seabed conditions. Based on the normal wave method using the seabed reflection coefficient, it analyzes the modal characteristics of normal waves under elastic seabed conditions, optimizes the sound source motion parameter estimation algorithm based on time-frequency interference fringes, and then proposes a sound source motion parameter estimation system and method based on sound field interference phenomena.

[0006] This invention is achieved through the following technical solution:

[0007] A method for estimating the motion parameters of an elastic seabed sound source based on acoustic field interference:

[0008] The method specifically includes the following steps:

[0009] Step 1: Analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment;

[0010] Step 2: Based on the interference structure characteristics of Step 1, perform a two-dimensional Fourier transform of the time-frequency interference structure to extract the spectral density peaks;

[0011] Step 3: Based on the normal mode characteristics of the elastic seabed waveguide in Step 1 and the coordinates of the peak spectral density points calculated in Step 2, the motion parameters of the sound source are estimated.

[0012] Furthermore, in step one,

[0013] The elastic marine waveguide environment includes a seawater layer, an elastic sedimentary layer, and a rocky substrate;

[0014] The medium in the seawater layer is an ideal fluid, and its sound speed at depth z is c(z), and its density ρ is independent of depth and distance.

[0015] The parameters of the elastic deposition layer include the longitudinal wave velocity c. p1 transverse wave speed of sound c s1 Density ρ1, longitudinal wave propagation loss α p1 transverse wave sound speed α s1 With the sedimentary depth H1;

[0016] The parameters of the rocky basement include the longitudinal wave velocity c. p2 transverse wave speed of sound c s2 Density ρ2, longitudinal wave propagation loss α p2 transverse wave sound speed α s2 .

[0017] Furthermore, in step one,

[0018] The expression for normal modes in a general waveguide can be rewritten as an expression for the amplitude of normal modes and the exponent e in an interference structure:

[0019]

[0020] Where the normal mode amplitude A n The expression for (ω,r) is:

[0021]

[0022] Write the normal wave interference amplitude I based on the sound pressure field expression. mn (ω,r) expression:

[0023]

[0024] Based on the normal mode interference amplitudes of each mode in formulas (1), (2), and (3), the normal mode that plays a major role in the acoustic field interference structure in the current elastic seabed waveguide is analyzed.

[0025] Furthermore, in step two,

[0026] A two-dimensional Fourier transform is performed based on the sound intensity LOFAR spectrum I(ω,t) of the hydrophone received signal, as shown in equation (4):

[0027]

[0028] Where v r ω is the radial velocity of the sound source, Δt is the observation time, ω0 is the center angular frequency of the observation band, Δω is the width of the observation band, ν is the frequency shift in the Fourier transform factor, and τ is the time delay in the Fourier transform factor.

[0029] Furthermore, in step two,

[0030] The amplitude distribution of the spectral density F(τ,ν) in different angular coefficient directions in the frequency shift-time delay domain is written as:

[0031]

[0032] Based on the maximum value of the angle coefficient obtained by the detector in the formula, determine the coordinates (τ1,ν1) of the peak point (τ1,ν1) that is closest to the origin and formed by adjacent modal normal waves in that direction.

[0033] Furthermore, in step three,

[0034] Based on the coordinates of the spectral density peak point obtained in step two and the normal mode characteristics of the elastic seabed, an expression for estimating the sound source motion parameters is given:

[0035]

[0036] Where N f The first normal mode order carrying the main energy is obtained from step one, N is the normal mode order carrying the main energy obtained from the analysis, and the β value is the interference invariant value calculated based on the formula:

[0037]

[0038] Where k n(n+1) It is the difference between the eigenvalues ​​of adjacent modal normals.

[0039] A system for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena:

[0040] The system includes: a feature analysis module, a Fourier transform module, and a motion parameter estimation module;

[0041] The feature analysis module is used to analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment.

[0042] The Fourier transform module performs a two-dimensional Fourier transform of the time-frequency interference structure based on the analyzed interference structure characteristics, and extracts the spectral density peak value.

[0043] The motion parameter estimation module estimates the motion parameters of the sound source based on the analyzed normal mode characteristics of the elastic seabed waveguide and the calculated coordinates of the peak points of the spectral density.

[0044] Furthermore, the elastic marine waveguide environment analyzed by the feature analysis module includes a seawater layer, an elastic sedimentary layer, and a rocky substrate;

[0045] The medium in the seawater layer is an ideal fluid, and its sound speed at depth z is c(z), and its density ρ is independent of depth and distance.

[0046] The parameters of the elastic deposition layer include the longitudinal wave velocity c. p1 transverse wave speed of sound c s1 Density ρ1, longitudinal wave propagation loss α p1 transverse wave sound speed α s1 With the sedimentary depth H1;

[0047] The parameters of the rocky basement include the longitudinal wave velocity c. p2 transverse wave speed of sound c s2 Density ρ2, longitudinal wave propagation loss α p2 transverse wave sound speed α s2 .

[0048] An electronic device includes a memory and a processor, the memory storing a computer program, the processor executing the computer program to implement the steps of any of the methods described above.

[0049] A computer-readable storage medium for storing computer instructions that, when executed by a processor, implement the steps of any of the methods described above.

[0050] Beneficial effects of the invention

[0051] Based on the normal mode characteristics of elastic seabed waveguides, this invention analyzes the surface wave characteristics of this ubiquitous marine environment, identifies the normal mode that plays a major role in the acoustic field interference structure in elastic seabed waveguides, optimizes the original two-dimensional Fourier transform algorithm for estimating sound source motion parameters, and accurately realizes the estimation of sound source motion parameters in this marine environment waveguide. Attached Figure Description

[0052] Figure 1 The images show the eigenfunction diagrams of different modes at different frequencies in an elastic seabed waveguide, obtained by a normal mode calculation program based on seabed parameter coefficients. (a) is the eigenfunction diagram for frequency a, and (b) is the eigenfunction diagram for frequency b.

[0053] Figure 2 The diagram shows the interference amplitude of different normal modes in the marine environmental waveguide in the embodiment.

[0054] Figure 3 The images are video interferograms, where (a) is the LOFAR spectrum of the received signal near the hydrophone and (b) is the LOFAR spectrum of the received signal far from the hydrophone.

[0055] Figure 4 The images show the spectral density results extracted after processing the video interference image, where (a) is the spectral density result close to the hydrophone and (b) is the spectral density result far from the hydrophone. Detailed Implementation

[0056] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0057] A method for estimating the motion parameters of an elastic seabed sound source based on acoustic field interference:

[0058] The method specifically includes the following steps:

[0059] Step 1: Analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment;

[0060] Step 2: Based on the interference structure characteristics of Step 1, perform a two-dimensional Fourier transform of the time-frequency interference structure to extract the spectral density peaks;

[0061] Step 3: Based on the normal mode characteristics of the elastic seabed waveguide in Step 1 and the coordinates of the peak spectral density points calculated in Step 2, the motion parameters of the sound source are estimated.

[0062] In step one,

[0063] The elastic marine waveguide environment includes a seawater layer, an elastic sedimentary layer, and a rocky substrate;

[0064] The medium in the seawater layer is an ideal fluid, and its sound speed at depth z is c(z), and its density ρ is independent of depth and distance.

[0065] The parameters of the elastic deposition layer include the longitudinal wave velocity c. p1 transverse wave speed of sound c s1 Density ρ1, longitudinal wave propagation loss α p1 transverse wave sound speed α s1 With the sedimentary depth H1;

[0066] The parameters of the rocky basement include the longitudinal wave velocity c. p2 transverse wave speed of sound c s2 Density ρ2, longitudinal wave propagation loss α p2 transverse wave sound speed α s2 .

[0067] In step one,

[0068] The expression for normal modes in a general waveguide can be rewritten as an expression for the amplitude of normal modes and the exponent e in an interference structure:

[0069]

[0070] Where the normal mode amplitude A n The expression for (ω,r) is:

[0071]

[0072] Write the normal wave interference amplitude I based on the sound pressure field expression. mn (ω,r) expression:

[0073]

[0074] The eigenfunctions of different modes at different frequencies in an elastic seabed waveguide obtained by the normal mode calculation program based on seabed parameter coefficients are as follows: Figure 1 Analysis revealed that, under different frequency conditions within the frequency band, the surface wave modes appeared at different normal mode numbers.

[0075] Based on the amplitude images of the normal mode interference of each order of formulas (1), (2), and (3), the normal mode that plays a major role in the acoustic field interference structure in the current elastic seabed waveguide is analyzed.

[0076] In step two,

[0077] A two-dimensional Fourier transform is performed based on the sound intensity LOFAR spectrum I(ω,t) of the hydrophone received signal, as shown in equation (4):

[0078]

[0079] Where v rω is the radial velocity of the sound source, Δt is the observation time, ω0 is the center angular frequency of the observation band, Δω is the width of the observation band, ν is the frequency shift in the Fourier transform factor, and τ is the time delay in the Fourier transform factor.

[0080] In step two,

[0081] The amplitude distribution of the spectral density F(τ,ν) in different angular coefficient directions in the frequency shift-time delay domain is written as:

[0082]

[0083] Based on the maximum value of the angle coefficient obtained by the detector in the formula, determine the coordinates (τ1,ν1) of the peak point (τ1,ν1) that is closest to the origin and formed by adjacent modal normal waves in that direction.

[0084] In step three,

[0085] Based on the coordinates of the peak spectral density point obtained in step two and the normal mode characteristics of the elastic seabed, an expression for estimating the motion parameters of the sound source under this environment is given:

[0086]

[0087] Where N f The first normal mode order carrying the main energy is obtained from step one, N is the normal mode order carrying the main energy obtained from the analysis, and the β value is the interference invariant value calculated based on the formula:

[0088]

[0089] Where k n(n+1) It is the difference between the eigenvalues ​​of adjacent modal normals.

[0090] A system for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena:

[0091] The system includes: a feature analysis module, a Fourier transform module, and a motion parameter estimation module;

[0092] The feature analysis module is used to analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment.

[0093] The Fourier transform module performs a two-dimensional Fourier transform of the time-frequency interference structure based on the analyzed interference structure characteristics, and extracts the spectral density peak value.

[0094] The motion parameter estimation module estimates the motion parameters of the sound source based on the analyzed normal mode characteristics of the elastic seabed waveguide and the calculated coordinates of the peak points of the spectral density.

[0095] The elastic marine waveguide environment analyzed by the feature analysis module includes a seawater layer, an elastic sedimentary layer, and a rocky substrate.

[0096] The medium in the seawater layer is an ideal fluid, and its sound speed at depth z is c(z), and its density ρ is independent of depth and distance.

[0097] The parameters of the elastic deposition layer include the longitudinal wave velocity c. p1 transverse wave speed of sound c s1 Density ρ1, longitudinal wave propagation loss α p1 transverse wave sound speed α s1 With the sedimentary depth H1;

[0098] The parameters of the rocky basement include the longitudinal wave velocity c. p2 transverse wave speed of sound c s2 Density ρ2, longitudinal wave propagation loss α p2 transverse wave sound speed α s2 .

[0099] Example:

[0100] The sound velocity in the seawater layer is set to be around 1500 m / s, and the seawater depth is 150 m; the sediment layer thickness H1 is 50 m, and the P-wave velocity is c. p1 The transverse wave speed changes from 1700 m / s to 1730 m / s, and the sound velocity c... s1 The velocity is 0 m / s, and the density ρ1 is 2 g / cm³. 3 Longitudinal wave propagation loss and transverse wave sound velocity α s1 Set to 0; P-wave velocity c in the rocky basement p2 The transverse wave speed is 3500 m / s, and the sound velocity is c. s2 The speed is 2200 m / s, and the density ρ2 is 2.5 g / cm³. 3 Longitudinal wave propagation loss α p2 With transverse wave sound speed α s2 Set to 0. According to the formula, the interference amplitudes of each normal mode of this marine environmental waveguide are as follows: Figure 2 As shown;

[0101] based on Figure 2 It can be observed that the main factors affecting the acoustic field interference structure in this waveguide are the normal modes of interference between (4,5), (5,6), (6,7) and (7,8), with the third normal mode being the surface wave. Their effect on the acoustic field at the receiving point far from the seabed is negligible.

[0102] In the simulation, the sound source depth and the hydrophone depth were both 20m. The sound sources were set to be close to and far from the hydrophone. The initial distance of the sound source close to the hydrophone was 10km, with a radial velocity of 4m / s; the initial distance of the sound source far from the hydrophone was 6.8km, with a radial velocity of 4m / s. The resulting LOFAR spectrum of the received signal is shown below.

[0103] Based on the two-dimensional Fourier transform of the formula and the extraction of the spectral density peak point direction, for Figure 3 The spectral density results extracted after video interference image processing are as follows: Figure 4 As shown:

[0104] according to Figure 4 The estimated coordinates of the spectral density peaks generated by adjacent modes, combined with the normal mode characteristics of the elastic seabed waveguide, yield the motion parameters of the sound source in the two examples as follows:

[0105]

[0106] Table 1. Estimation results of sound source motion parameters in the example of elastic seabed waveguide.

[0107] An electronic device includes a memory and a processor, the memory storing a computer program, the processor executing the computer program to implement the steps of any of the methods described above.

[0108] A computer-readable storage medium for storing computer instructions that, when executed by a processor, implement the steps of any of the methods described above.

[0109] The memory in this application embodiment can be volatile memory or non-volatile memory, or it can include both volatile and non-volatile memory. The non-volatile memory can be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), or flash memory. The volatile memory can be random access memory (RAM), which is used as an external cache. By way of example, but not limitation, many forms of RAM are available, such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDRSDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous linked dynamic random access memory (SLDRAM), and direct rambus RAM (DR RAM). It should be noted that the memory used in the methods described in this invention is intended to include, but is not limited to, these and any other suitable types of memory.

[0110] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented in software, it can be implemented, in whole or in part, as a computer program product. The computer program product includes one or more computer instructions. When the computer instructions are loaded and executed on a computer, all or part of the processes or functions described in the embodiments of this application are generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another via wired (e.g., coaxial cable, fiber optic, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium accessible to a computer or a data storage device such as a server or data center that integrates one or more available media. The available media may be magnetic media (e.g., floppy disks, hard disks, magnetic tapes), optical media (e.g., high-density digital video discs (DVDs)), or semiconductor media (e.g., solid-state disks (SSDs)).

[0111] In implementation, each step of the above method can be completed by integrated logic circuits in the processor's hardware or by instructions in software. The steps of the method disclosed in the embodiments of this application can be directly implemented by a hardware processor, or by a combination of hardware and software modules in the processor. The software modules can reside in random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or other mature storage media in the art. This storage medium is located in memory, and the processor reads information from the memory and, in conjunction with its hardware, completes the steps of the above method. To avoid repetition, detailed descriptions are omitted here.

[0112] It should be noted that the processor in the embodiments of this application can be an integrated circuit chip with signal processing capabilities. During implementation, each step of the above method embodiments can be completed by the integrated logic circuits in the processor's hardware or by instructions in software form. The processor can be a general-purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. It can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of this application. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the methods disclosed in the embodiments of this application can be directly embodied as being executed by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software modules can be located in random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or other mature storage media in the art. This storage medium is located in memory, and the processor reads the information in the memory and, in conjunction with its hardware, completes the steps of the above methods.

[0113] The foregoing has provided a detailed description of the elastic seabed sound source motion parameter estimation system and method based on acoustic field interference phenomena proposed in this invention. The principles and implementation methods of this invention have been explained. The above description of the embodiments is only for the purpose of helping to understand the method and core ideas of this invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this invention. Therefore, the content of this specification should not be construed as a limitation of this invention.

Claims

1. A method for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena, characterized in that: The method specifically includes the following steps: Step 1: Analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment; The elastic marine waveguide environment includes a seawater layer, an elastic sedimentary layer, and a rocky substrate; the normal mode modes that play a major role in the acoustic field interference structure in the current elastic seabed waveguide are analyzed. determining the order of the first order normal mode carrying the main energy and the order of the normal mode carrying the main energy ; Step 2: Based on the interference structure characteristics of Step 1, perform a two-dimensional Fourier transform of the time-frequency interference structure to extract the spectral density peaks; Step 3: Based on the normal mode characteristics of the elastic seabed waveguide in Step 1 and the coordinates of the peak points of the spectral density calculated in Step 2, the motion parameters of the sound source are estimated. Based on the coordinates of the spectral density peak point obtained in step two and the normal mode characteristics of the elastic seabed, an expression for estimating the sound source motion parameters is given: (6) The value is the interference invariant value calculated based on equation (7): It is the center angular frequency of the observation band. The coordinates of the peak points formed by adjacent modal normal waves; (7) in It is the difference between the eigenvalues ​​of adjacent modal normals.

2. The method according to claim 1, characterized in that: In step one, The seawater layer medium is an ideal fluid, at a depth of Its speed of sound is Its density It is independent of depth and distance; The parameters of the elastic deposition layer include the longitudinal wave velocity. transverse wave speed of sound ,density Longitudinal wave propagation loss transverse wave speed of sound With the depth of the sedimentary layer ; The parameters of the rocky basement include longitudinal wave velocity. transverse wave speed of sound ,density Longitudinal wave propagation loss transverse wave speed of sound .

3. The method according to claim 1, characterized in that: In step one, The expression for normal modes in a general waveguide can be rewritten as an expression for the amplitude of normal modes and the exponent e in an interference structure: (1) The amplitude of the normal mode The expression is: (2) Write the normal wave interference amplitude based on the sound pressure field expression. expression: (3) According to formulas (1), (2), (3), (3), the interference amplitude of each modal normal wave is given.

4. The method according to claim 2, characterized in that: In step two, Based on the LOFAR spectrum of the sound intensity of the signal received by the hydrophone Perform a two-dimensional Fourier transform, as shown in equation (4): (4) in It is the radial velocity of the sound source's motion. It is the length of the observation period. It is the observation frequency band width. It is the frequency shift in the Fourier transform factor. It is the time delay in the Fourier transform factor.

5. The method according to claim 4, characterized in that: In step two, Spectral density in different angular coefficient directions in the frequency shift-time delay domain Amplitude distribution detector Written as: (5) Based on the maximum angle coefficient obtained by the detector in equation (5), determine the coordinates of the peak point closest to the origin in that direction, formed by adjacent modal normal waves. .

6. A system for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena, characterized in that: The system is used to perform the method for estimating the motion parameters of an elastic seabed acoustic source based on acoustic field interference phenomena as described in any one of claims 1 to 5; The system includes: a feature analysis module, a Fourier transform module, and a motion parameter estimation module; The feature analysis module is used to analyze the normal mode characteristics and interference structure characteristics of the elastic seabed waveguide environment. This elastic seabed waveguide environment includes a seawater layer, an elastic sedimentary layer, and a rocky substrate. The analysis identifies the normal modes that play a major role in the acoustic field interference structure within the current elastic seabed waveguide, and determines the order of the first normal mode carrying the main energy. And the order of the normal mode that carries the main energy. ; The Fourier transform module performs a two-dimensional Fourier transform of the time-frequency interference structure based on the analyzed interference structure characteristics, and extracts the spectral density peak value. The motion parameter estimation module estimates the motion parameters of the sound source based on the analyzed normal mode characteristics of the elastic seabed waveguide and the calculated coordinates of the peak points of the spectral density. Based on the coordinates of the peak spectral density points calculated by the Fourier transform module and the normal mode characteristics of the elastic seabed, expressions for estimating the motion parameters of the sound source are given: (6) The value is the interference invariant value calculated based on equation (7): It is the center angular frequency of the observation band. The coordinates of the peak points formed by adjacent modal normal waves; (7) in It is the difference between the eigenvalues ​​of adjacent modal normals.

7. The system according to claim 6, characterized in that: The elastic marine waveguide environment analyzed by the feature analysis module includes a seawater layer, an elastic sedimentary layer, and a rocky substrate. The seawater layer medium is an ideal fluid, at a depth of Its speed of sound is Its density It is independent of depth and distance; The parameters of the elastic deposition layer include the longitudinal wave velocity. transverse wave speed of sound ,density Longitudinal wave propagation loss transverse wave speed of sound With the depth of the sedimentary layer ; The parameters of the rocky basement include longitudinal wave velocity. transverse wave speed of sound ,density Longitudinal wave propagation loss transverse wave speed of sound .

8. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 5.

9. A computer-readable storage medium for storing computer instructions, characterized in that, When the computer instructions are executed by the processor, they implement the steps of the method according to any one of claims 1 to 5.

Images