CN115408832A - River bank ecological slope protection anti-impact flow velocity rechecking method - Google Patents
River bank ecological slope protection anti-impact flow velocity rechecking method Download PDFInfo
- Publication number
- CN115408832A CN115408832A CN202210966734.XA CN202210966734A CN115408832A CN 115408832 A CN115408832 A CN 115408832A CN 202210966734 A CN202210966734 A CN 202210966734A CN 115408832 A CN115408832 A CN 115408832A
- Authority
- CN
- China
- Prior art keywords
- boundary
- bank
- flow velocity
- grid
- coordinates
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 238000011160 research Methods 0.000 claims abstract description 9
- 238000012163 sequencing technique Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 10
- 238000011144 upstream manufacturing Methods 0.000 claims description 7
- 239000002352 surface water Substances 0.000 claims description 6
- 238000013461 design Methods 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000000605 extraction Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002787 reinforcement Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000007103 stamina Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
Technical Field
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.
Background
River shoreline is regarded as the important protective screen of flood control, the important component of ecosystem, and its system governance receives government's more and more attention. The new era has new requirements for the regulation of the shoreline, and the existing regulation of the shoreline is not limited by the traditional bank protection and protection engineering, but is transferred to the construction of a coastal river ecological landscape zone through a series of engineering measures. Compared with traditional hard revetments such as grouted stone revetments, dry stone revetments, concrete revetments and the like, ecological revetments such as flexible ecological water and soil protection blankets, ecological vegetation net mats or Reynolds protection mats and the like are increasingly adopted in the embankment slope protection engineering due to the ecological friendliness.
Compared with the traditional hard slope protection method, the ecological slope protection method has certain requirements on the near-shore flow velocity, so that the water flow velocity of the bank slope edge needs to be calculated by adopting a two-dimensional hydrodynamic model and is compared with the anti-impact flow velocity of the ecological slope protection, and the engineering form and the applicability of the ecological slope protection are determined. MIKE21 is a two-dimensional hydrodynamic numerical calculation software researched and developed by the research institute of water conservancy in denmark, and a two-dimensional shallow water equation is discretely solved by adopting a finite element method, so that water depth and flow velocity results of a two-dimensional flow field can be obtained after boundary conditions are determined. MIKE21 is widely applied in China, and MIKE21 is applied to typical great rivers in China to solve important engineering problems, such as the simulation of hydrodynamics and salinity at the mouth of a Yangtze river, the simulation of water flow and silt at the mouth of a Zhujiang river, the simulation of water flow and waves at a Bohai Bay and the like. It can be said that the process of calculating two-dimensional hydrodynamic force by MIKE21 has become one of the industry standards for hydrodynamic force calculation in the water conservancy industry.
MIKE21 was designed over the last century and provided a graphical interface that was old, difficult to use, and particularly cumbersome in extracting the distribution of river bank flow velocities.
Disclosure of Invention
The invention aims to provide a flow velocity resistant rechecking method for ecological river bank revetment based on the defects of the prior art, which is based on an MIKE21 program expansion package (MIKE SDK) provided by Danish water conservancy research institute and combines with spatial geographic analysis software (QGIS), so that the rapid extraction and analysis of the two-dimensional flow field bank slope flow velocity calculated by the MIKE21 are realized, and the design efficiency of ecological revetment is greatly improved.
The invention discloses an anti-impact flow velocity rechecking method for ecological river bank revetment, which comprises the following steps:
generating a grid file according to a topographic map by using a surface water simulation system, indicating boundary conditions of a research area in the surface water simulation system, leading the grid file into an MIKE21, determining upstream and downstream boundaries, and calculating in the MIKE21 to obtain a dfsu file containing flow speed and water depth information;
reading grid coordinates from the calculation result dfsu file, obtaining boundary grid coordinates from the calculation result dfsu file by using boundary type coding, and generating a CSV file;
importing the generated coordinate CSV file into a QGIS, adding left and right bank attributes in an attribute table, and exporting the generated attribute table as an xlsx file after sorting according to the left and right bank attributes;
reading the terrain grid coordinates, matching the terrain grid coordinates with the left and right bank coordinates to obtain index information of boundary grid coordinates, reading the speed information of a calculation result, and obtaining boundary flow speed information by using the index information of the boundary grid coordinates;
and sequencing the boundary grid coordinates, calculating the accumulated length of the sequenced coordinate points along the axial direction of the river bank, corresponding to the pile number section of the ecological slope protection to be laid, comparing the maximum flow velocity of the bank side of the corresponding pile number section with the impact flow velocity of the ecological slope protection, and evaluating whether the ecological slope protection is applicable to the pile number section.
Preferably, the boundary conditions of the research area include a non-slip boundary, a flow boundary and a water level boundary.
Preferably, the importing MIKE21 determines an upstream boundary and a downstream boundary, and the calculating is performed in the MIKE21 to obtain a dfsu file containing information of flow velocity and water depth, where the dfsu file includes:
newly building a MIKE21FM project in MIKE ZERO software, importing a grid file, and determining an upstream boundary and a downstream boundary;
setting a water level flow boundary condition in the MIKE21, and calculating flow field flow velocity water level results under different working conditions;
and selecting dfsu for output, and checking flow rate and water level information in output contents.
Preferably, the boundary conditions of the given water level flow in MIKE21 include water level flow conditions of a river bank and an entrance and exit of a river.
Preferably, the reading the grid coordinates from the calculation result dfsu file includes reading the grid coordinates from the calculation result dfsu file by using a MATLAB mzReadMesh function provided by MIKE SDK.
Preferably, the adding left and right bank attributes in the attribute table includes:
and selecting scatter points belonging to a left bank from the QGIS to mark the left bank in the attribute table, and selecting scatter points belonging to a right bank to mark the right bank in the attribute table.
Preferably, the reading the terrain grid coordinate, matching the terrain grid coordinate with the left and right bank coordinates to obtain index information of the boundary grid coordinate, reading the speed information of the calculation result, and obtaining the boundary flow speed information by using the index information of the boundary grid coordinate includes:
reading a boundary coordinate attribute table to obtain coordinate information of a left bank and a right bank;
reading coordinates in the grid file, and matching the coordinates with the left and right bank coordinates to obtain index information of grid points positioned at the boundary in the calculation result dfsu file;
and obtaining boundary flow velocity information from the calculation result dfsu file by using the index information.
Preferably, the sorting the boundary grid coordinates includes:
if the river flow direction is from west to east, sorting coordinate points and corresponding boundary flow velocity according to the sequence that the X coordinate of the left and right bank boundaries is from small to large;
and if the river flows from east to west, sorting the coordinate points and the corresponding boundary flow velocity according to the sequence of the X coordinates of the left and right bank boundaries from large to small.
It is comparatively preferred, if the maximum velocity of flow of corresponding pile number section bank is not more than the anti speed of flow of dashing of ecological bank protection, then ecological bank protection is suitable for in this pile number section, if the maximum velocity of flow of corresponding pile number section bank is greater than the anti speed of flow of dashing of ecological bank protection, then ecological bank protection is not suitable for in this pile number section.
Preferably, the CSV file is filled with grid node coordinates in a form.
The invention has the beneficial effects that:
1. according to the method, the MIKE21 program expansion package MIKE SDK provided by Danish water conservancy research is combined with the QGIS (spatial geographic analysis software), so that the rapid extraction and analysis of the two-dimensional flow field bank slope flow velocity calculated by the MIKE21 are realized, and the design efficiency of ecological bank protection is greatly improved.
2. The method omits a complex post-processing method in the traditional bank flow velocity calculation, directly extracts the bank flow velocity by utilizing MATLAB API provided by MIKE SDK, and directly corresponds to the pile number in the engineering design practice, thereby greatly facilitating the rechecking of the anti-impact flow velocity of the ecological slope protection.
3. By specially encoding the boundary nodes during the construction of the grid, the boundary nodes are distinguished from the nodes in the river, so that post-processing is facilitated.
4. The boundary grid coordinates are led into the QGIS for post-processing, so that left and right bank parts in the boundary grid coordinates are effectively distinguished, and the flow velocity of the left bank and the right bank is conveniently and respectively calculated.
5. And the boundary grid coordinates are post-processed through a QGIS (Quadrature information System), so that the continuous coding of boundary grid nodes along a river levee axis is ensured, and the calculation of the boundary flow velocity pile number is facilitated.
Drawings
FIG. 1 is a schematic flow diagram of the present invention;
FIG. 2 is a topographical grid map of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantageous effects to be solved by the present application clearer, the present application is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of 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 that depart from 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 will be understood that the terms "comprises" and/or "comprising," when used in this specification and the appended claims, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It should also be understood that the term "and/or" as used in this specification and the appended claims refers to any and all possible combinations of one or more of the associated listed items and includes such combinations.
As used in this specification and the appended claims, the term "if" may be interpreted contextually as "when", "upon" or "in response to" determining "or" in response to detecting ". Similarly, the phrase "if it is determined" or "if a [ described condition or event ] is detected" may be interpreted contextually to mean "upon determining" or "in response to determining" or "upon detecting [ described condition or event ]" or "in response to detecting [ described condition or event ]".
Furthermore, in the description of the present application and the appended claims, the terms "first," "second," "third," and the like are used for distinguishing between descriptions and not necessarily for describing or implying relative importance.
Reference throughout this specification to "one embodiment" or "some embodiments," or the like, 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," or the like, in various places throughout this specification are not necessarily all referring to the same embodiment, but rather mean "one or more but not all embodiments" unless specifically stated otherwise. The terms "comprising," "including," "having," and variations thereof mean "including, but not limited to," unless expressly specified otherwise. "plurality" means "two or more".
Example one
Fig. 1 shows a schematic structural diagram of an impact-resistant flow velocity review method for ecological slope protection of river bank provided in a preferred embodiment of the present application (fig. 1 shows a first embodiment of the present application), and for convenience of description, only the parts related to the present embodiment are shown, and detailed descriptions are as follows:
the invention discloses an anti-impact flow velocity rechecking method for ecological river bank revetment, which comprises the following steps:
step 1: generating a grid file according to a topographic map by using a Surface water simulation System (SMS), indicating boundary conditions of a research area in the Surface water simulation System, including a non-slip boundary, a flow boundary and a water level boundary, newly building a MIKE21FM project in MIKE ZERO software, importing the grid file, determining an upstream boundary and a downstream boundary, setting a water level flow boundary condition in the MIKE21, and calculating to obtain a dfsu file containing flow speed and water depth information. Calculating flow rate and water level results of a flow field under different working conditions, selecting dfsu for output, and checking information such as flow rate, water level and the like in output contents;
step 2: reading grid coordinates from the calculation result dfsu file by using an MATLAB mzReadMesh function provided by the MIKE SDK, obtaining boundary grid coordinates from the calculation result dfsu file by using boundary type codes (the code of the river channel boundary in the grid file of the MIKE21 is 1), and generating a CSV file; the format of the CSV file is shown in table 1 (where the grid node coordinates are filled out, see table 3 for a specific example).
Table 1: boundary coordinate format table
| X | Y |
| … | … |
And 3, step 3: and importing the generated coordinate CSV file into a QGIS, and adding left and right bank attributes into an attribute table, wherein the format of the attribute table is shown in a table 2. And selecting scatter points belonging to the left bank from the QGIS to mark the left bank (0) in the attribute table, and then selecting scatter points belonging to the right bank to mark the right bank (1) in the attribute table. And sorting the generated attribute table according to the attributes of the left bank and the right bank, and respectively integrating the scatter points positioned on the left bank and the right bank together and exporting the scatter points into an xlsx file (format).
Table 2: boundary coordinate attribute table
| X | Y | Left and right bank |
| … | … | 0/1 |
*0 represents the left bank, 1 represents the right bank
And 4, step 4: reading a terrain grid coordinate, matching the terrain grid coordinate with left and right bank coordinates to obtain index information of a boundary grid coordinate, reading speed information of a calculation result, and obtaining boundary flow speed information by using the index information of the boundary grid coordinate; the method specifically comprises the following steps:
step 401: reading a result file calculated by MIKE21 by using a DFSU MATLAB API provided by the MIKE SDK, firstly checking the storage position of speed information (Current speed) in the DFSU file, and then reading the speed information in the calculation result.
Step 402: and reading the boundary coordinate attribute table by using MATLAB to obtain the coordinate information of the left bank and the right bank. And reading the coordinates in the grid file by using an MATLAB mzReadMesh function provided by the MIKE SDK, and matching the coordinates with the left and right bank coordinates to obtain the index information of the grid points positioned at the boundary in the calculation result dfsu file. And obtaining boundary flow velocity information from a calculation result dfsu file by using the index information.
And 5: and sequencing the boundary grid coordinates, calculating the accumulated length of the sequenced coordinate points along the axial direction of the river bank, corresponding to the pile number section of the ecological slope protection to be laid, comparing the maximum flow velocity of the bank side of the corresponding pile number section with the impact flow velocity of the ecological slope protection, and evaluating whether the ecological slope protection is applicable to the pile number section. If the maximum flow velocity on the corresponding pile number section bank side is not more than the anti-impact flow velocity of the ecological slope protection, the ecological slope protection is applicable to the pile number section, and if the maximum flow velocity on the corresponding pile number section bank side is more than the anti-impact flow velocity of the ecological slope protection, the ecological slope protection is not applicable to the pile number section.
The method for sorting the boundary grid coordinates comprises the following steps:
if the river flows from west to east, sorting the coordinate points and the corresponding boundary flow velocity according to the sequence that the X coordinate of the left and right bank boundaries is from small to large; and if the river flows from east to west, sorting the coordinate points and the corresponding boundary flow velocity according to the sequence of the X coordinates of the left and right bank boundaries from large to small.
Example two
The invention takes the anti-impact rechecking calculation of the ecological slope protection in the reinforcement and treatment engineering of the dyke of the new covered house flood diversion tunnel in the new male security area as an embodiment for detailed description, and has guiding significance for the anti-impact rechecking calculation of other ecological slope protection.
Step 1: generating a grid file from landform scattered points by using SMS (short message service), as shown in figure 2, newly building a MIKE21FM project in MIKE ZERO software, importing the grid file, determining boundary conditions (river bank and river channel inlet and outlet water level flow conditions), calculating flow rate and water level results of a flow field under different working conditions, selecting dfsu for output, and checking information such as flow rate, water level and the like in output contents.
Step 2: and reading the coordinate grid boundary coordinates by using a MATLAB mzReadMesh function provided by the MIKE SDK, and generating a CSV file. The coordinate grid boundary coordinates CSV file 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 |
And 3, step 3: and importing the generated coordinate CSV file into a QGIS, and adding left and right bank attributes into an attribute table. And selecting scatter points belonging to a left bank from the QGIS to mark the left bank (0) in the attribute table, and then selecting scatter points belonging to a right bank to mark the right bank (1) in the attribute table. And sorting the generated attribute table according to the attributes of the left bank and the right bank and exporting the attribute table as an xlsx file. The coordinate attribute file is shown in table 4.
Table 4: boundary left and right bank attribute table
| X | Y | Left and right bank |
| 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: and reading a result file calculated by MIKE21 by using a DFSU MATLAB API, wherein speed information (Current speed) in the model simulation is positioned at the 5 th of an output data set in a DFSU file, and reading all speed information.
Step 402: and reading the xlsx file of the attribute table by using MATLAB to obtain coordinate files of the left bank and the right bank. And reading the grid coordinates of the calculation result dfsu file by using an MATLAB mzReadMesh function provided by the MIKE SDK, and matching the grid coordinates with the left and right bank coordinates by using the coordinate position to obtain the index information of the boundary coordinate points in the grid of the calculation result dfsu file. And obtaining boundary flow velocity information by using the index information.
And 5: the new covered room flood diversion tunnel in the new male security area flows from west to east, so that the boundary coordinate points and the corresponding boundary flow rates are sequenced according to the sequence that the X coordinates of the left and right bank boundaries are from small to large. And calculating the accumulated length of the sorted riverbank coordinate points along the riverbank axis direction. This time new covering room flood diversion tunnel dyke reinforcement of new area of peace of stamina and improvement engineering need be at the whole section ecological bank protection of application, so 500 meters are the biggest velocity of flow in an interval statistics interval, and the speed of flow is resisted to the contrast ecological bank protection, and whether the evaluation ecological bank protection is suitable for in this stake number section. The accumulated length of the river bank coordinate points is the pile number corresponding to the grid point, the maximum flow velocity in the pile number is counted by taking the interval of 500 meters from the head of the river bank as an interval, and tables of the maximum flow velocity in the pile number are shown in tables 5 and 6.
Table 5: left bank boundary flow velocity pile number meter
| Initial pile number | Number of terminal pile | Maximum flow velocity in the segment |
| 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: right bank boundary flow velocity pile number table
| Initial pile number | Number of terminal pile | Maximum flow velocity in the segment |
| 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 should be 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 presented.
In the foregoing detailed description, various features are grouped together in a single embodiment for the purpose of streamlining 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, 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 the invention.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. To those skilled in the art; various modifications to these embodiments will be readily apparent, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the disclosure. Thus, the present disclosure is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
What has been described above includes examples of one or more embodiments. It is, of course, not possible to describe every conceivable combination of components or methodologies for purposes of describing the aforementioned embodiments, but one of ordinary skill in the art may recognize that many further combinations and permutations of various embodiments are possible. 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 "includes" is used in either the detailed description or the claims, such term is intended to be inclusive in a manner similar to the term "comprising" as "comprising" is interpreted when employed as a transitional word in a claim. Furthermore, any use of the term "or" in the specification of the claims is intended to mean a "non-exclusive or".
Those of skill in the art will further appreciate that the various illustrative logical blocks, units, and steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate the interchangeability of hardware and software, various illustrative components, elements, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design requirements of the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present embodiments.
The various illustrative logical blocks or elements described in connection with the embodiments disclosed herein may be implemented or performed with a general purpose processor, a digital signal processor, an Application Specific Integrated Circuit (ASIC), a field programmable gate array or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described herein. A general-purpose processor may be a microprocessor, but, in the alternative, the processor may be any conventional processor, controller, microcontroller, or state machine. A processor may also be implemented as a combination of computing devices, e.g., a digital signal processor and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a digital signal processor core, or any other similar configuration.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may be stored in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. For example, a storage medium may be coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC, which may be located in a user terminal. In the alternative, the processor and the storage medium may reside in different components in a user terminal.
In one or more exemplary designs, the 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 may be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both 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 purpose computer. For example, such computer-readable media can include, but is not limited to, RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to carry or store program code in the form of instructions or data structures and which can be read by a general-purpose or special-purpose computer, or a general-purpose or special-purpose processor. In addition, any connection is properly termed a computer-readable medium, and thus is included if the software is transmitted from a website, server, or other remote source via a coaxial cable, fiber optic cable, twisted pair, digital Subscriber Line (DSL), or wirelessly, e.g., infrared, radio, and microwave. Such discs (disk) and disks (disc) include compact disks, laser disks, optical disks, DVDs, floppy disks and blu-ray disks where disks usually reproduce data magnetically, while disks usually reproduce data optically with lasers. Combinations of the above may also be included in the computer-readable medium.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the embodiments of the present application, and they should be construed as being included in the present application.
Claims (10)
1. An anti-impact flow velocity rechecking method for ecological slope protection of river banks is characterized by comprising the following steps:
generating a grid file according to a topographic map by using a surface water simulation system, indicating boundary conditions of a research area in the surface water simulation system, leading the grid file into an MIKE21, determining upstream and downstream boundaries, and calculating in the MIKE21 to obtain a dfsu file containing flow speed and water depth information;
reading grid coordinates from the calculation result dfsu file, obtaining boundary grid coordinates from the calculation result dfsu file by using boundary type coding, and generating a CSV file;
importing the generated CSV file into a QGIS, adding left and right bank attributes in an attribute table, and exporting the generated attribute table into an xlsx file after sorting according to the left and right bank attributes;
reading the terrain grid coordinates, matching the terrain grid coordinates with the left and right bank coordinates to obtain index information of boundary grid coordinates, reading the speed information of a calculation result, and obtaining boundary flow speed information by using the index information of the boundary grid coordinates;
and sequencing the boundary grid coordinates, calculating the accumulated length of the sequenced coordinate points along the axis direction of the river bank, corresponding to the pile number section of the ecological protection slope to be laid, comparing the maximum flow velocity of the bank side of the corresponding pile number section with the impact flow velocity of the ecological protection slope, and evaluating whether the ecological protection slope is applicable to the pile number section.
2. The method for rechecking the flow resistance of the ecological slope protection of the river bank according to claim 1, wherein the boundary conditions of the research area include a non-slip boundary, a flow boundary and a water level boundary.
3. The method for rechecking the flow velocity for impact resistance of the ecological slope protection of the river bank according to claim 1, wherein the step of guiding into MIKE21, determining the upstream and downstream boundaries, and calculating in MIKE21 to obtain a dfsu file containing information of flow velocity and water depth comprises the following steps:
newly building a MIKE21FM project in MIKE ZERO software, importing a grid file, and determining an upstream boundary and a downstream boundary;
setting a water level flow boundary condition in the MIKE21, and calculating flow field flow velocity water level results under different working conditions;
and selecting dfsu output, and checking flow rate and water level information in the output content.
4. The method for rechecking the flow resistance of the ecological slope protection of the river bank according to claim 3, wherein the boundary conditions of the given water level flow in the MIKE21 include water level flow conditions of the river bank and the river channel inlet and outlet.
5. The method for rechecking the impact flow rate of the ecological river bank revetment according to claim 1, wherein the reading of the grid coordinates from the calculation result dfsu file comprises reading the grid coordinates from the calculation result dfsu file using a MATLAB mzReadMesh function provided by MIKE SDK.
6. The method for rechecking the impact flow rate of the ecological slope protection of the river bank as claimed in claim 1, wherein the adding left and right bank attributes in the attribute table comprises:
and selecting scattered points belonging to the left bank from the QGIS to mark the left bank in the attribute table, and selecting scattered points belonging to the right bank to mark the right bank in the attribute table.
7. The method for rechecking the flow velocity of the river bank ecological slope protection impact resistance of the claim 1, wherein the reading of the terrain grid coordinates, the matching of the terrain grid coordinates with the left and right bank coordinates to obtain the index information of the boundary grid coordinates, the reading of the speed information of the calculation result, and the obtaining of the boundary flow velocity information by using the index information of the boundary grid coordinates comprise:
reading a boundary coordinate attribute table to obtain coordinate information of a left bank and a right bank;
reading coordinates in the grid file, and matching the coordinates with the left and right bank coordinates to obtain index information of grid points positioned at the boundary in the dfsu file;
and obtaining boundary flow velocity information from the calculation result dfsu file by using the index information.
8. The method for rechecking the flow speed of the river bank ecological slope protection impact resistance of the river bank as claimed in claim 1, wherein the sorting of the boundary grid coordinates comprises:
if the river flows from west to east, sorting the coordinate points and the corresponding boundary flow velocity according to the sequence that the X coordinate of the left and right bank boundaries is from small to large;
and if the river flows from east to west, sorting the coordinate points and the corresponding boundary flow velocity according to the sequence of the X coordinates of the left and right bank boundaries from large to small.
9. The method for rechecking the flow speed of the river bank ecological slope protection impact resistance of the river bank according to claim 1, wherein the method comprises the following steps: if the maximum flow velocity of corresponding stake number section bank is not more than the anti-impact flow velocity of ecological bank protection, then ecological bank protection is suitable for at this stake number section, if the maximum flow velocity of corresponding stake number section bank is greater than the anti-impact flow velocity of ecological bank protection, then ecological bank protection is not suitable for at this stake number section.
10. The method for rechecking the flow speed of the river bank ecological slope protection impact resistance of the river bank according to claim 1, wherein the method comprises the following steps: and grid node coordinates are filled in the CSV file in a form.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210966734.XA CN115408832B (en) | 2022-08-12 | 2022-08-12 | River bank ecological slope protection impact-resistant flow velocity rechecking method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210966734.XA CN115408832B (en) | 2022-08-12 | 2022-08-12 | River bank ecological slope protection impact-resistant flow velocity rechecking method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115408832A true CN115408832A (en) | 2022-11-29 |
| CN115408832B CN115408832B (en) | 2023-12-01 |
Family
ID=84159622
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210966734.XA Active CN115408832B (en) | 2022-08-12 | 2022-08-12 | River bank ecological slope protection impact-resistant flow velocity rechecking method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115408832B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120101795A1 (en) * | 2010-10-26 | 2012-04-26 | The Government Of The United States, As Represented By The Secretary Of The | Tool for rapid configuration of a river model using imagery-based information |
| CN206815280U (en) * | 2016-12-21 | 2017-12-29 | 中国水利水电科学研究院 | A kind of ecological revetment based on tortuous channel |
| CN111898303A (en) * | 2020-08-05 | 2020-11-06 | 苏州大圜科技有限公司 | River basin water level and waterlogging forecasting method based on weather forecasting and hydrodynamic simulation |
| CN113642057A (en) * | 2021-07-13 | 2021-11-12 | 长江勘测规划设计研究有限责任公司 | Method for automatically generating flow velocity vector diagram in Autocad format |
| US20220074157A1 (en) * | 2019-01-05 | 2022-03-10 | Fudan University | Novel ecological slope protection with efficient water purification function and ecological improvment method |
-
2022
- 2022-08-12 CN CN202210966734.XA patent/CN115408832B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120101795A1 (en) * | 2010-10-26 | 2012-04-26 | The Government Of The United States, As Represented By The Secretary Of The | Tool for rapid configuration of a river model using imagery-based information |
| CN206815280U (en) * | 2016-12-21 | 2017-12-29 | 中国水利水电科学研究院 | A kind of ecological revetment based on tortuous channel |
| US20220074157A1 (en) * | 2019-01-05 | 2022-03-10 | Fudan University | Novel ecological slope protection with efficient water purification function and ecological improvment method |
| CN111898303A (en) * | 2020-08-05 | 2020-11-06 | 苏州大圜科技有限公司 | River basin water level and waterlogging forecasting method based on weather forecasting and hydrodynamic simulation |
| CN113642057A (en) * | 2021-07-13 | 2021-11-12 | 长江勘测规划设计研究有限责任公司 | Method for automatically generating flow velocity vector diagram in Autocad format |
Non-Patent Citations (7)
| Title |
|---|
| 刘启翔: "清泥河流域综合治理工程护坡防护形式选择", 陕西水利, pages 2 - 3 * |
| 张志林 等: "基于MIKE21 FM模型的河道流场图绘制", 东北水利水电, pages 2 - 5 * |
| 李杰;: "某中小河流岸坡整治及护岸工程选型分析", 广东水利水电, no. 01 * |
| 罗日洪;黄锦林;王立华;张令;张志伟;: "中小河流治理中雷诺护垫设计方法探讨", 广东水利水电, no. 01 * |
| 艾来提江・谢热甫丁;: "西白杨沟综合治理工程水面线推求及护坡方案比选", 陕西水利, no. 09 * |
| 赵通阳;符强;邰肇悦;: "城市防洪河道生态护岸形式选择与设计", 黄河水利职业技术学院学报, no. 01 * |
| 韩志聪;樊彦国;李祥昌;: "基于QGIS的黄河三角洲岸线演变预测***设计与实现", 测绘与空间地理信息, no. 06 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115408832B (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107239607B (en) | Model element and calculation scheme management method | |
| CN108984962B (en) | Large-particle damming dam burst flow calculation method under impact of burst type viscous debris flow | |
| CN110008608B (en) | Calculation method and application of earthquake-region burst type debris flow | |
| Mehta et al. | Analysis of scour depth in the case of parallel bridges using HEC-RAS | |
| Zolghadr et al. | Migration of sand mining pit in rivers: an experimental, numerical and case study | |
| CN102864754B (en) | The investigating method of a kind of bridge location section trraning pole scour depth | |
| Donald et al. | Hydraulic method to evaluate the performance of ditch check practices and products | |
| CN113221335B (en) | Wave height calculation method and device after landslide surge attenuation in curved river channel and storage medium thereof | |
| CN115408832B (en) | River bank ecological slope protection impact-resistant flow velocity rechecking method | |
| Norbury et al. | Combining hydraulic modelling with partnership working: towards practical natural flood management | |
| Fenocchi et al. | Using numerical and physical modeling to evaluate total scour at bridge piers | |
| Zhang et al. | Influence of bridge piers shapes on the flow of the lower Yellow River | |
| Liu et al. | Urban flood modelling in Qiqihar city based on MIKE flood | |
| CN109615195B (en) | Grading evaluation method for hydrological landform of mountain river | |
| Zhang et al. | Advances in Water Resources & Hydraulic Engineering: Proceedings of 16th IAHR-APD Congress and 3rd Symposium of IAHR-ISHS | |
| Ahmed | Influence of wetlands on black-creek hydraulics | |
| Moghadam et al. | Numerical study of energy dissipation and block barriers in stepped spillways | |
| Donald et al. | Comparative evaluation of wattle ditch checks composed of differing materials and properties | |
| CN115310322B (en) | Ecological base stream calculation method, system, device and storage medium | |
| Mo et al. | Simulation of one-dimensional dam-break flood routing based on HEC-RAS | |
| Vaughn et al. | Floating woody debris: Blocking sensitivity of labyrinth weirs in channel and reservoir applications | |
| Islam et al. | Impact analysis of Sand dredging from alluvial tidal river | |
| Fitri et al. | Determining the Effectiveness of Harapan Lake as Flood Retention Pond in Flood Mitigation Effort | |
| Zhang et al. | Study on the effect of morphological changes of bridge piers on water movement properties | |
| CN113312690A (en) | Estuary leaf vein-like river network tide boundary determination method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |