Ocean vortex three-dimensional flow field reconstruction method based on physical constraint

Abstract

The invention discloses an ocean vortex three-dimensional flow field reconstruction method based on physical constraints, and belongs to the technical field of physical oceanography and data modeling, and the method comprises the steps: calculating sea surface height abnormal data according to climate state data and sea surface observation data, recognizing vortex features, and dividing a vortex into grid points; according to the climate state data and the vortex vertical mode characteristic equation, constructing a vortex vertical mode function, and obtaining a vortex horizontal flow velocity vertical structure function; by taking the sea surface flow velocity abnormal condition and the deep layer flow velocity abnormal condition as constraints, determining a vortex vertical modal coefficient according to a vortex horizontal flow velocity vertical structure function, and reconstructing a vortex horizontal flow velocity vertical structure at each grid point; and superposing the vortex horizontal flow velocity vertical structure and the climate state horizontal flow velocity vertical structure to obtain a seawater horizontal flow velocity vertical structure to construct a vortex three-dimensional flow field structure. The problems that the prior art is limited to the sea surface, and the vertical structure of the vortex is difficult to accurately describe are solved.

Description

Technical Field

[0001] The present invention relates to a method for reconstructing the three-dimensional flow field of ocean vortices based on physical constraints, belonging to the technical fields of physical oceanography and data modeling. Background Art

[0002] In the study of ocean dynamics, vortices, as key factors affecting climate change, material transport, and ecosystem balance, have always been the focus of attention of scientific researchers. At present, the research on the three-dimensional structure of vortices mainly relies on in-situ measurement data, which are usually collected by equipment such as shipborne instruments, buoys, and moored buoys. In-situ measurement can directly obtain key parameters such as velocity, temperature, and salinity inside the vortices, providing valuable data support for case studies of vortices. However, this method is limited by factors such as the deployment range, operation time, and cost of observation equipment, and usually can only conduct short-term observations on vortices in specific regions, making it difficult to cover large-scale vortex activities or conduct long-term continuous monitoring. Therefore, the existing in-situ measurement technology has obvious limitations in studying large-scale and long-term vortex processes.

[0003] Although in-situ measurement data play an irreplaceable role in the study of the three-dimensional structure of vortices, their inherent limitations make it extremely difficult to comprehensively and deeply understand the characteristics of vortices. First of all, in-situ measurement is usually time-consuming, laborious, and costly, which greatly limits the frequency and scope of observations, making it a very challenging task to obtain the full three-dimensional structure of ocean vortices. Secondly, due to the limitations of observation means, existing research often can only conduct case studies on individual vortices, making it difficult to form a profound understanding of the general laws of vortices. In addition, with the rapid development of satellite remote sensing technology, although it provides a large amount of real-time and quasi-real-time information with a wide coverage range, high spatial resolution, and strong time continuity on the sea surface, most of these data are limited to the sea surface and are difficult to accurately depict the vertical structure of vortices, thus limiting the accuracy and reliability of the reconstruction of the three-dimensional flow field of vortices. Summary of the Invention

[0004] The purpose of the present invention is to provide a method for reconstructing the three-dimensional flow field of ocean vortices based on physical constraints, by reconstructing the three-dimensional flow field structure of vortices with the dynamic characteristics of vortices as physical constraint conditions, so as to solve the problem that the existing technology is limited to the sea surface and is difficult to accurately depict the vertical structure of vortices.

[0005] To solve the above technical problems, the present invention is implemented by adopting the following technical solutions:

[0006] The present invention provides a method for reconstructing the three-dimensional flow field of ocean vortices based on physical constraints, including:

[0007] Collecting the climatological data and sea surface observation data of the target sea area;

[0008] Calculate sea surface height anomaly data based on the climatological data and sea surface observation data;

[0009] Identify vortex features based on the sea surface height anomaly data, and divide the vortices into grids based on the spatial resolution of the sea surface height observation data according to the vortex features;

[0010] Construct a vertical mode function of the vortex based on the climatological data and the vertical mode characteristic equation of the vortex, and extract the first two vertical modes for superposition to obtain the vertical structure function of the horizontal velocity of the vortex;

[0011] Determine the vertical mode coefficients of the vortex by taking the sea surface velocity anomaly condition and the deep velocity anomaly condition as constraints respectively, according to the vertical structure function of the horizontal velocity of the vortex;

[0012] Reconstruct the vertical structure of the horizontal velocity of the vortex at each grid point according to the vertical mode function of the vortex and the vertical mode coefficients of the vortex;

[0013] Calculate the vertical structure of the climatological horizontal velocity according to the climatological data and the thermal wind relationship;

[0014] Superpose the vertical structure of the horizontal velocity of the vortex at each grid point and the vertical structure of the climatological horizontal velocity to obtain the vertical structure of the seawater horizontal velocity at each grid point;

[0015] Construct the three-dimensional flow field structure of the vortex according to the vertical structure of the seawater horizontal velocity at each grid point.

[0016] Further, the climatological data includes climatological temperature, climatological salinity, climatological flow field and climatological sea surface height;

[0017] The sea surface observation data includes observed sea surface height data;

