Method and system for constructing the acoustic channel of the ocean mesoscale eddy affecting the surface and subsurface

Abstract

The present invention discloses a method and system for constructing a surface and subsurface sound channel affected by a mesoscale eddy in the ocean. The method comprises obtaining hydrological parameters of a specified area within a specified time; calculating the Argo sound velocity profile and the sound velocity background field of the specified area and performing quality control to obtain an Argo sound velocity profile data set; extracting the characteristics of the mesoscale eddy, and generating a mesoscale eddy disturbance sound velocity field by temporal and spatial matching and merging of the Argo sound velocity profile data set and the mesoscale eddy trajectory data set; interpolating the underwater sound field and the sound velocity background field of the mesoscale eddy to separate the sound velocity background field and extract the underwater sound velocity disturbance structure of the mesoscale eddy, and obtaining the abnormal data of sound propagation loss under the influence of the mesoscale eddy in combination with the Bellhop acoustic toolbox to extract the characteristics of the surface and subsurface sound channels. The present invention aims to obtain the influence of the mesoscale eddy on the surface sound channel and accurately obtain the variability characteristics of the surface and subsurface sound channels corresponding to the mesoscale eddy.

Description

Technical Field

[0001] The present invention relates to the technical field of constructing underwater surface and subsurface sound channels of ocean mesoscale vortices, and in particular to a method and system for constructing surface and subsurface sound channels affected by ocean mesoscale vortices. Background Art

[0002] Sound waves are the only energy form known to mankind that is suitable for long-distance propagation in seawater. Marine acoustics is an irreplaceable technical means for the development and utilization of the ocean. Its role in underwater detection, navigation, positioning, communication and other fields is similar to that of radio waves in the air. In civilian use, water sound can be used for seabed topography mapping, natural resource exploration, etc.; in military use, water sound can be used for underwater communication, water sound listening and detection, water sound interference, target positioning and navigation, etc., providing necessary guarantees for underwater military activities. The propagation characteristics of water sound are closely related to the characteristics of the marine environment. The surface sound channel mainly exists in a mixed layer of a certain thickness under the sea surface, which is usually located within a depth of 100 m below the sea surface. The sound speed in the layer is distributed in a weak positive gradient. The sound rays with small grazing angles near the sea surface are reversed due to refraction in the mixed layer and reflected at the sea surface. This process is repeated continuously, so that the sound energy is confined to the surface of the ocean, thus forming a waveguide-type propagation. Therefore, this sound channel is called the surface sound channel (or mixed layer sound channel). The subsurface sound channel mainly exists at a depth of hundreds of meters below the sea surface. The formation mechanism of the subsurface sound channel is similar to that of the deep sea sound channel. Due to refraction, the sound rays are bounded in a water layer of a certain thickness above and below the sound channel axis at the depth where the sound velocity is the minimum, forming a convergence area of ​​sound energy. The characteristic of this sound channel axis is that the sound velocity in the upper layer is distributed in a negative gradient with the increase of pressure, while the lower layer is distributed in a negative gradient.

[0003] Mesoscale eddies are circulations with high energy, high velocity, and time-varying currents. The spatial horizontal scale is hundreds of kilometers, and its vertical structure can extend to one thousand meters underwater. It gathers the kinetic energy of seawater over a horizontal distance of tens to hundreds of kilometers. Its significant nonlinear characteristics enable it to transport heat and salt in the ocean. Studies have shown that the existence of mesoscale eddies affects the propagation of sound waves in the surface sound channel by changing the depth and thickness of the mixing layer. Due to the large vertical temperature difference in the ocean, the sound velocity structure of the upper ocean is significantly affected by mesoscale eddies. Therefore, hydroacoustic scholars have conducted in-depth research on the influence of mesoscale eddies on sound wave propagation. In terms of physical mechanism, the huge kinetic energy carried by mesoscale eddies can fully mix the seawater in the mesoscale eddy region, deepening the mixing layer. The different rotation directions can produce cold cores and warm cores in the center of the mesoscale eddies, and their disturbances also weaken as the distance from the center of the mesoscale eddies increases, thus producing different acoustic effects.

[0004] Research has shown that the influence of mesoscale eddies on the underwater sound velocity field can extend to a depth of 1000 m underwater, thereby changing the distribution of acoustic properties in the upper ocean. This depth covers the entire range of submarine activities and has a significant impact on military activities such as submarine concealment and sonar detection. On the other hand, because mesoscale eddies change the distribution pattern of acoustic properties in the upper ocean, surface and subsurface channels can be formed under certain conditions, thereby achieving long-distance propagation of sound waves.

[0005] Therefore, the impact of mesoscale vortices on the surface acoustic channel has high research value and is challenging and difficult. It is of great significance to improve the utilization of the acoustic effects of mesoscale vortices in the ocean and enhance the ability to ensure the underwater battlefield environment. Summary of the invention

[0006] Technical problem to be solved by the present invention: In view of the above-mentioned problems in the prior art, a method and system for constructing the surface and subsurface sound channels affected by mesoscale eddies in the ocean are provided. The present invention aims to obtain the influence of mesoscale eddies on the surface sound channels and accurately obtain the variability characteristics of the surface and subsurface sound channels corresponding to the mesoscale eddies.

[0007] In order to solve the above technical problems, the technical solution adopted by the present invention is:

[0008] A method for constructing an ocean mesoscale eddy impact surface and subsurface sound channel comprises the following steps:

