CN117541031B - Full-period proper ecological flow interval deducing method - Google Patents
Full-period proper ecological flow interval deducing method Download PDFInfo
- Publication number
- CN117541031B CN117541031B CN202410030698.5A CN202410030698A CN117541031B CN 117541031 B CN117541031 B CN 117541031B CN 202410030698 A CN202410030698 A CN 202410030698A CN 117541031 B CN117541031 B CN 117541031B
- Authority
- CN
- China
- Prior art keywords
- period
- flow
- ecological
- spawning
- control section
- 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 45
- 241000251468 Actinopterygii Species 0.000 claims abstract description 75
- 230000001186 cumulative effect Effects 0.000 claims abstract description 10
- 230000009467 reduction Effects 0.000 claims abstract description 10
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 230000004044 response Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 25
- 230000005526 G1 to G0 transition Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000017448 oviposition Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 230000007547 defect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Operations Research (AREA)
- Mining & Mineral Resources (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Sewage (AREA)
Abstract
The invention provides a method for calculating a full-period proper ecological flow interval; based on the historical natural daily runoff data of the river basin ecological control section, a cumulative flow proportion difference method is provided, a year-by-year difference sequence is obtained, the runoff period is divided into three stages by combining the main fish spawning and propagation period data of the river basin ecological control section, and for the annual stationary period, the flow corresponding to the historical daily runoff p and 1-p quantiles is calculated as the upper and lower limit ranges of the proper ecological flow of the period; for the spawning sensitive period of fish, calculating the upper and lower limit ranges of the proper ecological flow in the period by adopting a biological hydrological response model; for the intra-year fluctuation period, a suitable ecological flow interval range calculating method based on a scene recourse reduction technology is provided, and the upper and lower limit ranges of the suitable ecological flow in the period are calculated; and finally, summarizing and splicing the sections with the proper ecological flow rate in each stage to form a full-period range with the proper ecological flow rate.
Description
Technical Field
The invention relates to the field of hydrology ecological flow calculation, in particular to a full-period suitable ecological flow interval calculating method.
Background
Along with the construction and the operation of the water engineering, the development and the utilization degree of human beings on river water resources are continuously improved, and the health of the original river ecological system is affected to a certain extent. In order to maintain the structure and the function of the ecological system of the river and the lake and improve the quality and the stability of the ecological system of the river and the lake, the concept of ecological flow is provided and is widely applied to the fields of reservoir dispatching operation, river basin planning management and the like.
According to incomplete statistics, the current ecological flow calculation method is more than 200, and can be divided into four types of hydrology methods, hydraulics methods, habitat simulation methods and integral methods. The hydrologic method only needs flow data, so that the calculation is simple and convenient, and the application is the most extensive. However, the existing hydrology method has two defects: firstly, combining the hydrologic method with the historical runoff data, adopting a statistical method to calculate the ecological flow process, and failing to consider that the spawning, propagation and the like of aquatic organisms such as river fish and the like are ecological water requirements; second, the full-cycle natural runoff process of rivers presents a staged character over time, but the hydrologically deduced ecological flow process generally does not have an obvious time staged character.
In view of the above, a full-period suitable ecological flow interval calculating method is provided, and the data of historical runoff, river basin ecology control section fish spawning, propagation and the like are combined to calculate a suitable ecological flow interval range with a stage characteristic, so that an important reference can be provided for reservoir ecological flow scheduling.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a full-period suitable ecological flow interval calculating method, provides a full-period suitable ecological flow interval calculating method and provides technical support for further expanding a hydrologic ecological flow calculating method.
In order to achieve the above purpose, the invention adopts the following technical scheme:
the invention provides a full-period proper ecological flow interval deducing method, which comprises the following steps:
s1, collecting historical natural daily runoff data of ecological control sections of drainage basins, providing a cumulative flow proportion difference method, calculating the ratio of the daily cumulative flow of each year to the average flow of many years and a second-order backward difference value, and obtaining a year-by-year difference sequence;
s2, calculating an average value of inflection points of the year-by-year differential sequence, and dividing the runoff period into three stages according to the data of the river basin ecological control section, wherein the three stages are an annual stationary period, a fish spawning sensitive period and an annual fluctuation period respectively;
s3, calculating the annual stationary phase before sorting from small to large year by yearFractional number of bits, whereinAnd before orderingThe flow corresponding to the quantile is used as the upper and lower limit ranges of the proper ecological flow in the annual stationary period;
s4, for the fish spawning sensitive period, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish spawning sensitive period by adopting a biological hydrological response model based on the spawning and flow of the fish in the river basin ecological control section;
s5, for the annual fluctuation period, providing a calculation method of a suitable ecological flow interval range based on a scene recourse reduction technology, and calculating the ranges of the upper limit and the lower limit of the suitable ecological flow of the annual fluctuation period;
and S6, according to the S1-S5, summarizing and splicing the intervals with proper ecological flow in each stage, and obtaining the interval range with proper ecological flow in the whole period of the ecological control section of the river basin.
Further, in the step S1, the basin ecological control section history natural daily runoff data specifically includes:
;
wherein,representing a matrix formed by natural daily runoffs of the river basin ecological control section history;represent the firstYear of lifeFlow value per day, unit;Andrepresenting the total number of days and total years of a year, respectively.
Further, in the step S1, the cumulative flow ratio difference method specifically includes:
s101, calculating the accumulated flow ratio according to natural daily runoffs of the ecological control section history, wherein the accumulated flow ratio is specifically as follows:
;
wherein,represent the firstYear of lifeCumulative flow ratio of days;represents the average flow rate of the river basin ecological control section history for years, and the unit;
S102, calculating to obtain a year-by-year differential sequence by adopting a second-order backward difference method, wherein the method specifically comprises the following steps of:
;
wherein,representing a year-by-year differential sequence matrix;represent the firstYear of lifeSecond-order backward differential values of days;
the calculation formula of (2) is as follows:
。
further, in S2, the year-by-year differential sequence inflection point is specifically:
;
wherein,represent the firstA annual differential sequence inflection point group;andrespectively represent the firstThe 1 st and 2 nd inflection points of year areEqual to 0;
in the step S2, the average value of inflection points of the year-by-year differential sequence is calculated according to the following formula:
;
;
wherein,andmean values of inflection points of the 1 st and 2 nd year-by-year differential sequences are respectively represented;representing an upward integer function;
in the S2, the data of the river basin ecological control section is that the spawning period and the propagation period of the fish are the firstDay to dayThe day, wherein,;
based onAndand (b)Andthe runoff process of the ecological control section is divided into three stages, namely:
the annual stationary period is from day 1 to dayTianhe (Chinese character)Day to dayA day;
the spawning sensitive period of fish, the firstDay to dayDay, annual fluctuation period is the firstDay to dayDay, the (b)Day to dayAnd (3) days.
Further, in the step S3, the historical natural runoff in the annual stationary period is calculated to form a runoff matrix, which specifically includes:
;
wherein,a runoff matrix which represents the historical natural runoffs in the annual stationary period;
selecting the data of each row of runoffs in the historical natural runoff matrix in the annual stationary period before sorting according to the sorting from small to largeWhereinAndthe upper limit and the lower limit of the proper ecological flow in the annual stationary phase are formed by the following steps:
;
;
wherein,andrespectively representing the upper limit and the lower limit of the proper ecological flow in the annual stationary phase;andrespectively represent the stationary phase in the yearBefore the sorting from small to large year by yearBefore the quantile and the sortingQuantile flow value, unit。
Further, the fish spawning and flow rate of the river basin ecological control section of the S4 is specifically:
;
;
wherein,a flow matrix for representing spawning sensitive period of fish in the river basin ecological control section;an oviposition amount matrix for representing the oviposition sensitive period of the fish in the river basin ecological control section;represent the firstYear of lifeEcologically controlling spawning amount of sectional fishes in a tenna river basin, wherein the spawning amount is ten thousand per ten thousand grains;
the biological hydrological response model is adopted, and specifically comprises the following steps:
;
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinA sequence of day history traffic;representing the spawning sensitive period of the fish with the ecological control section of the river basinA daily historical spawning amount sequence;representing the spawning sensitive period of the fish with the ecological control section of the river basinFitting relation between spawning amount and corresponding flow of fish in the day.
Further, in the step S4, the upper limit and the lower limit of the suitable ecological flow rate of the spawning sensitive period of the fish are calculated, specifically:
s401, statistics of historical basin ecological control sectionThe maximum and minimum values of spawning amount of the main fish in the day are specifically as follows:
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinThe maximum spawning amount of the fish is mainly the highest in every ten thousand grains;representing the spawning sensitive period of the fish with the ecological control section of the river basinMinimum value of spawning quantity of fish in daily life, unit;
S402, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish in the river basin ecological control section in the spawning sensitive period, wherein the range is specifically as follows:
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinUpper limit of ecological flow rate suitable for daily use, unit;Representing the spawning sensitive period of the fish with the ecological control section of the river basinLower limit of ecological flow suitable for daily use, unit。
Further, in S5, the flow data of the annual fluctuation period specifically includes:
for the purpose ofThe scene back-substitution reduction technology is adopted, and specifically comprises the following steps:
;
;
;
wherein,representing a scene back-substitution reduction model;representing the reduced scene corresponding probability sequence;representing the post-reduction thScene correspondence probability;representing the reduced traffic scene matrix;representing a reduced scene number;
select scene after downscalingThe maximum value and the minimum value of the daily flow are used as the upper limit and the lower limit range of the proper ecological flow in the annual fluctuation period, and specifically:
;
;
wherein,andrespectively represent the annual fluctuation periodUpper and lower limits of suitable ecological flow per day, units。
Further, in the step S6, the interval range of the full period suitable ecological flow of the river basin ecological control section is specifically:
;
;
wherein,andrespectively represent the upper limit and the lower limit of the full period suitable ecological flow of the ecological control section of the river basin.
The beneficial effects of the invention are as follows: the method for calculating the full-period suitable ecological flow interval is provided for the first time, the historical runoff data of the river basin ecological control section are combined, the spawning and propagation periods of fish with the ecological control section are considered, the annual runoff process is divided in stages, the upper limit range and the lower limit range of the suitable ecological flow are calculated by adopting various methods, the defect of a hydrologic ecological flow calculation method can be overcome, and important references are provided for the development of ecological flow scheduling of a reservoir.
Drawings
FIG. 1 is a flow chart of a method for calculating a full-period suitable ecological flow interval according to the present invention;
FIG. 2 is a schematic diagram of a runoff staging of an embodiment of the present invention;
fig. 3 is a schematic diagram of a full-period suitable ecological flow interval of a river basin ecological control section according to an embodiment of the invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Referring to fig. 1, taking a river basin ecological control section of the Yangtze river as an example, a method for calculating a full-period suitable ecological flow interval is specifically described, and includes the following steps:
s1, collecting historical natural daily runoff data of ecological control sections of drainage basins, providing a cumulative flow proportion difference method, calculating the ratio of the daily cumulative flow of each year to the average flow of many years and a second-order backward difference value, and obtaining a year-by-year difference sequence;
the natural daily runoff data of the river basin ecological control section history specifically comprises the following steps:
;
wherein,representing a matrix formed by natural daily runoffs of the river basin ecological control section history;represent the firstYear of lifeFlow value per day, unit;Andrespectively represent the total day of one yearNumber and total years; for the situation of year wetting, deleting the flow data of 29 days of 2 months, and ensuring the regularity of the matrix;
the cumulative flow ratio difference method specifically comprises the following steps:
s101, calculating the accumulated flow ratio according to natural daily runoffs of the ecological control section history, wherein the accumulated flow ratio is specifically as follows:
;
wherein,represent the firstYear of lifeCumulative flow ratio of days;represents the average flow rate of the river basin ecological control section history for years,
s102, calculating to obtain a year-by-year differential sequence by adopting a second-order backward difference method, wherein the method specifically comprises the following steps of:
;
wherein,representing a year-by-year differential sequence matrix;represent the firstYear of lifeSecond-order backward differential values of days;
the calculation formula of (2) is as follows:
;
s2, calculating an average value of inflection points of the year-by-year differential sequence, and dividing the runoff period into three stages according to the data of the river basin ecological control section, wherein the three stages are an annual stationary period, a fish spawning sensitive period and an annual fluctuation period respectively;
in the step S2, the year-by-year differential sequence inflection points are specifically:
;
wherein,represent the firstA annual differential sequence inflection point group;andrespectively represent the firstThe 1 st and 2 nd inflection points of year areEqual to 0;
the average value of inflection points of the year-by-year differential sequence is calculated as follows:
;
;
wherein, based on the year-by-year differential sequence inflection point value obtained by S1, the average value of the inflection points of the differential sequence for many years is calculated, namely=97 sum=301 represents the average value of the inflection points of the 1 st and 2 nd year-by-year differential sequences, respectively; this is similar to the flood season of the basin in this example, 4 months to 9 months, and the non-flood season divided by 10 months to 3 months of the next year;representing an upward integer function;
in the S2, the spawning and propagation period of main fishes in the river basin ecological control section is the firstDay to dayIn the days, the total weight of the product,
wherein, according to the data such as the comprehensive planning of the river basin and the comprehensive dispatching report of the reservoir group of the river basin, the spawning sensitive period of the main fishes in the ecological control section of the river basin is 5 late to 7 early days,Based onAndand (b)Andcan control the section of the ecologyIs divided into three stages, namely: the annual stationary phase is from day 1 to dayTianhe (first)Day to dayThe spawning sensitive period of the fish is dayDay to dayDay, annual fluctuation period is the firstDay to dayDay, the (b)Day to dayAnd (3) days.
S3, for the annual stationary phase, calculating the historical natural runoffs before sorting from small to large year by yearFractional number of bits, whereinAnd before orderingThe flow corresponding to the quantile is used as the upper and lower limit ranges of the proper ecological flow in the annual stationary period;
in the step S3, calculating the historical natural runoff in the annual stationary period to form a runoff matrix, wherein the method specifically comprises the following steps:
;
wherein,a runoff matrix which represents the historical natural runoffs in the annual stationary period;
selecting the data of each row of runoffs in the historical natural runoff matrix in the annual stationary period before sorting according to the sorting from small to largeWhereinAndthe upper limit and the lower limit of the proper ecological flow in the annual stationary phase are formed by the following steps:
;
;
wherein,andrespectively representing the upper limit and the lower limit of the proper ecological flow in the annual stationary phase;andrespectively represent the stationary phase in the yearBefore the sorting from small to large year by yearBefore the quantile and the sortingQuantile flow value, unit。
S4, for the fish spawning sensitive period, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish spawning sensitive period by adopting a biological hydrological response model based on the spawning and flow of the fish in the river basin ecological control section;
the fish spawning and flow of the river basin ecological control section of the S4 are specifically:
;
;
wherein,a flow matrix for representing spawning sensitive period of fish in the river basin ecological control section;an oviposition amount matrix for representing the oviposition sensitive period of the fish in the river basin ecological control section;represent the firstYear of lifeEcologically controlling spawning amount of sectional fishes in a tenna river basin, wherein the spawning amount is ten thousand per ten thousand grains;
the biological hydrological response model is adopted, and specifically comprises the following steps:
;
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinA sequence of day history traffic;representing the spawning sensitive period of the fish with the ecological control section of the river basinA daily historical spawning amount sequence;representing the spawning sensitive period of the fish with the ecological control section of the river basinFitting relation between spawning amount and corresponding flow of fish in the day.
In the step S4, calculating the range of the upper limit and the lower limit of the proper ecological flow of the spawning sensitive period of the fish, specifically:
s401, statistics of historical basin ecological control sectionThe maximum and minimum values of spawning amount of the main fish in the day are specifically as follows:
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinThe maximum spawning amount of the fish is mainly the highest in every ten thousand grains;representing the spawning sensitive period of the fish with the ecological control section of the river basinMinimum value of spawning quantity of fish in daily life, unit;
S402, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish in the river basin ecological control section in the spawning sensitive period, wherein the range is specifically as follows:
;
;
wherein,representing the spawning sensitive period of the fish with the ecological control section of the river basinUpper limit of ecological flow rate suitable for daily use, unit;Representing the spawning sensitive period of the fish with the ecological control section of the river basinLower limit of ecological flow suitable for daily use, unit。
S5, for the annual fluctuation period, providing a calculation method of a suitable ecological flow interval range based on a scene recourse reduction technology, and calculating the ranges of the upper limit and the lower limit of the suitable ecological flow of the annual fluctuation period;
in the step S5, flow data in the annual fluctuation period specifically includes:
for the purpose ofThe scene back-substitution reduction technology is adopted, and specifically comprises the following steps:
;
;
;
wherein,representing a scene back-substitution reduction model;representing the reduced scene corresponding probability sequence;representing the post-reduction thScene correspondence probability;representing the reduced traffic scene matrix;representing a reduced scene number; here, the;
Select scene after downscalingThe maximum value and the minimum value of the daily flow are used as the upper limit and the lower limit range of the proper ecological flow in the annual fluctuation period, and specifically:
;
;
wherein,andrespectively represent the annual fluctuation periodUpper and lower limits of suitable ecological flow per day, units。
S6, according to the S1-S5, summarizing and splicing the sections with proper ecological flow in each stage to obtain the section range with proper ecological flow in the whole period of the ecological control section of the river basin;
in the step S6, the interval range of the full period suitable ecological flow of the river basin ecological control section is specifically:
;
;
wherein,andrespectively represent the upper limit and the lower limit of the full period suitable ecological flow of the ecological control section of the river basin.
The foregoing examples merely illustrate embodiments of the invention and are described in more detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention. Accordingly, the scope of protection of the present patent is to be determined by the appended claims.
Claims (4)
1. The method for calculating the full-period proper ecological flow interval is characterized by comprising the following steps of:
s1, collecting historical natural daily runoff data of ecological control sections of drainage basins, providing a cumulative flow proportion difference method, calculating the ratio of the daily cumulative flow of each year to the average flow of many years and a second-order backward difference value, and obtaining a year-by-year difference sequence;
s2, calculating an average value of inflection points of the year-by-year differential sequence, and dividing the runoff period into three stages according to the data of the river basin ecological control section, wherein the three stages are an annual stationary period, a fish spawning sensitive period and an annual fluctuation period respectively;
s3, for the annual stationary phase, calculating the historical natural runoffs before sorting from small to large year by yearFractional number of bits, whereinAnd before ordering->The flow corresponding to the quantile is used as the upper and lower limit ranges of the proper ecological flow in the annual stationary period;
s4, for the fish spawning sensitive period, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish spawning sensitive period by adopting a biological hydrological response model based on the spawning and flow of the fish in the river basin ecological control section;
s5, for the annual fluctuation period, providing a calculation method of a suitable ecological flow interval range based on a scene recourse reduction technology, and calculating the ranges of the upper limit and the lower limit of the suitable ecological flow of the annual fluctuation period;
s6, according to the S1-S5, summarizing and splicing the sections with proper ecological flow in each stage to obtain the section range with proper ecological flow in the whole period of the ecological control section of the river basin;
in the S1, the natural daily runoff data of the river basin ecological control section history is specifically:
;
wherein,representing a matrix formed by natural daily runoffs of the river basin ecological control section history; />Indicate->Year of lifeDay flow value, unit->;/>And->Representing the total number of days and total years of a year, respectively;
in the step S1, the accumulated flow ratio difference method specifically comprises the following steps:
s101, calculating the accumulated flow ratio according to natural daily runoffs of the ecological control section history, wherein the accumulated flow ratio is specifically as follows:
;
wherein,indicate->Age->Cumulative flow ratio of days; />Represents average flow rate of river basin ecological control section history for years, unit +.>;
S102, calculating to obtain a year-by-year differential sequence by adopting a second-order backward difference method, wherein the method specifically comprises the following steps of:
;
wherein,representing a year-by-year differential sequence matrix; />Indicate->Age->Second-order backward differential values of days;
the calculation formula of (2) is as follows:
;
in the step S5, flow data in the annual fluctuation period specifically includes:
;
for the purpose ofThe scene back-substitution reduction technology is adopted, and specifically comprises the following steps:
;
;
;
wherein,representing a scene back-substitution reduction model; />Representing the reduced scene corresponding probability sequence; />Representing the reduced->Scene correspondence probability; />Representing the reduced traffic scene matrix; />Representing a reduced scene number;
select scene after downscalingThe maximum value and the minimum value of the daily flow are used as the upper limit and the lower limit range of the proper ecological flow in the annual fluctuation period, and specifically:
;
;
wherein,and->Respectively represent the annual fluctuation period/>Upper and lower limits of suitable ecological flow per day, units;
S6, the interval range of the full period suitable ecological flow of the river basin ecological control section is specifically:
;
;
wherein,
and->Respectively representing the upper limit and the lower limit of the full period suitable ecological flow of the ecological control section of the river basin;
in the step S2, the year-by-year differential sequence inflection points are specifically:
;
wherein,indicate->A annual differential sequence inflection point group; />And->Respectively represent +.>1 st and 2 nd inflection points of +.>Equal to 0;
in the step S2, the average value of inflection points of the year-by-year differential sequence is calculated according to the following formula:
;
;
wherein,and->Mean values of inflection points of the 1 st and 2 nd year-by-year differential sequences are respectively represented; />Representing an upward integer function;
in the S2, the data of the river basin ecological control section is that the spawning period and the propagation period of the fish are the firstDay to->The day, wherein,,/>;
based onAnd->And->And->The runoff process of the ecological control section is divided into three stages, namely:
the annual stationary period is from day 1 to dayTianhe->Day to->A day;
the spawning sensitive period of fish, the firstDay to->Day, intra-annual fluctuation period is +.>Day to->Day, the (b)Day to->And (3) days.
2. The method for calculating the full-period suitable ecological flow interval according to claim 1, wherein in S3, the historical natural runoff in the annual stationary period is calculated to form a runoff matrix, and specifically:
;
wherein,a runoff matrix which represents the historical natural runoffs in the annual stationary period;
selecting the data of each row of runoffs in the historical natural runoff matrix in the annual stationary period before sorting according to the sorting from small to largeFractional number, wherein->And before ordering->The flow of the fractional digits constitutes the upper limit and the lower limit of the proper ecological flow in the annual stationary phase, and specifically comprises the following steps:
;
;
wherein,and->Respectively representing the upper limit and the lower limit of the proper ecological flow in the annual stationary phase; />Andrespectively represent the +.>Before the daily sorting from small to large year->Quantile and before ordering->Quantile flow value, unit->。
3. The method for calculating the full-period suitable ecological flow interval according to claim 2, wherein the fish spawning and flow of the river basin ecological control section of S4 is specifically:
;
;
wherein,a flow matrix for representing spawning sensitive period of fish in the river basin ecological control section; />An oviposition amount matrix for representing the oviposition sensitive period of the fish in the river basin ecological control section; />Indicate->Age->Ecologically controlling spawning amount of sectional fishes in a tenna river basin, wherein the spawning amount is ten thousand per ten thousand grains;
the biological hydrological response model is adopted, and specifically comprises the following steps:
;
;
;
wherein,represents the sensitive spawning period of the fish with the ecological control section in the river basin>A sequence of day history traffic; />Represents the sensitive spawning period of the fish with the ecological control section in the river basin>A daily historical spawning amount sequence; />Represents the sensitive spawning period of the fish with the ecological control section in the river basin>Fitting relation between spawning amount and corresponding flow of fish in the day.
4. A method for estimating a full-cycle fit ecological flow interval according to claim 3, characterized in that: in the step S4, calculating the range of the upper limit and the lower limit of the proper ecological flow of the spawning sensitive period of the fish, specifically:
s401, statistics of historical basin ecological control sectionThe maximum and minimum values of spawning amount of the main fish in the day are specifically as follows:
;
;
wherein,represents the sensitive spawning period of the fish with the ecological control section in the river basin>The maximum spawning amount of the fish is mainly the highest in every ten thousand grains; />Represents the sensitive spawning period of the fish with the ecological control section in the river basin>Minimum value of spawning quantity of fish in day main, unit +.>;
S402, calculating the range of the upper limit and the lower limit of the proper ecological flow of the fish in the river basin ecological control section in the spawning sensitive period, wherein the range is specifically as follows:
;
;
wherein,represents the sensitive spawning period of the fish with the ecological control section in the river basin>Upper limit of ecological flow rate suitable for daily use, unit;/>Represents the sensitive spawning period of the fish with the ecological control section in the river basin>Lower limit of ecological flow suitable for daily use, unit。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410030698.5A CN117541031B (en) | 2024-01-09 | 2024-01-09 | Full-period proper ecological flow interval deducing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410030698.5A CN117541031B (en) | 2024-01-09 | 2024-01-09 | Full-period proper ecological flow interval deducing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117541031A CN117541031A (en) | 2024-02-09 |
CN117541031B true CN117541031B (en) | 2024-03-29 |
Family
ID=89788477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410030698.5A Active CN117541031B (en) | 2024-01-09 | 2024-01-09 | Full-period proper ecological flow interval deducing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117541031B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348083A (en) * | 2019-06-26 | 2019-10-18 | 长江水利委员会长江科学院 | A kind of base flow adds the fish egg-laying season ecological flow design method of pulse |
CN113688497A (en) * | 2021-07-05 | 2021-11-23 | 郑州大学 | Ecological flow quantification method based on ecological-habitat-flow grouping response mechanism |
CN116245042A (en) * | 2022-12-22 | 2023-06-09 | 河海大学 | Ecological water demand calculating method based on hydrodynamic force-water environment-aquatic organisms |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111079085B (en) * | 2019-12-05 | 2023-03-24 | 西安理工大学 | Ecological flow determination method considering lifting amount |
-
2024
- 2024-01-09 CN CN202410030698.5A patent/CN117541031B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348083A (en) * | 2019-06-26 | 2019-10-18 | 长江水利委员会长江科学院 | A kind of base flow adds the fish egg-laying season ecological flow design method of pulse |
CN113688497A (en) * | 2021-07-05 | 2021-11-23 | 郑州大学 | Ecological flow quantification method based on ecological-habitat-flow grouping response mechanism |
CN116245042A (en) * | 2022-12-22 | 2023-06-09 | 河海大学 | Ecological water demand calculating method based on hydrodynamic force-water environment-aquatic organisms |
Non-Patent Citations (3)
Title |
---|
Type and timing of stream flow changes in urbanizing watersheds in the eastern U.S.;Kristina G. Hopkins, 等;《Elementa: Science of the Anthropocene》;20150611;第3卷(第000056期);第1-14页 * |
基于河流鱼类适宜生境控制的梯级水库优化调度方法研究;王琲;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20180615(第6期);第C037-15页 * |
综合物理栖息地与流量脉冲的二松干流生态需水过程研究;杨泽凡;《中国优秀硕士学位论文全文数据库基础科学辑》;20160315(第3期);第A006-4页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117541031A (en) | 2024-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Berhe et al. | MODSIM-based water allocation modeling of Awash River Basin, Ethiopia | |
Yilmaz et al. | An indicator based assessment for water resources management in Gediz River Basin, Turkey | |
Xia et al. | Water resources vulnerability and adaptive management in the Huang, Huai and Hai river basins of China | |
Zhang et al. | Using a hierarchical model framework to assess climate change and hydropower operation impacts on the habitat of an imperiled fish in the Jinsha River, China | |
CN107092786B (en) | Ecological base flow calculation method and system considering different states of river | |
Ely et al. | Dam-induced hydrologic alterations in the rivers feeding the Pantanal | |
Wang et al. | Reservoir ecological operation considering outflow variations across different time scales | |
Chaudhary et al. | Integrated 1D and 2D numerical model simulations for flushing of sediment from reservoirs | |
Sreelash et al. | Changes in the Hydrological Characteristics of Cauvery River draining the eastern side of southern Western Ghats, India | |
Loinaz et al. | Modeling ecohydrological impacts of land management and water use in the Silver Creek basin, Idaho | |
CN117541031B (en) | Full-period proper ecological flow interval deducing method | |
Guo et al. | Evaluation of hydrological regime alteration and ecological effects in the middle and lower of the Yangtze River, China | |
CN113688497B (en) | Ecological flow quantification method based on ecological-habitat-flow grouping response mechanism | |
Liang et al. | Evaluating future water security in the upper Yangtze River Basin under a changing environment | |
Li et al. | Incorporating ecological adaptation in a multi-objective optimization for the Three Gorges Reservoir | |
CN116542459A (en) | Ecological scheduling method for remodeling spawning ground habitat of fishes in variable water return area | |
Park et al. | Estimation of least ecological streamflow through long-term habitat evaluation using stream hydrodynamics and water quality factors | |
Zhou et al. | Uncertainty analysis for evaluating flow regime alteration of Jinsha River based on Indicators of Hydrologic Alteration | |
CN114037248A (en) | River ecological risk assessment method based on ecological deficit index | |
Wang et al. | Multi-scale flow regimes and driving forces analysis based on different models: a case study of the Wu River basin | |
Nguyen | Harmonizing multi-sectorial water management with minimum flow requirements in an anthropogenically impacted river basin. The case of Vu Gia–Thu Bon, Central Viet Nam | |
CN116993184A (en) | Water resource shortage assessment method, system and computer readable storage medium | |
CN118095581A (en) | Step small water power water inflow prediction method and system based on distributed hydrologic model | |
Frazier et al. | Effects of flow regulation in flow regime on the Murrumbidgee River, south eastern Australia: an assessment using a daily estimation hydrological model | |
CN115841192A (en) | Flood season staging method and system based on information entropy theory |
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 |