A rapid calculation method for the wave dissipation capacity of castle-shaped oyster reefs

Abstract

This invention belongs to the field of hydraulic engineering and discloses a rapid calculation method for the wave dissipation capacity of castle-shaped oyster reefs. It primarily targets closely arranged and submerged castle-shaped oyster reefs, using relative submergence degree of the reef, relative reef crest width, incident wave steepness, and oyster density as independent variables. A closed-form empirical formula for the transmission coefficient is obtained through multivariate nonlinear regression. The steps include data acquisition and conversion, parameter domain checking, rapid calculation, and output recording. The applicable conditions are ds / Hi>0 and B / Hi≤25. Its advantages are: this invention achieves rapid calculation of wave dissipation capacity through closed-form substitution calculation, with a single calculation completed in milliseconds, eliminating the need for complex numerical simulations or physical model experiments, and completely solving the pain points of high cost and long cycle of traditional assessment methods.

Description

Technical Field

[0001] This invention belongs to the field of water conservancy engineering, and relates to the field of coastal ecological restoration and slope protection engineering, especially a rapid calculation method for the wave dissipation capacity of castle-shaped oyster reefs. Background Technology

[0002] Oyster reefs are a key component of marine ecosystems, not only purifying water, promoting fishery resource enhancement, increasing biodiversity and ecosystem stability, but also playing important roles in wave and current mitigation and shoreline erosion reduction. They are a core means of achieving synergistic effects between coastal zone ecological restoration and disaster prevention and mitigation. Currently, over one-third of coastal zone protection and restoration projects include oyster reef restoration projects.

[0003] However, existing assessments of the wave-damping effectiveness of castle-shaped oyster reefs have significant shortcomings:

[0004] The evaluation methods are costly and time-consuming: they rely on numerical simulation (which requires specialized software and computing power) or physical model testing (which requires building a water tank and debugging waves), and cannot meet the needs of rapid comparison of multiple options in the early stages of engineering.

[0005] Existing formulas have poor applicability: Traditional wave-dissipating empirical formulas are derived from regular submerged breakwaters or non-ecological components (such as concrete blocks), and do not take into account the ecological characteristics of oyster reefs (such as oyster attachment density), making it difficult to accurately reflect their wave-dissipating mechanisms.

[0006] Therefore, there is an urgent need to develop a rapid calculation method for the wave dissipation capacity of castle-shaped oyster reefs that does not rely on complex simulations and incorporates key factors such as the relative submergence degree of the reef and oyster density, so as to support the early-stage comparison of multiple schemes and real-time on-site evaluation of the project. Summary of the Invention

[0007] The purpose of this invention is to solve the above-mentioned problems existing in the prior art and to provide a rapid calculation method for the wave dissipation capacity of castle-shaped oyster reefs.

[0008] To achieve the above objectives, the technical solution adopted by the present invention is as follows:

[0009] A rapid calculation method for the wave-dissipating capacity of castle-shaped oyster reefs includes the following steps:

[0010] S1. Collect incident wave and reef parameters: including incident wave height Hi, incident wavelength Li, and distance from the reef top to the still water surface. d s reef top width B Oyster attachment density ;

[0011] S2. Convert the parameters into four dimensionless independent variables: relative submergence of the reef ds / Hi, relative reef crest width B / Hi, incident wave steepness Hi / Li, and oyster density. 0;

[0012] S3. Perform parameter domain checks on the dimensionless independent variables ds / Hi and B / Hi of the reef body relative submergence.

[0013] S4. If the parameter domain check passes, substitute the four dimensionless independent variables into the closed empirical formula for the transmission coefficient obtained by multivariate nonlinear regression to obtain the predicted value of the transmission coefficient Kt.

[0014] S5. Output and record the calculation results, including the set of input parameters, the value of Kt, and the conclusion of the parameter domain check.

[0015] Preferably, the condition for the parameter domain check in step S3 is: the relative submergence degree of the reef ds / Hi>0 and the relative width of the reef top B / Hi≤25; if the condition is not met, an "uncertainty prompt" is output.

[0016] Preferably, the empirical formula for the closed-loop transmission coefficient in step S4 is:

[0017]

[0018] The formula is derived from regression analysis of 73 physical model experiments, with a goodness of fit R² = 0.768 and a significance level < 0.05, and can stably predict Kt within the parameter domain.

[0019] Preferably, the oyster attachment density in step S5 This is an ecological parameter, measured in units per square meter, used to reflect the ecological wave-dissipating characteristics of castle-shaped oyster reefs.

[0020] Preferably, the predicted transmission coefficient Kt ranges from 0 to 1, and the smaller Kt is, the better the wave-damping effect of the oyster reef.

[0021] Preferably, the rapid calculation is a closed-loop substitution operation, and the completion time is independent of the input size, which is suitable for sensitivity scanning and rolling comparison of multiple schemes in the early stage of engineering.

[0022] Another aspect of the present invention is to provide a rapid calculation system for the wave-dissipating capacity of castle-shaped oyster reefs.

[0023] To achieve the above objectives, the technical solution adopted by the present invention is as follows:

[0024] A rapid calculation system for the wave-damping capacity of castle-shaped oyster reefs includes a processor and a memory, the memory storing instructions that, when executed on the processor, implement all the steps of the method.

[0025] Due to the adoption of the above technical solution, the beneficial effects obtained by the present invention include:

[0026] 1. This invention achieves rapid calculation of wave dissipation capacity through closed-loop substitution calculation, with a single calculation completed in milliseconds. It does not rely on complex numerical simulations or physical model tests, completely solving the pain points of high cost and long cycle of traditional evaluation methods. This method is suitable for sensitivity scanning and rolling comparison of multiple schemes in the early stage of engineering. The completion time is independent of the input scale, which greatly shortens the scheme optimization cycle and provides efficient decision support for coastal ecological restoration projects.

[0027] 2. The formula of this invention is based on regression calibration of 73 sets of physical model experiments, with a goodness of fit R²=0.768 and statistical significance (P<0.05), and the results are reliable and verifiable. Parameter domain verification (relative submergence degree of reef ds / Hi>0, relative reef top width B / Hi≤25) avoids extrapolation of working conditions and reduces the risk of misuse. The calculation results conform to physical trends (e.g., as the reef top width increases → the transmission coefficient Kt decreases), ensuring consistency with actual engineering laws and providing accurate reference for scheme design.

[0028] 3. This invention is the first to measure oyster attachment density. Incorporating it as a core independent variable into wave dissipation calculations accurately reflects the ecological wave dissipation mechanism of castle-shaped oyster reefs, distinguishing it from the limitations of traditional non-ecological component formulas. The variable selection covers key controlling factors such as submersion degree, reef width, wave steepness, and density, which are highly consistent with engineering physical trends, facilitating rapid directional judgment and parameter fine-tuning, and filling the technical gap in rapid assessment of oyster reef wave dissipation capacity in ecological revetment engineering. Attached Figure Description

[0029] Figure 1 This is a flowchart illustrating the operation of the rapid calculation method for the wave-dissipating capacity of castle-shaped oyster reefs in this invention.

[0030] Figure 2 This is a schematic diagram of an embodiment of the castle-shaped oyster reef in this invention.

[0031] Figure 3 This is a parameter definition and geometric diagram of the present invention.

[0032] Figure 4 This is a comparison chart of the measured and predicted values ​​of the transmission coefficient of the closely arranged castle-shaped oyster reefs of this invention.

[0033] Figure 5 The present invention relates to different oyster attachment densities on reefs. Schematic diagram. Detailed Implementation

[0034] like Figure 1-5As shown, this invention proposes a rapid calculation method for the transmission coefficient Kt (Kt ranges from 0 to 1, with a smaller Kt indicating a better wave-dissipating effect) of closely arranged castle-like oyster reefs under submerged conditions. This method is applicable to closely arranged working conditions. It achieves rapid engineering estimation of input and output through a closed empirical formula for the transmission coefficient obtained by multivariate nonlinear regression, avoiding the high cost and time delay caused by frequent numerical or physical model evaluations during the scheme comparison stage.

[0035] This method uses a closed-form calculation of the transmission coefficient Kt based on a multivariate nonlinear regression empirical formula obtained from physical model experiments. The method incorporates the reef's relative submergence degree ds / Hi, relative reef crest width B / Hi, incident wave steepness Hi / Li, and oyster density. 0 is used as the independent variable.

[0036] The empirical formula fitting form for determining the transmission coefficient of closely packed castle-like oyster reefs is shown in Equation (1).

[0037]

[0038] In the formula: a1, b1, c1, d1, e1, f1 are parameters; Kt is the transmission coefficient; ds is the distance from the top of the oyster reef to the still water surface, which is taken as a positive value when submerged, in meters (m); B is the width of the reef top, in meters (m); Hi is the incident wave height, in meters (m); Li is the incident wavelength, in meters (m). Oyster attachment density, in units per square meter (ind / m2).

[0039] Based on 73 sets of physical model experiments, regression analysis was performed, and the parameter values ​​and 95% confidence intervals are shown in the table below. The goodness of fit R² = 0.768, and the significance < 0.05, indicating that Kt can be stably predicted within the given parameter range.

[0040]

[0041] Finally, the prediction model for the transmission coefficient of closely packed castle-like oyster reefs, determined based on multivariate nonlinear regression analysis, is shown in Equation (2).

[0042] When ds / Hi>0 and B / Hi≤25:

[0043]

[0044] Where a1=0.087, b1=-0.006, c1=0.078, d1=-0.001, e1=-1.107, f1=0.814; the empirical formula for multivariate nonlinear regression is ds / Hi, B / Hi, Hi / Li, 0 is the independent variable, corresponding to the transmission coefficient Kt when the castle-shaped oyster reefs are closely arranged and submerged;

[0045] The physical meaning and dimensions of the independent variables are as follows: ds is the distance from the top of the oyster reef to the still water surface, which is taken as a positive value when submerged, in meters (m); B is the width of the reef top, in meters (m); Hi is the incident wave height, in meters (m); Li is the incident wavelength, in meters (m). Oyster attachment density, in units per square meter (ind / m2).

[0046] Specifically, the method includes the following steps:

[0047] Collect incident wave and reef parameters, including incident wave height Hi, incident wavelength L, distance ds from reef top to still water surface (positive for submersion), reef top width B, and oyster attachment density. According to the model variable definitions, four dimensionless independent variables were obtained: relative submergence degree of the reef ds / Hi, relative reef crest width B / Hi, incident wave steepness Hi / L, and oyster density. 0. The physical meaning and geometric representation of the above variables. Figure 2 As shown; the dimensions of the variables are kept consistent to ensure consistency and verifiability in the calculations.

[0048] Parameter domain check, for ds / Hi, B / Hi, Hi / L, 0. Perform an applicability check to confirm whether it is within the calibrated parameter range; if it meets the parameter range conditions, set ds / Hi, B / Hi, Hi / L, Substituting 0 into the transmission coefficient calculation formula yields the predicted transmission coefficient value Kt; the preset applicable conditions selected by this formula include... d s / H i >0, B / H i ≤25; the calculation formula is used when this condition is met; when any independent variable exceeds the parameter range, an "uncertainty warning" is given in the output to avoid extrapolating the fast result to inapplicable working conditions.

[0049] Quickly calculate ds / Hi, B / Hi, Hi / L, Substituting these values ​​into the transmission coefficient calculation formula yields the predicted transmission coefficient value Kt. This rapid calculation model is a closed-form expression, independent of numerical simulation and real-time physical model solving, and can be directly used for early-stage multi-scheme comparison (suitable for early-stage engineering sensitivity scanning and rolling comparison of multiple schemes).

[0050] Output and Record: Generate a calculation record, including: the set of input parameters and Kt; this record is archived with the design file for easy review and traceability.

[0051] The implementation results of this method satisfy the following trend consistency check: when other variables remain unchanged, an increase in B / Hi or an increase in Hi / L causes a consistent response direction in which Kt decreases; under the condition corresponding to an increase in water depth, The effect of 0 on Kt diminishes marginally; this method is limited to the case of tightly packed castle-like oyster reefs in a submerged state, and does not involve other layouts.

