A method, system, device and medium for evaluating health of rivers and lakes

Abstract

The application discloses a kind of river and lake health evaluation method, system, equipment and medium, wherein the method obtains the health index data and multi-year average daily flow sequence set of river and lake;The multi-year average daily flow sequence set is divided into annual cycle, and a plurality of annual cycle data are obtained, and the empirical distribution of adjacent annual cycle data in time is different;All the annual cycle data are processed to generate a grade evaluation table, and a cycle evaluation table corresponding to each annual cycle data is obtained;According to the health index data, all the cycle evaluation tables are evaluated by health mapping, and the health evaluation result of the river and lake is obtained.The method can effectively improve the objectivity and accuracy of river and lake health evaluation.The application relates to the field of hydrological science and technology.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of hydrology science, and in particular, to a river and lake health evaluation method, system, device and medium. BACKGROUND

[0002] River and lake health evaluation refers to a comprehensive and systematic evaluation of the health status of a river or lake to understand its ecological, environmental and functional status and provide a scientific basis for its protection and management. It is of great significance to improve the level of river and lake governance and protection and promote the construction of water ecological civilization.

[0003] At present, the existing river and lake health evaluation technology mainly realizes the health evaluation of rivers and lakes through biological index evaluation method, physical and chemical index evaluation method or comprehensive evaluation method. However, these methods are heavily dependent on expert experience and do not consider the distribution characteristics of the evaluation index, so the objectivity and accuracy of the health evaluation results of rivers and lakes are not satisfactory.

[0004] Therefore, the problems of the prior art still need to be solved and optimized. SUMMARY

[0005] The purpose of the present application is to at least partially solve one of the technical problems in the related art.

[0006] To this end, one purpose of the embodiments of the present application is to provide a river and lake health evaluation method, system, device and medium, wherein the method can effectively improve the objectivity and accuracy of river and lake health evaluation.

[0007] In order to achieve the above technical purpose, the technical solutions adopted by the embodiments of the present application include:

[0008] In a first aspect, the embodiments of the present application provide a river and lake health evaluation method, comprising:

[0009] obtaining health index data and a multi-year average daily flow sequence set of a river and lake;

[0010] dividing the multi-year average daily flow sequence set into intra-annual cycles to obtain a plurality of intra-annual cycle data, the empirical distribution of time-adjacent intra-annual cycle data being different;

[0011] generating a grade evaluation table for all the intra-annual cycle data to obtain a cycle evaluation table corresponding to each intra-annual cycle data;

[0012] performing health mapping evaluation on all the cycle evaluation tables according to the health index data to obtain a health evaluation result of the river and lake.

[0013] In addition, the method according to the above embodiments of the present application can also have the following additional technical features:

[0014] Further, in one embodiment of the present application, the multi-year average daily flow sequence set is divided into intra-annual cycles to obtain a plurality of intra-annual cycle data, including:

[0015] Obtaining a division threshold value;

[0016] Obtaining a division starting point, a division inspection point, and a division ending point;

[0017] According to the division threshold value, the division starting point, the division inspection point, and the division ending point, the multi-year average daily flow sequence set is divided into intra-annual cycles to obtain a plurality of intra-annual cycle data.

[0018] Further, in one embodiment of the present application, the multi-year average daily flow sequence set is divided into intra-annual cycles to obtain a plurality of intra-annual cycle data according to the division threshold value, the division starting point, the division inspection point, and the division ending point, including:

[0019] According to the division starting point and the division inspection point, the multi-year average daily flow sequence set is divided into a first sequence to obtain a first flow sequence;

[0020] According to the division inspection point and the division ending point, the multi-year average daily flow sequence set is divided into a second sequence to obtain a second flow sequence;

[0021] According to the first flow sequence, the second flow sequence is verified for empirical distribution consistency to obtain an empirical distribution verification result;

[0022] According to the division threshold value, the division ending point is verified for a first threshold value to obtain a first threshold value verification result;

[0023] If the empirical distribution verification result is that the empirical distributions are different and the first threshold value verification result is that the division ending point is less than the division threshold value, the first flow sequence is retained, the division starting point is updated according to the division inspection point, and the updated division starting point is used to update the division inspection point and the division ending point, and then the step of obtaining the division starting point, the division inspection point, and the division ending point is executed again; or, if the empirical distribution verification result is that the empirical distributions are different and the first threshold value verification result is that the division ending point is equal to the division threshold value, a plurality of intra-annual cycle data is obtained according to the current second flow sequence and all the first flow sequences.

[0024] Further, in one embodiment of the present application, the method further includes:

[0025] If the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is less than the division threshold, the division termination point is updated, and then the steps of obtaining the division start point, the division test point, and the division termination point are executed again.

[0026] Or,

[0027] If the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is equal to the division threshold, the second threshold verification is performed on the division test point according to the division threshold, and a second threshold verification result is obtained.

[0028] Further, in an embodiment of the present application, the method further comprises:

[0029] If the second threshold verification result is that the division test point is less than the division threshold, the division test point is updated, and the division termination point is updated according to the updated division test point, and then the steps of obtaining the division start point, the division test point, and the division termination point are executed again.

