Sea surface wave digital twin modeling method, device, electronic device and storage medium

Abstract

The present invention relates to the technical field of sea surface modeling and provides a sea surface surge digital twin modeling method, apparatus, electronic device, and storage medium. The method comprises: determining surge parameters based on a surge factor and a hyperbolic tangent function, wherein the independent variables of the hyperbolic tangent function are determined by the real-time wave value and peak wave value of a target sea spectrum; determining a surge factor function based on the surge parameters and the wave number direction angle of the target sea spectrum; determining a target direction spread function based on a basic directional spread function, the surge factor function, and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target direction spread function based on the wave number direction angle is 1; and constructing a sea surface surge model using the target sea spectrum as a base spectrum and in combination with the target directional spread function. The present invention can construct a sea surface surge model with high fidelity while requiring relatively low computational effort.

Description

Technical Field

[0001] The present invention relates to the technical field of sea surface modeling, and in particular to a method, device, electronic device and storage medium for digital twin modeling of sea surface surges. Background Art

[0002] Sea surface modeling can be used for radar signal simulation and sea clutter analysis. However, in challenging scenarios, such as high sea conditions, collecting radar sea clutter data is expensive, and comprehensive measurement and analysis of radar electromagnetic echoes at varying wind speeds, directions, or grazing angles is impossible. To reduce costs, a sea surface surge model can be used as the basis for radar electromagnetic echo simulation for radar signal simulation and sea clutter analysis.

[0003] Currently, sea surge modeling methods are primarily categorized into physical modeling and mathematical modeling. Physical modeling yields models with high fidelity and conformity to physical laws, but due to the enormous computational complexity, it is difficult to apply to large-scale sea surge simulations. Mathematical modeling offers high computational speed, but the resulting sea surge models lack fidelity.

[0004] In summary, existing sea surface surge modeling methods cannot achieve a high-fidelity sea surface surge model with low computational complexity. Summary of the Invention

[0005] The present invention provides a sea surface wave digital twin modeling method, device, electronic device and storage medium to solve the problem in the prior art that it is impossible to achieve a high-fidelity sea surface wave model with low computational complexity.

[0006] The present invention provides a sea surface wave digital twin modeling method, comprising the following steps:

[0007] Determining surge parameters based on the surge factor and the hyperbolic tangent function, wherein the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum;

[0008] determining a surge factor function based on the surge parameter and the wave number direction angle of the target ocean spectrum;

[0009] Determining a target directional spread function based on the basic directional spread function, the surge factor function, and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1;

[0010] The target sea spectrum is used as the base spectrum and combined with the target direction spread function to construct a sea surface surge model.

[0011] According to a sea surface surge digital twin modeling method provided by the present invention, surge parameters are determined based on surge factors and hyperbolic tangent functions, including: determining surge parameters according to the following formula :

[0012] ;

[0013] Where tanh(·) represents the hyperbolic tangent function, e Indicates the surge factor, the value range is [0,1], k Indicates the real-time wave value of the target sea spectrum, k p Indicates the peak wave value of the target ocean spectrum.

[0014] According to a sea surface surge digital twin modeling method provided by the present invention, a surge factor function is determined based on the surge parameters and the wave number direction angle of the target sea spectrum, including: determining the surge factor function according to the following formula :

[0015] ;

[0016] in, k Indicates the real-time wave value of the target sea spectrum, θ represents the wave number direction angle of the target sea spectrum, Indicates surge parameters.

[0017] According to a sea surface surge digital twin modeling method provided by the present invention, a target directional expansion function is determined based on a basic directional expansion function, a surge factor function and a first regularization function, including: determining the target directional expansion function according to the following formula :

[0018] ;

[0019] ;

[0020] in, represents the basic direction expansion function, represents the surge factor function, represents the first regularization function, k Indicates the real-time wave value of the target sea spectrum, θ Indicates the wave number direction angle of the target ocean spectrum.

[0021] According to a sea surge digital twin modeling method provided by the present invention, the basic direction spread function for:

[0022] ;

[0023] ;

[0024] ;

[0025] ;

[0026] in, k p represents the peak wave value of the target sea spectrum, U represents the sea surface wind speed, represents the acceleration due to gravity, represents the second regularization coefficient, represents the gamma function, s and s p is an intermediate variable.

