A method, device, computer equipment and readable storage medium for shaping estuary and coastal ecological shoals by coordinating water and sand regulation with habitat creation

Abstract

The present invention discloses a method, apparatus, computer device, and readable storage medium for shaping estuarine and coastal ecological shoals that coordinate water and sediment regulation with habitat creation. The method includes: first, determining the shaping area, selecting a reference shoal, and obtaining relevant data; dividing the spatial units and interpolating them to construct a correlation model between shoal wetland elevation, sediment particle size, and ecological characteristics; generating a current terrain elevation field based on a randomized algorithm, and combining sediment stability and skeleton rule correction to obtain a design elevation field; obtaining erosion and deposition balance data through tidal and sediment numerical simulation, and inputting this data into the model to obtain candidate ecological characteristic data; and repeating the process to select the optimal design elevation field as the target habitat for shaping. This method integrates multiple factors to achieve coordinated water and sediment regulation and habitat creation, promoting the restoration and improvement of estuarine and coastal ecosystems.

Description

Technical Field

[0001] The present invention relates to the field of ecological management, and in particular to a method, device, computer equipment and readable storage medium for shaping estuary and coastal ecological shoals in a coordinated manner of water and sand regulation and habitat creation. Background Art

[0002] Estuarine and coastal ecosystems are crucial for maintaining biodiversity and protecting coastlines. However, human activities and changes in the natural environment have led to the degradation of the ecological functions of estuarine and coastal wetlands. Traditional ecological shoal formation methods often focus solely on topography or a single ecological factor, lacking the synergy between water and sediment regulation and habitat creation, making it difficult to effectively restore and enhance the functions of estuarine and coastal ecosystems. Summary of the Invention

[0003] The purpose of the present invention is to provide a method, device, computer equipment and readable storage medium for shaping estuary and coastal ecological shoals that coordinate water and sand regulation and habitat creation.

[0004] In a first aspect, an embodiment of the present invention provides a method for shaping an estuary and coastal ecological shoal by coordinating water and sediment regulation with habitat creation, comprising:

[0005] Determine the spatial scope of the area where ecological shoal shaping is to be carried out, wherein the area where ecological shoal shaping is to be carried out includes a plurality of candidate shoals;

[0006] Selecting at least two reference tidal flats from the plurality of candidate tidal flats, and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats;

[0007] Dividing the at least two reference tidal flats into a plurality of spatial units, and interpolating each of the spatial units based on the tidal flat wetland elevation information, the sediment median particle size, and the ecological characteristic data, to obtain the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units;

[0008] Constructing a target structural equation model representing a first correlation relationship included in each of the spatial units, wherein the first correlation relationship is a correlation relationship between reference tidal flat wetland elevation information, target sediment median particle size, and target ecological characteristic data corresponding to each of the spatial units;

[0009] Constructing a second association relationship for each of the spatial units based on a statistical regression analysis strategy, wherein the second association relationship is an association relationship between the reference tidal flat wetland elevation information and the target sediment median particle size corresponding to each of the spatial units;

[0010] Determining a current terrain elevation field of the area to be subjected to ecological shoal shaping based on a random generation algorithm, and modifying the current terrain elevation field using the second association relationship and a preset rule to obtain a design elevation field of the area to be subjected to ecological shoal shaping;

[0011] Conducting numerical simulation of tidal current and sediment on the designed elevation field to obtain the scouring and deposition equilibrium elevation field and the median particle size of scouring and deposition equilibrium sediment in the area to be developed for ecological shoal shaping under the regional scouring and deposition equilibrium state;

[0012] Inputting the scour-sedimentation balance elevation field and the scour-sedimentation balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data;

[0013] Repeating the step of determining the current terrain elevation field of the area to be ecological shoal shaping based on the random generation algorithm for a preset number of times, until the step of inputting the erosion-deposition balance elevation field and the erosion-deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, to obtain a plurality of candidate ecological characteristic data;

[0014] A target design elevation field corresponding to the target ecological characteristic data is selected from the plurality of candidate ecological characteristic data according to preset ecological indicators to be used as the target habitat for the ecological shallows shaping area to be developed.

[0015] In a second aspect, an embodiment of the present invention provides an estuary and coastal ecological shoal shaping device that coordinates water and sand regulation with habitat creation, comprising:

[0016] An acquisition module is used to determine the spatial scope of the ecological shoal shaping area to be carried out, and the surrounding area of ​​the ecological shoal shaping area to be carried out includes multiple candidate shoals; at least two reference shoals are selected from the multiple candidate shoals, and the tidal flat wetland elevation information, sediment median particle size and ecological characteristic data of the at least two reference shoals are obtained; the at least two reference shoals are divided into multiple spatial units, and each of the spatial units is interpolated based on the tidal flat wetland elevation information, the sediment median particle size and the ecological characteristic data to obtain the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each spatial unit; a target structural equation model is constructed to characterize the first correlation relationship included in each of the spatial units, and the first correlation relationship is the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each of the spatial units. The correlation between particle size and target ecological characteristic data; constructing a second correlation relationship included in each of the spatial units based on a statistical regression analysis strategy, wherein the second correlation relationship is the correlation relationship between the reference tidal flat wetland elevation information corresponding to each of the spatial units and the target sediment median particle size; determining the current terrain elevation field of the ecological shoal shaping area to be carried out based on a random generation algorithm, and correcting the current terrain elevation field using the second correlation relationship and preset rules to obtain the design elevation field of the ecological shoal shaping area to be carried out; performing tidal sediment numerical simulation on the design elevation field to obtain the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment in the ecological shoal shaping area to be carried out under the regional scouring and sedimentation balance state; inputting the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment into the target structural equation model to obtain candidate ecological characteristic data;

[0017] A shaping module is used to repeatedly execute the step of determining the current terrain elevation field of the ecological shallows shaping area to be carried out based on the random generation algorithm for a preset number of times, to the step of inputting the erosion and deposition balance elevation field and the erosion and deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, and obtain multiple candidate ecological characteristic data; from the multiple candidate ecological characteristic data, a target design elevation field corresponding to the target ecological characteristic data is screened out according to preset ecological indicators as the target habitat of the ecological shallows shaping area to be carried out for shaping.

[0018] An acquisition module is used to determine the spatial scope of the ecological shoal shaping area to be carried out, and the surrounding area of ​​the ecological shoal shaping area to be carried out includes multiple candidate shoals; at least two reference shoals are selected from the multiple candidate shoals, and the tidal flat wetland elevation information, sediment median particle size and ecological characteristic data of the at least two reference shoals are obtained; the at least two reference shoals are divided into multiple spatial units, and each of the spatial units is interpolated based on the tidal flat wetland elevation information, the sediment median particle size and the ecological characteristic data to obtain the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each spatial unit; a target structural equation model is constructed to characterize the first correlation relationship included in each of the spatial units, and the first correlation relationship is the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each of the spatial units. The correlation between particle size and target ecological characteristic data; constructing a second correlation relationship included in each of the spatial units based on a statistical regression analysis strategy, wherein the second correlation relationship is the correlation relationship between the reference tidal flat wetland elevation information corresponding to each of the spatial units and the target sediment median particle size; determining the current terrain elevation field of the ecological shoal shaping area to be carried out based on a random generation algorithm, and correcting the current terrain elevation field using the second correlation relationship and preset rules to obtain the design elevation field of the ecological shoal shaping area to be carried out; performing tidal sediment numerical simulation on the design elevation field to obtain the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment in the ecological shoal shaping area to be carried out under the regional scouring and sedimentation balance state; inputting the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment into the target structural equation model to obtain candidate ecological characteristic data;

[0019] A shaping module is used to repeatedly execute the step of determining the current terrain elevation field of the ecological shallows shaping area to be carried out based on the random generation algorithm for a preset number of times, to the step of inputting the erosion and deposition balance elevation field and the erosion and deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, and obtain multiple candidate ecological characteristic data; from the multiple candidate ecological characteristic data, a target design elevation field corresponding to the target ecological characteristic data is screened out according to preset ecological indicators as the target habitat of the ecological shallows shaping area to be carried out for shaping.