[0009] S1, obtains the hydrological parameters of the specified area within the specified time according to the input Argo float dataset and SODA3.4.2 climate state dataset;

[0010] S2, the hydrological parameters of the Argo float data set are used to calculate the Argo sound velocity profile of the specified area using the Mackenzie sound velocity empirical formula, and the hydrological parameters of the SODA3.4.2 climate state data set are used to calculate the sound velocity background field of the specified area using the Mackenzie sound velocity empirical formula. The Argo sound velocity profile calculated by the Argo float data set is quality controlled to remove bad data to obtain the Argo sound velocity profile data set;

[0011] S3, extracting the characteristics of the mesoscale eddy according to the input mesoscale eddy trajectory dataset;

[0012] S4, based on the characteristics of the mesoscale vortex and the position and time of the Argo float, the temporal and spatial matching and merging of the Argo sound velocity profile data set and the mesoscale vortex trajectory data set are completed, and vertical and horizontal smoothing is performed to obtain the underwater sound velocity characteristics of the mesoscale vortex and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex;

[0013] S5, interpolating the underwater sound field and the sound velocity background field of the mesoscale vortex formed by the characteristics of the mesoscale vortex to the specified depth respectively, subtracting the sound velocity background field from the underwater sound field of the mesoscale vortex to separate the sound velocity background field, and then extracting the sound velocity anomaly characteristics in four seasons of the year, thereby obtaining the underwater sound velocity disturbance structure of the mesoscale vortex;

[0014] S6, input the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, input the sound velocity background field into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field; subtract the climatic underwater sound pressure field from the mesoscale vortex underwater sound pressure field to obtain the abnormal data of sound propagation loss under the influence of the mesoscale vortex;

[0015] S7, combining the abnormal data of sound propagation loss under the influence of mesoscale eddies and the underwater sound velocity disturbance structure of mesoscale eddies, extracts the characteristics of the surface and subsurface sound channels affected by mesoscale eddies in the ocean.

[0016] Optionally, the hydrological parameters in step S1 include temperature, salinity and depth of seawater.

[0017] Optionally, the function expression of the Mackenzie sound speed empirical formula in step S2 is:

[0018] ,

[0019] In the above formula, represents the speed of sound, The temperature of the sea water, The salinity of seawater. Indicates the depth of sea water.

[0020] Optionally, in step S2, quality control is performed on the Argo sound speed profile calculated by the Argo float data set to remove bad data to obtain the Argo sound speed profile data set, including: calculating the sound speed standard deviation of the Argo sound speed profile calculated by the Argo float data set, and deleting the Argo sound speed profile outside three times the sound speed standard deviation, thereby achieving quality control to remove bad data to obtain the Argo sound speed profile data set.

[0021] Optionally, step S4 includes:

[0022] S4.1, based on the characteristics of the mesoscale vortex including the central coordinates, appearance time, radius, polarity and amplitude of the mesoscale vortex, the influence range and time of each mesoscale vortex are calculated using the central coordinates, appearance time and radius of the mesoscale vortex, and then the Argo sound velocity profile dataset is spatially and temporally matched with the mesoscale vortex dataset based on the position and time recorded by the Argo float;

[0023] S4.2, based on the center coordinates of the mesoscale vortex and the coordinates of the Argo float within twice its radius, obtain the radial normalized coordinates of the underwater sound speed profile of each mesoscale vortex relative to the center of the mesoscale vortex, and then integrate the underwater sound speed profiles of mesoscale vortices with the same characteristics into the same mesoscale vortex based on the normalized coordinates;

[0024] S4.3, set the sliding window size multiple times, find the optimal sliding window size according to the smoothing results of each vertical and horizontal smoothing process, and use the optimal sliding window sizes in the vertical and horizontal directions to perform vertical and horizontal smoothing processes respectively, so as to obtain the underwater sound velocity characteristics of mesoscale vortices with different amplitudes under different sound velocity background fields, and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of mesoscale vortices.

[0025] Optionally, step S4.2 includes:

[0026] S4.2.1, obtain the Argo sound velocity profile for each mesoscale eddy pair based on the coordinates of the center coordinates of the mesoscale eddy and the coordinates of the Argo float within the range of twice its radius on the same day of observation;

[0027] S4.2.2, project all paired Argo sound velocity profiles into a coordinate system with the vortex center of the corresponding mesoscale vortex as the origin Below, among them To fall in The first Argo sound velocity profiles and vortex centers The distance is the z-axis coordinate;

[0028] S4.2.3, calculate the The radial normalized distance of each Argo sound speed profile within a mesoscale vortex:

[0029] ,

[0030] In the above formula, is the radial normalized distance of the kth Argo sound velocity profile, , For the The radius of the mesoscale vortex is , To fall to The number of mesoscale eddies;

[0031] S4.2.4, the radial normalized distance of each Argo sound velocity profile Arrange them in the vortex center coordinate system according to their size, and then according to the radius of the synthetic vortex integrated into the same mesoscale eddy.

[0032] Optionally, in step S6, the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex is input into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, and the sound velocity background field is input into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field, including judging whether there is a surface sound channel in the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex, and if the mesoscale vortex underwater sound velocity characteristics constitute the surface sound channel, the surface sound channel is calculated. If there is a surface sound channel in the mesoscale eddy disturbance sound velocity field formed by the underwater sound velocity characteristics of the mesoscale vortex, the sound source is set at the sea surface. Otherwise, if there is a subsurface sound channel in the mesoscale eddy disturbance sound velocity field formed by the underwater sound velocity characteristics of the mesoscale vortex, the sound source is set at the depth of the subsurface sound channel axis, and the sound ray simulation is performed based on the operation type R and the sound propagation loss is calculated based on the operation type C in the ray model of the Bellhop acoustic toolbox, so as to obtain the mesoscale vortex underwater sound channel sound ray model and the mesoscale vortex underwater sound pressure field or the climatic underwater sound pressure field.

