Construction method and system for wave height time sequence process of long-period sea waves at in-port berths

Abstract

The invention discloses a construction method and system for a long-period sea wave height time sequence process at a berth in a port. The method comprises the following steps: constructing a wave propagation model; on the basis of actual measurement or reanalysis of wind wave length time series data, a typical-order wave spectrum process curve at the boundary of the model is obtained so as to drive the wave propagation model; on the basis of the driven wave propagation model, wave parameters at the berth are obtained, and the nonlinear fitting relation between the wave parameters at the boundary and the wave height at the berth is obtained through calculation; and acquiring long-time-sequence wave parameters at the boundary, and based on the long-time-sequence wave parameters at the boundary, obtaining gravity wave short-wave and low-frequency long-wave processes of the long-time-sequence sea waves at the berth through the nonlinear fitting relationship. According to the method, the long-time-sequence process of low-frequency long waves and gravity short waves at any berth of a harbor area can be rapidly achieved, dynamic load conditions influencing ship movement at the berth are provided, and basic wave conditions are provided for evaluating the number of days of ship operation loss.

Description

Technical Field

[0001] The invention belongs to the technical field of ocean wave simulation, and in particular relates to a method and a system for constructing a time series process of ocean wave height of a long period at a berth in a port. Background Art

[0002] Wave load is the main dynamic factor affecting the stability of ships in the port. Long-period waves in the port area are a necessary condition for correctly evaluating the ship's operating window and berthing and unberthing window. Especially in sea areas affected by long-period waves, gravity waves and low-frequency waves have different effects on the movement of ships. Low-frequency waves have stronger penetrating power. When encountering breakwaters, part of the wave energy may pass through the breakwater to produce transmitted waves, which in turn affects the berthing conditions in the port and causes a significant increase in the movement of ships. Therefore, quickly obtaining the wave time sequence process at the berth is helpful to evaluate the safety and feasibility of ship mooring under different wave conditions.

[0003] At present, when evaluating the safety of ship operations, it is usually based on wave data observed in short-term engineering areas or reanalyzed wave data, and the time-series wave data process at the berth is obtained using the relative wave height relationship provided in the propagation model. The existing technology makes it difficult to comprehensively consider the different effects of gravity short waves and low-frequency long waves on ship operations. Summary of the invention

[0004] The present invention proposes a method and system for constructing a time series process of long-period sea wave height at a berth in a port, so as to solve the technical problems existing in the above-mentioned prior art.

[0005] To achieve the above object, the present invention provides a method for constructing a time series process of long-period sea wave height at a berth in a port, comprising:

[0006] Construct wave propagation models;

[0007] Based on the wave propagation model, a wave spectrum process curve of a typical group at the boundary is obtained, and based on the wave spectrum process curve, a wave parameter at the boundary is obtained;

[0008] Based on the wave propagation model, wave parameters at the berth are obtained, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth is calculated, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth;

[0009] The long-period wave parameters at the boundary are obtained, and based on the long-time series wave parameters at the boundary, the gravity short waves and low-frequency long waves of the long-period waves at the berth are obtained through the nonlinear fitting relationship.

[0010] Preferably, the wave propagation model is a time domain wave propagation model of wave refraction and port wall reflection.

[0011] Preferably, before obtaining the wave spectrum process curve of the typical group at the boundary, the method further includes:

[0012] Based on the wave propagation model, a typical group with wave height-period-wave direction as characteristic values ​​is obtained;

[0013] Based on the typical groups, a wave spectrum process curve of the typical groups at the boundary is obtained.

[0014] Preferably, the wave spectrum process curve includes: a low-frequency wave spectrum distribution curve and a gravity wave spectrum distribution curve.

[0015] Preferably, the wave parameters at the boundary include: wave height, period and wave direction.

[0016] Preferably, the method further comprises acquiring wave parameters at the detection point, and obtaining a third nonlinear fitting relationship based on the wave parameters at the detection point and the wave parameters at the boundary;

