SUMMARYOne proposed geothermal heat extraction scheme relies on water circulation in a large vertical crack created by hydraulic fracturing in a hot dry impermeable rock mass. Water flow, heat convection and crack opening widths are analysed by finite elements. Governing field equations of the problem are first set up rigorously and then various small terms are identified and neglected, retaining the effects of pressure gradient, buoyancy, velocity head (kinetic energy) and head loss due to viscous friction in the water flow equation, and the effects of heat convection in water and heat conduction in rock in the heat transfer equation. The finite element scheme for water flow is based on a variational principle that is typical for diffusion problems, and for heat transfer it is based on the method of least-square residuals. The system of differential equations is highly non-linear. The non-linear terms and coefficients are treated in the finite element analysis as constant; the finite element analysis of the steady-state pressures, fluxes and temperatures is then iterated, evaluating all non-linear terms and coefficients on the basis of the solution obtained in the previous iteration. Numerically calculated fields at various times after the start of cooling are presented. They indicate some features favourable for the geothermal scheme, such as formation of eddy currents, and downward flux of water toward hotter rock. However, other important questions would have to be solved to gain full understanding of this proposed geothermal scheme.
INTRODUCflONIn one proposed scheme of heat extraction from hot dry rock geothermal reservoirs, a large vertical crack is to be induced by hydraulic fracturing from a borehole (Figure 1). The upper part of the crack is then intersected by a second borehole, and heat is extracted by circulating water through the crack. Influx of cold water is expected to cause secondary fracturing and large expansion of the region from which the heat can be withdrawn.Considerable interest has been recently shown in this scheme and various aspects of the problem have been studied. 1-3 One problem of particular interest is the circulation of water, heat convection and opening width of a single, main crack created by hydraulic fracture. This problem was analysed by McFarland 4 using a standard finite difference program for diffusion flows. The present paper deals with finite element analysis of the problem. Variations of the crack thickness due to temperature changes are taken into account, and differential equations describing the flow and heat convection in the crack more accurately than those used before are developed.