[0033] In addition, the present invention also provides a system for constructing the sound channels of the surface and subsurface layers affected by mesoscale eddies in the ocean, comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to execute the method for constructing the sound channels of the surface and subsurface layers affected by mesoscale eddies in the ocean.

[0034] In addition, the present invention also provides a computer-readable storage medium, which stores a computer program or instruction, and the computer program or instruction is programmed or configured to execute the method of constructing the ocean mesoscale eddy-affected surface and subsurface sound channels through a processor.

[0035] In addition, the present invention also provides a computer program product, including a computer program or instructions, which are programmed or configured to execute the method of constructing the ocean mesoscale eddy-affected surface and subsurface sound channels through a processor.

[0036] Compared with the prior art, the present invention mainly has the following advantages: the present invention comprises obtaining the hydrological parameters of a specified area within a specified time according to the input Argo float data set and the SODA3.4.2 climate state data set; calculating the hydrological parameters by using the Mackenzie sound speed empirical formula to obtain the corresponding sound speed profile, and filtering out the profile containing bad points by quality control to obtain the Argo sound speed profile data set; extracting the characteristics of the mesoscale vortex according to the input mesoscale vortex trajectory data set, completing the time-space matching of the Argo sound speed profile data set and the mesoscale vortex trajectory data set based on the characteristics of the mesoscale vortex and the Argo float position and time, and obtaining the underwater sound speed characteristics of the mesoscale vortex; the present invention separates the sound speed background field and eliminates interference by sliding smoothing to obtain the accurate sound speed disturbance characteristics corresponding to the mesoscale vortex; based on the underwater sound speed field of the mesoscale vortex, the present invention uses the Bellhop acoustic ray model to construct the characteristics of the surface and subsurface sound channels under the occurrence of the ocean mesoscale vortex, provides a reference for the effect of the mesoscale vortex on the surface and subsurface, and provides a theoretical support for studying the variability characteristics of the underwater surface and subsurface sound channels of different mesoscale vortices. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 Schematic diagram of the process of the embodiment of the present invention.

[0038] Figure 2 1 is the average amplitude distribution diagram of the anticyclonic vortex found in the Kuroshio extension zone during the study period in the embodiment of the present invention, wherein (a) is the average amplitude distribution diagram in spring, (b) is the average amplitude distribution diagram in summer, (c) is the average amplitude distribution diagram in autumn, and (d) is the average amplitude distribution diagram in winter.

[0039] Figure 3 3 is a density probability diagram of the radius and amplitude of all mesoscale vortices found in the Kuroshio extension area during the study in the embodiment of the present invention.

[0040] Figure 4 Schematic diagram of the matching synthesis of mesoscale vortices and Argo sound velocity profile data in an embodiment of the present invention.

[0041] Figure 5 Schematic diagram of the synthetic mesoscale vortex sound velocity structure of all amplitudes in the Kuroshio extension area in each season in the embodiment of the present invention, where (a) to (d) are the underwater sound velocity structures of anticyclonic vortices in each season, (e) to (h) are the underwater sound velocity structures of cyclonic vortices in each season, is the radius of the mesoscale vortex. The shadow in the figure represents the sound speed anomaly, and the black solid line represents the sound speed contour.

[0042] Figure 6: The channel sound pressure field formed by the 0.4 m amplitude mesoscale vortex in the Kuroshio extension area in each season in the embodiment of the present invention, where (a) is the channel sound pressure field of the surface channel in spring with a sound source depth of 10 m, (b) is the channel sound pressure field of the subsurface channel in summer with a sound source depth of 90 m, and (c) is the channel sound pressure field of the double channel in autumn with a sound source depth of 10 (120) m (when the sound source is set at 10 m, it is the surface channel, and when the sound source is set at 120 m, it is the subsurface channel). m), (d) is the sound pressure field of the surface sound channel in winter and the sound source depth of 10m, (e) the sound pressure field of the subsurface sound channel in spring and the sound source depth of 496m, (f) the sound pressure field of the subsurface sound channel in summer and the sound source depth of 481m, (g) the sound pressure field of the subsurface sound channel in autumn and the sound source depth of 498m, (h) the sound pressure field of the subsurface sound channel in winter and the sound source depth of 514m. DETAILED DESCRIPTION

[0043] The present invention aims to extract the underwater sound velocity structure of the mesoscale vortex by means of synthetic analysis, and use the Bellhop acoustic toolbox to simulate the sound line and calculate the sound pressure field to obtain the accurate surface and subsurface sound channel variability characteristics corresponding to the mesoscale vortex. The specific implementation methods of the present invention will be further described in detail below in conjunction with the accompanying drawings.

[0044] like Figure 1 As shown, the method for constructing the ocean mesoscale eddy impact surface and subsurface sound channels in this embodiment includes the following steps:

[0045] S1, obtains the hydrological parameters of the specified area within the specified time according to the input Argo float dataset and SODA3.4.2 climate state dataset;