[0018] Among them, calculating the sea surface height anomaly data based on the climatological data and the sea surface observation data includes:

[0019] Subtract the climatological sea surface height data from the observed sea surface height data to obtain the sea surface height anomaly data.

[0020] Further, the vortex features include the vortex center and the vortex boundary;

[0021] The vortex center includes the potential center of the anticyclonic vortex and the potential center of the cyclonic vortex;

[0022] Among them, identifying the vortex features based on the sea surface height anomaly data includes:

[0023] If the spatial resolution of the sea surface height anomaly data is 1 / 4º×1 / 4º, a 5×5 grid point window is scanned to find the extreme points, where the maximum point is the potential center of the anticyclonic vortex and the minimum point is the potential center of the cyclonic vortex;

[0024] Taking the sea surface height anomaly data at the vortex center as the starting value, the sea surface height anomaly change value is gradually increased or decreased outward until the outermost isoline only contains one vortex center, and the outermost isoline is the vortex boundary.

[0025] Furthermore, a vortex vertical mode function is constructed according to the climatological data and the vortex vertical mode characteristic equation, including:

[0026] Calculate the buoyancy frequency of the target sea area according to the climatological temperature and climatological salinity;

[0027] Construct a vortex vertical mode function according to the Coriolis parameter, the buoyancy frequency of the target sea area, the vertical coordinates of the grid points, and the vortex vertical mode characteristic equation.

[0028] Furthermore, the vortex vertical mode characteristic equation is expressed as:

[0029] ;

[0030] In the formula, represents differentiation, represents differentiation with respect to the vertical coordinate of the grid point, represents the eigenvalue, represents the Coriolis parameter, represents the buoyancy frequency of the target sea area, represents the th vertical mode of the vortex vertical mode function, represents the vertical coordinate of the grid point, that is, the seawater depth where the climatological temperature and climatological salinity of the target sea area are located.

[0031] Furthermore, the vortex horizontal velocity vertical structure function includes a vortex east-west velocity function and a vortex north-south velocity function;

[0032] Among them, the vortex east-west velocity function and the vortex north-south velocity function are respectively expressed as:

[0033] ;

[0034] In the formula, and respectively represent the vortex east-west velocity function and the vortex north-south velocity function, represents the vertical coordinate of the grid point, , , and represent the vertical mode coefficients of the vortex. Among them, and represent the east - west velocity coefficient and the north - south velocity coefficient of the first vertical mode of the vertical structure function of the vortex horizontal flow velocity respectively, of and represent the east - west velocity coefficient and the north - south velocity coefficient of the second vertical mode of the vertical structure function of the vortex horizontal flow velocity respectively, of

[0035] Furthermore, the calculation method of the sea - surface flow velocity anomaly condition includes:

[0036] Calculate through observing sea - surface height data and geostrophic relationship to obtain the sea - surface geostrophic current;

[0037] Calculate through climatological sea - surface height data and geostrophic relationship, the climatological sea - surface geostrophic current;

[0038] Subtract the climatological sea - surface geostrophic current from the sea - surface geostrophic current to obtain the sea - surface flow velocity anomaly condition;

[0039] The calculation method of the deep - layer flow velocity anomaly condition includes:

[0040] Set the horizontal flow velocity of the vortex at a depth of 2000 meters in the target sea area to 0 to obtain the deep - layer flow velocity anomaly condition.

[0041] Furthermore, under the sea - surface flow velocity anomaly condition, the vortex east - west velocity function and the vortex north - south velocity function are respectively expressed as:

[0042] ;

[0043] In the formula, , represent the vortex east - west velocity function and the vortex north - south velocity function respectively under the sea - surface flow velocity anomaly condition, , represent the first vertical mode and the second vertical mode of the vertical structure function of the vortex horizontal flow velocity respectively under the sea - surface flow velocity anomaly condition, , , and represent the vertical mode coefficients of the vortex. Among them, and represent respectively of the east - west velocity coefficient and the north - south velocity coefficient, and represent respectively The east-west flow velocity coefficient and the north-south flow velocity coefficient;

[0044] Under the condition of abnormal deep-layer flow velocity, the vortex east-west flow velocity function and the vortex north-south flow velocity function are respectively expressed as:

[0045] ;

[0046] In the formula, , respectively represent the vortex east-west flow velocity function and the vortex north-south flow velocity function under the condition of abnormal deep-layer flow velocity, , respectively represent the first vertical mode and the second vertical mode of the vertical structure function of the vortex horizontal flow velocity under the condition of abnormal deep-layer flow velocity.

[0047] Furthermore, the calculation method of the vortex vertical mode coefficient includes:

[0048] Calculating the vortex vertical mode coefficient according to the vortex east-west flow velocity function and the vortex north-south flow velocity function under the condition of abnormal sea surface flow velocity and the vortex east-west flow velocity function and the vortex north-south flow velocity function under the condition of abnormal deep-layer flow velocity.

[0049] Furthermore, the vertical structure of the seawater horizontal flow velocity includes the seawater east-west horizontal flow velocity and the seawater north-south horizontal flow velocity;

[0050] The vertical structure of the vortex horizontal flow velocity includes the vortex east-west horizontal flow velocity and the vortex north-south horizontal flow velocity;