[0020] In a third aspect, an embodiment of the present invention provides a computer device, comprising a processor and a non-volatile memory storing computer instructions, wherein when the computer instructions are executed by the processor, the computer device executes the method described in the first aspect.

[0021] In a fourth aspect, an embodiment of the present invention provides a readable storage medium, characterized in that the readable storage medium includes a computer program, and when the computer program is running, it controls the computer device where the readable storage medium is located to execute the method described in the first aspect.

[0022] Compared to existing technologies, the present invention offers the following advantages: A method, device, computer equipment, and readable storage medium for shaping estuarine and coastal ecological shoals that coordinate water and sediment regulation with habitat creation, disclosed herein, include: first, determining the scope of the shaping area, selecting a reference shoal, and acquiring relevant data; dividing spatial units and interpolating to construct a correlation model between shoal wetland elevation, sediment particle size, and ecological characteristics; generating a current terrain elevation field based on a randomized algorithm, and combining sediment stability and skeleton rule correction to obtain a design elevation field; obtaining scouring and deposition balance data through tidal and sediment numerical simulation, and inputting this into the model to obtain candidate ecological characteristic data; and repeating the process to select the optimal design elevation field as the target habitat for shaping. This method integrates multiple factors to achieve coordinated water and sediment regulation and habitat creation, promoting the restoration and improvement of estuarine and coastal ecosystems. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly describes the drawings required for use in the embodiments. It should be understood that the following drawings illustrate only certain embodiments of the present invention and should not be construed as limiting the scope of the present invention. Those skilled in the art can, without inventive effort, derive other relevant drawings from these drawings.

[0024] Figure 1 A schematic diagram of the steps of a method for shaping an estuary and coastal ecological shoal by coordinating water and sediment regulation and habitat creation according to an embodiment of the present invention;

[0025] Figure 2 Schematic diagram of the technical framework for shaping estuarine and coastal ecological shoals by coordinating water and sediment regulation and habitat creation according to an embodiment of the present invention;

[0026] Figure 3 A schematic diagram of a structural equation model analysis framework provided by an embodiment of the present invention;

[0027] Figure 4 A schematic diagram of an area where ecological shoal shaping is to be carried out according to an embodiment of the present invention;

[0028] Figure 5 A schematic block diagram of the structure of an estuary and coastal ecological shoal shaping device for coordinated water and sediment regulation and habitat creation provided by an embodiment of the present invention;

[0029] Figure 6 A schematic block diagram of the structure of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings of the embodiments of the present invention. It should be understood that the described embodiments are only a portion of the embodiments of the present invention, not all of them. Generally, the components of the embodiments of the present invention described and illustrated in the drawings herein may be arranged and designed in a variety of different configurations.

[0031] The specific embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0032] In order to solve the technical problems in the above background technology, Figure 1 This is a flow chart of a method for shaping an estuary and coastal ecological shoal by coordinating water and sand regulation with habitat creation provided in an embodiment of the present disclosure. The following is a detailed introduction to the method for shaping an estuary and coastal ecological shoal by coordinating water and sand regulation with habitat creation.

[0033] Step S201, determining the spatial scope of an area to be developed for ecological shoal shaping, wherein the area to be developed for ecological shoal shaping includes a plurality of candidate shoals;

[0034] Step S202: selecting at least two reference tidal flats from the plurality of candidate tidal flats, and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats;

[0035] Step S203: Divide the at least two reference tidal flats into a plurality of spatial units, and interpolate each of the spatial units based on the tidal flat wetland elevation information, the sediment median particle size, and the ecological characteristic data to obtain the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units;

[0036] Step S204: constructing a target structural equation model representing a first correlation relationship included in each of the spatial units, wherein the first correlation relationship is a correlation relationship between the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units;

[0037] Step S205: constructing a second association relationship for each of the spatial units based on a statistical regression analysis strategy, wherein the second association relationship is an association relationship between the reference tidal flat wetland elevation information and the target sediment median particle size corresponding to each of the spatial units;

[0038] Step S206: determining the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm, and modifying the current terrain elevation field using the second association relationship and a preset rule to obtain a design elevation field of the area to be ecological shoal shaping;

[0039] Step S207, performing a tidal sediment numerical simulation on the designed elevation field to obtain an erosion-deposition equilibrium elevation field and an erosion-deposition equilibrium sediment median particle size in the area to be developed for ecological shoal shaping under a regional erosion-deposition equilibrium state;

[0040] Step S208, inputting the erosion-deposition balance elevation field and the erosion-deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data;

[0041] Step S209, repeatedly executing the step of determining the current terrain elevation field of the area to be ecological shoal shaping based on the random generation algorithm for a preset number of times, until the step of inputting the erosion-deposition balance elevation field and the erosion-deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, thereby obtaining a plurality of candidate ecological characteristic data;

[0042] Step S210 , selecting a target design elevation field corresponding to target ecological characteristic data from the plurality of candidate ecological characteristic data according to preset ecological indicators as the target habitat for the ecological shoal shaping area to be carried out.

[0043] In an embodiment of the present invention, for example, the server is connected to the geographic information database of the estuary and coastal area and the relevant monitoring data platform. By collecting and analyzing data on the developed and developing mudflats in the estuary and coastal areas in the area, combined with regional hydrological sediment data, riverbed evolution simulation results, and national land space planning and other information. For example, in a certain estuary and coastal area, there are multiple developing mudflats and some mudflats that have been formed but whose ecological functions need to be improved. Based on these data, the server determines a larger area as the ecological shoal shaping area to be carried out. The area is defined by horizontal and vertical coordinates on the plane. For example, the horizontal coordinate range is from arrive , the vertical axis range is from arrive , in the vertical direction by arrive Determine the elevation range. Within this area, there are multiple candidate tidal flats with different topography, sediment conditions and preliminary ecological conditions. The server stores this spatial range information in the form of a three-dimensional terrain field array. It is stored in the form of , providing a basic framework for subsequent work.

[0044] The server then evaluated multiple candidate tidal flats surrounding the area identified in the first step. By analyzing previously collected satellite remote sensing imagery (such as GF-1 satellite data and QuickBird imagery) and LiDAR data, it made a preliminary assessment of the morphological stability of each candidate. It also consulted existing biological monitoring data to understand the ecological status of each tidal flat, such as the coverage of wetland plants and the number of benthic species. For example, from the numerous candidate tidal flats, the server selected two reference tidal flats, labeled Tidal Flat A and Tidal Flat B.

[0045] For Tidal Flat A, the server reads its tidal flat wetland elevation information from the elevation monitoring data file. These data are obtained through long-term GPS elevation measurement and LiDAR elevation scanning, and accurately record the elevation values ​​of different locations on the tidal flat. In terms of the median particle size of sediment, the server calls the analysis data of the surface sediment samples of Tidal Flat A. These samples were collected from different areas of the tidal flat, and the median particle size of the sediment was measured in the laboratory through screening and other methods. In terms of ecological characteristic data, the server obtains information on the species of wetland plants on the tidal flat, such as the presence of plants such as reeds and alkali sedge, as well as their distribution range; in terms of benthic organisms, there are species, numbers and distribution of organisms such as sandworms and clams; there are also activity records and species information of swimming organisms and birds on the tidal flat.

[0046] Similarly, for tidal flat B, the server also obtains the corresponding tidal flat wetland elevation information, sediment median particle size, and ecological characteristics data in a similar manner. This data is organized and stored in a dedicated database table for subsequent processing.

[0047] The server divides the selected mudflats A and B into spatial units. Assume that each mudflat is divided into 100m x 100m rectangular spatial units. For mudflat A, the server uses an interpolation algorithm (such as Kriging) from the stored mudflat wetland elevation information data to interpolate the discrete elevation measurement point data onto each spatial unit based on the center coordinate position of each spatial unit, obtaining the reference mudflat wetland elevation information corresponding to each spatial unit.

[0048] When processing sediment median particle size data, the server also distributes the sediment median particle size data obtained from different sampling points to each spatial unit through interpolation method based on the location of the spatial unit to obtain the target sediment median particle size.

