Correlation between the rock mass properties and maximum horizontal stress:A case study of overcoring stress measurements

Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stress,the three-dimensional(3D)stress tensors at 89 measuring points determined using an improved overcoring technique in nine mines in China were adopted,a newly defined characteristic parameter C_(ERP)was proposed as an indicator for evaluating the structural properties of rock masses,and a fuzzy relation matrix was established using the information distribution method.The results indicate that both the vertical stress and horizontal stress exhibit a good linear growth relationship with depth.There is no remarkable correlation between the elastic modulus,Poisson's ratio and depth,and the distribution of data points is scattered and messy.Moreover,there is no obvious relationship between the rock quality designation(RQD)and depth.The maximum horizontal stress σ_(H) is a function of rock properties,showing a certain linear relationship with the C_(ERP)at the same depth.In addition,the overall change trend of σ_(H) determined by the established fuzzy identification method is to increase with the increase of C_(ERP).The fuzzy identification method also demonstrates a relatively detailed local relationship betweenσ_H and C_(ERP),and the predicted curve rises in a fluctuating way,which is in accord well with the measured stress data.

Monitoring the change in horizontal stress with multi-wave time-lapse seismic response based on nonlinear elasticity theory

Monitoring the change in horizontal stress from the geophysical data is a tough challenge, and it has a crucial impact on broad practical scenarios which involve reservoir exploration and development, carbon dioxide (CO_(2)) injection and storage, shallow surface prospecting and deep-earth structure description. The change in in-situ stress induced by hydrocarbon production and localized tectonic movements causes the changes in rock mechanic properties (e.g. wave velocities, density and anisotropy) and further causes the changes in seismic amplitudes, phases and travel times. In this study, the nonlinear elasticity theory that regards the rock skeleton (solid phase) and pore fluid as an effective whole is used to characterize the effect of horizontal principal stress on rock overall elastic properties and the stress-dependent anisotropy parameters are therefore formulated. Then the approximate P-wave, SV-wave and SH-wave angle-dependent reflection coefficient equations for the horizontal-stress-induced anisotropic media are proposed. It is shown that, on the different reflectors, the stress-induced relative changes in reflectivities (i.e., relative difference) of elastic parameters (i.e., P- and S-wave velocities and density) are much less than the changes in contrasts of anisotropy parameters. Therefore, the effects of stress change on the reflectivities of three elastic parameters are reasonably neglected to further propose an AVO inversion approach incorporating P-, SH- and SV-wave information to estimate the change in horizontal principal stress from the corresponding time-lapse seismic data. Compared with the existing methods, our method eliminates the need for man-made rock-physical or fitting parameters, providing more stable predictive power. 1D test illustrates that the estimated result from time-lapse P-wave reflection data shows the most reasonable agreement with the real model, while the estimated result from SH-wave reflection data shows the largest bias. 2D test illustrates the feasibility of the proposed inversion method for estimating the change in horizontal stress from P-wave time-lapse seismic data.

Gradient principle of horizontal stress inducing rock burst in coal mine

Based on the stress distribution characteristics of rock burst multiple sites, the criterion of horizontal stress inducing layer dislocation rock burst was established. Accordingly, the influencing factors were analyzed. The analysis results indicate that the stress condition, edge of etastic zone depth, supporting strength, and the friction angle and cohesion among coal stratum, roof and floor are sensitive factors. By introducing double-couple model, the layer dislocation rock burst was explained and the energy radiation characteristics were analyzed. The SOS micro-seismic monitoring system was applied to observe the rock burst hazards about a mining face. The results show that P- and S-wave energy radiations produced by rock burst have directional characteristics. The energy radiation characteristics of the 22 rock bursts occurring on 79Z6 long-wall face are basically the same as theoretical results, that is, the ratio of S-wave energy of sensor 4 to 6 is about 1.5 and that of P-wave is smaller than 0.5. The consistency of the monitored characteristics of the energy radiation theoretically increases with the total energy increasing.

Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

Analytical investigations of in situ stress inversion from borehole breakout geometries