[0051] The vertical structure of the climatological mean horizontal flow velocity includes the climatological mean east-west horizontal flow velocity and the climatological mean north-south horizontal flow velocity;

[0052] The seawater east-west horizontal flow velocity and the seawater north-south horizontal flow velocity are respectively expressed as:

[0053] ;

[0054] In the formula, , respectively represent the seawater east-west horizontal flow velocity and the seawater north-south horizontal flow velocity, , respectively represent the vortex east-west horizontal flow velocity and the vortex north-south horizontal flow velocity, and respectively represent the climatological mean east-west horizontal flow velocity and the climatological mean north-south horizontal flow velocity, represents the vertical coordinate of the grid point.

[0055] Compared with the prior art, the beneficial effects achieved by the present invention:

[0056] 1. The present invention realizes the reconstruction of the three-dimensional flow field of ocean vortices efficiently and accurately by integrating climatological data, sea surface observation data, and using the vortex dynamics characteristics as physical constraint conditions. The present invention not only makes full use of the advantages of high spatio-temporal resolution of satellite remote sensing data but also combines climatological historical data, effectively making up for the deficiencies of a single data source. By constructing the vortex vertical mode function and extracting the first two vertical modes for superposition, and determining the vortex vertical mode coefficients with the sea surface velocity anomaly condition and the deep velocity anomaly condition as constraints, only the observed sea surface height data is required to quickly realize vortex feature recognition, vortex vertical structure construction, and high-precision three-dimensional flow field reconstruction of vortices, solving the problem that the prior art is limited to the sea surface and difficult to accurately depict the vertical structure of vortices.

[0057] 2. The present invention integrates climatological data including temperature, salinity, flow field, and sea surface height with sea surface observation data, and calculates sea surface height anomaly data to identify vortex features, and then constructs the vortex vertical mode function. It not only makes full use of the complementary advantages of multi-source data but also ensures the accuracy of the three-dimensional flow field reconstruction of vortices through physical constraints such as the Coriolis parameter and buoyancy frequency. The present invention can accurately depict the vertical structure characteristics of vortices, providing new details and accuracy for ocean vortex research.

[0058] 3. The present invention determines the vortex center by scanning the sea surface height anomaly data to find the extreme points, including anticyclonic vortices and cyclonic vortices, and defines the vortex boundary based on the sea surface height anomaly change value. It not only improves the efficiency and accuracy of vortex identification but also helps to deeply understand the dynamic behavior and influence range of vortices. By accurately defining the vortex boundary, the present invention provides a more reliable basis for the reconstruction of the three-dimensional flow field of vortices.

[0059] 4. The present invention innovatively combines the sea surface velocity anomaly condition and the deep velocity anomaly condition as physical constraints to determine the vortex vertical mode coefficients. The sea surface geostrophic flow is calculated by observing the sea surface height data and subtracted from the climatological sea surface geostrophic flow to obtain the sea surface velocity anomaly condition. At the same time, the horizontal velocity of the vortex at a depth of 2000 meters in the target sea area is set to 0 as the deep velocity anomaly condition. It not only considers the velocity change characteristics of vortices in different depth layers but also ensures the consistency and reliability of the three-dimensional flow field reconstruction of vortices through physical constraints. By integrating the sea surface velocity anomaly condition and the deep velocity anomaly condition, the present invention can more accurately reconstruct the three-dimensional flow field structure of vortices, providing strong support for fields such as ocean environmental monitoring, resource development, and disaster warning. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] Figure 1 is a schematic flowchart of a method for reconstructing the three-dimensional flow field of ocean vortices based on physical constraints provided by an embodiment of the present invention;

[0061] Figure 2 It is a schematic plan view of the vortex features identified based on sea surface height anomaly and each grid point constructed in the embodiment of the present invention;

[0062] Figure 3 It is a schematic plan view of the horizontal velocity field at a depth of 100 meters in the vortex three-dimensional flow field reconstruction structure provided by the embodiment of the present invention;

[0063] Figure 4 It is a schematic plan view of the horizontal velocity field at a depth of 100 meters in the vortex of the ARMOR3D data provided by the embodiment of the present invention;

[0064] Figure 5 It is a schematic plan view of the horizontal velocity field at a depth of 300 meters in the vortex three-dimensional flow field reconstruction structure provided by the embodiment of the present invention;

[0065] Figure 6 It is a schematic plan view of the horizontal velocity field at a depth of 300 meters in the vortex of the ARMOR3D data provided by the embodiment of the present invention. Detailed implementation manners

[0066] The technical solution of the present invention will be described in detail below through the accompanying drawings and specific embodiments. It should be understood that the embodiments of the present invention and the specific features in the embodiments are detailed descriptions of the technical solution of the present invention, rather than limitations on the technical solution of the present invention. Without conflict, the technical features in the embodiments of the present invention and the embodiments can be combined with each other.

[0067] The term "and / or" only describes the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone. In addition, the character " / " generally represents an "or" relationship between the associated objects before and after.

[0068] Embodiment 1

[0069] As Figure 1 shown, this embodiment introduces a method for reconstructing the three-dimensional flow field of ocean vortices based on physical constraints, including:

