CN113627095B - Ecological flow process construction method and system based on daily flow process change trend - Google Patents
Ecological flow process construction method and system based on daily flow process change trend Download PDFInfo
- Publication number
- CN113627095B CN113627095B CN202110787674.0A CN202110787674A CN113627095B CN 113627095 B CN113627095 B CN 113627095B CN 202110787674 A CN202110787674 A CN 202110787674A CN 113627095 B CN113627095 B CN 113627095B
- Authority
- CN
- China
- Prior art keywords
- daily
- flow process
- trend
- flow
- ecological
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 163
- 230000008569 process Effects 0.000 title claims abstract description 139
- 230000008859 change Effects 0.000 title claims abstract description 72
- 238000010276 construction Methods 0.000 title claims abstract description 29
- 241000894007 species Species 0.000 claims abstract description 24
- 238000011160 research Methods 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims abstract description 9
- 238000011835 investigation Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 31
- 230000009467 reduction Effects 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 230000004044 response Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides an ecological flow process construction method and system based on daily flow process variation trend, wherein the method comprises the following steps: step 1, investigation is conducted on the current situation of biological resources in a research area, and indication species and the lowest proper ecological flow of the research area are determined; step 2, selecting long-series daily flow data with small influence on human activities in early stage, and calculating a mean daily flow process for years; step 3, calculating the slope and deflection angle of each point according to the average daily flow process for many years; step 4, determining the upper and lower trend change points and trend change lines in the process of the average daily flow for years according to the slope; linear interpolation is carried out on the upper and lower trend line points to obtain upper and lower trend lines, and average values of the upper and lower trend lines are taken from overlapping sections of the upper and lower trend lines to obtain average daily flow process variation trend lines for a plurality of years; step 5, calculating the fluctuation rate of the variation trend according to the variation trend line of the average daily flow process for many years; and step 6, obtaining a daily ecological flow process line.
Description
Technical Field
The invention belongs to the technical field of ecological flow construction, and particularly relates to an ecological flow process construction method and system based on a daily flow process change trend.
Technical Field
Ecological water demand refers to the necessary flow (water volume, water level) to be maintained in a river or in a lake and its process, also called ecological flow in a river, which maintains the structure and function of the water ecosystem. With regard to ecological water demand, domestic research is relatively late, and a plurality of research methods are developed on the basis of introducing a large amount of foreign theories and methods. In recent years, the ecological water demand is being studied deeply, and the water ecological system is being emphasized more and more, but the meaning awareness of the ecological water demand is still different, the calculation method is also various, and the calculation of the ecological water demand is mainly divided into four major categories, namely a hydrology method, a hydraulics method, a habitat simulation method and an overall method. At present, the unified view is that ecological water needs to consider the requirement of a protection target, space-time variation and dynamic variation of water quantity.
In the ecological water demand calculation method, the data of the hydrology method is easy to obtain, the calculation is simple and quick, and the Tennant method and the improvement method thereof are most widely applied. The hydrodynamics method mainly considers the river section condition, wherein the R2CROSS method based on the Manning formula is more used. Habitat simulation considers living requirements of indicated species, and the result is most fit with river ecosystem requirements, but is complex in calculation. The whole method mainly comprises BBM method of south Africa and Australian whole evaluation method. The calculation methods of the ecological water demand are difficult to obtain the daily-scale ecological flow process, and also difficult to quickly and accurately obtain the actual condition data of the composite ecological flow process, and the calculation result cannot meet the actual demand.
Disclosure of Invention
The invention aims to solve the problems, and aims to provide an ecological flow process construction method and system based on the daily flow process variation trend, which can quickly and accurately obtain the data of the actual situation of the composite ecological flow process.
In order to achieve the above object, the present invention adopts the following scheme:
< method >
The invention provides an ecological flow process construction method based on a daily flow process change trend, which is characterized by comprising the following steps of:
step 1, investigation is conducted on the current situation of biological resources in a research area, and indication species in the research area are determined; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 ;
Step 2, selecting long-series daily flow data with small influence on human activities in early stage, and calculating a mean daily flow process for years;
step 3, calculating the slope k of the ith and (i+1) th day flow folding line segments of each point according to the average day flow process for years i And calculates the deflection angle theta of the ith day relative to the ith-1 th day flow folding line segment i ;
Step 4, selecting theta i >Point of 80 DEG, according to slope k i And determining the change points and the change lines of the upper trend and the lower trend in the flow process of the average daily flow for many years. When k is i <When 0, the corresponding point is an upper trend line point; when k is i >The corresponding point is a lower trend line point when 0; linearly interpolating the upper and lower trend line points to obtain upper and lower trend lines, and taking the average value of upper and lower trend lines at the overlapping segment (same date range) to obtain multiple trend linesAnnual average daily flow rate process change trend line;
step 5, calculating the fluctuation rate delta q of the daily flow variation trend according to the variation trend line of the average daily flow process for years J :
Non-coincident segments:
overlapping:
in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J ;
Step 6. At Δq J >The interval of 0 is the daily flow rate increase interval (water rising), and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is a solar flow rate reduction interval (falling into water), and the interval deltaq is taken J The average value of (1) is subjected to flow drop; the minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first change interval, carrying out fluctuation change according to the rule; finally, starting from the last water falling interval (generally the water rising interval), the daily flow rate is continuously reduced until the flow rate falls to the minimum proper ecological flow rate Q e0 The method comprises the steps of carrying out a first treatment on the surface of the Thus, the daily ecological flow process line can be obtained.
Preferably, the ecological flow process construction method based on the daily flow process variation trend provided by the invention can also have the following characteristics: in step 3, the slope k is calculated using the formula i And deflection angle theta i :
Wherein Q (i) is the daily flow rate on the ith day, x (i) is the date, θ i The value range of (2) is [0, pi ]]。
Preferably, the ecological flow process construction method based on the daily flow process variation trend provided by the invention can also have the following characteristics: in step 6, from the first flood interval, the flood change is performed according to the above-mentioned rule.
Preferably, the ecological flow process construction method based on the daily flow process variation trend provided by the invention further comprises the following steps: and 7, calculating the daily scale month or year ecological flow process through the daily ecological flow process.
< System >
Furthermore, the invention also provides an ecological flow process construction system based on the daily flow process change trend, which is characterized by comprising the following steps:
a determining unit for determining the species indicated in the investigation region by investigating the current situation of the biological resource in the investigation region; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 ;
An average daily flow rate calculation unit for selecting long-series daily flow rate data less affected by human activities in the early stage and calculating an average daily flow rate process for years;
slope and deflection angle calculation unit for calculating slope k of ith and (i+1) th day flow folding line segments of each point according to years average day flow process i And calculates the deflection angle theta of the ith day relative to the ith-1 th day flow folding line segment i ;
A change trend line acquisition unit for selecting θ i >Point of 80 DEG, according to slope k i Determining how muchThe change points and the change lines of the upward trend and the downward trend of the annual average daily flow process. When k is i <When 0, the corresponding point is an upper trend line point; when k is i >The corresponding point is a lower trend line point when 0; linear interpolation is carried out on the upper and lower trend line points to obtain upper and lower trend lines, the average value of the upper and lower trend lines is taken from the overlapping section of the upper and lower trend lines, and the change trend line of the average daily flow process for a plurality of years is obtained;
a change trend fluctuation rate calculating part for calculating change trend fluctuation rate deltaq according to the change trend line of the average daily flow process for years J :
Non-coincident segments:
overlapping:
in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J ;
A daily ecological flow process line acquisition unit for acquiring a daily ecological flow process line at Δq J >The interval of 0 is the daily flow rate increase interval, and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is the daily flow rate reduction interval, and the interval deltaq is taken J The average value of (1) is subjected to flow drop; the minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first change interval, carrying out fluctuation change according to the rule; finally, starting from the last water falling interval, the daily flow is continuously reduced until the flow falls to the minimum proper ecological flow Q e0 The method comprises the steps of carrying out a first treatment on the surface of the Thus obtaining a daily ecological flow process line;
an input display part which is communicated with the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part and the daily ecological flow process line acquiring part, and allows a user to input a control instruction and to display correspondingly according to the control instruction; and
and the control part is communicated with the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part, the daily ecological flow process line acquiring part and the input display part and controls the operation of the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part, the daily ecological flow process line acquiring part and the input display part.
Preferably, the ecological flow process construction method based on the daily flow process variation trend provided by the invention can also have the following characteristics: the input display unit can display the daily ecological flow rate process line obtained by the daily ecological flow rate process line obtaining unit according to the operation instruction.
Preferably, the gate dam group low-impact emergency ecological dispatching system based on the adjustable capacity provided by the invention can also have the following characteristics: the slope and deflection angle calculation unit calculates the slope k using the following formula i And deflection angle theta i :
Wherein Q (i) is the daily flow rate on the ith day, x (i) is the date, θ i The value range of (2) is [0, pi ]]。
Preferably, the gate dam group low-impact emergency ecological dispatching system based on the adjustable capacity provided by the invention can also have the following characteristics: in the daily ecological flow process line acquisition unit, fluctuation is performed according to the above-described rule from the first fluctuation interval.
Preferably, the low-impact emergency ecological dispatching system for the gate dam group based on the adjustable capacity provided by the invention further comprises: and a daily scale flow rate calculation unit for calculating a daily scale month or year ecological flow rate process by the daily ecological flow rate process.
Preferably, the gate dam group low-impact emergency ecological dispatching system based on the adjustable capacity provided by the invention can also have the following characteristics: the input display unit also displays the month or year ecological flow course obtained by the daily scale flow calculation unit according to the operation instruction.
Preferably, the gate dam group low-impact emergency ecological dispatching system based on the adjustable capacity provided by the invention also has the following characteristics: the input display unit can indicate species and minimum proper ecological flow rate Q of the indicated species to the research area determined by the determination unit according to the corresponding operation instruction e0 The average daily flow rate process for many years calculated by the average daily flow rate calculation unit can be displayed, and the slope k calculated by the slope and deflection angle calculation unit can be displayed i And deflection angle theta i The trend line acquisition unit may be configured to acquire a trend line of the average daily flow rate of the years, and display the trend line of the average daily flow rate of the years, and calculate the trend fluctuation rate Δq of the change calculated by the trend fluctuation rate calculation unit J And displaying.
Effects and effects of the invention
According to the ecological flow process construction method and system based on the daily flow process variation trend, the defect that daily flow process calculation is lacking in an ecological flow process calculation method is overcome, the data of calculation requirements are easy to obtain, the method is simple and convenient, the realization is easy, the calculation result can meet the flow requirements of indicated species, the water quantity variation of a natural flow process is met, the constructed ecological flow process is more accurate and meets the actual condition, the rationality of ecological flow result analysis is improved, and the river channel ecological system is kept stable.
Drawings
FIG. 1 is a schematic diagram of a yearly average daily flow rate process and a constructed daily ecological flow rate process according to an embodiment of the invention;
FIG. 2 is a schematic diagram of upper and lower trend lines and a change trend line according to an embodiment of the present invention;
FIG. 3 shows the fluctuation rate Δq according to the embodiment of the present invention J Schematic diagram.
Detailed Description
The following describes in detail the specific embodiments of the ecological flow process construction method and system based on the daily flow process variation trend according to the present invention with reference to the accompanying drawings.
< example >
The ecological flow process construction method based on the daily flow process change trend provided by the embodiment comprises the following steps:
firstly, investigating the current situation of biological resources in a research area, determining the indication species of the research area, and selecting a certain fish as the indication species for the river ecological flow process; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 。
And step two, selecting long-series daily flow data with small influence on human activities in early stage, and calculating a years average daily flow process, as shown in the years average daily flow process in figure 1.
Step three, calculating the slope k of each point according to the average daily flow process for many years i And deflection angle theta i :
Wherein Q (i) is the daily flow rate on the ith day, x (i) is the date, k i Slope, θ, of the flow trace for day i and day i+1 i Is the corner between two adjacent broken line sections, theta i The value range of (2) is [0, pi ]]。
Step four, selecting theta i >Point of 80 DEG, according to slope k i And determining the change points and the change lines of the upper trend and the lower trend in the flow process of the average daily flow for many years. When k is i <When 0, the corresponding point is an upper trend line point; when k is i >And the point corresponding to 0 is a lower trend line point. And (3) linearly interpolating the upper and lower trend line points to obtain upper and lower trend lines, and taking the average value of the upper and lower trend lines at the overlapping section (same date range) of the upper and lower trend lines to obtain the average daily flow process variation trend line for a plurality of years, as shown in figure 2.
Step five, calculating the fluctuation rate delta q of the change trend according to the change trend line of the average daily flow process for years J :
Non-coincident segments:
overlapping:in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J As shown in fig. 3.
Step six, at Δq J >The interval of 0 is the daily flow rate increase interval (water rising), and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is a solar flow rate reduction interval (falling into water), and the interval deltaq is taken J Is subjected to flow dropping.
The minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first variation section (generally, the fluctuation section), the fluctuation is performed according to the rule; finally, from the last water falling interval, the daily flow is continuously reduced,until the flow rate falls to the minimum proper ecological flow rate Q e0 . Thus, a daily ecological flow process line, such as the daily ecological flow process in FIG. 1, can be obtained.
And step seven, calculating the lunar ecological flow process through the daily ecological flow process.
Further, the embodiment also provides an ecological flow process construction system capable of automatically realizing the method, which comprises a determination part, an average daily flow calculation part, a slope and deflection angle calculation part, a change trend line acquisition part, a change trend fluctuation rate calculation part, a daily ecological flow process line acquisition part, an input display part and a control part.
The determining part is used for investigating the current situation of biological resources in the investigation region and determining the indication species of the investigation region; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 。
The average daily flow rate calculation unit selects long-series daily flow rate data which is less affected by human activities in the early stage, and calculates an average daily flow rate course for many years.
The slope and deflection angle calculation part calculates the slope k of the ith and (i+1) th day flow folding line segments of each point according to the average day flow process for years i And calculates the deflection angle theta of the ith day relative to the ith-1 th day flow folding line segment i 。
A change trend line acquisition unit: select θ i >Point of 80 DEG, according to slope k i And determining the change points and the change lines of the upper trend and the lower trend in the flow process of the average daily flow for many years. When k is i <When 0, the corresponding point is an upper trend line point; when k is i >The corresponding point is a lower trend line point when 0; and linearly interpolating the upper and lower trend line points to obtain upper and lower trend lines, and taking the average value of the upper and lower trend lines at the overlapping section of the upper and lower trend lines to obtain the variation trend line of the average daily flow process for a plurality of years.
The change trend fluctuation rate calculating part calculates the change trend fluctuation rate deltaq according to the change trend line of the average daily flow process for years J :
Non-coincident segments:
overlapping:
in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J ;
Daily ecological flow process line acquisition part: at Δq J >The interval of 0 is the daily flow rate increase interval, and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is the daily flow rate reduction interval, and the interval deltaq is taken J The average value of (1) is subjected to flow drop; the minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first change interval, carrying out fluctuation change according to the rule; finally, starting from the last water falling interval, the daily flow is continuously reduced until the flow falls to the minimum proper ecological flow Q e0 The method comprises the steps of carrying out a first treatment on the surface of the Thus, a daily ecological flow process line is obtained.
The input display part is communicated with the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part and the daily ecological flow process line acquiring part, so that a user inputs a control command and carries out corresponding display according to the control command.
The control part is communicated with the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part, the daily ecological flow process line acquiring part and the input display part, and controls the operation of the determining part, the average daily flow rate calculating part, the slope and deflection angle calculating part, the change trend line acquiring part, the change trend fluctuation rate calculating part, the daily ecological flow process line acquiring part and the input display part.
The above embodiments are merely illustrative of the technical solutions of the present invention. The ecological flow process construction method and system based on the daily flow process variation trend according to the present invention are not limited to the description of the above embodiments, but the scope of the invention is defined by the claims. Any modifications, additions or equivalent substitutions made by those skilled in the art based on this embodiment are within the scope of the invention as claimed in the claims.
Claims (10)
1. The ecological flow process construction method based on the daily flow process change trend is characterized by comprising the following steps of:
step 1, investigation is conducted on the current situation of biological resources in a research area, and indication species in the research area are determined; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 ;
Step 2, selecting long-series daily flow data with small influence on human activities in early stage, and calculating a mean daily flow process for years;
step 3, calculating the slope k of the ith and (i+1) th day flow folding line segments of each point according to the average day flow process for years i And calculates the deflection angle theta of the ith day relative to the ith-1 th day flow folding line segment i ;
Step 4, selecting theta i >Point of 80 DEG, according to slope k i Determining the upper and lower trend change points and trend change lines in the average daily flow process for many years; when k is i <When 0, the corresponding point is an upper trend line point; when k is i >The corresponding point is a lower trend line point when 0; linear interpolation is carried out on the upper and lower trend line points to obtain upper and lower trend lines, and average values of the upper and lower trend lines are taken from overlapping sections of the upper and lower trend lines to obtain average daily flow process variation trend lines for a plurality of years;
step 5, calculating the fluctuation rate delta q of the variation trend according to the variation trend line of the average daily flow process for many years J :
Non-coincident segments:
overlapping:
in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J ;
Step 6. At Δq J >The interval of 0 is the daily flow rate increase interval, and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is the daily flow rate reduction interval, and the interval deltaq is taken J The average value of (1) is subjected to flow drop; the minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first change interval, carrying out fluctuation change according to the rule; finally, starting from the last water falling interval, the daily flow is continuously reduced until the flow falls to the minimum proper ecological flow Q e0 The method comprises the steps of carrying out a first treatment on the surface of the Thus, the daily ecological flow process line can be obtained.
2. The ecological flow process construction method based on the daily flow process variation trend according to claim 1, wherein the method comprises the following steps:
wherein in step 3, the slope k is calculated using the formula i And deflection angle theta i :
Wherein Q (i) is the daily flow rate on the ith day, x (i) is the date, θ i The value range of (2) is [0, pi ]]。
3. The ecological flow process construction method based on the daily flow process variation trend according to claim 1, wherein the method comprises the following steps:
in step 6, the fluctuation is changed according to the rule from the first fluctuation interval.
4. The ecological flow process construction method based on the daily flow process variation trend according to claim 1, further comprising:
and 7, calculating the daily scale month or year ecological flow process through the daily ecological flow process.
5. Ecological flow process construction system based on daily flow process variation trend, which is characterized by comprising:
a determining unit for determining the species indicated in the investigation region by investigating the current situation of the biological resource in the investigation region; determining the minimum proper ecological flow Q of the indicated species according to the river cross section shape and the proper flow velocity and water depth curve of the indicated species e0 ;
An average daily flow rate calculation unit for selecting long-series daily flow rate data less affected by human activities in the early stage and calculating an average daily flow rate process for years;
slope and deflection angle calculation unit for calculating slope k of ith and (i+1) th day flow folding line segments of each point according to years average day flow process i And calculates the deflection angle theta of the ith day relative to the ith-1 th day flow folding line segment i ;
A change trend line acquisition unit for selecting θ i >Point of 80 DEG, according to slope k i Determining the upper and lower trend change points and trend change lines in the average daily flow process for many years; when (when)k i <When 0, the corresponding point is an upper trend line point; when k is i >The corresponding point is a lower trend line point when 0; linear interpolation is carried out on the upper and lower trend line points to obtain upper and lower trend lines, the average value of the upper and lower trend lines is taken from the overlapping section of the upper and lower trend lines, and the change trend line of the average daily flow process for a plurality of years is obtained;
a change trend fluctuation rate calculating part for calculating change trend fluctuation rate deltaq according to the change trend line of the average daily flow process for years J :
Non-coincident segments:
overlapping:
in which Q j ' is the j day flow of the upper/lower trend line of the non-coincident section, Q j "day-of-the-day flow of j for transition trend line of overlapping section, j 0 Days, j, corresponding to the first point starting with the upper/lower trend line 1 The number of days, j, corresponding to the first point of the upper and lower trend line overlapping section 2 The last point of the upper and lower trend line overlapping section is ended for the corresponding days; j (j) 0 -j 1 The segment is an upper and lower trend line non-coincident segment, and may appear before, on both sides of or after the coincident segment; the fluctuation and variation calculated by the overlapping section and the non-overlapping section jointly form the fluctuation and variation rate delta q J ;
A daily ecological flow process line acquisition unit for acquiring a daily ecological flow process line at Δq J >The interval of 0 is the daily flow rate increase interval, and the interval deltaq is taken J The average value of (2) is subjected to flow increase; Δq J <The interval of 0 is the daily flow rate reduction interval, and the interval deltaq is taken J The average value of (1) is subjected to flow drop; the minimum proper ecological flow Q is kept at the beginning of year e0 Starting from the first change interval, carrying out fluctuation change according to the rule; finally, starting from the last water falling interval, the daily flow is continuously reduced until the flow falls to the minimum proper ecological flow Q e0 The method comprises the steps of carrying out a first treatment on the surface of the Thereby obtaining the daily ecological flow process line;
An input display unit which is communicatively connected to the determination unit, the average daily flow rate calculation unit, the slope and deflection angle calculation unit, the change trend line acquisition unit, the change trend fluctuation rate calculation unit, and the daily ecological flow rate process line acquisition unit, and which allows a user to input a control instruction and to display the control instruction in response to the control instruction; and
and a control unit which is communicatively connected to the determination unit, the average daily flow rate calculation unit, the slope and deflection angle calculation unit, the change trend line acquisition unit, the change trend fluctuation rate calculation unit, the daily ecological flow rate process line acquisition unit, and the input display unit, and controls the operation of the determination unit, the average daily flow rate calculation unit, the slope and deflection angle calculation unit, the change trend line acquisition unit, the change trend fluctuation rate calculation unit, the daily ecological flow rate process line acquisition unit, and the input display unit.
6. The ecological flow process construction system based on the daily flow process variation trend according to claim 5, wherein:
the input display unit can display the daily ecological flow rate process line obtained by the daily ecological flow rate process line obtaining unit according to an operation instruction.
7. The ecological flow process construction method based on the daily flow process variation trend according to claim 5, wherein the ecological flow process construction method is characterized by:
wherein the slope and deflection angle calculation section calculates the slope k using the following formula i And deflection angle theta i :
Wherein Q (i) is the daily flow rate on the ith day, x (i) is the date, θ i The value range of (2) is [0, pi ]]。
8. The ecological flow process construction method based on the daily flow process variation trend according to claim 5, wherein the ecological flow process construction method is characterized by:
wherein, in the daily ecological flow process line acquisition part, the fluctuation change is performed according to the rule from the first fluctuation interval.
9. The ecological flow process construction method based on the daily flow process variation trend according to claim 5, further comprising:
and a daily scale flow rate calculation unit for calculating a daily scale month or year ecological flow rate process by the daily ecological flow rate process.
10. The ecological flow process construction system based on the daily flow process variation trend according to claim 9, wherein:
the input display part is used for displaying the month or year ecological flow course obtained by the daily scale flow calculation part according to an operation instruction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110787674.0A CN113627095B (en) | 2021-07-13 | 2021-07-13 | Ecological flow process construction method and system based on daily flow process change trend |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110787674.0A CN113627095B (en) | 2021-07-13 | 2021-07-13 | Ecological flow process construction method and system based on daily flow process change trend |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113627095A CN113627095A (en) | 2021-11-09 |
CN113627095B true CN113627095B (en) | 2023-10-24 |
Family
ID=78379594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110787674.0A Active CN113627095B (en) | 2021-07-13 | 2021-07-13 | Ecological flow process construction method and system based on daily flow process change trend |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113627095B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012017628A (en) * | 2010-07-09 | 2012-01-26 | Asia Air Survey Co Ltd | Flood flow estimation system and method |
CN107103188A (en) * | 2017-04-11 | 2017-08-29 | 山东大学 | A kind of determination method of urban river ecological flow threshold |
CN108168635A (en) * | 2018-01-31 | 2018-06-15 | 长江水利委员会水文局 | Discharge of river realtime on-line monitoring method and device |
CN111144621A (en) * | 2019-12-06 | 2020-05-12 | 河海大学 | Reservoir group flood control and power generation competition relation judging method based on trend judgment |
CN113065280A (en) * | 2021-03-17 | 2021-07-02 | 南京工业职业技术大学 | Method for predicting integral change trend of runoff in drainage basin |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9625497B2 (en) * | 2013-04-26 | 2017-04-18 | Telefonaktiebolaget Lm Ericsson (Publ) | Predicting a network performance measurement from historic and recent data |
-
2021
- 2021-07-13 CN CN202110787674.0A patent/CN113627095B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012017628A (en) * | 2010-07-09 | 2012-01-26 | Asia Air Survey Co Ltd | Flood flow estimation system and method |
CN107103188A (en) * | 2017-04-11 | 2017-08-29 | 山东大学 | A kind of determination method of urban river ecological flow threshold |
CN108168635A (en) * | 2018-01-31 | 2018-06-15 | 长江水利委员会水文局 | Discharge of river realtime on-line monitoring method and device |
CN111144621A (en) * | 2019-12-06 | 2020-05-12 | 河海大学 | Reservoir group flood control and power generation competition relation judging method based on trend judgment |
CN113065280A (en) * | 2021-03-17 | 2021-07-02 | 南京工业职业技术大学 | Method for predicting integral change trend of runoff in drainage basin |
Also Published As
Publication number | Publication date |
---|---|
CN113627095A (en) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Contribution of multiple climatic variables and human activities to streamflow changes across China | |
Aberle et al. | Flow resistance of emergent rigid and flexible floodplain vegetation | |
CN103940974B (en) | Based on the middling scale watershed soil erosion Dynamic Spatial And Temporal Analysis method of GIS | |
Gao et al. | Shallow groundwater plays an important role in enhancing irrigation water productivity in an arid area: The perspective from a regional agricultural hydrology simulation | |
CN107401140B (en) | The determination method of tributary lower reaches section channel forming disahcge | |
CN101944161A (en) | Calculation method of ecological water demand of wetland based on evaluation of disturbance degree of ecological system | |
CN110400014A (en) | A kind of coastal cities multi-source flood method for numerical simulation based on GIS grid operation | |
CN111802287A (en) | Spawning ground positioning inversion simulation method based on development time of drifting fish egg embryo body | |
CN110276150A (en) | A kind of Mountain Area river basal flow capacity system interpolation extension method based on Copula function | |
CN113627095B (en) | Ecological flow process construction method and system based on daily flow process change trend | |
Wu et al. | Modeling of mixed crop field water demand and a smart irrigation system | |
Xu et al. | Analysis of water balance in Poyang Lake basin and subsequent response to climate change | |
CN107621532A (en) | The method that the outstanding critical shear stress of shallow lake bed mud is determined based on mutation analysis | |
CN109145396B (en) | River course roughness zoning calibration method based on vegetation distribution | |
He et al. | A framework for attributing runoff changes based on a monthly water balance model: An assessment across China | |
He et al. | Diurnal variation characteristics of thermal structure in a deep reservoir and the effects of selective withdrawal | |
Paudel | Role of sediment in the design and management of irrigation canals: Sunsari Morang Irrigation Scheme, Nepal | |
CN105760665A (en) | Calculation method for ecological water requirement of affected plain area river network | |
Han et al. | Spatial-temporal coupling pattern between irrigation demand and soil moisture dynamics throughout wheat-maize rotation system in the North China Plain | |
CN109902266A (en) | A kind of streamflow routing method based on Copula function | |
CN108305244A (en) | A kind of division methods and system of the soft or hard region of variation of crop | |
CN112561205B (en) | Quantitative measurement method considering influence of silt dam and terrace on runoff | |
CN208523172U (en) | A kind of hidden pipe adds vertical shaft type salt-soda soil desalting system | |
Zhang et al. | The impact of plant morphology on flow structure: comparative analysis of two types of submerged flexible macrophyte | |
CN205373998U (en) | Sheet flow roll wave measuring device based on pressure sensor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |