Basin ecological flow accounting method and system

Abstract

The application discloses a method and system for calculating ecological flow of a river basin, and relates to the technical field of ecological flow calculation, which comprises the following steps: obtaining a scatter plot by screening ecological flow; dividing a year into a first ecological period, a second ecological period and a third ecological period based on the scatter plot; obtaining a functional relationship between ecological flow and months in the first ecological period, the second ecological period and the third ecological period, respectively, and marking the functional relationships as a first function, a second function and a third function; obtaining a first threshold function, a second threshold function and a third threshold function based on the ecological flow of the second number of historical months, the first function, the second function and the third function; and obtaining an ecological flow calculation result of a corresponding month based on the first threshold function, the second threshold function and the third threshold function; the application is used to solve the problem that the setting of minimum ecological flow in the prior art is relatively subjective, resulting in an inaccurate ecological flow calculation result.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ecological flow accounting, in particular to a watershed ecological flow accounting method and system. BACKGROUND

[0002] With the rapid development of social economy, the intensity of human development and utilization of water resources is increasing, and the allocation pattern of water resources in the watershed has changed significantly. As an important part of the ecological system, the ecological function of rivers depends on certain water flow conditions, and ecological flow is a key factor to ensure the health of river ecosystems.

[0003] With the continuous deepening of research in the fields of ecology, hydrology and water resources management, various ecological flow accounting methods have emerged, including hydrological methods, hydraulic methods, habitat simulation methods, and comprehensive methods. In the hydrological method, the flow threshold curve method takes the flow at 90% frequency as the minimum ecological flow in the river channel. The setting of 90% is subjective and cannot generate appropriate threshold values based on historical ecological flow. For example, the patent application with the application publication number CN113887965A discloses a watershed ecological flow accounting method. In this scheme, the conventional flow threshold curve method is used, which takes the flow at 90% as the minimum ecological flow in the river channel. The setting of 90% is subjective, and the setting of the minimum ecological flow in the existing ecological flow accounting technology is also subjective, resulting in inaccurate ecological flow accounting results. SUMMARY

[0004] The present application aims to at least solve one of the technical problems in the prior art. By obtaining ecological flow accounting results for corresponding months based on first, second and third threshold functions, the problem of subjective setting of minimum ecological flow in existing ecological flow accounting technology is solved, resulting in inaccurate ecological flow accounting results.

[0005] To achieve the above-mentioned purpose, the present application provides a watershed ecological flow accounting method, comprising the following steps:

[0006] Screening the ecological flow of the first number of historical months to obtain screened ecological flow;

[0007] Obtaining a scatter plot based on the screened ecological flow;

[0008] Dividing a year into a first ecological period, a second ecological period and a third ecological period based on the scatter plot;

[0009] Obtaining the functional relationship between ecological flow and month in the first ecological period, the second ecological period and the third ecological period, respectively, and marking them as first function, second function and third function;

[0010] The first threshold function, the second threshold function and the third threshold function are obtained based on the ecological flow of the second number of historical months, the first function, the second function and the third function;

[0011] The ecological flow accounting result of the corresponding month is obtained based on the first threshold function, the second threshold function and the third threshold function;

[0012] Further, the ecological flow of the first number of historical months is screened to obtain screened ecological flow, including the following sub-steps:

[0013] The historical ecological flow of the same month is sorted from small to large and marked as Ss1 to Ss i ;

[0014] It is judged whether 0.5*Xt is an integer, if yes, Ss （0.5*Xt） is taken as the intermediate position value; if not, the integers on the left and right sides of 0.5*Xt are obtained and marked as Zz and Zy, and the mean of Ss （Zz） and Ss （Zy） is taken as the intermediate position value; wherein Xt is the number of ecological flow of the same historical month;

[0015] It is judged whether a*Xt is an integer, if yes, Ss （a*Xt） is taken as the first proportion value; if not, the integers on the left and right sides of a*Xt are obtained and marked as Dz and Dy, and the mean of Ss （Dz） and Ss （Dy） is taken as the first proportion value; wherein a is an ecological flow coefficient, and the range of a is (0.5, 1).

[0016] Further, the ecological flow of the first number of historical months is screened to obtain screened ecological flow, including the following sub-steps:

