CN115408832A - River bank ecological slope protection anti-impact flow velocity rechecking method - Google Patents

Abstract

The invention relates to the technical field of hydraulic engineering, in particular to an anti-impact flow velocity rechecking method for ecological river bank revetment. Establishing a terrain grid with special codes of boundary nodes by using SMS, and importing the terrain grid into MIKE21 to calculate a river flow field to obtain a calculation result file; reading a terrain grid, extracting grid boundary coordinates according to codes, and importing the grid boundary coordinates into a QGIS; in the QGIS, left and right bank attributes of the boundary coordinates are marked, and a boundary grid coordinate attribute table is generated by sequencing according to the left and right bank attributes; extracting boundary flow velocity from the calculation result file by using a boundary coordinate attribute table; and (5) corresponding the boundary flow velocity with the number of the engineering pile, and rechecking whether the boundary flow velocity exceeds the impact flow velocity of the ecological slope protection of the river bank. By the aid of an MIKE21 program expansion package MIKE SDK provided by Danish water conservancy research and the combination of spatial geographic analysis software QGIS, rapid extraction and analysis of the bank slope flow velocity of the two-dimensional flow field calculated by the MIKE21 are achieved, and design efficiency of ecological bank protection is greatly improved.

Description

Review method for anti-scour flow velocity of river bank ecological slope protection

technical field

The invention relates to the technical field of water conservancy engineering, in particular to a method for reviewing anti-scour flow velocity of river bank ecological slope protection.

Background technique

As an important protective barrier for flood control and an important part of the ecosystem, the river shoreline has received more and more attention from the government. In the new era, there are new requirements for shoreline improvement. The existing coastline improvement is no longer limited to traditional slope protection and bank protection projects, but has been transformed into a series of engineering measures to create a riverside ecological landscape belt. Compared with traditional hard slope protection such as masonry slope protection, dry stone slope protection, concrete slope protection, etc., ecological revetment such as flexible ecological water and soil protection blankets, ecological vegetation mats or Renault pads are increasingly used due to their ecological friendliness. It has been adopted in embankment slope protection engineering.

Compared with the traditional hard slope protection method, the ecological revetment has certain requirements on the near-shore flow velocity, which requires the use of a two-dimensional hydrodynamic model to calculate the water flow velocity at the bank slope, and compare it with the anti-scouring flow rate of the ecological slope protection, so as to determine the ecological slope protection. Engineering form and suitability. MIKE21 is a two-dimensional hydrodynamic numerical calculation software researched and developed by the Danish Hydraulic Research Institute. It uses the finite element method to discretely solve the two-dimensional shallow water equation. After determining the boundary conditions, the water depth and flow velocity results of the two-dimensional flow field can be obtained. MIKE21 has been widely used in China. MIKE21 has been applied to major Chinese rivers to solve major engineering problems, such as the simulation of hydrodynamics and salinity in the Yangtze Estuary, the simulation of water flow and sediment in the Pearl River Estuary, and the simulation of water flow and waves in the Bohai Bay. . It can be said that the calculation of two-dimensional hydrodynamic process through MIKE21 has become one of the industry standards for hydrodynamic calculation in the water conservancy industry.

The design of MIKE21 was completed in the last century, and the graphical interface it provides is relatively old and difficult to use, especially in the extraction of river bank flow velocity distribution, which is very cumbersome.

Contents of the invention

The purpose of the present invention is to address the defects of the prior art, to provide a river bank ecological slope protection anti-scour flow velocity review method, based on the MIKE21 program expansion package MIKE SDK provided by the Danish Water Conservancy Research Institute and combined with the spatial geographic analysis software QGIS, to achieve the calculation of MIKE21 The rapid extraction and analysis of the two-dimensional flow field bank slope velocity greatly improves the design efficiency of ecological revetment.

A kind of river bank ecological slope protection anti-scour flow velocity review method of the present invention comprises:

Use the surface water simulation system to generate a grid file based on the topographic map, and indicate the boundary conditions of the study area in the surface water simulation system, import it into MIKE21, determine the upstream and downstream boundaries, perform calculations in MIKE21, and obtain dfsu including flow velocity and water depth information document;

Read the grid coordinates from the calculation result dfsu file, use the boundary type code to get the boundary grid coordinates from the calculation result dfsu file, and generate a CSV file;

Import the generated coordinate CSV file into QGIS, add the left and right bank attributes in the attribute table, and export the generated attribute table as an xlsx file after sorting according to the left and right bank attributes;

Read the terrain grid coordinates, match them with the left and right bank coordinates, obtain the index information of the boundary grid coordinates, read the velocity information of the calculation results, and use the index information of the boundary grid coordinates to obtain the boundary flow velocity information;

Sort the coordinates of the boundary grid, calculate the cumulative length of the sorted coordinate points along the river bank axis direction, and correspond to the pile number section that needs to be laid with ecological slope protection, compare the maximum flow velocity on the bank of the corresponding pile number section with the ecological slope protection resistance Evaluate the applicability of ecological slope protection in this pile number section.

More preferably, the boundary conditions of the study area include no-slip boundary, flow boundary and water level boundary.

More preferably, the import into MIKE21 determines the upstream and downstream boundaries, and calculates in MIKE21 to obtain a dfsu file containing flow velocity and water depth information including:

Create a new MIKE21 FM project in the MIKE ZERO software, import the grid file, and determine the upstream and downstream boundaries;

Given the boundary conditions of water level and flow in MIKE21, calculate the results of flow field velocity and water level under different working conditions;

Select the dfsu output, and check the flow rate and water level information in the output content.

More preferably, the given water level and flow boundary conditions in MIKE21 include river bank, river channel inlet and outlet water level and flow conditions.

More preferably, reading the grid coordinates from the calculation result dfsu file includes using the MATLAB mzReadMesh function provided by the MIKE SDK to read the grid coordinates from the calculation result dfsu file.

More preferably, adding the left and right bank attributes in the attribute table includes:

In QGIS, select the scatter points belonging to the left bank and mark the left bank in the attribute table, and select the scatter points belonging to the right bank to mark the right bank in the attribute table.

More preferably, the said read terrain grid coordinates, match them with the left and right bank coordinates, obtain the index information of the boundary grid coordinates, read the calculation result velocity information, and use the index information of the boundary grid coordinates to obtain the boundary flow velocity information include:

Read the boundary coordinate attribute table to obtain the coordinate information of the left and right banks;

Read the coordinates in the grid file, match them with the coordinates of the left and right banks, and obtain the index information of the grid points at the boundary in the calculation result dfsu file;

Use the index information to get the boundary flow velocity information from the calculation result dfsu file.

More preferably, sorting the boundary grid coordinates includes:

If the river flows from west to east, sort the coordinate points and corresponding boundary flow velocities according to the X coordinates of the left and right bank boundaries from small to large;

If the river flows from east to west, the coordinate points and corresponding boundary flow velocities are sorted according to the order of the X coordinates of the left and right bank boundaries from large to small.

More preferably, if the maximum flow velocity on the bank of the corresponding pile number section is not greater than the anti-scour velocity of the ecological slope protection, then the ecological slope protection is applicable to this pile number section; Ecological slope protection is not applicable in this chainage section.

More preferably, grid node coordinates are filled in the CSV file in the form of a table.

The beneficial effects of the present invention are:

1. Through the MIKE21 program expansion package MIKE SDK provided by the Danish Water Conservancy Research Institute and combined with the spatial geography analysis software QGIS, the present invention realizes the rapid extraction and analysis of the two-dimensional flow field bank slope velocity calculated by MIKE21, and greatly improves the ecological revetment design efficiency.

2. The present invention saves the cumbersome post-processing method in the calculation of the traditional shore flow velocity, and directly extracts the shore flow velocity by using the MATLAB API provided by the MIKE SDK, and directly corresponds to the pile number in the engineering design practice, which is more convenient Review of anti-scour velocity of ecological slope protection.

3. Through special encoding for the boundary nodes when constructing the grid, the boundary nodes are distinguished from the nodes in the river channel, which is convenient for post-processing.

4. By importing the boundary grid coordinates into QGIS for post-processing, the left and right banks in the boundary grid coordinates are effectively distinguished, which facilitates the calculation of the left and right bank velocity respectively.

5. The boundary grid coordinates are post-processed by QGIS to ensure that the boundary grid nodes are coded continuously along the embankment axis, which facilitates the calculation of the boundary velocity stake number.

Description of drawings

Fig. 1 is a schematic flow chart of the present invention;

Fig. 2 is a terrain grid diagram of the present invention.

Detailed ways

In order to make the technical problems, technical solutions and beneficial effects to be solved by the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present application, and are not intended to limit the present application.

In the following description, specific details such as specific system structures and technologies are presented for the purpose of illustration rather than limitation, so as to thoroughly understand the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.