[0070] Step 1: Collect the climatological data and sea surface observation data of the target sea area.

[0071] The climatological data provides the long-term average ocean environmental state of the target sea area, including temperature, salinity, flow field, and sea surface height. The climatological data is the basis for understanding the background environment of ocean vortices. The sea surface observation data provides real-time sea surface height information for capturing the dynamic changes of ocean phenomena such as vortices. The present invention prepares for analyzing and reconstructing the three-dimensional flow field of vortices by collecting the climatological data and sea surface observation data of the target sea area.

[0072] Step 2: Calculate the sea surface height anomaly data based on the climatological data and the sea surface observation data.

[0073] The sea surface height anomaly data is the deviation relative to the climatological sea surface height, reflecting the dynamic changes of the ocean surface, especially the sea surface undulations caused by ocean phenomena such as vortices. By calculating the sea surface height anomaly data, the present invention can more easily identify the position and boundary of vortices.

[0074] Step 3: Identify the vortex characteristics based on the sea surface height anomaly data, and divide the vortices into grids according to the spatial resolution of the sea surface height observation data based on the vortex characteristics.

[0075] The extreme points in the sea surface height anomaly data correspond to the vortex centers. The positions of the vortex centers can be determined by identifying the maximum and minimum points. Among them, the present invention uses the spatial resolution of the sea surface height observation data provided by the ocean satellite remote sensing data AVISO to divide the vortices into grids, and the ocean satellite remote sensing data AVISO determines the fineness of vortex identification. The present invention divides the vortices into grids in order to process and manage a large amount of data when analyzing and reconstructing the three-dimensional flow field of vortices subsequently.

[0076] Step 4: Construct the vortex vertical mode function according to the climatological data and the vortex vertical mode characteristic equation, and extract the first two vertical modes for superposition to obtain the vertical structure function of the vortex horizontal velocity.

[0077] The vortex vertical mode characteristic equation describes the structural characteristics of the vortex in the vertical direction, and the characteristic function can be derived from climatological data such as temperature and salinity and the vortex dynamics theory. The present invention extracts the first two vertical modes of the vortex vertical mode function for superposition to obtain the vertical structure function of the vortex horizontal velocity. Among them, the first two vertical modes can usually capture the main characteristics of the vortex vertical structure, thus simplifying the problem and improving the calculation efficiency. The vertical structure function of the vortex horizontal velocity is the key to determining the vortex vertical mode coefficients and reconstructing the three-dimensional flow field of the vortex subsequently.

[0078] Step 5: Determine the vortex vertical mode coefficients according to the vertical structure function of the vortex horizontal velocity, with the sea surface velocity anomaly condition and the deep sea velocity anomaly condition as constraints respectively.

[0079] The sea surface velocity anomaly condition and the deep sea velocity anomaly condition provide the velocity constraints of the vortex in the vertical direction. Taking the sea surface velocity anomaly condition and the deep sea velocity anomaly condition as physical constraints helps to more accurately determine the vortex vertical mode coefficients. The vortex vertical mode coefficients are the coefficients of the vortex vertical structure function, and the coefficients of the vortex vertical structure function determine the velocity distribution of the vortex in the vertical direction. By combining the sea surface velocity anomaly condition and the deep sea velocity anomaly condition, the present invention further improves the accuracy of the three-dimensional flow field reconstruction of the vortex.

[0080] Step Six: Reconstruct the vertical structure of the vortex horizontal velocity at each grid point according to the vortex vertical mode function and the vortex vertical mode coefficient.

[0081] After determining the vortex vertical mode function and the vortex vertical mode coefficient, the present invention transforms the three-dimensional vortex flow field from a theoretical model into actual data, and reconstructs the vertical structure of the vortex horizontal velocity at each grid point, providing a basis for subsequent superposition of the vertical structure of the climatological horizontal velocity and construction of the three-dimensional vortex flow field structure.

[0082] Step Seven: Superpose the vertical structure of the vortex horizontal velocity and the vertical structure of the climatological horizontal velocity at each grid point to obtain the vertical structure of the seawater horizontal velocity at each grid point.

[0083] The present invention calculates the vertical structure of the climatological horizontal velocity according to the climatological data and the thermal wind relationship.

[0084] The vertical structure of the climatological horizontal velocity represents the long-term average velocity distribution in the target sea area, while the vertical structure of the vortex horizontal velocity reflects the velocity changes caused by ocean phenomena such as vortices. By superposing the vertical structure of the vortex horizontal velocity and the vertical structure of the climatological horizontal velocity at each grid point, the present invention can obtain the actual vertical structure of the seawater horizontal velocity at each grid point.

[0085] Step Eight: Construct the three-dimensional vortex flow field structure according to the vertical structure of the seawater horizontal velocity at each grid point.

[0086] After obtaining the vertical structure of the seawater horizontal velocity at each grid point, the present invention can combine these data through methods such as spatial interpolation to construct a complete three-dimensional vortex flow field structure. The present invention provides a comprehensive understanding and visual representation of the three-dimensional vortex flow field, providing strong support for fields such as marine scientific research, environmental monitoring, and resource development.

[0087] Embodiment 2