[0030] Or,

[0031] If the second threshold verification result is that the division test point is equal to the division threshold, a third flow sequence is obtained according to the current first flow sequence and the current second flow sequence, and a plurality of intra-annual cycle data are obtained according to the third flow sequence and all previous first flow sequences.

[0032] Further, in an embodiment of the present application, the expression form of the empirical distribution function of the first flow sequence is:

[0033]

[0034] Among them, is an empirical distribution function; is a first flow sequence; i is a division start point; j is a division test point; t is a time point of the first flow sequence on a time sequence; is an indicator function of the first flow sequence.

[0035] Further, in an embodiment of the present application, the intra-annual cycle data are subjected to a grade evaluation table generation process to obtain a cycle evaluation table, which comprises:

[0036] A plurality of intra-annual cycle index sub-data of the intra-annual cycle data are obtained, and each intra-annual cycle index sub-data corresponds to a different lake health evaluation index type;

[0037] performing data analysis on all the intra-annual period index sub-data, obtaining an index mean and an index standard deviation of each intra-annual period index sub-data;

[0038] obtaining the period evaluation table according to all the index means and all the index standard deviations.

[0039] Further, in an embodiment of the present application, the health mapping evaluation is performed on all the period evaluation tables according to the health index data, and a health evaluation result of the river and lake is obtained.

[0040] According to the health index data, a target evaluation table corresponding to the health index data is obtained by performing table screening on all the period evaluation tables.

[0041] According to the target evaluation table, the health index data is evaluated by index scoring, and the health evaluation result of the river and lake is obtained.

[0042] In a second aspect, an embodiment of the present application provides a river and lake health evaluation system, comprising:

[0043] A first processing unit is configured to obtain health index data of a river and lake and a multi-year average daily flow sequence set;

[0044] A second processing unit is configured to divide the multi-year average daily flow sequence set into intra-annual periods, and obtain a plurality of intra-annual period data, wherein the empirical distribution of time-adjacent intra-annual period data is different;

[0045] A third processing unit is configured to generate a period evaluation table corresponding to each intra-annual period data by performing grade evaluation table generation processing on all the intra-annual period data;

[0046] A fourth processing unit is configured to perform health mapping evaluation on all the period evaluation tables according to the health index data, and obtain a health evaluation result of the river and lake.

[0047] In a third aspect, an embodiment of the present application further provides an electronic device, comprising:

[0048] at least one processor;

[0049] at least one memory configured to store at least one program;

[0050] When the at least one program is executed by the at least one processor, the at least one processor implements the method of the first aspect.

[0051] In a fourth aspect, an embodiment of the present application further provides a computer readable storage medium, wherein a processor executable program is stored, and the processor executable program is used to implement the method of the first aspect when executed by the processor.

[0052] The advantages and beneficial effects of the present application will be partially given in the following description, partially will become obvious from the following description, or will be learned by the practice of the present application:

[0053] The method, system, device and medium for river and lake health evaluation are disclosed in the embodiments of the present application, wherein the method obtains health index data and a multi-year average daily flow sequence set of a river and lake; the multi-year average daily flow sequence set is divided into intra-annual cycles to obtain a plurality of intra-annual cycle data, and the empirical distribution of intra-annual cycle data adjacent in time is different; all the intra-annual cycle data are processed to generate a grade evaluation table to obtain a cycle evaluation table corresponding to each intra-annual cycle data; and the health index data are used to perform health mapping evaluation on all the cycle evaluation tables to obtain a health evaluation result of the river and lake. The method divides the multi-year average daily flow sequence set into intra-annual cycles to obtain a plurality of intra-annual cycle data, and the empirical distribution of each intra-annual cycle data is different from that of other intra-annual cycle data adjacent in time sequence, and then the cycle evaluation table generated from the intra-annual cycle data is evaluated based on the health index data, which can effectively improve the objectivity and accuracy of the river and lake health evaluation. BRIEF DESCRIPTION OF DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following introduces the drawings of the related technical solutions in the embodiments of the present application or the prior art. It should be understood that the drawings in the following introduction are only for the convenience of clearly expressing part of the embodiments of the technical solutions 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.

[0055] Figure 1 A flowchart of a river and lake health evaluation method provided by the embodiments of the present application is shown in the figure.

[0056] Figure 2 A structure diagram of a river and lake health evaluation system provided by the embodiments of the present application is shown in the figure.

[0057] Figure 3 A structure diagram of an electronic device provided by the embodiments of the present application is shown in the figure. DETAILED DESCRIPTION

[0058] The embodiments of the present application are described below in detail, examples of which are shown in the drawings, wherein the same or similar notations represent the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the drawings are exemplary only, for the purpose of explanation, and are not to be understood as a limitation of the present application. For the step numbers in the following embodiments, they are only set for the convenience of explanation, and no limitation is made on the order between the steps, and the execution order of each step in the embodiments can be adaptively adjusted according to the understanding of those skilled in the art.

[0059] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application belongs. The terminology used in the specification herein is for describing the embodiments of the present application only and is not intended to limit the present application.

