Density layering flow velocity correction method and system based on linear multilayer non-static pressure theory

Abstract

The invention discloses a density stratification flow velocity correction method and system based on a linear multilayer non-static pressure theory, and the method comprises the steps: obtaining seawater density distribution and water level change distribution of a target sea area, and calculating a seawater dimensionless density profile based on the seawater density distribution; layering the target sea area along the longitudinal direction to obtain a plurality of seawater layers, and calculating the dimensionless density and the dimensionless thickness corresponding to each seawater layer based on the seawater dimensionless density profile; calculating a flow velocity attenuation coefficient based on the total number of seawater layers, the dimensionless density and the dimensionless thickness; performing spatial Fourier transform on the water level change distribution to obtain wavenumber domain distribution of the water level change; on the basis of a linear wave theory, calculating wave flow velocities under different wave numbers according to wave number domain distribution; and correcting the wave flow velocity based on the flow velocity attenuation coefficient, and outputting the corrected wave flow velocity.

Description

Technical Field

[0001] The invention belongs to the field of ocean wave velocity assessment and correction, and in particular relates to a density stratified velocity correction method and system based on linear multi-layer non-hydrostatic pressure theory. Background Art

[0002] Ocean fluctuations are of great significance to marine engineering, weather forecasting, marine ecological protection, etc. In marine engineering, such as offshore wind power, port construction and marine resource development, accurate prediction of ocean fluctuations is the prerequisite for designing safe and reliable engineering structures; when conducting marine meteorological forecasts, ocean fluctuations are often relied on to assess the dynamic changes of the ocean, enhance the timeliness and accuracy of forecasts, and provide more reliable meteorological services to society; ocean fluctuations not only affect the transfer of heat and momentum, but also have a direct impact on the habitat of marine organisms. Ocean fluctuations are manifested as wave velocity, which is a complex phenomenon with multiple levels and multiple factors interacting, involving multiple fields such as fluid dynamics, thermodynamics and environmental science.

[0003] As climate change, marine pollution and human activities intensify the impact on marine ecosystems, it becomes increasingly important to accurately predict and understand the behavior of ocean waves. In this context, researchers continue to explore more accurate wave models to improve the ability to predict ocean waves. However, traditional wave models are often based on the assumption of static pressure or uniform media, which makes it difficult to reflect the complex conditions in the ocean, such as density stratification caused by changes in temperature and salinity, and flow inhomogeneity, resulting in poor prediction accuracy of wave velocity.

[0004] Therefore, how to provide a density-stratified velocity correction method based on linear multi-layer non-static pressure to improve the prediction accuracy of wave velocity has become a technical problem that needs to be solved urgently. Summary of the invention

[0005] The purpose of the present invention is to overcome the deficiencies in the prior art and to provide a density stratified flow velocity correction method and system based on the linear multilayer non-static pressure theory.

[0006] In a first aspect, an embodiment of the present invention provides a density stratified flow velocity correction method based on linear multilayer non-hydrostatic pressure theory, the method comprising the following steps:

[0007] Obtain the seawater density distribution and water level change distribution in the target sea area, and calculate the dimensionless density profile of seawater based on the seawater density distribution;

[0008] The target sea area is layered longitudinally to obtain several seawater layers, and the dimensionless density and dimensionless thickness corresponding to each seawater layer are calculated based on the dimensionless density profile of seawater;

[0009] The velocity attenuation coefficient is calculated based on the total number of seawater layers, dimensionless density, and dimensionless thickness;

[0010] Perform spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, calculate the wave velocity under different wave numbers according to the wave number domain distribution;

[0011] The wave velocity is corrected based on the velocity attenuation coefficient and the corrected wave velocity is output.

[0012] In a second aspect, an embodiment of the present invention provides a density stratified flow rate correction system based on the linear multilayer non-static pressure theory to implement the above-mentioned density stratified flow rate correction method based on the linear multilayer non-static pressure theory, and the system includes:

