CN116432562B - Analysis method for influence degree of density jump layer element change on internal tide - Google Patents

Analysis method for influence degree of density jump layer element change on internal tide Download PDF

Info

Publication number
CN116432562B
CN116432562B CN202310691642.XA CN202310691642A CN116432562B CN 116432562 B CN116432562 B CN 116432562B CN 202310691642 A CN202310691642 A CN 202310691642A CN 116432562 B CN116432562 B CN 116432562B
Authority
CN
China
Prior art keywords
density
tide
sea
internal
depth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310691642.XA
Other languages
Chinese (zh)
Other versions
CN116432562A (en
Inventor
李秉天
杜立彬
张宇
张清华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202310691642.XA priority Critical patent/CN116432562B/en
Publication of CN116432562A publication Critical patent/CN116432562A/en
Application granted granted Critical
Publication of CN116432562B publication Critical patent/CN116432562B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides an analysis method for the influence degree of density jump layer element change on internal tide, and belongs to the technical field of internal tide simulation. The method comprises the following steps: s1, selecting the type of the topography of the ocean area, setting topography parameters and density jump layer parameters, wherein the density jump layer parameters comprise strength, depth and thickness, and calculating to obtain an initial density field; s2, applying trend velocity forcing in an open boundary, performing internal tide simulation, and outputting horizontal flow velocity, ocean temperature and ocean salinity distribution; s3, calculating the sea water density by using the sea temperature and the sea salinity; s4, changing the density jump layer parameters, and repeating the steps S2-S3; s5, analyzing the influence of the density jump layer parameter change on the internal tide intensity and/or the ocean density according to the distribution difference of the horizontal flow velocity and the sea water density. According to the invention, the influence degree of different elements of the density jump layer on the internal tide can be quantitatively analyzed by parameterizing the density jump layer lamination elements and analyzing the difference of the internal tide intensity and the ocean density under different density jump layer parameters.

Description

Analysis method for influence degree of density jump layer element change on internal tide
Technical Field
The invention relates to the technical field of internal tide simulation, in particular to an analysis method for the influence degree of density jump layer element change on internal tide.
Background
The internal tide is a low-frequency internal wave with tidal frequency, is an important medium for redistributing the mechanical energy in the ocean, and the intensity of the internal tide directly determines the fluctuation degree of the equal-density surface. Internal tide and its intensity change have important effects in many fields of production and life. For example, sea water mixing caused by internal tide or lifting of internal equal-density surface can cause great fluctuation in the sea, and acoustic signals propagated underwater and an internal interface are reflected and refracted at different angles, so that the propagation speed and propagation direction of the underwater signals are changed, the acoustic channels are finally changed, the acoustic signals are damaged, the underwater communication and the detection of targets are seriously hindered, and when a submarine in sailing encounters strong internal tide, the stability and the safety of the submarine face great challenges. Therefore, research on the influence factors of the internal tide intensity has very important value in the aspects of marine fishery production, naval vessel safety, port development, energy development, environmental protection and the like.
The sea water density layer junction variation, especially at the density jump layer, has a significant effect on the degree of fluctuation of internal tide such as internal tide intensity or constant density surface. The depth, thickness and strength of the skip-and-skip can all significantly affect the internal tide, and it is of interest how much these three main elements of the skip-and-skip-each affect the internal tide strength. Most of current models are based on ocean density layer knots obtained from real observation data and changes thereof simulate and forecast internal tide, however, depth, strength and thickness of density jump layers in different observation data are synchronously changed, and influence degree of each element (depth, thickness and strength) of the density jump layers on the internal tide cannot be quantitatively separated and extracted. Moreover, the current model can only forecast the influence of the change of the lower layer junction of the single topography on the intensity of the internal tide, and the topography of the internal tide source is complex and various, such as Liu Po, sea mountain, sea ridge and sea area all have obvious internal tide, the influence degree of the lower layer junction element of the complex topography on the internal tide can not be accurately captured by the single topography model, and the adaptability is required to be improved.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides an analysis method for the influence degree of density jump layer element change on internal tide.
In order to achieve the above object, the present invention is specifically achieved by the following techniques:
the invention provides an analysis method for the influence degree of density jump layer element change on internal tide, which comprises the following steps:
s1, selecting a type for simulating the topography of a marine area, and setting topography parameters and density jump parameters, wherein the density jump parameters comprise strength, depth and thickness, and calculating to obtain an initial density field based on the strength, the depth and the thickness data;
s2, applying trend velocity forcing in an open boundary, performing internal tide simulation, and outputting horizontal flow velocity, ocean temperature and ocean salinity distribution;
s3, calculating the distribution of the sea water density by using the output sea temperature and the sea salinity;
s4, changing the density jump layer parameters, and repeating the steps S2-S3;
s5, analyzing the influence of the density jump parameter change on the internal tide intensity and/or the ocean density according to the output horizontal flow velocity and the calculated distribution difference of the sea water density.
Further, in step S1, the types of the marine region topography include Liu Po topography, mountain topography and ridge topography.
Still further, the Liu Po terrain is represented by a first formula comprising:
；
wherein h is slope (x) At Liu Po sea depth, H is seabed depth, H s At Liu Po shallowest depth, L s Length of Liu Po, x s For the west boundary position, x is the east-west topographic coordinate and tanh represents the hyperbolic tangent function.
Still further, the sea mountain terrain and the ridge terrain are represented by a second formula comprising:
；
wherein H (x, y) is the depth of the terrain, H is the depth of the sea floor, H 0 For the maximum height of the topography, x 0 And y 0 The center coordinates of the terrains in the east-west direction and the north-south direction are respectively, x and y are terrains in the east-west direction and the north-south direction, and a and b are terrains in the east-west direction and the north-south direction; when a=b, for the sea mountain topography, when a<b, the sea ridge topography is obtained.
Further, in step S1, the initial density field is calculated according to a third formula, where the third formula includes:
；
wherein ρ is a (z) is the vertical density ρ 0 Is of a background densityDegree, z is the vertical direction topographic coordinates, Δρ a D and delta are the strength, depth and thickness of the dense jump layer, respectively.
Further, Δρ a Represented by the difference in density between the upper and lower interfaces of the skip-density layer.
Further, in step S2, the flow rate within the open boundary is calculated by a fourth formula, where the fourth formula includes:
；
wherein U is ob To open the boundary inner tide flow velocity, U am1 And U am2 Flow velocity amplitude, ω, of the first and second modes, respectively tide The inner tide frequency, t is time, pi is circumference rate, z is the vertical direction topographic coordinate, and h is the simulated area water depth.
Further, in step S2, the internal tide simulation is performed using MITgcm (MIT General Circulation Model) ocean numerical mode to output the horizontal flow rate, the ocean temperature and the ocean salinity profile.
Further, in step S3, a seawater kit in Matlab software is adopted, a sw_ptmp command is called, and the sea water density distribution with time and depth is calculated by using the sea temperature and the sea salinity.
Further, in step S5, the influence of the change of the density jump parameter on the internal tide intensity and/or the ocean density is analyzed according to a fifth formula and/or a sixth formula;
the fifth formula includes:
intra-tide intensity influence degree = V B -V S ；
Wherein V is B For horizontal flow rate of internal tide analog output under initial density jump layer parameter, V s The horizontal flow rate is outputted in an internal tide simulation mode after the density jump layer parameters are changed;
the sixth formula includes:
ocean density influence degree = ρ B -ρ S ；
Wherein ρ is B To calculate the sea water density at a certain depth ρ by internal tide simulation under the initial density jump parameter S The sea water density at the same depth is obtained by internal tide simulation calculation after changing the density jump layer parameter.
Compared with the prior art, the invention has the advantages that:
according to the invention, through parameterizing the junction elements (depth, thickness and strength) of each layer of the density jump layer and analyzing the difference of the internal tide intensity and the ocean density under different density jump layer parameters, the influence degree of different elements of the density jump layer on the internal tide intensity and the ocean density distribution can be quantitatively analyzed, and the method has important significance on exploring the specific influence of the change of the junction elements of the density jump layer on the internal tide.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for analyzing the effect of changes in the density jump factor on the internal tide in accordance with an embodiment of the present invention;
FIG. 2 is a graph showing the response of the horizontal flow rate and the internal tide intensity to the change of the density jump factor in the embodiment 1 of the present invention, wherein (a) and (d) are respectively an east-west horizontal flow rate graph and a north-south horizontal flow rate graph which are obtained by performing internal tide simulation under the initial density jump parameter, (b) and (e) are respectively an east-west horizontal flow rate graph and a north-south horizontal flow rate graph which are obtained by performing internal tide simulation after changing the density jump parameter, and (c) and (f) are respectively an influence degree graph of the change of the density jump parameter on the east-west direction internal tide intensity and an influence degree graph of the north-south internal tide intensity;
fig. 3 is a graph showing the results of the sea water density distribution and the response of the sea density to the change of the density jump factor in example 1 of the present invention, wherein (a) is the sea water density distribution obtained by performing the internal tide simulation calculation under the initial density jump parameter, (b) is the sea water density distribution obtained by performing the internal tide simulation calculation after changing the density jump parameter, and (c) is the influence degree graph of the change of the density jump parameter on the sea density.
Detailed Description
It should be noted that, without conflict, the embodiments of the present invention and features of the embodiments may be combined with each other. In addition, the terms "comprising," "including," "having," and "containing" are not limiting, as other steps and other ingredients may be added that do not affect the result. Materials, equipment, reagents are commercially available unless otherwise specified.
For a better understanding of the present invention, and not to limit its scope, all numbers expressing quantities, percentages, and other values used in the present invention are to be understood as being modified in all instances by the term "about". Accordingly, unless indicated otherwise, the numerical parameters set forth in the specification and claims are approximations that may vary depending upon the desired properties sought to be obtained. Each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.
In order that the above objects, features and advantages of the invention will be readily understood, a more particular description of the invention will be rendered by reference to specific embodiments thereof which are illustrated in the appended drawings.
The embodiment of the invention provides an analysis method for the influence degree of density jump layer element change on internal tide, which is shown in fig. 1 and comprises the following steps:
s1, selecting a type for simulating the topography of an ocean area, setting topography parameters and density jump parameters, wherein the density jump parameters comprise strength, depth and thickness, and calculating to obtain an initial density field based on the strength, depth and thickness data;
s2, applying trend velocity forcing in an open boundary, performing internal tide simulation, and outputting horizontal flow velocity, ocean temperature and ocean salinity distribution;
and (3) performing internal tide simulation by applying low-mode internal tide speed forcing at the open boundary for driving the simulation area to generate internal tide. The initial density field distribution and open boundary internal tide flow rate are the initial fields necessary to simulate internal tide. After the initial density field distribution is obtained through calculation, the internal tide is forced by applying the internal tide velocity of the open boundary, the internal tide is simulated by utilizing a MITgcm (MIT General Circulation Model) ocean numerical mode to output horizontal flow velocity, ocean temperature and ocean salinity distribution, wherein the horizontal flow velocity comprises an east-west horizontal flow velocity u and a north-south horizontal flow velocity v, the larger the absolute value of the flow velocity is, the stronger the internal tide is, and the internal tide can be used for representing the internal tide intensity.
S3, calculating the distribution of the sea water density by using the output sea temperature and sea salinity;
the specific operation comprises the following steps: and (3) invoking a sw_ptmp command by adopting a seawater kit in Matlab software, and calculating by utilizing ocean temperature and ocean salinity to obtain the distribution of the sea water density along with time and depth.
S4, changing the density jump layer parameters, and repeating the steps S2-S3;
s5, analyzing the influence of the density jump layer parameter change on the internal tide intensity and/or the ocean density according to the output horizontal flow velocity and the calculated distribution difference of the sea water density.
Specifically, the degree of influence on the internal tide intensity is calculated according to a fifth formula, which includes:
intra-tide intensity influence degree = V B -V S ；
Wherein V is B For horizontal flow rate of internal tide analog output under initial density jump layer parameter, V s To provide a horizontal flow rate for the endotide analog output after changing the skip-to-density parameters. The flow rate is one of means for indicating the intensity of the internal tide, and the larger the flow rate difference is, the larger the influence on the intensity of the internal tide is.
Calculating the degree of influence on the ocean density according to a sixth formula comprising:
ocean density influence degree = ρ B -ρ S ；
Wherein ρ is B To calculate the sea water density at a certain depth ρ by internal tide simulation under the initial density jump parameter S To change at the same timeAnd (5) carrying out internal tide simulation calculation after the density jump layer parameters to obtain the sea water density at the same depth. Seawater density is one of the means for representing the distribution of ocean density, and a larger difference in seawater density indicates a larger influence on the ocean density distribution.
The method comprises the steps of firstly setting initial values of depth, strength and thickness of a density jump layer, establishing an internal tide model for simulation, taking the initial values as a reference experiment, outputting simulation results including east-west horizontal flow velocity u, north-south horizontal flow velocity v, ocean temperature T, ocean salinity s and the like, and calculating the distribution of sea water density rho along with time and depth based on ocean temperature T and ocean salinity s data. And then, according to the needs of a user, resetting the density jump layer parameters (depth, thickness and strength) as a sensitive experiment to perform internal tide simulation, obtaining the simulation results of the other group of horizontal flow rates u and v and the ocean temperature T and the ocean salinity s, calculating to obtain the new distribution of the sea water density rho along with time and depth, and taking the difference between the two groups of results, so that the influence degree of the density jump layer parameter change on the internal tide intensity and the ocean density can be quantitatively analyzed. The internal tide intensity is represented by the horizontal flow rates u and v, and the ocean density is represented by the sea water density ρ distribution. The degree of influence of the internal tide intensity is quantified by the difference between the flow rates of the reference experiment and the sensitive experiment, and the degree of influence of the ocean density is quantified by the difference between the densities of the seawater at the same depth between the reference experiment and the sensitive experiment.
According to the invention, through parameterizing the density jump layer lamination element and analyzing the difference of the internal tide intensity and the ocean density under different density jump layer parameters, the influence degree of different elements (depth, thickness and intensity) of the density jump layer on the internal tide intensity and the ocean density distribution can be quantitatively analyzed, and the method has important significance on exploring the influence of the layer lamination change at the density jump layer on the internal tide.
In step S1, the types of marine region topography include Liu Po topography, mountain topography, and ridge topography. Generally, other areas than the aforementioned terrain are provided as flat-bottomed seas of uniform depth, the specific depth being specified by the user, for example, 1000 meters, 2000 meters.
Specifically, liu Po terrain is represented by a first formula comprising:
；
wherein h is slope (x) At Liu Po sea depth, H is seabed depth, H s At Liu Po shallowest depth, L s Length of Liu Po, x s For the west boundary position, x is the east-west topographic coordinate and tanh represents the hyperbolic tangent function. H. h is a s And L s Can be set autonomously by the user.
The terrain Liu Po simulated by the invention is three-dimensional space (x, y, z), and the depth of the width (the terrestrial coordinates in the north-south direction) y of Liu Po is not changed along with the change of the y coordinate, so that the width of the whole simulated area, namely the width of a sea basin, is not required to be obtained by the width y of Liu Po in the formula.
Specifically, the sea mountain terrain and the sea ridge terrain are expressed by a second formula including:
；
wherein H (x, y) is the depth of the terrain, H is the depth of the sea floor, H 0 For the maximum height of the topography, x 0 And y 0 The center coordinates of the terrains in the east-west direction and the north-south direction are respectively, x and y are terrains in the east-west direction and the north-south direction, and a and b are terrains in the east-west direction and the north-south direction; when a=b, the sea mountain topography is the sea mountain topography, when a<b, the sea ridge topography. h is a 0 The a and b can be set autonomously by the user.
The Liu Po topography is expressed by hyperbolic functions, the sea mountain topography or the sea ridge topography is expressed by two-dimensional high-silk functions, and the two-dimensional distribution (comprising the height and depth of the topography) of the simulated topography can be given for three-dimensional simulation of the internal tide; in addition, two different terrains (sea mountain: a=b, sea ridge: a < b) of sea mountain and sea ridge can be represented by setting parameters a and b, so that the internal tide of the sea area of Liu Po, sea mountain and sea ridge can be simulated respectively, and important parameters such as height and width of terrains such as Liu Po, sea mountain and sea ridge are set independently by a user according to actual conditions, so that the influence of the parameter change of the density jump layer on the internal tide intensity under different complex terrains can be predicted, and the adaptability is strong.
In step S1, the initial density field is calculated according to a third formula, where the third formula includes:
；
wherein ρ is a (z) is the vertical density ρ 0 For background density, z is the vertical direction topographic coordinate, Δρ a D and delta are the strength, depth and thickness of the dense jump layer, respectively; Δρ a Represented by the difference in density between the upper and lower interfaces of the density jump layer, is used to characterize the strength of the density jump layer.
It should be noted that, the user can set the intensity, depth and thickness of the densitometry layer according to specific needs, the intensity, depth and thickness can be changed simultaneously or only one or two parameters can be changed, the density of different depths is calculated by using the above formula, the data is used as the initial density field distribution in the internal tide simulation, and the junction elements of each densitometry layer are further parameterized.
Considering that the internal tide in the global ocean is mainly in low modes (first and second modes), the open boundary internal tide speed in step S2 is forced to mainly apply the first and second modes for driving the model simulation, and the specific open boundary internal tide speed is calculated by a fourth formula, which includes:
；
wherein U is ob To open the boundary inner tide flow velocity, U am1 And U am2 Flow velocity amplitude, ω, of the first and second modes, respectively tide The inner tide frequency, t is time, pi is circumference rate, z is the vertical direction topographic coordinate, and h is the simulated area water depth. U (U) am1 And U am2 Can be set autonomously by the user.
The invention will be further illustrated with reference to specific examples. It is to be understood that these examples are illustrative of the present invention and are not intended to limit the scope of the present invention. The experimental methods, which do not address specific conditions in the following examples, are generally in accordance with the conditions recommended by the manufacturer.
Example 1
The embodiment of the invention provides an analysis method for the influence degree of density jump layer element change on internal tide, which is shown in fig. 1 and comprises the following steps:
s1, selecting a type for simulating the topography of a marine area, and setting topography parameters and density jump layer parameters;
the topography of the sea bottom in this example is selected as the topography of the sea ridge; the other areas except the terrain are flat-bottomed seas with the water depth of 1000 m;
the ridge topography is expressed by the following formula:
；
wherein H (x, y) is the depth of the terrain, H is the depth of the sea bottom, and is specifically 1000m, H 0 For maximum height of topography, in particular 800m, x 0 And y 0 The center coordinates of the terrains in the east-west direction and the north-south direction are respectively, the specific values are (121.5 DEG E,20.5 DEG N), x and y are terrains in the east-west direction and the north-south direction respectively, and a and b are terrains in the east-west direction and the north-south direction respectively; a=100 km, b=550 km;
the density jump layer parameters comprise intensity, depth and thickness, and based on the intensity, depth and thickness data, an initial density field is calculated according to the following formula and is used for subsequent internal tide simulation;
；
wherein ρ is a (z) is the vertical density ρ 0 At a background density of, in particular, 1025kg/m 3 Z is the vertical direction topographic coordinates, Δρ a D and delta are the strength, depth and thickness of the density jump layer, respectively, specifically Δρ a =3 kg/m 3 、d=300m、δ=200m；
S2, applying flow speed forcing in an open boundary, performing internal tide simulation by using an MITgcm ocean numerical mode, and outputting horizontal flow velocity, ocean temperature and ocean salinity distribution, wherein the horizontal flow velocity comprises an east-west horizontal flow velocity u and a north-south horizontal flow velocity v;
the flow rate within the open boundary is calculated by the following formula:
；
wherein U is ob To open the boundary inner tide flow velocity, U am1 And U am2 The flow velocity amplitudes of the first mode and the second mode are respectively, and the flow velocity amplitudes of the first mode and the second mode are respectively in the east-west direction: u (U) am1 =0.2 m/s and U am2 =0.1 m/s for north-south direction: u (U) am1 =0.1 m/s and U am2 =0.05 m/s，ω tide For the intra-tide frequency, M is taken in this embodiment 2 An internal tide frequency of the internal tide;
s3, using a seawater toolkit in Matlab software to call a sw_ptmp command, and calculating the distribution of the sea water density along with time and depth by using the output sea temperature and sea salinity;
s4, changing the density jump layer parameters, and setting Deltaρ a =6 kg/m 3 Repeating steps S2-S3, d=200m and δ=300m;
s5, calculating the influence degree of the internal tide intensity and the influence degree of the sea density according to the following formula according to the internal tide flow velocity output before and after changing the density jump layer parameter and the calculated distribution difference of the sea water density;
intra-tide intensity influence degree = V B -V S ；
Wherein V is B For horizontal flow rate of internal tide analog output under initial density jump layer parameter, V s The horizontal flow rate is outputted in an internal tide simulation mode after the density jump layer parameters are changed;
ocean density influence degree = ρ B -ρ S ；
Wherein ρ is B To calculate the sea water density at a certain depth ρ by internal tide simulation under the initial density jump parameter S To be carried out after changing the parameters of the skip-density layerThe sea water density at the same depth was calculated by the internal tide simulation.
Fig. 2 shows the influence of the east-west direction inner tide intensity (see (c) in fig. 2) and the north-south direction inner tide intensity (see (f) in fig. 2) on the east-west direction inner tide intensity (see (a) in fig. 2) and the north-south direction inner tide intensity (see (f) in fig. 2) of the east-west direction inner tide simulation output under the initial density jump parameter, the east-west direction inner tide simulation output after changing the density jump parameter, and the north-south direction inner tide simulation output.
Fig. 3 shows the sea water density distribution obtained by performing the internal tide simulation calculation under the initial density jump parameter (see (a) in fig. 3), the sea water density distribution obtained by performing the internal tide simulation calculation after changing the density jump parameter (see (b) in fig. 3), and the degree of influence of the change of the density jump on the sea density (see (c) in fig. 3).
As can be seen from fig. 2-3, the method can quantitatively analyze the influence degree of different elements of the density jump layer on the internal tide intensity and the ocean density by parameterizing the junction elements (depth, thickness and intensity) of each layer of the density jump layer and analyzing the difference of the internal tide intensity and the ocean density under different density jump layer parameters, and can rapidly explore the specific influence of the junction element change on the internal tide.
Although the present disclosure is disclosed above, the scope of the present disclosure is not limited thereto. Various changes and modifications may be made by one skilled in the art without departing from the spirit and scope of the disclosure, and these changes and modifications will fall within the scope of the disclosure.

Claims (9)

1. A method for analyzing the degree of influence of a density jump factor change on an internal tide, comprising the steps of:
s1, selecting a type for simulating the topography of a marine area, and setting topography parameters and density jump parameters, wherein the density jump parameters comprise strength, depth and thickness, and calculating to obtain an initial density field based on the strength, the depth and the thickness data;
s2, applying trend velocity forcing in an open boundary, performing internal tide simulation, and outputting horizontal flow velocity, ocean temperature and ocean salinity distribution;
s3, calculating the distribution of the sea water density by using the output sea temperature and the sea salinity;
s4, changing the density jump layer parameters, and repeating the steps S2-S3;
s5, analyzing the influence of the density jump parameter change on the internal tide intensity and/or the ocean density according to the output horizontal flow velocity and the calculated distribution difference of the sea water density, wherein the method comprises the following steps:
analyzing the influence of the change of the density jump parameter on the internal tide intensity and/or the ocean density according to a fifth formula and/or a sixth formula;
the fifth formula includes:
intra-tide intensity influence degree = V B -V S ；
Wherein V is B For horizontal flow rate of internal tide analog output under initial density jump layer parameter, V s The horizontal flow rate is outputted in an internal tide simulation mode after the density jump layer parameters are changed;
the sixth formula includes:
ocean density influence degree = ρ B -ρ S ；
Wherein ρ is B To calculate the sea water density at a certain depth ρ by internal tide simulation under the initial density jump parameter S The sea water density at the same depth is obtained by internal tide simulation calculation after changing the density jump layer parameter.
2. The method according to claim 1, wherein in step S1, the types of the marine land topography include Liu Po topography, mountain topography and ridge topography.
3. The method of analyzing the degree of influence of changes in a density jump factor on an internal tide of claim 2, wherein the Liu Po topography is represented by a first formula comprising:
；
wherein h is slope (x) At Liu Po sea depth, H is seabed depth, H s At Liu Po shallowest depth, L s Length of Liu Po, x s For the west boundary position, x is east-west terrain coordinates.
4. The method of analyzing the degree of influence of changes in a density jump factor on an internal tide according to claim 2, wherein said sea mountain topography and said sea ridge topography are represented by a second formula comprising:
；
wherein H (x, y) is the depth of the terrain, H is the depth of the sea floor, H 0 For the maximum height of the topography, x 0 And y 0 The center coordinates of the terrains in the east-west direction and the north-south direction are respectively, x and y are terrains in the east-west direction and the north-south direction, and a and b are terrains in the east-west direction and the north-south direction; when a=b, for the sea mountain topography, when a<b, the sea ridge topography is obtained.
5. The method of claim 1, wherein in step S1, the initial density field is calculated according to a third formula, the third formula comprising:
；
wherein ρ is a (z) is the vertical density ρ 0 For background density, z is the vertical direction topographic coordinate, Δρ a D and delta are the strength, depth and thickness of the dense jump layer, respectively.
6. The method for analyzing the influence degree of a change in a density jump factor on an internal tide according to claim 5, wherein Δρ a Represented by the difference in density between the upper and lower interfaces of the skip-density layer.
7. The method according to claim 1, wherein in step S2, the open boundary internal flow rate is calculated by a fourth formula, the fourth formula including:
；
wherein U is ob To open the boundary inner tide flow velocity, U am1 And U am2 Flow velocity amplitude, ω, of the first and second modes, respectively tide The inner tide frequency, t is time, pi is circumference rate, z is the vertical direction topographic coordinate, and h is the simulated area water depth.
8. The method according to claim 1, wherein in step S2, the internal tide simulation is performed using MITgcm ocean numerical mode.
9. The method according to claim 1, wherein in step S3, a seawater kit in Matlab software is used to call a sw_ptmp command, and the sea water density distribution over time and depth is calculated using the sea temperature and the sea salinity.
CN202310691642.XA 2023-06-13 2023-06-13 Analysis method for influence degree of density jump layer element change on internal tide Active CN116432562B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310691642.XA CN116432562B (en) 2023-06-13 2023-06-13 Analysis method for influence degree of density jump layer element change on internal tide

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310691642.XA CN116432562B (en) 2023-06-13 2023-06-13 Analysis method for influence degree of density jump layer element change on internal tide

Publications (2)

Publication Number Publication Date
CN116432562A CN116432562A (en) 2023-07-14
CN116432562B true CN116432562B (en) 2023-08-22

Family

ID=87085838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310691642.XA Active CN116432562B (en) 2023-06-13 2023-06-13 Analysis method for influence degree of density jump layer element change on internal tide

Country Status (1)

Country Link
CN (1) CN116432562B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150068924A (en) * 2014-10-21 2015-06-22 주식회사 지앤지테크놀러지 Tidal current type generator
CN110008509A (en) * 2019-03-01 2019-07-12 中国海洋大学 A kind of interior estimates force characteristics analysis method under consideration Background Flow Field

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150068924A (en) * 2014-10-21 2015-06-22 주식회사 지앤지테크놀러지 Tidal current type generator
CN110008509A (en) * 2019-03-01 2019-07-12 中国海洋大学 A kind of interior estimates force characteristics analysis method under consideration Background Flow Field

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
浅强跃层大振幅内波研究;杜涛;地学前缘(S2);全文 *

Also Published As

Publication number Publication date
CN116432562A (en) 2023-07-14

Similar Documents

Publication Publication Date Title
CN107341341B (en) A kind of river mouth point source sudden water pollution event source tracing method
EP2235500B1 (en) Determining connectivity architecture in 2-d and 3-d heterogeneous data
CN103941290B (en) A kind of submarine cable movement locus analogy method and system
CN102109495A (en) Method for classifying types of mixed seabed sediment based on multi-beam sonar technology
CN117198330B (en) Sound source identification method and system and electronic equipment
CN104180873B (en) Single-wave-beam depth finder water depth gross error detection and correction method and system
CN109855688A (en) A kind of inland harbour marine exhaust discharge Measurement Method
CN113420440B (en) Vortex identification method based on ocean vertical structure
CN112270113A (en) Seabed acoustic scattering FEM analysis method in laminar flow layered structure mode
CN105718725B (en) The Quantification Methods of polluted bed mud and its pollutant
CN108594241B (en) A kind of AUV Sound stealth method for situation assessment
Zhang et al. An online path planning algorithm for autonomous marine geomorphological surveys based on AUV
CN116432562B (en) Analysis method for influence degree of density jump layer element change on internal tide
CN110456024A (en) A kind of method and system for analyzing gas hydrates stable region boundary carbon cycle process
CN114325851A (en) Method and system for detecting underwater ferromagnetic target based on multi-frequency radiation source
CN110046374A (en) A kind of calculation method of sound field based on Gaussian acoustic beam
CN104680583B (en) A kind of method that sea-floor relief automatically generates
Tang et al. Flow characteristics of open-channel flow with partial two-layered vegetation
EP2098887A1 (en) Determining connectivity architecture in 2-D and 3-D heterogeneous data
CN114325853A (en) Method for analyzing amplitude-phase-frequency characteristics of sequence time window of disturbance of power frequency electromagnetic field of underwater vehicle
CN110096745B (en) Wave abundance depth calculation method based on wave encounter period
CN113221651A (en) Seafloor sediment classification method using acoustic propagation data and unsupervised machine learning
Fan et al. Research on shipborne aided navigation system based on enhanced traffic environment perception
CN207318714U (en) A kind of marine traffic engineering Seaworthiness bathymetric survey integrated processing unit
CN111751882A (en) Method and device for generating sensitive rock physical parameters

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant