Stress-based tomography: potential, open-questions and future developments

Author(s):  
Peter Bayer ◽  
Mohammad Javad Afshari Moein ◽  
Márk Somogyvári ◽  
Lisa Ringel ◽  
Mohammadreza Jalali

<p>Fracture network characterization is critical for many subsurface engineering problems in petroleum, mining, nuclear waste disposal and Enhanced Geothermal Systems (EGS). Due to limited exposure, direct measurement of fracture network properties at great depth is not possible and geophysical imaging techniques cannot resolve the fractures. Therefore, tomographic imaging techniques have been proposed and applied to reconstruct the structural discontinuities of rock mass. Stress-based tomography is a novel concept aiming at probabilistic imaging of the fracture network using the stress perturbations along deep boreholes. Currently, this approach has only been successfully tested on two-dimensional fracture networks. However, its great potential to unravel the heterogeneous structure of fractured rocks at great depth motivates further scientific effort. Here, we present the potential, open questions, current challenges and necessary future developments in order to apply this methodology to image three-dimensional multiscale structure of the rock mass in the field. Other tomographic approaches such as tracer and hydraulic tomography invert tracer breakthrough curves (BTCs) and pressure response in an observational well. We suggest a joint and comparative tomographic analysis in a Bayesian inversion framework to reconstruct Discrete Fracture Networks (DFN). This is expected to provide a new view of the strengths of each tomographic variant.</p>

2020 ◽  
Author(s):  
Kyung Won Chang ◽  
Gungor Beskardes ◽  
Chester Weiss

<p>Hydraulic stimulation is the process of initiating fractures in a target reservoir for subsurface energy resource management with applications in unconventional oil/gas and enhanced geothermal systems. The fracture characteristics (i.e., number, size and orientation with respect to the wellbore) determines the modified permeability field of the host rock and thus, numerical simulations of flow in fractured media are essential for estimating the anticipated change in reservoir productivity. However, numerical modeling of fluid flow in highly fractured media is challenging due to the explosive computational cost imposed by the explicit discretization of fractures at multiple length scales. A common strategy for mitigating this extreme cost is to crudely simplify the geometry of fracture network, thereby neglecting the important contributions made by all elements of the complex fracture system.</p><p>The proposed “Hierarchical Finite Element Method” (Hi-FEM; Weiss, Geophysics, 2017) reduces the comparatively insignificant dimensions of planar- and curvilinear-like features by translating them into integrated hydraulic conductivities, thus enabling cost-effective simulations with requisite solutions at material discontinuities without defining ad-hoc, heuristic, or empirically-estimated boundary conditions between fractures and the surrounding formation. By representing geometrical and geostatistical features of a given fracture network through the Hi-FEM computational framework, geometrically- and geomechanically-dependent fluid flow properly can now be modeled economically both within fractures as well as the surrounding medium, with a natural “physics-informed” coupling between the two.</p><p>SNL is managed and operated by NTESS under DOE NNSA contract DE-NA0003525.</p>


Solid Earth ◽  
2020 ◽  
Vol 11 (6) ◽  
pp. 2119-2140
Author(s):  
Billy James Andrews ◽  
Zoe Kai Shipton ◽  
Richard Lord ◽  
Lucy McKay

Abstract. Fault architecture and fracture network evolution (and resulting bulk hydraulic properties) are highly dependent on the mechanical properties of the rocks at the time the structures developed. This paper investigates the role of mechanical layering and pre-existing structures on the evolution of strike–slip faults and fracture networks. Detailed mapping of exceptionally well exposed fluvial–deltaic lithologies at Spireslack Surface Coal Mine, Scotland, reveals two phases of faulting with an initial sinistral and later dextral sense of shear with ongoing pre-faulting, syn-faulting, and post-faulting joint sets. We find fault zone internal structure depends on whether the fault is self-juxtaposing or cuts multiple lithologies, the presence of shale layers that promote bed-rotation and fault-core lens formation, and the orientation of joints and coal cleats at the time of faulting. During ongoing deformation, cementation of fractures is concentrated where the fracture network is most connected. This leads to the counter-intuitive result that the highest-fracture-density part of the network often has the lowest open fracture connectivity. To evaluate the final bulk hydraulic properties of a deformed rock mass, it is crucial to appreciate the relative timing of deformation events, concurrent or subsequent cementation, and the interlinked effects on overall network connectivity.


Energies ◽  
2021 ◽  
Vol 14 (18) ◽  
pp. 5725
Author(s):  
Rafał Moska ◽  
Krzysztof Labus ◽  
Piotr Kasza

Hydraulic fracturing (HF) is a well-known stimulation method used to increase production from conventional and unconventional hydrocarbon reservoirs. In recent years, HF has been widely used in Enhanced Geothermal Systems (EGS). HF in EGS is used to create a geothermal collector in impermeable or poor-permeable hot rocks (HDR) at a depth formation. Artificially created fracture network in the collector allows for force the flow of technological fluid in a loop between at least two wells (injector and producer). Fluid heats up in the collector, then is pumped to the surface. Thermal energy is used to drive turbines generating electricity. This paper is a compilation of selected data from 10 major world’s EGS projects and provides an overview of the basic elements needed to design HF. Authors were focused on two types of data: geological, i.e., stratigraphy, lithology, target zone deposition depth and temperature; geophysical, i.e., the tectonic regime at the site, magnitudes of the principal stresses, elastic parameters of rocks and the seismic velocities. For each of the EGS areas, the scope of work related to HF processes was briefly presented. The most important HF parameters are cited, i.e., fracturing pressure, pumping rate and used fracking fluids and proppants. In a few cases, the dimensions of the modeled or created hydraulic fractures are also provided. Additionally, the current state of the conceptual work of EGS projects in Poland is also briefly presented.


2020 ◽  
Author(s):  
Taghi Shirzad ◽  
Stanisław Lasocki ◽  
Beata Orlecka‐Sikora

<p> In Enhanced Geothermal Systems pressurized injections play a role in developing fracture networks and enhancing the water transmissivity. However, the fractures may also coalesce into undesired pathways for fluid migration to enable the fluids to reach pre-existing faults. The properties of observed seismicity can shed some light on the fracture network development and from the standpoint of the possibility to form such undesired pathways. However, to reach this goal the seismic events should be well parameterized. In particular, the information on fault plane mechanisms is essential, which is often not readily accessible. In this study, we use the rupturing process with an accurate P-wave velocity model, which is obtained by the first arrival P-wave tomography approach, to compensate for an eventual lack of source mechanisms of micro-events. For this purpose, four characteristics of the sources (final/average displacement on the fault, the dimension of fault, rupture velocity and particle velocity) can be considered. A 3D model is defined around the hypocenter of each event, so that the size of this model directly depends on the event magnitude. After calculating the arrival time of the selected phase (e.g., P, S, p or s) for each station, all waveforms are then aligned, and stacked by different stacking (e.g., phase weight, N<sup>th</sup>-root) methods. By considering the maximum amplitude of the stacked waveform which is stimulated by each grid, the rupturing plane and the average velocity of rupturing can be obtained. This information of source can be replaced by the double-couple mechanism to investigate the fractures linking and tracking.</p><p>This work was supported under the <em>S4CE</em>: "Science for Clean Energy" project, which has received funding from the European Union’s Horizon 2020 research and innovation program, under grant agreement No 764810.</p>


2021 ◽  
Author(s):  
Amy Myers ◽  
Claire Harnett ◽  
Michael Heap ◽  
Eoghan Holohan ◽  
Thomas Walter

<p>Volcanic domes form when lava is too viscous to flow away from an active volcanic vent; instead, the lava accumulates into a mound consisting of a hotter, ductile core and a colder, brittle outer layer. An existing lava dome grows when new material is injected into the core of the dome, causing the  outer layer to stretch and develop tensile fractures. With continued dome growth, these weaknesses can propagate to form an extensive fracture network and the dome may fail. Collapse events often generate rock falls and debris avalanches, lahars, and high-speed pyroclastic flows, endangering populations residing at the base of a volcano. Since such fractures represent potential failure planes, in this project we aim to understand the role they have in destabilising lava domes.</p><p>This project will build on the work published by Harnett et al. (2018), which demonstrates the suitability of a discrete element modelling approach to simulate dome emplacement and evolution. Specifically, this project is designed to:</p><p>1. Use high-resolution photogrammetry to characterise the possible fracture states of a dome;</p><p>2. Establish up-scaled rock-mass properties by performing geomechanical experiments on both fractured and non-fractured samples of dome rock from prior collapses;</p><p>3. Develop a numerical model to investigate how the presence and properties of fracture networks influence dome stability.</p><p>The model, developed using PFC, will be used to identify critical fracture states that can signify a dome collapse is likely to occur. Under the current model, parallel bonds simulate the fluid magma core and flat joints simulate the solid talus material. This project will build on this original model by incorporating discrete fracture networks into the smooth-joint model to implement dome fracturing. The new model will look to investigate the effect  of a fracture network on a static dome that, when in its unfractured state, is stable under gravity. Additionally, the model will be designed such that inputs can include experimentally derived rock-mass properties. It is hoped that, by incorporating observational and experimental data into a more  complex model, the dynamic evolution of fractures in a growing lava dome can be investigated and the ongoing likelihood of a dome collapse event can be assessed.</p><p> </p><p>Harnett, C. E. et al., 2018. J. Volcanol. Geoth. Res., 359: 68-77.</p>


Sign in / Sign up

Export Citation Format

Share Document