[0027] According to a sea surface wave digital twin modeling method provided by the present invention, a sea surface wave model is constructed using the target sea spectrum as a base spectrum and in combination with the target direction spread function, including:

[0028] converting the target ocean spectrum into a two-dimensional ocean spectrum by using the target direction expansion function;

[0029] Constructing a sea surface height model in the frequency domain based on the two-dimensional sea spectrum;

[0030] Converting the sea surface height model in the frequency domain into the sea surface height model in the spatial domain;

[0031] The sea surface height model in the spatial domain is determined to be the sea surface surge model.

[0032] According to a sea surface wave digital twin modeling method provided by the present invention, the target sea spectrum adopts the Elfouhaily omnidirectional sea spectrum.

[0033] The present invention also provides a sea surface wave digital twin modeling device, comprising the following modules:

[0034] A surge parameter determination module is used to determine the surge parameters based on the surge factor and the hyperbolic tangent function, wherein the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum;

[0035] a surge factor function determination module, configured to determine a surge factor function based on the surge parameters and the wave number direction angle of the target sea spectrum;

[0036] a target directional spread function determination module, configured to determine a target directional spread function based on a basic directional spread function, a surge factor function, and a first regularization function, wherein the first regularization function is configured to ensure that an integral of the target directional spread function based on a wave number direction angle is 1;

[0037] The surge model construction module is used to construct a sea surface surge model based on the target sea spectrum as the base spectrum and in combination with the target direction expansion function.

[0038] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and running on the processor. When the processor executes the program, it implements the digital twin modeling method of sea surface surge as described in any one of the above.

[0039] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the digital twin modeling method of sea surge as described in any one of the above is implemented.

[0040] The digital twin modeling method, device, electronic device, and storage medium for sea surge provided by the present invention determine surge parameters by setting a surge factor and combining it with the hyperbolic tangent function to satisfy the asymptotically increasing property of surges. This allows for the subsequent determination of a surge factor function that superimposes the surge factor based on the surge parameters. Furthermore, the basic directional expansion function, the surge factor function, and the first regularization function are combined to obtain a target directional expansion function that reflects the properties of the surge. This allows for the construction of a highly realistic sea surge model using the target sea spectrum as the base spectrum and the target directional expansion function of the present invention. Furthermore, compared to current sea surge modeling methods based on fluid mechanics equations (physical modeling), the sea surge modeling method based on the target sea spectrum combined with the target directional expansion function proposed by the present invention has lower computational complexity and can effectively address the difficulty of modeling surges in large sea areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] In order to more clearly illustrate the technical solutions in the present invention or the prior art, a brief introduction is given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0042] Figure 1 It is a flow chart of the sea surface surge digital twin modeling method provided by the present invention.

[0043] Figure 2 It is a schematic diagram of the variance spectrum and curvature spectrum of the Elfouhaily omnidirectional sea spectrum in the sea surface wave digital twin modeling method provided by the present invention.

[0044] Figure 3 This is a graph showing the change of the target direction expansion function with the surge factor in the sea surface surge digital twin modeling method provided by the present invention.

[0045] Figure 4It is a sea surface morphology diagram obtained by modeling when the surge factor is 0 in the sea surface surge digital twin modeling method provided by the present invention.

[0046] Figure 5 It is a sea surface morphology diagram obtained by modeling when the surge factor is 0.5 in the sea surface surge digital twin modeling method provided by the present invention.

[0047] Figure 6 It is a sea surface morphology diagram obtained by modeling when the surge factor is 0.8 in the sea surface surge digital twin modeling method provided by the present invention.

[0048] Figure 7 It is a sea surface morphology diagram obtained by modeling when the surge factor is 1 in the sea surface surge digital twin modeling method provided by the present invention.

[0049] Figure 8 It is a sea surface morphology diagram obtained by modeling in the sea surface surge digital twin modeling method provided by the present invention when the surge factor is 1 and the wind speed is 2 meters per second.

[0050] Figure 9 It is a sea surface morphology diagram obtained by modeling in the sea surface surge digital twin modeling method provided by the present invention when the surge factor is 1 and the wind speed is 4 meters per second.

[0051] Figure 10 It is a sea surface morphology diagram obtained by modeling in the sea surface surge digital twin modeling method provided by the present invention when the surge factor is 1 and the wind speed is 6 meters per second.

