Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

基于卫星遥感影像和DEM融合的地质断裂带研究

应用遥感影像进行地表地质分析,具有数据量大、耗时短、费用低廉等优势,可以快速获取大范围区域的相关地质信息,因此,近年来,卫星遥感影像广泛应用于地质分析工作中。但是由于卫星影像在不同时期内的特征不同,且有大量不同的数据源,在不同影像上获取地质特征困难,因此,如何利用影像处技术增强影像中的地质特征,成为遥感影像地质判读的研究重点。地质断裂带调查一般借助各种直接与间接的调查方法以确认断裂带的存在与否和位置,了解断裂带的活动性,并提供点状或剖面数据材料。卫星影像在追踪研究贯穿地表且可为影像所解析的断裂带上,提供真实的影像证据,有助于直接分析影像中断裂带走向和延展。利用遥感影像处理计算,将卫星影像数据叠加到DEM上编制三维地形图,并利用计算机仿真,通过不同方向的投射光源,找出最适合的观察角及最能凸显出重点区域线状地类的状态。对于地质现象的分析及解释上,可提供较平面影像细致、准确的参考依据。     

Weakness and mechanical anisotropy of phyllosilicate-rich cataclasites developed after mylonites of a low-angle normal fault (Simplon Line,Western Alps)

The Simplon Fault Zone is a late-collisional low-angle normal fault (LANF) of the Western Alps. The hanging wall shows evidence of brittle deformation only, while the footwall is characterized by a c. 1 km-thick shear zone (the Simplon Fault Zone), which continuously evolved, during exhumation and cooling, from amphibolite facies conditions to brittle-cataclastic deformations. Due to progressive localization of the active section of the shear zone, the thermal-rheological evolution of the footwall resulted in a layered structure, with higher temperature mylonites preserved at the periphery of the shear zone, and cataclasites occurring at the core (indicated as the Simplon Line). In order to investigate the weakness of the Simplon Line, we studied the evolution of brittle/cataclastic fault rocks, from nucleation to the most mature ones. Cataclasites are superposed on greenschist facies mylonites, and their nucleation can be studied at the periphery of the brittle fault zone. This is characterized by fractures, micro-faults and foliated ultracataclasite seams that develop along the mylonitic SCC′ fabric, exploiting the weak phases mainly represented by muscovite and chlorite. Approaching the fault core, both the thickness and frequency of cataclasite horizons increase, and, as their thickness increases, they become less and less foliated. The fault core itself is represented by a thicker non-foliated cataclasite horizon. No Andersonian faults or fractures can be found in the footwall damage zone and core zone, whilst they are present in the hanging wall and in the footwall further from the fault. Applying a stress model based on slip tendency, we have been able to calculate that the friction coefficient of the Simplon Line cataclasites was <0.25, hence this fault zone is absolutely weak. In contrast with other fault zones, the weakening effect of fluids was of secondary importance, since they accessed the fault zone only after an interconnected fracture network developed exploiting the cataclasite network.     

Smoothing and re-roughening processes: The geometric evolution of a single fault zone

The geometry of a fault zone exerts a major control on earthquake rupture processes and source parameters. Observations previously compiled from multiple faults suggest that fault surface shape evolves with displacement, but the specific processes driving the evolution of fault geometry within a single fault zone are not well understood. Here, we characterize the deformation history and geometry of an extraordinarily well-exposed fault using maps of cross-sectional exposures constructed with the Structure from Motion photogrammetric method. The La Quinta Fault, located in southern California, experienced at least three phases of deformation. Multiple layers of ultracataclasite formed during the most recent phase. Crosscutting relations between the layers define the evolution of the structures and demonstrate that new layers formed successively during the deformation history. Wear processes such as grain plucking from one layer into a younger layer and truncation of asperities at layer edges indicate that the layers were slip zones and the contacts between them slip surfaces. Slip surfaces that were not reactivated or modified after they were abandoned exhibit self-affine geometry, preserving the fault roughness from different stages of faulting. Roughness varies little between surfaces, except the last slip zone to form in the fault, which is the smoothest. This layer contains a distinct mineral assemblage, indicating that the composition of the fault rock exerts a control on roughness. In contrast, the similar roughness of the older slip zones, which have comparable mineralogy but clearly crosscut one another, suggests that as the fault matured the roughness of the active slip surface stayed approximately constant. Wear processes affected these layers, so for roughness to stay constant the roughening and smoothing effects of fault slip must have been approximately balanced. These observations suggest fault surface evolution occurs by nucleation of new surfaces and wear by competing smoothing and re-roughening processes.