Multi-Objective Analysis of Bhakra Reservoir
Multi-Objective Analysis of Bhakra Reservoir
P S RAO* and S R A M A S E S H A N t
Indian Institute of Management, Bannerghatta Road, Bangalore 560076, India
t Civil Engineering Department, Indian Institute of Technology, Kanpur 208 016, India.
Abstract. The study outlines the methodology for working out a multiobjective reservoir
operation with conjunctive utilization of surfaceand groundwater. Irrigation and power have
been considered the two objectivesfor operation of the [Link] results of this study are
compared with those achieved with conventional operation in the past and show that there is
an increase in irrigation by 7 to 8 ~ due to conjunctive use of surface and groundwater. The
trade-offbetween irrigation and power has been developedand discussed. The transformation
surface between three objectives viz., power, irrigation and carry-over storage has been
developed and the merits of carry-over storage have been discussed.
Keywards. Reservoiroperation; multiobjectiveanalysis;conjunctive use of water.
1. Introduction
Irrigation and power are two of the most important purposes for which surface and
groundwater resources are used. The integrated operation of a complex system is
necessary for efficient management of limited resources to meet irrigation and power
demands. The importance of the conjunctive utilization of surface and groundwater in
an integrated framework has been realised (Irrigation Commission 1972) but has not
been put into practice perhaps due to institutional constraints and lack of coordination.
However, using a systems analysis procedure it seems possible to study the integrated
operation of a multiobjective reservoir system with conjunctive utilization of surface
and groundwater. A case study of the Bhakra reservoir, which has been in operation
since 1966 is described in this paper. The object of this study is to analyse the past
operation o f the reservoir within a multiobjective framework to evaluate the trade-offs
between irrigation and power demands in the past. The study also demonstrates the
methodology and procedures in multiobjective analysis within an integrated
framework.
2. System description
~'ro~,
fm-~.o~,/l/1 /
/ . Julun6ur
~ ~ ~\
\ \ ~'~"
% Power
.....
House No ](6angawa,)
car,ag-'-/{ / / ~ ~N, ~'---NangaL Hyde! Channel
9876 m cu m and a live storage capacity of 7814 m cu m above a dead storage level of
445.62 m. It covers an area of 168.35 km 2 when full. The total run-off at Bhakra for a
dependable year works out to 13,723 m cu m and that for a mean year to 16,441 m cu m.
A dependable year has been adopted by the planners as a hypothetical year in which the
mean discharge o f any 10-day period is equalled or exceeded by the corresponding 10-
day mean discharges in 67 % of the years for which data are considered.
The river flows are high for the four months from June to September due to
snowmelt and monsoon rains. Almost 47 % of the river flow of a dependable year
occurs during these four months. Water from Govindsagar (Bhakra reservoir) can be
passed through the turbines of two power houses, one on the right bank and the other
on the left bank at the foot of the Bhakra dam. Both the power houses have 5 turbines
each. Each of the generators in the right-bank power house has a maximum capacity of
120 MWwhen the head is 116 m or more, whereas each of the left-bank generators has a
maximum capacity o f 90MW when the head is 110m or more. These capacities fall
respectively to 70 MW and 53 MW when the head falls to about 80 m. Thus the power-
generating capacity varies from a maximum of 1050MW tO a minimum of 615 MW.
About 11 km downstream o f the Bhakra dam is the Nangal reservoir formed by the
28-95 m high Nangal dam. It serves as a head regulator for control of irrigation releases.
Part o f the water from Nangal is released to the Nangal hydel channel which is 64.48 km
long and has a carrying capacity o f 353"75 cumecs. The remaining water is reJeased to
the Sutlej. The Nangal hydel channel supplies water to two power houses in its path at
Ganguwal and Kotla with a total installed capacity of 154 MW. Water from the Nangal
hydel channel is then divided between the Bhakra main canal and the Sirhind canal. The
Bhakra main canal serves some o f the canals of the old western Jamuna canal system
which were merged with the Bhakra system and so receives some supplies from the
western Jamuna canal.
Downstream of Nangal, there are headworks at two places on the river Sutlej at
Ropar and Harike. At Ropar, water is diverted to the Bist Doab and Sirhind canals. The
Beas river joins the Sutlej at Harike. Water is diverted at Harike to the Rajasthan feeder
Study of Bhakra reservoir operation 181
and the Ferozepur feeder for the Eastern and Bikaner canals and the Sirhind feeder.
Water from all these canals is distributed through the network of the branch canals of
the irrigation distribution system. The details of the system are given in figure 1.
Punjab possesses large resources of groundwater which are generally suitable for
irrigation. Due to the presence of alluvium, a large-scale programme to install tubewells
has been undertaken. The State in 1974 had about 126,000 tubewells and another
75,000 were proposed under the fifth plan (1973-1978). An area of 0-6 to 0-8 m ha suffers
from waterlogging. The installation of tubewells would also solve the waterlogging
problem to some extent.
Daily
No. irrigation Total irrigation
of requirements requirements
Month Period days (in cumec days) (in cumec days)
and were modified only slightly since, although the irrigation demands have changed
due to the introduction of high-yielding varieties of crops which require more water.
The need to reassess irrigation requirements and revise the demands has been well
recognised.
Irrigation requirements are high from September to November owing to the water
required for the maturing of the kharif crops and the preparation and sowing of the
rabi crops. The requirements are low from December to April, except in the latter half
of February and in March when water is required for the maturing of the rabi crops.
The rabi crop requirements are smaller than the kharif crop requirements. In May and
June, the requirements are again high due to the preparation and sowing of the kharif
crops. During the monsoon season also the requirements are quite high as the areas
where such water is to be utilized have low rainfall.
The Bhakra system has been studied and reported by economists (Raj 1960; Ansari
1968; Minhas et a11972; Reidinger 1974) and by engineers (Singh 1964; Mehndiratta &
Hoon 1973). Raj concerned himself with some economic aspects of the Bhakra Nangal
project and a preliminary analysis of selected investment criteria. Ansari studied the
economics of irrigation rates and examined the basic arguments for incurring financial
deficits on irrigation supplies. Reidinger studied the institutional rationing of canal
water in an area in Hissar district which is served by one of the Bhakra canals.
Minhas et al (1972) made extensive simulation studies of the operations of the
Bhakra system (i) to determine for the depletion period the efficient combinations of
firm power with certainty and irrigation supplied with a given confidence level, (ii) to
evaluate the increase in availability of both irrigation and firm power due to conjunctive
utilization of surface and groundwater and to compare the economics of this scheme
with that of a thermal back-up. They observed that "in the absence of a satisfactory
economic measure of the relative worth of irrigation and power, one cannot determine
an ideal operating policy". But the important question of devising such a satisfactory
economic measure has been left unanswered in the study. As such the trade-off
possibilities delineated in a probabilistic sense cannot be used by decision-makers in
choosing the point of operation of the reservoir at the end of the filling period.
Singh (1964) put forward the idea of using seasonally available secondary power for
lifting groundwater by tubewells in the general framework of conjunctive utilization of
surface and ground waters. Minhas et al (1972) noted that the tubewells, when operated
solely to augment irrigation supplies, can utilize power which would otherwise go waste
and which therefore is 'free' in the economic sense of the term (since the economic value
of goods is based on their alternative uses). Operated solely to increase the level of firm
power output, tubewells can convert useless secondary power to valuable firm power.
Singh also pointed out the role of groundwater development as an antiwaterlogging
measure.
Mehndiratta & Hoon (1973) report the methodology adopted by B}aBto regulate
supplies from the Bhakra reservoir from 1966 to 1972. This report greatly helped in
understanding the system and serves as the source material for the present study. Water
power studies of BMB,given in tables 2 and 3, reveal that according to the planning
studies for the operation of the reservoir in 1969-70 and 1970-71, considerable
secondary or dump power was generated during the early part of the depletion period.
From the studies reviewed, it is clear that ti) the sudden peak in the irrigation
demands during the rabi sowing season, which is the early part of the depletion period
of the reservoir, necessitates the release of water with high potential for power
generation for irrigation; (ii) consumers of power who get used to the higher level of
firm power supplied during this period find it difficult to reduce their consumption of
OO
ta.
Gains or Total
losses requirements
Inflows between Contributions of canals Available
in river Bhakra from Sirhind to be fed Actual Actual Reservoir average
Sutlej & Ropar feeder from Bhakra release RF elevations power
Month Period (cumecs) (cumecs) (cumecs) (cumecs) (cumecs) (~0) (m) (MgO
1 2 3 4 5 6 7 8 9 10
OO
Table 3. Water power study for the year 1970-71
oo
o~
Gains or Total
Inflows losses Contributions requirement
in river between from of canal to Available
Sutlej Bhakra Sirhind be fed from Actual Reservoir average
at Bhakra & Ropar feeder Bhakra releases Actual elevation, power
Month Period (cumecs) (cumecs) (cumecs) (cumees) (cumees) RF (%) (m) (MW)
1 2 3 4 5 6 7 8 9 10 t,.,
October 1970 21-31 191.19 .. 96-45 654.64 410-60 77.5 492.86 439"00
November 1-10 170-89 .. 65'07 654"64 42ff75 74'9 491.34 447.8
11-20 152.12 .. 58.96 654.64 420'75 74"0 489"20 441.5
21-30 137.65 .. 55.87 654.64 420-75 73'5 487.07 434"5
December 1-10 130-06 + 2.83 53.58 611'75 396-44 74'0 484.94 399.9
11- 15 111.25 + 2.83 55.64 363'67 311'96 101-8 484.02 311"9
16-20 111.25 + 2'83 47.68 345.81 313.88 105"3 483.41 311'9
21-31 114.09 + 2.83 47.68 345.81 317-17 106.3 481"58 311.9
January 1971 I-I0 112.90 +5.66 29.08 345.81 322.00 103"I 479.76 311"9
11-15 114.66 + 5'66 41.40 345.81 325.70 107-8 478.84 311'9
16-20 114.66 + 5.66 39.04 363.67 328.27 102'5 477.92 311'9
21-31 106-98 + 5.66 53-32 363.67 332'72 107.7 475.79 311'8
February 1-10 105.37 + 11"33 63.94 390-46 339' 17 106" 1 473.35 311"8
11-20 105.90 + 11.33 48.99 545.52 345.92 74-4 474.22 311.9
21-28/29 109.21 + 11.33 70.79 545.52 352"45 79.6 468.09 311"9
March 1-10 118.76 + 5'66 84.95 545.52 359"54 82.5 466.65 311-9
11-20 116.19 + 5.66 84.95 545'52 368"20 84" 1 463'91 311.9
21-31 133.74 + 5.66 84.95 545.52 378-25 85'9 460.55 31 I-9
April 1-10 150.56 -8, 50 53.83 345.81 389'I 1 125.6 458'I I 311'9
11 -l 5 149.12 ~ 8.50 63.17 345.81 307-67 131 '6 456.59 311.9
16-20 149.12 -8.50 63'17 418'89 404.18 110-2 455.07 311'8
21-30 186.79 - 8"50 72.21 418.89 414"44 114"1 452.02 311 "9
May 1-10 243.67 - 16.99 96.45 703.50 427' 18 72"0 449"58 311 "9
11-20 274'56 -16"99 87.75 703'50 438'34 72'3 447'14 311'8
21-31 400-26 -,I 6-99 108-82 703' 50 446"30 76"4 446'23 " 311 "9
power subsequently; (iii) it has not been possible for Br~Bto satisfactorily coordinate
the releases for irrigation and power; and (iv) conjunctive utilization of surface and
groundwaters may lead to increased firm power level and reliability of irrigation.
4. Conjunctive utilization
where i = index for the subperiod i; subscript i denotes the value of the variable in
subperiod i, N = number of subperiods, and X = release from the reservoir in
cumec days.
4.3 Constraints
4.3a Power constraint: The reservoir release in each subperiod should generate the
firm power supplied in that subperiod.
X i >~Pi/Ki (for all i), (2)
188 P S Rao and S Ramaseshan
INFLOWS
@ BHAKRA
RESERVOIR
SECONDARY I I~RIG~.TION
POWER DEMAND
'r
I
l
I
i 1
IGROUNOWATER~
L ...... J
Figure 2. Conjunctiveuse model for the Bhakra reservoir.
Turbine efficiencies of the order of 0.85 and pump efficiencies of 0-7 are common. The
product of the turbine and pump efficiencies is assumed to be 0-6. Tail race elevation
(357 m) is taken from Singh (1965). The power conversion factors Ki are computed
from the data in tables 2 and 3. They are assumed constant during the subperiods. As
the length of the subperiod is short (8 to 10 days), the variation in the reservoir elevation
is small and the head is high, the effect of this approximation is negligible.
The depth to groundwater is assumed to be 15 m in October and 17 m in November.
Minhas et al (1972) mentioned that in Haryana and Punjab groundwater level rises to
15 to 25 m (50 to 80 ft) below ground surface even for deep tubewells. In many areas
groundwater is available much nearer the surface at depths of 3 to 4 m especially after
the rainy season in October and November. The assumption of 15 to 17 m depth to
groundwater is conservative.
The firm power level and reservoir factor adopted by aMa for the early depletion
period up to November 30 are respectively 552 MWand 0'85 in 1969-70 and 400 Mw and
0"75 in 1970-71. These values are used in this study.
5. Muitiobjective analysis
Releases
Secondary suggested
power in by Model Available Power for Water saved in reservoir
BMB Study Power Groundwater
Month Period (MW) cumecs cumec days (MW) (MW) cumecs cumec days
Releases
Secondary suggested Power Water saved
power in by model Available for in reservoir
BMB study power groundwater
Month Period (MW) cumecs cume¢ days (MW) (MW) cumecs cumec days
There has been sharp competition for and conflict in the use of the Bhakra reservoir
storage for irrigation and power (Mehndiratta & Hoon 1973; Reidinger 1974). In the
planning stage Bhakra was considered more of an irrigation project and the power
benefits were secondary. A firm power of 282 MW was expected from Bhakra power
houses in a dry year. With the rapid increase in power load in the region, the position
has completely changed. The power demand overrides other interests of the project
during most of the months and it has not been possible to satisfactorily coordinate the
releases. It is therefore worthwhile to consider irrigation and power as the two
objectives to operate the Bhakra system and to study the complementarity and conflict
between the two for the six years for which data are available.
5.3 Constraints
5.3a Irrigation constraint: The reservoir release in each subperiod should meet the
irrigation requirements to be supplied from the reservoir,
Xi >1Ai, with i -- 1, 2 . . . . . N, (6)
where N - - n u m b e r of subperiods during the depletion phase of the reservoir,
X = reservoir release in cumec days, and A = irrigation requirements to be supplied
from the reservoir in cumec days.
5.3b Power constraint: The reservoir release should generate the committed firm
power in each subperiod.
Xi/>--
P'flini with i-- 1,2 . . . . . N, (7)
gl
where fl~ = firm power coefficient in subperiod i, n~ = number of days in subperiod
i, and K = power conversion factor(Mw/cu mec).
5.3c Continuity constraint: The total volume of water released from the reservoir
during the depletion period should not exceed the total utilizable storage in the
reservoir at the beginning of the depletion period and the inflows in the depletion
period.
N N
x, v+ E 1. (8)
i=1 i=1
where V = volume of reservoir storage utilized during depletion period in cumec days,
and I = river inflow into the reservoir.
Table 6a. Availability of water in filling period
1966- 67 11114 913 382 4367 439.83 4-6-66 509.26 14-9-66 6463
1967-68 11397 1172 308 4305 443-80 3-6-67 512.23 21-9-67 6944 ta.
1968~9 9819 950 271 4317 454-74 31-5-68 506.76 12-9-68 6056
1969-70 12002 962 74 5230 451.03 27-5-69 513.30 21-9-69 7142
1970-71 7722 1011 247 3256 444.92 25-6-70 495-98 3-10-70 4527
1971-72 11077 999 358 3910 443-23 28-5-71 512.76 16-9-71 7043
Starling Minimum
reservoir reservoir
level in Pattern of level in Reservoir
Study No. meters inflow meters factor % Power
Starting Minimum
reservoir Pattern reservoir
Study level in of level in Reservoir
No. metres inflows metres factor % Power
(table 9). The fll values corresponding to the same were used in this study for the
respective years. The fli were found to be respectively 1"26 in 1969-70 and 1-282 in
1970-71. For 1966 to 1969 for which no data are available, a distribution similar to
1969-70 and 1970-71 with a fli value o f 1.27 is used.
5.4d Power conversion factors: The power conversion factors K i for 1969-70 and
1970-71 are deduced from the data in tables 2 and 3 respectively and the relationship
between K~ and reservoir elevation is derived. For the other four years, the
approximate reservoir elevations for each of the subperiods are estimated from the
anticipated and actual depletion curves in figures 3 (a to f) and the power conversion
factors corresponding to these elevations are used in this study. These depletion curves
generally correspond to RE values in the ranges of 0.75 to 0.85. For other RF values, the
depletion curves would be somewhat different and the reservoir elevations and hence
the power conversion factors would also be different. However, for simplicity, it is
assumed that the conversion factors are independent of the R~ value and the same
conversion factors as in the water power studies of BMa are used in each of the
subperiods. F r o m the solution thus obtained, it is possible to determine the reservoir
elevations implied and hence determine a new set o f values of K~ and repeat the
procedure till convergence is obtained. However, the irrigation requirement is
constraining generally only from 23 September to 10 December when the reservoir
elevation is very high and the depletion rate is low. Such refinement is not considered
here.
5.4e Sutlej river inflows: This study assumes the dependable year inflows for
planning as in the case of water power studies of aua. They are obtained from
Mehndiratta & H o o n (1973) and Lamba & Prem (1975) for all years except for 1971-72
for which there is a discrepancy between the two reported values. The dependable year
inflow for that year is corrected by interpolation. These are shown in table 10. For
actual operation, actual river inflows (table 6b) are used.
Study of Bhakra reservoir operation 195
Note:
(1) Study Nos. I to VIII are BMa Water Power Studies for 1971-72 (Mehndiratta & Hoon 1973).
(2) Study No. XI is one of actual operation.
5.4f Stored volume: For planning studies, live storage available above dead storage
level (DSL)is used in the different studies. In case carryover was planned for, carryover
storage was estimated from volume elevation curves and deducted from live storage
above DSLto give utilizable storage with carryover. It is presumed that the losses due to
evaporation and infiltration from the reservoir are considered by BMB in the
specification of demands or in specifying inflows. The total volume of water for
196 P S Rao and S Ramaseshan
utilization during the years is assumed to bc the sum of dependable year inflows and the
utilizablc storage from the rcservoir as in thc water powcr studics of [Link] studics of
actual operation the volume of water actually utilized during the dcplction period is
used. The details of the data are given in table 10.
i i i I ~ i i i J
510.54 -.Max. R.L.509.~6 on 14.9,~6" 510.54 .L. 512.23 on 2.9.67
495.3o 495.30 [ ~
z
480.06 . . . . \~
.J
Anticipated ~ \
depletion curve \ \
~ e a l Oeplefion Actual deplebon
~ 464.82
uJ curve 464.82 curve \\~
\\\
449.58
449.58 \
M'In.R.L.4S4.74 on 31 56~
441.95 441.96 - _. I r
PERIOO PERIOD iocTINOvIOEC~JM~[FEBIMAF~Iz,Pr~I~AY
a. 1966 - 57 b. 1967- 68
i i i i f
510.5~ x. R t 506.76 on 12 9.68 510 54 on 21.9.69-
v3
LO
E 495.3o 4~30
w
z
J
~ 480.06
Anticipated ..~ . Antiopated \
~ 464.B2
(~tplehon curve ~\
46~82
\
cur ve
uJ
L
nr
449 58 449.58
vim. RL 451.03 on 27 5.69 \ Mtn. R L 444.92 on 2
C- 1968 - 69 d. 1969- 70
Study of Bhakra reJervoir operation 197
5105a 9 71.
i l i l I
495.30 5 0 n 31070
~r 4 ~ 30
W
z
w 480.06
,,-i, ~8o 06
>
Z
W
g
~ 46~.82 Anbc~patcd
deplebon curve ~, 464 82 curve
Actual depletion ~\
449.58
441.96 " I I I t . . . . J-
441.96 I I ~ I I _
PER~D PERK)D 3C" i NC' 'IDFC , AN~ FEg!MaP'j a ,F~IMAV ]JL~,
¢ 1970-77 f 197]-72
oo
Stored Depletion Total supply Actual
supply period inflows planned for Stored depletion Total supply
above (dependable utilization supply period actually utilized
D.S.L. year) utilized inflows
Year (m cu m) (m cu m) (m cu m) (cumec days) (m cu m) (m cu m) (m cu m) (cume¢ days)
Case 1 Case 2 Case 1 Case 2 Case 1 Case 2 Case 1 Case 2 Case 1 Case 2 Case 1 Case 2
0"60 430-8 424.25 445.3 465.0 413"6 401"2 461.9 514-9 320"9 310-1 447.7 521.3
0"65 430-8 424-25 445'3 465"0 413.6 401"2 461"9 514-9 320"9 310.1 447-7 521"3
0"70 430"8 424-0 445"3 465.0 413"6 401.2 461-9 , 514'9 318"8 305"8 447.7 521"3
0'75 428"1 420"7 445.3 465"0 413.6 401~2 461-9 514"9 308'2 288.4 ~147.7 521"3
0"80 420'9 411-8 442"5 464.4 412"5 397-8 461"0 514"9 279-4 250"5 446"3 521-3
0"85 404.7 390"3 432'5 459.6 404.0 386-0 454"4 514"9 227-6 N.F. 437"6 521"3
0.90 353"0 339"0 414"5 442.6 383"5 350-8 439'8 512.0 N.F. N.F. 422"5 520'5
0"95 263"1 N.F. 367"8 397'9 334"2 299"3 409.0 503"7 N.F. N.F. 402.4 513-7
1.00 N.F. N.F. N.F. 339"7 N.F. N.F. 354"1 488"0 N.F. N.F. 352-0 499"9
N.F.: Not feasible; Case 1: Study based on stored supply planned for utilisation and dependable year inflows during the depletion period;
Case 2: Study based on utilized stored supply and the actual inflows during the depletion period.
Study of Bhakra reservoir operation 199
1.0
0.9
n,.
:\.\ \\ ';
~0.8 ,~-67--. t ~L\ ~1 I~ ,oo7-~ I ',
,.<
~ 0.7
n, , il Irl~ i,
, il II I! i i
0.6
Planned Operation
Actuat Operation
0.5 _1
200 300 400 500
FiRM POWER (MW)
IRRIGATION
i
o2 ~b
i," I
' C
r ;
,/ ls3 POWER
t31k / . . . .
D Oec~sio~ Polrff.
for Operation
Reservoir level at the beginning of 512.06 512.06 512.06 512.06 512.76 512.76 512.76
depletion period (m)
Reservoir level at the end of 445-92 451.41 456.90 452.32 452.32 445.62 445'62
depletion period (m)
Stored supply utilized during 6,871.61 6,521'69 6,162.52 6,462.68 6,634 7,043 7,043
depletion period (m cu m)
Dependable year 4,059 4,059 4,059 4,059 4,059 4,059 5,205*
inflows (m cu m)
Total supply utilized (cumec days) 10,931 10,581 10,222 10,522 10,693 11,102 12,248
Total supply utilizable (mcum) 126,475 122,476 118,270 121,743 123,726 128,458 141,758"*
Carryover storage (mourn) 19.43 349.92 709.08 409.93 408-93 0 0
Carryover storage as percent of 0.3 4.94 10.32 5.95 5.95 0 0
total storage of 6871.61 m c u m
Table 13. Maximum firm power for different retervoir factors and carryover storages
(water power studies for 1971-72)
Carryover
storage
(m cu m) 19.4 349.9 709.1 408.9 408.9 0 0
* Actual operation
in the below average years 1966-67 and 19684i9; slightly larger in the above average year
1967-68 and very much larger in the wet years 1969-70 and 1971-72. It shows that
using dependable year inflows for planning every year irrespective of how wet the year
has been and how high the reservoir has filled up led to very conservative planning in
wet years and the opposite of it in a dry year. It is not sound practice. It seems possible to
delineate an average year with a transformation curve between the above average and
RESERV(~R FACTOR
//
0.9
0.8 /
I0
O.O
~- 0.8
u_
o
0.7
0.5
0 0 ,5 I.o
CARRYO',/ER STORAGE (b ¢u m )
• RF 0 . 8 5
v
o~ A RF 0 905
o 4s0
==
U_
Z,00
0.5 ,b
CARRYOVER STORAGE ( b ¢u m )
below average years for which the assumption o f dependable year inflows in depletion
period would be correct. Forecasting models for the stochastic inflows (Krishnasami
1976) are required to predict depletion period inflows for use in planning.
The points o f operation D chosen in the planning studies in 1969-70 and 1971-72
(figure 4) lie inside the respective transformation curves indicating that the planned
operation is inferior and nonoptimal. The point D could be moved on to the curve at Dt or
D 2 and the operation would be Pareto-superior to that at D. A higher firm power level
would be attained for the same Rr o f 0L85 at D1 while a higher Rr would be obtained for
the same firm power level at D 2. RF planned was never higher than 0.85 and so it seems
reasonable to suppose that the points DI would be the decision points o f BMB for
operation in 1969-70 and 1971-72.
The RF planned for was 0.75 in the dry year 1970-71. The decision point chosen in
planning is not feasible with respect to the actual transformation curve. The firm power
Study of Bhakra reservoir operation 203
level planned for was stuck to and the adverse impact of the shortage must have been
borne by irrigation only. So the actual RF provided would have been lower than 0-75.
Irrigation supplies at RFof I were feasible in 3 out of 6 years and at RFof 0"9 tO 0"95 in
2 out of 6 years. In the dry year an RF of only 0-8 was feasible.
Firm power levels planned for varied from 312 MWin a dry year to 438 Mw in a wet
year. In actual operation it seemed that 510 to 520 raw of power was generated in the wet
years of 1969-70 and 1971-72 and 280 to 300Mw in the dry year of i970-71. The high
order of fluctuations in firm power from year to year can be reduced by carryover
storage from wet to dry years. For example, a carryover of about 1500mcu m storage
from 1969-70 to 1970-71 would have facilitated operation in both years around the
level of an average year (1968-69) with firm power of about 400 to 420 raw. The
recommendation of Minhas et al (1972) that a high dead storage level should be
maintained for high levels of firm power (which are possible in years of good inflows)
implies, in a way, the concept of carryover storage. Carryover storage was planned for
in 1971-72 possibly because of the experience of the previous two years, but in actual
operation this sound policy was not carried out.
Irrigation and power are almost wholly complementary up to an RFlevel of about 0.9
in years of good inflows but the complementarity could not be exploited because of the
uncertainty of depletion period inflows. This again points to using appropriate
streamflow forecasting models for planning reservoir operation.
The trade-offs betw¢cn irrigation and firm power in the planned and actual operation
for the six years are given in table 14. The trade-offs implied in planned operation vary
from 200 MW/RI: in average years to 360 MW/Rv in a dry year. Since the results
indicate a higher trade-off for irrigation in a dry year, the value system adopted in
planning is consistent.
In 1969-70 (figure 4), the trade-offs of 200 MW/RF at Dt and 500 MW/Rr at D 2
provide the bounds for the marginal rates of transformation implied in the choice of
decision point D for operation. The BraBhas used a maximum R~ of 0"85 for operation
(tables 7, 8 and 9). Assuming, therefore, an Rv of 0.85 corresponding to D~ the trade-off
implied is 200 MW/Rv. The trade-off on the transformation curve of actual operation
at the level of D~ corresponding to RF of 0-85 is 0-0 MW/Rr. This implies that, in actual
operation, there is complete complementarity between irrigation and firm power up to
Rr of 0"85 and that the irrigation supply could have been increased upt0 Rr levels of 0.9
or even upto 1 with a small loss of firm power. Even at an RVof 1"0, 50 MW more of firm
power than that planned for would have been produced (table 11).
In the dry year 1970-71 (figure 4) the trade-offbetween irrigation and power implied
by the decision is 360 MW/xF. The planned operation at point D is not feasible for
actual inflows. If an Rv level of 0"75 was maintained, a firm power level of 288-4 MW
would have been reached on the transformation curve for actual operation and it would
have meant a trade-off of 490 MW/Rv at that point. This is larger than the trade-off in
an average year. During the early part of the depletion period from 3 October 1970 to
about 15 December 1970 the operation was as per the plan, i.e. at point D with RVof 0"75
and firm power level of 400 MW (figure 3e). The adverse effect of the shortage seems to
have been borne by the latter part of the depletion period from 15 December 1970 to 31
May 1971. The average power generation for the year 1970-71 was 335 MW. It can be
deduced from these that the Rr level achieved was certainly less than 0-75 and must have
been somewhere between 0.65 and 0.7 with a trade-off of 160 MW/RF or less. This, in
turn, implies a higher value for power than irrigation compared to an average year. This
indicates not only a contradiction between planning and operation but also an
inconsistency in economic evaluation. If, as pointed out earlier, the reservoir was not
depleted as it was in 1969-70 but some storage was carried over to 1970-71, it would
have greatly cushioned the adverse impact due to'poor inflows in 1970-71 and enabled
the achievement of a higher RF for irrigation.
Trade-off (MW/XV)
Irrigation
level Planned Actual
Year (RF) operation operation
5.7b Transformation surface for 1971-72: The concept of the transformation surface
amons the three objectives has been illustrated in figure 5. ABC constitutes the surface
bounded by the three coordinate planes. D is the decision point for operation. If it is on
the surface, it is Pareto-optimal. Curves a2a3, b263 and c2c3 pass through the decision
point and lie respectively in planes ala2a 3 parallel to the coordinate plane BOC, blb2b 3
parallel to OAC and c~c2c3 parallel to OAB. The trade-off implied by decision point D
between any two of the objectives can be derived as the negative of the slope of the curve
in the plane parallel to the coordinate plane defined by the same objectives.
For 1971-72, water power studies for planning the reservoir operation were
conducted by SMBassuming different levels of irrigation and firm power supplies and
different carryover storages. They are listed in table 9. Four of the studies, viz., I, II, III
and VIII assumed an irrigation level at RF 0-85 and dependable year inflows during
depletion period but they differ in that they assume different carryover storages. Details
of these studies are given in table 12.
The maximum firm power levels for different reservoir factors and carryover storages
are given in table 13. The transformation surface among the three objectives, irrigation,
firm power and carryover storage for the four planning studies is shown in figure 6. The
eMB decided to operate the reservoir according to study No. VIII with a carryover
storage of 408.93 mcum, a firm power level of 420 MW and RF of 0"85. A firm power
level of 437-6 MW was possible for an RFof 0"85and carryover storage of 408-93 m cum
(table 13, figure 4). Hence the decision point D (figure 4) is inside the transformation
curve and though feasible is not Pareto-optimal. For the same carryover storage, any
operation between D 1 (RFof 0"85 and firm power of 437.6 MW) and Dz (RFof 0-905 and
firm power of 420 MW) will be better than the point D chosen and will be Pareto-
superior. The trade-offs implied by choice of DI or D2 as decision points can be
obtained from the transformation curves in figures 4, 7 and 8 and are respectively 230 to
350 MW/RF between irrigation and power, 4824 to 5976 m cu m/Rr between carryover
storage and irrigation, and 60-6 to 67.6 MW/b cu m (billion cubic metres) between firm
power and carryover storage. As the SVlevel of 0.85 was always adhered to as an upper
limit, it appears that D~ would have been the choice. Hence the trade-offs implied by DI
are relevant indicators of decision makers' values.
The planning studies of BMS (table 9) assumed that the reservoir would reach an
elevation of 512-06m on 16 September 1971. But the reservoir actually filled to an
elevation of 512.76 m on that date with 171 mcum more storage than planned for in
these studies. Study Nos. IX and X (table 9) are revisions to account for this increase.
Study No. IX assumes carryover storage of 408.9 m cure as in study No. VI and
corresponds to planned operation. Study No. X assumes no carryover storage as in
Study No. XI and corresponds to actual operation. They are reported in tables 12 and
13. The transformation curve for actual operation (figure 4) indicates that irrigation
could have been supplied at an RFof 0-9 instead of 0.85 without any loss of firm power
and this too at a power level of 521.30 MWas against 420 Mw planned for. An increase in
ar to 1 would have resulted in a loss of only about 20 MW.
If the reservoir was operated according to the decision on 16 September, 1971 the
carryover storage would have been much more than what was planned for because the
excess of actual inflows over the assumed dependable year inflows, would have been
stored, and carried over to the subsequent year. The reservoir elevation at the end of the
depletion period would have been about 8 m higher than the planned elevation of
452-32 m. In actual operation the reservoir was drawn down to about the dead storage
206 P S Rao and S Ramaseshan
6. Summary of results
The analysis of the results o f the study of the Bhakra reservoir operation from 1966 to
1972 indicates that (i) the depletion period inflows are correlated to the filling period
inflows and planning should be based on appropriate estimates o f depletion period
inflows; (ii) an RF o f 0-85 or above can be attained in an average year and an gF o f 1 can
be obtained at a loss of firm power of about 80 raw in an average year; (iii) planning and
operation should be based on consistent estimates of economic parameters; and
(iv) carryover storage will serve to reduce year to year fluctuations in firm power and
irrigation levels.
References
Ansari N 1968 Economics of irrigation rates--a study in Punjab and U ttar Pradesh (Bombay: Asia Publishing
House)
Irrigation Commission, Reports of the Irrigation Commission 1972 Vols. I, II, lII. Ministry of Irrigation and
Power, New Delhi
Krishnasami M 1976 Decoupled multivariate time series models for multisite streamflows Ph.D. thesis, IIT,
Kanpur
Lamba S S, Prem K S 1975 Integrated development of rivers Sutlej, Beas, and Ravi for optimum utility of
water in north-western areas of India, Presented at the IWRA Second World Congress on Water
Resources, New Delhi, Reprint Vol. 3, p. 79
Mehndiratta K R, Hoon R N 1973 J. Central Board oflrrigation and Power India 30:349
Minhas B S, Parikh K S, Srinivasan T N 1972 Sankhya, Indian J. Stat. B34:105
Raj K N 1960 Some economic aspects of the Bhakra Nangal Project-a preliminary analysis in terms of selected
investment criteria
Reidinger R B 1974 Econ. Dev. Cult. Chanqe 23:79
Singh H 1964 J. Central Board of lrrigation and Power, India 21:471
Singh J 1965 in Proceedings of the sixth regional conference on water resources development in Asia and thefar
East, Water Resources Series, United Nations, p. 281-295
Thomas H A Jr, Revelle R 1966 Manage. Sci. 12:13-296