[0049] For ecological characteristic data, taking wetland plants as an example, the server interpolates the actual distribution and density of wetland plants to calculate target ecological characteristic data such as the species and coverage area within each spatial unit. For benthic organisms, the server interpolates the species and abundance of benthic organisms within each spatial unit based on sampling data from different areas of the mudflat. Through this processing, the server assigns complete reference mudflat elevation information, target sediment median particle size, and target ecological characteristic data to each spatial unit of mudflat A.

[0050] For mudflat B, the server uses the same method to perform spatial unit division and data interpolation operations, ultimately obtaining the relevant data corresponding to each spatial unit of mudflat B. The server further organizes this processed data into a format that is easy to analyze and stores it in the database.

[0051] The server begins to construct a target structural equation model to analyze the correlation between various factors within each spatial unit. Taking a spatial unit of mudflat A as an example, the server uses the reference mudflat wetland elevation information, target sediment median particle size, and target ecological characteristic data corresponding to this spatial unit as variables. First, the server sets a hypothetical path, for example, assuming that terrain elevation affects the growth of wetland plants, thereby affecting the species and number of benthic organisms, and that the median sediment particle size may also affect the distribution of swimming organisms.

[0052] The server then processes the data from these numerous spatial units using statistical software (such as the lavaan package in R). During the model fitting process, the server continuously adjusts the relationship parameters between the variables to ensure that the model accurately reflects the associations in the actual data. By testing the significance of each path, the server retains the significant ones and removes the insignificant hypothetical paths. Ultimately, a target structural equation model is developed that effectively characterizes the relationships between the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data within that spatial unit.

[0053] The server also follows the same process to construct target structural equation models for each spatial unit of mudflat B. These models are stored by the server for subsequent analysis and prediction.

[0054] The server uses a statistical regression analysis strategy for each spatial unit: the reference tidal flat elevation and the target median sediment particle size. For example, for a spatial unit on tidal flat A, the server uses the reference tidal flat elevation as the independent variable and the target median sediment particle size as the dependent variable.

[0055] The server uses regression analysis methods such as the least squares method to fit these two sets of data for a large number of spatial units. By continuously adjusting the parameters of the regression equation, a function expression that can most accurately describe the relationship between the reference tidal flat wetland elevation information and the target sediment median particle size is found. ,in, is the median particle size of sediment, It is the elevation of tidal flat wetland.

[0056] Similarly, the server constructed corresponding secondary correlation functions for each spatial unit of mudflat B. These functions provided important quantitative evidence for determining sediment particle size based on elevation during subsequent landform shaping, or adjusting elevation based on sediment particle size requirements. They were recorded by the server in a dedicated model library.

[0057] The server first determines the current terrain elevation field of the ecological shallows shaping area based on a random generation algorithm (such as Perlin noise, Simplex noise, etc.). Taking a local area in the area as an example, the server randomly generates a series of elevation values ​​according to the set algorithm parameters to form a preliminary random elevation field. ,in It is a randomly generated elevation value within a certain range.

[0058] Then, the server uses the second association relationship constructed above, that is, the relationship function between the reference tidal flat wetland elevation information and the target sediment median particle size. , the server generates a random elevation value based on the unit , calculate the corresponding sediment median particle size through the second correlation function .

[0059] Then, the server calculates the sediment repose angle at the spatial unit according to the sediment repose angle formula (such as Zhang Hongwu formula, Cheng Niansheng formula, etc.) . Then according to the length of the space unit Calculate the maximum elevation difference between the spatial unit and the four adjacent spatial units .

[0060] The server checks whether the elevation difference between the spatial unit and the surrounding adjacent units satisfies the preset relationship (that is, the elevation difference between adjacent units cannot exceed If it is not satisfied, the server adjusts the elevation of the spatial unit to meet the requirements. By performing such checks and adjustments on all spatial units of the entire random elevation field, the server obtains a terrain elevation field corrected based on sediment stability.

[0061] In addition, the server also considers preset rules, such as setting some main and supporting frames (such as concrete or stone embankments) in the area. Based on these rules, the server further adjusts the modified terrain elevation field and finally obtains the design elevation field of the area to be developed for ecological shoal shaping. .

[0062] The server calls a special tidal sediment numerical simulation software (such as Delft3D, etc.) to convert the previously obtained design elevation field The initial terrain conditions are input into the simulation system. At the same time, the server also imports hydrological and sediment data such as water velocity, flow rate, and sediment content of the estuary and coastal area as boundary conditions.

[0063] During the simulation, the software calculates the flow of water over the designed elevation terrain, as well as the transport and deposition of sediment, based on principles of hydrodynamics and sediment kinematics. As the simulation progresses, the system gradually reaches a regional equilibrium between erosion and deposition.

[0064] The server extracts the erosion and deposition balance elevation field of the ecological shoal shaping area under the erosion and deposition balance state from the simulation results. The elevation field reflects the topographic changes after long-term water flow and sedimentation. At the same time, the server also obtains the median particle size of the sediment balance corresponding to each spatial unit. ,These data provide important topographic and sediment condition information for ,subsequent evaluation of ecological characteristics.

[0065] The server will get the erosion and siltation balance elevation field from the previous step. Median particle size of sediment in equilibrium with erosion and deposition This is then input as an independent variable into the previously constructed target structural equation model. Taking the target structural equation model corresponding to mudflat A as an example, the server calculates candidate ecological characteristic data, such as the species, abundance, and biomass of wetland plants, benthic organisms, nematodes, and birds in the area under the current erosion and deposition balance conditions, based on the relationships between the various variables set in the model. The Shannon-Wiener index is also used.

[0066] The server performs the same operation on the target structural equation model corresponding to mudflat B, obtaining the corresponding candidate ecological characteristic data. These candidate ecological characteristic data reflect the ecological status after erosion and deposition balance at the current design elevation field, providing a basis for subsequent screening of the optimal solution.

[0067] The server repeats all the steps from determining the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm to inputting the erosion-deposition equilibrium elevation field and the erosion-deposition equilibrium sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data a preset number of times (assuming 10 times).

[0068] With each iteration, the server generates a different random elevation field. After a series of corrections, simulations, and calculations, it obtains a different erosion-deposition equilibrium elevation field, median sediment particle size, and corresponding candidate ecological signature data. This ultimately yields multiple (here, 10) candidate ecological signature data sets, each corresponding to a different terrain design and ecological condition under erosion-deposition equilibrium conditions.

[0069] The server evaluates and screens the multiple candidate ecological characteristic data obtained previously based on preset ecological indicators (such as species richness, biomass size, Shannon-Wiener index, etc.).

[0070] The server compares the values ​​of various ecological indicators in each candidate ecological feature data set. For example, for species richness, the server calculates the total number of wetland plants, benthic organisms, nematodes, and bird species in each set; for biomass, the server calculates the total biomass of each type of organism; and for the Shannon-Wiener index, the server calculates the corresponding index value based on the abundance of each species.

[0071] By comprehensively comparing these ecological indicators, the server identifies the candidate ecological characteristic data set that performs best across all metrics and determines the corresponding design elevation field as the target design elevation field. This target design elevation field will serve as the target habitat for the ecological shoal shaping area to be implemented, aiming to achieve the best ecological results and promote ecological protection and restoration in estuary and coastal areas.

[0072] In an embodiment of the present invention, the step of selecting at least two reference tidal flats from the multiple candidate tidal flats and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats may be implemented through the following examples.

[0073] selecting at least two reference tidal flats from the plurality of candidate tidal flats according to preset tidal flat morphological stability indicators and good habitat indicators;

[0074] Acquiring remote sensing image data, DEM elevation data, wetland plant data, benthic organism data, swimming organism data, bird data, mudflat elevation data, and bed sand gradation data corresponding to the at least two reference mudflats;

[0075] Using ArcGIS to process the DEM elevation data to obtain the elevation information of the tidal flat wetland;

[0076] Determining the median particle size of the sediment based on the remote sensing image data, the tidal flat elevation data, and the bed sand gradation data;

[0077] The remote sensing image data is interpreted by a support vector machine combined with a deep learning strategy, and combined with the wetland plant data, benthic organism data, swimming organism data and bird data, the species, abundance and biomass of wetland plants, benthic organisms, swimming organisms and birds are obtained as the ecological characteristic data.

[0078] In the embodiment of the present invention, for example, in an ecological shoal shaping project in a certain estuary and coastal area, the server undertakes key data processing and analysis tasks.

[0079] The server first connects to a database that stores geographic information and ecological data about the estuary and coastal areas. The database records relevant information about numerous candidate mudflats. The server screens these candidate mudflats based on preset mudflat morphological stability indicators and good habitat indicators. For example, the mudflat morphological stability indicator is set as the mudflat boundary changing within a certain range in the past three years, and the good habitat indicator includes wetland plant coverage reaching a certain proportion and the benthic biodiversity index being higher than a certain threshold. By analyzing and comparing the data of each candidate mudflat, the server selected two reference mudflats from the many candidate mudflats, named Mudflat M and Mudflat N.

[0080] Next, the server begins acquiring data related to the two reference mudflats. It retrieves remote sensing imagery corresponding to mudflats M and N from the satellite data receiving system, including GF-1 satellite data and QuickBird imagery. It also obtains DEM elevation data from the elevation data storage module. It reads wetland plant data, benthic data, nematode data, and bird data from the biological monitoring data folder. It extracts mudflat elevation data from the topographic survey data file. Finally, it obtains bed sand gradation data from the sediment analysis report.

[0081] After acquiring the data, the server processed the DEM elevation data using ArcGIS software. For example, for Tidal Flat M, the server imported the DEM elevation data into ArcGIS and applied terrain analysis tools to the data, interpolating and smoothing it. This yielded accurate elevation information reflecting the elevation changes of Tidal Flat M, including detailed information such as altitude and slope at different locations. Similarly, the DEM elevation data for Tidal Flat N was processed in the same manner.

[0082] The server then determines the median sediment particle size based on remote sensing imagery, mudflat elevation data, and bed-sand gradation data. The server first analyzes the tonal and texture characteristics of different regions in the remote sensing imagery and, combined with the mudflat elevation data, determines the sediment deposition at various locations. Then, using a specific algorithm and empirical formula, the server references the bed-sand gradation data to calculate the median sediment particle size at each location on mudflats M and N.

[0083] Finally, the server interprets the remote sensing image data using a support vector machine combined with a deep learning strategy. The server inputs the remote sensing image data into the trained support vector machine and deep learning model, which then identifies and classifies the features in the image. Simultaneously, combining data on wetland plants, benthic organisms, nematodes, and birds, the server accurately determines the species, abundance, and biomass of these plants, benthic organisms, nematodes, and birds. For example, on mudflat M, wetland plants such as reeds and Suaeda salsa were identified. Their abundance was determined by calculating the number of plants per unit area, while their biomass was estimated based on the average weight and number of plants. A similar approach was used for benthic organisms, nematodes, and birds, ultimately resulting in complete ecological characteristic data. Through these steps, the server provides comprehensive and accurate data support for subsequent ecological shoal shaping work.

[0084] In an embodiment of the present invention, the current terrain elevation field of the area to be ecological shoal shaping is determined based on a random generation algorithm, and the current terrain elevation field is corrected using the second association relationship and preset rules to obtain the design elevation field of the area to be ecological shoal shaping, which can be implemented through the following examples.

[0085] Determine the current terrain elevation field of the area to be subjected to ecological shoal shaping based on a random generation algorithm;

[0086] Calculating the sediment repose angle of the spatial unit corresponding to the current terrain elevation field based on the second association relationship, and correcting the current terrain elevation field based on the sediment repose angle to obtain a corrected terrain elevation field;

[0087] The modified terrain elevation field is modified based on preset rules for the skeleton of the ecological shoal to be shaped, so as to obtain the final designed elevation field.

[0088] In the embodiment of the present invention, illustratively, in an estuary and coast ecological shoal shaping project, the server begins to determine a design elevation field of an area where ecological shoal shaping is to be carried out.

[0089] The server first determines the current terrain elevation field for the area where the ecological shoal shaping project is to be carried out, using a random generation algorithm, such as the Perlin noise algorithm. Taking a 500m x 500m sub-area within this area as an example, the server randomly generates a series of elevation values ​​within the sub-area according to the set algorithm parameters. These elevation values ​​fluctuate within a certain range, forming a preliminary random elevation field. The server stores this random elevation field as a three-dimensional array, where each element corresponds to the elevation value of a spatial unit (for example, 100m x 100m) within the sub-area.

[0090] Next, the server uses the second association relationship constructed previously, that is, the relationship between the reference tidal flat wetland elevation information and the target sediment median particle size, to calculate the sediment repose angle of the corresponding spatial unit in the current terrain elevation field. Taking a spatial unit U1 in the sub-area as an example, the server calculates the corresponding sediment median particle size through the second association relationship based on the elevation value of the U1 unit in the current terrain elevation field. Assuming that the calculated sediment median particle size is 0.2mm, the server then calculates the sediment repose angle formula (such as Zhang Hongwu's formula: β=35.3×( )^0.04, where (where is the median sediment particle size) and the sediment repose angle for unit U1 is calculated to be approximately 37.6 degrees. The server performs this calculation for all spatial units within the sub-area.

[0091] Based on the calculated sediment repose angle, the server corrects the current terrain elevation field. The server checks the elevation difference between each spatial unit and the four adjacent spatial units. For spatial unit U1, the server calculates the elevation difference between it and the surrounding units U2, U3, U4, and U5. Assume that the elevation of U1 is , the elevation of U2 is ,like Exceeds the maximum allowable elevation difference calculated based on the sediment repose angle (calculated by the formula , L is the length of the spatial unit (here, 100m). The server then adjusts the elevation of U1 or U2 to meet the stability requirements. By checking and adjusting all spatial units in the entire sub-region, the server obtains the corrected terrain elevation field.

[0092] Finally, the server further modifies the revised terrain elevation field based on the preset rules for the skeleton of the ecological shoal to be shaped. The project presets two main skeletons and three branch skeletons in this sub-area. Both the main skeleton and the branch skeleton are concrete embankments. The server adjusts the terrain elevation field at the location of the skeleton according to the preset position and elevation rules. For example, the top elevation of the main skeleton is preset to be 2 meters above the average high tide level, and the server adjusts the terrain elevation at the location of the main skeleton to the corresponding height. After adjusting the positions of all skeletons, the server obtained the final design elevation field. This design elevation field comprehensively considers the stability of the terrain and the engineering rules for the shaping of ecological shoals, and will provide an accurate terrain basis for subsequent tidal sediment numerical simulation and ecological assessment.

[0093] In an embodiment of the present invention, the determination of the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm can be implemented through the following examples.

[0094] Obtain the initial terrain elevation field of the area to be used for ecological shoal shaping, and process the initial terrain elevation field using a preset random noise algorithm to obtain a random terrain elevation field; by formula: , calculate the current terrain elevation field, where is the current terrain elevation field, is the initial terrain elevation field, is the random terrain elevation field, is a vector terrain field array, is the random terrain elevation adjustment coefficient.

[0095] In an embodiment of the present invention, illustratively, in an estuary and coastal ecological shoal shaping project, the server begins to execute a task of determining the current terrain elevation field of the area where the ecological shoal shaping is to be carried out.

[0096] The server first obtains the initial terrain elevation field of the ecological shallows shaping area to be carried out from the project's terrain data repository. Assuming that the area is a 1000m×1000m tidal flat area located at the mouth of a river, the initial terrain elevation field is represented by a three-dimensional vector terrain field array. The form of storage, where The value range of corresponds to the plane coordinate range of the mudflat area. Indicates the initial terrain elevation value.

[0097] Next, the server uses a preset random noise algorithm, such as the Simplex noise algorithm, to process the initial terrain elevation field. The server inputs the data of the initial terrain elevation field into the algorithm program, and the algorithm generates a series of random elevation disturbance values ​​in the area to form a random terrain elevation field. . Take a coordinate point in the area For example, in the initial terrain elevation field, the corresponding elevation value is After being processed by the random noise algorithm, an elevation value with random fluctuations within a certain range is generated. This value becomes the elevation value corresponding to that point in the random terrain elevation field. The server performs this process on all coordinate points in the entire tidal flat area to obtain a complete random terrain elevation field.

[0098] The server then passes the formula To calculate the current terrain elevation field. Here, the server sets the random terrain elevation adjustment coefficient α to 0.5 based on the actual needs of the project and the results of the previous test. , its elevation value in the initial terrain elevation field Elevation values ​​in random terrain elevation fields , calculated according to the formula: The server performs this calculation on each coordinate point in the mudflat area, thereby obtaining the current terrain elevation field of the entire area. .

[0099] Through these steps, the server successfully determined the current terrain elevation field for the area where ecological shoal formation will be carried out. This current terrain elevation field retains the basic characteristics of the initial terrain while introducing random terrain variations. This provides a basis for subsequent terrain elevation field corrections based on sediment stability and preset rules. This helps simulate more realistic and diverse terrain conditions to meet the topographic requirements of ecological shoal formation.

[0100] In an embodiment of the present invention, the calculating of the sediment repose angle of the spatial unit corresponding to the current terrain elevation field based on the second association relationship, and correcting the current terrain elevation field based on the sediment repose angle to obtain a corrected terrain elevation field, includes:

[0101] The second association relationship is obtained, where the second association relationship is: ,in, is the median particle size of sediment, is the elevation of the tidal flat wetland;

[0102] According to the second association relationship, through the formula: Calculate and obtain the median particle size of sediment in the spatial unit corresponding to the current terrain elevation field; wherein, is the coordinate position of the current spatial unit corresponding to the current terrain elevation field;

[0103] According to the median particle size of sediment, the sediment angle of repose formula is used: , where the sediment repose angle at the corresponding position is calculated, where is the sediment repose angle at the corresponding position;

[0104] For the four adjacent spatial units adjacent to the current spatial unit, the elevation relationship formula is satisfied: ;in, , is the maximum elevation difference between the current spatial unit and the four adjacent spatial units, is the length of the spatial unit, 、 、 as well as are the coordinate positions of the four adjacent spatial units respectively;

[0105] The current terrain elevation field is corrected using the elevation relationship formula to obtain a corrected terrain elevation field. .

[0106] In an embodiment of the present invention, exemplarily, in an estuary and coastal ecological shoal shaping project, the server performs correction work based on the second association relationship for the determined current terrain elevation field.

[0107] The server first obtains the second association from the database storing the model and relationship , which clarifies the median sediment particle size Elevation of tidal flats and wetlands The corresponding relationship between them.

[0108] Then, the server performs calculations for each spatial unit in the current terrain elevation field. Assume that the current terrain elevation field covers an area of ​​500m×500m, which is divided into 5×5 spatial units of 100m×100m. Taking the current spatial unit with coordinates (2,2) as an example, the server obtains the elevation value of the unit in the current terrain elevation field. According to the second correlation formula , calculate the median particle size of sediment corresponding to the spatial unit Assume that the calculation 0.3mm.

[0109] Next, the server uses the sediment repose angle formula ,Will Substitute into the formula to calculate the sediment repose angle at the corresponding position (2,2) After calculation, It is about 38.5 degrees. Knowing that the length of the space unit L is 100m, the server further calculates the maximum elevation difference between the current space unit and the four adjacent space units. .

[0110] Then, the server makes a judgment based on the elevation relationship formula for the four adjacent spatial units adjacent to the current spatial unit (2,2), namely the spatial units with coordinates (1,2), (3,2), (2,1), and (2,3). 、 、 、 Assume that the calculation is , exceeding , then the server will adjust the elevations of the two units to satisfy the elevation relationship formula.

[0111] The server performs the above operations on every spatial unit in the current terrain elevation field and its adjacent units. Through continuous judgment and adjustment, the entire current terrain elevation field is corrected using the elevation relationship formula, ultimately obtaining a revised terrain elevation field. This revised terrain elevation field takes into account sediment stability factors, ensuring that the terrain will not collapse due to excessive elevation differences in actual conditions, providing a stable and reliable terrain foundation for subsequent ecological shoal shaping work.

[0112] In an embodiment of the present invention, the modified terrain elevation field is corrected based on the preset rules for the skeleton of the ecological shoal to be shaped to obtain the final designed elevation field, which can be implemented through the following examples.

[0113] The main skeleton and branch skeletons are set for the ecological shallows shaping area to be carried out, wherein the position equation of the main skeleton is: ,in, , the position equation of the support skeleton is: ,in, ;

[0114] By formula: , calculate the design elevation field of the ecological shallows shaping area to be carried out, where This is the final designed elevation field.

[0115] In an embodiment of the present invention, for example, in an estuary and coastal ecological shoal shaping project, after the server completes the correction of the terrain elevation field based on sediment stability, it continues to further process the corrected terrain elevation field according to the preset rules of the skeleton for the ecological shoal shaping to be carried out, so as to obtain the final design elevation field.

[0116] The server first sets up the main skeleton and branch skeletons for the ecological shallows shaping area to be carried out according to the project plan. Assume that the ecological shallows shaping area to be carried out is an estuary mudflat with a length of 800m and a width of 600m. The server sets up two main skeletons and three branch skeletons. For the main skeleton, its position equation is ,in , the value range of x is (100,700). For example, the first main skeleton The top elevation z_1m is set to be 2.5 meters above the average high tide level; the second main frame The top elevation is also 2.5 meters above the average high tide level. For the support frame, the position equation is ,in , the value range of x is also (100,700). For example, the first branch skeleton , top elevation Set to 1.5 meters above the average high tide level; the second support frame The top elevation is 1.5 meters above the average high tide level; the third support frame The top elevation is 1.5 meters above the average high tide level.

[0117] Next, the server uses the formula This formula is used to calculate the design elevation field of the ecological shallows shaping area to be carried out. ,Right now When calculating, first determine an intermediate elevation field ,in The middle elevation value. and Represent the coordinate information of points P and S respectively, and They are the elevation at the corresponding position and the elevation at the corresponding position. will be and Coverage, based on this calculation result, the intermediate elevation field Processing is performed to obtain the final design elevation field .in, This is the terrain elevation field that was corrected based on sediment stability. Taking the coordinate point (300,100) as an example, in the corrected terrain elevation field The elevation value corresponding to this point is The server determines whether the point is located on the main frame or the branch frame. After calculation, if the point is not located on the main frame or the branch frame, then the design elevation field The elevation value of this point follows Taking the coordinate point (400,90) as an example, it is located on the second support frame, so in the design elevation field In the example, the elevation of this point is set to 1.5 meters above the top of the second support frame (i.e. value of ).

[0118] The server performs such judgment and assignment operations on all coordinate points in the entire ecological shallows shaping area to be carried out, and obtains the final design elevation field through calculation of the above formula. This designed elevation field comprehensively considers terrain stability and the skeleton rules for shaping ecological shoals, providing an accurate and engineering-compliant terrain foundation for subsequent tidal sediment numerical simulations and ecological assessments. It will guide the actual shaping of ecological shoals to achieve the desired ecological and engineering effects.

[0119] In order to more clearly describe the technical solution provided by the embodiment of the present invention, a relatively complete implementation method is provided below. Figure 2 , Figure 2Schematic diagram of the technical framework for shaping estuarine and coastal ecological shoals by coordinating water and sediment regulation and habitat creation provided by an embodiment of the present invention.

[0120] Step 1: Determine the spatial scope of the ecological shoal shaping area. (This symbol represents a three-dimensional terrain field array), x and y represent the horizontal and vertical coordinates, which are the spatial range of the shaping area, x∈( ,), y∈( ), z0 is the vertical coordinate, indicating the initial terrain elevation, ∈( The ecological shoal shaping area is divided into spatial units (which can be delineated into rectangular units with a unit length of L), where (x, y) is the center coordinate point of the unit.

[0121] The spatial extent of the ecological shoal to be reshaped can be determined by the following methods: First, survey and monitor existing and developing tidal flats along the estuary coast within the study area to clarify the spatial extent of existing tidal flats. Simulate regional hydrological sediment and riverbed evolution to identify mainstream deep troughs and the future development of tidal flats. Combined with local land and space planning and tidal flat utilization planning, the spatial extent of the ecological shoal to be reshaped can be determined. (The boundaries of developing tidal flats can be appropriately expanded based on relevant national and spatial planning to serve as the spatial extent of the ecological shoal to be reshaped. Here, the horizontal and vertical extents (i.e., x and y) are determined.)

[0122] Step 2: Conduct current ecological survey and analysis and assessment.

[0123] (1) Based on the well-developed tidal flats and wetlands around the area to be shaped identified in step 1, select several tidal flats with stable morphology and good habitats, conduct an assessment of the current status of the ecological environment, and provide target habitats for the subsequent shaping of ecological shallows.

[0124] (2) Collect and obtain remote sensing image data (such as GF-1 satellite data, QuickBird imagery, and LiDAR data) of the tidal flat wetlands in the study area, DEM elevation data, and survey and monitoring data on wetland plants, benthic organisms, swimming organisms, and birds in the tidal flat wetlands, as well as measured tidal flat elevation and bed sand gradation data. Use support vector machines and deep learning methods to interpret remote sensing impact data to obtain information on the types and biomass of wetland plants and benthic organisms. Use tools such as ArcGIS to process DEM elevation data and obtain tidal flat wetland elevation information. Based on the measured bed sand gradation data, remote sensing images, and DEM elevation data, more detailed and detailed information on the composition of tidal flat wetland sediments can be obtained.

[0125] (3) The selected tidal flat wetland is divided into spatial units of a certain size (e.g., 100 m × 100 m). The above-mentioned topographic, sediment, and ecological data are interpolated onto the demarcated units. The topographic data are elevation information, the sediment data are median particle size information (composition of bed sand particle size), and the ecological data are information on the species, abundance, and biomass of wetland plants, benthic organisms, nematodes, birds, and other species.

[0126] (4) Based on the above data, calculate the Shannon-Wiener index (which is used to characterize the biodiversity of the community) for each spatial unit.

[0127] (5) Based on the above steps, the topographic, sediment, and ecological data of each spatial unit were obtained, including: elevation, median particle size, species, abundance, and biomass of wetland plants, benthic organisms, nematodes, and birds, as well as information such as the Shannon-Wiener index.

[0128] Step 3: Analyze the relationship between tidal flat wetland elevation, particle size composition, and ecological characteristics. Using a structural equation model (SEM), we established relationships between the species, abundance, and biomass of wetland plants, benthic organisms, nematodes, and birds, as well as the Shannon-Wiener Index (dependent variables) and tidal flat wetland elevation and sediment median particle size (independent variables). The specific process is as follows:

[0129] (1) Establish hypothetical pathways, including the impact of terrain elevation and median sediment particle size on wetland plants, benthic organisms, nematodes, birds, etc.; the impact of wetland plants on benthic organisms, nematodes, birds, etc.; the impact of benthic organisms on nematodes, birds, etc.; the impact of nematodes on birds, etc.

[0130] (2) Based on the hypothesized paths, perform model fitting and path analysis using statistical software (such as the lavaan package or piecewiseSEM package in R), test the significance of each path, and adjust the model based on the results. Based on the model fitting results, remove insignificant hypothesized paths and add hypothesized paths that may have been missed. Evaluate the rationality of the model using model goodness-of-fit indicators (such as AIC value and R² value).

[0131] (3) The data information obtained can be divided into two parts, one for model calibration and the other for model verification. Figure 3 , Figure 3 Schematic diagram of the structural equation model analysis framework provided by an embodiment of the present invention.

[0132] Step 4: Based on the sediment median particle size and elevation information of each spatial unit, a statistical regression analysis method was used to establish the relationship between sediment median particle size (dependent variable) and tidal flat wetland elevation (independent variable), namely: ,in, represents the median particle size of sediment, and z represents the elevation of the tidal flat wetland.

[0133] Step 5: For areas where ecological shoal formation is required , where x∈( ,), y∈( ), ∈( ), a rule-based random generation method is used to generate the design elevation field , where, ,x∈( ,), y∈( ), ∈( ), the specific steps are as follows:

[0134] (1) First, based on the initial terrain elevation field, a random elevation field is generated using a random generation algorithm. For areas where ecological shoals need to be shaped , use Perlin noise, Simplex noise and other algorithms to generate elevation values, and normalize the randomly generated elevation values ​​to make them range between 0 and 1, and you can get ,in .

[0135] Therefore, the new terrain elevation field The following formula can be used for calculation: ,

[0136] in, represents the vector terrain field array. The above formula is the sum of terrain field vectors; z0-r is the terrain elevation after superimposing the random elevation field; α is the random terrain elevation adjustment coefficient, which can be determined based on the initial terrain field elevation range and the generated random elevation field range.

[0137] (2) Considering the stability characteristics of mudflat surface sediment, the randomly generated terrain elevation field Correction is performed to obtain the corrected terrain elevation field: The reason for the correction is that if the elevation difference between the randomly generated elevation fields of adjacent cells is too large, it will cause terrain collapse. Here, the sediment repose angle formula is used to consider the stability of the mudflat surface sediment and correct the terrain elevation field.

[0138] a) Based on the relationship between the tidal flat terrain elevation and the median sediment particle size established in step 4 above, , calculate and obtain the random elevation field Corresponding location The median particle size of sediment at is as follows: ,

[0139] b) Calculate the spatial unit according to the sediment repose angle formula (such as Zhang Hongwu's formula, Cheng Niansheng's formula, etc.) Sediment repose angle ,in, is the underwater angle of repose of sediment, in degrees; is the median particle size of sediment, in mm.

[0140] c) For spatial units The four adjacent space units 、 、 、 , which can be calculated separately from the spatial unit Maximum elevation difference: , where L is the length of the space unit.

[0141] Space Unit 4 space units adjacent to the surrounding area 、 、 、 The elevation relationship should satisfy the following relationship: ;

[0142] d) Assume that A spatial unit r starts iterative calculation of the elevation field of the shallow wetland Correction is performed to obtain the corrected morphologically stable shallow wetland elevation field .