This study aims to investigate the feasibility of deriving in situ horizontal stresses from the breakout width and depth using the analytical method.Twenty-three breakout data with different borehole sizes were collected and three failure criteria were studied.Based on the Kirsch equations,relatively accurate major horizontal stress(sH)estimations from known minor horizontal stress(sh)were achieved with percentage errors ranging from 0.33%to 44.08%using the breakout width.The Mogi-Coulomb failure criterion(average error:13.1%)outperformed modified Wiebols-Cook(average error:19.09%)and modified Lade(average error:18.09%)failure criteria.However,none of the tested constitutive models could yield reasonable sh predictions from known sH using the same approach due to the analytical expression of the redistributed stress and the nature of the constitutive models.In consideration of this issue,the horizontal stress ratio(sH/sh)is suggested as an alternative input,which could estimate both sH and sh with the same level of accuracy.Moreover,the estimation accuracies for both large-scale and laboratory-scale breakouts are comparable,suggesting the applicability of this approach across different breakout sizes.For breakout depth,conformal mapping and complex variable method were used to calculate the stress concentration around the breakout tip,allowing the expression of redistributed stresses using binomials composed of sH and sh.Nevertheless,analysis of the breakout depth stabilisation mechanism indicates that additional parameters are required to utilise normalised breakout depth for stress estimation compared to breakout width.These parameters are challenging to obtain,especially under field conditions,meaning utilising normalised breakout depth analytically in practical applications faces significant challenges and remains infeasible at this stage.Nonetheless,the normalised breakout depth should still be considered a critical input for any empirical and statistical stress estimation method given its significant correlation with horizontal stresses.The outcome of this paper is expected to contribute valuable insights into the breakout stabilisation mechanisms and estimation of in situ stress magnitudes based on borehole breakout geometries.

An improved gravity method to horizontal tectonic stresses and its applications in Tibet

It is significant for identifying mass movement patterns to invert horizontal tectonic stresses at different depths underneath Tibet.In recent years,a large number of achievements focusing on two-dimensional tectonic stresses have been obtained from gravity data.However,three-dimensional tectonic stresses in Tibet are still unknown or debatable.Therefore,in the present study an improved method to multilayer horizontal tectonic stresses using gravity observations is developed.The inverted multilayer horizontal tectonic stresses are in agreement with those from previous studies.In addition,rich tectonic structure and development can be revealed from the inverted multilayer horizontal tectonic stresses:(1)the distribution of horizontal tectonic stresses at various depths shows strong correlation with that of the tectonic elements,where major faults and earthquake epicenters are corresponding with stress highs and the stable basins are consistent with stress lows.(2)the mass movement patterns of whole Tibet present clockwise,and the material movement directions in the west and east are approximately southnorth and east-west,respectively.(3)in eastern Tibet,the eastward materials caused by the south-north extrusion between Indian and Eurasian plates are divided into two parts by the stable Sichuan Block,one flowing nearly southeast and the other moving almost northeast.The inverted multilayer horizontal tectonic stresses may provide direct evidences for mass movement patterns in Tibet.

Geostress measurements near fault areas using borehole stress-relief method

To minimize negative effects of geostress distribution on mining safety near the fault areas, the UPM40 triaxial geostress testing system was introduced to conduct in-situ geostress measurements at three sites and nine points by the borehole stress-relief method. The results of strain?confining pressure curves show that rock masses at the three measuring sites exhibit comprehensive linear elasticity in spite of various fissures or cracks within rocks. Horizontal and vertical stress components distribute discrepantly near the fault areas, and the maximum lateral pressure coefficient is as high as 6.15. The maximum principle stress ranges from 8.01 to 14.93 MPa, and stress directions are in the range of N78.07°W?N17.55°W. Geostresses near fault areas are dominated by the horizontal tectonic stresses, while the lower values, compared to those under similar geological conditions are due to stress release by the fault. Additionally, the fault and shear stress nearby are partially responsible for asymmetric elongation and southwesterly migration of orebodies.

Ground movement mechanism in tectonic stress metal mines with steep structure planes

When mining metal mines with steep structure planes by the caving method,there is a mechanical model in which the horizontal stress on the rock mass is simplified as a column before surface subsidence.The model is used to deduce critical support load and limiting column length for a given horizontal stress and support pressure.Considering the impact of the column effect,a method is proposed to determine the movement of the ground and caving area in a mine.After surface subsidence,the horizontal stress on a surrounding rock mass can be simplified to a cantilever beam mechanical model.Expressions for its bending fracture length are deduced,and a method is given to determine its stability.On this basis,an explanation for the large ground movement and subsidence scope was given.A case study shows that the damage effect of column and cantilever beam is significant for ground movement in metal-ore mine,and an appropriate correction value should be applied when designing for its angle of ground movements.

Characteristics of orthorhombic anisotropic seismic response induced by horizontal in situ stress in vertical transversely isotropic media