It should be understood that when used in this specification and the appended claims, the term "comprising" indicates the presence of described features, integers, steps, operations, elements and/or components, but does not exclude one or more other Presence or addition of features, wholes, steps, operations, elements, components and/or collections thereof.

It should also be understood that the term "and/or" used in the description of the present application and the appended claims refers to any combination and all possible combinations of one or more of the associated listed items, and includes these combinations.

As used in this specification and the appended claims, the term "if" may be construed, depending on the context, as "when" or "once" or "in response to determining" or "in response to detecting ". Similarly, the phrase "if determined" or "if [the described condition or event] is detected" may be construed, depending on the context, to mean "once determined" or "in response to the determination" or "once detected [the described condition or event] ]” or “in response to detection of [described condition or event]”.

In addition, in the description of the specification and appended claims of the present application, the terms "first", "second", "third" and so on are only used to distinguish descriptions, and should not be understood as indicating or implying relative importance.

Reference to "one embodiment" or "some embodiments" or the like in the specification of the present application means that a particular feature, structure, or characteristic described in connection with the embodiment is included in one or more embodiments of the present application. Thus, appearances of the phrases "in one embodiment," "in some embodiments," "in other embodiments," "in other embodiments," etc. in various places in this specification are not necessarily All refer to the same embodiment, but mean "one or more but not all embodiments" unless specifically stated otherwise. The terms "including", "comprising", "having" and variations thereof mean "including but not limited to", unless specifically stated otherwise. "Multiple" means "two or more".

Embodiment one

Fig. 1 shows the structural representation of a kind of riverbank ecological slope protection anti-scour velocity review method that the preferred embodiment of the application (Fig. 1 shows the first embodiment of the application) provides, for ease of description, only shows the same as the present application The relevant part of the embodiment is described in detail as follows:

A kind of river bank ecological slope protection anti-scour flow velocity review method of the present invention comprises:

Step 1: Use the Surface Water Modeling System (SMS, Surface Water Modeling System) to generate grid files based on topographic maps, and indicate the boundary conditions of the study area in the Surface Water Modeling System, including no-slip boundaries, flow boundaries, and water level boundaries. Create a new MIKE21FM project in the MIKE ZERO software, import the grid file, determine the upstream and downstream boundaries, set the water level and flow boundary conditions in MIKE21, and perform calculations to obtain the dfsu file including flow velocity and water depth information. Including calculating the results of flow field velocity and water level under different working conditions, selecting dfsu output, and checking the flow velocity, water level and other information in the output content;

Step 2: Use the MATLAB mzReadMesh function provided by MIKE SDK to read the grid coordinates from the calculation result dfsu file, and use the boundary type code to obtain the boundary grid coordinates from the calculation result dfsu file (the code of the river boundary in the grid file of MIKE21 1), and generate a CSV file; the format of the CSV file is shown in Table 1 (fill in the grid node coordinates, see Table 3 for specific examples).

Table 1: Boundary coordinate format table

x Y … …

Step 3: Import the generated coordinate CSV file into QGIS, and add the left and right bank attributes in the attribute table. The format of the attribute table is shown in Table 2. In QGIS, select the scatter points belonging to the left bank and mark the left bank (0) in the attribute table, and then select the scatter points belonging to the right bank and mark the right bank (1) in the attribute table. The generated attribute table is sorted according to the attributes of the left and right banks, and the scattered points located on the left and right banks are respectively rounded up and exported as an xlsx file (format).

Table 2: Boundary coordinate attribute table

x Y left and right banks … … 0/1

*0 represents the left bank, 1 represents the right bank

Step 4: Read the terrain grid coordinates, match them with the left and right bank coordinates, obtain the index information of the boundary grid coordinates, read the velocity information of the calculation results, and use the index information of the boundary grid coordinates to obtain the boundary flow velocity information; specifically include :

Step 401: Use the DFSU MATLAB API provided by MIKE SDK to read the result file of MIKE21 calculation, first check the storage location of the speed information (Current speed) in the dfsu file, and then read the speed information in the calculation result.

Step 402: Use MATLAB to read the boundary coordinate attribute table to obtain the coordinate information of the left and right banks. Use the MATLAB mzReadMesh function provided by MIKE SDK to read the coordinates in the grid file, match them with the coordinates of the left and right banks, and obtain the index information of the grid points on the boundary in the calculation result dfsu file. Using the index information, the boundary flow velocity information is obtained from the calculation result dfsu file.

Step 5: Sort the coordinates of the boundary grid, calculate the cumulative length of the sorted coordinate points along the axis of the river bank, and correspond to the pile number section where ecological slope protection needs to be laid, and compare the maximum flow velocity on the bank of the corresponding pile number section with the same The size of the anti-scour velocity of the ecological slope protection is used to evaluate whether the ecological slope protection is suitable for the pile number section. If the maximum flow velocity on the bank of the corresponding pile number section is not greater than the anti-scour velocity of the ecological slope protection, then the ecological slope protection is suitable for this pile number section; Station segments are not applicable.

Among them, the method of sorting the boundary grid coordinates is as follows:

If the river flows from west to east, sort the coordinate points and the corresponding boundary velocity according to the X coordinates of the left and right bank boundaries; Sort the coordinate points and their corresponding boundary flow velocities in descending order.

Embodiment two

The present invention is described in detail by taking the anti-scour recheck calculation of ecological slope protection in the embankment reinforcement and treatment project of newly built houses in Xiong'an New District as an example, and it also has guiding significance for the anti-scour recheck calculation of other ecological slope protection.

Step 1: Use SMS to generate grid files from terrain scatter points, as shown in Figure 2, create a new MIKE21FM project in MIKE ZERO software, import grid files, determine boundary conditions (river banks, river channel inlet and outlet water level and flow conditions), and calculate different working conditions For the results of flow velocity and water level in the flow field, select dfsu output, and check the flow velocity, water level and other information in the output content.

Step 2: Use the MATLAB mzReadMesh function provided by the MIKE SDK to read the boundary coordinates of the coordinate grid and generate a CSV file. The CSV file of coordinate grid boundary coordinates is shown in Table 3.

Table 3: Boundary grid coordinate table

x Y 526958.8 4322764 526873.1 4322715 526798.3 4322649 … … … … … … 504392.9 4328130 504307.8 4328182 504218.2 4328226

Step 3: Import the generated coordinate CSV file into QGIS, and add the left and right bank attributes in the attribute table. In QGIS, select the scatter points belonging to the left bank and mark the left bank (0) in the attribute table, and then select the scatter points belonging to the right bank and mark the right bank (1) in the attribute table. Export the generated attribute table as an xlsx file after sorting it according to the left and right bank attributes. The coordinate attribute file is shown in Table 4.

Table 4: Attribute table of the left and right banks of the boundary

x Y left and right banks 526958.8 4322764 0 526873.1 4322715 0 526798.3 4322649 0 … … … … … … … … … 503556.2 4326815 1 503502.5 4326899 1 503442.7 4326979 1

*0 represents the left bank, 1 represents the right bank

Step 401: Use DFSU MATLAB API to read the result file calculated by MIKE21. In this model simulation, the speed information (Current speed) is located in the fifth output data set in the dfsu file, and all speed information is read.

Step 402: Read the attribute table xlsx file with MATLAB to obtain the coordinate files of the left and right bank. Use the MATLAB mzReadMesh function provided by MIKE SDK to read the grid coordinates of the calculation result dfsu file, use the coordinate position to match it with the left and right bank coordinates, and obtain the index information of the boundary coordinate points in the calculation result dfsu file grid. Using the index information, the boundary flow velocity information is obtained.

Step 5: The flow direction of Xingaifang Floodway in Xiongan New Area is from west to east, so the boundary coordinate points and corresponding boundary flow velocities are sorted according to the X coordinates of the left and right bank boundaries from small to large. Calculate the cumulative length of the sorted river bank coordinate points along the river bank axis. The embankment reinforcement and treatment project of the newly built flood diversion channel in Xiong'an New Area needs to apply ecological slope protection in the whole section, so 500 meters is the maximum flow velocity in a statistical interval. section applies. The accumulative length of the coordinate points on the river bank is the corresponding chainage of the grid point. From the head of the river embankment, every 500 meters is used as an interval to count the maximum flow velocity in the chainage. See Table 5 and Table 6 for the maximum flow velocity in the chainage.

Table 5: Stake table of flow velocity at the boundary of the left bank

Starting stake Termination stake The maximum velocity in the section 0+000 0+500 0.18 0+500 1+000 0.17 1+000 1+500 0.15 … … … 30+000 30+500 0.71 30+500 31+000 0.67 31+000 31+500 0.68

Table 6: Stake table of flow velocity at the right bank boundary