[0143] (3) Topographic correction considering the rules (the skeleton of ecological shoal topography, used to reduce erosion and promote sedimentation).

[0144] a) For areas where ecological shoal formation is required, several main frameworks are set up and support skeleton The main frame and supporting frame are made of concrete or stone embankments, which will not be directly carried away by the tide. After the arrangement, a slow flow area will be formed to reduce erosion and promote sedimentation. Figure 4 , Figure 4 A schematic diagram of the area where ecological shoal shaping is to be carried out according to an embodiment of the present invention.

[0145] b) Conduct numerical simulations of tidal currents and sediment in the study area. Combined with the distribution of tidal currents and changes in sediment scouring and deposition, identify the main trough (which can be set as a strip-shaped area with a flow rate greater than 0.5 m / s and a tendency to scour). The main framework of the ecological shoal shaping area, i.e., the beach protection dike, can be set parallel to the main trough and at a certain distance (e.g., 200 m) to one side of the ecological shoal. The length and number of the dikes can be set based on actual conditions. The equation for setting the main framework position is: ,in, The top elevation of the main frame (beach protection embankment) can be set to several meters above the high tide level, such as , you can also set a different top elevation for each main frame (beach protection dike).

[0146] Several branch frames can be set in the direction perpendicular to the main frame (beach protection dike). The length and number of the branch frames can be set based on actual conditions. The equation for setting the branch frame position is: ,in, The top elevation of the support frame can be set to several meters above the high tide level, such as , you can also set a different top elevation for each support frame.

