Method for evaluating dynamic endurance capacity of coral to island reef water

Abstract

The invention relates to the technical field of coral tolerance, and discloses a method for evaluating the hydrodynamic tolerance of corals to island reefs, and the method comprises the steps: preparing a bionic coral model based on the morphological characteristics of a target coral species; different waves are simulated through a large-scale wave water tank, and wave height, flow velocity and stress data of the bionic coral model at different positions of the reef are measured; establishing a critical wave height-flow velocity-stress threshold curve of the damaged target coral species according to the wave height, the flow velocity and the stress data; and evaluating the endurance capacity of the target coral species according to the critical wave height-flow velocity-stress threshold curve of the damaged target coral species. According to the scheme, a research method can be provided for the tolerance of corals to seawater flow caused by waves under the real island terrain, a theoretical basis is provided for coral ecological protection, and a reliable technical means is provided for sustainable management of an island ecological system.

Description

Technical Field

[0001] The invention belongs to the technical field of coral tolerance, and in particular relates to a method for evaluating the tolerance of corals to island and reef hydrodynamics. Background Art

[0002] Coral reef ecosystems are important carriers of marine biodiversity and a key barrier to island and reef ecological protection. Currently, the assessment methods for coral tolerance to hydrodynamic environments mainly rely on field observations, laboratory simplified model experiments, and numerical simulations, but there are still significant deficiencies in technical paths and application effects.

[0003] In the existing technology, the field observation method can reflect the real environmental response by long-term monitoring of the correlation between the physical damage of corals and wave parameters in the natural environment. However, it is difficult to systematically reveal the dynamic coupling mechanism between coral morphological characteristics and hydrodynamic parameters due to the long observation period, uncontrollable environmental variables, and high data collection costs. Laboratory studies often use regular geometric bodies (such as cylinders and cubes) to simulate coral monomers and test their hydrodynamic responses in small water tanks, but such simplified models cannot truly reproduce the complex morphology of coral branching structures and their nonlinear responses to eddy currents and shear stresses, resulting in significant deviations between experimental results and actual coral damage thresholds. Numerical simulation technology can predict the flow field distribution around corals by establishing a fluid-structure interaction model, but its accuracy is highly dependent on boundary condition settings and material parameter assumptions, especially the over-simplification of the mechanical properties of the porous skeleton of corals, making it difficult to accurately quantify the stress accumulation effect of biomaterials under dynamic loads.

[0004] In addition, existing methods generally have the following technical bottlenecks: (1) There is a lack of systematic research on the coupling mechanism between coral morphological characteristics and multiple parameters of the hydrodynamic environment (wave height, flow velocity, stress), resulting in a single evaluation dimension; (2) There is no critical damage threshold curve based on biomechanical properties, which cannot provide a quantitative decision-making basis for coral ecological restoration; (3) The traditional physical model lacks geometric similarity with real corals, and experimental data are difficult to effectively extrapolate to complex island and reef terrain scenarios. These problems have seriously restricted the scientific evaluation of coral tolerance and made the formulation of ecological protection strategies lack accurate theoretical support.

[0005] Therefore, there is an urgent need to develop an assessment method that integrates bionic model construction, high-precision environmental simulation and multi-parameter threshold criteria to break through the limitations of existing technologies and provide reliable technical means for the sustainable management of island and reef ecosystems. Summary of the invention

[0006] In order to solve the above technical problems, the present invention provides a method for evaluating the tolerance of corals to the hydrodynamics of islands and reefs, which can provide a research method for the tolerance of corals to the flow of seawater caused by waves under real island and reef terrain, provide a theoretical basis for coral ecological protection, and provide reliable technical means for the sustainable management of island and reef ecosystems.

[0007] The present invention provides a method for evaluating the tolerance of corals to island and reef hydrodynamics, the method comprising the following steps: S1. Prepare a bionic coral model based on the morphological characteristics of the target coral species; S2. Use a large-scale wave tank to simulate different waves and measure the wave height, flow velocity and stress data of the bionic coral model at different locations on the reef flat; S3. Establishing the critical wave height-flow velocity-stress threshold curve for the destruction of target coral species based on wave height, flow velocity and stress data; S4. Evaluate the tolerance of the target coral species based on the critical wave height-flow velocity-stress threshold curve for damage to the target coral species.

[0008] Furthermore, in S1, preparing a bionic coral model based on the morphological characteristics of the target coral species includes: S11. determining an overall shape index and a branching shape index of the target coral species based on the morphological characteristics of the target coral species; S12, establishing a generalized model based on the overall shape index, branching shape index and target modeling height of the target coral species; S13. Determine the preparation material according to the generalized model, and prepare the bionic coral model to be tested by 3D printing; S14, performing a mechanical property test on the bionic coral model to be tested to obtain a mechanical property test result; S15, comparing the mechanical property test results with the actual values ​​of the mechanical properties of the target coral species; If the error between the mechanical property test result and the true value of the mechanical property is within the preset error range, the bionic coral model to be tested is used as the final bionic coral model; otherwise, the prepared material is adjusted and S13-S15 are executed again.

[0009] Furthermore, in S11, the overall shape index is the ratio of the projection height to the projection width, and the branch shape index is the ratio of the branch length to the branch diameter.

[0010] Further, in S14, a mechanical property test is performed on the bionic coral model to be tested, and the mechanical property test results obtained include: The bone strength of the bionic coral model to be tested was measured to obtain the elastic modulus, compressive strength and critical stress intensity factor; The matrix strength of the bionic coral model to be tested was measured to obtain the matrix bonding strength.

[0011] Furthermore, in S2, different waves were simulated by a large-scale wave tank, and the wave height, flow velocity and stress data of the bionic coral model at different locations on the reef flat were measured, including: S21. Setting the bionic coral model at different locations on the reef flat in a large-scale wave tank to simulate different wave conditions for testing; S22. During the test, the flow velocity distribution at different positions on the reef flat, the wave surface elevation and the stress time history curve of the bionic coral model under the action of waves and currents are recorded in real time. When the bionic coral model is destroyed, the test ends; Among them, the damage to the bionic coral model included matrix removal and branch breakage.

[0012] Further, in S22, when the base shear stress of the bionic coral model is greater than or equal to the matrix bonding strength, it is determined that the matrix is ​​removed; When the stress intensity factor of the bionic coral model is greater than or equal to the critical stress intensity factor, it is judged as a branch fracture.

[0013] Furthermore, in S3, the critical wave height-flow velocity-stress threshold curve for the destruction of the target coral species is established based on the wave height, flow velocity and stress data, including: S31. Determine the morphological parameter function through multiple regression and establish the wave height-velocity threshold equation for the destruction of the target coral species; The equation expression is as follows: H 0 =f(U 0 ,a,b,c); Among them, H 0 represents the critical wave height, U 0 represents the critical flow velocity, a, b, c represent the morphological parameter functions; S32. Draw the tolerance envelope of the target coral species based on the wave height-current threshold equation; S33, performing low-pass filtering on the stress time history curve to obtain the maximum principal stress; S34, calculating the stress coefficient according to the maximum principal stress; S35. Establish a correlation equation between the stress coefficient, wave height and flow velocity, and obtain a critical wave height-flow velocity-stress threshold curve for the destruction of the target coral species.

[0014] Furthermore, in S34, the stress coefficient is calculated according to the maximum principal stress, and the calculation formula is as follows: ; Among them, C σ represents the stress coefficient, σ maxrepresents the maximum principal stress, ρ w represents water density and U represents flow velocity.

[0015] The embodiments of the present invention have the following technical effects: This scheme prepares a bionic coral model by defining the overall shape index and the branch shape index, and ensures that the mechanical response of the model is consistent with that of the real coral through iterative optimization of the mechanical properties; simulates the real reef flat topography in a large-scale wave tank, and simultaneously collects wave height, flow velocity, and stress time-history curves. Through low-pass filtering and defining dimensionless stress coefficients, a three-dimensional correlation equation of wave height-flow velocity-stress is established. By quantifying the dynamic coupling relationship between the combined action of waves and flows and the internal stress response of corals, the coral tolerance mechanism of wave-induced flows under real topography is revealed. This can provide a research method for the tolerance of corals to wave-induced seawater flow under real island and reef topography, provide a theoretical basis for coral ecological protection, and provide reliable technical means for the sustainable management of island and reef ecosystems. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 This is a flow chart of a method for evaluating the hydrodynamic tolerance of corals to islands and reefs provided by an embodiment of the present invention; Figure 2 This is a generalized schematic diagram of a Staghorn coral provided by an embodiment of the present invention; Figure 2 (a) is a schematic diagram of the prototype coral of Acropora, Figure 2 (b) is a generalized coral diagram of Acropora; Figure 3 It is a schematic diagram of the layout of a coral wave and current tolerance test model provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0018] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be described clearly and completely below. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work belong to the scope of protection of the present invention.