[0046] S2, the hydrological parameters of the Argo float data set are used to calculate the Argo sound velocity profile of the specified area using the Mackenzie sound velocity empirical formula, and the hydrological parameters of the SODA3.4.2 climate state data set are used to calculate the sound velocity background field of the specified area using the Mackenzie sound velocity empirical formula. The Argo sound velocity profile calculated by the Argo float data set is quality controlled to remove bad data to obtain the Argo sound velocity profile data set;

[0047] S3, extracting the characteristics of mesoscale eddies based on the input META2.0 mesoscale eddy trajectory dataset;

[0048] S4, based on the characteristics of the mesoscale vortex and the position and time of the Argo float, the temporal and spatial matching and merging of the Argo sound velocity profile data set and the mesoscale vortex trajectory data set are completed, and vertical and horizontal smoothing is performed to obtain the underwater sound velocity characteristics of the mesoscale vortex and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex;

[0049] S5, interpolating the underwater sound field and the sound velocity background field of the mesoscale vortex formed by the characteristics of the mesoscale vortex to the specified depth respectively, subtracting the sound velocity background field from the underwater sound field of the mesoscale vortex to separate the sound velocity background field, and then extracting the sound velocity anomaly characteristics in four seasons of the year, thereby obtaining the underwater sound velocity disturbance structure of the mesoscale vortex;

[0050] S6, input the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, input the sound velocity background field into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field; subtract the climatic underwater sound pressure field from the mesoscale vortex underwater sound pressure field to obtain the abnormal data of sound propagation loss under the influence of the mesoscale vortex;

[0051] S7, combining the abnormal data of sound propagation loss under the influence of mesoscale eddies and the underwater sound velocity disturbance structure of mesoscale eddies, extracts the characteristics of the surface and subsurface sound channels affected by mesoscale eddies in the ocean.

[0052] The method for constructing the surface and subsurface sound channels of the ocean mesoscale eddies in this embodiment mainly includes seven steps, namely, data preprocessing (steps S1 to S3), matching synthesis (step S4), and feature extraction and analysis (steps S5 to S7).

[0053] In order to obtain the underwater sound velocity structure of the mesoscale vortex with the specified amplitude in the studied area and the studied time, it is necessary to use seawater temperature, seawater salinity and observation depth as input parameters of the Mackenzie sound velocity empirical formula. In step S1 of this embodiment, when obtaining the hydrological parameters of the specified area within the specified time according to the input Argo buoy data set and the SODA3.4.2 climate state data set, the hydrological parameters include the temperature, salinity and depth of the seawater, and the temperature, salinity and depth data of the seawater are screened out according to the input Argo buoy data set and the SODA3.4.2 climate state data set. The temperature, salinity and depth data can be directly obtained from the Argo buoy data set and the SODA3.4.2 climate state data set. The Argo buoy data set used in this embodiment is a data set released by the Coriolis Center after automatic quality control and processing, and the observation data of the Argo buoy recorded in the data set from 1998-05-03 to 2021-01-11 are selected, which mainly includes seawater pressure, temperature, salinity, conductivity, buoy number, coordinates and other related measurement information. This example selects the SODA3.4.2 climate state data from January 2000 to December 2019, which can provide three-dimensional monthly average ocean state variables of 50 vertical depth layers mapped to a regular 0.5° Mercator horizontal grid.

[0054] The functional expression of the Mackenzie sound velocity empirical formula in step S2 of this embodiment is:

[0055] ,

[0056] In the above formula, represents the speed of sound, The temperature of the sea water, The salinity of seawater. Indicates the depth of sea water.

[0057] In step S2 of this embodiment, the Argo sound velocity profile calculated by the Argo float data set is quality controlled to remove bad data to obtain the Argo sound velocity profile data set, including: calculating the sound velocity standard deviation of the Argo sound velocity profile calculated by the Argo float data set, and deleting the Argo sound velocity profile beyond three times the sound velocity standard deviation, so as to achieve quality control to remove bad data to obtain the Argo sound velocity profile data set. In addition, as an optional implementation, the quality control to remove bad data in this embodiment also includes: removing the Argo sound velocity profile with a temperature and salinity record of nan (not a number) at a certain depth (the temperature and salinity record of nan at a certain depth indicates that the Argo sound velocity profile has missing data at the depth, so it is considered that the Argo sound velocity profile has an erroneous record and should be removed), and removing the Argo sound velocity profile that does not cover the recording depth of 5 to 1000 m.

