Basin ecological flow accounting method and system

Abstract

The invention discloses a watershed ecological flow accounting method and system, and relates to the technical field of ecological flow accounting, and the method comprises the following steps: obtaining a decentralized point diagram based on the screened ecological flow; dividing a year into a first ecological period, a second ecological period and a third ecological period based on the scatter plot; function relations between ecological flow and months among the first ecological period, the second ecological period and the third ecological period are obtained and marked 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; obtaining an ecological flow accounting result of the corresponding month based on the first threshold function, the second threshold function and the third threshold function; the method is used for solving the problem that in an existing ecological flow accounting technology, setting of the minimum ecological flow is subjective, and consequently the ecological flow accounting result is not accurate.

Description

Technical Field

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

[0002] With the rapid development of society and economy, the intensity of human development and utilization of water resources has continued to increase, and the water resource distribution pattern within the basin has undergone significant changes. As an important part of the ecosystem, the maintenance of the ecological function of rivers depends on certain water flow conditions, and ecological flow is a key factor in ensuring the health of river ecosystems. With the continuous deepening of research in fields such as ecology, hydrology and water resources management, a variety of ecological flow accounting methods have emerged, including hydrological methods, hydraulic methods, habitat simulation methods, and comprehensive methods. In the flow threshold curve method of the hydrological method, the flow at a frequency of 90% is used as the minimum ecological flow in the river. The setting of 90% is relatively subjective and cannot generate a suitable threshold based on historical ecological flow. For example, a patent application with application publication number CN113887965A discloses a watershed ecological flow accounting method. This scheme uses a conventional flow threshold curve method, that is, the flow at 90% is used as the minimum ecological flow in the river. The setting of 90% is relatively subjective, that is, the setting of the minimum ecological flow in the existing ecological flow accounting technology is relatively subjective, resulting in inaccurate ecological flow accounting results. Summary of the Invention

[0003] The present invention aims to solve one of the technical problems in the prior art to at least a certain extent. By 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, the problem of the subjective setting of the minimum ecological flow in the existing ecological flow calculation technology, which leads to inaccurate ecological flow calculation results, can be solved.

[0004] To achieve the above objectives, this application provides a watershed ecological flow accounting method, which includes the following steps: 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; Furthermore, filtering the ecological flows of the first number of historical months to obtain the filtered 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).

[0005] Furthermore, filtering the ecological flows of the first number of historical months to obtain the filtered ecological flows further 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.

[0006] Furthermore, obtaining a partitioning scatter plot based on the 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.

[0007] Furthermore, 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.

[0008] Furthermore, 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.

[0009] Furthermore, 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 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 ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the first function predicted flows and the ecological flow of the corresponding historical months, 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 second function predicted flow and the ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the second function predicted flows and the ecological flow of the corresponding historical months, and mark it as the second difference threshold; Obtain the ecological flow corresponding to each serial number in the third function, and mark it as the third function predicted flow; calculate the difference between the third function predicted flow and the corresponding ecological flow of the historical month; obtain the maximum value of the difference between all the third function predicted flows and the corresponding ecological flow of the historical month, and mark it as the third difference threshold.

[0010] Furthermore, 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 includes the following sub-steps: Move the first function in the negative direction of the vertical axis by the distance of the first difference threshold to obtain the first threshold function; move the second function in the negative direction of the vertical axis by the distance of the second difference threshold to obtain the second threshold function; move the third function in the negative direction of the vertical axis by the distance of the third difference threshold to obtain the third threshold function.

[0011] Furthermore, obtaining the ecological flow accounting 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 as the ecological flow accounting result output of the corresponding month.

[0012] 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; 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.

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

[0014] Figure 1 It is a principle block diagram of the system of the present invention; Figure 2 A schematic diagram of a partitioned scatter plot of the present invention; Figure 3 Schematic diagram of the first month-division threshold and the second month-division threshold of the present invention; Figure 4 is a schematic diagram of the first function, the second function and the third function of the present invention; Figure 5 is a schematic diagram of the first threshold function, the second threshold function and the third threshold function of the present invention; Figure 6 Flow chart of the steps of the method of the present invention. DETAILED DESCRIPTION

[0015] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0016] Example 1, please refer to Figure 1 As shown, this application 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; The data screening module is used to screen the ecological flow of the first number of historical months to obtain the filtered ecological flow; select the average ecological flow of each month in the past ten years, and select the data of the past ten years, that is, the first number is 120, and the ecological flow of the historical months is the average ecological flow within the month; The data screening module is configured with a first ratio value and an intermediate position value acquisition strategy, which includes: 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); 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; 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.

[0017] The data filtering module is configured with data filtering strategies, which include: 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.

[0018] 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. In practical application, taking sequence number 6 as an example, the screening ecological threshold is: Cs=239m 3 / s.

[0019] The scatter plot drawing module is used to obtain the partition scatter plot based on the screening ecological flow; The scatter plot drawing module is configured with a scatter plot drawing strategy, which includes: 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 mean are used as division coordinate points to draw a division coordinate system to obtain a division scatter plot; In practical applications, please refer to Figure 2 As shown, draw a partition scatter plot; 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 ecological period division module is configured with ecological interval division strategies, which include: 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 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. In practical applications, please refer to Figure 3 As shown, the first and second division month thresholds obtained are 4 and 9 respectively. Therefore, the range of serial numbers less than 4 is marked as the first ecological period; the range of serial numbers greater than or equal to 4 and less than 9 is marked as the second ecological period; the range of serial numbers greater than or equal to 9 is marked as the third ecological period. The segmented analysis of the ecological flow accounting results is more accurate.

[0020] The function fitting module is used to obtain the functional relationship between ecological flow and month during the first ecological period, the second ecological period, and the third ecological period, which are marked as the first function, the second function, and the third function respectively; In practical applications, please refer to Figure 4 As shown, the first function, the second function and the third function are obtained; The function fitting module is configured with function fitting strategies, which include: During the first ecological period, the divided coordinate points are fitted with a function to obtain a first function; during the second ecological period, the divided coordinate points are fitted with a function to obtain a second function; during the third ecological period, the divided coordinate points are fitted with a function to obtain a 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 threshold acquisition module is configured with a threshold acquisition strategy, which includes: 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 ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the first function predicted flows and the ecological flow of the corresponding historical months, 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 second function predicted flow and the ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the second function predicted flows and the ecological flow of the corresponding historical months, and mark it as the second difference threshold; Obtain the ecological flow corresponding to each serial number in the third function, and mark it as the third function predicted flow; calculate the difference between the third function predicted flow and the ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the third function predicted flows and the ecological flow of the corresponding historical months, and mark it as the third difference threshold; In practical applications, please refer to Figure 5 As shown, taking the second function as an example, the second function predicts the flow rate to be 239m at sequence number 6. 3 / s, taking the ecological flow of a historical month as 230m 3 / s is an example, the calculated difference is 9m 3 / s, and then obtain the difference between the second function predicted flow and the corresponding historical month's ecological flow, the maximum value is 9m 3 / s, the second difference threshold is 9m 3 / s, similarly, get the second difference threshold 12m 3 / s, and the third difference threshold is 16m 3 / s.

[0021] The threshold acquisition module is configured with a threshold function acquisition strategy, which includes: Move the first function in the negative direction of the vertical axis by the distance of the first difference threshold to obtain the first threshold function; move the second function in the negative direction of the vertical axis by the distance of the second difference threshold to obtain the second threshold function; move the third function in the negative direction of the vertical axis by the distance of the third difference threshold to obtain the third threshold function.

[0022] The ecological flow accounting module is used to obtain the ecological flow accounting results of the corresponding month based on the first threshold function, the second threshold function and the third threshold function; The ecological flow accounting module is equipped with an ecological flow accounting strategy, which includes: Obtain the historical ecological flow corresponding to each serial number in the first threshold function, the second threshold function, and the third threshold function as the ecological flow accounting result of the corresponding month; In practical applications, please refer to Figure 5 As shown, the first threshold function, the second threshold function, and the third threshold function corresponding to each serial number are obtained as the ecological flow accounting result of the corresponding month. For example, when the serial number is 1, the first threshold function is obtained as 98m 3 / s, the month of the serial number corresponds to December, so 98m 3 / s, the ecological flow accounting result of the watershed in December is 98m 3 / s, and similarly obtain the basin ecological flow accounting results for other months.

[0023] Example 2, please refer to Figure 6 As shown, this application provides a watershed ecological flow accounting method, which includes the following steps: Step S1, filtering the ecological flows of a first number of historical months to obtain filtered ecological flows; Step S1 includes the following sub-steps: Step S101: sort the historical ecological flows of the same month from small to large and mark them as Ss1 to Ss2. i ; Step S102, 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; Step S103, 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); Step S104, obtaining the screening ecological threshold: ; Where Cs is the screening ecological threshold, Bz is the middle position value, and B1 is the first proportion value; Step S105 , deleting the historical ecological flow that is greater than the screening ecological threshold value to obtain the screening ecological flow.

[0024] Step S2, obtaining a partition scatter plot based on the filtered ecological flow; Step S2 includes the following sub-steps: Step S201, obtaining the filtered ecological flow of each historical month; Step S202, obtaining the mean of the filtered ecological flow for each historical month, and marking it as the historical ecological flow mean; Step S203: sort the historical ecological flow mean values ​​from small to large, and mark the sorted historical ecological flow mean values ​​from small to large as corresponding sequence numbers, where the sequence numbers are positive integers from 1 to n; Step S204: establish a coordinate system with the serial number as the horizontal axis data and the historical ecological flow mean as the vertical axis data, and mark it as the division coordinate system; Step S205: All serial numbers and corresponding historical ecological flow means are used as division coordinate points to draw a division coordinate system to obtain a division scatter plot.

[0025] Step S3, dividing the year into a first ecological period, a second ecological period, and a third ecological period based on the partition scatter plot; Step S3 includes the following sub-steps: Step S301: Connect 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; Step S302, calculate the 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; Step S303: 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; Step S304: In the divided scatter plot, the range of serial numbers that are less than the first division month threshold is marked as the first ecological period; the range of serial numbers that are 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; and the range of serial numbers that are greater than or equal to the second division month threshold is marked as the third ecological period.

[0026] Step S4, respectively 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, and marking them as the first function, the second function, and the third function, respectively; Step S4 includes the following sub-steps: Step S401, performing function fitting on the divided coordinate points in the first ecological period to obtain a first function; performing function fitting on the divided coordinate points in the second ecological period to obtain a second function; performing function fitting on the divided coordinate points in the third ecological period to obtain a third function.

[0027] Step S5, 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; Step S5 includes the following sub-steps: Step S501: Obtain the ecological flow corresponding to each serial number in the first function, marking it as the first function predicted flow; calculate the difference between the first function predicted flow and the ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the first function predicted flows and the ecological flow of the corresponding historical months, marking it as the first difference threshold; Step S502: Obtain the ecological flow corresponding to each serial number in the second function, marking it as the second function predicted flow; calculate the difference between the second function predicted flow and the ecological flow of the corresponding historical month; obtain the maximum value of the difference between all the second function predicted flows and the ecological flow of the corresponding historical months, marking it as the second difference threshold; Step S503, obtain the ecological flow corresponding to each serial number in the third function, mark it as the third function predicted flow; calculate the difference between the third function predicted flow and the corresponding historical month's ecological flow; obtain the maximum value of the difference between all third function predicted flows and the corresponding historical months' ecological flow, mark it as the third difference threshold.

[0028] Step S504, move the first function in the negative direction of the vertical axis by the distance of the first difference threshold to obtain the first threshold function; move the second function in the negative direction of the vertical axis by the distance of the second difference threshold to obtain the second threshold function; move the third function in the negative direction of the vertical axis by the distance of the third difference threshold to obtain the third threshold function.

[0029] Step S6, 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; Step S6 includes the following sub-steps: Step S601: Obtain the historical ecological flow corresponding to each serial number in the first threshold function, the second threshold function, and the third threshold function as the ecological flow accounting result of the corresponding month.

[0030] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media containing computer-usable program code. The storage medium may be implemented by any type of volatile or non-volatile storage device, or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0031] In the embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interface, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or 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; The ecological flow accounting results of the corresponding months are obtained based on the first threshold function, the second threshold function and the third threshold function.

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. The method for calculating watershed ecological flow according to claim 6, characterized in that: 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.

8. The method for calculating watershed ecological flow according to claim 7, characterized in that: 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 further includes the following sub-steps: Move the first function in the negative direction of the vertical axis by the distance of the first difference threshold to obtain the first threshold function; move the second function in the negative direction of the vertical axis by the distance of the second difference threshold to obtain the second threshold function; move the third function in the negative direction of the vertical axis by the distance of the third difference threshold to obtain the third threshold function.

9. The watershed ecological flow accounting method according to claim 8 is characterized in that: 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 as the ecological flow accounting result output of the corresponding month.

10. A watershed ecological flow accounting system, used to implement the watershed ecological flow accounting method according to any one of claims 1 to 9, 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