[0147] c) Correct the elevations of the main frame P and the support frame S and set them as 、 The final design elevation field is obtained. .

[0148] (4) Comprehensively consider the random elevation field generation, elevation correction based on the stability of the tidal flat morphology, and elevation correction based on rules to obtain the final design elevation field for the area where ecological shoal shaping is required. . ,in, Represents the vector terrain field array, and the above formula is the addition of terrain field vectors; is the final design elevation field; the symbol “⊇” here indicates the ecological shoal shaping area The elevations of the main frame P and the branch frame S are directly set or overwritten as 、 The symbol “∪” means that both the main skeleton P and the branch skeleton S are considered at the same time.

[0149] Step 6: Design elevation field for the area where ecological shallows need to be shaped , carry out numerical simulation of tidal sediment, calculate the regional erosion and deposition balance, and obtain the elevation field after the ecological shoal is shaped into regional erosion and deposition balance. and corresponding positions Median sediment particle size at .

[0150] Step 7: Elevation field of the ecological shoal shaping area after erosion and deposition balance , median particle size of sediment By inputting the species, abundance, biomass and Shannon-Wiener index of wetland plants, benthic organisms, nematodes and birds obtained in step 3 and the structural equation model of tidal flat wetland elevation and sediment median particle size, the ecological shoal shaping area can be obtained. The types, abundance, biomass and Shannon-Wiener index of wetland plants, benthic organisms, nematodes and birds are analyzed.