[0019] The embodiment of the present invention provides a method for evaluating the tolerance of corals to the hydrodynamics of islands and reefs. Figure 1is a flow chart of a method for evaluating the hydrodynamic tolerance of corals to islands and reefs provided by an embodiment of the present invention, see Figure 1 , the method comprises the following steps: S1. Prepare a bionic coral model based on the morphological characteristics of the target coral species.

[0020] In some embodiments, S1 specifically includes the following sub-steps: S11. Determine an overall shape index and a branching shape index of the target coral species based on the morphological characteristics of the target coral species.

[0021] S12. Establish a generalized model based on the overall shape index, branching shape index and target modeling height of the target coral species.

[0022] In this example, Acropora is used as the target coral species, the morphological data of Acropora is analyzed, and the accurate morphological data of coral is obtained by using 3D scanning technology. The overall shape index I is selected. g (Projection height H c With projection width W c The ratio of ) is 1:3, 1:1 and 3:1, and the branch shape index I b (Branch length L c With branch diameter D c The study was conducted on Acropora with a ratio of 8:1, 5:1, 3:1 and 1:1. The coral model was generalized to a height of 20 cm. Figure 2 is a schematic diagram of a generalized antler coral provided by an embodiment of the present invention. The generalized antler coral morphology is as follows Figure 2 As shown, Figure 2 (a) is a schematic diagram of the prototype coral of Acropora, Figure 2 (b) is a generalized coral diagram of Acropora, Figure 2 The overall morphological index of the generalized Acropora in (b) is 1:1, and the branching indices are 8:1, 5:1, 3:1, and 1:1.

[0023] S13. Determine the preparation materials according to the generalized model, and prepare the bionic coral model to be tested by 3D printing.

[0024] In some embodiments, the prepared material can use a light-curable resin to simulate the coral skeleton, and a coral sand-epoxy resin composite material to simulate the coral matrix.

[0025] S14. Perform a mechanical property test on the bionic coral model to be tested to obtain a mechanical property test result.

[0026] In some embodiments, the bone strength of the bionic coral model to be tested is measured, and the bone strength measurement includes: Uniaxial compression test: The loading rate of the testing machine head is 0.5 mm / min, the stress-strain curve is recorded, and the elastic modulus E and compressive strength σ are calculated. bc .

[0027] Fracture toughness test: Three-point bending method is used to calculate the critical stress intensity factor K IC .

[0028] The substrate strength of the bionic coral model to be tested is tested. The substrate strength test includes: Simulate the intertidal environment (temperature 25-30℃, salinity 30-35‰, pH 7.8-8.2), set up a dry-wet cycle (immersion / exposure every 6 hours), and conduct a shear test after 28 days to measure the substrate bond strength τ max =F / A, where F is the applied shear force and A is the projected area of ​​the matrix.

[0029] S15. Compare the mechanical property test results with the true values ​​of the mechanical properties of the target coral species.

[0030] In some embodiments, the true values ​​of the elastic modulus, compressive strength and matrix bonding strength of the coral are calculated by combining data investigation with previous field test data, and determined by limit destruction test. If the error between the mechanical property test result and the true value of the mechanical property is within the preset error range, the bionic coral model to be tested is used as the final bionic coral model; otherwise, the preparation material is adjusted and S13-S15 is re-executed; wherein the preset error range can be 5%.

[0031] S2. Use a large-scale wave tank to simulate different waves and measure the wave height, flow velocity and stress data of the bionic coral model at different locations on the reef flat.

[0032] S21. The bionic coral model was placed at different locations on the reef flat in a large-scale wave tank to simulate different wave conditions for testing.

[0033] Figure 3 is a schematic diagram of the layout of a coral wave and current tolerance test model provided by an embodiment of the present invention, see Figure 3 First, the island and reef terrain was generalized, using a 1:8 model scale. The model and instrument layout are as follows: Figure 3 The model is 5m wide and 40m long. The slope of the front slope of the island reef is considered as a 1:1 comprehensive slope. The distance between the foot of the steep slope of the island reef and the wave maker is 330m. Regular waves are used for wave simulation. The wave parameters mainly consider the extreme wave design with different recurrence periods. The wave height is 4m-12m, the period is 8-14s, the reef flat water depth is 0-3m, and the distance between the coral position and the reef edge is 0-300m. The main test parameters are shown in Table 1.

[0034] Table 1 Test parameters

[0035] The bionic coral models are arranged at the reef edge, reef crown, outer reef flat and inner reef flat.

[0036] S22. During the test, the velocity distribution at different locations on the reef flat, the wave surface elevation, and the stress time history curve of the bionic coral model under the action of waves and currents are recorded in real time. When the bionic coral model is destroyed, the test ends.

[0037] In some embodiments, an ADV current meter (accuracy ±1%) and a wave height meter are used to obtain the current velocity distribution (U(z)) and wave surface elevation (η(t)) at different locations on the reef flat.

[0038] Micro dynamic strain gauges (sampling frequency f=1000Hz) are embedded in the branches of the bionic coral model to synchronously record the stress time history curve (σ(t)) under the action of wave flow, and a high-definition camera system is used to observe the changes in coral morphology in real time.

[0039] Among them, the damage to the bionic coral model includes matrix removal and branch breakage: When the base shear stress of the bionic coral model is greater than or equal to the matrix bonding strength (τ ≥ τ max ), it was determined to be matrix removal; When the stress intensity factor of the bionic coral model is greater than or equal to the critical stress intensity factor (K I ≥K IC ), it is judged as a branch break.

[0040] S3. Establish a critical wave height-flow velocity-stress threshold curve for damage to the target coral species based on wave height, flow velocity and stress data.

[0041] In some embodiments, S3 includes the following sub-steps: S31. Determine the morphological parameter function through multivariate regression and establish the wave height-velocity threshold equation for the destruction of the target coral species.

[0042] In some embodiments, the equation is expressed as follows: H 0 =f(U 0 ,a,b,c); In some embodiments, f(U 0 ,a,b,c)=a×U 0 b +c; Among them, H 0 represents the critical wave height (i.e., the wave height corresponding to the destruction of the target coral species), U 0represents the critical flow velocity, i.e., the water velocity corresponding to the destruction of the target coral species), a, b, and c represent morphological parameter functions, which are related to the geometric characteristics of the target coral species and reflect the differences in the damage resistance of different coral structures. The experimental data can be substituted into the equation, and the specific expressions of a, b, and c can be fitted through statistical methods (such as the least squares method).

