A method, system, device and medium for quickly predicting the movement of underwater sound field convergence area caused by ocean mesoscale vortex
Through the acoustic field ray tracing theory and Gaussian vortex model, combined with satellite remote sensing and measured data, the location of the acoustic wave convergence zone under the influence of deep-sea vortex is quickly predicted, which solves the problem of change in the position of the acoustic wave convergence zone in the existing technology and improves the forecasting speed and accuracy.
Patent Information
- Application Number
- CN202210734264.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-06-27
- Publication Date
- 2025-09-05
- Estimated Expiration
- 2042-06-27
AI Technical Summary
The prior art lacks a fast forecast method for the location of acoustic wave convergence zones under the influence of mesoscale vortexes in the ocean, resulting in changes in the propagation path of the sound wave and changes in the position of the convergence zones, affecting underwater acoustic communication and marine biological detection.
The joint forecasting process of ocean vortex and sound field based on the acoustic field ray tracing theory is adopted, and the vortex is identified and tracked through satellite remote sensing data, the sound velocity field is reconstructed with actual measured data, and the characteristic parameters are extracted using the Gaussian vortex model to calculate the span of the whole sound line of the converging belt to quickly predict the location of the converging zone.
It realizes the rapid forecasting and convergence area in the deep-sea vortex frequency region, completes satellite remote sensing tracking, seawater temperature salt reconstruction and Gaussian vortex extraction, and improves the speed and accuracy of sound wave propagation forecasts.
Smart Images

Figure CN115292884B_ABST
Abstract
Description
Technical Field
[0001] The present invention belongs to the technical field of deep sea areas where mesoscale eddies frequently occur, and in particular relates to a method, system, device, and medium for rapidly predicting the position of underwater acoustic convergence zones caused by ocean mesoscale eddies. The method is used to rapidly predict the position of acoustic convergence zones passing through mesoscale eddies. Background Art
[0002] The global oceans are rich in dynamical phenomena. Ocean eddies, a common mesoscale dynamical phenomenon, are widely distributed within the depth range of the seawater mixing layer. Their horizontal spatial scales range from tens to hundreds of kilometers, and their temporal scales range from days to months. When acoustic methods are used in offshore resource development and scientific investigations, sound waves are inevitably affected by ocean eddies, resulting in direction-finding errors and even failure to locate underwater targets. The large number of mesoscale eddies in the deep sea cause spatiotemporal variations in hydrological environmental parameters, severely affecting the horizontal refraction and interlayer reflection of underwater sound waves. Mesoscale eddies alter the sound wave propagation path, causing drastic changes in sound propagation loss and shifting the location of convergence zones. The propagation of sound waves through ocean eddies has a significant impact on both human activities and marine animals. Ocean eddies not only affect underwater acoustic communication networks but also have an impact on the detection and communication of marine life.
[0003] The study of mesoscale eddies began in the 1930s and flourished in the 1970s. Hydroacousticians have long focused on underwater sound propagation, such as the location of acoustic energy convergence and the spatial and temporal distribution of sound intensity under the influence of deep-sea eddies. However, a suitable interface between ocean and hydroacoustic modeling and forecasting is currently lacking, resulting in slow resolution of acoustic fields in ocean environments. Therefore, there is an urgent need to establish a joint prediction process for ocean eddies and acoustic fields, enabling rapid prediction of the locations of acoustic convergence zones in eddy-prone waters. Summary of the Invention
[0004] This invention addresses the shifting of acoustic energy convergence locations under the influence of deep-sea mesoscale vortices. To address existing technical issues, it proposes a method, system, device, and medium for rapidly predicting the shifting of underwater acoustic convergence zones caused by ocean mesoscale vortices. This invention proposes a combined ocean vortex and acoustic field prediction process based on acoustic field ray tracing theory. This process tracks vortex locations, extracts Gaussian vortex characteristics of real vortices, reconstructs the spatial distribution of the sound velocity field, and develops an algorithm for solving the distance-dependent span of ambient sound rays. Ultimately, this method achieves the goal of rapidly predicting the shifting of acoustic convergence zones under the influence of vortices.
[0005] The present invention is achieved through the following technical solutions. The present invention proposes a method for rapidly predicting the position movement of the underwater sound field convergence zone caused by ocean mesoscale vortices, which specifically includes the following steps:
[0006] Step 1: Construct a mesoscale eddy tracking model: Use satellite remote sensing geostrophic velocity data to identify and track vortices, and determine the location of the vortex center, vortex boundary, and average radius using four constraints.
[0007] Step 2: Reconstructing measured data and satellite remote sensing seawater temperature and salinity data: A towed temperature chain and a drop-type temperature and depth gauge were used to observe the vortex. By analyzing the vortex's three-dimensional real-time measurement data, the temperature, salinity, and density of slices at different water depths were reconstructed. The sound speed was extended at depth using an empirical formula that shows how the sound speed changes with temperature, salinity, and depth.
[0008] Step 3: Using the least squares method to extract the characteristics of the vortex information based on the elliptical Gaussian vortex sound velocity model, the position and characteristic parameters of the vortex center are obtained, thereby processing the tracked vortex boundary;
[0009] Step 4: Calculate the span of the start and end sound lines of the convergence zone to make a rapid prediction.
[0010] Furthermore, the vortex boundary detection method is a detection algorithm for defining the vortex boundary by the contour lines of the stream function, the contour lines of relative vorticity, and the vortex boundary using the Okubo-Weiss criterion.
[0011] Furthermore, the four constraints are specifically:
[0012] (1) The geostrophic velocity v in the east-west direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center;
[0013] (2) The geostrophic velocity u in the north-south direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center;
[0014] (3) The velocity is at its minimum in the local area at the vortex center;
[0015] (4) The rotational characteristics of the mesoscale vortex velocity vector changes around the vortex center are consistent; the directions of two adjacent velocity vectors must be located in the same or two adjacent quadrants.
[0016] Furthermore, the empirical formula for the speed of sound is:
[0017] c=c0+(16.23+0.253T)z+(0.213-0.1T)z 2 +[0.016+0.0002(S-35)](S-35)Tz(1)
[0018] Where c0 = 1449.05 + 45.7T - 5.21T 2 +0.23T 3 +(1.333-0.126T+0.009T2 )(S-35), T is temperature, z is depth, and S is salinity.
[0019] Furthermore, the empirical formula for the speed of sound is limited to a temperature between 0 and 35°C and a salinity between 0‰ and 45‰.
[0020] Furthermore, in step three, the sound velocity distribution of the vortex is obtained based on the position and characteristic parameters of the vortex center:
[0021]
[0022] Where A c is the amplitude of the sound velocity change caused by the existence of the vortex, r0 and z0 are the positions of the vortex center, dr and dz are the characteristic parameters of the semi-major axis and semi-minor axis respectively; background sound velocity c0(z)=1500(1+0.0057[e -η -(1-η)]), which satisfies the Munk deep-sea sound velocity profile.
[0023] Furthermore, the step 4 is specifically as follows:
[0024] Due to the large size of the mesoscale vortex, the normal sound velocity gradient in the horizontal propagation direction of the sound ray is almost zero. When the horizontal refraction of the sound ray is ignored, it can be simplified to:
[0025]
[0026] Integrating the sound path in each direction gives:
[0027]
[0028] Further deducing the formula we get
[0029]
[0030]
[0031] When the critical angle |θ c |=arccos(c1 / c h ) will be reversed on the sea surface or seabed. A numerical method is used to solve the trajectory of the reversed sound line. When the sound source is located above the sound channel axis, the deep-sea sound velocity profile serves as the background environment and only two types of trajectories need to be calculated: the downward-bending sound line emitted at a 0° grazing angle and the upward-bending sound line with a small grazing angle and the sea surface inversion sound line. The position of the convergence zone can be determined.
[0032] The present invention proposes a rapid prediction system for the position movement of underwater sound field convergence area caused by ocean mesoscale vortices, the rapid prediction system comprising:
[0033] Model building module: Construct a mesoscale eddy tracking model: identify and track vortices using satellite remote sensing geostrophic velocity data, and determine the location of the vortex center, vortex boundary, and average radius using four constraints;
[0034] Reconstruction module: Reconstruction of measured data and satellite remote sensing seawater temperature and salinity data: Using a towed temperature chain and a drop-type temperature and depth gauge to observe vortices, the vortex is analyzed in real time in three dimensions to ultimately reconstruct the temperature, salinity, and density of slices at different water depths. The sound velocity is extended at depth using an empirical formula that expresses the change in sound velocity with temperature, salinity, and depth.
[0035] Feature extraction module: uses the least squares method to extract the features of vortex information based on the elliptical Gaussian vortex sound velocity model, obtains the position and characteristic parameters of the vortex center, and thus processes the tracked vortex boundary;
[0036] Rapid prediction module: Calculates the span of the start and end sound lines of the convergence zone to make rapid predictions.
[0037] The present invention proposes an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of a method for rapidly predicting the position movement of an underwater sound field convergence zone caused by an ocean mesoscale vortex are implemented.
[0038] The present invention proposes a computer-readable storage medium for storing computer instructions. When the computer instructions are executed by a processor, the steps of the method for quickly predicting the position movement of the underwater sound field convergence area caused by a mesoscale vortex in the ocean are realized.
[0039] The advantages of this invention are: it completely integrates satellite remote sensing tracking, seawater temperature and salinity reconstruction, Gaussian vortex extraction, solution of distance-dependent sound ray trajectories, and prediction of the start and end positions of the acoustic caustic zone under the influence of the vortex. By calculating the span of the start and end sound ray of the convergence zone, it can quickly predict the position movement of the convergence zone caused by deep-sea vortices. BRIEF DESCRIPTION OF THE DRAWINGS
[0040] Figure 1 This is a schematic diagram of the Northwest Pacific vortex tracking;
[0041] Figure 2 Schematic diagram of Gaussian vortex extraction;
[0042] Figure 3 Schematic diagram of ray tracing;
[0043] Figure 4 Schematic diagram of the sound speed profile of cyclonic vortex and anticyclonic vortex;
[0044] Figure 5 This is a schematic diagram of the acoustic ray trajectory at the beginning and end of the caustic zone under the influence of the vortex;
[0045] Figure 6 This is a flow chart of a method for rapidly predicting the position movement of underwater sound field convergence areas caused by ocean mesoscale vortices according to the present invention. DETAILED DESCRIPTION
[0046] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.
[0047] Combine Figures 1-6 The present invention proposes a method for rapidly predicting the position movement of the underwater sound field convergence area caused by the ocean mesoscale vortex, which specifically includes the following steps:
[0048] Step 1: Construct a mesoscale eddy tracking model: Use satellite remote sensing geostrophic velocity data to identify and track vortices, and determine the location of the vortex center, vortex boundary, and average radius using four constraints.
[0049] Step 2: Reconstruction of measured data and satellite remote sensing seawater temperature and salinity data: A towed conductivity temperature depth (CTD) and an expendable bathythermograph (XBT) were used to observe the vortex. By analyzing the three-dimensional real-time measurement data of the vortex, the temperature, salinity, and density of slices at different water depths were reconstructed. The sound speed was extended at depth using the empirical formula for the speed of sound given by Coppen, which shows how the speed of sound changes with temperature, salinity, and depth.
[0050] Step 3: Using the least squares method to extract the characteristics of the vortex information based on the elliptical Gaussian vortex sound velocity model, the position and characteristic parameters of the vortex center are obtained, thereby processing the tracked vortex boundary;
[0051] Step 4: Calculate the span of the start and end sound lines of the convergence zone to make a rapid prediction.
[0052] The vortex boundary detection method is a detection algorithm for the vortex boundary defined by the stream function contour line, the relative vorticity contour line and the Okubo-Weiss criterion.
[0053] The four constraints are specifically:
[0054] (1) The geostrophic velocity v in the east-west direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center;
[0055] (2) The geostrophic velocity u in the north-south direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center;
[0056] (3) The velocity is at its minimum in the local area at the vortex center;
[0057] (4) The rotational characteristics of the mesoscale vortex velocity vector changes around the vortex center are consistent; the directions of two adjacent velocity vectors must be located in the same or two adjacent quadrants.
[0058] The empirical formula for the speed of sound is:
[0059] c=c0+(16.23+0.253T)z+(0.213-0.1T)z 2 +[0.016+0.0002(S-35)](S-35)Tz(1)
[0060] Where c0 = 1449.05 + 45.7T - 5.21T 2 +0.23T 3 +(1.333-0.126T+0.009T 2 )(S-35), T is temperature, z is depth, and S is salinity.
[0061] The empirical formula for the speed of sound is limited to temperatures between 0 and 35°C and salinity between 0‰ and 45‰.
[0062] In step three, to model and study the underwater sound propagation issues through vortices, the elliptical Gaussian vortex sound velocity model has significant advantages. 80% of vortex types can be attributed to the Gaussian type. The sound velocity distribution of the vortex is obtained based on the position and characteristic parameters of the vortex center:
[0063]
[0064] Where A c is the amplitude of the sound velocity change caused by the existence of the vortex, r0 and z0 are the positions of the vortex center, dr and dz are the characteristic parameters of the semi-major axis and semi-minor axis respectively; background sound velocity c0(z)=1500(1+0.0057[e -η -(1-η)]), which satisfies the Munk deep-sea sound velocity profile.
[0065] The step 4 is specifically as follows:
[0066] Due to the large size of the mesoscale vortex, the normal sound velocity gradient in the horizontal propagation direction of the sound ray is almost zero. When the horizontal refraction of the sound ray is ignored, it can be simplified to:
[0067]
[0068] Integrating the sound path in each direction gives:
[0069]
[0070] Further deducing the formula we get
[0071]
[0072]
[0073] When the critical angle |θ c |=arccos(c1 / c h ) will be reversed on the sea surface or seabed. A numerical method is used to solve the trajectory of the reversed sound line. When the sound source is located above the sound channel axis, the deep-sea sound velocity profile serves as the background environment and only two types of trajectories need to be calculated: the downward-bending sound line emitted at a 0° grazing angle and the upward-bending sound line with a small grazing angle and the sea surface inversion sound line. The position of the convergence zone can be determined.
[0074] The present invention proposes a rapid prediction system for the position movement of underwater sound field convergence area caused by ocean mesoscale vortices, the rapid prediction system comprising:
[0075] Model building module: Construct a mesoscale eddy tracking model: identify and track vortices using satellite remote sensing geostrophic velocity data, and determine the location of the vortex center, vortex boundary, and average radius using four constraints;
[0076] Reconstruction module: Reconstruction of measured data and satellite remote sensing seawater temperature and salinity data: Using a towed temperature chain and a drop-type temperature and depth gauge to observe vortices, the vortex is analyzed in real time in three dimensions to ultimately reconstruct the temperature, salinity, and density of slices at different water depths. The sound velocity is extended at depth using an empirical formula that expresses the change in sound velocity with temperature, salinity, and depth.
[0077] Feature extraction module: uses the least squares method to extract the features of vortex information based on the elliptical Gaussian vortex sound velocity model, obtains the position and characteristic parameters of the vortex center, and thus processes the tracked vortex boundary;
[0078] Rapid prediction module: Calculates the span of the start and end sound lines of the convergence zone to make rapid predictions.
[0079] Example
[0080] by Figure 1 The northwestern Pacific Ocean is shown as an example. Two vortices with strong kinetic energy are selected. The colors in the figure represent the sea surface height anomalies, with blue representing areas below sea level and red representing areas above sea level. The black arrows represent the distribution of geostrophic current velocity on the sea surface.
[0081] Figure 2 The blue and red dashed lines represent the tracked vortex boundaries, while the green solid line represents the extracted Gaussian vortex boundary. Least squares methods were used to process the tracked vortex boundaries and extract the Gaussian vortex features of the actual vortex with an error of less than 25%. The radius of the cyclonic vortex is 79.42 km, with a vorticity of -20 m / s; the radius of the anticyclonic vortex is 69.26 km, with a vorticity of 25 m / s.
[0082] Ray tracing of sound ray trajectories in two-dimensional space is shown in the figure Figure 3 When the sound source is above the sound channel axis, the deep sea sound velocity profile as the background environment is as follows: Figure 4 shown.
[0083] Figure 5 As shown in Figure 1, under the influence of the mesoscale vortex, the inversion position of the initial and final caustics in the first convergence zone of the sound field will shift. It is only necessary to calculate two types of trajectories: the downward-bending sound line with a 0° grazing angle and the upward-bending sound line with a small grazing angle and a sea surface inversion. The final calculated convergence zone position movement position is shown in Table 1.
[0084] Table 1 Caustic zone position
[0085]
[0086] The present invention establishes a joint prediction model for ocean vortexes and acoustic fields, combines satellite remote sensing tracking, dynamic phenomenon extraction, and seawater temperature and salinity reconstruction, and realizes rapid prediction of the convergence zone position by calculating the start and end sound line span of the convergence zone. It completely unifies the prediction process of satellite remote sensing tracking, seawater temperature and salinity reconstruction, dynamic phenomenon extraction, and three-dimensional sound field propagation. It establishes a suitable interface for ocean and hydroacoustic modeling and prediction, uses the Gaussian vortex model to extract the characteristics of vortex information, and quickly predicts the start and end sound line span of the convergence zone. It tracks the vortex position, combines the measured data of the temperature and salinity profile, and proposes a method for reconstructing the hydrological environment parameters using the Gaussian vortex model; uses the extracted characteristic sound speed to solve the caustic line span of the sound field convergence zone under the distance-related environment, and simulates the sound energy distribution of the sound wave passing through the ocean vortex and the measured vortex. Compared with the background hydrological environment without vortexes, after the sound wave passes through the cyclonic vortex and anticyclonic vortex, the position of the convergence zone moves towards and away from the sound source and the sea surface respectively. It solves the problem of quickly predicting the location of the convergence zone in areas where deep-sea vortices frequently occur. In practical applications, its principle is easy to understand and highly operational.
[0087] The present invention proposes an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of a method for rapidly predicting the position movement of an underwater sound field convergence zone caused by an ocean mesoscale vortex are implemented.
[0088] The present invention proposes a computer-readable storage medium for storing computer instructions. When the computer instructions are executed by a processor, the steps of the method for quickly predicting the position movement of the underwater sound field convergence area caused by a mesoscale vortex in the ocean are realized.
[0089] The memory in the embodiments of the present application may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. The non-volatile memory may be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory may be a random access memory (RAM), which is used as an external cache. By way of example and not limitation, many forms of RAM are available, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), and direct RAM bus RAM (DR RAM). It should be noted that the memory of the methods described herein is intended to include, but is not limited to, these and any other suitable types of memory.
[0090] In the above embodiments, all or part of the embodiments may be implemented by software, hardware, firmware, or any combination thereof. When implemented using software, all or part of the embodiments may be implemented in the form of 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 the present application are generated. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via a wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) method. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated therein. The available medium may be a magnetic medium (eg, a floppy disk, a hard disk, a magnetic tape), an optical medium (eg, a high-density digital video disc (DVD)), or a semiconductor medium (eg, a solid state disc (SSD)).
[0091] During implementation, each step of the above method can be completed by an integrated logic circuit of the hardware in the processor or by instructions in the form of software. The steps of the method disclosed in conjunction with the embodiments of the present application can be directly embodied as being executed by a hardware processor, or can be executed by a combination of hardware and software modules in the processor. The software module can be located in a storage medium mature in the art such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in conjunction with its hardware. To avoid repetition, it will not be described in detail here.
[0092] It should be noted that the processor in the embodiments of the present application can be an integrated circuit chip with signal processing capabilities. During implementation, each step of the above method embodiment can be completed by an integrated logic circuit of the hardware in the processor or by instructions in the form of software. The above 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 device, a discrete gate or transistor logic device, or a discrete hardware component. The various methods, steps, and logic block diagrams disclosed in the embodiments of the present application can be implemented or executed. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in the embodiments of the present application can be directly embodied as being executed by a hardware decoding processor, or can be executed by a combination of hardware and software modules in the decoding processor. The software module can be located in a storage medium mature in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, or electrically erasable programmable memory, registers, etc. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.
[0093] The above is a detailed introduction to the method, system, equipment and medium for quickly predicting the position movement of the underwater sound field convergence area caused by the ocean mesoscale vortex proposed in the present invention. Specific examples are used in this article to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for general technical personnel in this field, according to the ideas of the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.
Claims
1. A method for rapidly predicting the position movement of underwater sound field convergence area caused by ocean mesoscale eddies, characterized in that: The specific steps include: Step 1: Construct a mesoscale eddy tracking model: Use satellite remote sensing geostrophic velocity data to identify and track vortices, and determine the location of the vortex center, vortex boundary, and average radius using four constraints. Step 2: Reconstructing measured data and satellite remote sensing seawater temperature and salinity data: A towed temperature chain and a drop-type temperature and depth gauge were used to observe the vortex. By analyzing the vortex's three-dimensional real-time measurement data, the temperature, salinity, and density of slices at different water depths were reconstructed. The sound speed was extended at depth using an empirical formula that shows how the sound speed changes with temperature, salinity, and depth. Step 3: Using the least squares method to extract the characteristics of the vortex information based on the elliptical Gaussian vortex sound velocity model, the position and characteristic parameters of the vortex center are obtained, thereby processing the tracked vortex boundary; Step 4: Calculate the span of the start and end sound lines of the convergence zone to make a rapid prediction; The four constraints are specifically: (1) The geostrophic velocity v in the east-west direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center; (2) The geostrophic velocity u in the north-south direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center; (3) The velocity in the local area at the vortex center is the minimum extreme value; (4) The rotation of the velocity vector of the mesoscale vortex around the vortex center is consistent; the directions of two adjacent velocity vectors must be located in the same or two adjacent quadrants.
2. The method according to claim 1, characterized in that The vortex boundary detection method is a detection algorithm for the vortex boundary defined by the stream function contour line, the relative vorticity contour line and the Okubo-Weiss criterion.
3. The method according to claim 1, characterized in that The empirical formula for the speed of sound is: (1) in, , T is temperature, z is depth, and S is salinity.
4. The method according to claim 3, characterized in that The empirical formula for the speed of sound is limited to temperatures between 0 and 35°C and salinity between 0‰ and 45‰.
5. The method according to claim 4, characterized in that In step 3, the sound velocity distribution of the vortex is obtained based on the position and characteristic parameters of the vortex center: (2) In the formula is the amplitude of the sound speed change caused by the existence of the vortex, and is the position of the vortex center, and are the characteristic parameters of the semi-major axis and semi-minor axis respectively; background sound speed , meeting the Munk deep-sea sound velocity profile.
6. The method according to claim 5, characterized in that The step 4 is specifically as follows: Due to the large size of the mesoscale vortex, the normal sound velocity gradient in the horizontal propagation direction of the sound ray is almost zero. When the horizontal refraction of the sound ray is ignored, it can be simplified to: (3) Integrating the sound path in each direction gives: (4) Further deducing the formula we get (5) (6) When the critical angle The outgoing sound line will be reversed at the sea surface or seabed. The numerical method is used to solve the trajectory of the reversed sound line. When the sound source is located above the sound channel axis, the deep-sea sound velocity profile is used as the background environment. Only two types of trajectories need to be calculated: the downward-bending sound line at a 0° grazing angle and the upward-bending sound line at a small grazing angle and the sea surface inversion. The location of the convergence zone can be determined.
7. A rapid prediction system for the movement of underwater sound field convergence area caused by ocean mesoscale eddies, characterized by: The rapid forecasting system comprises: Model building module: Construct a mesoscale eddy tracking model: identify and track vortices using satellite remote sensing geostrophic velocity data, and determine the location of the vortex center, vortex boundary, and average radius using four constraints; Reconstruction module: Reconstruction of measured data and satellite remote sensing seawater temperature and salinity data: Using a towed temperature chain and a drop-type temperature and depth gauge to observe vortices, the vortex is analyzed in real time in three dimensions to ultimately reconstruct the temperature, salinity, and density of slices at different water depths. The sound velocity is extended at depth using an empirical formula that expresses the change in sound velocity with temperature, salinity, and depth. Feature extraction module: uses the least squares method to extract the features of vortex information based on the elliptical Gaussian vortex sound velocity model, obtains the position and characteristic parameters of the vortex center, and thus processes the tracked vortex boundary; Rapid prediction module: calculates the span of the start and end sound lines of the convergence zone to make rapid predictions; The four constraints are specifically: (1) The geostrophic velocity v in the east-west direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center; (2) The geostrophic velocity u in the north-south direction has opposite signs on both sides of the vortex center, and its magnitude gradually increases as it moves away from the vortex center; (3) The velocity in the local area at the vortex center is the minimum extreme value; (4) The rotation of the velocity vector of the mesoscale vortex around the vortex center is consistent; the directions of two adjacent velocity vectors must be located in the same or two adjacent quadrants.
8. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.
9. A computer-readable storage medium for storing computer instructions, characterized in that: When the computer instructions are executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.
Citation Information
Patent Citations
Method for obtaining ocean mesoscale vortex wave impedance based on deep learning and processing terminal
CN113312830A
Method for constructing mesoscale vortex underwater sound velocity disturbance general structure
CN113984178A