The numerical simulation and inversion fitting of radon concentration distribution in homogeneous overburden above active fault zones

Based on the convection and diffusion mechanisms of radon migration, in this paper we deduce the two-dimensional differential equation for radon transportation in the overburden above active fault zones with an unlimited extension along the strike. Making use of the finite difference method, the radon concentration distribution in the overburden above active faults is calculated and modeled. The active fault zone parameters, such as the depth and the width of the fault zone, and the value of radon concentration, can be inverted from the measured radon concentration curve. These realize quantitative interpretation for radon concentration anomalies. The inversion results are in good agreement with the actual fault zone parameters.

Dominant pulse simulation of near fault ground motions

In this study, a new mathematical model is developed composed of two parts, including harmonic and polynomial expressions for simulating the dominant velocity pulse of near fault ground motions. Based on a proposed velocity function, the corresponding expressions for the ground acceleration and displacement time histories are also derived. The proposed model is then fitted using some selected pulse-like near fault ground motions in the Next Generation Attenuation （NGA） project library. The new model is not only simple in form but also simulates the long-period portion of actual velocity near fault records with a high level of precision. It is shown that the proposed model-based elastic response spectra are compatible with the near fault records in the neighborhood of the prevailing frequency of the pulse. The results indicate that the proposed model adequately simulates the components of the time histories. Finally, the energy of the proposed pulse was compared with the energy of the actual record to confirm the compatibility.

Numerical simulation analysis of covering rock strata as mining steep-inclined coal seam under fault movement

The fault is one important factor for the stability of overburden strata caused by steeply inclined coal seam. The stress and displacement change of overburden strata caused by steeply-inclined coal seam mining activity under faulting was simulated by FLAC2D finite differential program on the basis of Zhaogezhuang mining example belonging to Kailuan Mining Group. From the results, the stress and displacement clouding image after mining became complex because of the fault, that is, a kind of weak structural plane. The stress concentration region concentrated around the goaf, and also around the fault plane. As the mining depth increases, the stress and displacement within the fault zone change significantly. This movement and deformation characteristic of overburden strata can provide theoretical basis for the similar mining condition.

Numerical Simulation of Coal Floor Fault Activation Influenced by Mining

By means of the numerical simulation software ANSYS, the activation regularity of coal floor faults caused by mining is simulated. The results indicate that the variation in horizontal, vertical and shear stresses, as well as the horizontal and vertical displacements in the upper and the lower fault blocks at the workface are almost identical. Influ- enced by mining of the floor rock, there are stress releasing and stress rising areas at the upper part and at the footwall of the fault. The distribution of stress is influenced by the fault so that the stress isolines are staggered by the fault face and the stress is focused on the rock seam around the two ends of the fault. But the influence in fault activation on the upper or the lower fault blocks of the workface is markedly different. When the workface is on the footwall of the fault, there is a horizontal tension stress area on the upper part of the fault; when the workface is on the upper part of the fault, it has a horizontal compressive stress area on the lower fault block. When the workface is at the lower fault block, the maximum vertical displacement is 5 times larger then when the workface is on the upper fault block, which greatly in- creases the chance of a fatal inrush of water from the coal floor.

Numerical simulation of fault activity owing to hydraulic fracturing

We built a three-dimensional model to simulate the disturbance of the stress field near the reverse fault in Zhaziao, Leyi Township owing to hydraulic fracturing. The pore pressure, and shear and normal stresses during fracturing are analyzed in detail. Input rock mechanics parameters are taken from laboratory test data of shale samples from the study area. The simulation results suggest that after 16 hours of fluid injection, the pore-pressure variation can activate the reverse fault, i.e., we observe reverse slip, and the shear stress and displacement on the fault plane increase with time. The biggest stress–strain change occurs after one hour of fluid injection and the yield point appears about 0.5 h after injection. To observe the stress evolution in each section, the normal displacement on the boundary is constrained and the fault plane is set as nonpermeable. Thus, the sliding is limited and the shear displacement is only in the scale of millimeters, and the calculated magnitude of the induced earthquakes is between Mw-3.5 and Mw-0.2. The simulation results suggest that fluid water injection results in inhomogeneous fracturing. The main ruptured areas are around the injection positions, whereas the extent of rupturing and cracks in other areas are relatively small. Nevertheless, nonnegligible fault activation is recorded. Sensitivity analysis of the key parameters suggests that the pore pressure is most sensitive to the maximum unbalanced force and the internal friction angle strongly affects the fault slip. Finally, the comparison between the effective normal stress and the maximum and minimum principal stresses on the fault plane explains the fault instability, i.e., the Mohr circle moves towards the left with decreasing radius reduces and intersects the critical slip envelope, and causes the fault to slip.

