Fracture Size Emerges as Key to Long-Lasting Geothermal Heating Systems

Around 100 degrees Celsius, deep crystalline rock holds enough heat to warm entire cities, but only if engineers can pull that energy out of the ground without draining the reservoir prematurely. A new open-access study published in GeoEnergy Communications examines exactly this problem, using an analytical model to test how the design choices of a deep petrothermal multi-frac enhanced geothermal system shape its thermal output over decades of district heating service. The work, led by Shicheng Xing and colleagues, shifts attention away from the high-temperature reservoirs above 100 degrees Celsius that dominate electricity-focused geothermal research and toward the medium-temperature regime that suits heating applications in medium-sized cities.

The researchers built their analysis around an analytical solution for the temperature and energy evolution within fractures and the surrounding rock matrix, drawing on a Laplace-transform formulation adapted from earlier work by Akdas and Onur. In their framework, a horizontal or inclined well doublet intersects a set of parallel hydraulic fractures, and the total injection flow rate is divided evenly among them. Water injected at a set temperature flows along each fracture, absorbs heat from the hot rock, and returns to the surface, where the fracture outlet temperature and thermal power output can be calculated directly. The team verified this analytical scheme against the numerical distinct element code 3DEC, simulating a single fracture with two injection stages of hot and cooler water into rock at 30 degrees Celsius, and found nearly identical temperature responses from both approaches.

The baseline configuration reflects conditions typical of crystalline formations such as granite or gneiss, which make up nearly 70 percent of the upper three kilometres of the crust. System parameters were chosen so that output power falls between roughly 10 and 140 megawatts, a range matching the central heating demand of metropolitan areas, with a return flow temperature of 60 degrees Celsius. The model deliberately simplifies nature: rock properties are homogeneous and impermeable, fracture geometry is idealised, thermal properties and aperture remain fixed, and the natural heat flux from the Earth’s interior is neglected because its contribution, about 0.1 degrees Celsius per year, is marginal over the 30 to 50 year lifetimes typically demanded of such systems.

With this framework in place, the sensitivity analysis varied four key parameters of the underground heat exchanger: injection flow rate, fracture size, fracture spacing, and initial rock temperature. The results reveal a recurring tension between immediate performance and long-term sustainability. Higher flow rates deliver greater thermal power in the early years of operation, because more water circulates through the fractures and absorbs heat faster. Yet that same speed of extraction cools the rock near the fractures more quickly, shrinking the temperature difference between rock and fluid and accelerating the decline of both outlet temperature and power output. Lower flow rates start smaller but sustain their output for longer, benefiting from slower heat extraction and better thermal recharge from the surrounding rock.

Fracture size proved to be the decisive lever. When the team compared fracture areas ranging from 150 by 150 metres up to 400 by 400 metres, keeping flow rate, spacing, and initial rock temperature constant, larger fractures showed markedly slower declines in outlet temperature and maintained higher power output over the 50-year simulation window. All systems began with approximately the same initial power, since flow rate and rock temperature were identical, but the larger contact area between fluid and rock matrix allowed heat to be drawn from a wider region of the reservoir, so the temperature around the fractures decreased more gradually. This effect persisted even when the flow rate was raised by a factor of 1.4, although the higher flow rate still accelerated thermal depletion compared with the baseline case.

Fracture spacing, by contrast, offered a more modest benefit. Increasing spacing from 60 to 110 metres slowed the rate of decline in outlet temperature and power output, because wider spacing lets heat be extracted from a broader swath of the reservoir and reduces rapid localised cooling around individual fractures. However, the benefit diminished as spacing continued to grow, and the authors found that enlarging spacing was consistently less effective than enlarging fracture length. A larger spacing reduces localised cooling but does not directly increase the contact area between fluid and matrix the way a longer fracture does. The authors also caution that compressing the spacing too far makes the system inefficient.

The fourth parameter, initial rock temperature, which correlates with reservoir depth, produced the same double-edged pattern as flow rate. Across a temperature span from 75 to 155 degrees Celsius, hotter reservoirs delivered stronger initial power output, but they also cooled faster. The larger temperature difference between the injected fluid and the rock matrix accelerates heat extraction, causing outlet temperature to fall more rapidly over time. Deeper placement of the heat exchanger therefore buys a stronger start at the cost of a steeper decline, reinforcing the conclusion that no single parameter can be optimised in isolation.

The study’s authors are candid about the limits of their model. It assumes homogeneous rock properties and idealised parallel fractures, and it does not incorporate reservoir heterogeneity, fracture roughness, or non-parallel well configurations such as anti-parallel or oblique alignments, all of which may significantly influence heat extraction efficiency under real reservoir conditions. Validation so far rests on numerical comparisons rather than calibration against field data. The authors also note that the strongly non-linear hydro-thermal coupling of these systems demands specific simulations for optimisation, and that simple interpolation or extrapolation from their results should be used with caution.

Despite these caveats, the analysis points toward concrete design strategies. Because several considered configurations lead to significant thermal drawdown, particularly those with small fractures that restrict the area available for heat exchange, the authors recommend combining heat extraction with thermal storage technologies to slow or compensate the cooling process. Where the desired output temperature cannot be reached, they suggest evaluating whether heat pumps are necessary. To limit risk and build redundancy, they propose a multiple-doublet configuration, in which several well pairs share the load, restricting the number of fractures per doublet and thereby limiting borehole length and the risk of induced seismicity.

For cities weighing geothermal district heating, the message is that reservoir longevity is designed, not given. Enlarging fracture size stands out as the most effective single strategy for improving long-term thermal sustainability, while flow rate and depth trade early strength against faster exhaustion. Balancing heat extraction against long-term cooling of the heat exchanger, the authors conclude, is the most critical task in maximising thermal energy output while ensuring a sufficient system lifetime, and the parameter choices made before drilling will determine whether a multi-frac system heats a city for decades or fades within a fraction of that span.

Subject of Research: Parameter sensitivity of heat extraction efficiency in deep petrothermal multi-frac enhanced geothermal systems for district heating

Article Title: How different parameters influence the efficiency of heat extraction using deep petrothermal multi-frac systems for central heating purposes

Article References: Xing, S., Zhang, M., Tiedtke, F., Konietzky, H., Liu, Q., & Wang, F. (2025). How different parameters influence the efficiency of heat extraction using deep petrothermal multi-frac systems for central heating purposes. GeoEnergy Communications, 1(1), Article 8. https://doi.org/10.1007/s44421-025-00008-z

Image Credits: AI Generated

DOI: 10.1007/s44421-025-00008-z

Keywords: enhanced geothermal systems, district heating, petrothermal reservoirs, fracture size, flow rate, thermal drawdown, analytical model, sensitivity analysis, crystalline rock, induced seismicity, heat storage, multi-frac EGS