Breakthrough time of a geothermal reservoir

Abstract

This research identifies and provides a relative ranking of the parameters that control the thermal breakthrough time of a geothermal doublet. The ranking is based on simulations modelled by a three-dimensional model build in COMSOL Multiphysics with a simulation duration of 50 years. The model is a nonisothermal, isotropic, cuboid consisting of a sandstone aquifer surrounded by identical impermeable layers. The ranking is derived based on a solution space covering a minimum, maximum and mean value for each of the simulation parameters. The investigated simulation parameters are the density, porosity, permeability, thermal conductivity and specific heat capacity of both the surrounding rock formations and the aquifer, the depth and thickness of the aquifer, the flow rate, injection temperature and the well spacing. The mean values of the solution space form the base model, from which the parameters will divert separately to monitor the model’s behaviour on the varying parameters.

In this research, the breakthrough time is defined as the time at which the production temperature is decreased by 1% to 99% of its initial value. The ranking is based on comparing the proportional change of each parameter from the base model to the change in breakthrough time. The results show that the thickness, flow rate and the well spacing are the most crucial parameters influencing the thermal breakthrough time of a reservoir. Overall, the flow rate has the greatest impact, with a decrease of 16.7 years between a flow rate of 150 m3/hour and 250 m3/hour.

In addition, the surrounding rock parameters have a notably smaller impact on the thermal breakthrough time compared to the reservoir rock and process parameters. The surrounding rock parameter with the relatively largest impact on the breakthrough time is the specific heat capacity, which is only a change of 0.2 years between a specific heat capacity of 1150 and 1250 Jkg-1K-1.