Sedimentary strata typically exhibit the characteristics of transverse isotropy(VTI)with a vertical axis of symmetry.However,fractures in sedimentary strata tend to produce anisotropic closure due to horizontal in situ stress,resulting in pronounced orthorhombic anisotropy in VTI media under such stress conditions and influencing the propagation behavior of seismic waves.Previous studies have primarily focused on the elastic wave velocity anisotropy induced by applied stress in isotropic background media,neglecting the impact of VTI background media on the anisotropy induced by horizontal in situ stress and the response characteristics of seismic wave propagation.To address these gaps,we first establish the effective elastic stiffness tensor of VTI media under horizontal in situ stress using nonlinear acoustoelastic theory.Then,we derive the accurate and linearized approximate equations for P-wave seismic reflectivity of VTI media under horizontal in situ stress,based on wave equations and scattering theory,respectively.Finally,we compare and analyze the characteristics of orthorhombic anisotropic seismic response induced by horizontal in situ stress at various types of elastic reflection interfaces.Our results demonstrate that the linearized approximation of the seismic reflection response characteristics closely aligns with the accurate equations under conditions of small stress below 10 MPa,effectively capturing the azimuth-dependent orthorhombic anisotropy induced by horizontal in situ stress in VTI media.The results of this study also provide a novel theoretical approach and valuable insights into the seismic prediction of in situ stress.

Evolution of stress fields during the supercontinent cycle

We investigate the evolution of stress fields during the supercontinent cycle using the 2D Cartesian geometry model of thermochemical convection with the non-Newtonian rheology in the presence of floating deformable continents.In the course of the simulation,the supercontinent cycle is implemented several times.The number of continents considered in our model as a function of time oscillates around 3.The lifetime of a supercontinent depends on its dimension.Our results suggest that immediately before a supercontinent breakup,the over-lithostatic horizontal stresses in it(referring to the mean value by the computational area)are tensile and can reach-250 MPa.At the same time,a vast area beneath a supercontinent with an upward flow exhibits clearly the over-lithostatic compressive horizontal stresses of 50-100 МРа.The reason for the difference in stresses in the supercontinent and the underlying mantle is a sharp difference in their viscosity.In large parts of the mantle,the over-lithostatic horizontal stresses are in the range of±25 MPa,while the horizontal stresses along subduction zones and continental margins are significantly larger.During the process of continent-to-continent collisions,the compressive stresses can approximately reach 130 MPa,while within the subcontinental mantle,the tensile over-lithostatic stresses are about-50 MPa.The dynamic topography reflects the main features of the su-percontinent cycle and correlates with real ones.Before the breakup and immediately after the disin-tegration of the supercontinent,continents experience maximum uplift.During the supercontinent cycle,topographic heights of continents typically vary within the interval of about±1.5 km,relatively to a mean value.Topographic maxima of orogenic formations to about 2-4 km are detected along continent-to-continent collisions as well as when adjacent subduction zones interact with continental margins.

Stress states of thawed soil subgrade

The article presents the field measurement results of the stress states of roadbed thawed soil subgrade during the passage of trains. The dependences of the vertical and horizontal stresses on the velocity of the rolling stock motion, the axle load, and the distance from the sleeper sole have been obtained.

Hydraulic Fracturing Stress Measurements and Their Implication for the Design of the First Concrete Lined Pressure Shaft in Sri Lanka

The Urea Oya Multipurpose Development Project （UMDP） is a water transfer, hydropower and irrigation in the south-eastern part of the central highland region of Sri Lanka. During geotechnical site investigation program 42 hydraulic fracture tests and 42 impression packer tests were carried out in 3 boreholes to about 840 m depth. Based on the stress measurements the minimum and maximum horizontal stress ratios were calculated. In situ stress computations at all the tests were based on the assumption that the principal stress components were vertical （σv） and horizontal （σH and σh, the maximum and minimum, respectively）. The results of the measurements had a direct impact on the design of the major openings bearing a high overburden--the underground powerhouse and the transformer cavern--and revealed a significant optimization potential concerning the selection of the lining system of the pressure shaft.

Estimation of the horizontal in-situ stress magnitude and azimuth using previous drilling data