[0017] Based on the wave parameters at the boundary and the third nonlinear fitting relationship, the observation values ​​at the detection point are verified, wherein the observation values ​​at the detection point include: low-frequency wave height observation values ​​and gravity wave height observation values.

[0018] In order to achieve the above technical objectives, the present invention also provides a system for constructing a time series process of long-period sea wave height at a berth in a port, which is used to implement any of the methods for constructing a time series process of long-period sea wave height at a berth in a port, and the system comprises:

[0019] Model building module, used to build wave propagation models;

[0020] A parameter acquisition module, used to obtain a wave spectrum process curve of a typical group at the boundary based on the wave propagation model, and obtain wave parameters at the boundary based on the wave spectrum process curve;

[0021] A fitting calculation module is used to obtain wave parameters at the berth based on the wave propagation model, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; calculate and obtain a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth;

[0022] The parameter fitting module is used to obtain the long-period wave parameters at the boundary, and obtain the gravity short waves and low-frequency long waves of the long-period waves at the berth through the nonlinear fitting relationship based on the long-time series wave parameters at the boundary.

[0023] Compared with the prior art, the present invention has the following advantages and technical effects:

[0024] The invention provides a method for constructing a time series process of wave height of long-period sea waves at a berth in a port, comprising the following steps: firstly, constructing a wave propagation model; secondly, obtaining a wave spectrum process curve of a typical group at a boundary based on the wave propagation model, and obtaining wave parameters at the boundary based on the wave spectrum process curve; then obtaining the wave parameters at the berth, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; calculating and obtaining a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth; finally, obtaining the long-period wave parameters at the boundary, and obtaining the gravity short waves and low-frequency long waves of the long-period sea waves at the berth through the nonlinear fitting relationship based on the long time series wave parameters at the boundary.

[0025] The present invention can quickly realize the long time series process of low-frequency long waves and gravity short waves at any berth in the port area, provide the dynamic load conditions affecting the movement of ships at the berth, which is conducive to the correct evaluation of ship operation windows and berthing and unberthing windows, and provide important support for port design and operation. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] The drawings constituting a part of the present application are used to provide a further understanding of the present application. The illustrative embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0027] Figure 1 Schematic diagram of a wave propagation model from outside the port to the port area according to an embodiment of the present invention;

[0028] Figure 2 The wave height-period-wave direction scatter point distribution diagram and group setting schematic diagram of the embodiment of the present invention, wherein (a) is the effective wave height-effective period scatter point distribution diagram in the embodiment, wherein the black circles represent the typical group distribution schematic diagram; (b) is the effective wave height-wave direction scatter point distribution diagram in the embodiment, wherein the black circles represent the typical group distribution schematic diagram; (c) is the effective period-wave direction scatter point distribution diagram in the embodiment, wherein the black circles represent the typical group distribution schematic diagram;

[0029] Figure 3The wave height-period-wave direction scatter point distribution diagram and group setting schematic diagram of the embodiment of the present invention, wherein (a) is a typical group overall spectrum curve constructed in the embodiment, and (b) is a typical group low-frequency wave spectrum curve constructed in the embodiment;

[0030] Figure 4 A schematic diagram of the locations of inspection points and berths according to an embodiment of the present invention;

[0031] Figure 5 The wave height parameters obtained at the inspection point and the berth under different typical groups of the embodiments of the present invention are as follows: (a) the wave height of gravity waves at the inspection point under different typical groups (unit: cm), (b) the wave height of low-frequency waves at the inspection point under different typical groups (unit: cm), (c) the wave height of gravity waves at the berth under different typical groups (unit: cm), (d) the wave height of low-frequency waves at the berth under different typical groups (unit: cm);

