Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

DICE:An open-source MATLAB application for quantification and parametrization of digital outcrop model-based fracture datasets

An open-source MATLAB application(app)named Discontinuity Intensity Calculator and Estimator(DICE)was developed in order to quantitatively characterize the fractures,or in more general,discontinuities within a rocky outcrop in three-dimensional(3D)digital data,such as digital outcrop model(DOM).The workflow proposed for the parametrization of the discontinuities consists of the following steps:(1)Analysis and mapping of the fractures detected within the 3D DOMs;(2)Calculation of the orientation,position and dimensions of discontinuities that are represented by best-fit circular planes;(3)Determining the discontinuity parameters(dimension,distribution,spacing and intensity)by the DICE algorithm using different 3D oriented sampling techniques(3D oriented scanline,3D oriented circular scan window and spherical scan volume).Different sampling methods were bench tested with a synthetic,as well as a natural case study,and compared in order to understand the advantages and limitations of each technique.The 3D oriented circular scan window appears to be the most effective method for fracture intensity estimation with high accuracy(error 0.4%)and stability with variations in scan radius.

Stress tensors and quantification of fracture patterns to analyze connectivity and potential fluid flow in a mesa landform of the Northern Andes

This work applies stress tensors inversions and quantification of fracture patterns along the Mesa de Los Santos, in the Eastern Cordillera of Colombia, to better understand the potential fluid flow. It thus contributes to the conceptual hydrogeological model. The area was subdivided into three blocks, separated by the NW Potreros and the Los Santos faults, having minor inner faults of different orientations. This separation facilitates the analysis of the fractures measured in the field, which in general show high dip angles and a conjugate geometry in the northern block, tension fractures(Mode I) in the central block, and a random distribution in the southwestern block. WinTensor treatment of slickensides yielded a maximum horizontal stress(SHmax) of 111o, which coincides with the WNW-ESE tensor observed from the conjugate and tension joints. We then used Frac Pa Q to generate interpolation maps of fracture intensity and density. The maps show the highest values in the central block and the lowest in the northern block,where the precipitation is higher, causing intensive rock weathering and homogenization of the fracture planes. Although the highest values of connectivity by line are found to the south of the mesa, we suggest the possibility of greater flow from the recharge zone(NE) along bedding planes and open NW-SE fractures.

A Gel-Based Solidification Technology for Large Fracture Plugging

Fault fractures usually have large openings and considerable extension. Accordingly, cross-linked gel materials aregenerally considered more suitable plugging agents than water-based gels because the latter often undergo contaminationvia formation water, which prevents them from being effective over long times. Hence, in this study, aset of oil-based composite gels based on waste grease and epoxy resin has been developed. These materials havebeen observed to possess high compressive strength and resistance to the aforementioned contamination, therebyleading to notable increase in plugging success rate. The compressive strength, thickening time, and resistance toformation water pollution of these gels have been evaluated indoors. The results show that the compressivestrength of the gel can reach 11 MPa;additionally, the related gelation time can be controlled to be more than3 h, thereby providing a safe construction time;Invasion of formation water has a small effect on the gel strengthand does not shorten the thickening time. All considered performance indicators of the oil-based gel confirm itssuitability as a plugging agent for fault fractures.

Theoretical and experimental study of 3-D initial fracture and its significance to faulting

The experimental results of 3-D fracture under compression are introduced in brief and the theory of stress criterion of 3-D fracture is studied. Methods to imitate initial fractures are developed. It is pointed that there are important defects in the extreme value （EV） method ever proposed by Palaniswamy and Knauss. The major defect lies in that only two Euler angles （2EA） are considered, but another one is neglected. If the variation of all the three Euler angles （3EA） are considered, one can get better result which is consistent with the observation of faulting that extends on curved surfaces but not on planes. The method of evaluating maximal normal stress direction vector （NSDV） is proposed and further proved to be equivalent to the 3EA method. It is proved that the NSDV method can be further optimized to the method of composition of the first principal differential plane （CFPDP）. The results from CFPDP method can fit the curved surfaces of initial growth observed in the experiments of 3-D fracture. The CFPDP method can also be used to interpret the 3-D fractures of the slipping section between the asperities in the buried fault plane that is modeled as ellipse crack. The results of 3-D fracture can be applied to interpreting the related problems of faulting including the mechanism of a lot of shatter rocks with different dimensions, the cause of earthquakes occurred at the edge of plate under low shear stress, and the mechanism of anisotropy caused by the extensive dilatancy anisotropy （EDA） cracks.

Tensile Fractures and in situ Stress Measurement Data Constraints on Cretaceous-Present Tectonic Stress Field Evolution of the Tanlu Fault Zone in Shandong Province,North China Craton

Tectonic stress fields are the key drivers of tectonic events and the evolution of regional structures.The tectonic stress field evolution of the Tanlu fault zone in Shandong Province,located in the east of the North China Craton(NCC),may have preserved records of the NCC’s tectonic history.Borehole television survey and hydraulic fracturing were conducted to analyze the paleo and present tectonic stress fields.Three groups of tensile fractures were identified via borehole television,their azimuths being NNW-SSE,NW-SE and NE-SW,representing multiple stages of tectonic events.Hydraulic fracturing data indicates that the study region is experiencing NEE-SWW-oriented compression and nearly-N-Soriented extension,in accordance with strike-slip and compression.Since the Cretaceous,the orientation of the extensional stress has evolved counterclockwise and sequentially from nearly-NW-SE-oriented to NE-SW-oriented and even nearly N-S-oriented,the stress state having transitioned from strike-slip-extension to strike-slip-compression,in association with the rotating and oblique subduction of the Pacific Plate beneath the NCC,with the participation of the Indian Plate.

Quantitative multiparameter prediction of fault-related fractures： a case study of the second member of the Funing Formation in the Jinhu Sag, Subei Basin

In this paper, the analysis of faults with different scales and orientations reveals that the distribution of fractures always develops toward a higher degree of similarity with faults, and a method for calculating the multiscale areal fracture density is proposed using fault-fracture self-similarity theory. Based on the fracture parameters observed in cores and thin sections, the initial apertures of multiscale fractures are determined using the constraint method with a skewed distribution. Through calculations and statistical analyses of in situ stresses in combination with physical experiments on rocks, a numerical geomechanical model of the in situ stress field is established. The fracture opening ability under the in situ stress field is subsequently analyzed. Combining the fracture aperture data and areal fracture density at different scales, a calculation model is proposed for the prediction of multiscale and multiperiod fracture parameters, including the fracture porosity, the magnitude and direction of maximum permeability and the flow conductivity. Finally, based on the relationships among fracture aperture,density, and the relative values of fracture porosity and permeability, a fracture development pattern is determined.

Structural Traces of Secondary Faults(Fractures) along the Main Faults and Their Reliability as Kinematic Indicators

Secondary/minor structures occurring along the main fault surfaces are important indicators for judging the kinematic characteristics of faults.However,many factors can lead to the formation of these structures,which results in the difficulty for rapid judgment and application in the fields.A series of secondary faults/fractures developed due to the movement of main faults are the most important and widespread phenomena in the scope of brittle deformation.The morphology of the main fault surfaces is various,and former researchers mainly discussed the structures on the main even fault surfaces. However,the fluctuation of fault surfaces is the intrinsic character of the faults,and the intersection between the main fault and secondary faults/fractures can produce a series of kinematic indicators on the main fault surfaces.Based on previous studies and our observations,i.e.the structural traces of the P,R,R',T and X shears/faults along the main faults,some indicators which are rarely reported previously,are described in the paper.Furthermore,their reliabilities are also discussed,and more practical and reliable criteria are brought forward.We suggest that the simple application of congruous and incongruous steps without knowing their exact origins should be abandoned in the fields,and several types of indicators along one fault surface should be checked with each other as much as possible.Meanwhile,the origins of some other arcuate indicators on the fault surfaces are also discussed,and new models are brought forward.

Application of the third generation of coherent cube in recognizing faults and fractures

Three-dimensional coherent cube is an extremely effective new technique for interpreting seismic data. It has obvious advantages in many aspects compared with the conventional 3D data volume, such as recognizing faults and fractures, interpreting ancient channels, and edge detection of oil-gas reservoir. Coherent cube is to condense and extract information around a certain point in 3D data volume, and then highlight the original characteristics of the geologic body at this point. Therefore, in terms of its essence, coherent cube is a special seismic attribute cube and those points having rather small coherent value are related to the discontinuity of geologic body. In practical production, people often interpret horizontal slices or layer slices of coherent cube, and this provides advantageous foundations for resolving special problems in oil-gas exploration.……