[0013] The seawater data acquisition module is used to obtain the seawater density distribution and water level change distribution in the target sea area, and calculate the dimensionless density profile of seawater based on the seawater density distribution;

[0014] The seawater stratification module is used to stratify the target sea area longitudinally to obtain several seawater layers, and calculate the dimensionless density and dimensionless thickness corresponding to each seawater layer based on the dimensionless density profile of seawater;

[0015] A velocity attenuation coefficient calculation module, used to calculate the velocity attenuation coefficient based on the total number of seawater layers, dimensionless density and dimensionless thickness;

[0016] The wave velocity calculation module is used to perform spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, the wave velocity under different wave numbers is calculated according to the wave number domain distribution;

[0017] The wave velocity correction module is used to correct the wave velocity based on the velocity attenuation coefficient and output the corrected wave velocity.

[0018] In a third aspect, an embodiment of the present invention provides an electronic device comprising a memory and a processor, wherein the memory is coupled to the processor; wherein the memory is used to store program data, and the processor is used to execute the program data to implement the above-mentioned density stratified flow velocity correction method based on linear multilayer non-static pressure theory.

[0019] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned density stratified flow velocity correction method based on the linear multilayer non-hydrostatic pressure theory.

[0020] In a fifth aspect, an embodiment of the present invention provides a computer program product, including a computer program / instruction, which, when executed by a processor, implements the above-mentioned density stratified flow velocity correction method based on linear multilayer non-hydrostatic pressure theory.

[0021] Compared with the prior art, the present invention has the following beneficial effects:

[0022] The present invention provides a density stratified velocity correction method based on linear multi-layer non-hydrostatic theory, comprising: obtaining seawater density distribution and water level change distribution in a target sea area, and calculating a dimensionless density profile of seawater based on the seawater density distribution; stratifying the target sea area longitudinally to obtain a plurality of seawater layers, and calculating a dimensionless density and a dimensionless thickness corresponding to each seawater layer based on the dimensionless density profile of seawater; calculating a velocity attenuation coefficient based on the total number of layers, dimensionless density and dimensionless thickness of the seawater layer; performing spatial Fourier transformation on the water level change distribution to obtain a wave number domain distribution of the water level change; calculating wave velocities at different wave numbers based on the linear wave theory and the wave number domain distribution; and correcting the wave velocity based on the velocity attenuation coefficient, and outputting the corrected wave velocity. The present invention considers the motion characteristics of each seawater layer and their interactions respectively to comprehensively capture the propagation characteristics and energy distribution of ocean waves, reveal the dynamic relationship between different seawater layers, correct the shortcomings of traditional wave models, calculates wave velocities under different wave numbers based on the water level change distribution, and corrects the wave velocity by the velocity attenuation coefficient. Arbitrary accuracy can be achieved by continuously superimposing the number of seawater layers, thereby improving the prediction accuracy of wave velocity. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative labor.

[0024] Figure 1 It is a flow chart of a density stratified flow velocity correction method based on linear multi-layer non-static pressure provided by an embodiment of the present invention;

[0025] Figure 2 is a schematic diagram of a dimensionless density profile of seawater provided by an embodiment of the present invention;

[0026] Figure 3 is a parameter schematic diagram of the dimensionless density profile of eight layers of seawater provided by an embodiment of the present invention;

[0027] Figure 4 is a schematic diagram of wave velocity correction provided by an embodiment of the present invention;

[0028] Figure 5 It is a structural schematic diagram of a density stratified flow velocity correction system based on linear multi-layer non-static pressure provided by an embodiment of the present invention;

[0029] Figure 6 It is a schematic diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0030] The present invention is further described below in conjunction with embodiments. The description of the following embodiments is only used to help understand the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principles of the present invention, several improvements and modifications may be made to the present invention, and these improvements and modifications also fall within the scope of protection of the claims of the present invention.

[0031] In addition, the technical solutions between the various embodiments can be combined with each other, but they must be based on the fact that they can be implemented by ordinary technicians in this field. When the combination of technical solutions is contradictory or cannot be implemented, it should be deemed that such combination of technical solutions does not exist and is not within the scope of protection required by this application.

[0032] like Figure 1 As shown, an embodiment of the present invention provides a density stratified flow velocity correction method based on linear multilayer non-static pressure theory, and the method includes the following steps:

[0033] Step S1, obtaining the seawater density distribution and water level change distribution of the target sea area, and calculating the dimensionless density profile of the seawater based on the seawater density distribution.

[0034] The process of calculating the dimensionless density profile of seawater based on the seawater density distribution includes:

[0035] The dimensionless density profile of seawater is obtained by dividing the seawater density distribution by a characteristic density value; wherein the characteristic density value is the maximum value in the seawater density distribution or the seabed density value.

[0036] Step S2, stratifying the target sea area longitudinally to obtain a plurality of seawater layers, and calculating the dimensionless density and dimensionless thickness corresponding to each seawater layer based on the dimensionless density profile of the seawater.

[0037] Furthermore, the process of stratifying the target sea area longitudinally to obtain a plurality of seawater layers includes:

[0038] Set the number of sampling points, and stratify the target sea area longitudinally according to the number of sampling points, so that the total number of seawater layers is greater than twice the number of sampling points;

[0039] or,

[0040] The total number of seawater layers is set according to the empirical formula, and the expression is as follows:

[0041]

[0042] Where N represents the total number of seawater layers. Indicates the size of the error, represents the dimensionless density of the seafloor, Represents the dimensionless density value of the sea surface.

[0043] Step S3, calculating the velocity attenuation coefficient based on the total number of seawater layers, dimensionless density and dimensionless thickness.

[0044] Furthermore, the expression of the flow velocity attenuation coefficient is as follows:

[0045]

[0046] In the formula, represents the velocity attenuation coefficient, is the acceleration due to gravity, For water depth, is the wave number, Indicates The dimensionless density value corresponding to the seawater layer, Indicates The dimensionless thickness value corresponding to the seawater layer, It represents the pressure change of water under the action of sea surface waves. Represents the volume transport caused by gravity during wave propagation.

[0047] Step S4, performing spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, calculating the wave flow velocity at different wave numbers according to the wave number domain distribution.

[0048] Specifically, assuming that at the initial moment When , perform Fourier transform on the distribution in the spatial domain to obtain the wavenumber domain distribution of water level change, which is expressed as follows:

[0049]

[0050] In the formula, express The Fourier transform in the wavenumber domain is is an imaginary unit, is the wave number;

[0051] Based on the linear wave theory, The fluctuations of the sea surface are as follows:

[0052]

[0053] In the formula, is the wave speed. In linear wave theory, the expression of wave speed is as follows:

[0054]

[0055] In the formula, is the acceleration due to gravity, For water depth.

[0056] Step S5, correcting the wave velocity based on the velocity attenuation coefficient, and outputting the corrected wave velocity.

[0057] Furthermore, the linear wave theory does not consider the influence of density stratification, so the wave speed it calculates is too fast. This patent corrects the wave speed c by considering the stratification of seawater density, so that the model can obtain more accurate results. The expression of the corrected wave velocity is as follows:

[0058]

[0059] In the formula, is the wave speed, represents the velocity attenuation coefficient, is the acceleration due to gravity, For water depth, is the wave number.

[0060] Next, this example elaborates on a density stratified velocity correction algorithm based on linear multilayer non-hydrostatic pressure theory - the process of theoretical derivation.

[0061] In order to facilitate the demonstration of the detailed calculation process, this example assumes that the model is in two-dimensional space. Calculated in the coordinate system. Figure 2 As shown, is the sea surface height, is the seafloor depth, defined is the total water depth. In this example, the water flow is divided into The thickness of each water layer is defined as:

[0062]

