River ecology assessment method and device based on hydrodynamic force, medium and product

Abstract

The invention discloses a river ecology assessment method and device based on hydrodynamic force, a medium and a product, and relates to the field of ecological hydraulics, and the method comprises the steps: obtaining basic feature elements and historical flow velocity data of a river cross section; obtaining time sequence data of the flow velocity of the surface of the river by using non-contact measurement equipment; determining entropy parameters of the river cross section according to the historical flow velocity data; obtaining two-dimensional flow velocity field data of the cross section of the river channel by using an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameters and the time series data; according to the two-dimensional flow velocity field data, calculating two hydraulic complexity indexes including a kinetic energy gradient parameter and a normalized energy change rate; obtaining a two-dimensional distribution field of the hydraulic complexity index by using a two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; and performing quantitative evaluation on the river ecosystem according to the two-dimensional distribution field. According to the invention, non-contact measurement can be fully utilized, and rapid and accurate evaluation of the river ecosystem can be realized.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of ecological hydraulics, and in particular to a river ecological evaluation method based on water power, equipment, medium and product. BACKGROUND

[0002] River ecosystems are an important part of global biodiversity, and bear multiple key functions such as regulating climate, purifying water quality, maintaining species diversity, and providing ecological services. However, the exertion of its ecological functions is largely influenced by the significant effects of river hydrodynamic processes. Specifically, the flow structure in rivers, such as velocity gradient and secondary flow, provides diverse habitats for aquatic organisms, and the heterogeneity of these habitats is of great significance to the survival, reproduction and distribution of organisms. Considering that aquatic habitats are closely related to hydrodynamic characteristics (such as velocity gradient), the hydrodynamic complexity indicators (kinetic energy gradient parameter (M1) and normalized energy change rate (M2)) are proposed from the perspective of river hydrodynamic structure to analyze the characteristics of aquatic habitats.

[0003] However, traditional river ecosystem monitoring methods mostly rely on field sampling surveys, which are not only time-consuming and labor-intensive, but also difficult to efficiently and accurately evaluate the hydrodynamic characteristics of biological habitats, and there is a certain risk in the measurement process. Especially in large rivers, field measurement has high risk and difficulty, which makes the evaluation of river ecology more challenging.

[0004] Therefore, there is an urgent need for a low-interference, efficient and convenient monitoring and evaluation method for rivers to achieve rapid and accurate evaluation of river ecosystems. SUMMARY

[0005] The purpose of the present application is to provide a river ecological evaluation method based on water power, equipment, medium and product, which can realize rapid and accurate evaluation of river ecosystems.

[0006] To achieve the above purpose, the present application provides the following solutions: In a first aspect, the present application provides a river ecological evaluation method based on water power, which comprises: obtaining basic characteristic elements of a river cross section and corresponding historical flow rate data; the basic characteristic elements include: measurement data of riverbank surface topography and underwater topography, measurement point coordinates and corresponding relative elevations obtained by a certain distance in the river cross section; the historical flow rate data includes: average flow rate and maximum flow rate of the river cross section; obtaining surface flow rate of the river cross section and corresponding time series data by using a non-contact measurement device; the non-contact measurement device includes: a radar speed measuring instrument or a large particle image velocimeter; determine an entropy parameter of the river cross section according to the historical flow rate data; and obtain two-dimensional flow rate field data of the river cross section based on an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and time sequence data corresponding to the surface flow rate; determine a hydraulic complexity index according to the two-dimensional flow rate field data of the river cross section; and obtain a two-dimensional distribution field of the hydraulic complexity index of the river cross section by using a two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; the hydraulic complexity index includes a kinetic energy gradient parameter and a kinetic energy change parameter; quantitatively evaluate the river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0007] Optionally, the obtaining of the basic characteristic elements of the river cross section and the corresponding historical flow rate data specifically includes: measure the underwater topography at a set distance on the river cross section by using an ultrasonic detector or a boat-mounted Doppler profile flowmeter to obtain measurement data of the underwater topography; measure the relative elevations and distances of the two banks along the direction of the river cross section to fixed pillars on the two banks by using a level and a range finder with a relative elevation zero coordinate point to obtain measurement data of the bank surface topography; the relative elevation zero coordinate point is a fixed pillar on one side of the river cross section; convert the relative elevations by using the relative elevation zero coordinate point and determine the three-dimensional cross-sectional topography of the river cross section and topographic measurement points according to the measurement data of the underwater topography and the measurement data of the bank surface topography; obtain the corresponding historical flow rate data according to the three-dimensional cross-sectional topography of the river cross section and the topographic measurement points.