3D simulation of near-fault strong ground motion: comparison between surface rupture fault and buried fault

In this paper, near-fault strong ground motions caused by a surface rupture fault （SRF） and a buried fault （BF） are numerically simulated and compared by using a time-space-decoupled, explicit finite element method combined with a multi-transmitting formula （MTF） of an artificial boundary. Prior to the comparison, verification of the explicit element method and the MTF is conducted. The comparison results show that the final dislocation of the SRF is larger than the BF for the same stress drop on the fault plane. The maximum final dislocation occurs on the fault upper line for the SRF; however, for the BE the maximum final dislocation is located on the fault central part. Meanwhile, the PGA, PGV and PGD of long period ground motions （≤ 1 Hz） generated by the SRF are much higher than those of the BF in the near-fault region. The peak value of the velocity pulse generated by the SRF is also higher than the BE Furthermore, it is found that in a very narrow region along the fault trace, ground motions caused by the SRF are much higher than by the BF. These results may explain why SRFs almost always cause heavy damage in near-fault regions compared to buried faults.

Numerical simulation on fault water-inrush based on fluid-solid coupling theory

About 75% water-inrush accidents in China are caused by geological structure such as faults, therefore, it is necessary to investigate the water-inrush mechanism of faults to provide references for the mining activity above confined water. In this paper, based on the fluid-solid coupling theory, we built the stress-seepage coupling model for rock, then we combined with an example of water-inrush caused by fault, studied the water-inrush mechanism by using the numerical software COMSOL Mutiphysics, analyzed the change rule of shear stress, vertical stress, plastic area and water pressure for stope with a fault, and estimated the water-inrush risk at the different distances between working faces and the fault. The numerical simula- tion results indicate that： （1） the water-inrush risk will grow as the decrease of the distance between working face and the fault; （2） the failure mode of the rock in floor with fault is shear failure; （3） the rock between water-containing fault and working face failure is the reason for water-inrush.

Geomechanics model test and numerical simulation of 2G-NPR bolt support effect in an active fault tunnel

Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.

Quantification and assessment of fault uncertainty and risk using stochastic conditional simulations

The effect of geological uncertainty on the development and mining of underground coal deposits is a key issue for longwall mining, as the presence of faults generates substantial monetary losses. This paper develops a method for the conditional simulation of fault systems and uses the method to quantify and assess fault uncertainty. The method is based on the statistical modelling of fault attributes and the simulation of the locations of the centres of the fault traces. Fault locations are generated from the thinning of a Poisson process using a spatially correlated probability field. The proposed algorithm for simulating fault traces takes into account soft data such as geological interpretations and geomechanical data. The simulations generate realisations of fault populations that reproduce observed faults, honour the statistics of the fault attributes, and respect the constraints of soft data, providing the means to thereby model and assess the related fault uncertainty.

Numerical Simulation of Earthquake Nucleation Process and Seismic Precursors on Faults

To understand precursory phenomena before seismic fault sliP,this work focuses onearthquake nucleation process on a fault plane through numerical simulation.Rate and statedependent friction law with variable normal stress is employed in the analysis.The resultsshow that in the late stage of nucleation process:(1)The maximum slip velocity ismonotonically accelerating;(2)The slipping hot spot(where the slip rate is maximum)migrates spontaneously from a certain instant,and such migration is spatially continuous;(3)When the maximum velocity reaches a detectable order of magnitude(at least one orderof magnitude greater than the loading rate),the remaining time is 20 hours or longer,andthe temporal variation of slip velocity beyond this point may be used as a precursoryindicator;(4)The average slip velocity is related to the remaining time by a log-log linearrelation,which means that a similar relation between rate of microseismicity and remainingtime may also exist;(5)when normal stress variation is taken

Measuring the Qatar-Kazeron Fault Dip Using Random Finite Fault Simulation of September 27, 2010 Kazeron Earthquake and Analytical Signal Map of Satellite Magnetic Data

