Event-based method for calculating the repair weight of a water-damaged sea dike facing block

Abstract

The application discloses an event-based water-destroyed sea embankment facing block weight repair calculation method, which comprises the following steps: determining a typical typhoon according to tropical cyclone optimal path data, sea embankment completion time, facing block arrangement and embankment front tide level process, and collecting multi-element coordinated observation data during the typical typhoon; constructing an astronomical tide-storm surge mathematical model and a typhoon wave mathematical model according to the typical typhoon, and then correcting the two mathematical models by using the multi-element coordinated observation data; calculating the embankment front tide level and the typhoon wave corrected by scouring after correction; designing a facing block stable weight parameterization scheme, obtaining a stable calculation weight based on the embankment front tide level and the typhoon wave, and determining an optimal scheme in combination with an actual weight of the facing block; calculating a design wind speed according to meteorological station wind speed data, calculating a design tide level and a design wave corrected by scouring based on a relationship between the tide level data and the wind wave, and obtaining a facing block repair weight by using the optimal scheme. The method has a high guiding effect on coastal protection engineering design.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of coastal protection engineering, and particularly relates to a method for calculating the repair weight of a water-damaged sea embankment facing block based on events. BACKGROUND

[0002] Storm surge and typhoon wave are common marine disasters caused by typhoons, and storm surge can submerge coastal towns. As an important hydraulic structure for resisting storm surge and typhoon wave, the sea embankment plays a vital role in protecting the lives and property of coastal people. The facing block is covered on the outer layer of the embankment core material and is used to directly bear the impact of waves and water flow to protect other parts of the embankment from erosion and damage, and the stability of the facing block is crucial to the safety of the entire sea embankment.

[0003] In the design process, the stability of the facing block is mainly determined by the stability weight calculation formula, and important projects are verified by wave tank tests in the laboratory by making a scaled-down model of the sea embankment. There are many existing calculation formulas and research results for the stability of the facing block, but they are basically obtained through wave tank tests in the laboratory. In addition to the scale effect of the model test itself, the sea embankment model, water, beachfront topography, sea waves and wind speed used in the test are different from the actual situation, and cannot fully reflect the real process of the facing block of the sea embankment receiving the action of typhoon waves during a typhoon. SUMMARY

[0004] The purpose of the embodiments of the present application is to provide a method for calculating the repair weight of a water-damaged sea embankment facing block based on events, so as to solve the technical problems of inaccurate calculation of the stability weight of the facing block and poor applicability of the calculation method in actual sea embankment projects, difficult decision-making and deviation from the actual situation.

[0005] According to a first aspect of the embodiments of the present application, a method for calculating the repair weight of a water-damaged sea embankment facing block based on events is provided, comprising:

[0006] S1: collecting the design completion report, safety appraisal report, tidal station tidal data, meteorological station wind speed data, tropical cyclone best path data and satellite image of the water-damaged sea embankment, determining the sea embankment completion time, facing block arrangement, actual weight of the facing block, pre-embankment tidal process, sea embankment water damage and non-water damage time points;

[0007] S2: determining a typical typhoon according to the tropical cyclone best path data, sea embankment completion time, facing block arrangement and pre-embankment tidal process, collecting multi-element coordinated observation data during the typical typhoon, and classifying the typical typhoon according to the sea embankment water damage and non-water damage time points;

[0008] S3: constructing an astronomical tide-storm surge mathematical model and a typhoon wave mathematical model according to the typical typhoon, and then correcting the two mathematical models by using the multi-element collaborative observation data, and calculating the dike front tide level and the typhoon wave corrected by scouring after correction;

[0009] S4: designing a parameterized scheme of the armor block stable weight according to the typical typhoon type, obtaining the stable calculation weight based on the dike front tide level and the typhoon wave, and calculating the comprehensive decision factor of each scheme after combining the actual weight of the armor block to determine the optimal scheme;

[0010] S5: calculating the design wind speed according to the wind speed data of the weather station, calculating the design tide level and the design wave corrected by scouring based on the relationship between the tide level data and the wind wave, and obtaining the repair weight of the armor block by using the optimal scheme.

[0011] Optionally, the multi-element collaborative observation data includes weather station pressure and wind speed data, ocean buoy wave data, nearshore tide level station tide level data, and dike front geological data, topographic data and vegetation data.

[0012] Optionally, according to the tropical cyclone best track data, the sea dike completion time, the armor block arrangement and the dike front tide level process, a typical typhoon is determined, multi-element collaborative observation data during the typical typhoon is collected, and the typical typhoon is classified according to the water-destroyed and non-water-destroyed time points of the sea dike, including:

[0013] S21: determining all typhoons that occur since the completion of the sea dike by using the tropical cyclone best track data, and determining the typical typhoon and its time process in which the highest tide level in front of the dike is higher than the top elevation of the bottom armor block of the sea dike according to the armor block arrangement and the dike front tide level process;

[0014] S22: collecting the multi-element collaborative observation data during the typical typhoon;

[0015] S23: classifying the typical typhoon into a non-water-destroyed typhoon and a water-destroyed typhoon according to the water-destroyed and non-water-destroyed time points of the sea dike.

[0016] Optionally, an astronomical tide-storm surge mathematical model and a typhoon wave mathematical model are constructed according to the typical typhoon, and then the two mathematical models are corrected by using the multi-element collaborative observation data, and the dike front tide level and the typhoon wave corrected by scouring after correction are calculated, including:

[0017] S31: collecting ERA5 pressure field, ERA5 wind field and offshore astronomical tide prediction process corresponding to the time according to the typical typhoon, comparing the ERA5 pressure field and ERA5 wind field with the collected meteorological station pressure and wind speed data, respectively correcting the ERA5 pressure field and ERA5 wind field by using the theoretical pressure field and wind field to obtain the corrected ERA5 wind field and ERA5 pressure field;

[0018] S32: establishing an astronomical tide-storm surge mathematical model by using the topographic data, inputting the corrected ERA5 pressure field and ERA5 wind field into the astronomical tide-storm surge mathematical model, adopting the offshore astronomical tide prediction process as the model sea area boundary, calculating to obtain a storm surge level, comparing the storm surge level with the nearshore tide station tide level data to correct the offshore astronomical tide prediction process, and recalculating the astronomical tide-storm surge mathematical model until the difference between the calculated tide level and the nearshore tide station tide level data is less than a set threshold, and finally obtaining a dike front tide level in a typical typhoon process;