[0052] Figure 11 This is a comparison diagram of the echoes of the random sea surface and the sea surface waves under sea condition 1 in the sea surface wave digital twin modeling method provided by the present invention.

[0053] Figure 12 This is a comparison diagram of the echoes of the random sea surface and the sea surface surge in the case of sea state 3 in the sea surface surge digital twin modeling method provided by the present invention.

[0054] Figure 13 This is a comparison diagram of the echoes of random sea surface and sea surface waves in sea state 5 in the sea surface wave digital twin modeling method provided by the present invention.

[0055] Figure 14 It is a distribution simulation diagram of sea condition 1 in the sea surface surge digital twin modeling method provided by the present invention.

[0056] Figure 15 It is a distribution simulation diagram of sea state 3 in the sea surface wave digital twin modeling method provided by the present invention.

[0057] Figure 16It is a distribution simulation diagram of sea state 5 in the sea surface wave digital twin modeling method provided by the present invention.

[0058] Figure 17 This is a distribution simulation diagram of measured sea clutter in the sea surface wave digital twin modeling method provided by the present invention.

[0059] Figure 18 It is a structural schematic diagram of the sea surface surge digital twin modeling device provided by the present invention.

[0060] Figure 19 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0061] To make the objectives, technical solutions, and advantages of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0062] The digital twin modeling method of sea surge in the embodiment of the present invention is as follows: Figure 1 As shown, it includes steps S110 to S140.

[0063] Step S110: Determine surge parameters based on the surge factor and the hyperbolic tangent function, wherein the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum, wherein wave number = 2π / wavelength.

[0064] The target ocean spectrum is a pre-selected ocean spectrum used to generate the sea surface model. However, the sea surface model generated using the target ocean spectrum is typically a linear ocean surface model and cannot simulate the effects of swells. Therefore, it is necessary to add a swell factor and combine it with the hyperbolic tangent function to ensure that the swell formation (elongation) increases asymptotically, ultimately allowing the construction of a sea surface swell model.

[0065] Optionally, in this embodiment, the target ocean spectrum preferably adopts the Elfouhaily omnidirectional ocean spectrum, such as Figure 2The following diagrams show the variance spectrum (left) and curvature spectrum (right) of the Elfouhaily omnidirectional ocean spectrum at a wind speed of 10 m / s. The first peak is clearly visible in the variance spectrum. Since the second peak is not apparent at wave number k = 370 rad / m, multiplying the wave number cubed to convert the variance spectrum into a curvature spectrum clearly reveals the second peak at k = 370 rad / m. This corresponds to the double peak observed in real sea areas. Therefore, the Elfouhaily omnidirectional ocean spectrum was selected as the target ocean spectrum, resulting in a more realistic sea surge model.

[0066] Step S120: Determine a surge factor function based on the surge parameter and the wave number direction angle of the target ocean spectrum. For example, the surge factor function may be a cosine power function of an independent variable determined by the wave number direction angle of the target ocean spectrum. Specifically, the surge factor function is obtained by raising the cosine power function to the power of the surge parameter.

[0067] In this step, since the surge factor function includes surge parameters, and the surge parameters are obtained based on the surge factor combined with the hyperbolic tangent function, the surge factor function also satisfies the property of asymptotic increase in surge formation (elongation). Based on the surge factor function, the surge effects caused by ultra-far wind zones and gravity waves can be simulated. The surge factor function causes the waves to gradually elongate as the wavelength increases, asymptotically approaching the maximum elongation value.

[0068] Step S130: Determine a target directional spread function based on the base directional spread function, the surge factor function, and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1, so as to ensure that the target directional spread function does not affect the total amount of ocean spectrum energy.

[0069] In this step, the basic directional spread function is combined with the surge factor function and the first regularization function to obtain a target directional spread function, which also satisfies the property of asymptotic increase in surge shaping (elongation).

[0070] Step S140: Using the target ocean spectrum as the base spectrum and combining it with the target directional spread function, a sea surface swell model is constructed. For example, using the Elfouhaily omnidirectional ocean spectrum as the target ocean spectrum and combining it with a conventional basic directional spread function (such as the unilateral cosine directional spread function proposed by Longuet-Higgins) can only generate a linear sea surface. However, using the target ocean spectrum as the base spectrum and combining it with the target directional spread function of this embodiment can construct a more realistic sea surface swell model.

[0071] In the digital twin modeling method for sea surge waves of this embodiment, the surge parameters are determined by setting the surge factor and combining it with the hyperbolic tangent function to satisfy the asymptotically increasing property of the surge. This allows the surge factor function, which is superimposed on the surge factor, to be determined based on the surge parameters. Furthermore, the target directional spread function, the surge factor function, and the first regularization function are combined to obtain a target directional spread function that reflects the properties of the surge waves. This allows the target ocean spectrum to be used as the base spectrum, combined with the target directional spread function of this embodiment, to construct a highly realistic sea surge model. Furthermore, compared to current sea surge modeling methods based on fluid mechanics equations (physical modeling), the sea surge modeling method based on the target ocean spectrum combined with the target directional spread function proposed in this embodiment has lower computational complexity and can effectively address the difficulty of modeling surge waves over large sea areas.

[0072] In some embodiments, determining the surge parameter based on the surge factor and the hyperbolic tangent function includes: determining the surge parameter according to the following formula :

[0073] (1).

[0074] Where tanh(·) represents the hyperbolic tangent function, e Indicates the surge factor, the value range is [0,1], k Indicates the real-time wave value of the target sea spectrum, k p Indicates the peak wave value of the target ocean spectrum.

[0075] In this embodiment, the surge factor The square of produces the perceptual linearization of the surge parameters, which brings about the surge effect. The hyperbolic tangent function tanh(·) is introduced to construct the surge parameters. The hyperbolic tangent function tanh(·) has the following characteristics:

[0076] (1) Nonlinear saturation characteristics: simulates the gradual attenuation process of surge energy when it deviates from the main wave direction.

[0077] (2) Parameter adjustability: by introducing surge factor ∈[0,1] regulates the distribution of surge waves.

[0078] (3) Directional bias capability: Achieve asymmetric stretching of energy distribution along the main wave axis.

[0079] when = 0, it degenerates into the classical model, i.e., the linear sea surface model, corresponding to the fully developed wind and wave state; →1, the distribution of the directional spread function is significantly elongated, simulating the mature swells that propagate over long distances, and the intermediate value represents the mixed state of wind waves and swells. Figure 3As shown in the figure, the closer the surge factor is to 1, the more elongated and asymmetric the corresponding target direction spread function is. The sea surface surge model constructed in combination with the bottom spectrum can better demonstrate the ductility and asymmetry of the surge, that is, the constructed sea surface surge model is more realistic. On the other hand, the closer the surge factor is to 0, the more symmetric the corresponding target direction spread function is, indicating that the surge effect is worse. It should be noted that: generally, when e>0.5, the surge state can be seen. Not all surges are necessarily in the most complete form. As long as e is not equal to 0, the final constructed sea surface surge model is.

[0080] In this embodiment, since the hyperbolic tangent function tanh(·) has the above three characteristics, the target direction spread function also has the above three characteristics. Therefore, the sea surface surge model constructed by combining the target direction spread function with the target sea spectrum can simulate the sea surface surge more realistically.

[0081] In some embodiments, determining a surge factor function based on the surge parameter and the wave number direction angle of the target sea spectrum includes: determining the surge factor function according to the following formula :

[0082] (2).

[0083] in, k Indicates the real-time wave value of the target sea spectrum, θ represents the wave number direction angle of the target sea spectrum, Indicates surge parameters.

[0084] In this embodiment, the directional spread function of the Longuet-Higgins unilateral cosine form can be used as the basis, and the above-mentioned surge parameter is used as the power of the cosine function. Since the hyperbolic tangent function tanh(·) has the above-mentioned three characteristics, the corresponding surge factor function also has the above-mentioned three characteristics, so that the finally constructed sea surface surge model can simulate the sea surface surge more realistically.

[0085] In some embodiments, determining a target directional spread function based on the base directional spread function, the surge factor function, and the first regularization function includes: determining the target directional spread function according to the following formula: :

[0086] (3);

[0087] (4).

[0088] in, represents the basic direction expansion function, represents the surge factor function, represents the first regularization function, kIndicates the real-time wave value of the target sea spectrum, θ Indicates the wave number direction angle of the target ocean spectrum.