[0063] Among them, the coefficient is the water layer thickness coefficient, and its value range is 0-1. Indicates The water layer and The interface layer between the water layers, where the sea surface With the seabed Can be regarded as a special interface layer, namely:

[0064]

[0065] Therefore, each interface layer The position can be calculated:

[0066]

[0067] It should be noted that when the lower bound of the summation symbol in formula (4) is greater than the upper bound, the summation term is considered to be 0. Indicates The density of the water layer is assumed to be constant, Indicated at the interface layer The density of the seafloor and the sea surface is approximately , ; The density on the remaining interface layer is approximately:

[0068]

[0069] Based on the non-hydrostatic free surface flow model, the continuity equation and momentum equation are:

[0070]

[0071]

[0072]

[0073] The boundary conditions are:

[0074]

[0075]

[0076] in, represents density, Respectively represent the horizontal direction and vertical direction speed, Indicates time, For sea surface changes, is the acceleration due to gravity, and are the vertical velocities of the sea surface and seabed, respectively, and are the horizontal velocities of the sea surface and seabed, respectively, is the pressure term, and its expression is:

[0077]

[0078] in, is the characteristic density, which is a constant; is a non-static pressure term. The pressure term is on the sea surface The value is 0, that is .

[0079] Equations (6)-(8) along each layer arrive Integrate, substitute the density approximation formula (5), boundary conditions (9)-(10), and pressure decomposition formula (11), perform linear approximation and remove the nonlinear terms to obtain the following equation:

[0080]

[0081]

[0082]

[0083] in, The range is 1 to , Indicates The horizontal and vertical speeds of the layer. Indicates The vertical velocity of the boundary layer is Indicates The horizontal and vertical velocities of the boundary layer. Indicates The boundary layer and The non-static pressure of the boundary layer.

[0084] and Indicates and The dimensionless density value of the water layer is expressed as:

[0085]

[0086] Indicates The effect coefficient of the static pressure term of the water layer is expressed as:

[0087]

[0088] At this time, the linearized boundary conditions are written as:

[0089]

[0090]

[0091] Equations (12)-(14) represent the Movement in the water layer, so the whole system has a total of equation, the following example changes the variables Make the following transformations:

[0092]

[0093] in represents the new variable after transformation. At this time, equation (13) can be rewritten as:

[0094]

[0095] Transform equation (12) from the seafloor to the The vertical velocity of the interface layer can be obtained by accumulating the water layer and substituting it into the new variable (15). The expression is:

[0096]

[0097] in and Respectively about and The expression is:

[0098]

[0099]

[0100] In addition, When , according to (21) and boundary conditions (17), we can obtain:

[0101]

[0102] According to the gradient formula, the vertical velocity in the water layer It is approximately estimated as the average vertical velocity of the adjacent interface layer, that is:

[0103]

[0104] Transform equation (14) from The non-static pressure term of the interface layer can be obtained by accumulating the pressure from the water layer to the sea surface and substituting it into formula (25): The expression is:

[0105]

[0106] According to formula (26), the following expression can be calculated:

[0107]

[0108] in is the coefficient, and its expression is:

[0109]

[0110] Substituting (27) into (20), we can obtain the following expression:

[0111]

[0112] The coefficient The expression is:

[0113]

[0114] coefficient Layer thickness coefficient and layer density coefficient Once these two sets of parameters are determined, the coefficients can be directly given. unknown quantity Next, we calculate the wave velocity based on a small amplitude system, and let:

[0115]

[0116]

[0117] in and Respectively and The amplitude of represents the wave number, is the angular frequency, is the imaginary unit. Substituting (31) and (32) into (24) and (29), we can obtain and The linear equations are:

[0118]

[0119]

[0120] The coefficient matrix form of the linear equations formed by (33) and (34) is:

[0121]