[0017] The screened ecological threshold is calculated as: ; wherein Cs is the screened ecological threshold, Bz is the intermediate position value, and B1 is the first proportion value;

[0018] The historical ecological flow greater than the screened ecological threshold is deleted to obtain the screened ecological flow.

[0019] Further, the division of the scatter diagram is based on the screened ecological flow, including the following sub-steps:

[0020] The screened ecological flow of each historical month is obtained;

[0021] The mean of the screened ecological flow of each historical month is calculated and marked as the historical ecological flow mean;

[0022] Sort the historical ecological flow mean in ascending order, and mark the sorted result of the historical ecological flow mean in ascending order as the corresponding sequence number, wherein the sequence number is a positive integer from 1 to n;

[0023] Establish a coordinate system with the sequence number as the horizontal axis data and the historical ecological flow mean as the vertical axis data, and mark it as a division coordinate system;

[0024] Draw the division coordinate system with all sequence numbers and corresponding historical ecological flow means as division coordinate points to obtain a division scatter plot.

[0025] Further, based on the division scatter plot, the year is divided into a first ecological period, a second ecological period and a third ecological period, including the following sub-steps:

[0026] Connect adjacent division coordinate points from left to right according to the division scatter plot to obtain a line segment, mark it as a first connection line segment, and obtain the slope of the first connection line segment, mark it as F m , wherein m is a positive integer, and m represents the order of the connection line segment from small to large;

[0027] Calculate the difference between adjacent F m , which is FC m =F m -F m-1 : wherein FC m is the difference between adjacent F m , and m is greater than or equal to 2 at this time;

[0028] Sort FC m from large to small, and obtain the corresponding F m of the first two FC m , mark it as a division month threshold; mark the first division month threshold and the second division month threshold according to the small and large values of the division month threshold, respectively;

[0029] In the division scatter plot, mark the range of sequence numbers less than the first division month threshold as the first ecological period; mark the range of sequence numbers greater than or equal to the first division month threshold and less than the second division month threshold as the second ecological period; and mark the range of sequence numbers greater than or equal to the second division month threshold as the third ecological period.

[0030] Further, the functional relationship between the ecological flow and the month in the first ecological period, the second ecological period and the third ecological period is obtained respectively, and is marked as the first function, the second function and the third function, including the following sub-steps:

[0031] The divided coordinate points are functionally fitted to obtain a first function in a first ecological period; the divided coordinate points are functionally fitted to obtain a second function in a second ecological period; and the divided coordinate points are functionally fitted to obtain a third function in a third ecological period.

[0032] Further, obtaining the first threshold function, the second threshold function and the third threshold function based on the ecological flow of the second number of historical months, the first function, the second function and the third function comprises the following sub-steps:

[0033] The ecological flow corresponding to each serial number in the first function is obtained and marked as a first function predicted flow; the difference between the first function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of the difference between all first function predicted flows and the corresponding historical month ecological flow is obtained and marked as a first difference threshold value.

[0034] The ecological flow corresponding to each serial number in the second function is obtained and marked as a second function predicted flow; the difference between the second function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of the difference between all second function predicted flows and the corresponding historical month ecological flow is obtained and marked as a second difference threshold value.

[0035] The ecological flow corresponding to each serial number in the third function is obtained and marked as a third function predicted flow; the difference between the third function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of the difference between all third function predicted flows and the corresponding historical month ecological flow is obtained and marked as a third difference threshold value.

[0036] Further, obtaining the first threshold function, the second threshold function and the third threshold function based on the ecological flow of the second number of historical months, the first function, the second function and the third function further comprises the following sub-steps:

[0037] The first function is moved in the negative direction of the vertical axis by a distance of the first difference threshold value to obtain the first threshold function; the second function is moved in the negative direction of the vertical axis by a distance of the second difference threshold value to obtain the second threshold function; and the third function is moved in the negative direction of the vertical axis by a distance of the third difference threshold value to obtain the third threshold function.

[0038] Further, obtaining the ecological flow accounting result of the corresponding month based on the first threshold function, the second threshold function and the third threshold function comprises the following steps:

[0039] The historical ecological flow corresponding to each serial number in the first threshold function, the second threshold function and the third threshold function is obtained as the ecological flow accounting result of the corresponding month and output.

[0040] The present application also provides a watershed ecological flow accounting system, including: a data screening module, a scatter plot drawing module, an ecological period division module, a function fitting module, a threshold acquisition module and an ecological flow accounting module;

[0041] The data screening module is used to screen the ecological flows of a first number of historical months to obtain screened ecological flows;

[0042] The scatter plot drawing module is used to obtain a partition scatter plot based on the screening ecological flow;

[0043] The ecological period division module is used to divide a year into a first ecological period, a second ecological period and a third ecological period based on the division scatter plot;

[0044] The function fitting module is used to obtain the functional relationship between the ecological flow and the month during the first ecological period, the second ecological period, and the third ecological period, respectively marked as the first function, the second function, and the third function;

[0045] The threshold acquisition module is used to acquire a first threshold function, a second threshold function and a third threshold function based on the ecological flow of a second number of historical months, the first function, the second function and the third function;

[0046] The ecological flow accounting module is used to obtain the ecological flow accounting result of the corresponding month based on the first threshold function, the second threshold function and the third threshold function.

[0047] Beneficial effects of the present invention: The present invention obtains ecological flow accounting results for corresponding months based on the first threshold function, the second threshold function, and the third threshold function. The advantage is that the accounting results can be obtained based on historical ecological flow, and the obtained accounting results are more accurate than the conventional flow threshold curve method.

[0048] The present invention divides a year into the first ecological period, the second ecological period and the third ecological period based on a scatter plot. The advantage is that due to the seasonal changes each year, the ecological flow and the corresponding water demand of each month are different. Therefore, different ecological periods are divided to make the calculation results more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 is a functional block diagram of the system of the present invention;

[0050] Figure 2 A schematic diagram of a partitioned scatter plot of the present invention;

[0051] Figure 3 Schematic diagram of the first month-division threshold and the second month-division threshold of the present invention;

[0052] Figure 4 Schematic diagram of the first function, the second function and the third function of the present application;

[0053] Figure 5 Schematic diagram of the first threshold function, the second threshold function and the third threshold function of the present application;

[0054] Figure 6 Flow chart of the steps of the method of the present application. DETAILED DESCRIPTION

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

[0056] Embodiment 1, please refer to Figure 1 As shown in the figure, the present application provides a watershed ecological flow accounting system, which comprises a data screening module, a scatter plot drawing module, an ecological period division module, a function fitting module, a threshold acquisition module and an ecological flow accounting module;

[0057] The data screening module is configured to screen the ecological flow of the first number of historical months to obtain screened ecological flow; the average ecological flow of each month in the past ten years is selected, and the data in the past ten years is selected, i.e., the first number is 120, and the ecological flow of the historical months is the average ecological flow in the month;

[0058] The data screening module is configured with a first proportion value and an intermediate position value acquisition strategy, and the first proportion value and the intermediate position value acquisition strategy comprises:

[0059] The historical ecological flow of the same month is sorted from small to large and marked as Ss1 to Ss i ;

[0060] It is judged whether 0.5*Xt is an integer, if yes, Ss （0.5*Xt） is taken as the intermediate position value; if not, the integers on the left and right sides of 0.5*Xt are obtained and marked as Zz and Zy, and the mean of Ss （Zz） and Ss （Zy） is taken as the intermediate position value; wherein Xt is the number of ecological flow of the same historical month;

[0061] It is judged whether a*Xt is an integer, if yes, Ss （a*Xt） is taken as the first proportion value; if not, the integers on the left and right sides of a*Xt are obtained and marked as Dz and Dy, and the mean of Ss （Dz） and Ss （Dy）The mean of is taken as the first proportional value; where a is the ecological flow coefficient, and the range of a is (0.5, 1); an ecological flow with a large historical distribution is selected, so the range of a is (0.5, 1), and a is preferably set to 0.75;

[0062] In practical applications, since the data is obtained in the past ten years, that is, Xt is 10, taking the serial number 6 as an example, it is determined that 0.5*10 is an integer, and Ss （5） =225m 3 / s is used as the middle position value, and it is determined that 0.75*10 is not an integer. （7） =231m 3 / sSs （8） =233m 3 The average value of / s is 232m 3 / s is the first proportional value.