[0019] S33: establishing a typhoon wave mathematical model by using the topographic data and vegetation data, inputting the corrected ERA5 wind field and the storm surge level into the typhoon wave mathematical model to calculate a typhoon wave, comparing the typhoon wave with the ocean buoy wave data to correct the bottom friction, white cap and breaking coefficient of the typhoon wave mathematical model, and recalculating the typhoon wave mathematical model until the difference between the calculated typhoon wave element and the ocean buoy wave data is less than a set threshold, and finally obtaining a typhoon wave at a half-wavelength position in front of a dike toe in a typical typhoon process;

[0020] S34: calculating a dike toe scouring depth by using the typhoon wave according to the dike front geological data, correcting the topography in front of the sea dike by using the dike toe scouring depth, inputting the corrected topographic data into the typhoon wave mathematical model, repeating S33 and S34 until the absolute value of the difference between the typhoon wave heights at the half-wavelength position in front of the dike toe before and after is less than a set threshold, and finally obtaining a scouring corrected typhoon wave at the half-wavelength position in front of the dike toe in the typical typhoon process.

[0021] Optionally, the calculation method of the dike toe scouring depth is as follows:

[0022] If the dike front beach is a silt coast, the dike toe scouring depth calculation formula is as follows:

[0023] ;

[0024] In the formula, S is the scouring depth; H is the wave height; α is the slope angle of the seaward face of the sea dike; L is the wavelength; d is the water depth;

[0025] If the beach in front of the dike is sandy coast, the dike toe scour depth calculation formula is as follows:

[0026] .

[0027] Optionally, according to the typical typhoon type, a parameterization scheme of the armor block stable weight is designed, the stable calculation weight is obtained based on the tidal level in front of the dike and the typhoon wave, the optimal scheme is determined after the comprehensive decision factor of each scheme is calculated combined with the actual weight of the armor block, including:

[0028] S41: According to the typical typhoon type, one kind of armor block critical stable weight calculation method and one kind of value scheme of its multi-power influence factor are selected respectively to form an armor block stable weight parameterization scheme;

[0029] S42: The seawall is divided into several seawall intervals along the seawall axis direction, and for the identified typical typhoon process, the armor block stable weight parameterization scheme is used to calculate the armor block stable weight of each seawall interval according to the tidal level and the scour corrected typhoon wave of each seawall interval W i ;

[0030] S43: According to all the typical typhoon processes, the comprehensive decision factor of each armor block stable weight parameterization scheme is calculated by using the calculation stable weight W i and the actual weight G of each seawall interval armor block K ;

[0031] S44: The armor block stable weight parameterization scheme corresponding to the maximum comprehensive decision factor value is selected as the optimal parameterization scheme of the armor block stable weight of the corresponding seawall interval.

[0032] Optionally, the comprehensive decision factor K calculation step is:

[0033] S431: According to the water damage typhoon, the selected one kind of armor block stable weight parameterization scheme is used to calculate the result coincidence rate for the unstable segment and the stable segment of the seawall, specifically including: for a seawall interval of the seawall unstable segment, if the calculation stable weight W i of the armor block is greater than the actual weight G of the armor block, it is determined that the block is unstable, and the armor block unstable result and the actual armor block unstable result coincidence rate SW1 are calculated according to all the seawall intervals in the seawall unstable segment; for a seawall interval of the seawall stable segment, if the calculation stable weight W iWD1, if the calculated stability weight of the armor block is less than the actual weight G of the armor block, it is determined that the block is stable, and the stability result of the armor block is calculated according to all sea embankment intervals of the sea embankment stable section, and the coincidence rate of the actual stability result of the armor block is WD1;

[0034] WHL1 = WD1, if the calculated stability weight of the armor block is less than the actual weight G of the armor block, it is determined that the block is stable, and the stability result of the armor block is calculated according to all sea embankment intervals of the sea embankment stable section, and the coincidence rate of the actual stability result of the armor block is WD1; α 1 × SW1 + (1 - 1) × WD1, α 1) × WD1, α 1 is a weight coefficient, if WHL1 is greater than a set threshold value β 1, then the typhoon period process calculated by the armor block stability weight parameterization scheme is assigned 1, otherwise 0, and the sum of the assignments of all the water-destroyed typhoons is S instab ;

[0035] S432: according to the non-water-destroyed typhoon, a selected armor block stability weight parameterization scheme is used, and for a sea embankment interval, if the calculated stability weight of the armor block is less than the actual weight G of the armor block, it is determined that the block is stable, and the stability result of the armor block is calculated according to all sea embankment intervals of the sea embankment stable section, and the coincidence rate of the actual stability result of the armor block is WD2; W i WD2, if the calculated stability weight of the armor block is less than the actual weight G of the armor block, it is determined that the block is stable, and the stability result of the armor block is calculated according to all sea embankment intervals of the sea embankment stable section, and the coincidence rate of the actual stability result of the armor block is WD1;

[0036] WHL2 = WD2, if WHL2 is greater than a set threshold value β 2, then the typhoon period process calculated by the armor block stability weight parameterization scheme is assigned 1, otherwise 0; and the sum of the assignments of all the non-water-destroyed typhoons is S stab ;

[0037] S433: comprehensive decision factor of the stability weight parameterization scheme , N instab and N stab respectively, the number of water-destroyed typhoons and non-water-destroyed typhoons.

[0038] Optionally, the design wind speed is calculated according to the meteorological station wind speed data, the design tide level and the scour corrected design wave are calculated based on the tide level data and the wave-tide relationship, and the repair weight of the armor block is obtained by using the optimal scheme, including:

[0039] S51: according to the meteorological station wind speed data, a P III type distribution frequency analysis method is used to obtain a design wind speed in front of the sea embankment under the sea embankment repair standard;

[0040] S52: Using the tide station tide data, the harmonic analysis method is used to obtain the long-term tide station mean high tide level of the astronomical spring tide, the linear interpolation method is used to calculate the dike front astronomical spring tide average high tide level according to the dike front design wind speed, and the dike front design tide level is calculated based on the local storm surge water level increase relationship;

[0041] S53: According to the dike front design wind speed, the design wave element in the outer sea is determined based on the wave relationship, the nearshore wave propagation mathematical model is established by using the terrain data, and the design wave at the dike front half wavelength is calculated combined with the design tide level;

[0042] S54: According to the dike front geological data, the dike toe scouring depth is calculated by using the design tide level and the design wave, the dike front topographic data is corrected, the corrected topographic data is input into the nearshore wave propagation mathematical model for recalculation, and S53 and S54 are repeated until the absolute value of the difference between the effective wave heights at the dike toe half wavelength before and after two times is less than a set threshold value, and the scouring corrected design wave at the dike toe half wavelength is obtained.

[0043] S55: According to the optimal parameterization scheme of the armor block stable weight, the critical stable weight of the armor block under the sea dike repair standard is calculated by using the design tide level and the scouring corrected design wave.

[0044] According to a second aspect of the embodiments of the present application, an electronic device is provided, comprising:

[0045] One or more processors;

[0046] Memory for storing one or more programs;

[0047] When the one or more programs are executed by the one or more processors, the one or more processors implement the method of the first aspect.

[0048] According to a third aspect of the embodiments of the present application, a computer readable storage medium is provided, which stores computer instructions, and the instructions are executed by a processor to implement the method of the first aspect.

[0049] The technical scheme provided by the embodiments of the present application can include the following beneficial effects:

[0050] As can be seen from the above embodiments, the method for determining the armor block stable weight parameterization scheme based on the stable condition of the armor block in the field under the historical typical typhoon is adopted, the scale effect caused by the scale reduction in the indoor test is overcome, and then the sea dike armor block stable weight calculation is more in line with the actual marine dynamic conditions of the local sea area, and the accuracy of the armor block stable weight calculation is improved.

[0051] The method of using the scouring depth of the dike toe to correct the terrain of the mathematical model of the typhoon wave and the mathematical model of the nearshore wave propagation overcomes the problem of the smaller calculated wave in front of the dike in the previous wave mathematical model due to the non-consideration of the seabed scouring, and further makes the wave model be able to correct the water depth in front of the dike according to the seabed scouring depth, so as to make the calculated results of the typhoon wave and the design wave more consistent with the actual wave elements in front of the seawall during the typhoon period.

[0052] The method of using the design wind speed to calculate the design tide level and the design wave based on the wind wave tide relationship of the local sea area overcomes the problem of the incoordination of the values of the traditional design tide level and the design wave due to the different calculation methods, and further increases the possibility of the design tide level and the design wave occurring in the same typhoon event, and improves the scientificity and rationality of the hydrological parameters of the seawall repair design.

[0053] It should be understood that the foregoing general description and the following detailed description are only exemplary and explanatory, and cannot limit the present application. BRIEF DESCRIPTION OF DRAWINGS

[0054] The accompanying drawings incorporated in the specification hereof and forming a part thereof illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the application.

[0055] Figure 1 is a flow chart of an event-based water-damaged seawall facing block repair weight calculation method provided by the embodiment of the present application.

[0056] Figure 2 is a schematic diagram of a typical typhoon period tide wave process in front of a certain seawall provided by the embodiment of the present application.

[0057] Figure 3 is a wave field distribution diagram of a typical typhoon period in front of a certain seawall provided by the embodiment of the present application.

[0058] Figure 4 is a comparison diagram of the calculated quality and the actual quality of a certain seawall with different facing block stable weight parameterization schemes provided by the embodiment of the present application.

[0059] Figure 5 is a decision factor calculation schematic diagram of a certain seawall facing block stable weight parameterization scheme in different typhoon periods provided by the embodiment of the present application.

[0060] Figure 6 is a decision result of a certain seawall facing block stable weight parameterization scheme provided by the embodiment of the present application.

[0061] Figure 7 is a block diagram of an electronic device provided by the embodiment of the present application. DETAILED DESCRIPTION

[0062] The exemplary embodiments will be described in detail herein with reference to the attached drawings. The description herein relates to the drawings, in which the same numbers represent the same or similar elements, throughout the several views. The following exemplary embodiments are described in order to provide a thorough understanding into the present application. It will be apparent, however, to one skilled in the art that these specific embodiments are not the only ways, in which the present application can be practiced. Rather, they are provided as examples in order to convey the scope of the present application to the skilled man, as detailed in the appended claims.

[0063] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and / or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and / or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0064] It is to be understood that the singular forms "a", "an", and "the" include plural referents unless the context clearly dictates otherwise. It will be further understood that the terms "comprises" and / or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and / or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0065] Figure 1 is a flow chart of a method for calculating the repair weight of a water-damaged revetment block based on events according to an exemplary embodiment, as shown in Figure 1 The method can include the following steps:

[0066] S1: Collecting the design completion report of the water-damaged revetment, the safety appraisal report, the tidal data of the tidal station, the wind speed data of the weather station, the best path data of the tropical cyclone, and the satellite image, determining the completion time of the revetment, the arrangement of the revetment blocks, the actual weight of the revetment blocks, the tidal process in front of the revetment, the water-damaged and non-water-damaged time points of the revetment;

[0067] Specifically, the present embodiment takes a seawall in the eastern coastal area of Zhejiang Province as an example, which is directly facing the open sea area with poor wave conditions. The design completion report of the water-destroyed seawall, the safety appraisal report, the tidal data of the tidal station, the wind speed data of the weather station, the best path data of the tropical cyclone, and the satellite image are collected. According to the design completion report, the seawall completion time, the armor block arrangement, and the actual weight of the armor block 0.55 t are determined. According to the design completion report, the safety appraisal report, and the satellite image, the seawall water-destroyed and non-water-destroyed time points are determined. The seawall has been subjected to multiple typhoons since its completion, and combined with the wind speed data of the weather station and the best path data of the tropical cyclone, it is determined that the damage occurred mainly during the 1918 typhoon “Mina” and the 2016 typhoon “Yan Hua”. According to the tidal data of the tidal station, the pre-dyke tidal process is determined.

[0068] Through the collection of the above-mentioned data, the whole process of the construction and operation of the seawall can be comprehensively understood, and the process and time of the water-destroyed event of the seawall can be understood in detail, thereby providing a basis for the subsequent decision-making of the repair weight calculation method of the seawall armor block.

[0069] S2: According to the best path data of the tropical cyclone, the seawall completion time, the armor block arrangement, and the pre-dyke tidal process, the typical typhoon is determined, the multi-element coordinated observation data during the typical typhoon is collected, and the typical typhoon is classified according to the seawall water-destroyed and non-water-destroyed time points; this step includes the following sub-steps:

[0070] S21: Using the best path data of the tropical cyclone, all typhoons occurring since the seawall completion time are determined, and according to the armor block arrangement and the pre-dyke tidal process, the typical typhoon and its time process in which the highest pre-dyke tidal level is higher than the top elevation of the lowest armor block of the seawall are determined among all typhoons;

[0071] Specifically, the main typhoon process after the construction of the seawall is collected, the typhoon process which has a greater impact on the seawall armor block is identified, and at the same time, the pre-dyke tidal level during the typhoon period is greater than the top elevation of the lowest armor block by 2.50 m, Figure 2 For the typical typhoon period tidal wave process in front of a seawall, the seawall has mainly encountered 1617, 1709, 1718, 1720, 1808, 1814, 1825, 1909, 1918, 2004, 2009, 2106 typhoons since its completion in 2016 to the water-destroyed in 2021, and according to the armor block arrangement and the pre-dyke tidal process, only during the 1808, 1814, 1918, 2106 typhoon, the highest pre-dyke tidal level obviously exceeds the top elevation of the lowest armor block of the seawall, therefore, the 1918 typhoon “Mina”, the 2016 typhoon “Yan Hua”, the 1808 typhoon “Maria”, and the 1814 typhoon “Capricorn” are taken as the typical typhoon.