Oil exploration and production,well stability,sand production,geothermal drilling,waste-water or CO_(2) sequestration,geohazards assessment,and EOR processes such as hydraulic fracturing,require adequate information about in-situ stresses.There are several methods for analyzing the magnitude and direction of in-situ stresses.The evaluation of tensile fractures and shear fractures in vertical oil and gas wellbores using image logs is one of these methods.Furthermore,when image logs are run in boreholes,they can be extremely costly and possibly stop the drilling.The data for this study were gathered from seven directional wells drilled into a strike-slip fault reservoir in southern Iran.Vertical stress,minimum horizontal stress,pore pressure,Poisson's ratio of formations,and 233 mud loss reporting points make up the entire data.This is the first time maximum horizontal stress direction has been calculated without referring to image log data.In addition,the points of lost circulation were categorized into natural and induced fracture.The results revealed that,the maximum horizontal stress direction of the reservoir was calculated at 65northeast-southwest.The error rate is roughly 10when comparing the results of this investigation to those obtained from the image log.The maximum horizontal stress direction is calculated precisely.In terms of tensile fracture pressure,the in-situ stress ratio identifies the safest as well as the most critical inclination and azimuth for each well.

Numerical simulation study on hard-thick roof inducing rock burst in coal mine

In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words

Numerical simulation of hydraulic fracture propagation in tight oil reservoirs by volumetric fracturing

Volumetric fracturing is a primary stimulation technology for economical and effective exploitation of tight oil reservoirs. The main mechanism is to connect natural fractures to generate a fracture network system which can enhance the stimulated reservoir volume. By using the combined finite and discrete element method, a model was built to describe hydraulic fracture propagation in tight oil reservoirs. Considering the effect of horizontal stress difference, number and spacing of perforation clus- ters, injection rate, and the density of natural fractures on fracture propagation, we used this model to simulate the fracture propagation in a tight formation of a certain oil- field. Simulation results show that when the horizontal stress difference is lower than 5 MPa, it is beneficial to form a complex fracture network system. If the horizontal stress difference is higher than 6 MPa, it is easy to form a planar fracture system; with high horizontal stress differ- ence, increasing the number of perforation clusters is beneficial to open and connect more natural fractures, and to improve the complexity of fracture network and the stimulated reservoir volume （SRV）. As the injection rate increases, the effect of volumetric fracturing may be improved; the density of natural fractures may only have a great influence on the effect of volume stimulation in a low horizontal stress difference.

LiDAR mapping of ground damage in a heading re-orientation case study

The Subtropolis Mine is a room-and-pillar mine extracting the Vanport limestone near Petersburg,Ohio,U.S.In February of 2018,mine management began implementing a heading re-orientation to better control the negative effects of excessive levels of horizontal stress.The conditions in the headings improved,but as expected,stress-related damage concentrated within crosscuts.The mine operator has worked to lessen the impact of the instabilities in the outby crosscuts by implementing several engineering controls.With the implementation of each control,conditions were monitored and analyzed using observational and measurement techniques including 3D LiDAR surveys.Since the heading re-orientation,several 3D LiDAR surveys have been conducted and analyzed by researchers from the National Institute for Occupational Safety and Health(NIOSH).This study examines(1)the characteristics of each 3D LiDAR survey,(2)the change in the detailed strata conditions in response to stress concentrations,and(3)the change detection techniques between 3D LiDAR surveys to assess entry stability.Ultimately,the 3D LiDAR surveys proved to be a useful tool for characterizing ground instability and assessing the effectiveness of the engineering controls used in the heading re-orientation at the Subtropolis Mine.

Coal pillar mechanics of violent failure in U.S. Mines

