A water-power-based river ecological evaluation method, device, medium and product

Abstract

The application discloses a water-power-based river ecological evaluation method and device, medium and product, relates to the field of ecological hydraulics, and comprises the following steps: acquiring basic characteristic elements of a river cross section and historical flow rate data; acquiring time series data of river surface flow rate by using a non-contact measurement device; determining an entropy parameter of the river cross section according to the historical flow rate data; obtaining two-dimensional flow rate field data of the river cross section by using an entropy model and a two-dimensional inverse distance weighted interpolation algorithm according to the entropy parameter and the time series data; calculating two hydraulic complexity indexes, namely, a kinetic energy gradient parameter and a normalized energy change rate, according to the two-dimensional flow rate field data; obtaining a two-dimensional distribution field of the hydraulic complexity indexes by using the two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity indexes; and quantitatively evaluating a river ecological system according to the two-dimensional distribution field. The application can fully utilize non-contact measurement and realize rapid and accurate evaluation of the river ecological system.

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 velocity gradient and secondary flow in rivers provide 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 flow 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:

[0007] In a first aspect, the present application provides a river ecological evaluation method based on water power, which comprises:

[0008] 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, measurement point coordinates and corresponding relative elevations obtained by a certain distance in the river cross section; the historical flow velocity data includes: average flow velocity and maximum flow velocity of the river cross section;

[0009] The surface flow velocity of the river cross section and corresponding time sequence data are acquired by using a non-contact measurement device; the non-contact measurement device comprises a radar velocity measuring instrument or a large particle image velocimetry;

[0010] The entropy parameter of the river cross section is determined according to the historical flow velocity data; and the two-dimensional flow velocity field data of the river cross section are obtained 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;

[0011] The hydraulic complexity index is determined according to the two-dimensional flow velocity field data of the river cross section; and the two-dimensional distribution field of the hydraulic complexity index of the river cross section is obtained by using a two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; the hydraulic complexity index comprises a kinetic energy gradient parameter and a kinetic energy change parameter;

[0012] The river ecosystem is quantitatively evaluated according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0013] Optionally, the basic characteristic elements of the river cross section and corresponding historical flow velocity data are acquired, specifically comprising:

[0014] The underwater topography is measured at a set distance on the river cross section by using an ultrasonic wave detector or a walking Doppler profile flow velocity meter to acquire the measurement data of the underwater topography;

[0015] The relative elevations and distances of the two banks are measured along the direction of the river cross section by using a level and a range finder to the fixed columns on the two banks with a relative elevation zero coordinate point as the fixed column of one bank of the river cross section, to obtain the measurement data of the bank surface topography;

[0016] The relative elevation is converted by using the relative elevation zero coordinate point, and the three-dimensional cross-sectional topography of the river cross section and the topographic measurement points are determined according to the measurement data of the underwater topography and the measurement data of the bank surface topography;

[0017] The corresponding historical flow velocity data are acquired according to the three-dimensional cross-sectional topography of the river cross section and the topographic measurement points.

[0018] Optionally, the entropy parameter of the river cross section is determined according to the historical flow velocity data; and the two-dimensional flow velocity field data of the river cross section are obtained based on an entropy model and a two-dimensional inverse distance weighted IDW interpolation algorithm according to the entropy parameter and the time sequence data corresponding to the surface flow velocity, specifically comprising:

[0019] The entropy function is determined by using a formula ; wherein, is the entropy function, M is the entropy parameter, U m is the average flow velocity of the river cross section, U maxis the maximum flow velocity of the river cross section, e is a natural logarithm;

[0020] determining an entropy parameter of the river cross section according to an entropy function;

[0021] determining a section point flow velocity of the river cross section based on an entropy model according to the entropy parameter of the river cross section and time sequence data corresponding to the surface flow velocity;

[0022] obtaining two-dimensional flow velocity field data of the river cross section by using a two-dimensional inverse distance weighted interpolation algorithm according to the section point flow velocity.

[0023] Optionally, a hydraulic complexity index is determined according to the two-dimensional flow velocity field data of the river cross section, and specifically includes:

[0024] a kinetic energy gradient parameter M1 is determined by using a formula

[0025] a normalized energy change rate M2 is determined by using a formula

[0026] wherein V1 and V2 are flow velocities of adjacent measuring points, Δs is a measuring point interval, V ave is an average flow velocity of the adjacent points, V min is a smaller flow velocity value in the adjacent points.

[0027] Optionally, the river ecosystem is quantitatively evaluated according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section, and specifically includes:

[0028] an average value of each measuring line is determined according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section, and a plane distribution value of the hydraulic complexity index is obtained by using a two-dimensional inverse distance weighted interpolation algorithm according to the average value of the measuring line;

[0029] an average value of the hydraulic complexity index is determined according to the plane distribution value of the hydraulic complexity index;

[0030] a comparison result is obtained by comparing the hydraulic complexity index of the river cross section with the average value of the hydraulic complexity index;

[0031] an evaluation result of the river ecosystem is determined according to the comparison result.

[0032] Optionally, the average value of the hydraulic complexity index is determined according to the plane distribution value of the hydraulic complexity index, and specifically includes:

[0033] an average value m of kinetic energy gradient parameters M1 of measuring points (x i ,y i ) is determined by using a formula m = ∑M1(x i ,y i ) / k; ​​

[0034] The average value n of the normalized energy variation rate M2 of the measuring point (x i ,y i ) is determined by the formula n=∑M2(x i ,y i ) / k.

[0035] Wherein, k is the total number of measuring points in the river channel plane range.

[0036] In a second aspect, the present application provides a water-power-based river ecological evaluation device, which comprises:

[0037] A data acquisition module is configured to acquire 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, measuring 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.

[0038] A 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.

[0039] A 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.

[0040] A 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 variation parameter.

[0041] An evaluation module is configured to quantitatively evaluate a river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0042] In a third aspect, the present application provides a computer device, which comprises: a memory, a processor and a computer program stored in the memory and executable on the processor; the processor executes the computer program to realize the water-power-based river ecological evaluation method.

[0043] In a fourth aspect, the present application provides a computer readable storage medium, which stores a computer program; when the computer program is executed by a processor, the water-power-based river ecological evaluation method is realized.

[0044] In a fifth aspect, the present application provides a computer program product comprising a computer program which, when executed by a processor, implements the water-power-based river ecological assessment method.

[0045] According to the specific embodiments provided in the present application, the present application has the following technical effects:

[0046] The present application provides a water-power-based river ecological assessment method, device, medium and product. The surface flow velocity of a river cross section and the corresponding time sequence data are obtained by using a non-contact measurement device, which fully embodies the efficiency, safety and ecological friendliness of non-contact measurement, avoiding the interference of traditional contact measurement on the river ecological system. The method has high automation, strong adaptability to complex environments, high measurement safety, etc., can guarantee the measurement work during the high flood flow period, and avoid the safety hazards to the hydrological measurement personnel. The heterogeneity of the river ecological system (river biological habitat) is quantitatively evaluated by combining the entropy model and the hydraulic complexity index, and then the automatic and rapid river ecological system evaluation is realized. The present application provides a scientific basis for river ecological protection and biodiversity research, and can solve the problems of high cost and difficulty in real-time online monitoring and evaluation of river ecological assessment, which helps to improve the river fine control ability and the level of basin management and management ability. BRIEF DESCRIPTION OF DRAWINGS

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

[0048] Fig. 1 The figure is a flowchart of the water-power-based river ecological assessment method in an embodiment of the present application.

[0049] Fig. 2 The figure is a schematic diagram of the principle of the water-power-based river ecological assessment method in an embodiment of the present application.

[0050] Fig. 3 The figure is a technical route flowchart of the water-power-based river ecological assessment method in an embodiment of the present application.

[0051] Fig. 4 The figure is a distribution diagram of M1 and M2 sections when the river flow is 124578 m / s. 3

[0052] Fig. 5 The figure is a distribution diagram of M1 and M2 sections when the river flow is 151663 m / s. 3 ​M1 and M2 cross-section distribution diagrams at t=0 and t=1;