[0052] Figure 3 This diagram illustrates parameter definitions and geometric representations. In the diagram, Hr, Tr, Ht, Tt, and Ti are common terms in hydraulic engineering and wave dynamics. Hr is the wave height of the reflected wave, referring to the height of the wave that returns to the original medium (incident side) after the incident wave encounters the castle-shaped oyster reef without passing through the reef. Tr is the period of the reflected wave, referring to the time required for two adjacent wave crests or troughs to pass through the same fixed point. Ht is the wave height of the transmitted wave, referring to the wave height formed behind the reef (transmission side) after the incident wave passes through the castle-shaped oyster reef. The ratio of wave height Ht to incident wave height Hi is called the transmission coefficient Kt (Kt=Ht / Hi), which reflects the degree to which the reef weakens wave energy (the smaller Kt is, the better the wave dissipation effect); Tt is the period of the transmitted wave, which refers to the time required for two adjacent wave crests or troughs to pass through the same fixed point; Ti is the period of the incident wave, which refers to the time required for two adjacent wave crests or troughs to pass through the same fixed point; ds is the distance from the top of the oyster reef to the still water surface; hs is the distance from the top of the oyster reef to the bottom of the water; d is the distance from the still water surface to the bottom of the water.

[0053] Figure 4 The figure shows a comparison between the measured and predicted values ​​of the transmission coefficient of closely arranged castle-like oyster reefs. As can be seen from the figure, the measured and predicted values ​​fit well, and the formula used in this rapid calculation method can provide good guidance for engineering practice.

[0054] Specifically, here's an example of how the above formula can be applied:

[0055] Input parameters: Hi=0.5m, Li=5m, ds=0.2m, B=1m =100ind / m²

[0056] Step 1 Conversion: ds / Hi = 0.4, B / Hi = 2, Hi / Li = 0.1 0=100

[0057] Step 2 check: ds / Hi=0.4>0, B / Hi=2≤25 → Meets the conditions

[0058] Step 3 Calculation:

[0059] Output result:

[0060] Kt≈0.815 (moderate wave-damping effect);

[0061] Parameter domain check: Meets the requirements;

[0062] Archive: Records all input parameters and calculation results.

[0063] It should be noted that:

[0064] Parameter domain restrictions:

[0065] The following conditions must be met: ds / Hi>0 (submerged state) and B / Hi≤25 (relative reef width);

[0066] When the value exceeds the limit, output an "uncertainty warning" and prohibit the use of formula calculations.

[0067] Ecological parameter collection:

[0068] 0 (oyster density) should be the actual attachment density under tight arrangement (field sampling or design value), unit ind / m²;

[0069] like If 0 represents 0 (no oysters attached), the formula still applies, but the wave-damping effect will be significantly reduced (consistent with trend verification). 0 (Increase → Kt decreases).

[0070] Trend verification:

[0071] Verify whether the calculation results conform to physical laws:

[0072] B / Hi increases → Kt decreases (the wider the reef top, the better the wave dissipation);

[0073] As Hi / Li increases, Kt decreases (the steeper the wave, the better the wave reduction).

[0074] Water depth increases → The impact of 0 on Kt weakens (marginal effect of ecological characteristics).

[0075] Based on existing technologies, this application breaks through the bottleneck of "ecological reef wave dissipation calculation relying on simulation / experiment". By incorporating ecological parameters, constructing closed formulas, and limiting parameter domains, it fills the technical gap in the rapid calculation of wave dissipation capacity of castle-shaped oyster reefs, and realizes rapid and accurate assessment of the wave dissipation capacity of castle-shaped oyster reefs.

[0076] In addition, this invention also discloses a system for executing the above-described method, including a processor and a memory. The memory stores instructions that, when executed on the processor, implement all the steps of the method. The core function of this system is to provide an automated and efficient execution platform for the rapid calculation method of the wave-dissipating capacity of castle-shaped oyster reefs, specifically in the following aspects:

[0077] 1. Automated full-process execution: By running instructions in the memory through the processor, it automatically completes all steps of parameter acquisition and dimensionless conversion, parameter domain checking, closed-form formula calculation, result output and recording, replacing manual operation, avoiding human error and ensuring the standardization of the calculation process;

[0078] 2. Achieve rapid evaluation in milliseconds: Based on the efficient calculation logic of the empirical formula for closed transmission coefficient, the system can complete a single wave dissipation capability evaluation within milliseconds, supporting sensitivity scanning and rolling comparison of multiple schemes in the early stage of engineering (such as quickly comparing the wave dissipation effect by adjusting parameters such as reef top width and oyster density), which greatly shortens the scheme optimization cycle;

[0079] 3. Ensure computational reliability and security: Strictly perform parameter domain checks (verify whether ds / Hi>0 and B / Hi≤25 are satisfied), automatically output "uncertainty prompt" when exceeding the limit to avoid misuse of extrapolation; stably reproduce the regression formula based on 73 sets of physical experiments to ensure that the results are consistent with engineering physical laws (such as the transmission coefficient Kt decreases as the reef top width increases), providing accurate reference for scheme design;

[0080] 4. Support for engineering decision-making and traceability: Automatically generate calculation records containing input parameters, transmission coefficient Kt value, and parameter domain check conclusions, which facilitates design document archiving and result verification, meeting the requirements of "traceability and verifiability" in engineering management;

[0081] 5. Reduce assessment costs and technical barriers: Without relying on complex numerical simulations or physical model tests, the system can quickly assess and replace traditional high-cost methods, allowing engineers (even those without professional simulation skills) to easily apply ecological revetment wave reduction assessment technology, and helping oyster reef restoration projects to achieve rapid optimization of solutions that synergistically enhance "ecology + disaster prevention".

[0082] The foregoing descriptions and embodiments are provided to enable those skilled in the art to understand and apply the present invention. It will be apparent to those skilled in the art that various modifications can be easily made to these contents, and the general principles described herein can be applied to other embodiments without creative effort. Therefore, the present invention is not limited to the foregoing descriptions and embodiments. Improvements and modifications made by those skilled in the art based on the disclosure of the present invention without departing from its scope should be within the protection scope of the present invention.

Claims

1. A method for rapidly calculating the wave-dissipating capacity of castle-shaped oyster reefs, characterized in that, Includes the following steps: S1. Collect incident wave and reef parameters: including incident wave height Hi, incident wavelength Li, and distance from the reef top to the still water surface. d s reef top width B Oyster attachment density ; S2. Convert the parameters into four dimensionless independent variables: relative submergence of the reef ds / Hi, relative reef crest width B / Hi, incident wave steepness Hi / Li, and oyster density. 0; S3. Perform parameter domain checks on the dimensionless independent variables ds / Hi and B / Hi of the reef body relative submergence. S4. If the parameter domain check passes, substitute the four dimensionless independent variables into the empirical formula for the closed-loop transmission coefficient obtained from multivariate nonlinear regression to obtain the predicted transmission coefficient value Kt; wherein, the empirical formula for the closed-loop transmission coefficient is: ； The formula is derived from regression analysis of 73 physical model experiments, with a goodness of fit R² = 0.768 and a significance level < 0.05, and can stably predict Kt within the parameter domain. S5. Output and record the calculation results, including the set of input parameters, the value of Kt, and the conclusion of the parameter domain check.

2. The method according to claim 1, characterized in that, The condition for parameter domain checking in step S3 is: the relative submergence degree of the reef ds / Hi>0 and the relative width of the reef top B / Hi≤25; if the condition is not met, an "uncertainty prompt" will be output.

3. The method according to claim 1, characterized in that, The oyster attachment density mentioned in step S5 This is an ecological parameter, measured in units per square meter, used to reflect the ecological wave-dissipating characteristics of castle-shaped oyster reefs.

4. The method according to claim 1, characterized in that, The predicted transmission coefficient Kt ranges from 0 to 1, and the smaller Kt is, the better the wave-dissipating effect of the oyster reef.

5. The method according to claim 1, characterized in that, The rapid calculation is a closed-loop substitution operation, and the completion time is independent of the input size. It is suitable for sensitivity scanning and rolling comparison of multiple schemes in the early stage of engineering.

6. A system for performing the method according to any one of claims 1-5, characterized in that, It includes a processor and a memory, the memory storing instructions that, when executed on the processor, implement all the steps of the method.

Images