[0058] In step S3 of this embodiment, the characteristics of the mesoscale eddy are extracted according to the input META2.0 mesoscale eddy trajectory dataset. In addition, other mesoscale eddy trajectory datasets can also be used as needed. When the META2.0 mesoscale eddy trajectory dataset extracts the characteristics of the mesoscale eddy, it includes loading the center coordinates, appearance time, polarity, amplitude and radius of the mesoscale eddy. Mesoscale eddy trajectories atlas version 2.0 (META 2.0) is mainly used for special research on ocean mesoscale eddies. The original data of the product comes from the daily altimeter sampling of the two satellites from January 1, 1993 to March 7, 2020, and is reanalyzed and processed according to the mesoscale eddy identification and tracking scheme provided by Chelton. The META2.0 mesoscale eddy dataset contains the sea surface characteristics of mesoscale eddies in the global ocean from 1993-01-01 to 2020-03-07, mainly including the amplitude of the mesoscale eddy (the height difference between the SLA extreme value in the mesoscale eddy and the SLA that defines the contour of the mesoscale eddy circumference), the type of mesoscale eddy (cyclone is -1, anticyclone is +1), the coordinates of the center of the mesoscale eddy (latitude and longitude), the observation sequence number (the starting date of the first detection of the mesoscale eddy), the maximum average linear velocity, radius, observation date, trajectory identification number, etc. Figure 2 As shown in the figure, the large value area of ​​the anticyclonic vortex radius appears in 143°E~168°E, 34°N~39°N throughout the year, which is about 3° northward of the large value area of ​​the cyclonic vortex amplitude. The entire area is similar to the "shooting star" shape of the large value area of ​​the cyclonic vortex amplitude. The amplitude of the mesoscale vortex in the region is generally above 0.2 m, while most of the other regions are only below 0.1 m. In spring, the larger amplitude mesoscale vortices appear in the area of ​​145°E~148°E, 35°N~37°N, with a maximum value of 0.37 m. In summer, the area of ​​large amplitude of mesoscale vortices has no obvious change compared with spring, and the largest amplitude of mesoscale vortices appears in 146°E~148°E, 36°N~37°N, with a value of 0.38 m, which is close to that in spring. In autumn, the area of ​​large amplitude has no obvious change, and the largest amplitude of mesoscale vortices is 0.37 m, which appears in 147°E~148°E, 36°N~37°N. In winter, the area of ​​large amplitude has basically no change compared with autumn, and the larger amplitude of mesoscale vortices appears in 143°E~149°E, 35°N~37°N, with a maximum value of 0.36 m. In step S3, when the mesoscale eddy characteristics are statistically analyzed, it is found that the radius of the mesoscale eddy in the Kuroshio extension area is concentrated in 50-150 km, and the amplitude is concentrated in 0.05-0.1 m and 0.4 m. Figure 3 shown.

[0059] Step S4 of this embodiment includes:

[0060] S4.1, based on the characteristics of the mesoscale vortex including the central coordinates, appearance time, radius, polarity and amplitude of the mesoscale vortex, the influence range and time of each mesoscale vortex are calculated using the central coordinates, appearance time and radius of the mesoscale vortex, and then the Argo sound velocity profile dataset is spatially and temporally matched with the mesoscale vortex dataset based on the position and time recorded by the Argo float;

[0061] S4.2, based on the center coordinates of the mesoscale vortex and the coordinates of the Argo float within twice its radius, obtain the radial normalized coordinates of the underwater sound speed profile of each mesoscale vortex relative to the center of the mesoscale vortex, and then integrate the underwater sound speed profiles of mesoscale vortices with the same characteristics into the same mesoscale vortex based on the normalized coordinates;

[0062] S4.3, set the sliding window size multiple times, find the optimal sliding window size according to the smoothing results of each vertical and horizontal smoothing process, and use the optimal sliding window sizes in the vertical and horizontal directions to perform vertical and horizontal smoothing processes respectively, so as to obtain the underwater sound velocity characteristics of mesoscale vortices with different amplitudes under different sound velocity background fields, and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of mesoscale vortices.

[0063] The synthetic analysis method can utilize a large amount of asynchronous hydrological observation profile data and project it into a coordinate system with the center of the mesoscale vortex as the origin, thereby forming a mesoscale vortex structure image under statistical average. At this time, the size of the synthesized mesoscale vortex is dimensionless. Step S4.2 of this embodiment includes:

[0064] S4.2.1, obtain the Argo sound velocity profile for each mesoscale eddy pair based on the coordinates of the center coordinates of the mesoscale eddy and the coordinates of the Argo float within the range of twice its radius on the same day of observation;

[0065] S4.2.2, project all paired Argo sound velocity profiles into a coordinate system with the vortex center of the corresponding mesoscale vortex as the origin Below, among them To fall in The first Argo sound velocity profiles and vortex centers The distance is the z-axis coordinate;

[0066] S4.2.3, calculate the The radial normalized distance of each Argo sound speed profile within a mesoscale vortex:

[0067] ,

[0068] In the above formula, is the radial normalized distance of the kth Argo sound velocity profile, , For the The radius of the mesoscale vortex is , To fall to The number of mesoscale eddies;

[0069] S4.2.4, the radial normalized distance of each Argo sound velocity profile Arrange them in the vortex center coordinate system according to their size, and then according to the radius of the synthetic vortex integrated into the same mesoscale eddy.