[0060] At present, the existing river and lake health evaluation technology mainly realizes the health evaluation of rivers and lakes through biological index evaluation method, physical and chemical index evaluation method or comprehensive evaluation method. The biological index evaluation method takes aquatic organisms as indicator species, and reflects the health status of rivers and lakes by monitoring the changes of species, quantity and distribution of aquatic organisms. Specifically, the fish integrity index (F-IBI), the benthic animal integrity index (B-IBI), the Shannon-Wienner diversity index, the plankton evaluation method and the biological integrity index can be used to reflect the health status of rivers and lakes.

[0061] The physical and chemical index evaluation method reflects the health status of rivers and lakes by water quality, hydrological characteristics, river morphology or lake topography. The representative indexes of water quality evaluation include chemical oxygen demand (COD), biochemical oxygen demand (BOD), ammonia nitrogen, total phosphorus and dissolved oxygen. The hydrological characteristics evaluation takes flow, water level and flow rate as the starting point, and the representative index is IHA (Indicators of Hydrologic Alteration). The river morphology or lake topography evaluation includes indexes such as the degree of winding of the river, river width and water depth.

[0062] The comprehensive evaluation method analyzes the health status of rivers and lakes by selecting multiple indexes reflecting the health of rivers and lakes, following the path of weighting and scoring to build a comprehensive evaluation system. For example, by analyzing the mutual relationship among pressure, state and response through the pressure-state-response model (PSR), the health status of rivers and lakes is evaluated. Or, for example, by constructing a judgment matrix to determine the weight between the levels of the structure, the quantitative value of the comprehensive evaluation index is obtained to obtain the evaluation result of the health of rivers and lakes.

[0063] But since the biological index evaluation method and the physical and chemical index evaluation method basically adopt quantile method to divide the evaluation grade of the evaluation index, the commonly used quantiles include 10%, 25%, 50%, 75%, 80%, 90%, etc., the setting of the quantiles depends on the experience of experts, and cannot reflect the distribution law of the evaluation index itself, so that the biological index evaluation method and the physical and chemical index evaluation method cannot consider the distribution characteristics of the to-be-evaluated index, and the objectivity and accuracy of the health evaluation result of the river and lake are not satisfactory. In addition, the comprehensive evaluation method generally adopts expert experience method, subjective and objective combination method (objective weighting method + expert experience method) to determine the weight, and the evaluation index system and recursive hierarchical structure are strongly dependent on expert experience, so that the health evaluation result of the river and lake also lacks objectivity, and the accuracy of the health evaluation of the river and lake is not satisfactory.

[0064] In addition, the above method usually sets the lower standard of the river and lake in the dry season as the annual average standard, which does not consider the differences of different rivers and lakes in natural geographical conditions, climate characteristics, economic and social development and other factors, and the constraint degree of the health evaluation of the river and lake is small, which makes the reference value of the health evaluation result of the river and lake not high.

[0065] Therefore, the embodiments of the present application provide a river and lake health evaluation method, system, device and medium, wherein the method divides the multi-year average daily flow sequence set into several intra-annual period data, and the empirical distribution of each intra-annual period data and other intra-annual period data adjacent to it in time sequence is different, and then generates a period evaluation table based on the intra-annual period data, which can adaptively generate evaluation grades with higher adaptation degree for each intra-annual period of the river and lake while fully considering the distribution characteristics of the to-be-evaluated index, without relying on expert experience to determine the evaluation grade or weight, which is beneficial to improve the objectivity and accuracy of the health evaluation result of the river and lake.

[0066] In addition, the method constructs a period evaluation table through several intra-annual period index sub-data, which can fully consider the differences of the river and lake in natural geographical conditions, climate characteristics, economic and social development and other factors, which is beneficial to improve the constraint degree of the health evaluation of the river and lake, and further improve the reference value of the health evaluation result of the river and lake.

[0067] Reference Figure 1 In the embodiments of the present application, a river and lake health evaluation method comprises:

[0068] Step 110, obtaining health index data and a multi-year average daily flow sequence set of the river and lake;

[0069] In the embodiment of the present application, the health index data can be index data of the river and lake at a certain time date, and the index data can include at least one of a longitudinal connectivity index of the river and lake, an ecological flow satisfaction degree, a water plant community condition, and public satisfaction degree. The multi-year average daily flow sequence set includes a plurality of multi-year average daily flow data of the river and lake, each multi-year average daily flow data corresponds to a different time date within a year, and the number of multi-year average daily flow data in the multi-year average daily flow sequence set can be flexibly set according to actual conditions, and the embodiment of the present application will not be described here.

[0070] In step 120, the multi-year average daily flow sequence set is divided into intra-year cycles to obtain a plurality of intra-year cycle data, and the empirical distribution of time-adjacent intra-year cycle data is different.

[0071] In the embodiment of the present application, the multi-year average daily flow sequence set of the river and lake can be divided into intra-year cycles based on the time date sequence within a year, so as to obtain a plurality of intra-year cycle data, and the number of intra-year cycle data can be an integer greater than or equal to 1. In addition, in all obtained intra-year cycle data, two adjacent intra-year cycle data in the time sequence are subject to different empirical distributions, that is, there is no consistency between the two adjacent intra-year cycle data in the time sequence.