[0072] S22: Collecting multi-element coordinated observation data during the typical typhoon;

[0073] Specifically, the multi-element coordinated observation data includes meteorological station pressure and wind speed data, ocean buoy wave data, nearshore tide station tide level data, and dike front geological data, topographic data, and vegetation data.

[0074] S23: According to the sea dike water damage and non-water damage time point, the typical typhoon is divided into non-water damage typhoon and water damage typhoon;

[0075] Specifically, according to satellite photos, visits and data collection, the sea dike was intact on September 25, 2016, December 20, 2017, February 14, 2018, and August 28, 2019, and was water damaged on April 21, 2021, December 10, 2019, and July 31, 2021. According to the time of typhoon occurrence, No. 1918 typhoon "Mina" and No. 2016 typhoon "Yan Hua" are listed as water damage typhoons, and No. 1808 typhoon "Maria" and No. 1814 typhoon "Aquarius" are listed as non-water damage typhoons.

[0076] S3: According to the typical typhoon, an astronomical tide-storm surge mathematical model and a typhoon wave mathematical model are constructed, and the multi-element coordinated observation data is used to correct the two mathematical models, and the corrected dike front tide level and scour corrected typhoon wave are calculated; This step includes the following sub-steps:

[0077] S31: According to the typical typhoon, collect ERA5 pressure field, ERA5 wind field and open sea astronomical tide prediction process corresponding to the time, compare the ERA5 pressure field and ERA5 wind field with the collected meteorological station pressure and wind speed data, respectively. Utilize the theoretical pressure field and wind field to correct the ERA5 pressure field and ERA5 wind field to obtain the corrected ERA5 wind field and ERA5 pressure field;

[0078] Specifically, the ERA5 pressure field and wind field corresponding to the typhoon period are downloaded from the European Center for Medium-Range Weather Forecasts, and the Holland theoretical wind field and pressure field are used to correct the ERA5 pressure field and wind field based on the wind speed and pressure of Shipu Station. This makes the corrected ERA5 wind field and ERA5 pressure field more in line with the actual situation, making the subsequent sea dike front storm surge and typhoon wave calculation more accurate.

[0079] S32: Establish an astronomical storm surge mathematical model using the topographic data, input the corrected ERA5 pressure field and ERA5 wind field into the astronomical storm surge mathematical model, use the astronomical tide prediction process of the outer sea as the model sea area boundary, calculate the storm surge level, compare it with the nearshore tide station tide level data, and use it to correct the astronomical tide prediction process of the outer sea, and then recalculate the astronomical storm surge mathematical model until the difference between the calculated tide level and the nearshore tide station tide level data is less than a set threshold, and finally obtain the dike front tide level in a typical typhoon process;

[0080] Specifically, the MIKE 21 developed by the Demark Water Institute is selected as the basis for the coupling calculation of astronomical tide and storm surge, which can arrange the open boundary in the deep water sea area where the astronomical tide and the water level increase are basically linear with coarse grid, and can nest the fine grid in the local area to finely simulate the nearshore sea area with strong water level increase. The large and small calculation domains provide each other with water level and flow rate boundary conditions, and exchange water level and momentum information. It can also better handle the dynamic boundary problem. The astronomical tide level of the sea boundary is calculated by the global tide wave model TPXO.6, which contains 8 main components (M2, S2, N2, K2, K1, O1, P1, Q1) and 2 long-period components (Mf and Mm). The global grid number of the model is 1 440 × 721, and the resolution is 0.25°. The boundary tide wave, sea surface forcing mode and two-dimensional hydrodynamic mode can constitute the coupling mode of astronomical tide and storm surge. The corrected ERA5 pressure field and ERA5 wind field of Typhoon No. 1918 "Mina", Typhoon No. 2016 "Yan Hua", Typhoon No. 808 "Maria" and Typhoon No. 1814 "Capricorn" are input into the astronomical storm surge mathematical model. The calculation range of the astronomical storm surge mathematical model is 21.5°N-41°N, 116.5°E-127°E, and the calculation area is 224 million km 2 , which fully utilizes the advantages of triangular grid. The grid is arranged according to the principle of key water area grid being dense and other water area grid being sparse. The grid layout in the calculation domain takes into account the differences in water flow and topographic gradient. The model has a total of 98 153 units and 52 280 nodes, of which the grid resolution at the outer sea boundary is 50 km, the grid resolution in the coastal area of Zhejiang is 200-500 m, and the grid size in the coastal area in front of the dike is 10 m. By comparing the tide level data of Shipu Station with the calculated values of the astronomical storm surge mathematical model, and by using the outer sea boundary input by TPXO.6, the calculated values of the astronomical storm surge mathematical model are consistent with the measured tide level of Shipu Station, and the error is kept within 5%. Then the dike front tide level in a typical typhoon process is extracted from the astronomical storm surge mathematical model, thereby improving the accuracy of the dike front tide level.

[0081] S33: Establish a mathematical model of typhoon wave using the terrain data and vegetation data, input the corrected ERA5 wind field and the storm tide into the mathematical model of typhoon wave, calculate the typhoon wave, compare it with the ocean buoy wave data, and then correct the bottom friction, white cap and breaking coefficient of the mathematical model of typhoon wave, recalculate the mathematical model of typhoon wave until the difference between the calculated typhoon wave elements and the ocean buoy wave data is less than a set threshold, and finally obtain the typhoon wave at the half-wavelength position in front of the dike during the typical typhoon process;

[0082] Specifically, the third generation SWAN wave spectrum model is used to establish a mathematical model of typhoon wave, which can describe the evolution of wave field under specific wind, flow and underwater topography conditions in shallow water. The control equation is as follows:

[0083]

[0084] In the formula: S term is the source and sink term, , , , , represent energy input generated by wind, wave-wave interaction, breaking wave dissipation, bed loss respectively. The corrected ERA5 wind field of typhoon No. 1918 "Mina", typhoon No. 2016 "Yan Hua", typhoon No. 808 "Maria" and typhoon No. 1814 "Aquarius" is input into the mathematical model of typhoon wave. The calculation results of the mathematical model of typhoon wave are compared with the measured data of Dachen Ocean Station near Zhejiang coast, and the bottom friction, white cap and breaking coefficient of the mathematical model of typhoon wave are adjusted and corrected, so that the error between the calculation results of the mathematical model of typhoon wave and the measured data of Dachen Ocean Station is less than 10%, and finally the typhoon wave in front of the seawall is obtained, as shown in Figure 3 , thereby improving the accuracy of the wave parameters in front of the seawall.

[0085] S34: According to the dike front geological data, the typhoon wave is calculated to calculate the scour depth of the dike toe, the corrected terrain data is input into the mathematical model of typhoon wave, and S33 and S34 are repeated until the absolute value of the typhoon wave height difference at the half-wavelength position in front of the dike toe is less than a set threshold, and finally the corrected typhoon wave at the half-wavelength position in front of the dike toe during the typical typhoon process is obtained.

[0086] Specifically, the stratum distributed in the seawall and its vicinity range is mainly the Holocene plain fill, silty clay mixed silt, silt, silty clay and the upper Pleistocene silty clay, and the beach in front of the seawall is basically silty coast. The silty coast embankment toe scour depth calculation method is used to calculate the scour depth in front of the seawall under the action of typhoon No. 1918 “Mina”, typhoon No. 2016 “Yan Hua”, typhoon No. 808 “Maria” and typhoon No. 1814 “Makara” combined with the tide level and typhoon wave in front of the seawall during the period, and then the terrain in front of the seawall is corrected. The corrected seawall terrain data is input into the typhoon wave mathematical model, and S33 and S34 are repeated until the absolute value of the difference between the wave heights before and after the half wavelength in front of the embankment toe under two typhoons is less than 0.2m, and finally the scour corrected typhoon wave at the half wavelength in front of the embankment toe during the process of typhoon No. 1918 “Mina”, typhoon No. 2016 “Yan Hua”, typhoon No. 808 “Maria” and typhoon No. 1814 “Makara” is obtained, which further improves the accuracy of the typhoon wave in front of the embankment.

[0087] The calculation method of the embankment toe scour depth is as follows:

[0088] If the beach in front of the embankment is silty coast, the embankment toe scour depth calculation formula is as follows:

[0089] ;

[0090] In the formula, S is the scour depth; H is the wave height; α is the seawall water surface slope angle; L is the wavelength; d is the water depth;

[0091] If the beach in front of the embankment is sandy coast, the embankment toe scour depth calculation formula is as follows:

[0092] .

[0093] S4: According to the typical typhoon type, a parameterization scheme of the armor block stable weight is designed, the stable calculation weight is obtained based on the tide level in front of the embankment and the typhoon wave, and the comprehensive decision factor of each scheme is calculated after combining the actual weight of the armor block to determine the optimal scheme; this step includes the following sub-steps:

[0094] S41: According to the typical typhoon type, a kind of armor block critical stable weight calculation method and a kind of value scheme of its multi-power influence factor are selected respectively to form an armor block stable weight parameterization scheme;

[0095] Specifically, for Typhoon No. 1918 "Mina", Typhoon No. 2016 "Yan Hua", Typhoon No. 808 "Maria" and Typhoon No. 1814 "Capricorn", a kind of critical stable weight calculation method of facing block and a kind of value scheme of its multi-power influence factor are selected to form a parameterization scheme of facing block stable weight, so as to calculate the stable weight of facing block during each typhoon. The critical stable weight calculation method of facing block includes: Hudson method, Hudson wave direction correction method, Tanpo method and Soviet method. Among them:

[0096] Hudson method:

[0097] In the formula, W is the individual weight of the facing block of the facing layer, block stone; gamma b is the bulk density of prefabricated concrete special-shaped block or block stone; gamma is the water body bulk density; H is the design wave height; m is the slope; K D is the stability coefficient of facing block when the wave acts in the positive direction.

[0098] Hudson wave direction correction method:

[0099] In the formula, k is the test parameter; β is the angle between the wave direction line and the normal line of the longitudinal axis of the slope embankment, which is suitable for 22.5°~67.5°.

[0100] Tanpo method:

[0101] In the formula, K is the special-shaped block system coefficient; L is the wave length, which is obtained by period T According to the wave dispersion relationship.

[0102] Soviet method:

[0103] In the formula: mu is the coefficient, which is 0.0058 for tetrahedral pyramid and 0.0034 for three-column hexagonal block.

[0104] The multi-power influence factor includes wave height, period, wave direction and water density.

[0105] The wave height value method includes: one percent wave height H 1% , two percent wave height H 2% , four percent wave height H4% , 5% wave height H 5% , significant wave height H 13% , average H a .

[0106] The period value method includes: 1% period T 1% , 2% period T 2% , 4% period T 4% , 5% period T 5% , significant wave period T 13% , average period T a , spectral peak period T P .

[0107] The wave direction value method includes: average wave direction, most unfavorable wave direction, wave direction at maximum wave height moment.

[0108] The water body unit weight value method includes: fresh water unit weight gamma 淡水 , seawater unit weight gamma 海水 , typhoon period high sediment content seawater unit weight gamma 海水 + S 1x g . S 1 is the sediment content of the water body in front of the dike during the typhoon period.

[0109] S42: Divide the seawall along the seawall axis direction into several seawall intervals, for the identified typical typhoon process, according to the tide level and scour corrected typhoon wave of each seawall interval, adopt the armor block stable weight parameterization scheme to calculate the armor block stable weight of each seawall interval W i ;

[0110] Specifically, the seawall is divided every 10m, and the stable weight of each interval during typhoon No. 1918 "Mina", typhoon No. 2016 "Yan Hua", typhoon No. 808 "Maria" and typhoon No. 1814 "Capricorn" is calculated by using the armor block stable weight parameterization scheme W i , see Figure 4 .

[0111] S43: According to all the typical typhoon processes, the calculated stable weight of the armor block of each seawall interval is used Wi and the actual weight G Calculate the comprehensive decision factor of each armor block stable weight parameterization scheme K ;

[0112] Specifically, for Typhoon No. 1918 "Mina", Typhoon No. 2016 "Yan Hua", Typhoon No. 1808 "Maria" and Typhoon No. 1814 "Capricorn", the comprehensive decision factor of the armor block stable weight parameterization scheme is calculated, see Figure 5 .

[0113] The comprehensive decision factor K The calculation steps are:

[0114] S431: According to the water-destroyed typhoon, a selected armor block stable weight parameterization scheme is used to calculate the result fitting rate for the unstable section and the stable section of the seawall, specifically including: for a seawall section of the seawall unstable section, if the calculated stable weight of the armor block W i is greater than the actual weight G of the armor block, it is determined that the block is unstable, and according to all seawall sections in the seawall unstable section, the armor block unstable result and the actual armor block unstable result fitting rate SW1 are calculated; for a seawall section of the seawall stable section, if the calculated stable weight of the armor block W i is less than the actual weight G of the armor block, it is determined that the block is stable, and according to all seawall sections in the seawall stable section, the armor block stable result and the actual armor block stable result fitting rate WD1 are calculated;