[0053] Fig. 6 The river flow is 154474 m 3 M1 and M2 cross-section distribution diagrams at t=0 and t=1. DETAILED DESCRIPTION

[0054] The technical solutions in the embodiments of the present application will be described clearly and completely below with the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the present application.

[0055] The above purposes, features and advantages of the present application can be more obvious and easy to understand. The present application will be further described in detail below with the drawings and specific embodiments.

[0056] In an exemplary embodiment, as Figs. 1-3 shown, a water-based river ecological assessment method is provided, which comprises S101 to S105. Wherein:

[0057] S101, acquiring 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;

[0058] S101 specifically includes:

[0059] S11, using an ultrasonic detector or a walking Doppler profile flowmeter to measure underwater topography on the river cross-section at a set distance to obtain measurement data of the underwater topography;

[0060] 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 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;

[0061] S13, converting the relative elevations 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 based on the measurement data of the underwater topography and the measurement data of the riverbank surface topography;

[0062] S14, acquiring corresponding historical flow rate data based on the three-dimensional cross-section topography of the river cross-section and the topographic measurement point.

[0063] S102, acquiring surface flow velocity and corresponding time series data of the river cross section by using a non-contact measurement device; the non-contact measurement device includes but is not limited to a radar velocity measuring instrument or a large particle image velocimeter;

[0064] The selection of the non-contact measurement device is used to measure the surface flow velocity and water level of the river. The radar velocity measuring instrument is driven by an automatic control device to measure the surface flow velocity and water level at different positions of the river cross section at a certain interval, and record the position and flow velocity data (X i ,D i ,U i ) at time t;

[0065] S103, determining the entropy parameter of the river cross section according to the historical flow velocity data; and obtaining the two-dimensional flow 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 series data corresponding to the surface flow velocity;

[0066] S103 specifically includes:

[0067] S31, determining the entropy function by using the formula ; wherein, 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 of the river cross section, and e is the natural logarithm;

[0068] At least 20 groups of historical flow velocity data are used to determine the entropy function by using the formula in S31, and then the entropy parameter M is determined by inverse calculation.

[0069] S32, determining the entropy parameter of the river cross section according to the entropy function;

[0070] S33, determining the cross-sectional point flow velocity of the river cross section based on the entropy model according to the entropy parameter of the river cross section and the time series data corresponding to the surface flow velocity;

[0071] S32 specifically includes:

[0072] S3.1, determining the maximum flow velocity downward displacement value on the vertical line;

[0073] According to the distance x i from the left bank, considering the water depth D(x i ) at which the measuring line is located, the maximum flow velocity downward displacement value δ(x i ) on the vertical line is determined, and the formula is as follows:

[0074] ;

[0075] wherein, x iis the distance of the i-th measuring line position from the left bank, δ(x i ) is the maximum flow velocity downward displacement value at the measuring line position, is the water depth at the measuring line position, is determined through an iterative cycle.

[0076] S3.2, determination of the maximum vertical flow velocity;

[0077] According to the flow velocity formula of the entropy model, the maximum vertical flow velocity value at each vertical line x i of the cross section is obtained based on the surface flow velocity, the entropy parameter and the maximum flow velocity downward displacement value of the cross section of the river channel:

[0078] ;

[0079] wherein, V(x i ) is the surface flow velocity value at the measuring line.

[0080] S3.3, determination of the point flow velocity of the cross section;

[0081] The point flow velocity of the cross section is calculated according to the flow velocity formula of the entropy model using the entropy parameter and the maximum vertical flow velocity value of the cross section of the river channel, so as to determine the flow velocity data of the cross section of the river channel at the topographic survey point position (x i , y i ):

[0082] ;

[0083] wherein, h(x i ) is the distance of the maximum flow velocity from the water surface at the vertical line, is the point position in the vertical direction, is the total number of vertical lines of the cross section of the river channel.

[0084] S3.4, iterative cycle calculation of the point flow velocity of the cross section;

[0085] The maximum flow velocity downward displacement value is updated, and the iterative cycle calculation process of the flow field is carried out. An entropy parameter is obtained in each cycle process, and the calculation value of the entropy function and the measured value of the entropy function are compared. When the difference value meets the following formula, the calculation process is completed, and the point flow velocity of the cross section of this cycle is outputted;