[0072] In some embodiments, the step 120 of dividing the multi-year average daily flow sequence set into intra-year cycles to obtain a plurality of intra-year cycle data includes:

[0073] A1, obtaining a division threshold;

[0074] A2, obtaining a division starting point, a division test point, and a division termination point;

[0075] In the embodiment of the present application, the multi-year average daily flow sequence set can be divided into intra-year cycles in a loop manner, so as to divide each multi-year average daily flow data in the multi-year average daily flow sequence set into a corresponding intra-year cycle, and thus obtain a plurality of intra-year cycle data.

[0076] It can be understood that the division threshold can be set according to actual conditions. Specifically, if the multi-year average daily flow sequence set is a sequence set of multi-year average daily flow data of a common year, the corresponding division threshold can be 365; or if the multi-year average daily flow sequence set contains a sequence set of multi-year average daily flow data of a leap year, that is, the multi-year average daily flow sequence set contains multi-year average daily flow data of February 29, the corresponding division threshold can be 366.

[0077] It should be noted that the division starting point is used to indicate the starting point of the intra-annual period division of the multi-year average daily flow data in the multi-year average daily flow sequence set in the current loop process; and the division test point and the division termination point are similar to the division starting point, which can be simply analogized.

[0078] It is worth mentioning that the division starting point, the division test point and the division termination point in the embodiments of the present application can be points on a time sequence, and the points on the time sequence can be associated with the multi-year average daily flow data of the corresponding time date in the multi-year average daily flow sequence set. For example, if January 1 is the first day of the year, the multi-year average daily flow data of January 1 can correspond to the first point on the time sequence; or if January 3 is the third day of the year, the multi-year average daily flow data of January 1 can correspond to the third point on the time sequence.

[0079] In addition, the division test point in the embodiments of the present application can be greater than or equal to the division starting point, and the division termination point can be greater than or equal to the division test point. Exemplarily, if the current loop process is the first loop process, the division starting point, the division test point and the division termination point can be pre-set values, such as the division starting point and the division test point being 1, and the division termination point being 2; or if the current loop process is the second or more loop process, the division starting point is the division starting point at the end of the last loop process, and the division test point and the division termination point are the same.

[0080] A3, according to the division threshold, the division starting point, the division test point and the division termination point, the multi-year average daily flow sequence set is divided into intra-annual period, and a plurality of intra-annual period data is obtained.

[0081] Further, the step A3, according to the division threshold, the division starting point, the division test point and the division termination point, the multi-year average daily flow sequence set is divided into intra-annual period, and a plurality of intra-annual period data is obtained, comprising:

[0082] A31, according to the division starting point and the division test point, the multi-year average daily flow sequence set is divided into first sequence, and first flow sequence is obtained;

[0083] A32, according to the division test point and the division termination point, the multi-year average daily flow sequence set is divided into second sequence, and second flow sequence is obtained;

[0084] A33, according to the first flow sequence, the second flow sequence is verified for empirical distribution consistency, and empirical distribution verification result is obtained;

[0085] A34, performing a first threshold verification on the division end point according to the division threshold, to obtain a first threshold verification result;

[0086] A35, if the experience distribution verification result is that the experience distributions are different and the first threshold verification result is that the division end point is less than the division threshold, retaining the first flow sequence, updating the division start point according to the division check point, and updating the division check point and the division end point according to the updated division start point, and then returning to perform the steps of obtaining the division start point, the division check point and the division end point;

[0087] Alternatively, A36, if the experience distribution verification result is that the experience distributions are different and the first threshold verification result is that the division end point is equal to the division threshold, obtaining a plurality of intra-annual period data according to the current second flow sequence and all the first flow sequences.

[0088] In the embodiments of the present application, for a certain cycle process, the multi-year average daily flow sequence set can be divided into sequences according to the division start point and the division check point in the current cycle process, so as to obtain a plurality of multi-year average daily flow data between the division start point and the division check point. Specifically, if the point value of the division start point is 2, the point value of the division check point is 4, and the multi-year average daily flow sequence set is {Q1, Q2, Q3, Q4, Q5}, the corresponding first flow sequence includes the second multi-year average daily flow data Q2 in the time sequence, the third multi-year average daily flow data Q3 in the time sequence and the fourth multi-year average daily flow data Q4 in the time sequence.

[0089] It can be understood that the second flow sequence is similar to the content of the aforementioned first flow sequence, and can be simply analogized. The experience distribution consistency verification in step A33 can be based on K-S (Kolmogorov-Smirnov) test method to calculate the consistency between the first flow sequence and the second flow sequence, so as to obtain the experience distribution verification result. Exemplarily, the experience distribution function of the first flow sequence can be represented as:

[0090]

[0091] wherein, is the experience distribution function; is the first flow sequence; i is the division start point; j is the division check point; t is the time point of the first flow sequence in the time sequence; is the indicator function of the first flow sequence.