This paper seeks to enhance the understanding that the horizontal stresses build up and release during coal pillar loading and unloading(post-failure) drawing upon three decades of observations, geomechanical monitoring and numerical modeling in bump-prone U.S. mines. The focus is on induced horizontal stress in mine pillars and surrounding strata as highly stressed pillars punch into the roof and floor, causing shear failure and buckling of strata; under stiff stratigraphic units of some western US mines, these events could be accompanied by violent failure of pillar cores. Pillar punching eventually results in tensile stresses at the base of the pillar, facilitating transition into the post-failure regime; this transition will be nonviolent if certain conditions are met, notably the presence of interbedded mudstones with low shear strength properties and proper mine designs for controlling seismicity and dynamic loads. The study clearly shows high confining stress build-up in coal pillars resulting in up to twice higher peak vertical stress and high strain energy accumulations in some western US mines in comparison with peak stresses predicted using common empirical pillar design methods. It is the unstable release of this strain energy that can cause significant damage resulting from pillar dilation and ground movements. These forces are much greater than the capacity of most common internal support systems, resulting in horizontal stressinduced roof falls locally, in mines under unremarkable far-field horizontal stress. Attention should be placed on pillar designs as increasing support density may prove to be ineffective. This mechanism is analyzed using field measurements and generic finite-difference stress analyses. The study confirms the higher load carrying capacity of confinement-controlled coal seams in comparison with structurally controlled coal seams. Such significant differences in confining stresses are not taken into account when estimating peak pillar strength using most common empirical techniques such as those proposed by Bieniawski and Salamon. While using lower pillar strength estimates may be considered conservative,it underestimates the actual capacity of pillars in accumulating much higher stress and strain energies,misleading the designer and inadvertently diminishing mine safety. The role of induced horizontal stress in mine pillars and surrounding strata is emphasized in coal pillar mechanics of violent failure. The triggering mechanism for the violent events is sudden loss of pillar confinement due to dynamic loading resulting from failure of overlying stiff and strong strata. Evidence of such mechanism is noted in the field by observed red-dust at the coal-rock interfaces at the location of coal bumps and irregular, periodic caving in room-and-pillar mines quantified through direct pressure measurements in the gob.

Mechanical Properties of Asphalt Pavement Structure in Highway Tunnel

A linear full 3D finite element method (FEM) was performed in order to present the key design parameters of highway tunnel asphalt pavement under double-wheel load on rectangular loaded area considering horizontal contact stress induced by the acceleration/deceleration of vehicles.The key design parameters are the maximum horizontal tensile stresses at the surface of the asphalt layer,the maximum horizontal tensile stresses at the bottom of the asphalt layer and the maximum vertical shear stresses at the surface of the as- phalt layer were calculated.The influencing factors such as double-wheel weight;asphalt layer thickness;base course stiffness modulus and thickness;and the contact conditions among the structure layers on these key design parameters were also examined separately to propose construction procedures of highway tunnel asphalt pavement.

Modeling the Contemporary Stress Field and Deformation Pattern of Eastern Mediterranean

The contemporary stress field in the earth＇s crust is important and provides insights into mechanisms that drive plate motions. In this study, elastic plane stress finite element modeling incorporating realistic rock parameters was used to calculate the stress field, displacement field, and defor- mation of the plate interactions in the eastern Mediterranean. Modeled stress data for the African- Arabian-Anatolian plate interactions with fixed European platform correlate well with observed contemporary stress indicator from the world stress map （WSM） and focal mechanism of earthquakes; while displacement field agrees qualitatively well with GPS vectors and sense of motion indicated by focal mechanisms for large crustal earthquakes （Ms〉6） and plate motion models. Modeling result shows the direction of maximum horizontal compressive stress （σHmax） toward the direction of absolute motion of these plates. Large perturbations in σHmax orientations are shown to occur in and around tectonic boundaries between those plates. It is observed that, although the African plate acts mostly as indenter, which transmits the collisional motion from the Arabian plate to the Anatolian plate, in the current situation, the far-field stress, probably from the subduction in Aegean Arc, is needed to satisfy the contemporary stress field in Anatolia.

Statistical distribution of geomechanical properties and‘Sweet Spots’identification in part of the upper Bakken

Completions and Reservoir Quality are two key attributes that are used to characterize nonconventional hydrocarbon assets.This is because,for optimum exploitation of these unconventional assets,horizontal wells need to be drilled in“Sweet Spots”(i.e.,regions where Completions and Reservoir Quality are both superior).One way to quantify these qualities is to use reservoir and geomechanical properties.These properties can be estimated on a location basis from well logs,and then mapped over terrain using geostatistical modeling.This study presents a‘Sweet Spots’identification workflow based on three performance indexes(Storage Potential Index,Brittleness Index,and Horizontal Stress Index)that can be used to quantify CQ and RQ.The performance indexes are computed from petrophysical property volumes(of Young's Modulus,Bulk Modulus,Shear Modulus,Poisson's Ratio,Minimum Horizontal Stress,Volume of Shale,Total Organic Carbon,Thickness,and Porosity)which are in turn computed from well logs and geostatistical simulation.In the end,the study offers a method to compare the predicted“Sweet Spots”against available production data via their correlation coefficient.The resulting reasonable formation property maps,the successful identification of‘Sweet Spots’,and a correlation coefficient of 0.88(between the predicted“Sweet Spots”and well production data)point to the potential of the proposed effort.