[0115] The comprehensive fitting rate during the water-destroyed typhoon WHL1 = α 1×SW1+(1- α 1)×WD1, α 1 is the weight coefficient, if WHL1 is greater than the set threshold value β 1, the typhoon period process calculated by the armor block stable weight parameterization scheme is assigned 1, otherwise it is assigned 0, and the sum of the assignments of all the water-destroyed typhoons is S instab ;

[0116] S432: According to the non-water-destroyed typhoon, a selected armor block stable weight parameterization scheme is used, and for a seawall section, if the calculated stable weight of the armor block W i is less than the actual weight G of the armor block, it is determined that the block is stable, and according to all seawall sections, the armor block stable result and the actual armor block stable result fitting rate WD2 are calculated;

[0117] The comprehensive fitting rate during the non-water-destroyed typhoon WHL2 = WD2, if WHL2 is greater than the set threshold valueβ 2, then the typhoon period process calculated by the stable weight parameterization scheme of the kind of revetment block is assigned 1, otherwise 0; the sum of the typhoon period assignments of all the non-flooded typhoons is S stab ;

[0118] S433: Comprehensive decision factor of stable weight parameterization scheme , N instab and N stab are the number of flooded typhoons and non-flooded typhoons, respectively.

[0119] S44: Select the stable weight parameterization scheme of the revetment block corresponding to the maximum comprehensive decision factor value as the optimal parameterization scheme of the stable weight of the revetment block in the corresponding sea embankment section.

[0120] Specifically, the final comprehensive decision factor is sorted in descending order, and the decision factor K is the maximum of 84%, see Figure 6 , and the corresponding stable weight parameterization scheme of the revetment block is selected as the optimal stable weight parameterization scheme of the revetment block. After typhoon "Yan Hua" No. 2106, only a small part of the revetment block was damaged, and no large-scale damage occurred, so the stable weight of the revetment block is basically around the actual weight of the block. The calculation result of the optimal stable weight parameterization scheme of the revetment block is basically around the actual weight of the block, thereby proving the correctness of the decision method.

[0121] S5: Calculate the design wind speed according to the wind speed data of the meteorological station, calculate the design tide level and scour-corrected design wave based on the relationship between the tide level data and the wave, and obtain the repair weight of the revetment block using the optimal scheme; This step includes the following sub-steps:

[0122] S51: According to the wind speed data of the meteorological station, using P III type distribution frequency analysis method, the design wind speed in front of the dike under the repair standard of the dike is obtained;

[0123] Specifically, the wind speed data of Shipu Station around the dike since its establishment in 1956 is collected, and P III type distribution frequency analysis method is used, wherein the P-III type probability density function is:

[0124]

[0125] In the formula: ; ; . The design wind speed under the repair standard of the dike is calculated to be 25.9 m / s.

[0126] S52: Using the tide station tide data, the harmonic analysis method is used to obtain the long-term tide station astronomical spring average high tide level, and the linear interpolation method is used to calculate the dike front astronomical spring average high tide level according to the tide wave propagation characteristics, and then according to the design wind speed in front of the dike, the design tide level in front of the dike is calculated based on the wind speed increase water relationship in the local storm surge increase water;

[0127] Specifically, according to the location of the seawall, combined with the measured tide level of the upstream Pigeon Head Tide Station and the downstream Shipu Tide Station, according to the tide wave propagation characteristics, the linear interpolation method is used to obtain the astronomical spring average high tide level in front of the seawall, and according to the design wind speed of 25.9m / s under the repair standard of the seawall, based on the wind speed increase water relationship in the local storm surge increase water,

[0128]

[0129] In the formula: zeta is the storm surge increase water; tau a is the wind shear stress; l is the horizontal scale of the sea area; d is the water depth, rho a , rho w is the air and seawater density; C d is the wind drag coefficient; V is the wind speed; g is the gravitational acceleration; k is the empirical coefficient. The design tide level in front of the dike under the repair standard is 4.73m, which considers that the high tide level in front of the seawall is mainly caused by the wind stress during the typhoon period, which is the main reason for inducing shallow water storm surge, and can better coordinate with the wind speed and wave elements of typhoon.

[0130] S53: According to the design wind speed in front of the dike, the design wave elements in the outer sea are determined based on the wind wave relationship, a nearshore wave propagation mathematical model is established using the terrain data, and the design wave in front of the dike at half wavelength is calculated combined with the design tide level;

[0131] Specifically, the wind wave relationship is as follows:

[0132]

[0133]

[0134] In the formula: H 1 / 10 - one-tenth wave height; V - wind speed; F is the wind zone length; A and B - coefficient;C 、 D and E The coefficient is 0.25, 0.7, 0.7. According to the design wind speed of 25.9 m / s under the repair standard of the seawall, the design wave elements of the outer sea are determined. The nearshore wave propagation mathematical model is established by using the terrain data, and the design wave at the half-wavelength in front of the seawall is calculated by combining the design tide of 4.73 m under the repair standard in front of the seawall.

[0135] S54: According to the geological data in front of the seawall, the design tide and the design wave are used to calculate the scour depth of the seawall toe, and the terrain data in front of the seawall is corrected. The corrected terrain data is input into the nearshore wave propagation mathematical model for recalculation, and S53 and S54 are repeated until the absolute value of the difference between the effective wave heights at the half-wavelength in front of the seawall before and after is less than a set threshold value, and the scour corrected design wave at the half-wavelength in front of the seawall is obtained.

[0136] Specifically, the silt coast seawall toe scour depth calculation method is used, the design tide of 4.73 m under the repair standard in front of the seawall and the design wave are used to calculate the seawall toe scour depth to correct the model, until the absolute value of the difference between the effective wave heights before and after at the half-wavelength in front of the seawall is less than 0.05 m, and the scour corrected design wave H 13% at the half-wavelength in front of the seawall is 4.00 m. The design wave considering the scour of the seawall toe is more in line with the actual situation that the water depth in front of the seawall increases due to the scour of the seawall toe during typhoon period, and finally the wave height increases.

[0137] S55: According to the optimal parameterization scheme of the stable weight of the armor block, the design tide and the scour corrected design wave are used to calculate the critical stable weight of the armor block under the repair standard of the seawall;

