The behavior of low-ternperature geothermal reservoirs under exploitation is simulated using analytical lurnped-parameter models. These models consider the effects of fluid production and reinjection, as well as natural recharge, on the pressures (or water levels) of low-temperature, liquid-dominated geothermal systems. The cornputed responses for constant production/injection flow rates are given in the form of analytical expressions. Variable flow rate cases are modeled, based on the Duharnel's principle. Reservoir parameters are obtained by applying a weighted nonlinear least-squares estimation technique in which measured field data are history matched to the corresponding model response. By using history-matched models, the future performance of the reservoir can be predicted for different production/injection scenarios in order to optimize the management of a given geothermal system.