In this research the fault parameters causing the September 27, 2010 Kazeron Earthquake with a magnitude of MW = 5.8 (BHRC) were determined using the random finite fault method. The parameters were recorded by 27 accelerometer stations. Simulation of strong ground motion is very useful for areas about which little information and data are available. Considering the distribution of earthquake records and the existing relationships, for the fault plane causing the September 27, 2010 Kazeron Earthquake the length of the fault along the strike direction and the width of the fault along the dip direction were determined to be 10 km and 7 km, respectively. Moreover, 10 elements were assumed along the length and 7 were assumed along the width of the plane. Research results indicated that the epicenter of the earthquake had a geographic coordination of 29.88N - 51.77E, which complied with the results reported by the Institute of Geophysics Tehran University (IGTU). In addition, the strike and dip measured for the fault causing the Kazeron Earthquake were 27 and 50 degrees, respectively. Therefore, the causing fault was almost parallel to and coincident with the fault. There are magnetic discontinuities on the analytical signal map with a north-south strike followed by a northwest-southeast strike. The discontinuities are consistent with the trend of Kazeron fault but are several kilometers away from it. Therefore, they show the fault depth at a distance of 12 km from the fault surface.

Fractal Features and Computer Simulation of Fan-Shaped Sedimentary Bodies in Faulted Basins

Ｆａｎ－ｓｈａｐｅｄｓｅｄｉｍｅｎｔａｒｙｂｏｄｉｅｓａｒｅｔｈｏｓｅｍａｉｎｌｙｃｏｍ－ｐｏｓｅｄｏｆｔｅｒｉｇｅｎｏｕｓｃｏａｒｓｅｄｅｐｏｓｉｔｓａｎｄｂｅｉｎｇｆａｎ－ｓｈａｐｅｄｉｎｍｏｒｐｈｏｌｏｇｙ，ｉｎｃｌｕｄｉｎｇａｌｕｖｉａｌｆａ...

Geological Characteristics and Numerical Simulation of Badong Fault in TGP Reservoir Area

The new county-seat town of Badong in the reservoir area of the Three Gorges Project is located on a huge arcuate slope with a convex bank toward north. The slope is cut by a fault, Badong fault, trending in east-west in its back part. It is concerned if the huge arcuate slope is related to mass rock creep, and what is the role of the Badong fault in the formation of the huge arc slope? The Badong fault was put into main consideration in this paper. The data from field investigation were reviewed. Three main features of the Badong fault were summarized: a bedding fault between the Jialingjiang formation (T\-1j) and Badong formation (T\-2b), breccias with compound component, and multiple stages of activity. It was proposed that most of the breccias were formed by fracture-filling. To understand the state of stress and behavior of deformation of the fault during the incision of the Yangtze River as well as the initiation and development of the slope, numerical simulation was conducted. Results indicate that there is a tensional stress zone in the upper part of the fault, and that activity of the fault is dominated by bedding sliding. Opening was also noted in the upper part of the fault in the late periods. The results are consistent with the field observation of the fault. The displacement in the slope is small, which makes us conclude that there is no certain relation between the formation of the arcuate slope and the Badong fault.

Construction of Network Fault Simulation Platform and Event Samples Acquisition Techniques for Event Correlation

Event correlation is one key technique in network fault management. For the event sample acquisition problem in event correlation, a novel approach is proposed to collect the samples by constructing network simulation platform. The platform designed can set kinds of network faults according to user's demand and generate a lot of network fault events, which will benefit the research on efficient event correlation techniques.

Numerical simulation on the influences of Wenchuan earthquake on the surrounding faults

On 12 May 2008, the devastating Wenchuan earthquake struck the Longmenshan fault zone, which comprised the eastern margin of the Tibetan Plateau, and this fault zone was predominantly a convergent boundary with a right-lateral strike-slip component. After such a large-magnitude earthquake, it was crucial to analyze the influences of the earthquake on the surrounding faults and the potential seismic activity. In this paper, a complex viscoelastic model of western Sichuan and eastern Tibet regions was constructed including the topography. Based on the findings of co-seismic static slip distribution, we calculated the stress change caused by the Wenchuan earthquake with the post-seismic relaxation into consideration. Our preliminary results indicated that： （1） The tectonic stressing rate was relatively high in Kunlun mountain pass-Jiangcuo, Ganzi-Yushu, Xianshuihe and Zemuhe faults; while in the east Kunlun and Longriba was medium; also the value was less in the Minjiang, Longmenshan, Anninghe and Huya faults. As to the Longmenshan fault, the value was 0.28×10-3 MPa/a to 0.35×10-3 MPa/a, which is coincident with the previous long recurrence interval of Wenchuan earthquake; （2） The Wenchuan earthquake not only caused the Coulomb stress decrease in the source region, but also the stress increase in the two terminals, especially the northeastern segment, which is comparatively consistent with the aftershock distribution. Meanwhile, the high concentration areas of the static slip distribution were corresponding to the Coulomb stress reductions; （3） The Coulomb stress change caused by Wenchuan earthquake showed significant increase on five major faults, which were northwestern segment of Xianshuihe fault, eastern Kunlun fault, Longriba fault, Minjiang fault and Huya fault respectively; also the Coulomb stress on the fault plane of the Yushu earthquake was faintly increased; （4） We defined the recurrence interval as the time needed to accumulate the magnitude of the stress drop, and the recurrence interval of Wenchuan earthquake was estimated about 1 714 a to 2 143 a correspondingly.

Magnetic Force Simulation of Cables in Microgrid during Faults

In recent years, more and more electric utilities are using underground cables to distribute electric power rather than overhead transmission line. However, the cost of installation and maintenance of underground cables is very expensive. Thus, the proper design and damage prediction of cables are crucial. This paper is focused on the magnetic force waveforms simulation of cables under different types of faults using PSCAD and COMSOL. The results show that three-phase fault leads to the largest magnetic forces and the maximum magnitude of the forces in the x-direction is about 2.5 N. Also, the magnetic fields surrounding the cables are different depending on the arrangements of cables buried method. Although the magnitude is small, considering the long distance and long operating time of underground cables, the forces between cables can cause failures under some conditions. In the future, more types of faults such as high impedance fault and different protecting technologies can be studied.

Numerical simulation on the movement law of overlying strata in the stope with a fault and analysis of its influence on the ground gas drainage boreholes

In order to study the influence of a fault on the movement law of the overlying strata as well as its effect on the gas drainage boreholes, based on the practical situation of 1242（1） panel at Xieqiao Mine in Huainan, the Finite Element Method （FEM） model was built up, and the distribution of the stress field and the displacement field of the overlying strata in the stope with a fault were simulated by using the FEM software ANSYS. The results indicate that because of the existence of the fault, the horizontal displacement of overlying strata near the gas drainage borehole becomes larger than that in the stope without a fault, and the distribution of the stress field of the overlying strata changes greatly. When the working face is far away from the fault, the distribution of the stress field is approximately symmetrical. As the working face advances to the place 50 m away from the fault, the stress range at the right side goaf area is as twice as that at the left side. Here, the stress distribution area of goaf area and the fault plane run through, the fracture-connected-zone is formed. It can be presumed that the gas adsorbed in the coal and rock will flow into the fault zone along the fracture-connected-zone, which causes the quantity of gas drainage reduce remarkably.

Numerical simulation of the floor water-inrush in working face influenced by fault structure

Used numerical simulation method to study the floor water-inrush mechanism in working face which was influenced by fault structure, and set up many kinds of models and performs numerical calculation by fully using large finite element soft-ANSYS and element birth-death method. The results show that the more high the underground water pressure, the more big the floor displacement and possibility of water-inrush; the floor which has fault structure is more prone to water-inrush than the floor which not has fault structure, the floor which has multi-groups cracks is more prone to water-inrush than the floor which has single-group cracks. The numerical simulation result forecasts the water-inrush in working face preferably.

Fault simulation in ultrahigh voltage substation

Ultrahigh voltage(UHV)and extra-high voltage(EHV)have been widely used in power system,so the requirement for technology of operation management,as well as the accident analysis and disposal are more pressing.Based on the theory of Newton-Laphson power flow method,this simulation system imitates the primary device and secondary circuit,together with the normal,the abnormal and accidents in operation.Through simulation in various situations,it validates that this design is capable of simulating the complicated faults in UHV and EHV accurately and in real time.In addition,it can analyse and dispose them efficiently.

SEM Numerical Simulation of Vertical and Inclined Fault Zone Trapped Waves and Comparison of Their Wave Fields

Fault zone trapped waves （FZTWs） mainly travel along the fractured fault zone （FZ） which is of low velocity and high attenuation. FZTWs often carry significant information about a fault＇s internal structure, so it is important to understand their wave field characteristics for FZ structure inversion. Most previous simulations are based on vertical faults, while in this paper we implement the FZTW simulations on vertical or inclined faults and compare their wave fields in both time and frequency domains. The results show that the existence of fault zone and inclined angle of fault can significantly influence the features of waves near faults. In amplitude, a fault zone can generate a larger amplitude of waves. The velocity contrast between two wails of fault may lead to amplification of amplitudes in the low velocity fault wall. In frequency, a fault zone tends to influence the waves in the low frequency range. In a pattern of particle polarization of FZTWs, it tends to be single direction for vertical faults but fork to multiple directions for inclined faults, which might provide a new way to study the fault zone with FZTWs. These conclusions may be valuable for FZ structure inversion, and will enhance the knowledge on near-fault strong ground motions.