[0092] It should be noted that the empirical distribution function of the second flow sequence is similar to the empirical distribution function of the first flow sequence, which can be simply analogized. After obtaining the empirical distribution functions of the first flow sequence and the second flow sequence respectively, the confidence between the two empirical distribution functions can be calculated, and the significance level between the first flow sequence and the second flow sequence can be calculated based on the obtained confidence. If the significance level is less than a preset critical threshold, an empirical distribution verification result representing different empirical distributions is generated; or if the significance level is greater than or equal to the preset critical threshold, an empirical distribution verification result representing the same empirical distribution is generated. The preset critical threshold can be any one of 0.01, 0.02, 0.05, etc. The examples in the present application are only for illustration.

[0093] It is worth mentioning that the first threshold verification can be to compare the size relationship between the division threshold and the division end point, so as to obtain the first threshold verification result. Specifically, if the empirical distribution verification result is different empirical distributions and the first threshold verification result is that the division end point is less than the division threshold, it indicates that there is no consistency between the first flow sequence and the second flow sequence, and there is multi-year average daily flow data that has not been traversed in the multi-year average daily flow sequence set. At this time, the first flow sequence of the current loop process can be retained, and then the point value of the division test point is assigned to the division start point, so as to obtain an updated division start point. Then, the division test point and the division end point are updated based on the updated division start point. Specifically, the updated division test point can be equal to the updated division end point, and the updated division end point can be equal to the sum of the updated division start point and the value 1. Finally, step A2 is returned to be executed.

[0094] It should be added that if the empirical distribution verification result is different empirical distributions and the first threshold verification result is that the division end point is equal to the division threshold, it indicates that there is no consistency between the first flow sequence and the second flow sequence, and there is no multi-year average daily flow data that has not been traversed in the multi-year average daily flow sequence set. That is, the first flow sequence and the second flow sequence are the last two flow sequences generated by the multi-year average daily flow sequence set, and the first flow sequence and the second flow sequence are the sequence set of two adjacent multi-year average daily flow data in the time sequence. At this time, each first flow sequence retained in the previous loop process can be determined as an intra-annual period data, the first flow sequence obtained in the current loop process can be determined as an intra-annual period data, and the second flow sequence obtained in the current loop process can be determined as an intra-annual period data. The number of multi-year average daily flow data in each intra-annual period data belongs to the same empirical distribution.

[0095] In some embodiments, the method further comprises:

[0096] A37, if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is less than the division threshold, the division termination point is updated, and then the steps of obtaining the division start point, the division inspection point, and the division termination point are executed again;

[0097] or,

[0098] A38, if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is equal to the division threshold, the second threshold verification is performed on the division inspection point according to the division threshold, and a second threshold verification result is obtained.

[0099] In the embodiments of the present application, if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is less than the division threshold, it indicates that the first flow sequence and the second flow sequence have consistency and the multi-year average daily flow sequence set has multi-year average daily flow data that has not been traversed, at this time, the division termination point can be updated, specifically, the point value of the division termination point can be operated by 1 to obtain an updated division termination point, and then the step A2 is executed again.

[0100] It can be understood that if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division termination point is equal to the division threshold, the size relationship between the division threshold and the division inspection point in the current loop process can be compared to obtain a second threshold verification result.

[0101] In some embodiments, the method further comprises:

[0102] A381, if the second threshold verification result is that the division inspection point is less than the division threshold, the division inspection point is updated, and the division termination point is updated according to the updated division inspection point, and then the steps of obtaining the division start point, the division inspection point, and the division termination point are executed again;

[0103] or,

[0104] A382, if the second threshold verification result is that the division inspection point is equal to the division threshold, a third flow sequence is obtained according to the current first flow sequence and the current second flow sequence, and a plurality of intra-annual period data are obtained according to the third flow sequence and all previous first flow sequences.

[0105] In the embodiment of the present application, if the second threshold verification result is that the division test point is less than the division threshold, the point value of the division test point can be operated by 1 to obtain an updated division test point; meanwhile, the updated division test point is operated by 1 again, and the point value of the division test point after being operated by 1 again is assigned to the division termination point, and then the step A2 is returned to be executed.

[0106] It can be understood that, if the second threshold verification result is that the division test point is equal to the division threshold, it indicates that the division test point has reached the division termination point, the first flow sequence and the second flow sequence are the last two flow sequences of the multi-year average daily flow sequence set, and the first flow sequence and the second flow sequence have consistency, that is, the first flow sequence and the second flow sequence belong to the same intra-year period, at this time, the current first flow sequence and the second flow sequence can be integrated, and the obtained third flow sequence is determined as an intra-year period data, and each first flow sequence reserved in the previous loop process is determined as an intra-year period data.

[0107] It is worth mentioning that, the intra-year period division of the embodiment of the present application can divide the multi-year average daily flow data into intra-year period data according to the time date of the year, so as to obtain a plurality of intra-year period data, each intra-year period data is a set of multi-year average daily flow data distinguished according to time period, the plurality of intra-year period data not only can fully reflect the distribution characteristics of the to-be-evaluated index of the river and lake, but also is beneficial to improving the adaptation degree of the subsequent generated evaluation grade and the river and lake, and thus is beneficial to improving the objectivity and accuracy of the health evaluation result of the river and lake.

[0108] In step 130, a grade evaluation table generation process is performed on all the intra-year period data to obtain a period evaluation table corresponding to each intra-year period data.

[0109] In the embodiment of the present application, a period average table corresponding to each intra-year period data can be generated based on each multi-year average daily flow data in each intra-year period data, and the sum of the intra-year periods of all the period average tables is the entire time sequence.

[0110] In some embodiments, the grade evaluation table generation process is performed on the intra-year period data to obtain a period evaluation table, including:

[0111] B1, a plurality of intra-year period index sub-data of the intra-year period data are obtained, each intra-year period index sub-data corresponds to a different river and lake health evaluation index type;

[0112] B2, data analysis is performed on all the intra-year period index sub-data to obtain an index mean and an index standard deviation of each intra-year period index sub-data;

[0113] B3, obtaining the period evaluation table according to the mean value of all the indexes and the standard deviation of all the indexes.

[0114] In the embodiment of the present application, for the period data in a certain year, first, the period in the year corresponding to all the multi-year average daily flow data of the period data in the year can be obtained, and the period in the year can be constructed by the time date of all the multi-year average daily flow data of the period data in the year; then, based on the obtained period in the year, a plurality of period in the year index sub-data corresponding to the period data in the year is determined, and each period in the year index sub-data is used to represent the river and lake health evaluation index in the corresponding period in the year, which can be the river and lake longitudinal connectivity index, ecological flow sufficiency, water plant community condition or public satisfaction of the river and lake in the corresponding period in the year.

[0115] It can be understood that for a period in the year index sub-data of the period data in the year, the index mean μ n and the index standard deviation σ n of the period in the year index sub-data in the period in the year can be calculated respectively. n n n n n n n n n n n n n Then, based on all the obtained index thresholds, a sub-evaluation table corresponding to the period in the year index sub-data is constructed, and the sub-evaluation table can be shown in Table 1 as follows:

[0116] Table 1

[0117]

[0118] It should be noted that the sub-evaluation tables of the remaining period in the year index sub-data are similar to the foregoing, which can be simply analogized. After obtaining the sub-evaluation tables of all the period in the year index sub-data of the period data in the year, all the sub-evaluation tables can be integrated to obtain the period evaluation table corresponding to the period data in the year, and the remaining period data in the year is the same, and the present application will not be repeated here.

[0119] ​​​​​​​​​​​​​​Step 140, performing health mapping evaluation on all the period evaluation tables according to the health indicator data to obtain the health evaluation result of the river and lake.

[0120] In the embodiments of the present application, the health evaluation result of the river and lake at a certain time date can be obtained by table lookup scoring based on the health indicator data of the river and lake at the time date and all the period evaluation tables.

[0121] In some embodiments, the step 140 of performing health mapping evaluation on all the period evaluation tables according to the health indicator data to obtain the health evaluation result of the river and lake comprises:

[0122] C1, performing table screening on all the period evaluation tables according to the health indicator data to obtain a target evaluation table corresponding to the health indicator data;

[0123] C2, performing index scoring evaluation on the health indicator data according to the target evaluation table to obtain the health evaluation result of the river and lake.

[0124] In the embodiments of the present application, the period evaluation table corresponding to the health indicator data can be determined by judging whether the time date attribute in the health indicator data belongs to the in-year period of the period evaluation table, based on the time date attribute in the health indicator data and the in-year period corresponding to each period evaluation table, and the period evaluation table corresponding to the health indicator data is determined as the target evaluation table.

[0125] It can be understood that, in the embodiments of the present application, the health indicator data includes the river longitudinal connectivity index and the ecological flow sufficiency, and the step C2 can be to respectively obtain the grade score of the river longitudinal connectivity index on the corresponding sub evaluation table and the grade score of the ecological flow sufficiency on the corresponding sub evaluation table, and to perform accumulation operation on all the obtained grade scores to obtain the health evaluation result of the river and lake.

[0126] A river and lake health evaluation system according to the embodiments of the present application is described in detail below with reference to the accompanying drawings.

[0127] Referring to Figure 2 The river and lake health evaluation system according to the embodiments of the present application comprises:

[0128] A first processing unit 101 is configured to obtain health indicator data of a river and lake and a multi-year average daily flow sequence set;

[0129] A second processing unit 102 is configured to divide the multi-year average daily flow sequence set into in-year periods to obtain a plurality of in-year period data, and the experience distribution of the in-year period data adjacent in time is different;

[0130] The third processing unit 103 is configured to generate a periodicity evaluation table for each of the intra-year periodicity data.

[0131] The fourth processing unit 104 is configured to evaluate the health of the river or lake according to the health index data and the periodicity evaluation table.

[0132] It can be understood that the content in the above method embodiments is applicable to the system embodiments, the system embodiments specifically implement the functions of the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0133] With reference to Figure 3 The embodiment of the present application further provides an electronic device, which comprises:

