METHOD FOR OPERATION OF HYDROPOWER RESERVOIR WITH A 2-PARAMETER ELEVATION RULE CURVE
20170234292 · 2017-08-17
Inventors
Cpc classification
F05B2220/706
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
H02K7/1823
ELECTRICITY
Y04S10/50
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
F03B15/14
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
Y02E10/20
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
G06Q10/04
PHYSICS
F05B2270/20
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
F03B3/02
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
International classification
F03B15/14
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
F03B13/08
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
Abstract
Disclosed herein is a novel method for operating a hydropower reservoir which is an improvement over the existing single-parameter (the current month) USACE Rule Curve approach, the improvement comprising the consideration of a second parameter, namely the water level of the reservoir at the beginning of the month, in the decision-making process for operation of the reservoir.
Claims
1. A method for optimizing the energy produced by a hydropowered reservoir over a time wherein: generation of a first Rule Curve for the reservoir in accordance with a USACE operational strategy; generation of a second Rule Curve for the reservoir which modifies said first Rule Curve by considering a second parameter, namely, the water level of the reservoir at the beginning of a month; operating said reservoir in accordance with said second Rule Curve.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0013]
[0014]
[0015]
[0016]
[0017]
[0018]
[0019]
[0020]
[0021]
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
[0028]
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
DETAILED DESCRIPTION OF INVENTION
[0040] Analytical Description of the Reservoir Power Plant System
[0041] Operation Pattern for Continuous Power/Firm Energy
[0042]
[0043] a) CPp and CPb must be lower or equal than the power output at the minimum power pool level.
[0044] b) Total time, Tb+Tp, equals 24 hr.
[0045] c) The relationship between the base and peak continuous power, fB=CPb/CPp, lies between zero and one. A value of zero means the power-plant is dedicated to generate peak load, whereas a value of one means the plant generates in base load only.
[0046] Mass Balance of Reservoir
[0047] The mass balance of the reservoir is represented by the following recursive equations, for any given period:
Vf=Vi+ΔV (“EQ 1”)
ΔV=Qinf ΔT−Qtp ΔTp−Qtb ΔTb+(Kr R−Ke Ev)Area−Qeco ΔT−Qloss ΔT (“EQ 2”) [0048] Subject to the inequality constrain:
Vmin≦Vf≦Vmax (“EQ 3”)
[0049] In which: [0050] Vf: Volume stored in reservoir at the end of the month. [0051] Vi: Volume stored in reservoir at the beginning of the month, known value. [0052] ΔV: Variation on volume stored during the period, without restrictions imposed by the reservoir levels (maximum and minimum). [0053] Qinf: Average Inflow into the reservoir during the current period. [0054] R: Total monthly rainfall over the reservoir area during the current period. [0055] Ev: Total monthly evaporation over the reservoir area during the current period. [0056] Qtp: Average monthly flow into turbines during peak time, Tp. [0057] Qtb: Average monthly flow into turbines during base time, Tb. [0058] Qeco: Average ecological flow and other demands during the period which do not pass through the power house. [0059] Qloss: Average water losses during the period. [0060] Vmin: Volume stored at the minimum reservoir power pool level. [0061] Vmax: Volume stored at the maximum reservoir power pool level. [0062] ΔT: Time conversion factor from the units of flow into total monthly volume, for full day. [0063] ΔTp: Time conversion factor from the units of flow into total monthly volume, for peak hours. [0064] ΔTb: Time conversion factor from the units of flow into total monthly volume, for base hours. [0065] Area: Average surface area of the reservoir during the period. [0066] Kr: Effective rainfall and conversion constant. [0067] Ke: Effective evaporation and conversion constant.
[0068] The relationship between level-area storage of the reservoir is represented by the level-area-storage curve. If the water surface level falls below the minimum power pool level, the reservoir cannot yield the required power and energy output, and the available water must be allocated to the diverse demands on a predefined priority basis. The final level is then set to the minimum power pool level. If the water level goes above the maximum power pool level, it generates a spill, which is then calculated. The final level is set to the maximum power pool level. This procedure is repeated on a month-by-month basis, for the duration of the hydrological time series available.
[0069] Power Output of Power Plant
[0070] The equations representing the power output of the power-plant are the following:
P=γQt H ef.sub.tef.sub.g (“EQ 4”)
E=Pp ΔTp+Pb ΔTb (“EQ 5”) [0071] Subject to the inequality constrains on the continuous power to generate firm energy:
Pp≧CPp (EQ 6)
and
Pb≧CPb (“EQ7”)
[0072] In which: [0073] P: Power output of the power-plant, it is indexed for the generation on peak time, Pp; and base time, Pb. [0074] CPp: Continuous power during peak time. [0075] CPb: Continuous power during base time. [0076] E: Total energy generated during period. [0077] γ: specific weight of water. [0078] Qt: water flow into the turbines, indexed for peak and base time, Qtp and Qtb. [0079] H: average monthly net head on the turbines, indexed for peak and base time, Hp and Hb. [0080] ef.sub.t: average efficiency of the turbine units, indexed for generation at peak and base time. [0081] ef.sub.g: average combined efficiency of the generator and transformer units, it is presumed constant throughout the simulation.
[0082] Net Head
[0083] The net head is a function of the reservoir elevation, the energy losses through the conveyance systems (free-surface and pressure flow) from the reservoir into the turbines, and the tailwater elevation rating, as follows:
H=Elev−HL−TW (“EQ 8”)
[0084] In which: [0085] H: Average net head on the turbines, indexed for peak and base power. [0086] Elev: Average surface level of reservoir during period. [0087] HL: Average energy losses during the period, due to conveyance from reservoir into the power-plant, indexed for peak and base power. It is a function of total flow into the power-plant and unit flow into the turbines. [0088] TW: Average tailwater level during the period, indexed for peak and base power, it is a function of total flow into the power plant, and must take into account any spillway discharges that affect the power house discharge channel levels.
[0089] Efficiency
[0090] The efficiency of the turbines is a function of the net head and unit flow on the turbines. The relationship between efficiency-head-flow is part of the data set and can be read from the efficiency hill-diagrams of typical turbines or of the installed turbines. The combined efficiency of the generator and transformer units is generally presumed constant.
[0091] Estimation of the Composite Monthly Reservoir Parameter, PT
[0092] The backward simulation algorithm used to define the continuous power and firm energy can be represented by the equations:
Vi=Vf−ΔV (“EQ 9”)
Pp=CPp (“EQ 10”)
Pb=CPb (“EQ 11”)
[0093] Subject to the inequality constrain:
Vmin≦Vi≦Vmax (“EQ 12”)
[0094] In which the variables are defined in equations (EQ 1) and (EQ 2), with the initial condition of reservoir at minimum power pool level at the end of the time series. The equations are solved recursively backwards, beginning at the end of the time series data, up to the initial period.
[0095] Once the continuous power and firm energy are identified by the backward simulation procedure, the following parameter is calculated for the time-series:
PT(i,j)=ΔV(i,j)Hef(i,j) (“EQ 13”)
[0096] In which: [0097] i, j: indexes for the year and month of the time series. [0098] PT(i,j): monthly value of the variation on stored effective potential energy in the reservoir. [0099] ΔV(i,j): monthly value of the variation on storage during the period. [0100] Hef(i,j): monthly average product of net head and efficiency of the power-plant (turbines, generator and transformers combined).
[0101] Estimation of the Variation in Stored Volume of Reservoir
[0102] The calculation of ΔV(i,j) involves two variants, defined as: unrestricted variation (A10) and restricted variation (A20) of stored volume, as follows:
[0103] a) A10—Unrestricted variation of stored volume in reservoir. It is calculated prior to the imposition of the inequality constrains on the stored volume from (EQ 2), restated below for ease of reference:
ΔV=Qinf ΔT−Qtp ΔTp−Qtb ΔTb+(Kr R−Ke Ev)Area−Qeco ΔT−Qloss ΔT
[0104] b) A20—Restricted variation of stored volume in reservoir. It is calculated after the imposition of the inequality constrains on the stored volume (EQ 12):
Vmin≦Vi≦Vmax (“EQ 12”)
ΔV=Vf−Vi (“EQ 14”)
[0105] If the inequality constrains are not binding, the two calculated values, unrestricted (A10) and restricted (A20), are equal.
[0106] Estimation of the Average Net Effective Head
[0107] The net effective head is the product of net head and efficiency of the power-plant. The net head is defined in equation (EQ 8) restated below, and the efficiency is a function of the net head and unit flow in the turbines.
HEF=H ef.sub.tef.sub.g (“EQ 15”)
H=Elev−HL−TW (“EQ 8”)
[0108] In which: [0109] HEF: Net effective head, indexed for peak and base power. [0110] ef.sub.t: Efficiency of the turbines, indexed for peak and base power. [0111] ef.sub.g: Combined efficiency of the generators and transformers, it is presumed constant throughout simulation.
[0112] The average net effective head is a composite measure over the time period, weighted on the duration of peak and base load. The equation is the following:
HEFavg=(HEFp Tp+HEFb Tb)/24 (“EQ 16”)
[0113] In which: [0114] HEFavg: weighted average net effective head over the period. [0115] HEFp: net effective head on peak power. [0116] HEFb: net effective head on base power. [0117] Tp: duration of peak power. [0118] Tb: duration of base power.
[0119] Data Processing of the Composite Monthly Reservoir Parameter PT to Obtain the Maximum Monthly Values
[0120] The control value for the estimation of the Rule Curve is the maximum envelope of the monthly values of the PT-Parameter, that is:
PTX(j)=max(PT(ij)),i=1,n (“EQ 17”)
[0121] In which: [0122] i,j: indexes for the year and month of the time series. [0123] PTX(j): maximum value of PT parameter for month j, over the entire series. [0124] n: duration of the time series, years.
[0125] This process is done for the two variants of the PT parameter, the unrestricted and restricted case. Each case gives an envelope to be used in the estimation of the 2-Parameter Rule Curve
[0126] Estimation of the End-of-Month Level for the 2-Parameter Elevation Rule Curve
[0127] The estimation of the end-of-month level proceeds after the extreme maximum monthly values of the PT parameter is calculated. The calculations follow the inverse procedure to that which gave origin to the PT parameter, and are described below:
[0128] For any given month or time period:
[0129] a) Define initial reservoir level: Any value ranging between the minimum and maximum power pool levels. To obtain the full range of end-of-month storage, the calculations must be repeated with initial level varying between the extremes, with the adequate variation to obtain a good definition of the curve.
[0130] b) Define the extreme values of the effective rainfall and evaporation: The parameter RE(i,j)=+Kr R(i,j)−Ke Ev(i,j) is the effective contribution of the rainfall and evaporation to the mass balance equation. Use the time series available to estimate REmin(j), the minimum monthly value of this parameter to use in the mass-balance equations.
[0131] c) Estimate the variation in storage to calculate the required end-of-month storage for all the initial reservoir levels required, considering the generation of continuous power and firm energy only. The equations used are, for any given month and initial reservoir level:
ΔV=PTX/Hef (“EQ 18”)
Vf=Vi+ΔV−REmin Area+Qeco ΔT+Qloss ΔT (“EQ 19”)
Vmin≦Vf≦Vmax (EQ 3)
[0132] In which [0133] ΔV: variation in storage from the beginning to end of month. [0134] PTX: monthly extreme value of the PT parameter for the month. [0135] Hef: average effective net head for the period, considering the generation of continuous power and firm energy only. [0136] Vi: reservoir storage at the beginning of the month, associated with the initial level. [0137] Vf: reservoir storage at the end of the month. [0138] REmin: minimum value of the effective contribution of rainfall and evaporation, for the month. [0139] Area: average surface area of the reservoir during the period.
[0140] d) Use the level-area-storage curve to determine the end-of-month level for the Rule Curve. Once the storage at the end of the month is calculated, use the level-area-storage curve to estimate the reservoir level that defines the Elevation Rule Curve.
[0141] Sample Calculations
[0142] Provided here is the complete calculation of the 2-Parameter Rule Curve, in line with the calculations proposed by the USACE Reference to develop their Rule Curve, are discussed. First, there are the data requirements, then the processing of the data with the backward simulation to generate the control parameters, and, finally, the estimation of the values of the 2-Parameter Rule Curve for each period and range of values within the admissible levels in the reservoir. The experimental site is the Tayucai reservoir, a site under study in the Upper Caroni Basin, in Venezuela.
[0143] Data Requirements
[0144] For optimal results, the data should be be the most complete and up-to-date available, covering the hydrological aspects, the characteristics of the reservoir and power-plant, and other water demands to and losses from the reservoir.
[0145] Hydrology:
[0146] The hydrological data includes the total monthly streamflow into the reservoir, as well as the data for direct monthly rainfall to and evaporation from the reservoir. The length of the record period should be as long as possible. The data is taken from the appropriate sources or estimated from the hydrological modeling of the basin.
[0147] Reservoir:
[0148] The data required for the reservoir are the following:
[0149] a) Level-area-storage curve: It is taken from the appropriate topographical maps, or from bathymetric studies of existing reservoirs. The data can be presented in tabular and graphical forms. The interpolation between data points can be made with a line segment, or with the aid of a curve-fitting procedure. The level-area-storage curve is presented in
[0150] b) Representative Levels: It is required to determine the maximum and minimum levels of the power pool, EL.sub.max and EL.sub.min, to define the usable storage for regulation purposes. These levels define the maximum and minimum storage, Vmax and Vmin, and the usable storage as the difference between the maximum and minimum storage values. For this sample, the minimum power pool level is 340 masl, the maximum power pool level is 380 masl. Minimum storage is 1,456 mill. m.sup.3, and maximum storage 9,419 mill. m.sup.3, the usable storage is 7,963 mill. m.sup.3, which is roughly 10% of the mean annual inflow to the reservoir.
[0151] Power Plant:
[0152] a) Operating pattern: The operating pattern is fully defined with the following parameters: [0153] Duration of peak time, Tp: 4 hour [0154] Duration of base time, Tb: 20 hour [0155] Ratio between base and peak continuous power, fB: Case 1: 1.00 (base load), Case 2: 0.40
[0156] b) Tailwater level: For reaction turbines, the tailwater level is the rating curve (flow vs. level) of the discharge channel For impulse turbines, the tailwater level is the fixed elevation of the discharge jets into the turbines. In this example, the rating curve of the discharge channel is shown in
[0157] c) Turbines: The turbines are defined by their number, type, design head and flow, the usable range of net head and flow, the ratio of maximum power over rated power, and by the efficiency hill-diagram. The turbines can be analyzed jointly as a single unit or separated, in which case the number of operating turbines varies with the net head, flow and power produced by the power-plant. In any case, the operation of the turbines should be the most efficient.
[0158] The power-plant is located at the foot of the dam, with one penstock for each turbine. In this example, design head is fixed at 80 m. The machines are Francis turbines, with the efficiency hill diagrams for a machine of specific speed 344 (metric). Total power output is about 3,800 MW, with several turbines (in the range of 10 to 15 machines). Hence, it is expected that the most efficient operation of the power-plant will provide that the unit flow of the machines will be greater than or nearly equal to the maximum efficiency flow for any given head and total power output. The combined efficiency of the generator and transformer units is presumed constant, equal to 0.96.
[0159] d) Head Losses: The power-plant is located at the foot of the dam. The head losses in the penstock are considered constant, equal to 1.5 m.
[0160] Other Demands
[0161] Other demands include ecological flow (minimum required flow, fisheries, navigation, other uses), and other water demands such as irrigation or urban water supply. In this case, the only demand is maintaining a minimum flow downstream from the power-plant equal to 10% of mean annual inflow, some 250 m3/s, to maintain ecological flow. If the total turbined flow at peak or base times is greater than the minimum required flow of 250 m3/s, no action is taken; however if the turbined flow is smaller than the minimum required, the difference between the minimum required flow and turbined flow is released into the river channel without going through the power plant.
[0162] Estimation of the Composite Parameter PT and Maximum Values PTX
[0163] Presented here is the backward simulation algorithm, and the calculation of the required parameters, along with the estimation of the USACE Rule Curve, for comparison and combination purposes, as will be seen below. The calculations are made for the two (2) cases indicated: Case 1 (base load) with fB equal to 1.0; and Case 2 (peak load) with fB equal to 0.4.
[0164] Backward Simulation:
[0165] The result of the backward simulation is presented in
[0166] Composite Parameter PT and Maximum Parameter Value PTX:
[0167]
[0168] Estimation of the 2-Parameter Rule Curves
[0169]
[0170] Estimation of the USACE Curve
[0171]
[0172] Combined Use of the 2-Parameter and USACE Rule Curves
[0173] The principal purpose of both the 2-Parameter and USACE Rule Curves is to preserve the Finn Energy yield of the reservoir. Both elevation rule curves represent a recommended end-of-month level based on the current time of year and, for the 2-Parameter Rule Curve, on the reservoir level at the beginning of the month.
[0174] One form of exploration of the complementary solution space is by making a linear combination of level/storage from the two Rule Curves with the appropriate weights, plus the extreme values (maximum and minimum) of the range. Thus is it formed the Combined 2-Parameter-USACE Rule Curve, which can be tested by simulation to select the most appropriate curve for a reservoir and power-plant operating in a wide range of rated capacity.
[0175] Extreme Values
[0176] This combined extreme curve is represented by any of the following equations:
Maximum Value: L.sub.C=max(L.sub.2-P,L.sub.USACE) a)
Minimum Value: L.sub.C=min(L.sub.2-P,L.sub.USACE) b)
[0177] In which: [0178] L.sub.C: Combined Rule Curve Level [0179] L.sub.2-P: Level of the 2-Parameter Rule Curve [0180] L.sub.USACE: Level of the USACE Rule Curve
[0181]
[0182] Linear Combinations
[0183] The linear combination of the level/storage is represented by any the following equations:
Linear Combination on Level: L.sub.C=w L.sub.2-P+(1−w)L.sub.USACE a)
Linear Combination on Storage: S.sub.C=w S.sub.2-P+(1−w)S.sub.USACE. b)
[0184] The Combined Rule Curve level, L.sub.C, is then interpolated from the elevation-storage curve.
[0185] In which: [0186] L.sub.C: Combined Rule Curve Level [0187] L.sub.2-P: 2-Parameter Rule Curve Level [0188] L.sub.USACE: USACE Rule Curve Level [0189] S.sub.C: Combined Rule Curve Storage [0190] S.sub.2-P: 2-Parameter Rule Curve Storage [0191] S.sub.USACE: USACE Rule Curve Storage [0192] w: weight assigned to the 2-parameter rule curve in the linear combination, lies within the normalized range 0 to 1. The weight can be the same throughout the year, or can vary with the month-seasons
[0193]
[0194]
[0195]
[0196] The simulated results show a smooth, non-linear variation between the energy generated by the 2-Parameter Rule Curves and the USACE Rule Curve as the weight assigned to the 2-Parameter Rule Curve varies from 0.90 down to 0.1. The maximum annual energy is generated by the use of the 2-Parameter Unrestricted A10 Rule Curve.
[0197] Linear Combination on Storage
[0198]
[0199] The simulated results show a smooth, non-linear variation between the energy generated by the 2-Parameter Rule Curves and the USACE Rule Curve as the weight assigned to the 2-Parameter Rule Curve varies from 0.90 down to 0.1. The maximum annual energy is generated by the use of the 2-Parameter Unrestricted A10 Rule Curve.
[0200] In conclusion, the evaluation of the combination of the 2-Parameter and USACE Rule Curves indicates that the 2-Parameter Unrestricted A10 Rule Curve generates the maximum mean annual energy output.
[0201] Use of the 2-Parameter Rule Curve
[0202] The use of the 2-Parameter Rule Curve can increase total energy output from a single reservoir over the energy produced by using the USACE Rule Curve, while keeping the Firm Energy Yield of the reservoir, over a wide range of rated capacities in the power plant. The example is again Tayucai reservoir, in which the rated capacity of the power plant varies between 2,000 MW and 6,000 MW, with normal power pool level set at 380 masl and minimum power pool level set at 340 masl.
[0203] Characteristics of the Power Plant
[0204]
[0205] Energy Produced
[0206] USACE Rule Curve
[0207]
[0208] 2-Parameter Unrestricted A10
[0209]
[0210] The continuous power and firm energy generation does not change, however the average energy output increases by 2,240 GW-h/year which represents an increase of 15% over the maximum average energy and of 26% over the maximum secondary energy generated by the USACE Rule Curve alone.
[0211] 2-Parameter Restricted A20
[0212]
[0213] The continuous power and firm energy generation does not change, however the average energy output increases by 1,740 GW-h/year which represents an increase of 11% over the maximum average energy and of 20% over the maximum secondary energy generated by the USACE Rule Curve alone.
[0214] Comparison
[0215]
[0216] As can be seen, the Unrestricted A10 Rule Curve generates more energy from the reservoir for every rated capacity installed in the Tayucai Power Plant, from a low 2,000 MW (Plant Capacity 80%) to a high of 6,000 MW (plant capacity 30%).
[0217] The mean annual energy generated with the 2-Parameter A10 and A20 Rule Curves reaches a maximum value and remains fairly constant near the maximum value generated notwithstanding the increase in rated capacity of the power plant.