[0138] Specifically, the design tide under the repair standard is 4.73 m, the wave height H 4% is 4.43 m, the effective wave height H 13% is 4.00 m, and the average period T is 14 s. The repair weight of the seawall armor block is 2.24 t, which is 0.55 t higher than the original actual weight of the seawall armor block. The wind and wave resistance of the seawall armor block is greatly improved, and the seawall has not been damaged obviously since the repair.

[0139] ​From the above embodiments, it can be seen that the application adopts the historical typhoon to evaluate the calculation method of the repair weight of the water-destroyed seawall facing block, so that the calculation scheme of the stable weight of the facing block is more in line with the actual situation of the seawall facing block; at the same time, the technical scheme of calculating the design tide and the design wave by using the design wind speed makes the design tide and the design wave coordinated, which conforms to the actual situation of the wind-induced storm surge and the typhoon wave in the typhoon period of the coastal sea in the southeast coastal area of China; the calculation process of the actual typhoon wave and the design wave in front of the seawall considers the scouring of the seawall toe, which solves the problem of the previous design wave being too small to make the stable weight calculation of the facing block too small and unstable.

[0140] Correspondingly, the application also provides an electronic device, comprising: one or more processors; a memory for storing one or more programs; when the one or more programs are executed by the one or more processors, the one or more processors implement the event-based water-destroyed seawall facing block repair weight calculation method as described above. As shown in the Figure 7 The device based on the event-based water-destroyed seawall facing block repair weight calculation method provided by the embodiment of the application is a hardware structure diagram of any data processing capable device. In addition to the processor, the memory, the DMA controller, the disk, and the non-volatile memory shown in the Figure 7 In addition to the processor, the memory, the DMA controller, the disk, and the non-volatile memory shown in the

[0141] Correspondingly, the application also provides a computer readable storage medium having computer instructions stored thereon, which are executed by a processor to implement the event-based water-destroyed seawall facing block repair weight calculation method as described above. The computer readable storage medium can be an internal storage unit of any data processing capable device, such as a hard disk or a memory. The computer readable storage medium can also be an external storage device, such as a plug-in hard disk, a smart media card (SMC), an SD card, a flash card, etc. Further, the computer readable storage medium can include both the internal storage unit of any data processing capable device and the external storage device. The computer readable storage medium is used to store the computer program and other programs and data required by the data processing capable device, and can also be used to temporarily store data that has been output or will be output.

[0142] Other embodiments of the application will be apparent to those skilled in the art from consideration of the specification and practice of the application disclosed herein. It is intended that the application be construed as including any patern, process, method, technique, composition of matter, or device directly or indirectly

[0143] It is to be understood that the application is not limited to the precise construction hereinafter described and as shown in the attached drawings, and that various changes in the details and combinations of the application can be made without departing from the scope of the application or any of the embodiments thereof. The scope of the application is limited only by the claims appended hereto.

Claims

1. An event-based method for calculating the repair weight of a water-damaged sea dike revetment block, characterized by, The method comprises the following steps: S1: collecting the design completion report of the water-destroyed seawall, the safety appraisal report, the tidal level data of the tidal station, the wind speed data of the weather station, the best path data of the tropical cyclone, and the satellite image, determining the seawall completion time, the armor block arrangement, the actual weight of the armor block, the pre-dike tidal level process, the seawall water-destroying and non-water-destroying time point; S2: determining a typical typhoon according to the best path data of the tropical cyclone, the seawall completion time, the armor block arrangement and the pre-dike tidal level process, collecting multi-element coordinated observation data during the typical typhoon, and classifying the typical typhoon according to the seawall water-destroying and non-water-destroying time point; S3: constructing an astronomical tide-storm surge mathematical model and a typhoon wave mathematical model according to the typical typhoon, and then modifying the two mathematical models by using the multi-element coordinated observation data, and calculating the pre-dike tidal level and the scour-modified typhoon wave after the modification; S4: designing a parameterized scheme of the stable weight of the armor block according to the type of the typical typhoon, obtaining the stable calculation weight based on the pre-dike tidal level and the typhoon wave, and calculating the comprehensive decision factor of each scheme by combining the actual weight of the armor block to determine the optimal scheme; S5: calculating the design wind speed according to the wind speed data of the weather station, calculating the design tidal level and the scour-modified design wave based on the relationship between the tidal level data and the wind wave, and obtaining the repair weight of the armor block by using the optimal scheme; wherein the astronomical tide-storm surge mathematical model and the typhoon wave mathematical model are constructed according to the typical typhoon, and then the two mathematical models are modified by using the multi-element coordinated observation data, and the pre-dike tidal level and the scour-modified typhoon wave are calculated after the modification, which comprises the following steps: S31: collecting the ERA5 pressure field, the ERA5 wind field and the astronomical tide prediction process of the outer sea corresponding to the time according to the typical typhoon, comparing the ERA5 pressure field and the ERA5 wind field with the collected pressure and wind speed data of the weather station, respectively modifying the ERA5 pressure field and the ERA5 wind field by using the theoretical pressure field and the wind field, and obtaining the modified ERA5 wind field and the ERA5 pressure field; S32: establishing an astronomical tide-storm surge mathematical model by using the terrain data, inputting the modified ERA5 pressure field and the ERA5 wind field into the astronomical tide-storm surge mathematical model, using the astronomical tide prediction process of the outer sea as the model sea area boundary, calculating the storm tide level, comparing the storm tide level with the nearshore tidal station tidal level data, modifying the astronomical tide prediction process of the outer sea, recalculating the astronomical tide-storm surge mathematical model, and finally obtaining the pre-dike tidal level in the typical typhoon process until the difference between the calculated tidal level and the nearshore tidal station tidal level data is less than a set threshold value; S33: establishing a typhoon wave mathematical model by using the terrain data and the vegetation data, inputting the modified ERA5 wind field and the storm tide level into the typhoon wave mathematical model, calculating the typhoon wave, comparing the typhoon wave with the ocean buoy wave data, modifying the bottom friction, the white cap and the breaking coefficient of the typhoon wave mathematical model, recalculating the typhoon wave mathematical model, and finally obtaining the typhoon wave at the half-wavelength position in front of the seawall in the typical typhoon process until the difference between the calculated typhoon wave element and the ocean buoy wave data is less than a set threshold value. S34: According to the embankment front geological data, the embankment toe scour depth is calculated by using the typhoon wave, the embankment toe scour depth is used to correct the terrain in front of the seawall, and the corrected terrain data is input into the typhoon wave mathematical model, S33 and S34 are repeated until the absolute value of the difference between the two typhoon wave heights at the half wavelength in front of the embankment toe is less than a set threshold, and finally the scour corrected typhoon wave at the half wavelength in front of the embankment toe in the typical typhoon process is obtained.