Starting stake Termination stake The maximum velocity in the section 0+000 0+500 0.64 0+500 1+000 0.56 1+000 1+500 0.68 … … … 30+000 30+500 0.26 30+500 31+000 0.27 31+000 31+500 0.29

It is understood that the specific order or hierarchy of steps in the processes disclosed is an example of exemplary approaches. Based upon design preferences, it is understood that the specific order or hierarchy of steps in the processes may be rearranged without departing from the scope of the present disclosure. The accompanying method claims present elements of the various steps in a sample order, and are not meant to be limited to the specific order or hierarchy described.

In the foregoing Detailed Description, various features are grouped together in a single embodiment to simplify the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments of the subject matter require more features than are expressly recited in each claim. Rather, as the following claims reflect, the invention lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby expressly incorporated into the Detailed Description, with each claim standing on its own as a separate preferred embodiment of this invention.

The foregoing description of the disclosed embodiments was provided to enable any person skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may also be applied to other embodiments without departing from the spirit and scope of the disclosure. Thus, the present disclosure is not intended to be limited to the embodiments presented herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

The foregoing description includes illustrations of one or more embodiments. Of course, it is impossible to describe all possible combinations of components or methods to describe the above-mentioned embodiments, but those skilled in the art should recognize that various embodiments can be further combined and permuted. Accordingly, the embodiments described herein are intended to embrace all such alterations, modifications and variations that fall within the scope of the appended claims. Furthermore, to the extent that the term "comprises" is used in the specification or claims, the word is encompassed in a manner similar to the term "comprises" as interpreted when "comprises" is used as a link in the claims. Furthermore, any use of the term "or" in the specification of the claims is intended to mean a "non-exclusive or".

Those skilled in the art can also understand that various illustrative logical blocks, units, and steps listed in the embodiments of the present invention can be implemented by electronic hardware, computer software, or a combination of both. To clearly demonstrate the interchangeability of hardware and software, the various illustrative components, units and steps above have generally described their functions. Whether such functions are implemented by hardware or software depends on the specific application and overall system design requirements. Those skilled in the art may use various methods to implement the described functions for each specific application, but such implementation should not be understood as exceeding the protection scope of the embodiments of the present invention.

Various illustrative logic blocks or units described in the embodiments of the present invention can be discretely processed by a general-purpose processor, a digital signal processor, an application-specific integrated circuit (ASIC), a field programmable gate array or other programmable logic devices. Gate or transistor logic, discrete hardware components, or any combination of the above designed to implement or operate the described functions. The general-purpose processor may be a microprocessor, and optionally, the general-purpose processor may also be any conventional processor, controller, microcontroller or state machine. A processor may also be implemented by a combination of computing devices, such as a digital signal processor and a microprocessor, multiple microprocessors, one or more microprocessors combined with a digital signal processor core, or any other similar configuration to accomplish.

The steps of the method or algorithm described in the embodiments of the present invention may be directly embedded in hardware, a software module executed by a processor, or a combination of both. The software modules may be stored in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, removable disk, CD-ROM or any other storage medium in the art. Exemplarily, the storage medium can be connected to the processor, so that the processor can read information from the storage medium, and can write information to the storage medium. Optionally, the storage medium can also be integrated into the processor. The processor and the storage medium can be set in the ASIC, and the ASIC can be set in the user terminal. Optionally, the processor and the storage medium may also be set in different components in the user terminal.

In one or more exemplary designs, the above functions described in the embodiments of the present invention may be implemented in hardware, software, firmware or any combination of the three. If implemented in software, the functions can be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes computer storage media and communication media that facilitate transfer of a computer program from one place to another. Storage media may be any available media that can be accessed by a general purpose or special computer. For example, such computer-readable media may include, but are not limited to, RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other device that can be used to carry or store instructions or data structures and Other medium of program code in a form readable by a general-purpose or special-purpose computer, or a general-purpose or special-purpose processor. In addition, any connection is properly defined as a computer-readable medium, for example, if the software is transmitted from a website site, server, or other remote source via a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL) Or transmitted by wireless means such as infrared, wireless and microwave are also included in the definition of computer readable media. Disks and discs include compact discs, laser discs, optical discs, DVDs, floppy discs, and Blu-ray discs. Disks usually reproduce data magnetically, while discs usually reproduce data optically with lasers. Combinations of the above can also be contained on a computer readable medium.

The above-described embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still implement the foregoing embodiments Modifications to the technical solutions described in the examples, or equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the application, and should be included in the Within the protection scope of this application.