[0088] Based on the same inventive concept as Embodiment 1, this embodiment introduces the implementation steps of a method for reconstructing the three-dimensional vortex flow field in the ocean based on physical constraints, including:

[0089] Step 1: Collect the climatological data and sea surface observation data of the target sea area.

[0090] In some embodiments, the climatological data includes climatological temperature, climatological salinity, climatological flow field, and climatological sea surface height; the climatological data can be extracted from existing climatological databases such as the World Ocean Atlas WOA, the three-dimensional ocean temperature and salinity dataset ARMOR3D, and the ocean satellite remote sensing data AVISO, or can be calculated from other reanalysis data through monthly or annual averages.

[0091] In some embodiments, the sea surface observation data includes sea surface height observation data, and in this embodiment, the sea surface observation data is extracted from AVISO data.

[0092] Step 2: Calculate the sea surface height anomaly data according to the climatological data and the sea surface observation data.

[0093] In some embodiments, calculating the sea surface height anomaly data according to the climatological data and the sea surface observation data includes:

[0094] Subtract the climatological sea surface height data from the observed sea surface height data to obtain the sea surface height anomaly data.

[0095] Step 3: Identify the vortex features according to the sea surface height anomaly data, and divide the vortices into grids based on the spatial resolution of the sea surface height observation data according to the vortex features.

[0096] The vortex features include the vortex center and the vortex boundary; the vortex center includes the potential center of the anticyclonic vortex and the potential center of the cyclonic vortex.

[0097] In some embodiments, identifying the vortex features according to the sea surface height anomaly data includes:

[0098] If the spatial resolution of the sea surface height anomaly data is 1 / 4º×1 / 4º, scan with a 5×5 grid point window to find the extreme points, where the maximum point is the potential center of the anticyclonic vortex and the minimum point is the potential center of the cyclonic vortex;

[0099] Taking the sea surface height anomaly data of the vortex center as the starting value, gradually increase or decrease the sea surface height anomaly change value outward until the outermost contour line only contains one vortex center, and the outermost contour line is the vortex boundary. Wherein, in this embodiment, the change step of gradually increasing or decreasing the sea surface height anomaly change value outward is 0.001 m.

[0100] Step 4: Construct a vortex vertical mode function according to the climatological data and the vortex vertical mode characteristic equation, and extract the first two vertical modes for superposition to obtain the vertical structure function of the vortex horizontal flow velocity;

[0101] In some embodiments, constructing the vortex vertical mode function according to the climatological data and the vortex vertical mode characteristic equation includes:

[0102] Calculate the buoyancy frequency of the target sea area according to the climatological temperature and the climatological salinity;

[0103] Construct a vortex vertical mode function according to the Coriolis parameter, the buoyancy frequency of the target sea area, the vertical coordinate of the grid point, and the vortex vertical mode characteristic equation.

[0104] In this embodiment, the vertical mode characteristic equation of the vortex is expressed as:

[0105] ;

[0106] In the formula, represents taking the derivative, represents taking the derivative with respect to the vertical coordinate of the grid point, represents the eigenvalue, represents the Coriolis parameter, represents the buoyancy frequency of the target sea area, represents the th vertical mode of the vertical mode function of the vortex, represents the vertical coordinate of the grid point, that is, the seawater depth where the climatological temperature and climatological salinity of the target sea area are located.

[0107] In this embodiment, the vertical structure function of the vortex horizontal velocity includes the east - west velocity function of the vortex and the north - south velocity function of the vortex. The east - west velocity function of the vortex and the north - south velocity function of the vortex are respectively expressed as:

[0108] ;

[0109] In the formula, , respectively represent the east - west velocity function of the vortex and the north - south velocity function of the vortex, represents the vertical coordinate of the grid point, , , and represent the vertical mode coefficients of the vortex. Among them, and respectively represent the east - west velocity coefficient and the north - south velocity coefficient of the first vertical mode of the vertical structure function of the vortex horizontal velocity, and respectively represent the east - west velocity coefficient and the north - south velocity coefficient of the second vertical mode of the vertical structure function of the vortex horizontal velocity.

[0110] Step 5: Respectively, with the sea - surface velocity anomaly condition and the deep - layer velocity anomaly condition as constraints, determine the vertical mode coefficients of the vortex according to the vertical structure function of the vortex horizontal velocity.

[0111] In some embodiments, the calculation method of the sea - surface velocity anomaly condition includes:

[0112] Calculate the sea - surface geostrophic current through observing sea - surface height data and the geostrophic relationship.

[0113] Calculate the climatological sea surface geostrophic current through climatological sea surface height data and geostrophic relationships.

[0114] Subtract the climatological sea surface geostrophic current from the sea surface geostrophic current to obtain the sea surface velocity anomaly condition.

[0115] In some embodiments, the calculation method of the deep layer velocity anomaly condition includes:

[0116] Set the horizontal velocity of the vortex at a depth of 2000 meters in the target sea area to 0 to obtain the deep layer velocity anomaly condition.

[0117] In this embodiment, under the sea surface velocity anomaly condition, the east-west velocity function of the vortex and the north-south velocity function of the vortex are respectively expressed as:

[0118] ;