[0122] If the system of equations is to have a non-zero solution, it is necessary to ensure that the determinant of the matrix (35) is 0. Calculate the determinant of (35) and put Substituting in, first the first column is proposed , the remaining columns are presented ; then multiply the first row by , the first column is divided by ,get:

[0123]

[0124] If the value of formula (36) is 0, it can be deduced that:

[0125]

[0126] in, Representation Matrix The algebraic co-factor of the element in the first row and first column of is:

[0127]

[0128] All horizontal velocity terms in the momentum equation (29) are included. The coefficient of pressure, which represents the change in pressure on the water body due to sea surface waves. Indicates row 1, The algebraic co-factor of the column elements can be written as:

[0129]

[0130] In the formula, Characterizing sea surface spatial gradients and horizontal velocities interaction, which describes the volume transport caused by gravity during wave propagation.

[0131] in, Indicates Row, No. The algebraic cofactor of the column elements. So the wave speed You can solve:

[0132]

[0133] The attenuation coefficient It can be expressed as:

[0134]

[0135] Formula (41) gives the velocity attenuation coefficient when the number of layers N is determined, which can be applied to the correction calculation of the velocity. However, how to determine the value of N still needs to be discussed. An approximate estimation method will be given below.

[0136] when , the wave speed (40) can be written as:

[0137]

[0138] If we assume that each water layer is equally divided,

[0139]

[0140] Calculate now The situation tends to infinity:

[0141]

[0142] If the water depth is mapped to [0,1], (44) can be rewritten as:

[0143]

[0144] Formula (45) can be used to calculate any continuous density change. Assuming that the density changes linearly from the seabed to the sea surface, the flow velocity is:

[0145]

[0146] The flow rate of discrete layers is:

[0147]

[0148] If we want to ensure that the results of (46) and (47) are close, we can calculate them by relative error, which can be concluded:

[0149]

[0150] in, Indicates the error size, usually 10 -3 Formula (48) provides a simple estimation method, which can make a quick estimate of the number of layers for different density profiles.

[0151] In order to demonstrate the calculation process, the above conclusions are applied to the generation and propagation of tsunami waves in two-dimensional space. In the above figure, it is assumed that the initial seafloor topography is flat and the seafloor deformation for:

[0152]

[0153] in, is the Heaviside function, is the rupture velocity, the calculation formula for the sea surface is:

[0154]

[0155] in, This is the attenuation coefficient (i.e., correction coefficient) calculated by formula (41). is the Laplace transform operator, is the angular frequency.

[0156] To demonstrate the calculation process, the above conclusions can be applied to simulate the generation and propagation of tsunami waves. In the above example, it is assumed that the seafloor deformation for:

[0157]

[0158] in represents the horizontal width of the seafloor deformation, represents the Heavside function, and its formula is as follows:

[0159]

[0160] The corresponding initial moment of generation The sea surface at that time is:

[0161]

[0162] Therefore When , the sea surface fluctuation is

[0163]

[0164] At this time, for formula (47) Corrected, setting the sea surface density to 1020 kg / m 3 , the seafloor density is 1050 kg / m 3 The density changes linearly from the bottom to the bottom. =10 -3 According to formula (48), when the seawater layer is greater than or equal to eight layers, it is sufficient to ensure that the error is less than the given requirement. Therefore, this example takes the eight-layer seawater layer model as an example. Figure 3 Set the horizontal width of the seabed deformation kilometers, the calculation result is Figure 4 As shown, the dotted line is the result without density, and the solid line is the result with density. It can be found that the water level of the result with density is slightly higher and the propagation speed is slightly slower. Therefore, this method provides a fast flow velocity correction algorithm.

[0165] On the other hand, Figure 5 As shown, the present invention provides a density stratified flow rate correction system based on linear multilayer non-static pressure theory to implement the above-mentioned density stratified flow rate correction method based on linear multilayer non-static pressure theory, and the system includes:

[0166] The seawater data acquisition module is used to obtain the seawater density distribution and water level change distribution in the target sea area, and calculate the dimensionless density profile of seawater based on the seawater density distribution;

[0167] The seawater stratification module is used to stratify the target sea area longitudinally to obtain several seawater layers, and calculate the dimensionless density and dimensionless thickness corresponding to each seawater layer based on the dimensionless density profile of seawater;

[0168] A velocity attenuation coefficient calculation module, used to calculate the velocity attenuation coefficient based on the total number of seawater layers, dimensionless density and dimensionless thickness;

[0169] The wave velocity calculation module is used to perform spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, the wave velocity under different wave numbers is calculated according to the wave number domain distribution;

[0170] The wave velocity correction module is used to correct the wave velocity based on the velocity attenuation coefficient and output the corrected wave velocity.

[0171] Regarding the system in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0172] For the system embodiment, since it basically corresponds to the method embodiment, the relevant parts can refer to the partial description of the method embodiment. The system embodiment described above is only schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the present application. A person of ordinary skill in the art can understand and implement it without creative work.

[0173] Accordingly, the present application also provides an electronic device, comprising: one or more processors; a memory for storing one or more programs; when the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned density stratified flow velocity correction method based on the linear multilayer non-static pressure theory. Figure 6 As shown in FIG. 1 , a hardware structure diagram of a device having data processing capability is provided for the density stratified flow velocity correction method based on the linear multilayer non-static pressure theory provided by an embodiment of the present invention, except Figure 6 In addition to the processor, memory and network interface shown, any device with data processing capability in which the apparatus in the embodiment is located may also include other hardware, generally based on the actual functions of the device with data processing capability, which will not be described in detail.

[0174] Accordingly, the present application also provides a computer-readable storage medium on which computer instructions are stored, and when the instructions are executed by the processor, the density stratified flow rate correction method based on the linear multilayer non-static pressure theory as described above is implemented. The computer-readable storage medium can be an internal storage unit of any device with data processing capabilities described in any of the aforementioned embodiments, such as a hard disk or a memory. The computer-readable storage medium can also be an external storage device, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), an SD card, a flash card (Flash Card), etc. equipped on the device. Furthermore, the computer-readable storage medium can also include both an internal storage unit and an external storage device of any device with data processing capabilities. The computer-readable storage medium is used to store the computer program and other programs and data required by any device with data processing capabilities, and can also be used to temporarily store data that has been output or is to be output.

[0175] Those skilled in the art will readily appreciate other embodiments of the present application after considering the description and practicing the contents disclosed herein. The present application is intended to cover any variations, uses or adaptations of the present application, which follow the general principles of the present application and include common knowledge or customary technical means in the art that are not disclosed in the present application. The description and examples are intended to be exemplary only.

[0176] It should be understood that the present application is not limited to the exact construction that has been described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof.

Claims

1. A density stratified flow velocity correction method based on linear multilayer non-static pressure theory, characterized in that: The method comprises the following steps: Obtain the seawater density distribution and water level change distribution in the target sea area, and calculate the dimensionless density profile of seawater based on the seawater density distribution; The target sea area is layered longitudinally to obtain several seawater layers, and the dimensionless density and dimensionless thickness corresponding to each seawater layer are calculated based on the dimensionless density profile of seawater; The velocity attenuation coefficient is calculated based on the total number of seawater layers, dimensionless density, and dimensionless thickness; Perform spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, calculate the wave velocity under different wave numbers according to the wave number domain distribution; The wave velocity is corrected based on the velocity attenuation coefficient and the corrected wave velocity is output.

2. The density stratified flow velocity correction method based on linear multilayer non-static pressure theory according to claim 1 is characterized in that: The process of calculating the dimensionless density profile of seawater based on the seawater density distribution includes: The dimensionless density profile of seawater is obtained by dividing the seawater density distribution by a characteristic density value; wherein the characteristic density value is the maximum value in the seawater density distribution or the seabed density value.