[0151] Step 8: Repeat steps 5 to 7 to obtain several design elevation fields. (where i represents several scenarios), considering the design elevation field with the best ecological indicators (including more species, more abundance, larger biomass, higher Shannon-Wiener index, etc.) Shaped as target habitat.

[0152] Please refer to Figure 5 , Figure 5 An embodiment of the present invention provides an estuary and coastal ecological shoal shaping device 110 for coordinated water and sediment regulation and habitat creation, comprising:

[0153] The acquisition module 1101 is used to determine the spatial scope of the ecological shoal shaping area to be carried out, and the surrounding area of ​​the ecological shoal shaping area to be carried out includes multiple candidate shoals; select at least two reference shoals from the multiple candidate shoals, and obtain the tidal flat wetland elevation information, sediment median particle size and ecological characteristic data of the at least two reference shoals; divide the at least two reference shoals into multiple spatial units, and interpolate each of the spatial units based on the tidal flat wetland elevation information, the sediment median particle size and the ecological characteristic data to obtain the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each of the spatial units; construct a target structural equation model that characterizes the first correlation relationship included in each of the spatial units, and the first correlation relationship is the reference tidal flat wetland elevation information, target sediment median particle size and target ecological characteristic data corresponding to each of the spatial units. The correlation between the median particle size and the target ecological characteristic data; constructing a second correlation relationship included in each of the spatial units based on a statistical regression analysis strategy, wherein the second correlation relationship is the correlation relationship between the reference tidal flat wetland elevation information corresponding to each of the spatial units and the target sediment median particle size; determining the current terrain elevation field of the ecological shoal shaping area to be carried out based on a random generation algorithm, and correcting the current terrain elevation field using the second correlation relationship and preset rules to obtain the design elevation field of the ecological shoal shaping area to be carried out; performing tidal sediment numerical simulation on the design elevation field to obtain the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment in the regional scouring and sedimentation balance state of the ecological shoal shaping area to be carried out; inputting the scouring and sedimentation balance elevation field and the median particle size of the scouring and sedimentation balance sediment into the target structural equation model to obtain candidate ecological characteristic data;

[0154] The shaping module 1102 is used to repeatedly execute the step of determining the current terrain elevation field of the ecological shallows shaping area to be carried out based on the random generation algorithm for a preset number of times, to the step of inputting the erosion and deposition balance elevation field and the erosion and deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, and obtain multiple candidate ecological characteristic data; from the multiple candidate ecological characteristic data, the target design elevation field corresponding to the target ecological characteristic data is screened out according to the preset ecological indicators as the target habitat of the ecological shallows shaping area to be carried out for shaping.

[0155] It should be noted that the implementation principles of the aforementioned estuarine and coastal ecological shallows shaping device 110 for coordinated water and sediment control and habitat creation can be referenced to the implementation principles of the aforementioned estuarine and coastal ecological shallows shaping method for coordinated water and sediment control and habitat creation, and will not be elaborated upon here. It should be understood that the division of the various modules of the above device is merely a division of logical functions. In actual implementation, they can be fully or partially integrated into a single physical entity, or physically separated. Furthermore, these modules can be implemented entirely as software invoked by processing elements; entirely as hardware; or partially as software invoked by processing elements, while others can be implemented in hardware. For example, the estuarine and coastal ecological shallows shaping device 110 for coordinated water and sediment control and habitat creation can be a separate processing element, or integrated into a chip of the aforementioned device. Furthermore, it can be stored in the form of program code in the memory of the aforementioned device, with a processing element of the aforementioned device invoking and executing the functions of the aforementioned estuarine and coastal ecological shallows shaping device 110 for coordinated water and sediment control and habitat creation. The implementation of the other modules is similar. Furthermore, these modules can be fully or partially integrated together, or implemented independently. The processing element described herein may be an integrated circuit having signal processing capabilities. In the implementation process, each step of the above method or each module above may be completed by hardware integrated logic circuits in the processor element or software instructions.

[0156] For example, the above modules may be one or more integrated circuits configured to implement the above methods, such as one or more application specific integrated circuits (ASICs), one or more digital signal processors (DSPs), or one or more field programmable gate arrays (FPGAs). For another example, when a module is implemented by scheduling program code on a processing element, the processing element may be a general-purpose processor, such as a central processing unit (CPU) or other processor that can call program code. For another example, these modules may be integrated together and implemented in the form of a system-on-a-chip (SOC).

[0157] The embodiment of the present invention provides a computer device 100, which includes a processor and a non-volatile memory storing computer instructions. When the computer instructions are executed by the processor, the computer device 100 executes the aforementioned estuary and coastal ecological shoal shaping device 110 for coordinated water and sand regulation and habitat creation. Figure 6 As shown, Figure 6This is a block diagram of a computer device 100 according to an embodiment of the present invention. The computer device 100 includes an estuary and coastal ecological shoal shaping device 110 for coordinated water and sediment regulation and habitat creation, a memory 111 , a processor 112 , and a communication unit 113 .

[0158] In order to achieve data transmission or interaction, the memory 111, the processor 112 and the communication unit 113 are electrically connected to each other directly or indirectly. For example, the electrical connection between these components can be achieved through one or more communication buses or signal lines. The estuary and coastal ecological shallows shaping device 110 for coordinated water and sand regulation and habitat creation includes at least one software function module that can be stored in the memory 111 in the form of software or firmware or solidified in the operating system (OS) of the computer device 100. The processor 112 is used to execute the estuary and coastal ecological shallows shaping device 110 for coordinated water and sand regulation and habitat creation stored in the memory 111, such as the software function modules and computer programs included in the estuary and coastal ecological shallows shaping device 110 for coordinated water and sand regulation and habitat creation.

[0159] An embodiment of the present invention provides a readable storage medium, which includes a computer program. When the computer program is running, it controls the computer device where the readable storage medium is located to execute the aforementioned estuary and coastal ecological shoal shaping device 110 that coordinates water and sand regulation and habitat creation.

[0160] For illustrative purposes, the foregoing description has been made with reference to specific embodiments. However, the above illustrative discussion is not intended to be exhaustive or to limit the present disclosure to the precise forms disclosed. Numerous modifications and variations are possible in light of the above teachings. These embodiments have been selected and described in order to best illustrate the principles of the present disclosure and its practical application, thereby enabling those skilled in the art to best utilize the present disclosure and to utilize various embodiments with various modifications as appropriate for the specific application contemplated.

Claims