2. The method of claim 1, wherein, The multi-element cooperative observation data includes meteorological station pressure and wind speed data, ocean buoy wave data, nearshore tide station tide level data, and embankment front geological data, terrain data and vegetation data.

3. The method of claim 1, wherein, According to the tropical cyclone optimal path data, the seawall completion time, the armor block arrangement and the embankment front tide level process, a typical typhoon is determined, multi-element cooperative observation data during the typical typhoon is collected, and the typical typhoon is classified according to the seawall water damage and non-water damage time points, including: S21: Using the tropical cyclone optimal path data, all typhoons occurring since the completion of the seawall are determined, and according to the armor block arrangement and the embankment front tide level process, the typical typhoon and its time process are determined among all typhoons whose highest embankment front tide level is higher than the top elevation of the bottom armor block of the seawall; S22: Collecting multi-element cooperative observation data during the typical typhoon; S23: According to the seawall water damage and non-water damage time points, the typical typhoon is divided into non-water damage typhoon and water damage typhoon.

4. The method of claim 1, wherein, The calculation method of the embankment toe scour depth is: If the embankment front beach is a silt coast, the embankment toe scour depth calculation formula is as follows: ； wherein S is the scour depth; H is the wave height; α is the seawall water face slope angle; L is the wave length; d is the water depth; If the embankment front beach is a sandy coast, the embankment toe scour depth calculation formula is as follows: 。 5. The method of claim 1, wherein, According to the typical typhoon type, a parameterized scheme of armor block stable weight is designed, the stable calculation weight is obtained based on the embankment front tide level and the typhoon wave, and the optimal scheme is determined after calculating the comprehensive decision factor of each scheme combined with the actual weight of the armor block, including: S41: According to the typical typhoon type, a kind of armor block critical stable weight calculation method and a kind of value selection scheme of its multi-power influence factor are selected respectively to form an armor block stable weight parameterized scheme; S42: divide the seawall into several seawall sections along the seawall axis direction, for the identified typical typhoon process, according to the tide level and the scour corrected typhoon wave of each seawall section, the stable weight parameterization scheme of the armor block is used to calculate the stable weight of the armor block of each seawall section W i ; S43: Calculate the computed stable weight of each revetment block for each sea wall section according to all the typical typhoon processes W i and the actual weight G Calculate the overall decision factor for each revetment block stable weight parameterization scheme K ; S44: The armor block stable weight parameterized scheme corresponding to the maximum comprehensive decision factor value is selected as the optimal parameterized scheme of the armor block stable weight corresponding to the seawall interval.

6. The method of claim 5, wherein, The integrated decision factor K The calculation step is: S431: According to the water-destroyed typhoon, a selected one of the stable weight parameterization schemes of the facing block is adopted to calculate the coincidence rate of the results for the unstable section and the non-unstable section of the seawall, specifically including: for a seawall section of the seawall unstable section, if the calculated stable weight of the facing block is greater than the actual weight G of the facing block, it is determined that the block is unstable, and according to all seawall sections in the seawall unstable section, the coincidence rate of the unstable results of the facing block and the actual unstable results of the facing block is calculated as SW1; for a seawall section of the seawall stable section, if the calculated stable weight of the facing block is less than the actual weight G of the facing block, it is determined that the block is stable, and according to all seawall sections in the seawall stable section, the coincidence rate of the stable results of the facing block and the actual stable results of the facing block is calculated as WD1. W i W i ​​ The comprehensive coincidence rate of the water-destroyed typhoon period is WHL1 α 1×SW1+(1- α 1)×WD1, α 1 is a weight coefficient, if WHL1 is greater than a set threshold β 1, then the typhoon period process calculated by the stable weight parameterization scheme of the kind of armor block is assigned 1, otherwise it is assigned 0, and the sum of the assignments of all the water-destroyed typhoons is S instab ; S432: According to the non-water destructive typhoon, a selected one of the armor block stable weight parameterization schemes is used. For a seawall section, if the calculated stable weight of the armor block is less than the actual weight G of the armor block W i The block is determined to be stable. According to all seawall sections, the calculated stable result of the armor block and the actual stable result of the armor block are calculated. The coincidence rate is WD2. The non-water-destroyed typhoon period comprehensive coincidence rate WHL2 = WD2, if WHL2 is greater than a set threshold β 2, then the typhoon period process of the kind of armor block stable weight parameterization scheme is assigned 1, otherwise 0; the sum of all the typhoon period assignments of the non-water-destroyed typhoon is S stab ; S433: Comprehensive decision factor for stabilizing the weight parameterization scheme , N instab and N stab Nw and Nn are the numbers of water- and non-water-destroyed typhoons, respectively.

7. The method of claim 2, wherein, According to the meteorological station wind speed data, the design wind speed is calculated, the design tide level and scour corrected design wave are calculated based on the relationship between the tide level data and the wind wave, and the repair weight of the armor block is obtained by using the optimal scheme, including: S51: According to the meteorological station wind speed data, P III type distribution frequency analysis method is used to obtain the design wind speed in front of the embankment under the seawall repair standard; S52: Using the tide level data of the tide station, the harmonic analysis method is used to obtain the long-term average high tide level of the astronomical spring tide of the tide station, the linear interpolation method is used to calculate the average high tide level of the astronomical spring tide in front of the embankment according to the tide wave propagation characteristics, and then the design tide level in front of the embankment is calculated based on the relationship between the local storm surge and the wind speed surge according to the design wind speed in front of the embankment. S53: according to the design wind speed in front of the dike, determine the design wave element in the outer sea based on the relationship between wind and wave, establish a nearshore wave propagation mathematical model using the topographic data, and calculate the design wave at the half wavelength in front of the dike in combination with the design tide level; S54: according to the geological data in front of the dike, calculate the scour depth at the dike toe using the design tide level and the design wave, correct the topographic data in front of the dike, input the corrected topographic data into the nearshore wave propagation mathematical model for recalculation, repeat S53 and S54 until the absolute value of the difference between the effective wave heights at the half wavelength in front of the dike toe in two consecutive times is less than a set threshold value, and obtain the scour-corrected design wave at the half wavelength in front of the dike toe; S55: according to the optimal parameterization scheme of the stable weight of the armor block, calculate the critical stable weight of the armor block under the repair standard of the seawall using the design tide level and the scour-corrected design wave.

8. An electronic device, comprising: Comprise: one or more processors; a memory for storing one or more programs; when the one or more programs are executed by the one or more processors, so that the one or more processors implement the method of any one of claims 1-7.

9. A computer readable storage medium having stored thereon computer instructions, wherein, The instruction is executed by the processor to implement the method of any one of claims 1-7.

Images