[0063] The data filtering module is configured with data filtering strategies, which include:

[0064] The screening ecological threshold is obtained as: ; Where Cs is the screening ecological threshold, Bz is the middle position value, and B1 is the first proportion value; the data is regarded as uniformly distributed, and the maximum average ecological flow under the uniform distribution is obtained based on the middle position value and the first proportion value. If it is greater than this value.

[0065] The filtered ecological flow is obtained by deleting the historical ecological flow that is greater than the screening ecological threshold. Since the ecological flow accounting is to obtain the minimum ecological flow to maintain the normal operation of the ecology, the larger part of the data is deleted.

[0066] In practical application, taking sequence number 6 as an example, the screening ecological threshold is: Cs=239m 3 / s.

[0067] The scatter plot drawing module is used to obtain the partition scatter plot based on the screening ecological flow;

[0068] The scatter plot drawing module is configured with a scatter plot drawing strategy, which includes:

[0069] Get the filtered ecological flow for each historical month;

[0070] Calculate the mean of the filtered ecological flow for each historical month and mark it as the historical ecological flow mean;

[0071] Sort the historical ecological flow means from small to large, and mark the results of the sorting of the historical ecological flow means from small to large as corresponding serial numbers, where the serial numbers are positive integers from 1 to n;

[0072] With the serial number as the horizontal axis data and the historical ecological flow mean as the vertical axis data, a coordinate system is established and marked as the division coordinate system;

[0073] All serial numbers and the corresponding historical ecological flow mean are used as division coordinate points to draw a division coordinate system to obtain a division scatter plot;

[0074] In practical applications, please refer to Figure 2 As shown, draw a partition scatter plot;

[0075] The ecological period division module is used to divide a year into a first ecological period, a second ecological period and a third ecological period based on the division scatter plot;

[0076] The ecological period division module is configured with ecological interval division strategies, which include:

[0077] Connect the adjacent division coordinate points from left to right according to the division scatter plot to obtain a line segment, marked as the first connection line segment, and obtain the slope of the first connection line segment, marked as F m , where m is a positive integer, and m represents the order of connecting line segments from small to large;

[0078] Calculate adjacent F m The difference is: FC m =F m -F m-1 :FC m For adjacent F m The difference between , at this time m is greater than or equal to 2;

[0079] FC m Sort from largest to smallest and get the first two FCs m Corresponding F m , marked as the monthly division threshold; the monthly division threshold is marked as the first monthly division threshold and the second monthly division threshold according to the smaller and larger values;

[0080] In the partitioned scatter plot, the range of serial numbers less than the first partition month threshold is marked as the first ecological period; the range of serial numbers greater than or equal to the first partition month threshold and less than the second partition month threshold is marked as the second ecological period; the range of serial numbers greater than or equal to the second partition month threshold is marked as the third ecological period. This is because the ecological flow in different seasons of each year is different, for example, the ecological flow is higher in summer and lower in winter, and the ecological flow demand obtained in different seasons is also different.

[0081] In practical applications, please refer to Figure 3As shown, the obtained first division month threshold and second division month threshold are 4 and 9 respectively, so the range of the serial number less than 4 is marked as the first ecological period; the range of the serial number greater than or equal to 4 and less than 9 is marked as the second ecological period; the range of the serial number greater than or equal to 9 is marked as the third ecological period, and the segmented analysis of the ecological flow calculation result is more accurate.

[0082] The function fitting module is configured to obtain the function relationship between the ecological flow and the month in the first ecological period, the second ecological period and the third ecological period, respectively, and mark them as the first function, the second function and the third function.

[0083] In practical application, please refer to Figure 4 As shown, the first function, the second function and the third function are obtained.

[0084] The function fitting module is configured with a function fitting strategy, and the function fitting strategy includes:

[0085] The first function is obtained by function fitting of the division coordinate points in the first ecological period; the second function is obtained by function fitting of the division coordinate points in the second ecological period; and the third function is obtained by function fitting of the division coordinate points in the third ecological period.

[0086] The threshold acquisition module is configured to obtain the first threshold function, the second threshold function and the third threshold function based on the ecological flow of the second number of historical months, the first function, the second function and the third function.

[0087] The threshold acquisition module is configured with a threshold acquisition strategy, and the threshold acquisition strategy includes:

[0088] The ecological flow corresponding to each serial number in the first function is obtained and marked as the first function predicted flow; the difference between the first function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of all the first function predicted flow and the corresponding historical month ecological flow difference is obtained and marked as the first difference threshold.

[0089] The ecological flow corresponding to each serial number in the second function is obtained and marked as the second function predicted flow; the difference between the second function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of all the second function predicted flow and the corresponding historical month ecological flow difference is obtained and marked as the second difference threshold.

[0090] The ecological flow corresponding to each serial number in the third function is obtained and marked as the third function predicted flow; the difference between the third function predicted flow and the corresponding historical month ecological flow is calculated; and the maximum value of all the third function predicted flow and the corresponding historical month ecological flow difference is obtained and marked as the third difference threshold.

[0091] In practical application, please refer to Figure 5 As shown in the second function, the second function predicts the flow as 239m 3 / s, and the ecological flow of a historical month is 230m 3 / s, the difference is 9m 3 / s, and the difference between the second function prediction flow and the corresponding historical month ecological flow is 9m 3 / s, the second difference threshold is 9m 3 / s, and the second difference threshold is 12m 3 / s, and the third difference threshold is 16m 3 / s.

[0092] The threshold value acquisition module is configured with a threshold function acquisition strategy, and the threshold function acquisition strategy includes:

[0093] The first function is moved in the negative direction of the vertical axis by a distance of the first difference threshold to obtain a first threshold function; the second function is moved in the negative direction of the vertical axis by a distance of the second difference threshold to obtain a second threshold function; and the third function is moved in the negative direction of the vertical axis by a distance of the third difference threshold to obtain a third threshold function.

[0094] The ecological flow accounting module is configured to obtain the ecological flow accounting result of the corresponding month based on the first threshold function, the second threshold function and the third threshold function.

[0095] The ecological flow accounting module is configured with an ecological flow accounting strategy, and the ecological flow accounting strategy includes:

[0096] The historical ecological flow corresponding to each sequence number in the first threshold function, the second threshold function and the third threshold function is obtained as the ecological flow accounting result of the corresponding month.

[0097] In practical application, please refer to Figure 5 As shown in the figure, the historical ecological flow corresponding to each sequence number in the first threshold function, the second threshold function and the third threshold function is obtained as the ecological flow accounting result of the corresponding month, for example, when the sequence number is 1, the first threshold function corresponds to 98m 3 / s, and the month corresponding to the sequence number is December, so the ecological flow of December is 98m 3 / s, and the ecological flow accounting result of the basin in December is 98m 3 / s, and the ecological flow accounting result of the basin in December is 98m

[0098] Embodiment 2, please refer to Figure 6 As shown in the figure, the present application provides a method for accounting the ecological flow of a basin, which includes the following steps:

[0099] Step S1, screening the ecological flow of the first number of historical months to obtain a screened ecological flow; Step S1 includes the following sub-steps:

[0100] Step S101, sorting the historical ecological flow of the same month from small to large and marking them as Ss1 to Ss i ;

[0101] Step S102, determining whether 0.5*Xt is an integer, if yes, taking Ss （0.5*Xt） as the intermediate position value; if not, obtaining the integers on the left and right sides of 0.5*Xt, respectively marked as Zz and Zy, and taking the mean of Ss （Zz） and Ss （Zy） as the intermediate position value; wherein Xt is the number of ecological flows of the same historical month;

[0102] Step S103, determining whether a*Xt is an integer, if yes, taking Ss （a*Xt） as the first proportion value; if not, obtaining the integers on the left and right sides of a*Xt, respectively marked as Dz and Dy, and taking the mean of Ss （Dz） and Ss （Dy） as the first proportion value; wherein a is the ecological flow coefficient, and the range of a is (0.5, 1);

[0103] Step S104, taking the screened ecological threshold as: ; wherein Cs is the screened ecological threshold, Bz is the intermediate position value, and B1 is the first proportion value;

[0104] Step S105, deleting the historical ecological flow greater than the screened ecological threshold to obtain the screened ecological flow.

