This paper describes the theory and algorithm for a practical, analytical model for simulation of geothermal reservoirs. The model can be applied to any type of geothermal reservoir: all liquid, all steam, or two-phase. The model assumes radial geometry with three distinct zones: (1) an innermost, circular production zone; (2) an intermediate, concentric zone subjected to fluid flow and heat transfer but no production or injection; and (3) an outermost, radially infinite (or finite) aquifer zone. The reservoir is assumed to be producing from the innermost zone, as dictated by the capacity of a power plant. The innermost zone is treated as a tank from the viewpoint of mass and energy production, using a lumped-parameter model, but one in which the pressure distribution is calculated by an analytical formulation of pressure behavior of off-centered wells enclosed in a constant pressure boundary circle. The surrounding intermediate zone from which fluids migrate into the production zone is assumed to have a radial temperature profile, T(r), depending upon the reservoir. This profile can be either continuous or a step profile. The injection occurs in the outermost zone and induces a flow of fluid from the outermost zone through the middle zone into the central tank. There is no constraint on the relative size of various zones. Injection in the central production zone (such as in a five-spot development) is possible. If there is no injection, natural transient water influx from the aquifer is considered.
For a given step of depletion in the tank model, mass and energy produced are calculated and the pressure distribution in the central part of the reservoir is computed. The mass of fluid injected and/or natural water influx from the aquifer zone is estimated, and corresponding fluid flow from the intermediate zone into the central zone is calculated. The new temperature profile in the intermediate zone is then computed, and the average enthalpy of the influx water into the central zone is then estimated. The average temperature and pressure in the central zone are then calculated and a new step of depletion is initiated using the lumped-parameter model. If either constant enthalpy or constant available work is the specified constraint on production, an iterative approach is used to assure satisfaction of the constraint. This process is continued until the abandonment temperature or pressure is reached.
Application of this model to the East-Mesa reservoir shows a reasonable agreement with the more complex simulations done with an Intercomp model by Republic Geothermal, Inc. A sensitivity study of this reservoir has been made.