3. The density stratified flow velocity correction method based on linear multilayer non-static pressure theory according to claim 1 is characterized in that: The process of stratifying the target sea area longitudinally to obtain several seawater layers includes: Set the number of sampling points, and stratify the target sea area longitudinally according to the number of sampling points, so that the total number of seawater layers is greater than twice the number of sampling points; or, The total number of seawater layers is set according to the empirical formula, and the expression is as follows: ； Where N represents the total number of seawater layers. Indicates the size of the error, represents the dimensionless density of the seafloor, Represents the dimensionless density value of the sea surface.

4. The density stratified flow velocity correction method based on linear multilayer non-static pressure theory according to claim 1 is characterized in that: The process of calculating the velocity attenuation coefficient based on the total number of seawater layers, dimensionless density, and dimensionless thickness includes: ； In the formula, represents the velocity attenuation coefficient, is the acceleration due to gravity, For water depth, is the wave number, Indicates The dimensionless density value corresponding to the seawater layer, Indicates The dimensionless thickness value corresponding to the seawater layer, It represents the pressure change of water under the action of sea surface waves. It represents the volume transport caused by gravity during the propagation of waves.

5. The density stratified flow velocity correction method based on linear multilayer non-static pressure theory according to claim 1 is characterized in that: The water level change distribution is spatially Fourier transformed to obtain the wave number domain distribution of the water level change; based on the linear wave theory, the process of calculating the wave velocity under different wave numbers according to the wave number domain distribution includes: Assume that at the initial moment When , perform Fourier transform on the distribution in the spatial domain to obtain the wavenumber domain distribution of water level change, which is expressed as follows: ； In the formula, express The Fourier transform in the wavenumber domain is is an imaginary unit, is the wave number; Based on the linear wave theory, The sea surface fluctuations are as follows: ； In the formula, is the wave speed. In linear wave theory, the expression of wave speed is as follows: ； In the formula, is the acceleration due to gravity, For water depth.

6. A density stratified flow velocity correction method based on linear multilayer non-static pressure theory according to claim 1 or 4, characterized in that: The wave velocity is corrected based on the velocity attenuation coefficient to obtain the corrected wave velocity, which is expressed as follows: ； In the formula, is the wave speed, represents the velocity attenuation coefficient, is the acceleration due to gravity, For water depth, is the wave number.

7. A density stratified flow rate correction system based on linear multilayer non-static pressure theory, characterized in that: To implement the density stratified flow velocity correction method based on linear multilayer non-static pressure theory as described in any one of claims 1 to 6, the system comprises: The seawater data acquisition module is used to obtain the seawater density distribution and water level change distribution in the target sea area, and calculate the dimensionless density profile of seawater based on the seawater density distribution; The seawater stratification module is used to stratify the target sea area longitudinally to obtain several seawater layers, and calculate the dimensionless density and dimensionless thickness corresponding to each seawater layer based on the dimensionless density profile of seawater; A velocity attenuation coefficient calculation module, used to calculate the velocity attenuation coefficient based on the total number of seawater layers, dimensionless density and dimensionless thickness; The wave velocity calculation module is used to perform spatial Fourier transform on the water level change distribution to obtain the wave number domain distribution of the water level change; based on the linear wave theory, the wave velocity under different wave numbers is calculated according to the wave number domain distribution; The wave velocity correction module is used to correct the wave velocity based on the velocity attenuation coefficient and output the corrected wave velocity.

8. An electronic device, comprising a memory and a processor, characterized in that: The memory is coupled to the processor; wherein the memory is used to store program data, and the processor is used to execute the program data to implement the density stratified flow velocity correction method based on linear multilayer non-static pressure theory as described in any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the density stratified flow velocity correction method based on the linear multilayer non-hydrostatic pressure theory as described in any one of claims 1 to 6 is implemented.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instruction is executed by a processor, the density stratified flow velocity correction method based on the linear multilayer non-hydrostatic pressure theory described in any one of claims 1 to 6 is implemented.

Images