[0086] ;

[0087] wherein, Φ(M p ) is the calculation value of the entropy function, Φ(M obs ) is the measured value of the entropy function, and ε represents the difference value, which is generally 0.01.

[0088] ​​​​S34, according to the cross-section point flow rate, using two-dimensional inverse distance weighted interpolation (IDW) algorithm, the river cross-section two-dimensional flow field data is obtained.

[0089] The two-dimensional inverse distance weighted interpolation algorithm is to take the reciprocal of the distance as the weight for weighted average, so as to obtain the value of the unknown point, and the specific process of the two-dimensional inverse distance weighted interpolation algorithm is as follows:

[0090] S1, the coordinates and range of the position to be interpolated are determined, and the coordinates and values of the known sample points are determined.

[0091] S2, for each position to be interpolated, the distance between it and all known sample points is calculated, and the distance is converted into weight. The weight is proportional to the reciprocal of the distance, and the following formula is used to calculate:

[0092] ;

[0093] Wherein, ω 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, and p is an adjustable parameter, generally taking the value of 2 (Euclidean distance) or 3 (Manhattan distance).

[0094] S3, the function value of each sample point is weighted and averaged using the weight of each sample point to obtain the interpolation result:

[0095] ;

[0096] Wherein, Z(x0) is the estimated value of the position x0 to be interpolated, Z(x i ) is the observed value of the known position x i .

[0097] S104, according to the two-dimensional flow field data of the river cross-section, the hydraulic complexity index is determined; and according to the hydraulic complexity index, the two-dimensional distribution field of the hydraulic complexity index of the river cross-section is obtained by using the two-dimensional inverse distance weighted interpolation algorithm, and as shown in Figs. 4-6 ; The hydraulic complexity index includes kinetic energy gradient parameter and kinetic energy change parameter;

[0098] S104 specifically includes:

[0099] The kinetic energy gradient parameter M1 is determined by the formula ;

[0100] The kinetic energy change parameter M2 is determined by the formula ;

[0101] Wherein, V1 and V2 are the flow rates of adjacent measuring points, Δs is the distance between measuring points, V ave is the average flow rate of adjacent points, and V minFor the smaller flow rate value in adjacent points, V is the flow rate, and s is the distance.

[0102] S105, quantitatively evaluate the river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

[0103] S105 specifically includes:

[0104] S51, determine the average value of each measuring line according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section; and determine the planar distribution value of the hydraulic complexity index according to the average value of the measuring line;

[0105] S52, determine the threshold value of the hydraulic complexity index according to the planar distribution value of the hydraulic complexity index;

[0106] The threshold value m of the kinetic energy gradient parameter M1 of the measuring point (x i ,y i ) is determined by the formula m = ∑M1(x i ,y i ) / k;

[0107] The threshold value n of the kinetic energy change parameter M2 of the measuring point (x i ,y i ) is determined by the formula n = ∑M2(x ,y

[0001] ) / k;

[0108] Wherein, k is the total number of calculation measuring points in the river plane range.

[0109] The hydraulic complexity index of the river cross section is compared with the threshold value of the hydraulic complexity index to obtain a comparison result;

[0110] Specifically, according to the relative size of the hydraulic complexity index (M1, M2) of the river cross section, 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.

[0111] S53, determine the evaluation result of the river ecosystem according to the comparison result.

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

[0113] The area with higher M2 is usually an energy trap, and these areas are usually energy-intensive habitats, which are of great significance to filter-feeding macroinvertebrates and juvenile fish and other organisms. These areas provide abundant food resources and suitable habitats, which have a significant impact on the survival and reproduction of organisms.

[0114] The application can solve the problems of high cost, difficulty in real-time online monitoring and evaluation of river ecological evaluation. The use of non-contact measurement equipment and built-in entropy model calculation can realize real-time online evaluation of river biological habitat status; the application is based on non-contact measurement of the river channel, has the advantages of high automation, strong adaptability to complex environment, high measurement safety, can ensure the measurement work during high flood flow period, and avoids the safety hazards to the hydrological measurement personnel. At the same time, it also has the characteristics of low disturbance to the water environment, etc., which can obtain the flow data of the river surface without contacting the water, thereby providing more abundant and efficient information for the evaluation of the river ecosystem. The application makes full use of the surface flow rate data, not only realizes the reconstruction of the flow field of the whole cross section of the river channel, but also combines ecological hydraulics and is applied to the water ecological evaluation direction to analyze and evaluate the behavior of aquatic organisms and their habitats. The 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 application has strong operability, can be programmed, has portability and nesting, can realize automatic and rapid river ecosystem evaluation through the entropy model algorithm, provides technical support for river intelligent perception and comprehensive detection of the river basin, provides support for biodiversity protection, and helps to improve the fine control ability of the river and the management and management ability level of the river basin.

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

[0116] In one exemplary embodiment, a water power-based river ecological evaluation device is provided, comprising:

[0117] A data acquisition module is 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 obtained by a certain distance in the river cross section; and the historical flow rate data includes average flow rate and maximum flow rate of the river cross section;

[0118] a real-time data acquisition module configured to acquire surface flow velocities and corresponding time series data of a river cross section by using a non-contact measurement device, wherein the non-contact measurement device comprises a radar velocity measurement instrument or a large-scale particle image velocimetry;

[0119] a two-dimensional flow velocity field data determination module configured to determine an entropy parameter of the river cross section according to 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 series data corresponding to the surface flow velocities;

[0120] a hydraulic complexity index determination module 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, wherein the hydraulic complexity index comprises a kinetic energy gradient parameter and a kinetic energy variation parameter;

[0121] 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.

[0122] In an example embodiment, a computer device, which can be a server or a terminal, is provided. The computer device includes a processor, a memory, an input / output interface (I / O) 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 includes 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 a network connection. The computer program is executed by the processor to implement a water-power-based river ecosystem evaluation method.

[0123] In an example embodiment, a computer device is provided, which includes 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.

[0124] 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.

[0125] In an example embodiment, a computer program product is provided, including a computer program which, when executed by a processor, implements the steps of any of the above method embodiments.