[0105] Step S2, obtaining a scatter plot based on the screened ecological flow; Step S2 includes the following sub-steps:

[0106] Step S201, obtaining the screened ecological flow of each historical month;

[0107] Step S202, taking the mean of the screened ecological flow of each historical month, marked as the historical ecological flow mean;

[0108] Step S203, sorting the historical ecological flow mean from small to large, and marking the sorting result of the historical ecological flow mean from small to large as the corresponding sequence number, wherein the sequence number is a positive integer from 1 to n;

[0109] Step S204, taking the sequence number as the horizontal axis data and the historical ecological flow mean as the vertical axis data to establish a coordinate system, marked as a division coordinate system;

[0110] Step S205, all sequence numbers and corresponding historical ecological flow averages are plotted as division coordinate points in a division coordinate system to obtain a division scatter plot.

[0111] Step S3, dividing a year into a first ecological period, a second ecological period and a third ecological period based on the division scatter plot; step S3 includes the following sub-steps:

[0112] Step S301, connecting adjacent division coordinate points from left to right according to the division scatter plot to obtain a line segment, marked as a first connection line segment, obtaining a slope of the first connection line segment, marked as F m , wherein m is a positive integer, and m from small to large represents the front and back order of the connection line segments;

[0113] Step S302, calculating the difference between adjacent F m is: FC m = F m - F m-1 , wherein FC m is the difference between adjacent F m , and m takes a value greater than or equal to 2 at this time;

[0114] Step S303, sorting FC m from large to small, obtaining the corresponding F m of the first two FC m , marked as a division month threshold; the division month threshold is marked as a first division month threshold and a second division month threshold according to the small and large values respectively;

[0115] Step S304, in the division scatter plot, the range of sequence numbers less than the first division month threshold is marked as the first ecological period; the range of sequence numbers greater than or equal to the first division month threshold and less than the second division month threshold is marked as the second ecological period; the range of sequence numbers greater than or equal to the second division month threshold is marked as the third ecological period.

[0116] Step S4, obtaining the functional relationship between ecological flow and month in the first ecological period, the second ecological period and the third ecological period respectively, marked as a first function, a second function and a third function respectively; step S4 includes the following sub-steps:

[0117] Step S401, function fitting is performed on the division coordinate points in the first ecological period to obtain the first function; function fitting is performed on the division coordinate points in the second ecological period to obtain the second function; function fitting is performed on the division coordinate points in the third ecological period to obtain the third function.

[0118] In step S5, the first threshold function, the second threshold function and the third threshold function are obtained based on the ecological flow of the second number of historical months, the first function, the second function and the third function. Step S5 includes the following sub-steps:

[0119] In step S501, the ecological flow corresponding to each serial number in the first function is obtained and marked as the first function predicted flow. The difference between the first function predicted flow and the corresponding historical month ecological flow is calculated. The maximum value of the difference between all first function predicted flows and the corresponding historical month ecological flows is obtained and marked as the first difference threshold value.

[0120] In step S502, the ecological flow corresponding to each serial number in the second function is obtained and marked as the second function predicted flow. The difference between the second function predicted flow and the corresponding historical month ecological flow is calculated. The maximum value of the difference between all second function predicted flows and the corresponding historical month ecological flows is obtained and marked as the second difference threshold value.

[0121] In step S503, the ecological flow corresponding to each serial number in the third function is obtained and marked as the third function predicted flow. The difference between the third function predicted flow and the corresponding historical month ecological flow is calculated. The maximum value of the difference between all third function predicted flows and the corresponding historical month ecological flows is obtained and marked as the third difference threshold value.

[0122] In step S504, the first function is moved in the negative direction of the vertical axis by a distance of the first difference threshold value to obtain the first threshold function. The second function is moved in the negative direction of the vertical axis by a distance of the second difference threshold value to obtain the second threshold function. The third function is moved in the negative direction of the vertical axis by a distance of the third difference threshold value to obtain the third threshold function.

[0123] In step S6, the ecological flow accounting result of the corresponding month is obtained based on the first threshold function, the second threshold function and the third threshold function. Step S6 includes the following sub-steps:

[0124] In step S601, the historical ecological flow corresponding to each serial number in the first threshold function, the second threshold function and the third threshold function is obtained as the ecological flow accounting result of the corresponding month and output.

[0125] Those skilled in the art will appreciate that embodiments of the present application can be readily used as a method, a system or a computer program product. Accordingly, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application can take the form of a computer program product on one or more computer-usable storage media (or computer- readable storage media) having computer-usable program code embodied in the medium. The medium can be any available storage media that can be accessed by a computer. By way of example, and not limitation, such computer-usable storage media can include RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other storage medium that can be used to carry or store desired computer program code in the form of instructions or data structures and that can be accessed by a computer. Also, the present application can be embodied in a computer program product that can be traded as goods or merchandise, through the storage medium described above or any other suitable medium. When the computer program code is executed by a computer, a series of instructions or a series of operations is performed. Also, the present application can be embodied in a computer program product that can be traded as goods or merchandise, through the storage medium described above or any other suitable medium. When the computer program code is executed by a computer, a series of instructions or a series of operations is performed. Figure 1 one or more flows and / or blocks Figure 1 one or more flows and / or blocks

[0126] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, another division manner can be used. For example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units or components shown or discussed can be indirect coupling or communication connection through some communication interfaces. The coupling or communication connection can be electrical, mechanical or in other forms.

Claims

1. The watershed ecological flow accounting method is characterized by: The steps include: Filter the ecological flows of the first number of historical months to obtain the filtered ecological flows; Based on the screening ecological flow, a scatter plot of the partition is obtained; Based on the partition scatter plot, a year is divided into the first ecological period, the second ecological period, and the third ecological period; Obtain the functional relationship between ecological flow and month during the first ecological period, the second ecological period, and the third ecological period, respectively, and mark them as the first function, the second function, and the third function; Obtaining a first threshold function, a second threshold function, and a third threshold function based on the ecological flow of a second number of historical months, the first function, the second function, and the third function; Obtaining ecological flow accounting results for corresponding months based on the first threshold function, the second threshold function, and the third threshold function; Acquiring the first threshold function, the second threshold function, and the third threshold function based on the ecological flow of the second number of historical months, the first function, the second function, and the third function includes the following sub-steps: Obtain the ecological flow corresponding to each serial number in the first function, and mark it as the first function predicted flow; calculate the difference between the first function predicted flow and the corresponding historical month's ecological flow; Obtain the maximum value of the difference between the flow predicted by the first function and the ecological flow of the corresponding historical month, and mark it as the first difference threshold; Obtain the ecological flow corresponding to each serial number in the second function and mark it as the second function predicted flow; Calculate the difference between the flow predicted by the second function and the ecological flow of the corresponding historical month; Obtain the maximum value of the difference between the flow predicted by the second function and the ecological flow of the corresponding historical month, and mark it as the second difference threshold; Get the ecological flow corresponding to each serial number in the third function and mark it as the predicted flow of the third function; Calculate the difference between the flow predicted by the third function and the ecological flow of the corresponding historical month; Obtain the maximum value of the difference between the flow predicted by the third function and the ecological flow of the corresponding historical month, and mark it as the third difference threshold; The first function is moved in the negative direction of the vertical axis by the distance of the first difference threshold to obtain a first threshold function; the second function is moved in the negative direction of the vertical axis by the distance of the second difference threshold to obtain a second threshold function; the third function is moved in the negative direction of the vertical axis by the distance of the third difference threshold to obtain a third threshold function; Obtaining the ecological flow calculation results of the corresponding month based on the first threshold function, the second threshold function, and the third threshold function includes the following steps: The historical ecological flow corresponding to each serial number in the first threshold function, the second threshold function, and the third threshold function is obtained and output as the ecological flow accounting result of the corresponding month.

2. The watershed ecological flow accounting method according to claim 1 is characterized in that: Screening the ecological flows of the first number of historical months to obtain the screened ecological flows includes the following sub-steps: The historical ecological flows of the same month are sorted from small to large and marked as Ss1 to Ss i ; Determine whether 0.5*Xt is an integer. If so, set Ss （0.5*Xt） As the middle position value; if not, get the integers on the left and right sides of 0.5*Xt, mark them as Zz and Zy respectively, and calculate Ss （Zz） With Ss （Zy） The mean of is taken as the middle position value; where Xt is the number of ecological flows in the same historical month; Determine whether a*Xt is an integer. If so, set Ss （a*Xt） As the first ratio value; if not, get the integers on the left and right sides of a*Xt, mark them as Dz and Dy respectively, and calculate Ss （Dz） With Ss （Dy） The mean of is taken as the first proportional value; where a is the ecological flow coefficient, and the range of a is (0.5, 1).

3. The watershed ecological flow accounting method according to claim 2 is characterized in that: Screening the ecological flows of the first number of historical months to obtain the screened ecological flows also includes the following sub-steps: The screening ecological threshold is obtained as: ; Where Cs is the screening ecological threshold, Bz is the middle position value, and B1 is the first proportion value; The screening ecological flow is obtained by deleting the historical ecological flow that is greater than the screening ecological threshold.

4. The method for calculating watershed ecological flow according to claim 3, characterized in that: Obtaining a partitioning scatter plot based on filtered ecological flow includes the following sub-steps: Get the filtered ecological flow for each historical month; Obtain the mean of the filtered ecological flow for each historical month and mark it as the historical ecological flow mean; Sort the historical ecological flow means from small to large, and mark the results of the sorting of the historical ecological flow means from small to large as corresponding serial numbers, where the serial numbers are positive integers from 1 to n; With the serial number as the horizontal axis data and the historical ecological flow mean as the vertical axis data, a coordinate system is established and marked as the division coordinate system; All serial numbers and the corresponding historical ecological flow means are used as division coordinate points to draw a division coordinate system to obtain a division scatter plot.

5. The watershed ecological flow accounting method according to claim 4 is characterized in that: Dividing a year into a first ecological period, a second ecological period, and a third ecological period based on the partition scatter plot includes the following sub-steps: Connect the adjacent division coordinate points from left to right according to the division scatter plot to obtain a line segment, marked as the first connection line segment, and obtain the slope of the first connection line segment, marked as F m , where m is a positive integer, and m represents the order of connecting line segments from small to large; Calculate adjacent F m The difference is: FC m =F m -F m-1 :FC m For adjacent F m The difference between , at this time m is greater than or equal to 2; FC m Sort from largest to smallest and get the first two FCs m Corresponding F m , marked as the monthly division threshold; the monthly division threshold is marked as the first monthly division threshold and the second monthly division threshold according to the smaller and larger values; In the partitioned scatter plot, the range of serial numbers that are less than the first partition month threshold is marked as the first ecological period; the range of serial numbers that are greater than or equal to the first partition month threshold and less than the second partition month threshold is marked as the second ecological period; the range of serial numbers that are greater than or equal to the second partition month threshold is marked as the third ecological period.

6. The method for calculating watershed ecological flow according to claim 5, characterized in that: Obtaining the functional relationship between the ecological flow and the month during the first ecological period, the second ecological period, and the third ecological period, respectively labeled as the first function, the second function, and the third function, includes the following sub-steps: During the first ecological period, the divided coordinate points are fitted with functions to obtain a first function; during the second ecological period, the divided coordinate points are fitted with functions to obtain a second function; during the third ecological period, the divided coordinate points are fitted with functions to obtain a third function.

7. A watershed ecological flow accounting system, used to implement the watershed ecological flow accounting method according to any one of claims 1 to 6, characterized in that: It includes data screening module, scatter plot drawing module, ecological period division module, function fitting module, threshold acquisition module and ecological flow accounting module; The data screening module is used to screen the ecological flows of a first number of historical months to obtain screened ecological flows; The scatter plot drawing module is used to obtain a partition scatter plot based on the screening ecological flow; The ecological period division module is used to divide a year into a first ecological period, a second ecological period and a third ecological period based on the division scatter plot; The function fitting module is used to obtain the functional relationship between the ecological flow and the month during the first ecological period, the second ecological period, and the third ecological period, respectively marked as the first function, the second function, and the third function; The threshold acquisition module is used to acquire a first threshold function, a second threshold function and a third threshold function based on the ecological flow of a second number of historical months, the first function, the second function and the third function; The ecological flow accounting module is used to obtain the ecological flow accounting result of the corresponding month based on the first threshold function, the second threshold function and the third threshold function.

Images