Claims (10)

1. A river bank ecological slope protection anti-scour velocity review method, is characterized in that, comprising: Use the surface water simulation system to generate a grid file based on the topographic map, and indicate the boundary conditions of the study area in the surface water simulation system, import it into MIKE21, determine the upstream and downstream boundaries, perform calculations in MIKE21, and obtain dfsu including flow velocity and water depth information document; Read the grid coordinates from the calculation result dfsu file, use the boundary type code to get the boundary grid coordinates from the calculation result dfsu file, and generate a CSV file; Import the generated coordinate CSV file into QGIS, add the left and right bank attributes in the attribute table, and export the generated attribute table as an xlsx file after sorting according to the left and right bank attributes; Read the terrain grid coordinates, match them with the left and right bank coordinates, obtain the index information of the boundary grid coordinates, read the velocity information of the calculation results, and use the index information of the boundary grid coordinates to obtain the boundary flow velocity information; Sort the coordinates of the boundary grid, calculate the cumulative length of the sorted coordinate points along the river bank axis direction, and correspond to the pile number section that needs to be laid with ecological slope protection, compare the maximum flow velocity on the bank of the corresponding pile number section with the ecological slope protection resistance Evaluate the applicability of ecological slope protection in this pile number section. 2. The anti-scour flow velocity review method of river bank ecological slope protection according to claim 1, wherein the boundary conditions of the study area include a no-slip boundary, a flow boundary and a water level boundary. 3. the riverbank ecological slope protection anti-scour flow velocity review method according to claim 1, is characterized in that, described imports in MIKE21, determines upstream and downstream boundary, calculates in MIKE21, obtains the dfsu file that comprises flow velocity, depth of water information and comprises: Create a new MIKE21 FM project in the MIKE ZERO software, import the grid file, and determine the upstream and downstream boundaries; Given the boundary conditions of water level and flow in MIKE21, calculate the results of flow field velocity and water level under different working conditions; Select the dfsu output, and check the flow rate and water level information in the output content. 4. The river bank ecological slope protection anti-scour flow velocity review method according to claim 3, characterized in that, the given water level and flow boundary conditions in MIKE21 include river bank, river channel inlet and outlet water level and flow conditions. 5. the river bank ecological slope protection anti-scour velocity review method according to claim 1, is characterized in that, described read grid coordinates from calculation result dfsu file and comprise utilizing the MATLAB mzReadMesh function that MIKE SDK provides from calculation result dfsu file Read grid coordinates. 6. the method for reviewing the anti-scour velocity of river bank ecological slope protection according to claim 1, is characterized in that, described in attribute table, increases left and right banks attribute and comprises: In QGIS, select the scatter points belonging to the left bank and mark the left bank in the attribute table, and select the scatter points belonging to the right bank to mark the right bank in the attribute table. 7. the river bank ecological slope protection anti-scour velocity review method according to claim 1, is characterized in that, described read topographic grid coordinates, it is matched with left and right bank coordinates, obtains the index information of boundary grid coordinates, reads Calculate the velocity information of the result, and use the index information of the boundary grid coordinates to obtain the boundary flow velocity information including: Read the boundary coordinate attribute table to obtain the coordinate information of the left and right banks; Read the coordinates in the grid file, match them with the coordinates of the left and right banks, and obtain the index information of the grid points at the boundary in the calculation result dfsu file; Use the index information to get the boundary flow velocity information from the calculation result dfsu file. 8. the river bank ecological slope protection anti-scour velocity review method according to claim 1, is characterized in that, sorting the boundary grid coordinates comprises: If the river flows from west to east, sort the coordinate points and corresponding boundary flow velocities according to the X coordinates of the left and right bank boundaries from small to large; If the river flows from east to west, the coordinate points and corresponding boundary flow velocities are sorted according to the order of the X coordinates of the left and right bank boundaries from large to small. 9. The river bank ecological slope protection anti-scour velocity review method according to claim 1 is characterized in that: if the maximum flow velocity on the bank of the corresponding stake number section is not greater than the ecological slope protection anti-scour velocity, then the ecological slope protection is applicable in this stake number section, If the maximum flow velocity on the bank of the corresponding pile number section is greater than the anti-scour velocity of the ecological slope protection, the ecological slope protection is not applicable to the pile number section. 10. The method for reviewing anti-scour velocity of river bank ecological slope protection according to claim 1, characterized in that: grid node coordinates are filled in the CSV file in the form of a table.

Images