[0008] Optionally, the determining of the entropy parameter of the river cross section according to the historical flow rate data and the obtaining of the two-dimensional flow rate field data of the river cross section based on the entropy model and the two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and the time sequence data corresponding to the surface flow rate specifically include: determine an entropy function by using a formula determine the entropy function; wherein, the entropy function, M is the entropy parameter, U m the average flow rate of the river cross section, U max the maximum flow rate of the river cross section, and e is a natural logarithm; determine the entropy parameter of the river cross section according to the entropy function; determine the cross-sectional point flow rate of the river cross section based on the entropy model according to the entropy parameter of the river cross section and the time sequence data corresponding to the surface flow rate; obtain the two-dimensional flow rate field data of the river cross section by using a two-dimensional inverse distance weighted interpolation algorithm according to the cross-sectional point flow rate.

[0009] Optionally, a hydraulic complexity index is determined based on the two-dimensional velocity field data of the river cross section, specifically including: Using the formula Determine the kinetic energy gradient parameter M1; Using the formula Determine the normalized energy change rate M2; Where V1 and V2 are the flow velocities at adjacent measuring points, Δs is the distance between measuring points, and V ave is the average velocity of adjacent points, V min is the smaller flow velocity value among adjacent points.

[0010] Optionally, the quantitative assessment of the river ecosystem based on the two-dimensional distribution field of the hydraulic complexity index of the river cross section specifically includes: The average value of each survey line is determined based on the two-dimensional distribution field of the hydraulic complexity index of the river cross section; and the plane distribution value of the hydraulic complexity index is obtained using the two-dimensional inverse distance weighted interpolation algorithm based on the average value of the survey line; Determine the average value of the hydraulic complexity index according to the plane distribution value of the hydraulic complexity index; The hydraulic complexity index of the river cross section is compared with the average value of the hydraulic complexity index to obtain the comparison result; Based on the comparison results, the assessment results of the river ecosystem are determined.

[0011] Optionally, determining the average value of the hydraulic complexity index according to the plane distribution value of the hydraulic complexity index specifically includes: Using the formula m=∑M1(x i ,y i ) / k determines the measuring point (x i ,y i ) the average value m of the kinetic energy gradient parameter M1; Using the formula n=∑M2(x i ,y i ) / k determines the measuring point (x i ,y i ) The average value n of the normalized energy change rate M2; Where k is the total number of measurement points within the river plane.

[0012] In a second aspect, the present application provides a river ecology assessment device based on hydrodynamics, the river ecology assessment device based on hydrodynamics comprising: The data acquisition module is configured to acquire basic characteristic elements of the river cross section and corresponding historical flow velocity data; the basic characteristic elements include: measurement data of riverbank surface topography and underwater topography, measurement point coordinates and corresponding relative elevations acquired at a certain distance on the river cross section; and the historical flow velocity data includes: average flow velocity and maximum flow velocity of the river cross section; The real-time data acquisition module is configured to acquire surface flow velocity of the river cross section and corresponding time sequence data by using a non-contact measurement device; the non-contact measurement device includes: a radar velocity measuring instrument or a large-scale particle image velocimetry; The two-dimensional flow velocity field data determination module is configured to determine an entropy parameter of the river cross section according to the historical flow velocity data; and obtain two-dimensional flow velocity field data of the river cross section based on an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and the time sequence data corresponding to the surface flow velocity; The hydraulic complexity index determination module is configured to determine a hydraulic complexity index according to the two-dimensional flow velocity field data of the river cross section; and obtain a two-dimensional distribution field of the hydraulic complexity index of the river cross section by using a two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; the hydraulic complexity index includes: a kinetic energy gradient parameter and a kinetic energy change parameter; The evaluation module is configured to quantitatively evaluate the river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0013] In a third aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to implement the water-power-based river ecosystem evaluation method.

[0014] In a fourth aspect, the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the water-power-based river ecosystem evaluation method.

[0015] In a fifth aspect, the present application provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to implement the water-power-based river ecosystem evaluation method.

[0016] According to the specific embodiments provided by the present application, the present application has the following technical effects: The application provides a river ecological evaluation method, device, medium and product based on water power, which fully embodies the efficiency, safety and ecological friendliness of non-contact measurement by using a non-contact measurement device to obtain surface flow velocity and corresponding time sequence data of a river cross section, avoids interference of traditional contact measurement on a river ecological system, has high automation, strong complex environment adaptability, high measurement safety and the like, can guarantee measurement work in a high flood flow period, and avoids safety hazards to hydrological measurement personnel, realizes automatic and rapid river ecological system evaluation by combining an entropy model and a water power complexity index to quantitatively evaluate heterogeneity of a river ecological system (a river biological habitat), provides a scientific basis for river ecological protection and biodiversity research, solves problems of great difficulty, high cost and difficulty in real-time online monitoring and evaluation of river ecological evaluation, and helps to improve river fine control ability and basin management and administration level. BRIEF DESCRIPTION OF DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments will be briefly introduced below. Obviously, the drawings in the following description only constitute some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0018] Fig. 1 A flowchart of a river ecological evaluation method based on water power in an embodiment of the present application; Fig. 2 A principle diagram of a river ecological evaluation method based on water power in an embodiment of the present application; Fig. 3 A technical route flowchart of a river ecological evaluation method based on water power in an embodiment of the present application; Fig. 4 M1 and M2 cross section distribution diagrams when the river flow is 124578 m 3 / s; Fig. 5 M1 and M2 cross section distribution diagrams when the river flow is 151663 m 3 / s; Fig. 6 M1 and M2 cross section distribution diagrams when the river flow is 154474 m 3 / s. DETAILED DESCRIPTION

[0019] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative effort belong to the scope of protection of the present application.

[0020] The above purposes, features and advantages of the present application will be more apparent and understandable. The present application will be further described in detail below with reference to the drawings and specific embodiments.

[0021] In an exemplary embodiment, as shown in Figs. 1-3 A water-based river ecological assessment method is provided, which comprises S101 to S105. Wherein: S101, obtaining basic characteristic elements of a river cross section and corresponding historical flow rate data; the basic characteristic elements include: measurement data of riverbank surface topography and underwater topography obtained at a certain distance on the river cross section, measurement point coordinates and corresponding relative elevations; the historical flow rate data includes: average flow rate and maximum flow rate of the river cross section; S101 specifically includes: S11, using an ultrasonic detector or a walking Doppler profile flowmeter to measure the underwater topography on the river cross section at a set distance to obtain the measurement data of the underwater topography; S12, using a level and a range finder to measure the relative elevations and distances of both sides of the riverbank along the direction of the river cross section to the fixed column on both banks with a relative elevation zero coordinate point to obtain the measurement data of the riverbank surface topography; the relative elevation zero coordinate point is a fixed column on one side of the river cross section; S13, converting the relative elevations by using the relative elevation zero coordinate point and determining the three-dimensional cross section topography of the river cross section and the topographic measurement point position according to the measurement data of the underwater topography and the measurement data of the riverbank surface topography; S14, obtaining the corresponding historical flow rate data according to the three-dimensional cross section topography of the river cross section and the topographic measurement point position.

[0022] S102, using a non-contact measurement device to obtain the surface flow rate of the river cross section and corresponding time sequence data; the non-contact measurement device includes but is not limited to: a radar speed measuring instrument or a large particle image velocimeter; Wherein, the selection of the non-contact measurement device is to measure the surface flow rate and water level of the river. An automatic control device is used to drive the radar speed measuring instrument to measure the surface flow rate and water level at different positions of the river cross section at a certain distance interval, and record the position and flow rate data (X i ,D iU i ); S103, determining an entropy parameter of the river cross section according to historical flow rate data; and obtaining two-dimensional flow rate field data of the river cross section based on an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and time sequence data corresponding to the surface flow rate; S103 specifically includes: S31, determining an entropy function by using a formula ; wherein, is the entropy function, M is the entropy parameter, U m is the average flow rate of the river cross section, U max is the maximum flow rate of the river cross section, and e is a natural logarithm; At least 20 groups of historical flow rate data are used to determine the entropy function by using the formula in S31, and the entropy parameter M is determined by inverse calculation.