1. A method for shaping estuarine and coastal ecological shoals by coordinating water and sediment regulation with habitat creation, characterized in that: include: Determine the spatial scope of the area where ecological shoal shaping is to be carried out, wherein the area where ecological shoal shaping is to be carried out includes a plurality of candidate shoals; Selecting at least two reference tidal flats from the plurality of candidate tidal flats, and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats; Dividing the at least two reference tidal flats into a plurality of spatial units, and interpolating each of the spatial units based on the tidal flat wetland elevation information, the sediment median particle size, and the ecological characteristic data, to obtain the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units; Constructing a target structural equation model representing a first correlation relationship included in each of the spatial units, wherein the first correlation relationship is a correlation relationship between reference tidal flat wetland elevation information, target sediment median particle size, and target ecological characteristic data corresponding to each of the spatial units; Constructing a second association relationship for each of the spatial units based on a statistical regression analysis strategy, wherein the second association relationship is an association relationship between the reference tidal flat wetland elevation information and the target sediment median particle size corresponding to each of the spatial units; Determining a current terrain elevation field of the area to be subjected to ecological shoal shaping based on a random generation algorithm, and modifying the current terrain elevation field using the second association relationship and a preset rule to obtain a design elevation field of the area to be subjected to ecological shoal shaping; Conducting numerical simulation of tidal current and sediment on the designed elevation field to obtain the scouring and deposition equilibrium elevation field and the median particle size of scouring and deposition equilibrium sediment in the area to be developed for ecological shoal shaping under the regional scouring and deposition equilibrium state; Inputting the scour-sediment balance elevation field and the scour-sediment balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data; Repeating the step of determining the current terrain elevation field of the area to be ecological shoal shaping based on the random generation algorithm for a preset number of times, until the step of inputting the erosion-deposition balance elevation field and the erosion-deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, to obtain a plurality of candidate ecological characteristic data; A target design elevation field corresponding to the target ecological characteristic data is selected from the plurality of candidate ecological characteristic data according to preset ecological indicators to be used as the target habitat for the ecological shallows shaping area to be developed.

2. The method according to claim 1, characterized in that The step of selecting at least two reference tidal flats from the plurality of candidate tidal flats and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats includes: selecting at least two reference tidal flats from the plurality of candidate tidal flats according to preset tidal flat morphological stability indicators and good habitat indicators; Acquiring remote sensing image data, DEM elevation data, wetland plant data, benthic organism data, swimming organism data, bird data, mudflat elevation data, and bed sand gradation data corresponding to the at least two reference mudflats; Using ArcGIS to process the DEM elevation data to obtain the elevation information of the tidal flat wetland; Determining the median particle size of the sediment based on the remote sensing image data, the tidal flat elevation data, and the bed sand gradation data; The remote sensing image data is interpreted by a support vector machine combined with a deep learning strategy, and combined with the wetland plant data, benthic organism data, swimming organism data and bird data, the species, abundance and biomass of wetland plants, benthic organisms, swimming organisms and birds are obtained as the ecological characteristic data.

3. The method according to claim 2, characterized in that The target structural equation model representing the first correlation relationship included in each of the spatial units is constructed, wherein the first correlation relationship is the correlation relationship between the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units, including: The Shannon-Wiener index corresponding to each of the spatial units is calculated based on the species, abundance, and biomass of the wetland plants, benthic organisms, nematodes, and birds, and is used to characterize the biodiversity of the community; Constructing a hypothetical path, the hypothetical path including the relationship between the elevation information of the tidal flat wetland, the median particle size of the sediment, the ecological characteristic data, and the biodiversity; Using preset statistical software to perform model fitting and path analysis on the hypothesized path, and testing the significance of each path to obtain a model fitting result; Model tuning is performed based on the model fitting results and the significance of each path to obtain the target structural equation model.

4. The method according to claim 1, wherein The method of determining the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm, and correcting the current terrain elevation field using the second association relationship and a preset rule to obtain a design elevation field of the area to be ecological shoal shaping includes: Determine the current terrain elevation field of the area to be subjected to ecological shoal shaping based on a random generation algorithm; Calculating the sediment repose angle of the spatial unit corresponding to the current terrain elevation field based on the second association relationship, and correcting the current terrain elevation field based on the sediment repose angle to obtain a corrected terrain elevation field; The modified terrain elevation field is modified based on preset rules for the skeleton of the ecological shoal to be shaped, so as to obtain the final designed elevation field.

5. The method according to claim 4, characterized in that The determining of the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm includes: Acquire an initial terrain elevation field of the area where the ecological shoal shaping is to be carried out, and process the initial terrain elevation field using a preset random noise algorithm to obtain a random terrain elevation field; By formula: , calculate the current terrain elevation field, where is the current terrain elevation field, is the initial terrain elevation field, is the random terrain elevation field, is a vector terrain field array, is the random terrain elevation adjustment coefficient.

6. The method according to claim 5, characterized in that The calculating of the sediment repose angle of the spatial unit corresponding to the current terrain elevation field based on the second association relationship, and correcting the current terrain elevation field based on the sediment repose angle to obtain a corrected terrain elevation field, includes: The second association relationship is obtained, where the second association relationship is: ,in, is the median particle size of sediment, is the elevation of the tidal flat wetland; According to the second association relationship, through the formula: Calculate and obtain the median particle size of sediment in the spatial unit corresponding to the current terrain elevation field; wherein, is the coordinate position of the current spatial unit corresponding to the current terrain elevation field; According to the median particle size of sediment, the sediment angle of repose formula is used: , where the sediment repose angle at the corresponding position is calculated, where is the sediment repose angle at the corresponding position; For the four adjacent spatial units adjacent to the current spatial unit, the elevation relationship formula is satisfied: ;in, is the maximum elevation difference between the current spatial unit and the four adjacent spatial units, is the length of the spatial unit, as well as are the coordinate positions of the four adjacent spatial units respectively; The current terrain elevation field is corrected using the elevation relationship formula to obtain a corrected terrain elevation field. .

7. The method according to claim 6, characterized in that The modified terrain elevation field is modified based on the preset rules for the skeleton of the ecological shoal to be shaped to obtain the final designed elevation field, including: The main skeleton and branch skeletons are set for the ecological shallows shaping area to be carried out, wherein the position equation of the main skeleton is: ,in, , the position equation of the support skeleton is: ,in, ; By formula: , calculate the design elevation field of the ecological shallows shaping area to be carried out, where This is the final designed elevation field.

8. A device for shaping estuary and coastal ecological shoals that coordinates water and sand regulation with habitat creation, characterized in that: include: An acquisition module is used to determine the spatial scope of the area to be developed for ecological shoal shaping, wherein the area to be developed for ecological shoal shaping includes a plurality of candidate shoals; Selecting at least two reference tidal flats from the plurality of candidate tidal flats, and obtaining tidal flat wetland elevation information, sediment median particle size, and ecological characteristic data of the at least two reference tidal flats; Dividing the at least two reference tidal flats into a plurality of spatial units, and interpolating each of the spatial units based on the tidal flat wetland elevation information, the sediment median particle size, and the ecological characteristic data, to obtain the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units; Constructing a target structural equation model representing a first correlation relationship included in each of the spatial units, wherein the first correlation relationship is a correlation relationship between the reference tidal flat wetland elevation information, the target sediment median particle size, and the target ecological characteristic data corresponding to each of the spatial units; constructing a second correlation relationship included in each of the spatial units based on a statistical regression analysis strategy, wherein the second correlation relationship is a correlation relationship between the reference tidal flat wetland elevation information and the target sediment median particle size corresponding to each of the spatial units; Determine the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm, and modify the current terrain elevation field using the second association relationship and a preset rule to obtain a design elevation field of the area to be ecological shoal shaping; perform tidal sediment numerical simulation on the design elevation field to obtain an erosion-deposition equilibrium elevation field and an erosion-deposition equilibrium sediment median particle size of the area to be ecological shoal shaping under a regional erosion-deposition equilibrium state; input the erosion-deposition equilibrium elevation field and the erosion-deposition equilibrium sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data; a shaping module, configured to repeatedly execute the steps of determining the current terrain elevation field of the area to be ecological shoal shaping based on a random generation algorithm a preset number of times, to the step of inputting the erosion-deposition balance elevation field and the erosion-deposition balance sediment median particle size into the target structural equation model to obtain candidate ecological characteristic data, thereby obtaining a plurality of candidate ecological characteristic data; A target design elevation field corresponding to the target ecological characteristic data is selected from the plurality of candidate ecological characteristic data according to preset ecological indicators to be used as the target habitat for the ecological shallows shaping area to be developed.

9. A computer device, characterized in that: The computer device includes a processor and a non-volatile memory storing computer instructions. When the computer instructions are executed by the processor, the computer device executes the method according to any one of claims 1 to 7.

10. A readable storage medium, characterized in that: The readable storage medium includes a computer program, and when the computer program is executed, the computer device where the readable storage medium is located is controlled to execute the method according to any one of claims 1 to 7.

Images