[0032] Figure 6 Schematic diagram of the test results of the simulated wave height and the fitted wave height at the test point of the embodiment of the present invention, wherein (a) the low-frequency wave test results of the test point are obtained under different typical groups, and (b) the gravity wave test results of the test point are obtained under different typical groups;

[0033] Figure 7 Schematic diagram of the comparison and verification process of the fitted wave height and the observed wave height at the check point of the embodiment of the present invention, wherein (a) is the verification process of the low-frequency wave sequence at the check point, and (b) is the verification process of the gravity wave sequence at the check point;

[0034] Figure 8 It is a schematic diagram of the fitting wave height time series process at the berth according to an embodiment of the present invention; wherein (a) is the gravity wave sequence process at the berth, and (b) is the low-frequency wave sequence process at the berth;

[0035] Fig. 9 The following is a flow chart of a construction method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0036] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0037] It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and that, although a logical order is shown in the flowcharts, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0038] Embodiment 1

[0039] like Fig. 9As shown, in this embodiment, a method for constructing a time series process of long-period sea wave height at a berth in a port is provided, comprising:

[0040] Construct wave propagation models;

[0041] Specifically, a wave numerical model of waves propagating from outside the port to inside the port is established, and the corresponding wave conditions of the port area under different incident wave conditions can be obtained from the model. The model needs to consider the diffraction, refraction, and reflection of port buildings and quay walls. The time domain wave propagation model uses open source software FUNWAVE, SWASH or MIKE-BW software, etc., and no specific requirements are made here.

[0042] When establishing a wave propagation model from outside the port to the port area, a refined shoreline file, water depth file and input spectrum file are prepared. Specifically, the simulated port area plan layout shoreline file, water depth file and input parameter conditions are obtained in advance. The port area shoreline file must correctly describe the port area plan layout characteristics and the outline of the shoreline around the port area. The grid generator and file generator are used to generate the wave calculation grid and the input spectrum file in the required format respectively. These files are then input into the port area propagation wave mathematical model to calculate the output wave parameters and water surface line process, where the wave parameter data at least includes wave height and wave period. Figure 1 Shown is an example of the simulation range of a wave propagation model for a project in Peru.

[0043] Based on the wave propagation model, a wave spectrum process curve of a typical group at the boundary is obtained, and based on the wave spectrum process curve, a wave parameter at the boundary is obtained;

[0044] Before obtaining the wave spectrum process curve of the typical group at the boundary, it also includes:

[0045] Based on the wave propagation model, a typical group with wave height-period-wave direction as characteristic values ​​is obtained;

[0046] Based on the typical groups, a wave spectrum process curve of the typical groups at the boundary is obtained.

[0047] Specifically, a scatter density analysis of the three parameters of wave-period-wave direction at the model boundary is carried out. It contains three sub-graphs, wave height-period scatter distribution, wave height-wave direction scatter distribution and wave direction-period scatter distribution. The selection of typical groups from the scatter distribution graph needs to meet the principle of full coverage of wave height-period-wave direction, the period and wave direction corresponding to the maximum wave height, the period and wave direction corresponding to the minimum wave height; the wave height and wave direction corresponding to the maximum period, the wave height and wave direction corresponding to the minimum period; the full wave direction in the boundary wave sequence; the wave height-period scatter density concentration area, the wave height-wave direction scatter density concentration area, and the wave period-wave direction scatter density concentration area.

[0048] Figure 2In the examples shown, the main characteristic values ​​of effective wave height are between 1.00m and 3.50m, the main characteristic values ​​of short wave effective period are between 8s and 20s, and the main characteristic values ​​of wave direction are between 210° and 250°. Combined with the layout characteristics of the project plan, the selected short wave height characteristic values ​​are 0.6m, 1m, 1.2m, 1.4m, 1.5m, 2m, 2.25m, 2.5m, 2.75m, 3m, 3.25m and 3.5m; the short wave period characteristic values ​​are 6s, 9s, 10s, 12s, 14s, 15s, 16s, 18s and 20s; the direction characteristic values ​​are 210°, 230°, 250° and 270°, which form a total of 432 combinations, from which 116 groups are selected as typical groups for the calculation of the wave model propagating from outside the port to inside the port.

[0049] Obtain the spectral shape of the incident wave spectrum under typical groups. The input spectrum includes the spectrum distribution of low-frequency waves and gravity waves. The gravity wave adopts the JONSWAP wave generation method. When the wave frequency is less than 0.03Hz, the low-frequency wave part adopts a uniformly distributed white noise spectrum. According to the spectral characteristics of the sea area where the project area is located, a suitable spectral shape and a corresponding spectral peak factor of 3.3 are adopted.

[0050] Figure 3 The incident wave spectrum curve for the typical group (input Hs = 1.5m, Tm = 12s) is shown. The wave energy distribution when the frequency is in the range of 0.03 to 1Hz is obtained by using the JONSWAP theoretical spectrum calculation method. The spectrum peak factor is selected as 3.3. When the frequency is less than 0.03Hz, the wave energy distribution is obtained by using the white noise spectrum method, as shown in Figure 3 If the wave frequency is greater than 0.03 Hz and the wave energy is lower than the low-frequency wave energy, the low-frequency wave energy is used. Figure 3 As shown in (a) in .

[0051] The construction method also includes obtaining wave parameters at the detection point, and obtaining a third nonlinear fitting relationship based on the wave parameters at the detection point and the wave parameters at the boundary; based on the wave parameters at the boundary and the third nonlinear fitting relationship, verifying the observation values ​​at the detection point, wherein the observation values ​​at the detection point include: low-frequency wave height observation values ​​and gravity wave height observation values.

[0052] Specifically, the wave dynamic propagation process in the port is calculated, and the wave parameters at the inspection point and berth under different calculation groups are obtained. The locations of the inspection point and berth, where wave observations should be carried out at the inspection point to verify and compare the wave height parameters, and the berth should be located in the protected berth area in the port. Multiple berths can be selected for analysis. This demonstration only takes one point for example analysis, where the coordinates of the inspection point are (11.58°S, 77.30°E) and the water depth is -30m, and the coordinates of the berth are (11.59°S, 77.28°E) and the water depth is -15m. Figure 4 shown.

[0053] When carrying out the wave dynamic propagation process in the port for the selected typical groups, the wave height parameter values ​​at the inspection point and berth under different calculation group conditions are obtained. The wave height parameter values ​​include the gravity wave height value and the low-frequency wave height value.

[0054] Figure 5 In the demonstration example, the typical group of incident wave height ranges from 0.6m to 3.5m, the incident period ranges from 6s to 20s, the gravity wave height ranges from 0.2 to 3.5m at the inspection point, the low-frequency wave height ranges from 0 to 0.6m at the inspection point, and is positively correlated with the wave height and period; the gravity wave height ranges from 0 to 0.9m at the berth point, the low-frequency wave height ranges from 0 to 0.5m at the inspection point, and is nonlinearly correlated with the wave height and period.

[0055] The fitting relationship between the input wave parameters and the low-frequency wave height and gravity wave height at the test point is established. First, the wave height H, period T and incident wave direction D at the boundary in the above typical calculation group are used as dependent variables, and the low-frequency wave height (H S_IG ), gravity wave height (H S_SW ) The time series data is used as the independent variable and the nonlinear regression equation is used to establish the third nonlinear fitting relationship. S_IG With H, T, D, H 2 、T 2 and HT multi-parameter low-frequency wave nonlinear regression equation, and establish H S_SW With H, T, D, H 2 、T 2 The multi-parameter gravity wave nonlinear regression equations of HT and HT are shown in the following equations (1) and (2). The goodness of fit is 0.99, the P test value of each parameter fitting is close to 0, and the average value of the standard error of the three parameters is 0.11. S_IG Compared with the predicted low-frequency wave height test results, wave height H S_SW The test results of predicted gravity wave values ​​are shown in Figure 6 (a)~ Figure 6 (b).

[0056] H S_SW=-12.14H+19.49T+3.34D-1.33H 2 -0.49T 2 +2.85HT-490.7 (1)

[0057] H S_IG =0.28H-1.18T+0.26D+0.37H 2 -0.02T 2 +0.92HT-26.8 (2)

[0058] The wave height time series process at the check point is verified and analyzed. The wave height time series process parameters of the boundary input are input into the third nonlinear fitting relationship model mentioned above, and the low-frequency wave height and gravity wave height process are calculated and compared with the observed low-frequency wave height and gravity wave height for verification. Figure 7 The observed wave height at the test point shown is from 2018 / 07 / 07 to 2018 / 10 / 10, and the wave parameter sequence values ​​in this period at the boundary are obtained. The wave height period is input as an input parameter into the third nonlinear fitting relationship model, and the low-frequency wave height value and gravity wave height value at the test point in this period are obtained. The correlation coefficient in the low-frequency wave height verification is 0.93, and the relative error is 8.1%; the correlation coefficient in the gravity wave height verification is 0.92, and the relative error is 10.6%. The verification indicators meet the verification requirements.

[0059] Based on the wave propagation model, wave parameters at the berth are obtained, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth is calculated, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth;

[0060] Specifically, the wave height H and period T at the boundary in the above typical calculation group are used as dependent variables, and the low-frequency wave height (H S_IG ), gravity wave height (H S_SW ) The time series data is used as the independent variable and a nonlinear regression equation is used to establish a nonlinear fitting relationship. S_IG With H, T, D, H 2 The second nonlinear regression equation of the multi-parameter low-frequency wave of HT and HD was established. S_SW With H, T, D, H 2 The first nonlinear regression equation of multi-parameter gravity waves of HT and HD has a goodness of fit of 0.99, the P test value of each parameter fitting is close to 0, and the average value of the standard error of the three parameters is 0.11. S_IGCompared with the predicted low-frequency wave height test results, wave height H S_SW The results are verified against the predicted gravity wave values.

[0061] H S_SW =-63.17H-11.95T-6.42D-0.88H 2 +0.95HT+0.28HD+805.3 (3)

[0062] H S_IG =-8.58H-2.37T-0.88D+0.19H 2 +0.64HT+0.02HD+116.2 (4)

[0063] The long-period wave parameters at the boundary are obtained, and based on the long-time series wave parameters at the boundary, the gravity short waves and low-frequency long waves of the long-period waves at the berth are obtained through the nonlinear fitting relationship.

[0064] Specifically, the wave height time series process parameters inputted from the boundary are inputted into the first nonlinear fitting relationship and the second nonlinear fitting relationship, and the low-frequency wave height and gravity wave height processes at the berth are obtained by calculation. Figure 8 The wave height observation time at the berth is from 1979 / 01 / 01 to 2019 / 12 / 21. The wave parameter sequence values ​​within this period at the boundary are obtained, and the wave height period is input as an input parameter into the first nonlinear fitting relationship and the second nonlinear fitting relationship to obtain the low-frequency wave height value and gravity wave height value at the berth within this time period.

[0065] Compared with the prior art, this embodiment has the following advantages and technical effects:

[0066] This embodiment can quickly realize the long time series process of low-frequency long waves and gravity short waves at any berth in the port area, provide the dynamic load conditions affecting the movement of ships at the berth, and provide basic wave conditions for evaluating the number of days lost in ship operations.

[0067] Embodiment 2

[0068] This embodiment provides a system for constructing a time series process of long-period sea wave height at a berth in a port, which is used to execute the method for constructing a time series process of long-period sea wave height at a berth in a port described in Embodiment 1. The system includes:

[0069] Model building module, used to build wave propagation models;