[0070] The specific steps of step S4.2 can be divided into the following steps: (1) Pairing of Argo sound speed profiles with mesoscale vortices. Find all Argo sound speed profiles that fall into the selected area, and for each Argo sound speed profile, find the mesoscale vortex center with the shortest time interval and the closest spatial distance to pair them. The standard set in this embodiment for pairing the mesoscale vortex with the Argo sound speed profile is that the observation time is on the same day and the Argo falls within 2 times the radius of the mesoscale vortex. (2) Project the Argo sound speed profile onto the normalized vortex center coordinate system That is, the normalization of the radial distance. Figure 4 As shown in the figure, firstly, all paired Argo sound velocity profiles are projected onto the coordinate system of the corresponding mesoscale vortex with the center of the mesoscale vortex as the origin. Below, among them It falls on The first Argo sound velocity profiles and mesoscale vortex centers The distance, and then the radius of the mesoscale eddy Normalized Argo sound velocity profile and the distance from the center of the mesoscale vortex , and obtain the normalized distance of the Argo sound velocity profile ,in ( is the total number of Argo sound velocity profiles screened out), and is used as the radial coordinate during synthesis to eliminate the influence of the horizontal size difference of the mesoscale vortex. Here, the normalized distance The range is After obtaining the normalized radial distance of each Argo sound velocity profile, they can be arranged in the vortex center coordinate system according to their size. The radius of the synthetic mesoscale vortex is defined as Finally, the sliding window size is set multiple times, and the optimal sliding window size is found according to each smoothing result, and the result obtained by smoothing the optimal sliding window is used as the final result. In this embodiment, the optimal vertical sliding window is 1 and the optimal horizontal sliding window is 1000, thereby constructing the characteristics of the underwater two-dimensional sound velocity field of the mesoscale vortex. Figure 4 As shown in the figure, vortex 1 to vortex 3 merge to form a composite vortex, r1 is the radial distance from the Argo sound velocity profile matching vortex 1 to the center of the vortex, 2R1 is twice the radius of vortex 1, E1 is the vortex center of vortex 1, r1 * is the normalized radial distance of r1 (r1 * =r1 / R1), r2 is the radial distance from the center of the Argo sound velocity profile matching vortex 2 to the center of the vortex, 2R2 is twice the radius of vortex 2, E2 is the vortex center of vortex 2, r2 * is the normalized radial distance of r2 (r2* =r2 / R2), r3 is the radial distance from the Argo sound velocity profile matching vortex 3 to the center of the vortex, 2R3 is twice the radius of vortex 3, E3 is the vortex center of vortex 3, r3 * is the normalized radial distance of r3 (r3 * =r3 / R3). 2R is twice the radius of the synthetic vortex.

[0071] In step S5 of this embodiment, the underwater sound field of the mesoscale vortex and the sound velocity background field formed by the characteristics of the mesoscale vortex are interpolated to the specified depth, specifically to the depth of 5:1:1000 m, and the sound velocity background field is separated from the underwater sound field of the mesoscale vortex, and when the obtained underwater sound field of the mesoscale vortex and the sound velocity background field are subtracted, they should be in the same time period (month or season). In this embodiment, March to May are selected as spring, June to August as summer, September to November as autumn, and December to February of the following year as winter, and the sound velocity anomaly feature extraction of the four seasons is performed. The results are as follows: Figure 5 The disturbance of the vortex mainly exists in the sea water above 1000m deep, but the data in the dataset is too shallow (<5m) and the data is missing, and forced interpolation is likely to lead to false data. At the same time, in order to ensure the precision of the data, we choose to interpolate to a depth of 5:1:1000m.

[0072] In step S6 of this embodiment, when the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex is input into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, and the sound velocity background field is input into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field, it includes judging whether there is a surface sound channel in the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex. If the mesoscale vortex underwater sound velocity characteristics constitute the surface sound channel, the surface sound channel is calculated. If there is a surface sound channel in the mesoscale eddy disturbance sound velocity field formed by the underwater sound velocity characteristics of the mesoscale vortex, the sound source is set at the sea surface. Otherwise, if there is a subsurface sound channel in the mesoscale eddy disturbance sound velocity field formed by the underwater sound velocity characteristics of the mesoscale vortex, the sound source is set at the depth of the subsurface sound channel axis, and the sound ray simulation is performed based on the operation type R and the sound propagation loss is calculated based on the operation type C in the ray model of the Bellhop acoustic toolbox, so as to obtain the mesoscale vortex underwater sound channel sound ray model and the mesoscale vortex underwater sound pressure field or the climatic underwater sound pressure field. In step S6, the mesoscale vortex disturbance sound velocity field and sound velocity background field data are input into the Bellhop acoustic toolbox for sound line simulation and sound propagation loss calculation, including: constructing corresponding .env and .ssp files according to the above-mentioned mesoscale vortex disturbance sound velocity field and sound velocity background field, setting the sound source frequency to 3000Hz, and the top option to "QVW"; according to the obtained mesoscale vortex underwater sound velocity field characteristics, setting the sound source at the sea surface (when there is a surface sound channel) or the depth of the subsurface sound channel axis (when there is a subsurface sound channel), and performing sound line simulation (operation type R) and sound propagation loss calculation (operation type C). Subtracting the obtained mesoscale vortex disturbance sound pressure field from the background sound pressure field to obtain the sound propagation loss anomaly, thereby extracting the characteristics of sound waves propagating in the ocean surface and subsurface sound channels affected by the mesoscale vortex. In this embodiment, when studying the sound channel generated by the mesoscale vortex with an amplitude of 0.4 m, it is found that the anticyclonic vortex will produce a deeper surface sound channel in spring and winter, a subsurface sound channel in summer, and a double sound channel in autumn; the cyclonic vortex will lift the deep sea sound channel and transform it into a subsurface sound channel, and at the same time disintegrate the original surface sound channel of the background field in winter, such as Figure 6 shown.

[0073] In this embodiment, step S7 is used to synthesize the results obtained in step S5 and step S6, analyze the influence of mesoscale vortices on the surface and subsurface sound channels in the Kuroshio extension area, and find that the cyclonic vortex has a lifting effect on the maximum sound speed, and at the same time, by compressing the thickness of the mixed layer, the surface sound channel becomes shallower and thinner. When the intensity of the mesoscale vortex is large enough, the surface sound channel can disappear. That is, the cyclonic vortex has a reverse effect on the formation of the surface sound channel. In particular, when the cyclonic vortex intensity is large enough, the minimum sound speed can be produced at the depth of the center of the mesoscale vortex, providing conditions for the formation of the subsurface sound channel. The anticyclonic vortex has a downward pressure on the maximum sound speed, and it increases the width of the surface sound channel by deepening the thickness of the mixed layer, which has a positive effect on the formation of the surface sound channel. In particular, when the sound speed of the sea surface is very large and the sound speed of the upper layer is distributed in a strong negative gradient, such as in summer and autumn, under certain conditions, the anticyclonic vortex can produce a minimum sound speed in the subsurface layer, thereby forming a subsurface sound channel.

