Impact of injection rate ramp-up on nucleation and arrest of dynamic fault slip

Fluid injection into underground formations reactivates preexisting geological discontinuities such as faults or fractures. In this work, we investigate the impact of injection rate ramp-up present in many standard injection protocols on the nucleation and potential arrest of dynamic slip along a planar pressurized fault. We assume a linear increasing function of injection rate with time, up to a given time $$t_c$$ t c after which a maximum value $$Q_m$$ Q m is achieved. Under the assumption of negligible shear-induced dilatancy and impermeable host medium, we solve numerically the coupled hydro-mechanical model and explore the different slip regimes identified via scaling analysis. We show that in the limit when fluid diffusion time scale $$t_w$$ t w is much larger than the ramp-up time scale $$t_c$$ t c , slip on an ultimately stable fault is essentially driven by pressurization at constant rate. Vice versa, in the limit when $$t_c/t_w \gg 1$$ t c / t w ≫ 1 , the pressurization rate, quantified by the dimensionless ratio $$\dfrac{Q_m t_w}{t_c Q^*}$$ Q m t w t c Q ∗ with $$Q^*$$ Q ∗ being a characteristic injection rate scale, does impact both nucleation time and arrest distance of dynamic slip. Indeed, for a given initial fault loading condition and frictional weakening property, lower pressurization rates delay the nucleation of a finite-sized dynamic event and increase the corresponding run-out distance approximately proportional to $$\propto \left( \dfrac{Q_m t_w}{t_c Q^*}\right) ^{-0.472}$$ ∝ Q m t w t c Q ∗ - 0.472 . On critically stressed faults, instead, the ramp-up of injection rate activates quasi-static slip which quickly turn into a run-away dynamic rupture. Its nucleation time decreases non-linearly with increasing value of $$\dfrac{Q_m t_w}{t_c Q^*}$$ Q m t w t c Q ∗ and it may precede (or not) the one associated with fault pressurization at constant rate only.

SEISMIC RESONANCE: WAVELET-FREE REFLECTIVITY RETRIEVAL VIA MODIFIED CEPSTRAL DECOMPOSITION

Spectral decomposition is a proven tool in seismic interpretation, aiding interpreters to highlight channels, map temporal bed thickness and other geological discontinuities. Once seismic data is spectrally decomposed, notch patterns in the amplitude spectra are indicative of the reservoir layer’s thickness and/or its interval velocity. Additional cepstral decomposition will allow direct extraction of bed time-thickness or arrival time under particular reflectivity series setup. We build on these observations to establish a more generalized workflow for reflectivity retrieval method without the need to understand the details of the wavelet, provided the starting seismic is stably phased via phase correction during processing. We demonstrate reflector time and its ‘apparent strength’ can be identified in a transformed seismic resonance domain resulted from a modified cepstrum analysis. In this domain, each reflector can be characterized from obvious linear hot spots. The timing and strength of those linear hot spots will reveal reflector times and scaled reflectivity coefficients. This new method is subsequently applied for thickness prediction of a target reservoir in a complex geological setting, with large thickness variations and weak impedance contrast with underlying lithology previously complicating identification of base-reservoir. In a deep-water field blind test, the sand thicknesses evaluated from this method are found to be close to true vertical thickness found in wells.

Determination of Rock Slope Stability Using Stereographic Projection between Sulav and Amadiyah Resorts

Geological discontinuities play a significant role in the assessment of rock slope stability. Rock slope stability has been studied on the main road between Sulav and Amadiya resorts in Duhok governorate on the southern limb of Mateen anticline, to determine the expected rock slides on this road. Five (5) stations were chosen to study these rock slides that may occur on these steep slopes. All these stations within Pila Spi Formation that consists of hard dolomitic limestone and covering the areas from Sulav resort towards Amadiya district with a length of up to 2.5 Kms. The Stereographic analysis was used to study and classify the stability of these slopes. The analysis showed in all stations the possibility of plane sliding to happen on the bedding plane, and the wedge sliding between the bedding plane and planes of all joint sets, as well as the occurrence of rockfall on some stations.

Effects of External Dynamic Disturbances And Structural Plane On Rock Fracturing Around Deep Underground Cavern

The occurrence of disasters in deep mining engineering has been confirmed to be closely related to the external dynamic disturbances and geological discontinuities. Thus, a combined finite-element approach was employed to simulate the failure process of an underground cavern, which provided insights into the failure mechanism of deep hard rock affected by factors such as the dynamic stress-wave amplitudes, disturbance direction, and dip angles of the structural plane. The crack-propagation process, stress-field distribution, displacement, velocity of failed rock, and failure zone around the circular cavern were analyzed to identify the dynamic response and failure properties of the underground structures. The simulation results indicated that the dynamic disturbance direction had less influence on the dynamic response for the constant in situ stress state, while the failure intensity and damage range around the cavern always exhibited a monotonically increasing trend with an increase in the dynamic load (stress-wave amplitudes). The crack distribution around the circular cavern exhibited an asymmetric pattern, possibly owing to the stress-wave reflection behavior and attenuation effect along the propagation route. Geological discontinuities significantly affected the stability of nearby caverns subjected to dynamic disturbances, during which the failure intensity exhibited the pattern of an initial increase followed by a decrease with an increase in the dip angle of the structural plane. Additionally, the dynamic disturbance direction led to variations in the crack distribution for specific structural planes and stress states. These results indicate that the failure behavior should be the integrated response of the excavation unloading effect, geological conditions, and external dynamic disturbances.

GeomInt: geomechanical integrity of host and barrier rocks–experiments, models and analysis of discontinuities

The present paper gives an overview of the GeomInt project “Geomechanical integrity of host and barrier rocks—experiment, modelling and analysis of discontinuities” which has been conducted from 2017–2020 within the framework of the “Geo:N Geosciences for Sustainability” program. The research concept of the collaborative project is briefly introduced followed by a summary of the most important outcomes. The research concept puts geological discontinuities into the centre of investigations—as these belong to the most interesting and critical elements for any subsurface utilisation. Thus, while research questions are specific, they bear relevance to a wide range of applications. The specific research is thus integrated into a generic concept in order to make the results more generally applicable and transferable. The generic part includes a variety of conceptual approaches and their numerical realisations for describing the evolution of discontinuities in the most important types of barrier rocks. An explicit validation concept for the generic framework was developed and realised by specific “model-experiment-exercises” (MEX) which combined experiments and models in a systematic way from the very beginning. 16 MEX have been developed which cover a wide range of fundamental fracturing mechanisms, i.e. swelling/shrinkage, fluid percolation, and stress redistribution processes. The progress in model development is also demonstrated by field-scale applications, e.g. in the analysis and design of experiments in underground research laboratories in Opalinus Clay (URL Mont Terri, Switzerland) and salt rock (research mine Springen, Germany).

Impact of injection rate ramp-up on nucleation and arrest of dynamic fault slip

Fluid injection into underground formations reactivates preexisting geological discontinuities such as faults or fractures. In this work, we investigate the impact of injection rate ramp-up present in many standard injection protocols on the nucleation and potential arrest of dynamic slip along a planar pressurized fault. We assume a linear increasing function of injection rate with time, up to a given time tc after which a maximum value Qm is achieved. Under the assumption of negligible shear-induced dilatancy and impermeable host medium, we solve numerically the coupled hydro-mechanical model and explore the different slip regimes identified via scaling analysis. We show that in the limit when fluid diffusion time scale tw is much larger than the ramp-up time scale tc, slip on an ultimately stable fault is essentially driven by pressurization at constant rate. Vice versa, in the limit when tc/tw ≫ 1, the pressurization rate, quantified by the dimensionless ratio (Qm tw / tc Q∗), does impact both nucleation time and arrest distance of dynamic slip. Indeed, for a given initial fault loading condition and frictional weakening property, lower pressurization rates delay the nucleation of a finite-sized dynamic event and increase the corresponding run-out distance approximately proportional to (Qm tw / tc Q∗)^(-0.472). On critically stressed faults, instead, the ramp-up of injection rate activates quasi-static slip which quickly turn into a run-away dynamic rupture. Its nucleation time decreases non-linearly with increasing value of (Qm tw / tc Q∗) and it may precede (or not) the one associated with fault pressurization at constant rate only.

Fall of Ground Management Through Underground Joint Mapping: Shallow Chrome Mining Case Study

Joint properties and useful mapping techniques are key to the fall of ground management in underground mining. The study analysed the fall of ground management techniques at the mine with a view to identify the causal factors of the falls of ground. This paper practically demonstrates how two mapping techniques were used to obtain joint data. A brief description of geological discontinuities at the study area is given in the paper. Joint mapping was carried out in both the North and South sections of the mine. Procedures followed when collecting joint data are also provided. The collected joint data was used to evaluate rock fall probability. Rockfall probabilistic analysis carried out in the study indicates that about 80% of all key blocks formed are 1m3 in size. Results show that larger blocks are more likely to fail through the rotation. Furthermore, to prevent small blocks from falling between support units, areal coverage is suggested in heavily jointed rock masses. Probabilistic analysis can be used to evaluate the probability of rock falls, and support design for stability enhancement. The research noted that at the centre of effective falls of ground management are accurate and precise structural/joint mapping. This research is part of an MSc Engineering study.

Review: Hydrogeology of weathered crystalline/hard-rock aquifers—guidelines for the operational survey and management of their groundwater resources

Hard rocks or crystalline rocks (i.e., plutonic and metamorphic rocks) constitute the basement of all continents, and are particularly exposed at the surface in the large shields of Africa, India, North and South America, Australia and Europe. They were, and are still in some cases, exposed to deep weathering processes. The storativity and hydraulic conductivity of hard rocks, and thus their groundwater resources, are controlled by these weathering processes, which created weathering profiles. Hard-rock aquifers then develop mainly within the first 100 m below ground surface, within these weathering profiles. Where partially or noneroded, these weathering profiles comprise: (1) a capacitive but generally low-permeability unconsolidated layer (the saprolite), located immediately above (2) the permeable stratiform fractured layer (SFL). The development of the SFL’s fracture network is the consequence of the stress induced by the swelling of some minerals, notably biotite. To a much lesser extent, further weathering, and thus hydraulic conductivity, also develops deeper below the SFL, at the periphery of or within preexisting geological discontinuities (joints, dykes, veins, lithological contacts, etc.). The demonstration and recognition of this conceptual model have enabled understanding of the functioning of such aquifers. Moreover, this conceptual model has facilitated a comprehensive corpus of applied methodologies in hydrogeology and geology, which are described in this review paper such as water-well siting, mapping hydrogeological potentialities from local to country scale, quantitative management, hydrodynamical modeling, protection of hard-rock groundwater resources (even in thermal and mineral aquifers), computing the drainage discharge of tunnels, quarrying, etc.

Improvement of operational exploration and mining safety using the shear-wave reflection method of mine seismology

During the whole life of mines, the geological survey departments are faced with the critical objectives of operational exploration. These objectives are sometimes impossible to be met adequately without geophysics. Neutrality of geophysical studies is often below the desired level. The shear-wave reflection method of mine seismology, with separation of reflections, which is developed by the present authors, allows: investigating the shape, thickness and dip angles of geologic beds and ore bodies, detecting geological discontinuities such as faults, joints and cavities, as well as determining physical and mechanical properties of rocks. The capacity of the shear-wave reflection method with separation of reflections is described in a case-study of the Upper Kama Potash Salt Deposit. The seismic tests were carried out in a few underground excavations. The developed method in the test conditions is accurate to the first meters at a spacing more than 300 m. Accordingly, it is possible to study the whole thickness of the salt strata in the Upper Kama Deposit, from the salt table to the anhydrous clay roof, specifically, to assess the structure and composition of the waterproof strata, to identify mine-threatening geological discontinuities, and to build the model of physical and mechanical properties of rock mass. Introduction of the newly developed method can greatly enhance efficiency of actual operational exploration. The method has yielded the best results during or immediately after heading of permanent and development openings. Considering similarity of surface and underground mineral mining in terms of access to solid rocks, it is expected to have the same outcomes in open pit mines.

Multidirectional eigenvalue-based coherence attribute for discontinuity detection

The conventional coherence attribute is typically applied to migrated full-stacked seismic data volumes to detect geological discontinuities. Recently, multispectral, multiazimuth, and multioffset coherence attributes have been proposed and implemented with different seismic data volumes of specific frequencies, azimuths, and offsets to enhance discontinuities. Generally, geological anomalies, such as faults and channels, will be better illuminated by a perpendicular rather than a parallel direction for computation. Therefore, we propose a multidirectional eigenvalue-based coherence attribute by establishing multiple covariance matrices along certain different directions on a single post-stack volume. We adopt two methods to compute multidirectional coherence attribute. One is to compute multiple coherence volumes in different directions and to define the minimum as the final multidirectional coherence. This method is time-consuming, but could provide partial and overall discontinuity simultaneously. The other method obtains one coherence volume by summing covariance matrices in different directions, which is computationally efficient, but only provides overall discontinuity. The performance of 3D physical model and field data volumes demonstrates that multidirectional coherence can highlight subtle geologic structures with a higher resolution than conventional coherence. This suggests that multidirectional coherence attribute may serve as an effective tool for detecting the distribution of geologic discontinuities in seismic interpretation.