[0126] It should be noted that the user information (including but not limited to user equipment 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 related data need to comply with relevant regulations.

[0127] It can be understood by those skilled in the art that all or part of the processes in the above-mentioned embodiments can be completed by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer readable storage medium, and when executed, can include the processes of the above-mentioned embodiments. Any reference to memory, database or other medium used in the embodiments provided by 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.

[0128] The database involved in the embodiments provided by 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 blockchain, etc., without being limited thereto. The processor involved in the embodiments provided by the present application can be a general-purpose 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.

[0129] 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 obtaining the authorization given by the owner of the corresponding device.

[0130] The technical features of the above embodiments can be combined in any manner. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described, however, as long as the combinations of the technical features do not exist contradictions, they should be considered as the scope of the present disclosure.

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

Claims

1. A water power-based river ecological assessment method, characterized by, The water-power-based river ecological evaluation method comprises the following steps: obtaining basic characteristic elements of a river cross section and corresponding historical flow velocity data; the basic characteristic elements comprise: measurement data of riverbank surface topography and underwater topography obtained by a certain distance on the river cross section, measurement point coordinates and corresponding relative elevations; the historical flow velocity data comprise: average flow velocity and maximum flow velocity of the river cross section; obtaining 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 comprises: a radar velocity measuring instrument or a large particle image velocimetry; determining an entropy parameter of the river cross section according to the historical flow velocity data; and obtaining 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; determining a hydraulic complexity index according to the two-dimensional flow velocity field data of the river cross section; and obtaining 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 comprises: a kinetic energy gradient parameter and a normalized energy change rate; quantitatively evaluating a river ecological system according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

2. The water power-based river ecological assessment method according to claim 1, characterized in that, The obtaining of the basic characteristic elements of the river cross section and the corresponding historical flow velocity data specifically comprises the following steps: measuring underwater topography on the river cross section at a set distance by using an ultrasonic wave detector or a walking type Doppler profile flow velocity meter to obtain measurement data of the underwater topography; measuring relative elevations and distances of both sides of the riverbank along the direction of the river cross section by using a level and a range finder to the fixed column on both banks to obtain 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; converting the relative elevation by using the relative elevation zero coordinate point and determining three-dimensional cross section 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 riverbank surface topography; obtaining corresponding historical flow velocity data according to the three-dimensional cross section topography of the river cross section and the topographic measurement points.

3. The water power-based river ecological assessment method according to claim 1, characterized in that, The determining of the entropy parameter of the river cross section according to the historical flow velocity data and the obtaining of the two-dimensional flow velocity 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 velocity specifically comprise the following steps: Using the formula determining the entropy function; wherein, 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 of the river cross section, e is the natural logarithm; determining the entropy parameter of the river cross section according to an entropy function; determining cross section point flow velocity 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 velocity; obtaining the two-dimensional flow velocity field data of the river cross section by using the two-dimensional inverse distance weighted interpolation algorithm according to the cross section point flow velocity.

4. The water power-based river ecological assessment method according to claim 1, wherein, The determining of the hydraulic complexity index according to the two-dimensional flow velocity field data of the river cross section specifically comprises the following steps: Using the formula determining the kinetic energy gradient parameter M1; Using the formula determining the normalized rate of energy change M2; where V1and V2are the flow rates of adjacent measuring points, Δs is the distance between measuring points, V ave is the average flow rate of adjacent points, V min is the smaller flow rate value of adjacent points.

5. The water power-based river ecological assessment method according to claim 1, wherein, The quantitatively evaluating of the river ecological system according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section specifically comprises the following steps: The average value of each measuring line is determined according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section; and the two-dimensional inverse distance weighted interpolation algorithm is used to obtain the planar distribution value of the hydraulic complexity index according to the average value of the measuring line; The average value of the hydraulic complexity index is determined according to the planar distribution value of the hydraulic complexity index; The comparison result is obtained by comparing the hydraulic complexity index of the river cross section with the average value of the hydraulic complexity index; The evaluation result of the river ecosystem is determined according to the comparison result.

6. The water power-based river ecological assessment method according to claim 5, wherein, The average value of the hydraulic complexity index is determined according to the planar distribution value of the hydraulic complexity index, and 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; The average value n of the normalized energy change rate M2 of the measuring point (x i ,y i ) is determined by the formula n=∑M2(x i ,y i ) / k; Wherein, k is the total number of calculation measuring points in the river plane range.

7. A water power based river ecological assessment apparatus, characterized by, The river ecological evaluation device based on water power includes: The data acquisition module is used to acquire the basic characteristic elements of the river cross section and the corresponding historical flow rate data; the basic characteristic elements include: the measurement data of the riverbank surface topography and the underwater topography obtained by a certain distance in the river cross section, the measuring point coordinates and the corresponding relative elevation; the historical flow rate data includes: the average flow rate and the maximum flow rate of the river cross section; The real-time data acquisition module is used to acquire the surface flow rate of the river cross section and the corresponding time sequence 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; The two-dimensional flow rate field data determination module is used to determine the entropy parameter of the river cross section according to the historical flow rate data; and the two-dimensional flow rate field data of the river cross section is obtained 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; The hydraulic complexity index determination module is used to determine the hydraulic complexity index according to the two-dimensional flow rate field data of the river cross section; and the two-dimensional distribution field of the hydraulic complexity index of the river cross section is obtained by using the two-dimensional inverse distance weighted interpolation algorithm according to the hydraulic complexity index; the hydraulic complexity index includes: kinetic energy gradient parameter and kinetic energy change parameter; The evaluation module is used to quantitatively evaluate the river ecosystem according to the two-dimensional distribution field of the hydraulic complexity index of the river cross section.

8. A computer device comprising: The memory, the processor and the computer program stored in the memory and executable on the processor are characterized in that the processor executes the computer program to realize the river ecological evaluation method based on water power in any one of claims 1-6.

9. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the river ecological evaluation method based on water power in any one of claims 1-6.

10. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to realize the river ecological evaluation method based on water power in any one of claims 1-6.

Images