[0134] at least one processor 201;

[0135] at least one memory 202, configured to store at least one program;

[0136] When the at least one program is executed by the at least one processor 201, the at least one processor 201 implements the above method embodiments.

[0137] Similarly, it can be understood that the content in the above method embodiments is applicable to the device embodiments, the device embodiments specifically implement the functions of the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0138] The embodiment of the present application further provides a computer readable storage medium, which stores a program executable by the processor 201, and the program executable by the processor 201 is used to implement the above method embodiments when executed by the processor 201.

[0139] Similarly, the content in the above method embodiments is applicable to the computer readable storage medium embodiments, the computer readable storage medium embodiments specifically implement the functions of the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0140] In alternative embodiments, the functions / operations in the flow diagrams can occur in sequences other than those depicted. For example, two operations shown in succession can in fact be executed substantially concurrently or the operations sometimes can be executed in the reverse order depending upon the functionality / operations involved. Also, embodiments presented and described in this application are provided by way of example only. The

[0141] Moreover, while this application has been described in the context of functional modules, it is to be understood that one or more of the functions and / or features can be integrated in a single physical device and / or software module or one or more functions and / or features can be implemented in separate physical devices or software modules. It will also be appreciated that detailed discussion of the actual implementation of each module is not necessary for an understanding of the application. Rather, the actual implementation is a matter of choice apart from the conception and the disclosure provided that the attributes of each module remain intact. Accordingly, the disclosure is to be understood in a manner that

[0142] If the functions are implemented in software, the functions can be stored in or implemented as one or more software modules on a computer-readable medium. Based on the understanding thus far, the technical solutions of the present application, in essence, or the parts that make contributions to the prior art, or parts of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a number of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the embodiments of the method. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk, and various other media that can store program codes.

[0143] The logic and / or steps represented in the flow diagrams or otherwise described herein, for example, can be embodied in non-transitory computer- readable media, executed by an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions, or in conjunction with which the instructions can be executed. In the context of this specification, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The computer-readable medium can be, for example but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, device, or propagation medium.

[0144] More specific examples (a non-exhaustive list) of the computer-readable medium include the following: an electrical connection (electronic) having one or more wires, a portable computer diskette (magnetic), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium can also be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, for example, via optical scanning of the paper or other medium, then compiled, interpreted, or otherwise processed in a suitable manner, if necessary, and stored in a computer memory.

[0145] It should be understood that aspects of the present application can be implemented in hardware, software, firmware, or combinations thereof. In the above embodiments, various steps or methods can be implemented in software or firmware that is stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any of the following technologies, known in the art, or combinations thereof, can be used: a discrete logic circuit having logic gates for implementing logic functions upon data signals, an application specific integrated circuit having appropriate combinational logic gates, a programmable gate array (PGA), a field programmable gate array (FPGA), and / or the like.

[0146] In the above description of the present specification, reference is made to the descriptions of the terms "one embodiment / one example", "another embodiment / another example", or "some embodiments / some examples" and the like, meaning that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. The illustrative representations of the above terms in the present specification do not necessarily refer to the same embodiment or example. Also, the described particular features, structures, materials, or characteristics can be combined in any appropriate manner in one or more embodiments or examples.

[0147] While the embodiments of the application have been shown and described, it is to be understood that the embodiments can be varied, modified, substituted and changed by those skilled in the art without departing from the principles and spirit of the application, the scope of which is defined by the claims and their equivalents.

[0148] The above is a specific description of the preferred embodiments of the application, but the application is not limited to the embodiments, and those skilled in the art can make various equivalent modifications or substitutions without departing from the spirit of the application, and these equivalent modifications or substitutions are all included in the scope defined by the claims of the application.

Claims

1. A method for evaluating the health of a river or lake, characterized by, The method comprises the following steps: obtaining health index data and a multi-year average daily flow sequence set of a river or lake; dividing the multi-year average daily flow sequence set into intra-annual cycles to obtain a plurality of intra-annual cycle data, and the empirical distribution of time-adjacent intra-annual cycle data is different; generating a ranking evaluation table for all the intra-annual cycle data to obtain a cycle evaluation table corresponding to each intra-annual cycle data; performing health mapping evaluation on all the cycle evaluation tables according to the health index data to obtain a health evaluation result of the river or lake; the step of dividing the multi-year average daily flow sequence set into intra-annual cycles to obtain a plurality of intra-annual cycle data comprises the following steps: obtaining a division threshold; obtaining a division starting point, a division test point and a division ending point; dividing the multi-year average daily flow sequence set into intra-annual cycles according to the division threshold, the division starting point, the division test point and the division ending point to obtain a plurality of intra-annual cycle data; the step of dividing the multi-year average daily flow sequence set into intra-annual cycles according to the division threshold, the division starting point, the division test point and the division ending point to obtain a plurality of intra-annual cycle data comprises the following steps: performing first sequence division on the multi-year average daily flow sequence set according to the division starting point and the division test point to obtain a first flow sequence; performing second sequence division on the multi-year average daily flow sequence set according to the division test point and the division ending point to obtain a second flow sequence; performing empirical distribution consistency verification on the second flow sequence according to the first flow sequence to obtain an empirical distribution verification result; performing first threshold verification on the division ending point according to the division threshold to obtain a first threshold verification result; if the empirical distribution verification result is that the empirical distributions are different and the first threshold verification result is that the division ending point is less than the division threshold, the first flow sequence is retained, the division starting point is updated according to the division test point, and the updated division starting point is used to update the division test point and the division ending point, and then the step of obtaining the division starting point, the division test point and the division ending point is performed again; or, if the empirical distribution verification result is that the empirical distributions are different and the first threshold verification result is that the division ending point is equal to the division threshold, a plurality of intra-annual cycle data are obtained according to the current second flow sequence and all the first flow sequences.

2. The method of claim 1, wherein, The method further comprises the following steps: if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division ending point is less than the division threshold, the division ending point is updated, and then the step of obtaining the division starting point, the division test point and the division ending point is performed again; or, if the empirical distribution verification result is that the empirical distributions are the same and the first threshold verification result is that the division ending point is equal to the division threshold, second threshold verification is performed on the division test point according to the division threshold to obtain a second threshold verification result.

3. The method of claim 2, wherein, The method further comprises the following steps: If the second threshold verification result is that the division test point is less than the division threshold, the division test point is updated, and the division end point is updated according to the updated division test point, and then the division start point, the division test point and the division end point are obtained again; Or, If the second threshold verification result is that the division test point is equal to the division threshold, a third flow sequence is obtained according to the current first flow sequence and the current second flow sequence, and a plurality of intra-annual cycle data are obtained according to the third flow sequence and all previous first flow sequences.

4. The method of claim 1, wherein, A grade evaluation table generation process is performed on the intra-annual cycle data to obtain a cycle evaluation table, including: Obtaining a plurality of intra-annual cycle index sub-data of the intra-annual cycle data, each of which corresponds to a different river and lake health evaluation index type; Data analysis is performed on all the intra-annual cycle index sub-data to obtain the index mean and index standard deviation of each intra-annual cycle index sub-data; According to all the index mean and all the index standard deviation, the cycle evaluation table is obtained.

5. The method of claim 1, wherein, According to the health index data, the health mapping evaluation of all the cycle evaluation tables is performed to obtain the health evaluation result of the river and lake, including: According to the health index data, a table screening is performed on all the cycle evaluation tables to obtain a target evaluation table corresponding to the health index data; According to the target evaluation table, an index scoring evaluation is performed on the health index data to obtain the health evaluation result of the river and lake.

6. A river and lake health evaluation system characterized by, Including: A first processing unit is configured to obtain health index data of a river and lake and a multi-year average daily flow sequence set; A second processing unit is configured to perform intra-annual cycle division on the multi-year average daily flow sequence set to obtain a plurality of intra-annual cycle data, the experience distribution of which is different between adjacent intra-annual cycle data in time; A third processing unit is configured to perform a grade evaluation table generation process on all the intra-annual cycle data to obtain a cycle evaluation table corresponding to each intra-annual cycle data; A fourth processing unit is configured to perform health mapping evaluation on all the cycle evaluation tables according to the health index data to obtain a health evaluation result of the river and lake; The intra-annual cycle division on the multi-year average daily flow sequence set to obtain a plurality of intra-annual cycle data, including: Obtaining a division threshold; Obtaining a division start point, a division test point and a division end point; According to the division threshold, the division start point, the division test point and the division end point, the intra-annual cycle division is performed on the multi-year average daily flow sequence set to obtain a plurality of intra-annual cycle data; The intra-annual cycle division on the multi-year average daily flow sequence set to obtain a plurality of intra-annual cycle data according to the division threshold, the division start point, the division test point and the division end point, including: According to the division start point and the division test point, a first sequence division is performed on the multi-year average daily flow sequence set to obtain a first flow sequence; According to the division test point and the division termination point, the second flow sequence set is subjected to second sequence division to obtain a second flow sequence; According to the first flow sequence, the second flow sequence is subjected to empirical distribution consistency verification to obtain an empirical distribution verification result; According to the division threshold value, the division termination point is subjected to first threshold value verification to obtain a first threshold value verification result; If the empirical distribution verification result is that the empirical distributions are different and the first threshold value verification result is that the division termination point is less than the division threshold value, the first flow sequence is retained, the division starting point is updated according to the division test point, and the updated division starting point is used to update the division test point and the division termination point, and then the step of obtaining the division starting point, the division test point and the division termination point is executed again; or, if the empirical distribution verification result is that the empirical distributions are different and the first threshold value verification result is that the division termination point is equal to the division threshold value, a plurality of annual periodic data are obtained according to the current second flow sequence and all the first flow sequences.

7. An electronic device, comprising: Comprise: At least one processor; At least one memory for storing at least one program; When the at least one program is executed by the at least one processor, the at least one processor implements the method of any one of claims 1-5.

8. A computer readable storage medium having stored therein a program which is executable by a processor, characterized in that, The program executable by the processor when executed by the processor is used to implement the method of any one of claims 1-5.

Images