[0074] In summary, the method of this embodiment includes loading the Argo buoy data set, selecting hydrological parameters including temperature, salinity, and depth, and calculating the Argo sound velocity profile through the Mackenzie sound velocity empirical formula; eliminating the profiles outside the three times standard deviation by calculating the sound velocity standard deviation to achieve data quality control; loading the SODA3.4.2 climate state data set, selecting hydrological parameters including temperature, salinity, and depth, and calculating the sound velocity background field through the Mackenzie sound velocity empirical formula; loading the (META2.0) mesoscale eddy trajectory atlas to extract the surface features of the mesoscale eddy; matching and synthesizing the mesoscale eddy with the sound velocity disturbance to construct oceanographic maps with different amplitudes in different seasons. Underwater sound velocity field of mesoscale eddies in the ocean; subtract the background sound velocity from the obtained sound velocity field of mesoscale eddy disturbance to obtain the abnormal variation characteristics of sound velocity; load the mesoscale eddy data set to extract the occurrence time, central longitude and latitude coordinates, polarity, radius, amplitude and other parameters of the mesoscale eddy; use the Bellhop acoustic toolbox to simulate the sound line and calculate the sound propagation loss to obtain the sound propagation loss of sound waves under the sound velocity background field and the sound velocity field of mesoscale eddy disturbance; subtract the sound propagation loss under the interference of mesoscale eddies with different amplitudes in each season from the sound propagation loss in the absence of interference (background field) to obtain the sound propagation loss anomaly, and analyze the propagation characteristics of sound waves in the surface and subsurface sound channels and the influence of mesoscale eddies on sound wave propagation. The present invention can construct universal characteristics of underwater surface and subsurface acoustic channels of mesoscale vortices, provide theoretical support for the variability of surface and subsurface acoustic channels under the occurrence of mesoscale vortices with different characteristics, provide a reference for the study of sound propagation effects in a mesoscale vortex environment, and provide theoretical support for the study of the variability characteristics of underwater surface and subsurface acoustic channels of different mesoscale vortices.

[0075] In addition, this embodiment also provides a system for constructing the sound channels of the surface and subsurface layers affected by mesoscale eddies in the ocean, comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to execute the method of constructing the sound channels of the surface and subsurface layers affected by mesoscale eddies in the ocean.

[0076] In addition, this embodiment also provides a computer-readable storage medium, which stores a computer program / instructions, and the computer program or instructions are programmed or configured to execute the method of constructing the ocean mesoscale eddy-affected surface and subsurface sound channels through a processor.

[0077] In addition, this embodiment also provides a computer program product, including a computer program / instruction, which is programmed or configured to execute the method of constructing the ocean mesoscale eddy-affected surface and subsurface sound channels through a processor.

[0078] Those skilled in the art should understand that the technical solutions provided by the embodiments of the present application may be in the form of methods, systems, or computer program products. Therefore, the present application may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may take the form of a computer program product implemented on one or more computer-readable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the process Figure 1 A process or multiple processes and / or boxes Figure 1 These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing device to work in a specific way, so that the instructions stored in the computer-readable memory produce a product including an instruction device, which implements the functions specified in the process. Figure 1 A process or multiple processes and / or boxes Figure 1These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide for implementing the process in the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0079] The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should also be regarded as the protection scope of the present invention.

Claims

1. A method for constructing an ocean mesoscale eddy impact surface and subsurface sound channel, characterized in that: The steps include: S1, obtains the hydrological parameters of the specified area within the specified time according to the input Argo float dataset and SODA3.4.2 climate state dataset; S2, the hydrological parameters of the Argo float data set are used to calculate the Argo sound velocity profile of the specified area using the Mackenzie sound velocity empirical formula, and the hydrological parameters of the SODA3.4.2 climate state data set are used to calculate the sound velocity background field of the specified area using the Mackenzie sound velocity empirical formula. The Argo sound velocity profile calculated by the Argo float data set is quality controlled to remove bad data to obtain the Argo sound velocity profile data set; S3, extracting the characteristics of the mesoscale eddy according to the input mesoscale eddy trajectory dataset; S4, based on the characteristics of the mesoscale vortex and the position and time of the Argo float, the temporal and spatial matching and merging of the Argo sound velocity profile data set and the mesoscale vortex trajectory data set are completed, and vertical and horizontal smoothing is performed to obtain the underwater sound velocity characteristics of the mesoscale vortex and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex; S5, interpolating the underwater sound field and the sound velocity background field of the mesoscale vortex formed by the characteristics of the mesoscale vortex to the specified depth respectively, subtracting the sound velocity background field from the underwater sound field of the mesoscale vortex to separate the sound velocity background field, and then extracting the sound velocity anomaly characteristics in four seasons of the year, thereby obtaining the underwater sound velocity disturbance structure of the mesoscale vortex; S6, input the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, input the sound velocity background field into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field; subtract the climatic underwater sound pressure field from the mesoscale vortex underwater sound pressure field to obtain the abnormal data of sound propagation loss under the influence of the mesoscale vortex; S7, combining the abnormal data of sound propagation loss under the influence of mesoscale eddies and the underwater sound velocity disturbance structure of mesoscale eddies, extracts the characteristics of the surface and subsurface sound channels affected by mesoscale eddies in the ocean.

2. The method for constructing the ocean mesoscale eddy-induced surface and subsurface acoustic channels according to claim 1, characterized in that: The hydrological parameters in step S1 include the temperature, salinity and depth of the seawater.

3. The method for constructing the ocean mesoscale eddy-induced surface and subsurface sound channel according to claim 2, characterized in that: The functional expression of the Mackenzie sound velocity empirical formula in step S2 is: ， In the above formula, represents the speed of sound, The temperature of the sea water, The salinity of seawater. Indicates the depth of sea water.

4. The method for constructing the ocean mesoscale eddy-induced surface and subsurface acoustic channels according to claim 1, characterized in that: In step S2, quality control is performed on the Argo sound velocity profile calculated by the Argo float data set to remove bad data to obtain the Argo sound velocity profile data set, including: calculating the sound velocity standard deviation of the Argo sound velocity profile calculated by the Argo float data set, and deleting the Argo sound velocity profile outside three times the sound velocity standard deviation, thereby achieving quality control to remove bad data to obtain the Argo sound velocity profile data set.

5. The method for constructing the ocean mesoscale eddy-induced surface and subsurface acoustic channels according to claim 1, characterized in that: Step S4 includes: S4.1, based on the characteristics of the mesoscale vortex including the central coordinates, appearance time, radius, polarity and amplitude of the mesoscale vortex, the influence range and time of each mesoscale vortex are calculated using the central coordinates, appearance time and radius of the mesoscale vortex, and then the Argo sound velocity profile dataset is spatially and temporally matched with the mesoscale vortex dataset based on the position and time recorded by the Argo float; S4.2, based on the center coordinates of the mesoscale vortex and the coordinates of the Argo float within twice its radius, obtain the radial normalized coordinates of the underwater sound speed profile of each mesoscale vortex relative to the center of the mesoscale vortex, and then integrate the underwater sound speed profiles of mesoscale vortices with the same characteristics into the same mesoscale vortex based on the normalized coordinates; S4.3, set the sliding window size multiple times, find the optimal sliding window size according to the smoothing results of each vertical and horizontal smoothing process, and use the optimal sliding window sizes in the vertical and horizontal directions to perform vertical and horizontal smoothing processes respectively, so as to obtain the underwater sound velocity characteristics of mesoscale vortices with different amplitudes under different sound velocity background fields, and the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of mesoscale vortices.

6. The method for constructing the ocean mesoscale eddy-induced surface and subsurface acoustic channels according to claim 5, characterized in that: Step S4.2 includes: S4.2.1, obtain the Argo sound velocity profile for each mesoscale eddy pair based on the coordinates of the center coordinates of the mesoscale eddy and the coordinates of the Argo float within the range of twice its radius on the same day of observation; S4.2.2, project all paired Argo sound velocity profiles into a coordinate system with the vortex center of the corresponding mesoscale vortex as the origin Below, among them To fall in The first Argo sound velocity profiles and vortex centers The distance is the z-axis coordinate; S4.2.3, calculate the The radial normalized distance of each Argo sound speed profile within a mesoscale vortex: ， In the above formula, is the radial normalized distance of the kth Argo sound velocity profile, , For the The radius of the mesoscale vortex is , To fall to The number of mesoscale eddies; S4.2.4, the radial normalized distance of each Argo sound velocity profile Arrange them in the vortex center coordinate system according to their size, and then according to the radius of the synthetic vortex integrated into the same mesoscale eddy.

7. The method for constructing the ocean mesoscale eddy-induced surface and subsurface sound channel according to claim 1, characterized in that: In step S6, the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex is input into the ray model of the Bellhop acoustic toolbox to calculate the mesoscale vortex underwater sound channel sound line model and the mesoscale vortex underwater sound pressure field, and the sound velocity background field is input into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field, including judging whether there is a surface sound channel in the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex. If the mesoscale vortex underwater sound velocity characteristics constitute the surface sound channel, the sound velocity background field is input into the ray model of the Bellhop acoustic toolbox to calculate the climatic underwater sound pressure field. If there is a surface sound channel in the sound velocity field of the mesoscale vortex disturbance, the sound source is set on the sea surface. Otherwise, if there is a subsurface sound channel in the mesoscale vortex disturbance sound velocity field composed of the underwater sound velocity characteristics of the mesoscale vortex, the sound source is set at the depth of the subsurface sound channel axis, and the sound ray simulation is performed based on the operation type R and the sound propagation loss is calculated based on the operation type C in the ray model of the Bellhop acoustic toolbox, so as to obtain the mesoscale vortex underwater sound channel sound ray model and the mesoscale vortex underwater sound pressure field or the climatic underwater sound pressure field.

8. A system for constructing an ocean mesoscale eddy impact surface and subsurface sound channel, comprising a microprocessor and a memory connected to each other, characterized in that: The microprocessor is programmed or configured to execute the method for constructing the ocean mesoscale eddy-affected surface and subsurface sound channels as described in any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program or instruction stored therein, characterized in that: The computer program or instruction is programmed or configured to execute the method for constructing the ocean mesoscale eddy-affected surface and subsurface sound channels as described in any one of claims 1 to 7 through a processor.

10. A computer program product comprising a computer program or instructions, characterized in that The computer program or instruction is programmed or configured to execute the method for constructing the ocean mesoscale eddy-affected surface and subsurface sound channels as described in any one of claims 1 to 7 through a processor.

Images