[0089] Furthermore, the basic directional spread function can adopt the directional spread function of the Longuet-Higgins single-sided cosine form. The directional spread function of the single-sided cosine form has good scalability, fast integral calculation and power calculation, and is easy to calculate. Specifically, the basic directional spread function for:

[0090] (5);

[0091] (6);

[0092] (7);

[0093] (8).

[0094] in, k p represents the peak wave value of the target sea spectrum, U represents the sea surface wind speed, represents the acceleration due to gravity, represents the second regularization coefficient, represents the gamma function, s and s p is an intermediate variable.

[0095] The Longuet-Higgins model uses a cosine power form. Although it can describe the concentrated characteristics of the main wave energy, its symmetrical bell-shaped distribution makes it difficult to describe the directional extension characteristics unique to swells. By analyzing the measured swell data, it is found that when the wind zone is fully developed, the wave direction distribution presents a significant asymmetric elongated shape, and the energy attenuation presents a nonlinear asymptotic characteristic. Therefore, the traditional Longuet-Higgins model cannot show the characteristics of swells. In this embodiment, by using formula (3), the surge factor function and the first regularization function are superimposed on the basic directional extension function of Longuet-Higgins to obtain a target directional extension function that can describe the directional extension characteristics unique to swells, and ensure that the target directional extension function does not affect the total size of the ocean spectrum energy.

[0096] In some embodiments, the target sea spectrum is used as a base spectrum, and a sea surface surge model is constructed in combination with the target directional spread function, including:

[0097] The target ocean spectrum is converted into a two-dimensional ocean spectrum by the target direction expansion function. For example, the Elfouhaily omnidirectional ocean spectrum is used as the target ocean spectrum. , the target spectrum is converted into (i.e. Elfouhaily omnidirectional ocean spectrum) is converted into a two-dimensional ocean spectrum.

[0098] (9).

[0099] in, is a function representing the two-dimensional ocean spectrum, represents the wave number vector, , , , n and m Sea surface x peacekeeping y The number of sampling points of the dimension, is the wave number direction angle.

[0100] The sea surface height model in the frequency domain is constructed based on the two-dimensional sea spectrum. .

[0101] (10).

[0102] in, t Indicates the moment, represents a complex Gaussian random sequence with a mean of 0 and a variance of 1. and are the length and width of the constructed sea surface, represents the complex conjugate, and Both represent inverse Fourier transform.

[0103] Convert the sea surface height model in the frequency domain to the sea surface height model in the spatial domain .

[0104] (11).

[0105] Where i represents the imaginary unit, represents the Fourier transform, v Represents the (x,y) coordinates in the sea level coordinate system.

[0106] The sea surface height model in the spatial domain is determined to be the sea surface surge model, that is, The final sea surge model was constructed.

[0107] In this embodiment, the Elfouhaily omnidirectional sea spectrum is used as the target sea spectrum, and a sea surface surge model that can truly reflect the surge characteristics is constructed through the above steps.

[0108] like Figure 4 、 Figure 5 、 Figure 6 and Figure 7 As shown in the figure, based on the target direction expansion function, under the conditions of sea surface wind speed of 4m / s, wind direction of 0°, sea surface length and width of 50 meters, and the number of sampling points in the x and y dimensions of 100, different surge factors e The sea surface simulation effect. It can be seen that the surge factor e When it is 0 (i.e. Figure 4 ), the sea surface presents a random undulating shape, and with the surge factor e As the surge factor increases, the sea surface gradually shows the shape of swells. e When it is 1 (i.e. Figure 7 ), the waves are strongly elongated, showing the maximum swell effect.

[0109] like Figure 8 、 Figure 9 and Figure 10 As shown, they are respectively at the same surge factor e ( e =1), the sea surface shapes at wind speeds of 2m / s, 4m / s, and 6m / s. It can be seen that the swell effect can be well reproduced even at low wind speeds. However, at higher wind speeds, the swell effect is driven by the roll effect and is less obvious. However, the realistic sea surface shape is still presented.

[0110] The normalized sea clutter radar cross section (RCS) is calculated based on specific radar parameters (such as operating frequency, polarization mode, and pulse repetition rate) and environmental conditions (such as wind speed, wind direction, and sea state). These normalized RCS values ​​are then substituted into a far-field scattering model (TSC model), which takes into account factors such as the distance between the radar and the sea surface, the radar altitude, and the angle of incidence to generate sea clutter echo signals under different sea conditions.

[0111] like Figure 11 、 Figure 12 and Figure 13 , respectively, compare the radar sea clutter echo signal from a random sea surface (left) with the radar sea clutter echo signal from a sea swell model (right) in sea states 1, 3, and 5. The specific results for sea states 1, 3, and 5 are shown in Table 1. It is clear that, while the sea state still influences the echo, the radar echo from the sea surface in the presence of swells has a more striped and compressed shape, better reflecting the swell effect. This means that the radar echo from sea swells better reflects the actual sea surface radar clutter, eliminating the need for electromagnetic simulation, resulting in faster simulation speed and lower simulation costs.

[0112] Table 1 Description of sea state 1, sea state 3 and sea state 5

[0113]

[0114] like Figure 14 、 Figure 15 and Figure 16 Figures 2 and 3 show the Cumulative Distribution Function (CDF) statistical plots of the sea swell models generated by the above embodiments for sea conditions 1, 3, and 5. Each figure demonstrates the CDF's fitting of three distributions: the K distribution, the Weibull distribution, and the lognormal distribution. It can be seen that under different sea conditions, the echoes generated by the sea swell models constructed by the above embodiments are relatively consistent with the typical sea clutter amplitude distributions previously described.

[0115] At the same time, if Figure 17 As shown, the measured sea clutter ( Figure 17 ) were compared, Figure 17 The fitting effects of the three distributions above are used for the measured sea clutter collected in 2022. It can be seen that the process of generating the sea surge model and performing electromagnetic scattering simulation using the above embodiments can well simulate the real sea clutter situation.

[0116] As shown in Table 2, the KS test primarily reflects the maximum difference between the empirical distribution and the theoretical distribution; smaller values ​​indicate better fit. The root mean square error (RMSE) measures the accuracy of the probability density function (PDF) fit; smaller values ​​indicate smaller fit errors. It can be seen that, according to the KS test, the sea surge model constructed in this invention performs better than the measured sea clutter in fitting the CDF distribution and is consistent with the classical sea clutter distribution. However, with respect to the RMSE, it performs worse than the measured sea clutter. The conflict between the RMSE and KS test results may stem from the contradiction between local and global fit (a low RMSE reflects a better fit in the midsection of the PDF, while a high KS value indicates differences in the CDF tail) as well as data heterogeneity.

[0117] Table 2 Distribution fit test

[0118]

[0119] The digital twin modeling device for sea surface waves provided by the present invention is described below. The digital twin modeling device for sea surface waves described below and the digital twin modeling method for sea surface waves described above can be referenced to each other.

[0120] The digital twin modeling device for sea surge in an embodiment of the present invention is as follows: Figure 18 Shown, including:

[0121] The surge parameter determination module 1810 is used to determine the surge parameters based on the surge factor and the hyperbolic tangent function, where the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum.

[0122] The surge factor function determination module 1820 is configured to determine a surge factor function based on the surge parameters and the wave number direction angle of the target sea spectrum.

[0123] The target directional spread function determination module 1830 is used to determine the target directional spread function based on the basic directional spread function, the surge factor function and the first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1.

[0124] The surge model construction module 1840 is configured to construct a sea surface surge model using the target sea spectrum as a base spectrum and in combination with the target directional spread function.

[0125] The digital twin modeling device for sea surge waves in this embodiment determines surge parameters by setting a surge factor and combining it with the hyperbolic tangent function to satisfy the asymptotically increasing nature of surge waves. This allows for the subsequent determination of a surge factor function based on the surge parameters, which superimposes the surge factor. This device then combines the basic directional spread function, the surge factor function, and a first regularization function to obtain a target directional spread function that reflects the properties of the surge waves. This allows the target ocean spectrum to be used as the base spectrum, combined with the target directional spread function of this embodiment, to construct a highly realistic sea surge model. Furthermore, compared to current sea surge modeling methods based on fluid mechanics equations (physical modeling), the sea surge modeling method based on the target ocean spectrum combined with the target directional spread function proposed in this embodiment has lower computational complexity and can effectively address the difficulty of modeling surge waves over large sea areas.

[0126] In some embodiments, the surge parameter determination module 1810 is specifically configured to determine the surge parameter according to the following formula: :