[0119] In the formula, , respectively represent the east-west velocity function of the vortex and the north-south velocity function of the vortex under the sea surface velocity anomaly condition, , respectively represent the first vertical mode and the second vertical mode of the vertical structure function of the vortex horizontal velocity under the sea surface velocity anomaly condition;

[0120] Under the deep layer velocity anomaly condition, the east-west velocity function of the vortex and the north-south velocity function of the vortex are respectively expressed as:

[0121] ;

[0122] In the formula, , respectively represent the east-west velocity function of the vortex and the north-south velocity function of the vortex under the deep layer velocity anomaly condition, , respectively represent the first vertical mode and the second vertical mode of the vertical structure function of the vortex horizontal velocity under the deep layer velocity anomaly condition.

[0123] Step 6: Reconstruct the vertical structure of the vortex horizontal velocity at each grid point according to the vortex vertical mode function and the vortex vertical mode coefficient.

[0124] In this embodiment, the calculation method of the vortex vertical mode coefficient includes:

[0125] Calculate the vortex vertical mode coefficient according to the east-west velocity function of the vortex and the north-south velocity function of the vortex under the sea surface velocity anomaly condition and the east-west velocity function of the vortex and the north-south velocity function of the vortex under the deep layer velocity anomaly condition.

[0126] Step 7: Superimpose the vertical structure of the vortex horizontal flow velocity at each grid point and the vertical structure of the climatological horizontal flow velocity to obtain the vertical structure of the seawater horizontal flow velocity at each grid point, as Figure 2 shown.

[0127] In some embodiments, the vertical structure of the climatological horizontal flow velocity is calculated based on climatological temperature, salinity, sea surface height data, and the thermal wind relationship.

[0128] In some embodiments, the vertical structure of the seawater horizontal flow velocity includes the east - west horizontal flow velocity of seawater and the north - south horizontal flow velocity of seawater; the vertical structure of the vortex horizontal flow velocity includes the east - west horizontal flow velocity of the vortex and the north - south horizontal flow velocity of the vortex; the vertical structure of the climatological horizontal flow velocity includes the east - west horizontal flow velocity of the climatology and the north - south horizontal flow velocity of the climatology.

[0129] In this embodiment, the east - west horizontal flow velocity of seawater and the north - south horizontal flow velocity of seawater are respectively expressed as:

[0130] ;

[0131] wherein, , respectively represent the east - west horizontal flow velocity of seawater and the north - south horizontal flow velocity of seawater, , respectively represent the east - west horizontal flow velocity of the vortex and the north - south horizontal flow velocity of the vortex, and respectively represent the east - west horizontal flow velocity of the climatology and the north - south horizontal flow velocity of the climatology, represents the vertical coordinate of the grid point.

[0132] Step 8: Construct a three - dimensional flow field structure of the vortex according to the vertical structure of the seawater horizontal flow velocity at each grid point.

[0133] As Figure 3 shown, it shows a plan view of the horizontal flow velocity field at a depth of 100 m of the three - dimensional flow field reconstruction structure of the vortex provided by the embodiment of the present invention.

[0134] As Figure 4 shown, it shows a plan view of the horizontal flow velocity field at a depth of 100 m of the vortex of the ARMOR3D data.

[0135] As Figure 5 shown, it shows a plan view of the horizontal flow velocity field at a depth of 300 m of the three - dimensional flow field reconstruction structure of the vortex provided by the embodiment of the present invention.

[0136] As Figure 6As shown, it presents a planar schematic diagram of the horizontal velocity field of the vortex at a depth of 300 meters of the ARMOR3D data.

[0137] Among them, Figure 3 and Figure 5 are respectively the planar schematic diagrams of the horizontal velocity fields at depths of 100 meters and 300 meters of the reconstructed structure of the three-dimensional flow field of the vortex obtained by the method provided by the present invention. Figure 4 and Figure 6 are respectively the planar schematic diagrams of the horizontal velocity fields at depths of 100 meters and 300 meters of the vortex extracted from the existing ARMOR3D data.

[0138] By comparing the horizontal velocity fields at depths of 100 meters and 300 meters, it can be seen that the method provided by the present invention can accurately construct the three-dimensional flow field characteristics of the vortex, solving the problem that the prior art is limited to the sea surface and it is difficult to accurately depict the vertical structure of the vortex.

[0139] Embodiment 3

[0140] Based on the same inventive concept as other embodiments, this embodiment introduces a computer-readable storage medium, on which computer instructions are stored, and when the computer instructions are executed by a processor, the steps of the method in the above Embodiment 1 or 2 are implemented.

[0141] Embodiment 4

[0142] Based on the same inventive concept as other embodiments, this embodiment introduces a computer program product, including computer instructions, and when the computer instructions are executed by a processor, the steps of the method in the above Embodiment 1 or 2 are implemented. In summary of the above embodiments, the present invention efficiently and accurately realizes the reconstruction of the three-dimensional flow field of ocean vortices by integrating climatological data, sea surface observation data, and using the dynamic characteristics of vortices as physical constraint conditions. The present invention not only makes full use of the advantages of high spatio-temporal resolution of satellite remote sensing data, but also combines climatological historical observation data, effectively making up for the deficiencies of a single data source. The present invention constructs the vertical mode function of the vortex and extracts the first two vertical modes for superposition, and determines the vertical mode coefficients of the vortex with the sea surface velocity anomaly condition and the deep velocity anomaly condition as constraints. Only the observed sea surface height data is required to quickly realize vortex feature recognition, construction of the vertical structure of the vortex, and high-precision three-dimensional flow field reconstruction of the vortex, solving the problem that the prior art is limited to the sea surface and it is difficult to accurately depict the vertical structure of the vortex.

