The problem of multiple reservoirs storage control

For each reservoir i and at each time step t, system dynamics can be expressed as a set of water balance equations of the form:

Si(t + 1) = Si(t) + Ii(t) - Li(t) - Ri(t)

where Si(t) is the actual storage of the i reservoir, Ii(t) are the hydrological inflows, Li(t) are losses due to evaporation and leakage and Ri(t) is the controlled release rate.

The problem is to define either the unknown storages Si(t + 1) or, equivalently, the unknown releases Ri(t), in order to maximize the performance of the system. Note that the above magnitudes are under several physical constraints, due to the storage capacity of reservoirs and the discharge capacity of aqueducts, respectively.

Assuming a system consisting of N reservoirs and a simulation period of length T, the mathematical model has N x T parameters (control variables) and a similar large number of mathematical constraints.

ΝΑΜΑ Consulting Engineers and Planners S.A.

NTUA - Water Resources, Hdraulic & Maritime Engineering

Water Supply and Sewerage Company of Karditsa Municipality

AEPHORIKE of Dodecanesus Islands S.A.

Marathon Data Systems (MDS)