[0127] ;

[0128] Where tanh(·) represents the hyperbolic tangent function, e Indicates the surge factor, the value range is [0,1], k Indicates the real-time wave value of the target sea spectrum, k p Indicates the peak wave value of the target ocean spectrum.

[0129] In some embodiments, the surge factor function determination module 1820 is specifically configured to determine the surge factor function according to the following formula: :

[0130] ;

[0131] in, k Indicates the real-time wave value of the target sea spectrum, θ represents the wave number direction angle of the target sea spectrum, Indicates surge parameters.

[0132] In some embodiments, the target direction spread function determination module 1830 is specifically configured to determine the target direction spread function according to the following formula: :

[0133] ;

[0134] ;

[0135] in, represents the basic direction expansion function, represents the surge factor function, represents the first regularization function, k Indicates the real-time wave value of the target sea spectrum, θ Indicates the wave number direction angle of the target ocean spectrum.

[0136] In some embodiments, the base direction spread function for:

[0137] ;

[0138] ;

[0139] ;

[0140] ;

[0141] in, k p represents the peak wave value of the target sea spectrum, U represents the sea surface wind speed, represents the acceleration due to gravity, represents the second regularization coefficient, represents the gamma function, s and s p is an intermediate variable.

[0142] In some embodiments, the surge model building module 1840 is specifically configured to:

[0143] converting the target ocean spectrum into a two-dimensional ocean spectrum by using the target direction expansion function;

[0144] Constructing a sea surface height model in the frequency domain based on the two-dimensional sea spectrum;

[0145] Converting the sea surface height model in the frequency domain into the sea surface height model in the spatial domain;

[0146] The sea surface height model in the spatial domain is determined to be the sea surface surge model.

[0147] In some embodiments, the target sea spectrum adopts the Elfouhaily omnidirectional sea spectrum.

[0148] Figure 19 An example of a physical structure diagram of an electronic device is shown below. Figure 19 As shown, the electronic device may include: a processor 1910, a communications interface 1920, a memory 1930, and a communication bus 1940, wherein the processor 1910, the communications interface 1920, and the memory 1930 communicate with each other via the communication bus 1940. The processor 1910 may call the logic instructions in the memory 1930 to execute the sea surge digital twin modeling method, which includes:

[0149] The surge parameters are determined based on the surge factor and the hyperbolic tangent function, wherein the independent variables of the hyperbolic tangent function are determined by the real-time wave value and the peak wave value of the target sea spectrum.

[0150] A surge factor function is determined based on the surge parameters and the wave number direction angle of the target sea spectrum.

[0151] A target directional spread function is determined based on the basic directional spread function, the surge factor function and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1.

[0152] The target sea spectrum is used as the base spectrum and combined with the target direction spread function to construct a sea surface surge model.

[0153] Furthermore, the logic instructions in the aforementioned memory 1930 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product, stored in a storage medium, includes instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0154] On the other hand, the present invention further provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the sea surge digital twin modeling method provided by the above methods, which includes:

[0155] The surge parameters are determined based on the surge factor and the hyperbolic tangent function, wherein the independent variables of the hyperbolic tangent function are determined by the real-time wave value and the peak wave value of the target sea spectrum.

[0156] A surge factor function is determined based on the surge parameters and the wave number direction angle of the target sea spectrum.

[0157] A target directional spread function is determined based on the basic directional spread function, the surge factor function and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1.

[0158] The target sea spectrum is used as the base spectrum and combined with the target direction spread function to construct a sea surface surge model.

[0159] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the sea surge digital twin modeling method provided by the above methods, the method comprising:

[0160] The surge parameters are determined based on the surge factor and the hyperbolic tangent function, wherein the independent variables of the hyperbolic tangent function are determined by the real-time wave value and the peak wave value of the target sea spectrum.

[0161] A surge factor function is determined based on the surge parameters and the wave number direction angle of the target sea spectrum.

[0162] A target directional spread function is determined based on the basic directional spread function, the surge factor function and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1.

[0163] The target sea spectrum is used as the base spectrum and combined with the target direction spread function to construct a sea surface surge model.

[0164] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of the present embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0165] Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented using software plus a necessary general-purpose hardware platform, or of course, hardware. Based on this understanding, the essence of the above technical solution, or the portion that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, or an optical disk, and includes a number of instructions for causing a computer device (such as a personal computer, server, or network device) to execute the methods described in each embodiment or certain portions of the embodiments.