[0143] The present invention identifies vortex features by integrating climatological data including temperature, salinity, flow field, and sea surface height with sea surface observation data, and calculating sea surface height anomaly data, and then constructs a vertical mode function of the vortex. It not only makes full use of the complementary advantages of multi-source data, but also ensures the accuracy of the reconstruction of the three-dimensional flow field of the vortex through physical constraints such as the Coriolis parameter and buoyancy frequency. The present invention can accurately depict the vertical structural characteristics of the vortex, providing new details and accuracy for ocean vortex research.

[0144] The present invention determines the vortex center by scanning the sea surface height anomaly data to find the extreme points, including anticyclonic vortices and cyclonic vortices, and defines the vortex boundary based on the sea surface height anomaly change value. It not only improves the efficiency and accuracy of vortex identification, but also helps to deeply understand the dynamic behavior and influence range of the vortex. By accurately defining the vortex boundary, the present invention provides a more reliable basis for the reconstruction of the three-dimensional flow field of the vortex.

[0145] The present invention innovatively combines the sea surface velocity anomaly condition and the deep velocity anomaly condition as physical constraints to determine the vertical mode coefficients of the vortex. The sea surface geostrophic current is calculated by observing the sea surface height data and subtracted from the climatological sea surface geostrophic current to obtain the sea surface velocity anomaly condition; at the same time, the horizontal velocity of the vortex at a depth of 2000 meters in the target sea area is set to 0 as the deep velocity anomaly condition. It not only considers the velocity change characteristics of the vortex at different depth layers, but also ensures the consistency and reliability of the reconstruction of the three-dimensional flow field of the vortex through physical constraints. By integrating the sea surface velocity anomaly condition and the deep velocity anomaly condition, the present invention can more accurately reconstruct the three-dimensional flow field structure of the vortex, providing strong support for fields such as ocean environmental monitoring, resource development, and disaster warning.

[0146] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0147] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to the embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate for implementation in the processFigure 1 one or more processes and / or blocks Figure 1 means for the functions specified in one or more blocks

[0148] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including an instruction means that implements the functions specified in one process Figure 1 one or more processes and / or blocks Figure 1 or more blocks

[0149] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operational steps are performed 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 steps for implementing the functions specified in one process Figure 1 one or more processes and / or blocks Figure 1 or more blocks

[0150] The embodiments of the present invention have been described above in conjunction with the accompanying drawings. However, the present invention is not limited to the above specific embodiments. The above specific embodiments are merely illustrative and not restrictive. Under the inspiration of the present invention, those of ordinary skill in the art can also make many forms without departing from the spirit and scope protected by the claims of the present invention. All of these are within the protection scope of the present invention.

Claims

1. A method for reconstructing three-dimensional flow field of ocean vortex based on physical constraints, characterized in that: include: Collect climate data and sea surface observation data of the target sea area; Calculating sea surface height anomaly data based on the climatological data and sea surface observation data; Identify vortex characteristics based on sea surface height anomaly data and divide the vortex into grids based on the vortex characteristics according to the spatial resolution of the sea surface height observation data; The vortex vertical modal function is constructed according to the climatological data and the vortex vertical modal characteristic equation, and the first two vertical modes are extracted and superimposed to obtain the vortex horizontal flow velocity vertical structure function; Taking the abnormal conditions of sea surface velocity and deep-layer velocity as constraints, the vertical modal coefficients of the vortex are determined according to the vertical structure function of the vortex horizontal velocity. According to the vortex vertical modal function and vortex vertical modal coefficient, the vertical structure of the vortex horizontal flow velocity at each grid point is reconstructed; The vertical structure of climatological horizontal velocity is calculated based on climatological data and thermal wind relationship; The vertical structure of the vortex horizontal velocity at each grid point and the vertical structure of the climatological horizontal velocity are superimposed to obtain the vertical structure of the seawater horizontal velocity at each grid point. The vortex three-dimensional flow field structure is constructed according to the vertical structure of the horizontal flow velocity of seawater at each grid point.

2. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 1 is characterized in that: The climatological data include climatological temperature, climatological salinity, climatological flow field and climatological sea level height; The sea surface observation data includes observed sea surface height data; The method of calculating the sea surface height anomaly data based on the climatological data and the sea surface observation data includes: The sea surface height anomaly data is obtained by subtracting the climatological sea surface height data from the observed sea surface height data.

3. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 1 is characterized in that: The vortex characteristics include a vortex center and a vortex boundary; The vortex center includes an anticyclonic vortex potential center and a cyclonic vortex potential center; Wherein, identifying vortex characteristics according to abnormal sea level height data includes: If the spatial resolution of the sea surface height anomaly data is 1 / 4º×1 / 4º, a 5×5 grid point window is scanned to find the extreme points, where the maximum point is the potential center of the anticyclonic vortex and the minimum point is the potential center of the cyclonic vortex; Taking the sea level height anomaly data of the vortex center as the starting value, the sea level height anomaly change value is gradually increased or decreased outward until the outermost contour line only contains one vortex center, and the outermost contour line is the vortex boundary.

4. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 2 is characterized in that: The vortex vertical mode function is constructed based on the climatological data and the vortex vertical mode characteristic equation, including: Calculate the buoyancy frequency of the target sea area according to the climatological temperature and climatological salinity; The vortex vertical modal function is constructed based on the Coriolis parameters, the buoyancy frequency of the target sea area, the vertical coordinates of the grid points and the vortex vertical modal characteristic equation.

5. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 4 is characterized in that: The vortex vertical modal characteristic equation is expressed as: ； In the formula, represents the derivation, Represents the vertical coordinates of the grid points Seeking guidance, represents the eigenvalue, represents the Coriolis parameter, Indicates the buoyancy frequency of the target sea area, represents the vertical mode function of the vortex Vertical Mode Seeking guidance, Represents the vertical coordinates of the grid points The vertical coordinate is the sea water depth where the climatological temperature and climatological salinity of the target sea area are located.

6. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 4, characterized in that: The vortex horizontal velocity vertical structure function includes a vortex east-west velocity function and a vortex north-south velocity function; The east-west flow velocity function of the vortex and the north-south flow velocity function of the vortex are respectively expressed as: ； In the formula, , They represent the vortex east-west velocity function and the vortex north-south velocity function, respectively. represents the vertical coordinate of the grid point, , , and represents the vortex vertical modal coefficient, where and They represent the first vertical mode of the vertical structure function of the vortex horizontal flow velocity. The east-west velocity coefficient and the north-south velocity coefficient are and They represent the second vertical mode of the vertical structure function of the vortex horizontal flow velocity. The east-west flow velocity coefficient and the north-south flow velocity coefficient.

7. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 1, characterized in that: Said The calculation methods of abnormal sea surface velocity conditions include: The sea surface geostrophic current is obtained by calculating the sea surface height data and geostrophic relationship; The climatological sea surface geostrophic current is obtained by calculation through climatological sea surface height data and geostrophic relationship. The abnormal conditions of sea surface velocity are obtained by subtracting the climatological sea surface geostrophic current from the sea surface geostrophic current. The calculation method of the abnormal deep velocity condition includes: The horizontal flow velocity of the vortex at a depth of 2000 meters in the target sea area is set to 0, and the abnormal deep flow velocity condition is obtained.

8. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 7, characterized in that: Under the condition of abnormal sea surface current velocity, the east-west velocity function of the vortex and the north-south velocity function of the vortex are respectively expressed as: ； In the formula, , They represent the vortex east-west velocity function and the vortex north-south velocity function under abnormal sea surface velocity conditions, respectively. , They represent the first and second vertical modes of the vertical structure function of the vortex horizontal velocity under abnormal sea surface velocity conditions, respectively. , , and represents the vortex vertical modal coefficient, where and Respectively The east-west velocity coefficient and the north-south velocity coefficient are and Respectively The east-west velocity coefficient and the north-south velocity coefficient; Under the condition of abnormal deep flow velocity, the east-west flow velocity function of the vortex and the north-south flow velocity function of the vortex are respectively expressed as: ； In the formula, , They represent the vortex east-west velocity function and the vortex north-south velocity function under the condition of abnormal deep velocity. , They respectively represent the first vertical mode and the second vertical mode of the vertical structure function of the vortex horizontal velocity under the condition of abnormal deep flow velocity.

9. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 8, characterized in that: The method for calculating the vortex vertical modal coefficient comprises: The vortex vertical modal coefficients are calculated based on the vortex east-west flow velocity function and the vortex north-south flow velocity function under abnormal sea surface flow velocity conditions and the vortex east-west flow velocity function and the vortex north-south flow velocity function under abnormal deep flow velocity conditions.

10. The method for reconstructing the three-dimensional flow field of ocean vortex based on physical constraints according to claim 1, characterized in that: The vertical structure of the seawater horizontal flow velocity includes the seawater east-west horizontal flow velocity and the seawater north-south horizontal flow velocity; The vertical structure of the vortex horizontal flow velocity includes the vortex east-west horizontal flow velocity and the vortex north-south horizontal flow velocity; The climatic horizontal flow velocity vertical structure includes the climatic east-west horizontal flow velocity and the climatic north-south horizontal flow velocity; The east-west horizontal flow velocity of the seawater and the north-south horizontal flow velocity of the seawater are respectively expressed as: ； In the formula, , They represent the east-west horizontal flow velocity and the north-south horizontal flow velocity of seawater, respectively. , They represent the east-west horizontal flow velocity of the vortex and the north-south horizontal flow velocity of the vortex, and They represent the climatic east-west horizontal flow velocity and the climatic north-south horizontal flow velocity, respectively. Represents the vertical coordinate of the grid point.

Images