[0023] S32, determining the entropy parameter of the river cross section according to the entropy function; S33, determining the cross-sectional point flow rate of the river cross section based on the entropy model according to the entropy parameter of the river cross section and time sequence data corresponding to the surface flow rate; S32 specifically includes: S3.1, determining the maximum flow rate downward shift value on the vertical line; According to the distance of the vertical line x i from the left bank, the water depth D(x i ) at which the vertical line is located is considered to determine the maximum flow rate downward shift value δ(x i ) on the vertical line, and the formula is as follows: ; wherein, x i is the i-th measuring line position from the left bank, δ(x i ) is the maximum flow rate downward shift value at the measuring line position, is the water depth at the measuring line position, is determined through an iterative cycle.

[0024] S3.2, determining the maximum flow rate on the vertical line; According to the flow rate formula of the entropy model, the surface flow rate of the river cross section, the entropy parameter and the maximum flow rate downward shift value on the vertical line are used as basic data to obtain the maximum flow rate value V(x i ) on the vertical line at each vertical line x : ; wherein, is the surface flow rate value at the x i measuring line.

[0025] S3.3, determining the cross-sectional point flow rate; Using the entropy parameter of the river cross section and the maximum vertical velocity value, the point velocity of the cross section is calculated according to the velocity formula of the entropy model to determine the topographic measurement point (x i ,y i ) at the river cross section : ; in, is x i The distance from the water surface where the maximum velocity on the vertical line is, is the point in the vertical direction, is the total number of perpendicular lines in the river cross section.

[0026] S3.4, iterative loop calculation of flow velocity at cross-section points; Update the downward shift value of the maximum flow velocity and expand the iterative loop calculation process of the flow field. Each loop process obtains an entropy parameter. When the difference between the calculated value of the entropy function and the measured value of the entropy function meets the following formula, the calculation process is completed and the cross-sectional flow velocity of this loop is output; ; Among them, Φ(M p ) is the calculated value of the entropy function, Φ(M obs ) is the measured value of the entropy function, ε represents the difference, which is generally taken as 0.01.

[0027] S34, based on the flow velocity at the cross-section points, the two-dimensional flow velocity field data of the river cross section is obtained using the two-dimensional inverse distance weighted interpolation (IDW) algorithm.

[0028] Among them, the two-dimensional inverse distance weighted interpolation algorithm is based on the inverse of the distance as the weight for weighted averaging to obtain the value of the unknown point. The specific process of the two-dimensional inverse distance weighted interpolation algorithm is as follows: S1, determine the coordinates and range of the position to be interpolated, as well as the coordinates and values ​​of the known sample points.

[0029] S2, for each position to be interpolated, calculates its distance to all known sample points and converts the distance into a weight. The weight is proportional to the inverse of the distance and is calculated using the following formula: ; Among them, ω i is the weight of the i-th sample point, d i is the distance between the sample point and the position to be interpolated. p is an adjustable parameter, which is generally set to 2 (Euclidean distance) or 3 (Manhattan distance).

[0030] S3, use the weight of each sample point to perform weighted averaging of its function value to obtain the interpolation result: ; Among them, Z(x0) is the estimated value of the position x0 to be interpolated, Z(x i ) is the known position x i Observed value of .

[0031] S104, determining a hydraulic complexity index based on the two-dimensional velocity field data of the river cross section; and using a two-dimensional inverse distance weighted interpolation algorithm based on the hydraulic complexity index to obtain a two-dimensional distribution field of the hydraulic complexity index of the river cross section, and Figs. 4-6 As shown; the hydraulic complexity index includes: kinetic energy gradient parameter and kinetic energy change parameter; S104 specifically includes: Using the formula Determine the kinetic energy gradient parameter M1; Using the formula Determine the kinetic energy change parameter M2; Where V1 and V2 are the flow velocities at adjacent measuring points, Δs is the distance between measuring points, and V ave is the average velocity of adjacent points, V min is the smaller velocity value among adjacent points, V is the velocity, and s is the distance.

[0032] S105, quantitatively assess the river ecosystem based on the two-dimensional distribution field of hydraulic complexity indicators in the river cross section.