[0166] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A sea surface wave digital twin modeling method, characterized in that: include: Determining surge parameters based on the surge factor and the hyperbolic tangent function, wherein the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum; determining a surge factor function based on the surge parameter and the wave number direction angle of the target sea spectrum; Determining a target directional spread function based on the basic directional spread function, the surge factor function, and a first regularization function, wherein the first regularization function is used to ensure that the integral of the target directional spread function based on the wave number direction angle is 1; Taking the target sea spectrum as the base spectrum and combining it with the target direction spread function, a sea surface surge model is constructed; Wherein, based on the basic directional spread function, the surge factor function and the first regularization function, the target directional spread function is determined, including: determining the target directional spread function according to the following formula : ； ； in, represents the basic direction expansion function, represents the surge factor function, represents the first regularization function, k Indicates the real-time wave value of the target sea spectrum, θ Indicates the wave number direction angle of the target sea spectrum; Among them, the basic direction expansion function for: ； ； ； ； in, k p represents the peak wave value of the target sea spectrum, U represents the sea surface wind speed, represents the acceleration due to gravity, represents the second regularization coefficient, represents the gamma function, s and s p is an intermediate variable.

2. The sea surface wave digital twin modeling method according to claim 1, characterized in that: Based on the surge factor and the hyperbolic tangent function, the surge parameters are determined, including: determining the surge parameters according to the following formula : ； Where tanh(·) represents the hyperbolic tangent function, e Indicates the surge factor, the value range is [0,1], k Indicates the real-time wave value of the target sea spectrum, k p Indicates the peak wave value of the target ocean spectrum.

3. The sea surge digital twin modeling method according to claim 1, characterized in that: Determining a surge factor function based on the surge parameter and the wave number direction angle of the target sea spectrum includes: determining the surge factor function according to the following formula : ； in, k Indicates the real-time wave value of the target sea spectrum, θ represents the wave number direction angle of the target sea spectrum, Indicates surge parameters.

4. The sea surge digital twin modeling method according to claim 1, characterized in that: The target sea spectrum is used as the base spectrum, and the target direction spread function is combined to construct a sea surface surge model, including: converting the target ocean spectrum into a two-dimensional ocean spectrum by using the target direction expansion function; Constructing a sea surface height model in the frequency domain based on the two-dimensional sea spectrum; Converting the sea surface height model in the frequency domain into the sea surface height model in the spatial domain; The sea surface height model in the spatial domain is determined to be the sea surface surge model.

5. The sea surge digital twin modeling method according to any one of claims 1 to 4, characterized in that: The target sea spectrum adopts the Elfouhaily omnidirectional sea spectrum.

6. A digital twin modeling device for sea surge, characterized in that: include: A surge parameter determination module is used to determine the surge parameters based on the surge factor and the hyperbolic tangent function, wherein the independent variable of the hyperbolic tangent function is determined by the real-time wave value and the peak wave value of the target sea spectrum; a surge factor function determination module, configured to determine a surge factor function based on the surge parameters and the wave number direction angle of the target sea spectrum; a target directional spread function determination module, configured to determine a target directional spread function based on a basic directional spread function, a surge factor function, and a first regularization function, wherein the first regularization function is configured to ensure that an integral of the target directional spread function based on a wave number direction angle is 1; a surge model construction module, configured to construct a sea surface surge model using the target sea spectrum as a base spectrum and in combination with the target direction expansion function; The target direction spread function determination module is specifically configured to determine the target direction spread function according to the following formula: : ； ； in, represents the basic direction expansion function, represents the surge factor function, represents the first regularization function, k Indicates the real-time wave value of the target sea spectrum, θ Indicates the wave number direction angle of the target sea spectrum; The basic direction spread function for: ； ； ； ； in, k p represents the peak wave value of the target sea spectrum, U represents the sea surface wind speed, represents the acceleration due to gravity, represents the second regularization coefficient, represents the gamma function, s and s p is an intermediate variable.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the computer program, the sea surge digital twin modeling method according to any one of claims 1 to 5 is implemented.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the sea surge digital twin modeling method according to any one of claims 1 to 5 is implemented.

Images