[0043] S32. Draw the tolerance envelope of the target coral species based on the wave height-current velocity threshold equation.

[0044] In some embodiments, the fitted equation is visualized in the HU coordinate system to form a dividing line, where the area below the envelope is a safe area (corals are not damaged) and the area above the envelope is a dangerous area (corals may be damaged).

[0045] S33. Perform low-pass filtering on the stress time history curve to obtain the maximum principal stress.

[0046] In some embodiments, the stress time history curve σ(t) is low-pass filtered (cut-off frequency fc=50 Hz) to extract the maximum principal stress σ max .

[0047] S34. Calculate the stress coefficient based on the maximum principal stress.

[0048] In some embodiments, the calculation formula of the stress coefficient is as follows: ; Among them, C σ represents the stress coefficient, σ max represents the maximum principal stress, ρ w represents water density and U represents flow velocity.

[0049] S35. Establish a correlation equation between the stress coefficient, wave height and flow velocity, and obtain a critical wave height-flow velocity-stress threshold curve for the destruction of the target coral species.

[0050] Since the stresses on corals are caused by the combined effects of waves and currents, their physical nature can be broken down into: Wave force contribution: wave height H is transmitted to the coral structure through inertial force and drag force, and the inertial force is positively correlated with wave height H; Contribution of water flow force: The drag force generated by the flow velocity U is proportional to the square of the flow velocity; Therefore, the internal stress of coral can be regarded as the relationship between H and U 2 The joint function of H / U 2 It comprehensively reflects the wave height (H) and flow velocity (U 2 ) is closer to the actual mechanism of coral stress, so the stress coefficient and H / U 2 The correlation equation of the critical wave height-flow velocity-stress threshold curve for the destruction of the target coral species was obtained.

[0051] S4. Evaluate the tolerance of the target coral species based on the critical wave height-flow velocity-stress threshold curve for damage to the target coral species.

[0052] This scheme prepares a bionic coral model by defining the overall shape index and the branch shape index, and ensures that the mechanical response of the model is consistent with that of the real coral through iterative optimization of the mechanical properties; simulates the real reef flat topography in a large-scale wave tank, and simultaneously collects wave height, flow velocity, and stress time-history curves. Through low-pass filtering and defining dimensionless stress coefficients, a three-dimensional correlation equation of wave height-flow velocity-stress is established. By quantifying the dynamic coupling relationship between the combined action of waves and flows and the internal stress response of corals, the coral tolerance mechanism of wave-induced flows under real topography is revealed. This can provide a research method for the tolerance of corals to wave-induced seawater flow under real island and reef topography, provide a theoretical basis for coral ecological protection, and provide reliable technical means for the sustainable management of island and reef ecosystems.

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