[0033] S105 specifically includes: S51, determining an average value of each measuring line based on the two-dimensional distribution field of the hydraulic complexity index of the river cross section; and determining a planar distribution value of the hydraulic complexity index based on the average value of the measuring line; S52, determining a threshold value of the hydraulic complexity index according to the plane distribution value of the hydraulic complexity index; Using the formula m=∑M1(x i ,y i ) / k determines the measuring point (x i ,y i ) The threshold m of the kinetic energy gradient parameter M1; Using the formula n=∑M2(x i ,y i ) / k determine the measuring point (x i ,y i ) The threshold value n of the kinetic energy change parameter M2; Where k is the total number of measurement points within the river plane.

[0034] Comparing the hydraulic complexity index of the river cross section with the threshold of the hydraulic complexity index to obtain a comparison result; Specifically, according to the relative size of the river cross-section hydraulic complexity index (M1, M2), the river is divided into several sub-regions on the plane, wherein the regions corresponding to M1 and M2 exceeding m and n are analyzed as high resistance regions and energy intensive regions respectively, and the river aquatic habitat is comprehensively evaluated.

[0035] S53, determining the evaluation result of the river ecosystem according to the comparison result.

[0036] The region with higher M1 is a high resistance region, which may exert greater movement restriction on weak swimming species, while the rheophilic organisms adapted to high resistance environment are more suitable, and has an important influence on the energy consumption and survival strategy of organisms.

[0037] The region with higher M2 is usually an energy trap, which is usually an energy intensive habitat, and has important significance for filter-feeding macroinvertebrates and juvenile fish and other organisms. These regions provide rich food resources and suitable habitats, which have an important influence on the survival and reproduction of organisms.

[0038] The present application can solve the problems of high cost, difficulty in real-time online monitoring and evaluation of river ecological evaluation. Using non-contact measurement equipment and built-in entropy model calculation, real-time online evaluation of river biological habitat state can be realized; the present application is based on non-contact measurement of river channel, has the advantages of high automation, strong adaptability to complex environment, high measurement safety, etc., can guarantee the measurement work in high flood flow period, and avoid the safety hidden danger of hydrological measurement personnel. At the same time, it also has the characteristics of low disturbance to water environment, etc., which can obtain the flow data of the surface of the river without contacting the water, thereby providing more abundant and efficient information for the evaluation of river ecosystem. The present application makes full use of surface flow velocity data, not only realizes the reconstruction of flow velocity field of river channel cross-section, but also combines ecological hydraulics and applies to water ecological evaluation direction to analyze and evaluate the behavior of aquatic organisms and their habitat. The present application has strong applicability and can be applied to different river types of large rivers, medium and small rivers, especially effectively cope with the challenges of difficult actual measurement and ecological evaluation in large rivers. The present application has strong operability, can be programmed, has portability and nesting, can realize automatic and rapid river ecosystem evaluation through entropy model algorithm, provides technical support for river intelligent perception and comprehensive detection of river basin, provides support for biodiversity protection, and helps to improve the river fine control ability and the level of river basin management ability.

[0039] Based on the same inventive concept, the embodiments of the present application also provide a water-power-based river ecology evaluation device for implementing the water-power-based river ecology evaluation method as described above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme described in the above method, and therefore the specific limitations in one or more water-power-based river ecology evaluation device embodiments provided below can refer to the limitations of the water-power-based river ecology evaluation method described above, which will not be repeated here.

[0040] In one exemplary embodiment, a water-power-based river ecology evaluation device is provided, comprising: a data acquisition module configured to acquire basic characteristic elements of a river cross section and corresponding historical flow rate data; the basic characteristic elements include measurement data of riverbank surface topography and underwater topography, measurement point coordinates and corresponding relative elevations acquired at a certain distance on the river cross section; and the historical flow rate data includes average flow rate and maximum flow rate of the river cross section; a real-time data acquisition module configured to acquire surface flow rate of the river cross section and corresponding time series data by using a non-contact measurement device; the non-contact measurement device includes a radar speed measuring instrument or a large particle image velocimetry; a two-dimensional flow rate field data determination module configured to determine an entropy parameter of the river cross section according to the historical flow rate data; and obtain two-dimensional flow rate field data of the river cross section based on an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and the time series data corresponding to the surface flow rate; a hydraulic complexity index determination module configured to determine a hydraulic complexity index according to the two-dimensional flow rate field data of the river cross section; and obtain a two-dimensional distribution field of the hydraulic complexity index of the river cross section by using a two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; the hydraulic complexity index includes kinetic energy gradient parameters and kinetic energy variation parameters; an evaluation module configured to quantitatively evaluate a river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0041] In an example embodiment, a computer device is provided, which can be a server or a terminal. The computer device comprises a processor, a memory, an input / output interface and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for running the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is configured to exchange information between the processor and external devices. The communication interface of the computer device is configured to communicate with external terminals through network connection. The computer program is executed by the processor to implement a water-power-based river ecological assessment method.