On the Orientation of Fractures with Transpressional and Transtensional Wrenches in Pre-Existing Faults

The orientation of fractures with transpressional and transtensional wrenches in pre-existing faults has not been quantitatively determined. Based on Coulomb failure criterion and Byerlee’s frictional sliding criterion, this paper has indicated quantitative geometric relationships between the pre-existing fault and the local induced principal stress axes caused by the rejuvenation of the pre-existing fault. For a hidden pre-existing fault with some cohesion, the angles between the local induced principal stress axes and the pre-existing fault quantitatively vary with the applied stress and the cohesion coefficient, the ratio of the thickness of the cover layer to the thickness of the whole wrench body, whether transpressional or transtensional wrenches occur. For a surface pre-existing fault with zero cohesion, the angles between the pre-existing fault and the local induced principal stress axes are related to the rock inner frictional angle regardless of both the applied stress and the cohesion coefficient where transpressional wrenches occur, and the local induced maximum principal stress axis is identical with the applied maximum principal stress axis where transtensional wrenches occur. Therefore, the geometric relationships between the pre-existing faults and their related fractures are defined, because the local induced principal stress axes determine the directions of the related fractures. The results can be applied to pre-existing weak fabrics. They can help to understand and analyze wrench structures in outcrops or subsurface areas. They are of significance in petroleum exploration.

Response of underground pipeline through fault fracture zone to random ground motion

It is assumed that a pipeline is laid through a vertical fault fracture zone, and is excited by seismic ground motion modelled as stationary stochastic process. For horizontal incidence of waves, the cross-PSD （Power Spectral Density） function is developed using wave propagation theory, while for vertical incidence of waves the cross-PSD function is composed by auto-PSD model, coherence model and site response model. As the seismic input, the eross-PSD function is used to calculate the the axial and lateral seismic responses of underground pipeline through the fracture zone. The results show that the incident directions of seismic waves, width and soil property of the fracture zone have great influence on underground pipeline. It is suggested that the flexible joints with appropriate stiffness should be added into the pipeline near the interfaces between the fracture zone and the surrounded media.

Geomechanical Modeling of Stress and Fracture Distribution during Contractional Fault-Related Folding

Understanding and predicting the distribution of fractures in the deep tight sandstone reservoir are important for both gas exploration and exploitation activities in Kuqa Depression. We analyzed the characteristics of regional structural evolution and paleotectonic stress setting based on acoustic emission tests and structural feature analysis. Several suites of geomechanical models and experiments were developed to analyze how the geological factors influenced and controlled the development and distribution of fractures during folding. The multilayer model used elasto-plastic finite element method to capture the stress variations and slip along bedding surfaces, and allowed large deformation. The simulated results demonstrate that this novel Quasi-Binary Method coupling composite failure criterion and geomechanical model can effectively quantitatively predict the developed area of fracture parameters in fault-related folds. High-density regions of fractures are mainly located in the fold limbs during initial folding stage, then gradually migrate from forelimb to backlimb, from limbs to hinge, from deep to shallow along with the fold uplift. Among these factors, the fold uplift and slip displacement along fault have the most important influence on distributions of fractures and stress field, meanwhile the lithology and distance to fault have also has certain influences. When the uplift height exceeds approximately 55 percent of the total height of fold the facture density reaches a peak, which conforms to typical top-graben fold type with large amplitude and high-density factures in the top. The overall simulated results match well with core observation and FMI results both in the whole geometry and fracture distribution.

Influences of pre-existing fracture on ground deformation induced by normal faulting in mixed ground conditions

Physical model tests have been conducted by various researchers to investigate fault rupture propagation and ground deformation induced by bedrock faulting. However, the effects of pre-existing fracture on ground deformation are not fully understood. In this work, six centrifuge tests are reported to investigate the influence of pre-existing fracture on ground deformation induced by normal faulting in sand, clay and nine-layered soil with interbedded sand and clay layers. Shear box tests were conducted to develop a filter paper technique, which was adopted in soil model preparation to simulate the effects of pre-existing fracture in centrifuge tests. Centrifuge test results show that ground deformation mechanism in clay, sand and nine-layered soil strata is classified as a stationary zone, a shearing zone and a rigid body zone. Inclination of the strain localization is governed by the dilatancy of soil material. Moreover, the pre-existing fracture provides a preferential path for ground deformation and results in a scarp at the ground surface in sand. On the contrary, fault ruptures are observed at the ground surface in clay and nine-layered soil strata.

Using structure restoration maps to comprehensively identify potential faults and fractures in compressional structures

Faults and fractures of multiple scales are frequently induced and generated in compressional structural system. Comprehensive identification of these potential faults and fractures that cannot be distinguished directly from seismic profile of the complex structures is still an unanswered problem. Based on the compressional structural geometry and kinematics theories as well as the structural interpretation from seismic data, a set of techniques is established for the identification of potential faults and fractures in compressional structures. Firstly, three-dimensional(3D) patterns and characteristics of the faults directly interpreted from seismic profile were illustrated by 3D structural model. Then, the unfolding index maps, the principal structural curvature maps, and tectonic stress field maps were obtained from structural restoration. Moreover, potential faults and fractures in compressional structures were quantitatively identified relying on comprehensive analysis of these three maps. Successful identification of the potential faults and fractures in Mishrif limestone formation and in Asmari dolomite formation of Buzurgan anticline in Iraq demonstrates the applicability and reliability of these techniques.

Fracture Characters and Seismogenic Fault of the Yajiang Earthquake Sequence with M_S6.0 in Sichuan in 2001

The Yajiang earthquake sequence in 2001, with the major events of M S 5.1 on Feb. 14 and of M S 6.0 on Feb.23, are significant events in the Sichuan region during the last 13 years. Eighty-eight earthquakes in the sequence with at least 5 distinct onset parameters for each recorded by the Sichuan Seismic Network in the period of Jan. 1 through June 30, 2001 were chosen for this study. The events are relocated and the focal mechanism is derived from P-wave onsets for 13 events with relatively larger magnitudes. The focal depth of all earthquakes fall between a range of 2km to 16km, with dominant distribution between 9km to 11km. The foreshocks, the M S5.1 earthquake and the M S6.0 earthquake and their aftershocks are all located close to the Zihe fault and the dominant epicentral distribution is in NW direction, identical to that of the fault. The fracture surface of the focal mechanism is determined in accordance to the mass transfer orientation in the recent earth deformation field in the Yajiang region. The P axes of the principal compressive stress in focal mechanism solutions of the 13 events show bigger vertical components, and the horizontal projection trending SE. The earthquakes are of left-lateral, strike-slip normal, and normal strike-slip types. The rupture surface of most earthquakes strike NW-SE, dipping SW. Based on the above information, we conclude that the Zihe fault that crosses the earthquake area, striking NW and dipping SW, is the seismogenic fault for the Yajiang earthquake sequence.

Controls of strike-slip fault on fractures:Insight from 3D discrete element simulation

The fracture-cave reservoirs controlled by strike-slip faults are the main targets for oil and gas exploration of ultradeep carbonate in the Tarim Basin.It is of great significance to clarify the distribution rules of fractures related to strike-slip faults for guiding the exploration and development of ultra-deep oil and gas.In this study,six groups of strike-slip fault 3D models based on discrete element numerical simulation method have been created to investigate characteristics of fault-related fracture development and distribution law.In addition,we compared the modeling results to the measurement of fractures from the outcrop of a strike-slip fault in the Northern Tarim Basin to verify their validity.The results show that the stress environment is stable in the simple strike-slip section,and fractures intersecting with the strike-slip direction at a small angle are the principal fracture.In the releasing stepover and double-bend sections,the maximum principal stress changes from horizontal to vertical during the formation of pull-apart zones,where the principal fractures intersect the strike-slip direction at a large angle.The maximum principal stress in the restraining stepover and double-bend section remains horizontal,but their strikes change significantly with the increment of fault displacement.Thus,fractures intersecting the strike-slip direction at a small angle will become principal fractures early on,while those parallel to or anti-intersecting the strike-slip direction at a small angle will become principal fractures later.There are obvious differences in the development of fractures in different structural positions of strike-slip faults.Fractures are mainly concentrated in the fault tips,connections,and fault plane,and the magnitude of the fault damage zone is relatively larger in the first two.Compared with fault displacement,the principal damage zone(PDZ)shows stronger control on the distribution and development intensity of fractures.With the increment of fault displacement,the width of the fault damage zone and the fracture density first rapidly increases before the formation of PDZ and then slows down.Moreover,the formation time of PDZ in the restraining double-bend and stepover section is earlier than the simple strike-slip,releasing double-bend,and releasing stepover sections,and absorb more strain before the formation of the principal displacement zone.Thus,the restraining sections have the highest fracture intensity,followed by the pull-apart sections,then the simple strikeslip section.The results play an important role in understanding the development law of fractures related to strike-slip faults in different arrangements and move modes.

Integration of Electrical Resistivity and Electromagnetic Radiation Methods for Fracture Flow System Detection

An electrical resistivity and electromagnetic emission survey was carried out involving the use of vertical electrical soundings (VES) and natural pulse electromagnetic field of the earth (NPEMFE). The use of this new methodology managed to detect the fracture flow system rupture zones in the underground, also answered the questions about the deferent subsurface water bodies. The present study focuses on Marsaba-Feshcha sub-basin in the northeast of the Dead Sea. Due to the scarcity of boreholes in the study area, several geophysical methods were implanted. The combination of these two methods (VES and NPEMFE) with the field observations and East-West transversal faults with the coordination (624437/242888) was determined, cutting through the anticlines with their mainly impervious cores with fracture length of >400 m. These transversal faults saddle inside Nabi Musa syncline (Boqea syncline), leading to a hydraulic connection between the Lower and the Upper Aquifer. Due to the identified transversal fault, the water of the Upper and Lower Aquifer mixed and emerged as springs at Ein Feshcha group.

Seismic and Tectonic Correspondence of Major Earthquake Regions in Southern Ghana with Mid-Atlantic Transform-Fracture Zones

For four centuries now, southern Ghana has been known to be seismically active, and there is no clear geological explanation for the cause of the seismicity. By evaluating new field data and information with re-interpreted historical earthquake data of southern Ghana, the nature of the seismicity of southern Ghana has been elucidated. The mutual connection between the earthquake epicentres and the remote causes by Mid-Atlantic transform faults and fracture zones has been established. The seismic regions of southern Ghana have been linked separately to tectonic faults and activities of the St. Paul’s and Romanche transform-fracture zone systems offshore in the Gulf of Guinea to onshore. It is concluded that the seismicity of southern Ghana is due to tectonic activities of the St. Paul’s and Romanche transform-fracture systems. The Accra region earthquakes originate from reactivation of faults in the Romanche transform-fracture zone, and propagate onshore through Accra and environs. The Axim region earthquakes come from reactivated faults linked to the St Paul’s fracture zone, which go through southern Cote D’Ivoire to Ghana. Seismotectonic movements along the St Paul’s transform and fracture zones have quieted since 1879. But movement along the Romanche Transform fault and Fracture zone is active, causing ongoing seismicity of southern Ghana.

Failure responses of rock tunnel faces during excavation through the fault-fracture zone

It is essential to cast light on the construction risks in tunnel excavations through the fault-fracture zone(FFZ).This study adopts the material point method(MPM)to simulate the failure responses of a rock tunnel face during excavation through the FFZ.A numerical study was conducted to compare a physical model test and validate the feasibility of using the MPM in simulating tunnel face failure.One hundred ninety numerical simulation cases were constructed to represent a rock tunnel excavation project with different site con-figurations.The simulation results suggest that the cohesion and the friction angle significantly influence failure responses.The tunnel cover depth can magnify the failure responses,and the FFZ thickness significantly affects the mobilized rock mass volume when the FFZ consists of a weak rock mass.The numerical simulation results suggest three deformation patterns:face bulge,partial failure,and slide collapse.The failure responses can be characterized by stress arch,slip surface,angle of reposing,and influence range.The insights suggested by the face failure responses during excavation through the FFZ can aid field engineers in determining the scope of possible damage,and in establishing emergency measures to minimize losses if such failure occurs.