CN110578317A - hydrological model reservoir discharge capacity simulation method - Google Patents
hydrological model reservoir discharge capacity simulation method Download PDFInfo
- Publication number
- CN110578317A CN110578317A CN201910860489.2A CN201910860489A CN110578317A CN 110578317 A CN110578317 A CN 110578317A CN 201910860489 A CN201910860489 A CN 201910860489A CN 110578317 A CN110578317 A CN 110578317A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- capacity
- discharge
- initial
- time period
- 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.)
- Granted
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
- E02B1/02—Hydraulic models
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a hydrological model reservoir discharge capacity simulation method, which belongs to the technical field of distributed hydrological model simulation application and comprises the following steps: collecting basic information of the reservoir, simulating the inflow rate and the evaporation capacity of the reservoir in the current time period by adopting a hydrological model, and calculating the initial reservoir capacity and the over-limit lower discharge capacity exceeding the safety limit reservoir capacity in the current time period; correcting the discharge amount exceeding the safety limit storage capacity to obtain the initial discharge amount at the end of the time period, and further calculating the initial storage capacity of the reservoir at the end of the time period; and correcting the initial reservoir capacity and the lower discharge of the reservoir obtained at the end of the time period by adopting the reservoir capacity information to obtain the final reservoir capacity and the lower discharge of the reservoir at the end of the current time period, and recording relevant state information to obtain a final simulation result. The method adopts limited reservoir parameter information to simulate the reservoir water storage and drainage processes, and improves the simulation precision of the hydrological model; the whole reservoir dispatching simulation process is simple, and the program improvement is easy to realize.
Description
Technical Field
The invention relates to the technical field of distributed hydrological model simulation application, in particular to a hydrological model reservoir discharge capacity simulation method.
background
the distributed hydrological model is an effective means for exploring and recognizing complex hydrological cycle processes and mechanisms and is also an effective tool for solving a plurality of hydrological practical problems. With the increase of human activities, the influence on natural water circulation is increased gradually. The natural water circulation process can be expressed by a mathematical function and has certain regularity. However, the influence of human activities is not mathematical and has a certain randomness. Reservoir construction is one of the most important influencing factors influencing the runoff process of the river channel, and has the influence of enlarging and increasing the runoff of the river channel. In the course of simulating river course flow, if the influence of reservoir dispatching is not considered, the simulation result is far from the actual situation. In the process of simulating reservoir dispatching, information such as specific dispatching rules of each reservoir, water level-reservoir capacity-discharge capacity curves of the reservoirs and the like needs to be known. However, in practical application, the information collection is difficult, and even no relevant information exists, so that the complete reservoir dispatching simulation process is difficult. When the hydrological model is used for basin water circulation process simulation, the annual scale and even the perennial scale are generally adopted, so if a simpler annual scale reservoir dispatching simulation rule is provided to enable the annual runoff process simulation to be relatively accurate, the model simulation efficiency can be effectively improved, and the accuracy is improved.
Disclosure of Invention
the invention aims to provide a hydrological model reservoir discharge amount simulation method, so that the problems in the prior art are solved.
in order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a hydrological model reservoir discharge amount simulation method comprises the following steps:
S1, collecting basic reservoir information including 5 pieces of reservoir capacity information and/or 3 pieces of control water level maximum discharge capacity and/or 2 month-by-month average discharge capacities and/or economic social water intake;
s2, simulating the inflow rate and the evaporation capacity of the reservoir in the current time period by adopting the existing any hydrological model, and calculating the initial reservoir capacity and the over-limit discharge capacity exceeding the safety limit reservoir capacity in the current time period;
s3, correcting the discharge amount exceeding the safety limit storage capacity to ensure that the maximum discharge amount does not exceed the safety limit storage capacity, obtaining the initial discharge amount at the end of the time interval, and further calculating the initial storage capacity of the reservoir at the end of the time interval;
And S4, correcting the calculated initial reservoir capacity and the calculated lower discharge capacity of the reservoir at the end of the time period by using the reservoir capacity information to obtain the final reservoir capacity and the calculated lower discharge capacity of the reservoir at the end of the current time period, and recording related state information to obtain a final simulation result.
preferably, in step S1:
The used storage capacity is corresponding to 5 control water levels of the reservoir, and the arrangement sequence of the storage capacity in the front completely comprises the storage capacity in the rear according to the relation that: total reservoir capacity, design flood level reservoir capacity, normal water storage level reservoir capacity, flood control limit level reservoir capacity and dead reservoir capacity;
The maximum water level discharge control flow is used for controlling the maximum discharge flow when the water level is positioned in the corresponding water level line interval, and comprises checking the maximum discharge flow below the flood level, designing the maximum discharge flow below the flood level and the maximum discharge flow below the normal water storage level; wherein the maximum discharge below the check flood level is a maximum discharge above a design flood level; the maximum discharge below the design flood level is the maximum discharge of the water level below the design flood level and above the normal water storage level; the maximum discharge below the normal water storage level is the maximum discharge above the dead water level when the water level is lower than the normal water storage level; the maximum discharge rate below the dead water level is 0; and checking the maximum discharge flow below the flood level > the maximum discharge flow below the design flood level > the maximum discharge flow below the normal water storage level.
if the information of the water level-reservoir capacity-discharge capacity curve of the reservoir can be collected, the reservoir capacity and discharge capacity of a specific water level can be directly calculated from the curve, and the generalized processing is replaced.
Preferably, the monthly average discharge flow in the step S1 includes a power generation discharge flow or a minimum ecological flow, and the power generation discharge flow and the minimum ecological flow are monthly average values, which represent the conditions of monthly power generation discharge and ecological water demand;
The water intake quantity of the economic society is the water quantity for industry, agriculture and life which is obtained from a reservoir; the part of water is generally year-scale, needs to be scaled down, and meets the requirement of the minimum time scale of the model.
preferably, S2 includes the steps of:
s201, the simulation method of the inflow rate of the reservoir in the current time period comprises the following steps: simulating and calculating all river flows imported into the current reservoir on the same day by adopting a hydrological model or a statistical equation to serve as reservoir inflow;
The simulation method for the large evaporation capacity of the reservoir comprises the following steps: and calculating the water surface evaporation capacity by using a Peneman formula to serve as the maximum value of reservoir evaporation.
S202, calculating the initial reservoir capacity in the current time period, wherein the calculation formula is as follows: the initial reservoir capacity of the current time interval is the last reservoir capacity of the last time interval plus the river inflow amount of the current time interval, the evaporation amount of the reservoir of the current time interval and the economic and social water intake amount of the current time interval;
S203, safely limiting the storage capacity to be the storage capacity corresponding to the limited water level specified by the reservoir safety, wherein the limited storage capacity in the flood season is the flood control limited water level, and the limited water level in the non-flood season is the normal water storage level; and comparing the limited storage capacity with the initial storage capacity of the initial reservoir at the end of the time period, if the safety of the initial storage capacity of the initial reservoir at the end of the time period exceeds the limited storage capacity, the overrun is equal to the initial storage capacity of the initial reservoir at the end of the time period and the limited storage capacity, and if the overrun is not equal to 0.
Preferably, step S3 includes the steps of:
s301, firstly, correcting the overrun lower discharge amount calculated in step S2 by using the minimum ecological flow amount, so that the lower reservoir discharge amount is not less than the minimum ecological discharge amount in the month, that is, if the overrun lower discharge amount is less than the minimum ecological flow amount, the lower reservoir discharge amount in the current step is the minimum ecological discharge amount, otherwise, the lower reservoir discharge amount in the current step is the overrun lower discharge amount in step S2;
S302, then, the power generation discharge amount is used to correct the lower reservoir discharge amount calculated in step S301, so that the lower reservoir discharge amount is not less than the power generation discharge amount in the month, that is, if the lower reservoir discharge amount calculated in step S301 is less than the power generation discharge amount, the current step calculates the lower reservoir discharge amount as the power generation discharge amount, otherwise, the current step calculates the lower discharge amount as the power generation discharge amount in step S301; it should be noted that if the power generation drainage amount is not provided, the present step is not executed;
and S303, finally, calculating to obtain the initial downward discharge amount at the end of the time period.
Preferably, step S303 includes the steps of:
Comparing the initial reservoir capacity with 5 reservoir characteristic storage capacities at the end of the time interval to obtain a storage capacity interval in which the current storage capacity is positioned, and calculating to obtain the maximum discharge control quantity of the corresponding storage capacity; if a curve of the water level, the reservoir capacity and the discharging capacity of the reservoir is collected, the maximum discharge control quantity can be directly calculated according to the initial reservoir capacity at the end of the time period through the curve; and correcting the initial lower discharge amount at the end of the reservoir period calculated in the step S302 by using the maximum discharge control amount, so that the lower discharge amount of the reservoir is not more than the maximum discharge control amount.
Preferably, step S4 includes the steps of:
S401, when the calculated initial reservoir capacity at the end of the time interval is larger than the total reservoir capacity, correcting the lower reservoir discharge amount to be the lower reservoir discharge amount + the initial reservoir capacity at the end of the time interval to be the total reservoir capacity, setting the initial reservoir capacity at the end of the time interval to be the total reservoir capacity, recording the information of the overflow state of the reservoir, and ending the step S4; otherwise, entering step S402;
s402, when the initial reservoir capacity is larger than the full reservoir capacity at the end of the calculated time period, ending the step S4; otherwise, go to step S403;
s403, when the initial reservoir lower discharge capacity is larger than the difference between the dead reservoir capacity and the initial reservoir capacity at the end of the time period, reducing the reservoir lower discharge capacity, complementing the difference for correction, namely correcting the reservoir lower discharge capacity as the reservoir lower discharge capacity + the initial reservoir capacity at the end of the time period-the dead reservoir capacity, setting the final reservoir capacity at the end of the time period as the dead reservoir capacity, and ending the step S4; otherwise, correcting the initial reservoir capacity at the end of the time interval as the initial reservoir capacity at the end of the time interval and the lower reservoir discharge, setting the final lower reservoir discharge as 0, and entering step S404;
S404, when the calculated final initial reservoir capacity at the end of the correction period is greater than 0, the final reservoir capacity at the end of the correction period is equal to the initial reservoir capacity at the end of the correction period, and the step S4 is ended; otherwise, go to step S405;
s405, when the evaporation capacity of the reservoir can meet the negative shortage of the storage capacity of the reservoir, correcting the evaporation capacity of the reservoir to be the evaporation capacity of the reservoir + the initial storage capacity of the reservoir at the end of the time period, setting the final storage capacity of the reservoir at the end of the time period to be 0, recording the dry and anhydrous state information of the reservoir, and ending the step S4; otherwise, correcting the final reservoir capacity at the end of the time interval as the initial reservoir capacity at the end of the time interval and the evaporation capacity of the reservoir, setting the evaporation capacity of the reservoir as 0, and entering the step S406;
s406, when the economic social water intake quantity can meet the negative shortage of the reservoir capacity, correcting the economic social water intake quantity to be the economic social water intake quantity plus the initial reservoir capacity at the end of the time period, setting the final reservoir capacity at the end of the time period to be 0, recording the economic social water intake shortage quantity to be the absolute value of the initial reservoir capacity at the end of the time period, recording the information of the dry and anhydrous state of the reservoir, and ending the step S4; otherwise, correcting the final reservoir capacity at the end of the time interval, namely the initial reservoir capacity at the end of the time interval and the water intake of the economic society, correcting the water intake shortage of the economic society, namely the water intake of the economic society, setting the water intake of the economic society to be 0, and recording the negative reservoir capacity state information of the reservoir.
the invention has the beneficial effects that:
The invention provides a hydrological model reservoir dispatching simulation method, which adopts limited reservoir parameter information to simulate reservoir water storage and drainage processes and improves hydrological model simulation precision; the whole reservoir dispatching simulation process is simple, and the program improvement is easy to realize.
drawings
FIG. 1 is a flow chart of a hydrological model reservoir scheduling simulation method provided by the invention;
FIG. 2 is the reservoir parameter information that the present invention needs to collect;
FIG. 3 is a logic diagram for determining the lower reservoir discharge capacity and the actual discharge capacity of the reservoir;
fig. 4 is a flow chart of a reservoir scheduling simulation program in embodiment 2.
Detailed Description
in order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention.
Examples
as shown in fig. 1, an embodiment of the present invention provides a method for simulating a reservoir discharge amount of a hydrological model, including the following steps:
s1, collecting basic information of the reservoir according to information such as design files or actual operation reports of the reservoir, wherein the basic information mainly comprises 5 pieces of water level reservoir capacity information (total reservoir capacity, design flood water level reservoir capacity, normal water storage level reservoir capacity, flood control limit water level reservoir capacity and dead reservoir capacity), 3 pieces of control water level maximum discharge capacity (maximum discharge capacity below check flood level, maximum discharge capacity below design flood level and maximum discharge capacity below normal water storage level), 2 pieces of monthly average discharge capacity (power generation discharge capacity and minimum ecological discharge capacity) and economic and social water intake capacity;
S2, simulating the inflow rate and the evaporation capacity of the reservoir in the current time period by adopting a model, and calculating the initial reservoir capacity and the lower discharge capacity exceeding the safety limit capacity in the current time period;
s3, correcting the leakage exceeding the safety limit storage capacity according to the monthly power generation leakage flow, the minimum ecological flow and the control water level maximum leakage flow to obtain the initial leakage flow at the end of the time period;
and S4, correcting the calculated reservoir capacity and the lower discharge capacity of the reservoir at the end of the time period by adopting the information such as the total reservoir capacity and the dead reservoir capacity of the reservoir, obtaining the reservoir capacity and the lower discharge capacity of the reservoir at the end of the current time period, calculating relevant state information and obtaining a final simulation result.
first, according to step S1, basic information of reservoirs is collected, wherein the correlation between 5 reservoir levels is shown in fig. 2, where the total reservoir capacity > design flood level reservoir capacity > normal reservoir level reservoir capacity ≧ flood control limit level reservoir capacity > dead reservoir capacity. The 3 maximum discharge control amounts Qmax1 are the maximum discharge control amounts for which the storage capacity is between the total storage capacity and the design flood level storage capacity, Qmax2 is the maximum discharge control amount for which the storage capacity is between the design flood level storage capacity and the normal storage level storage capacity, and Qmax3 is the maximum discharge control amount for which the storage capacity is between the normal storage level storage capacity and the dead storage capacity. The monthly average power generation drainage amount is collected, and no data can be provided, which indicates that no power generation drainage is performed. And collecting the minimum monthly ecological flow, wherein the flow can be calculated by methods such as Tennet and the like. And collecting the economic and social water, and performing down-scale spreading of the adult water consumption according to the calculated time scale.
Next, the initial reservoir capacity and the amount of water at the overrun level are calculated according to step S2. Here, the time interval initial reservoir capacity is the last time interval final reservoir capacity + the river inflow in the time interval, the evaporation capacity of the reservoir in the time interval, and the economic and social water intake in the time interval; the overrun discharge amount is the storage capacity exceeding the limited water level storage capacity of the reservoir, namely, if the initial time period storage capacity is smaller than the limited storage capacity, the overrun discharge amount is equal to 0, otherwise, the overrun discharge amount is equal to the initial time period storage capacity-the limited storage capacity. It should be noted that, in the flood season, the flood control limited reservoir capacity is used as the limited reservoir capacity, and in the non-flood season, the normal reservoir capacity is used as the limited reservoir capacity
and thirdly, calculating initial reservoir lower discharge amount by adopting the overrun discharge amount, the power generation discharge amount, the minimum ecological base flow and the maximum discharge control amount, and correcting the initial time interval end reservoir capacity according to the time interval initial reservoir capacity. Wherein, the maximum value of the overrun discharge, the power generation discharge and the minimum ecological flow is taken, and then the minimum value is taken together with the maximum discharge control amount to obtain the initial reservoir discharge amount, if shown in fig. 3. It should be noted here that if there is no power generation leakage flow, it is not considered. If there is no maximum bleed-down control amount, maximum bleed-down control is not performed. Judging the maximum discharge control quantity according to the initial reservoir capacity of the time interval, namely if the maximum discharge control quantity is less than the dead reservoir capacity, the maximum discharge control quantity is equal to 0; if the water level is less than the normal water level storage capacity, the water level is equal to Qmax 3; in addition, the range in which the reservoir capacity of the final reservoir and the total reservoir capacity, the design flood level reservoir capacity and the normal water storage level reservoir capacity are judged to be located is adopted, and then the range is obtained by adopting thread interpolation.
and fourthly, correcting the calculated initial reservoir lower discharge capacity and the initial period final reservoir storage capacity according to whether the total storage capacity is exceeded or not and whether the total storage capacity is smaller than the dead storage capacity or not. If the total storage capacity is exceeded, adding the exceeded part into the initial lower discharge of the reservoir to obtain the final lower discharge; if the water discharge capacity is smaller than the dead storage capacity, correcting by reducing the water discharge capacity of the reservoir; and if the corrected storage capacity is less than 0, correcting by reducing the evaporation capacity and the economic and social water intake quantity to make the storage capacity of the reservoir not negative as much as possible.
Example 2
In this embodiment, the method is embedded in the hydrological model based on the method in embodiment 1, so as to realize the simulation of the reservoir scheduling drainage process.
The embodiment mainly lists a simulation flow chart of the reservoir water discharge scheduling process, and as shown in fig. 4, the specific program implementation depends on the language of the embedded model. Wherein V has the unit m3q has the unit m3And the unit of the T is s, the value is the number of seconds in the basic time unit of the model, and the value is used for calculating the conversion of the flow and the volume in a time period, for example, the time T is 86400s in daily simulation.
Firstly, river inflow, reservoir evaporation and economic and social water intake in the current time period are calculated by adopting a hydrological model, and initial reservoir capacity V0 in the current time period is obtained by accumulating the final reservoir capacity in the last time period.
secondly, in order to ensure the reservoir safety, the calculated initial reservoir capacity V0 is used for calculating the part needing to be discharged when exceeding the limited reservoir capacity as the initial water discharge quantity Q0 according to the limited reservoir capacity. If V0 exceeds the restricted capacity, Q0 ═ V0-restricted capacity)/T; if not, Q0 is 0. The non-flood season uses the normal water storage level as the limiting water level, and the flood season uses the flood control limiting water level as the limiting water level.
thirdly, the minimum ecological flow in the month of the period is used as a limiting condition for correction, so that the lower reservoir discharge capacity is not less than the current minimum ecological flow, namely Q1 is max (Q0, minimum ecological flow).
fourthly, the Q1 is corrected according to whether the power generation water release amount is provided or not, and the water release amount Q2 is obtained. If the generated water amount is not provided, Q2 is Q1; if provided, Q2 is max (Q1, generated current).
Fifthly, the drainage quantity Q3 is obtained by correcting the Q2 according to whether the maximum drainage control flow is provided or not. If the maximum bleed flow control is not provided, then Q3 is Q2. If so, judging according to the initial storage capacity V0, and when V0 is lower than the dead storage capacity, Qmax is 0; when V0 is lower than the normal storage level capacity, Qmax is the maximum discharge capacity below the normal storage level; when the V0 is larger than the normal storage level capacity, calculating to obtain Qmax by adopting a linear interpolation method according to the initial reservoir capacity V0, the different characteristic reservoir capacities and the corresponding maximum discharge capacity; then, the correction Q3 is min (Q2, Qmax). It should be noted here that if the water level-reservoir capacity-discharge capacity curve is taken directly when data is collected, the curve is directly used to calculate the maximum discharge capacity Qmax corresponding to the reservoir capacity of V1.
Sixthly, calculating the initial reservoir capacity V1-initial reservoir capacity V0-lower reservoir discharge Q3T at the end of the time interval.
Seventhly, judging whether the reservoir capacity V1 is larger than the total capacity, if so, calculating the lower reservoir discharge Q at the end of the time period as the lower reservoir discharge Q3+ (the reservoir capacity V1-the total capacity)/T at the end of the time period, calculating the reservoir capacity V at the end of the time period as the total capacity, recording the information of the full overflow state of the reservoir, and ending the simulation at the current time period; if not, continuing the subsequent steps.
Eighth, judging whether the reservoir capacity V1 is smaller than the dead reservoir capacity, if not, calculating the final reservoir capacity V which is V1 at the end of the time period, calculating the lower discharge Q which is Q3 at the end of the time period, and ending the current time period simulation; if so, continuing the subsequent steps.
ninth, since the reservoir capacity V1 is smaller than the dead capacity, the dead capacity requirement is satisfied first by reducing the letdown. Normally the dead volume is below that without a leakage flow. If the lower reservoir discharge Q3 can meet the difference of subtracting V1 from the dead storage capacity, calculating the lower reservoir discharge Q at the end of the time interval as Q3+ (V1-dead storage capacity)/T, calculating the reservoir storage capacity V at the end of the time interval as dead storage capacity, and ending the simulation at the current time interval; and if the lower discharge quantity of the reservoir cannot meet the difference, calculating the corrected reservoir capacity V2 (V1 + Q3) T at the end of the time interval, calculating the lower discharge quantity Q of the reservoir (0), and continuing the subsequent steps.
tenth, judging whether the reservoir capacity V2 of the correction reservoir at the end of the time interval is less than 0, if not, calculating the reservoir capacity V at the end of the time interval as V2, and ending the simulation at the current time interval; if still less than 0, the subsequent steps are continued.
Eleventh, the shortage that the reservoir capacity is negative is met by reducing the evaporation of the reservoir. If the evaporation capacity can meet the shortage of the reservoir capacity, correcting the evaporation capacity of the reservoir to be the evaporation capacity of the reservoir and the reservoir capacity V2 at the end of the time period, calculating the reservoir capacity V at the end of the time period to be 0, and recording the dry and anhydrous state information of the reservoir; and if the evaporation capacity can not meet the deficit, calculating the reservoir capacity V3 of the reservoir at the end of the time period as V2+ the evaporation capacity of the reservoir, setting the evaporation capacity of the reservoir equal to 0, and continuing the subsequent steps.
Twelfth, the shortage that the reservoir is negative is met by adopting a mode of reducing the economic and social water intake. If the economic social water intake can meet the shortage of the reservoir capacity, correcting the economic social water intake amount to be the economic social water intake amount plus the reservoir capacity V3 at the end of the time period, calculating the reservoir capacity V at the end of the time period to be 0, recording the shortage of the economic social water intake amount V3, and recording the water consumption and water-free state information of the reservoir; and if the water intake of the economic society cannot meet the shortage, calculating the time interval end reservoir capacity V which is the time interval end reservoir capacity V3 plus the water intake of the economic society, recording the shortage of the water intake of the economic society which is the water intake of the economic society, setting the water intake of the economic society which is 0, and recording the negative reservoir capacity state information of the reservoir. The current session simulation is ended.
Claims (7)
1. A hydrological model reservoir discharge simulation method is characterized by comprising the following steps:
S1, collecting basic reservoir information including 5 pieces of reservoir capacity information and/or 3 pieces of control water level maximum discharge capacity and/or 2 month-by-month average discharge capacities and/or economic social water intake;
S2, simulating the inflow rate and the evaporation capacity of the reservoir in the current time period by adopting a hydrological model, and calculating the initial reservoir capacity and the over-limit discharge capacity exceeding the safety limit reservoir capacity in the current time period;
s3, correcting the discharge amount exceeding the safety limit storage capacity to obtain the initial discharge amount at the end of the time period, and further calculating the initial reservoir storage capacity at the end of the time period;
and S4, correcting the initial reservoir capacity and the lower discharge capacity of the reservoir obtained at the end of the time period by adopting the reservoir capacity information to obtain the final reservoir capacity and the lower discharge capacity at the end of the current time period, and recording relevant state information to obtain a final reservoir discharge capacity simulation result.
2. The method for simulating the reservoir discharge amount of a hydrological model according to claim 1, wherein in step S1:
the used storage capacity is corresponding to 5 control water levels of the reservoir, and the arrangement sequence of the storage capacity in the front completely comprises the storage capacity in the rear according to the relation that: total reservoir capacity, design flood level reservoir capacity, normal water storage level reservoir capacity, flood control limit level reservoir capacity and dead reservoir capacity;
the maximum water level discharge control flow is used for controlling the maximum discharge flow when the water level is positioned in the corresponding water level line interval, and comprises checking the maximum discharge flow below the flood level, designing the maximum discharge flow below the flood level and the maximum discharge flow below the normal water storage level; wherein the maximum discharge below the check flood level is a maximum discharge above a design flood level; the maximum discharge below the design flood level is the maximum discharge of the water level below the design flood level and above the normal water storage level; the maximum discharge below the normal water storage level is the maximum discharge above the dead water level when the water level is lower than the normal water storage level; the maximum discharge rate below the dead water level is 0; and checking the maximum discharge capacity below the flood level > the maximum discharge capacity below the design flood level > the maximum discharge capacity below the normal water storage level;
if the information of the water level-reservoir capacity-discharge capacity curve of the reservoir can be collected, the reservoir capacity and discharge capacity of a specific water level can be directly calculated from the curve, and the generalized processing is replaced.
3. The method for simulating the discharged water amount of the hydrological model reservoir according to claim 1, wherein the monthly average discharge rate in step S1 includes a power generation discharge rate or a minimum ecological flow rate, and the power generation discharge rate and the minimum ecological flow rate are monthly average values, which represent monthly power generation discharge and ecological water demand conditions;
the water intake quantity of the economic society is the water quantity for industry, agriculture and life which is obtained from a reservoir; the part of water is generally year-scale, needs to be scaled down, and meets the requirement of the minimum time scale of the model.
4. The method for simulating the discharged water amount of the hydrological model reservoir according to claim 1, wherein the step S2 comprises the steps of:
S201, the simulation method of the inflow rate of the reservoir in the current time period comprises the following steps: simulating and calculating all river flows imported into the current reservoir on the same day by adopting a hydrological model or a statistical equation to serve as reservoir inflow;
s202, calculating the initial reservoir capacity in the current time period, wherein the calculation formula is as follows: the initial reservoir capacity of the current time interval is the last reservoir capacity of the last time interval plus the river inflow amount of the current time interval, the evaporation amount of the reservoir of the current time interval and the economic and social water intake amount of the current time interval;
S203, safely limiting the storage capacity to be the storage capacity corresponding to the limited water level specified by the reservoir safety, wherein the limited storage capacity in the flood season is the flood control limited water level, and the limited water level in the non-flood season is the normal water storage level; and comparing the limited storage capacity with the initial storage capacity of the initial reservoir at the end of the time period, if the safety of the initial storage capacity of the initial reservoir at the end of the time period exceeds the limited storage capacity, the overrun is equal to the initial storage capacity of the initial reservoir at the end of the time period and the limited storage capacity, and if the overrun is not equal to 0.
5. The method for simulating the discharged water amount of the hydrological model reservoir according to claim 1, wherein the step S3 comprises the steps of:
S301, firstly, correcting the overrun lower discharge amount calculated in step S2 by using the minimum ecological flow amount, so that the lower reservoir discharge amount is not less than the minimum ecological discharge amount in the month, that is, if the overrun lower discharge amount is less than the minimum ecological flow amount, the lower reservoir discharge amount in the current step is the minimum ecological discharge amount, otherwise, the lower reservoir discharge amount in the current step is the overrun lower discharge amount in step S2;
s302, then, the power generation discharge amount is used to correct the lower reservoir discharge amount calculated in step S301, so that the lower reservoir discharge amount is not less than the power generation discharge amount in the month, that is, if the lower reservoir discharge amount calculated in step S301 is less than the power generation discharge amount, the current step calculates the lower reservoir discharge amount as the power generation discharge amount, otherwise, the current step calculates the lower discharge amount as the power generation discharge amount in step S301; it should be noted that if the power generation drainage amount is not provided, the present step is not executed;
And S303, finally, calculating to obtain the initial downward discharge amount at the end of the time period.
6. the method for simulating the reservoir discharge amount of a hydrological model according to claim 5, wherein step S303 comprises the steps of:
Comparing the initial reservoir capacity with 5 reservoir characteristic storage capacities at the end of the time interval to obtain a storage capacity interval in which the current storage capacity is positioned, and calculating to obtain the maximum discharge control quantity of the corresponding storage capacity; if a curve of the water level, the reservoir capacity and the discharging capacity of the reservoir is collected, the maximum discharge control quantity can be directly calculated according to the initial reservoir capacity at the end of the time period through the curve; and correcting the initial lower discharge amount at the end of the reservoir period calculated in the step S302 by using the maximum discharge control amount, so that the lower discharge amount of the reservoir is not more than the maximum discharge control amount.
7. The method for simulating the discharged water amount of the hydrological model reservoir according to claim 1, wherein the step S4 comprises the steps of:
S401, when the calculated initial reservoir capacity at the end of the time interval is larger than the total reservoir capacity, correcting the lower reservoir discharge amount to be the lower reservoir discharge amount + the initial reservoir capacity at the end of the time interval to be the total reservoir capacity, setting the initial reservoir capacity at the end of the time interval to be the total reservoir capacity, recording the information of the overflow state of the reservoir, and ending the step S4; otherwise, entering step S402;
S402, when the initial reservoir storage capacity is larger than or equal to the dead storage capacity at the end of the calculated time period, ending the step S4; otherwise, go to step S403;
s403, when the initial reservoir lower discharge capacity is larger than the difference between the dead reservoir capacity and the initial reservoir capacity at the end of the time period, reducing the reservoir lower discharge capacity, complementing the difference for correction, namely correcting the reservoir lower discharge capacity as the reservoir lower discharge capacity + the initial reservoir capacity at the end of the time period-the dead reservoir capacity, setting the final reservoir capacity at the end of the time period as the dead reservoir capacity, and ending the step S4; otherwise, correcting the initial reservoir capacity at the end of the time interval as the initial reservoir capacity at the end of the time interval and the lower reservoir discharge, setting the final lower reservoir discharge as 0, and entering step S404;
s404, when the calculated final initial reservoir capacity at the end of the correction period is greater than 0, the final reservoir capacity at the end of the correction period is equal to the initial reservoir capacity at the end of the correction period, and the step S4 is ended; otherwise, go to step S405;
S405, when the evaporation capacity of the reservoir can meet the negative shortage of the storage capacity of the reservoir, correcting the evaporation capacity of the reservoir to be the evaporation capacity of the reservoir + the initial storage capacity of the reservoir at the end of the time period, setting the final storage capacity of the reservoir at the end of the time period to be 0, recording the dry and anhydrous state information of the reservoir, and ending the step S4; otherwise, correcting the final reservoir capacity at the end of the time interval as the initial reservoir capacity at the end of the time interval and the evaporation capacity of the reservoir, setting the evaporation capacity of the reservoir as 0, and entering the step S406;
s406, when the economic social water intake quantity can meet the negative shortage of the reservoir capacity, correcting the economic social water intake quantity to be the economic social water intake quantity plus the initial reservoir capacity at the end of the time period, setting the final reservoir capacity at the end of the time period to be 0, recording the economic social water intake shortage quantity to be the absolute value of the initial reservoir capacity at the end of the time period, recording the information of the dry and anhydrous state of the reservoir, and ending the step S4; otherwise, correcting the final reservoir capacity at the end of the time interval, namely the initial reservoir capacity at the end of the time interval and the water intake of the economic society, correcting the water intake shortage of the economic society, namely the water intake of the economic society, setting the water intake of the economic society to be 0, and recording the negative reservoir capacity state information of the reservoir.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910860489.2A CN110578317B (en) | 2019-09-11 | 2019-09-11 | Hydrological model reservoir discharge capacity simulation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910860489.2A CN110578317B (en) | 2019-09-11 | 2019-09-11 | Hydrological model reservoir discharge capacity simulation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110578317A true CN110578317A (en) | 2019-12-17 |
| CN110578317B CN110578317B (en) | 2021-01-15 |
Family
ID=68812132
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910860489.2A Active CN110578317B (en) | 2019-09-11 | 2019-09-11 | Hydrological model reservoir discharge capacity simulation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110578317B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111328671A (en) * | 2020-03-05 | 2020-06-26 | 红河哈尼族彝族自治州水利水电工程地质勘察咨询规划研究院 | Reservoir photovoltaic pumping irrigation control system and method for realizing automatic frequency adjustment |
| CN111597612A (en) * | 2020-05-08 | 2020-08-28 | 河海大学 | Structure recurrence period calculation method based on nonlinear flood storage and discharge model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09138708A (en) * | 1995-11-14 | 1997-05-27 | Hitachi Ltd | Method and device for controlling constant water level |
| CN105089003A (en) * | 2015-08-07 | 2015-11-25 | 中国电建集团成都勘测设计研究院有限公司 | Reservoir flood routing calculating method |
| CN105243458A (en) * | 2015-11-10 | 2016-01-13 | 河海大学 | Reservoir dispatching method based on multi-target shuffled frog leaping and differential algorithms |
| CN108269199A (en) * | 2017-12-25 | 2018-07-10 | 河海大学 | A kind of small hydro project spatial and temporal distributions formula of object-oriented goes out flow calculation methodologies |
-
2019
- 2019-09-11 CN CN201910860489.2A patent/CN110578317B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09138708A (en) * | 1995-11-14 | 1997-05-27 | Hitachi Ltd | Method and device for controlling constant water level |
| CN105089003A (en) * | 2015-08-07 | 2015-11-25 | 中国电建集团成都勘测设计研究院有限公司 | Reservoir flood routing calculating method |
| CN105243458A (en) * | 2015-11-10 | 2016-01-13 | 河海大学 | Reservoir dispatching method based on multi-target shuffled frog leaping and differential algorithms |
| CN108269199A (en) * | 2017-12-25 | 2018-07-10 | 河海大学 | A kind of small hydro project spatial and temporal distributions formula of object-oriented goes out flow calculation methodologies |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111328671A (en) * | 2020-03-05 | 2020-06-26 | 红河哈尼族彝族自治州水利水电工程地质勘察咨询规划研究院 | Reservoir photovoltaic pumping irrigation control system and method for realizing automatic frequency adjustment |
| CN111597612A (en) * | 2020-05-08 | 2020-08-28 | 河海大学 | Structure recurrence period calculation method based on nonlinear flood storage and discharge model |
| CN111597612B (en) * | 2020-05-08 | 2021-05-11 | 河海大学 | Structure recurrence period calculation method based on nonlinear flood storage and discharge model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110578317B (en) | 2021-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Si et al. | Revealing the water-energy-food nexus in the Upper Yellow River Basin through multi-objective optimization for reservoir system | |
| CN110851977B (en) | Water supply-power generation-ecological multi-target scheduling graph optimization method based on ecological flow | |
| Li et al. | Dynamic control of flood limited water level for reservoir operation by considering inflow uncertainty | |
| CN111080157B (en) | Method and system for scheduling phosphorus discharge amount of cascade hydropower station | |
| CN104376410B (en) | A kind of planing method of Distributed Generation in Distribution System | |
| CN105243438A (en) | Multi-year regulating storage reservoir optimal scheduling method considering runoff uncertainty | |
| US11620422B2 (en) | Water supply simulation method and tool for interlaced system of river system and canal system based on groundwater model | |
| CN101705671A (en) | Yellow River upstream cascade hydroelectric station operation design and optimized dispatching method as well as equipment | |
| CN105243502A (en) | Hydropower station scheduling risk assessment method and system based on runoff interval prediction | |
| CN110578317B (en) | Hydrological model reservoir discharge capacity simulation method | |
| CN111898253B (en) | Reservoir dispatching and downstream river ecological environment protection cooperation value evaluation method | |
| Diao et al. | Risk analysis of flood control operation mode with forecast information based on a combination of risk sources | |
| CN112036632B (en) | Optimal scheduling method based on cascade reservoir ecological power generation multi-target medium and long term random scheduling model | |
| CN115099477B (en) | Reservoir drought limit water level optimization and drought-resisting scheduling method | |
| CN113469528B (en) | Reservoir group multi-target flood control scheduling risk analysis method considering space-time correlation multi-dimensional uncertainty | |
| CN114819322B (en) | Forecasting method for flow of lake entering lake | |
| CN114358379A (en) | Optimization method, device and system for dry water limit level dispatching line of cascade reservoir | |
| CN101789053A (en) | Ecology scheduling technology for water-supply reservoir | |
| CN106909798A (en) | A kind of Daily rainfall multi-mode collection approach based on cumulative probability curve | |
| CN113592281A (en) | Power station performance analysis method, device, equipment and storage medium | |
| CN113078677A (en) | Energy consumption risk eliminating method considering uncertainty of renewable energy | |
| CN111328671A (en) | Reservoir photovoltaic pumping irrigation control system and method for realizing automatic frequency adjustment | |
| CN106557614B (en) | Flood frequency analysis method based on Halphen IB distribution | |
| CN114169590B (en) | Reservoir-in-storage runoff forecast correction method and system based on multi-scenario division | |
| CN115860421A (en) | Ecological flow dynamic calculation method and system adaptive to target guarantee rate |
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 |