Large-scale dynamic hydrofracturing, healing and fracture network characterization

Permeability of a rock is a dynamic property that varies spatially and temporally. Fractures provide the most efficient channels for fluid flow and thus directly contribute to the permeability of the system. Fractures usually form as a result of a combination of tectonic stresses, gravity (i.e. lithostatic pressure) and fluid pressures. High pressure gradients alone can cause fracturing, the process which is termed as hydrofracturing that can determine caprock (seal) stability or reservoir integrity. Fluids also transport mass and heat, and are responsible for the formation of veins by precipitating minerals within open fractures. Veining (healing) thus directly influences the rock’s permeability. Upon deformation these closed factures (veins) can refracture and the cycle starts again. This fracturing-healing-refacturing cycle is a fundamental part in studying the deformation dynamics and permeability evolution of rock systems. This is generally accompanied by fracture network characterization focusing on network topology that determines network connectivity. Fracture characterization allows to acquire quantitative and qualitative data on fractures and forms an important part of reservoir modeling. This thesis highlights the importance of fracture-healing and veins’ mechanical properties on the deformation dynamics. It shows that permeability varies spatially and temporally, and that healed systems (veined rocks) should not be treated as fractured systems (rocks without veins). Field observations also demonstrate the influence of contrasting mechanical properties, in addition to the complexities of vein microstructures that can form in low-porosity and permeability layered sequences. The thesis also presents graph theory as a characterization method to obtain statistical measures on evolving network connectivity. It also proposes what measures a good reservoir should have to exhibit potentially large permeability and robustness against healing. The results presented in the thesis can have applications for hydrocarbon and geothermal reservoir exploration, mining industry, underground waste disposal, CO2 injection or groundwater modeling.

Geothermal district heating networks: modelling novel operational strategies incorporating heat storage

The value of integrating a heat storage into a geothermal district heating system has been investigated. The behaviour of the system under a novel operational strategy has been simulated focusing on the energetic, economic and environmental effects of the new strategy of incorporation of the heat storage within the system. A typical geothermal district heating system consists of several production wells, a system of pipelines for the transportation of the hot water to end-users, one or more re-injection wells and peak-up devices (usually fossil-fuel boilers). Traditionally in these systems, the production wells change their production rate throughout the day according to heat demand, and if their maximum capacity is exceeded the peak-up devices are used to meet the balance of the heat demand. In this study, it is proposed to maintain a constant geothermal production and add heat storage into the network. Subsequently, hot water will be stored when heat demand is lower than the production and the stored hot water will be released into the system to cover the peak demands (or part of these). It is not intended to totally phase-out the peak-up devices, but to decrease their use, as these will often be installed anyway for back-up purposes. Both the integration of a heat storage in such a system as well as the novel operational strategy are the main novelties of this thesis. A robust algorithm for the sizing of these systems has been developed. The main inputs are the geothermal production data, the heat demand data throughout one year or more and the topology of the installation. The outputs are the sizing of the whole system, including the necessary number of production wells, the size of the heat storage and the dimensions of the pipelines amongst others. The results provide several useful insights into the initial design considerations for these systems, emphasizing particularly the importance of heat losses. Simulations are carried out for three different cases of sizing of the installation (small, medium and large) to examine the influence of system scale. In the second phase of work, two algorithms are developed which study in detail the operation of the installation throughout a random day and a whole year, respectively. The first algorithm can be a potentially powerful tool for the operators of the installation, who can know a priori how to operate the installation on a random day given the heat demand. The second algorithm is used to obtain the amount of electricity used by the pumps as well as the amount of fuel used by the peak-up boilers over a whole year. These comprise the main operational costs of the installation and are among the main inputs of the third part of the study. In the third part of the study, an integrated energetic, economic and environmental analysis of the studied installation is carried out together with a comparison with the traditional case. The results show that by implementing heat storage under the novel operational strategy, heat is generated more cheaply as all the financial indices improve, more geothermal energy is utilised and less fuel is used in the peak-up boilers, with subsequent environmental benefits, when compared to the traditional case. Furthermore, it is shown that the most attractive case of sizing is the large one, although the addition of the heat storage most greatly impacts the medium case of sizing. In other words, the geothermal component of the installation should be sized as large as possible. This analysis indicates that the proposed solution is beneficial from energetic, economic, and environmental perspectives. Therefore, it can be stated that the aim of this study is achieved in its full potential. Furthermore, the new models for the sizing, operation and economic/energetic/environmental analyses of these kind of systems can be used with few adaptations for real cases, making the practical applicability of this study evident. Having this study as a starting point, further work could include the integration of these systems with end-user demands, further analysis of component parts of the installation (such as the heat exchangers) and the integration of a heat pump to maximise utilisation of geothermal energy.

Quantifying uplift and denudation of a thermally heterogeneous crust: a detailed multi-thermochronometeric study of central west Britain

Topography is often thought as exclusively linked to mountain ranges formed by plates collision. It is now, however, known that apart from compression, uplift and denudation of rocks may be triggered by rifting, like it happens at elevated passive margins, and away from plate boundaries by both intra-plate stress causing reactivation of older structures, and by epeirogenic movements driven by mantle dynamics and initiating long-wavelength uplift. In the Cenozoic, central west Britain and other parts of the North Atlantic margins experienced multiple episodes of rock uplift and denudation that have been variable both at spatial and temporal scales. The origin of topography in central west Britain is enigmatic, and because of its location, it may be related to any of the processes mentioned above. In this study, three low temperature thermochronometers, the apatite fission track (AFT) and apatite and zircon (U-Th-Sm)/He (AHe and ZHe, respectively) methods were used to establish the rock cooling history from 200◦C to 30◦C. The samples were collected from the intrusive rocks in the high elevation, high relief regions of the Lake District (NW England), southern Scotland and northern Wales. AFT ages from the region are youngest (55–70 Ma) in the Lake District and increase northwards into southern Scotland and southwards in north Wales (>200 Ma). AHe and ZHe ages show no systematic pattern; the former range from 50 to 80 Ma and the latter tend to record the post-emplacement cooling of the intrusions (200–400 Ma). The complex, multi-thermochronometric inverse modelling suggests a ubiquitous, rapid Late Cretaceous/early Palaeogene cooling event that is particularly marked in Lake District and Criffell. The timing and rate of cooling in southern Scotland and in northern Wales is poorly resolved as the amount of cooling was less than 60◦C. The Lake District plutons were at >110◦C prior to the early Palaeogene; cooling due to a combined effect of high heat flow, from the heat producing granite batholith, and the blanketing effect of the overlying low conductivity Late Mesozoic limestones and mudstones. Modelling of the heat transfer suggests that this combination produced an elevated geothermal gradient within the sedimentary rocks (50–70◦C/km) that was about two times higher than at the present day. Inverse modelling of the AFT and AHe data taking the crustal structure into consideration suggests that denudation was the highest, 2.0–2.5 km, in the coastal areas of the Lake District and southern Scotland, gradually decreasing to less than 1 km in the northern Southern Uplands and northern Wales. Both the rift-related uplift and the intra-plate compression poorly correlate with the timing, location and spatial distribution of the early Palaeogene denudation. The pattern of early Palaeogene denudation correlates with the thickness of magmatic underplating, if the changes of mean topography, Late Cretaceous water depth and eroded rock density are taken into consideration. However, the uplift due to underplating alone cannot fully justify the total early Palaeogene denudation. The amount that is not ex- plained by underplating is, however, roughly spatially constant across the study area and can be referred to the transient thermal uplift induced by the mantle plume arrival. No other mechanisms are required to explain the observed pattern of denudation. The onset of denudation across the region is not uniform. Denudation started at 70–75 Ma in the central part of the Lake District whereas the coastal areas the rapid erosion appears to have initiated later (65–60 Ma). This is ~10 Ma earlier than the first vol- canic manifestation of the proto-Iceland plume and favours the hypothesis of the short period of plume incubation below the lithosphere before the volcanism. In most of the localities, the rocks had cooled to temperatures lower than 30◦C by the end of the Palaeogene, suggesting that the total Neogene denudation was, at a maximum, several hundreds of metres. Rapid cooling in the last 3 million years is resolved in some places in southern Scotland, where it could be explained by glacial erosion and post-glacial isostatic uplift.