[0070] A parameter acquisition module, used to obtain a wave spectrum process curve of a typical group at the boundary based on the wave propagation model, and obtain wave parameters at the boundary based on the wave spectrum process curve;

[0071] A fitting calculation module is used to obtain wave parameters at the berth based on the wave propagation model, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; calculate and obtain a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth;

[0072] The parameter fitting module is used to obtain the long-period wave parameters at the boundary, and obtain the gravity short waves and low-frequency long waves of the long-period waves at the berth through the nonlinear fitting relationship based on the long-time series wave parameters at the boundary.

[0073] This embodiment has all the advantages of the method for constructing the time series of long-period sea wave height at a berth in a harbor described in the first embodiment.

[0074] The above is only a preferred specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed in the present application should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for constructing a time series of long-period sea wave height at a berth in a port, characterized in that: include: Construct wave propagation models; Based on the wave propagation model, a wave spectrum process curve of a typical group at the boundary is obtained, and based on the wave spectrum process curve, a wave parameter at the boundary is obtained; Based on the wave propagation model, wave parameters at the berth are obtained, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth is calculated, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth; The long-period wave parameters at the boundary are obtained, and based on the long-period wave parameters at the boundary, the gravity short waves and low-frequency long waves of the long-period waves at the berth are obtained through the nonlinear fitting relationship.

2. The method for constructing the time series of long-period sea wave height at a berth in a port according to claim 1 is characterized in that: The wave propagation model is a time domain wave propagation model of wave refraction and port wall reflection.

3. The method for constructing the time series of long-period sea wave height at a berth in a port according to claim 1 is characterized in that: Before obtaining the wave spectrum process curve of the typical group at the boundary, it also includes: Based on the wave propagation model, a typical group with wave height-period-wave direction as characteristic values ​​is obtained; Based on the typical groups, a wave spectrum process curve of the typical groups at the boundary is obtained.

4. The method for constructing the time series of long-period sea wave height at a berth in a port according to claim 1 is characterized in that: The wave spectrum process curve includes: a low-frequency wave spectrum distribution curve and a gravity wave spectrum distribution curve.

5. The method for constructing the time series of long-period sea wave height at a berth in a port according to claim 4 is characterized in that: The wave parameters at the boundary include: wave height, wave period and wave direction.

6. The method for constructing the time series of long-period sea wave height at a berth in a port according to claim 1 is characterized in that: It also includes obtaining wave parameters at the detection point, and obtaining a third nonlinear fitting relationship based on the wave parameters at the detection point and the wave parameters at the boundary; Based on the wave parameters at the boundary and the third nonlinear fitting relationship, the observation values ​​at the detection point are verified, wherein the observation values ​​at the detection point include: low-frequency wave height observation values ​​and gravity wave height observation values.

7. A system for constructing the time series of long-period sea wave height at a berth in a port, characterized in that: A method for constructing a time series process of long-period sea wave height at a berth in a port according to any one of claims 1 to 6, the system comprising: Model building module, used to build wave propagation models; A parameter acquisition module, used to obtain a wave spectrum process curve of a typical group at the boundary based on the wave propagation model, and obtain wave parameters at the boundary based on the wave spectrum process curve; A fitting calculation module is used to obtain wave parameters at the berth based on the wave propagation model, wherein the wave parameters at the berth include: gravity short waves at the berth and low-frequency long waves at the berth; calculate and obtain a nonlinear fitting relationship between the wave parameters at the boundary and the wave parameters at the berth, wherein the nonlinear fitting relationship includes: a first nonlinear fitting relationship between the wave parameters at the boundary and the gravity short waves at the berth and a second nonlinear fitting relationship between the wave parameters at the boundary and the low-frequency long waves at the berth; The parameter fitting module is used to obtain the long-period wave parameters at the boundary, and obtain the gravity short waves and low-frequency long waves of the long-period waves at the berth through the nonlinear fitting relationship based on the long-time series wave parameters at the boundary.

Images