[0042] In an example embodiment, a computer device is provided, which comprises a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps in the above method embodiments.

[0043] In an example embodiment, a computer readable storage medium is provided, which stores a computer program. The computer program is executed by a processor to implement the steps in the above method embodiments.

[0044] In an example embodiment, a computer program product is provided, which comprises a computer program. The computer program is executed by a processor to implement the steps in the above method embodiments.

[0045] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of the related data need to comply with relevant regulations.

[0046] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical storage, high-density embedded non-volatile memory, resistive memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0047] The database involved in the embodiments provided in the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a blockchain, etc., without being limited thereto. The processor involved in the embodiments provided in the present application can be a general processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0048] In the present application, all actions of obtaining signals, information or data are performed under the premise of complying with the corresponding data protection regulations and policies of the country where the device is located, and under the premise of obtaining authorization from the owner of the corresponding device.

[0049] The technical features of the above embodiments can be combined arbitrarily. In order to make the description simple, all possible combinations of the technical features in the above embodiments are not described, but as long as the combinations of the technical features do not exist contradictory, they should be considered as the scope of the present application.

[0050] The principles and implementation manners of the present application are described herein by using specific examples, and the above examples are only used to help understand the method of the present application and its core idea; meanwhile, for those skilled in the art, according to the idea of the present application, the specific implementation manners and application ranges will have changes. In conclusion, the content of the specification should not be understood as a limitation of the present application.

Claims

1. A river ecological assessment method based on hydrodynamics, characterized in that: The river ecological assessment method based on hydrodynamics includes: Obtaining basic characteristic elements of a river cross section and corresponding historical flow velocity data; the basic characteristic elements include: measurement data of riverbank surface topography and underwater topography obtained at certain intervals in the river cross section, coordinates of measurement points, and corresponding relative elevations; the historical flow velocity data includes: average flow velocity and maximum flow velocity in the river cross section; Using non-contact measurement equipment to obtain the surface flow velocity of the river cross section and the corresponding time series data; the non-contact measurement equipment includes: a radar velocity measuring instrument or a large particle image velocimeter; Determine the entropy parameter of the river cross section based on historical flow velocity data; and obtain the two-dimensional flow velocity field data of the river cross section based on the entropy parameter and the time series data corresponding to the surface flow velocity, based on the entropy model and the two-dimensional inverse distance weighted interpolation algorithm; Determining a hydraulic complexity index based on two-dimensional velocity field data of a river cross section; and obtaining a two-dimensional distribution field of the hydraulic complexity index of the river cross section using a two-dimensional inverse distance weighted interpolation algorithm based on the hydraulic complexity index; the hydraulic complexity index includes a kinetic energy gradient parameter and a normalized energy change rate; The river ecosystem is quantitatively assessed based on the two-dimensional distribution field of hydraulic complexity indicators in the river cross section.

2. The river ecological assessment method based on hydrodynamics according to claim 1, characterized in that: The acquisition of the basic characteristic elements of the river cross section and the corresponding historical flow velocity data specifically includes: Using an ultrasonic sounder or a traveling Doppler profiler, underwater topography is measured at set distances on the cross section of the river to obtain underwater topography measurement data; Using a level and a distance meter, measure the relative elevation and distance of the riverbanks on both sides along the cross section of the river, starting from a relative elevation zero coordinate point, to fixed columns on both sides, to obtain measurement data of the riverbank surface topography; the relative elevation zero coordinate point is a fixed column on the riverbank on one side of the river cross section; The relative elevation is converted using the relative elevation zero coordinate point, and the three-dimensional cross-sectional topography of the river channel and the topographic measurement points are determined based on the measurement data of the underwater topography and the river bank surface topography; According to the three-dimensional cross-section topography of the river channel and the topographic measurement points, the corresponding historical flow velocity data is obtained.

3. The river ecological assessment method based on hydrodynamics according to claim 1, characterized in that: The entropy parameter of the river cross section is determined based on the historical flow velocity data; and based on the time series data corresponding to the entropy parameter and the surface flow velocity, the two-dimensional flow velocity field data of the river cross section is obtained based on the entropy model and the two-dimensional inverse distance weighted IDW interpolation algorithm, specifically including: Using the formula Determine the entropy function; where, is the entropy function, M is the entropy parameter, U m is the average flow velocity of the river cross section, U max is the maximum flow velocity in the river cross section, and e is the natural logarithm; Determine the entropy parameter of the river cross section based on the entropy function; According to the entropy parameters of the river cross section and the time series data corresponding to the surface flow velocity, the flow velocity at the cross section of the river is determined based on the entropy model; According to the flow velocity at the cross-section points, the two-dimensional flow velocity field data of the river cross section is obtained using the two-dimensional inverse distance weighted interpolation algorithm.

4. The river ecological assessment method based on hydrodynamics according to claim 1, characterized in that: Based on the two-dimensional velocity field data of the river cross section, the hydraulic complexity index is determined, including: Using the formula Determine the kinetic energy gradient parameter M1; Using the formula Determine the normalized energy change rate M2; Where V1 and V2 are the flow velocities at adjacent measuring points, Δs is the distance between measuring points, and V ave is the average velocity of adjacent points, V min is the smaller flow velocity value among adjacent points.

5. The river ecological assessment method based on hydrodynamics according to claim 1, characterized in that: The quantitative assessment of the river ecosystem based on the two-dimensional distribution field of the hydraulic complexity index of the river cross section specifically includes: The average value of each survey line is determined based on the two-dimensional distribution field of the hydraulic complexity index of the river cross section; and the plane distribution value of the hydraulic complexity index is obtained using the two-dimensional inverse distance weighted interpolation algorithm based on the average value of the survey line; Determine the average value of the hydraulic complexity index according to the plane distribution value of the hydraulic complexity index; The hydraulic complexity index of the river cross section is compared with the average value of the hydraulic complexity index to obtain the comparison result; Based on the comparison results, the assessment results of the river ecosystem are determined.

6. The method for river ecological assessment based on hydrodynamics according to claim 5, characterized in that: Determining the average value of the hydraulic complexity index according to the plane distribution value of the hydraulic complexity index specifically includes: Using the formula m=∑M1(x i ,y i ) / k determine the measuring point (x i ,y i ) the average value m of the kinetic energy gradient parameter M1; Using the formula n=∑M2(x i ,y i ) / k determine the measuring point (x i ,y i ) The average value n of the normalized energy change rate M2; Where k is the total number of measurement points within the river plane.

7. A river ecological assessment device based on hydrodynamics, characterized in that: The river ecology assessment equipment based on hydrodynamics includes: A data acquisition module is used to acquire basic characteristic elements of a river cross section and corresponding historical flow velocity data; the basic characteristic elements include: measurement data of river bank surface topography and underwater topography acquired at a certain distance in the river cross section, coordinates of measurement points, and corresponding relative elevations; the historical flow velocity data includes: average flow velocity and maximum flow velocity in the river cross section; A real-time data acquisition module is used to obtain the surface flow velocity of the river cross section and the corresponding time series data using a non-contact measurement device; the non-contact measurement device includes: a radar velocity measuring instrument or a large particle image velocimeter; A two-dimensional velocity field data determination module is used to determine the entropy parameter of the river cross section based on the historical velocity data; and based on the time series data corresponding to the entropy parameter and the surface velocity, based on the entropy model and the two-dimensional inverse distance weighted interpolation algorithm, the two-dimensional velocity field data of the river cross section is obtained; A hydraulic complexity index determination module is configured to determine a hydraulic complexity index based on two-dimensional velocity field data of a river cross section; and to obtain a two-dimensional distribution field of the hydraulic complexity index of the river cross section using a two-dimensional inverse distance weighted interpolation algorithm based on the hydraulic complexity index; the hydraulic complexity index includes a kinetic energy gradient parameter and a kinetic energy change parameter; The evaluation module is used to quantitatively evaluate the river ecosystem based on the two-dimensional distribution field of hydraulic complexity indicators in the river cross section.

8. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the hydrodynamic-based river ecological assessment method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the river ecological assessment method based on hydrodynamics according to any one of claims 1 to 6 is implemented.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the river ecological assessment method based on hydrodynamics according to any one of claims 1 to 6 is implemented.

Images