Claims

1. A method for evaluating the tolerance of corals to the hydrodynamic forces of islands and reefs, characterized in that: The method comprises the following steps: S1. Prepare a bionic coral model based on the morphological characteristics of the target coral species; S2. Simulating different waves through a large-scale wave tank to measure the wave height, flow velocity and stress data of the bionic coral model at different locations on the reef flat; S3, establishing a critical wave height-flow velocity-stress threshold curve for the target coral species to be destroyed according to the wave height, flow velocity and stress data; S4. Evaluate the tolerance of the target coral species according to the critical wave height-flow velocity-stress threshold curve for damage to the target coral species.

2. The method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 1, characterized in that: In S1, preparing a bionic coral model based on the morphological characteristics of the target coral species includes: S11. determining an overall shape index and a branching shape index of the target coral species based on the morphological characteristics of the target coral species; S12, establishing a generalized model according to the overall shape index, branch shape index and target modeling height of the target coral species; S13, determining the preparation material according to the generalized model, and preparing the bionic coral model to be tested by 3D printing; S14, performing a mechanical property test on the bionic coral model to be tested to obtain a mechanical property test result; S15, comparing the mechanical property test result with the true value of the mechanical property of the target coral species; If the error between the mechanical property test result and the true value of the mechanical property is within a preset error range, the bionic coral model to be tested is used as the final bionic coral model; otherwise, the prepared material is adjusted and S13-S15 are re-executed.

3. A method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 2, characterized in that: In S11, the overall shape index is the ratio of the projection height to the projection width, and the branch shape index is the ratio of the branch length to the branch diameter.

4. The method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 2, characterized in that: In S14, the mechanical properties of the bionic coral model to be tested are tested, and the mechanical properties test results obtained include: Performing bone strength measurement on the bionic coral model to be tested to obtain elastic modulus, compressive strength and critical stress intensity factor; The matrix strength of the bionic coral model to be tested is measured to obtain the matrix bonding strength.

5. A method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 4, characterized in that: In S2, different waves are simulated by a large-scale wave tank, and wave height, flow velocity and stress data of the bionic coral model at different locations on the reef flat are measured, including: S21, placing the bionic coral model at different locations on the reef flat in a large-scale wave tank, simulating different wave conditions for testing; S22. During the test, the flow velocity distribution at different positions on the reef flat, the wave surface elevation, and the stress time history curve of the bionic coral model under the action of waves and currents are recorded in real time. When the bionic coral model is destroyed, the test ends; The destruction of the bionic coral model includes matrix removal and branch breakage.

6. A method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 5, characterized in that: In S22, when the base shear stress of the bionic coral model is greater than or equal to the matrix bonding strength, it is determined that the matrix is ​​removed; When the stress intensity factor of the bionic coral model is greater than or equal to the critical stress intensity factor, it is determined to be a branch fracture.

7. A method for evaluating the ability of corals to tolerate hydrodynamic forces on islands and reefs according to claim 5, characterized in that: In S3, establishing a critical wave height-flow velocity-stress threshold curve for damage to the target coral species based on the wave height, flow velocity and stress data includes: S31, determining the morphological parameter function through multiple regression, and establishing a wave height-flow velocity threshold equation for the destruction of the target coral species; The equation expression is as follows: H0=f(U0,a,b,c); Among them, H0 represents the critical wave height, U0 represents the critical flow velocity, a, b, and c represent the morphological parameter functions; S32, drawing a tolerance envelope of the target coral species according to a wave height-current velocity threshold equation; S33, performing low-pass filtering on the stress time history curve to obtain the maximum principal stress; S34, calculating a stress coefficient according to the maximum principal stress; S35, establishing a correlation equation between the stress coefficient and the wave height and flow velocity, and obtaining a critical wave height-flow velocity-stress threshold curve for the target coral species to be destroyed.

8. A method for evaluating the hydrodynamic tolerance of corals to islands and reefs according to claim 7, characterized in that: In S34, the stress coefficient is calculated according to the maximum principal stress, and the calculation formula is as follows: ； Among them, C σ represents the stress coefficient, σ max represents the maximum principal stress, ρ w represents water density and U represents flow velocity.

Images