Granite deformation and behavior of acoustic emission sequence under the temperature and pressure condition at different crust depths

IntroductionFocal depths of shallow strong intraplate earthquakes are mostly distributed in the highstrength range of lithosphere controlled by the rheology property of granite and diorite (Brace,Kohlstedt, 1980; Sibson, 1982; Meissner, Strehlau, 1982; CHEN, Molnap 1983). This is obviouslyrelated to the change of rock deformation characteriStics at different crust depths. So, for earthquake stUdy, both of the rock failure types (fractUre or rock flow) and its mechanical instabilityforms (…     

Displacement capacity of masonry piers: parametric numerical analyses versus international building codes

The nonlinear behaviour of masonry piers loaded in their plane is investigated by parametric numerical simulations. Each pier has a cantilever scheme, is loaded by a constant axial load and is subjected to an increasing horizontal displacement at the top. The macro-modelling approach is used to perform numerical analyses, adopting two different constitutive laws: a total strain crack model and a plastic model. The numerical model is calibrated on a block-masonry type for which experimental tests are available in literature. Parametric numerical simulations are performed by varying the aspect-ratio and the compression level, in order to assess the influence of such parameters on both shear strength and displacement capacity. By comparing numerical results with formulas of international codes, a good agreement for the shear strength is obtained, while significant differences are observed for the displacement capacity, which is influenced by both parameters. The authors propose a simple empirical formula for the displacement capacity, obtained by fitting the numerical results. The expression can be useful in the practical design for considering the influence of aspect-ratio and compression level, currently neglected by building codes.     

岩石声发射的损伤模式及其在地震研究中的初步应用

本文运用损伤理论和微元强度统计理论导出了声发射参数与损伤参量及应变的关系模型。依据这一模型,分析了单轴压缩时岩石破坏阶段声发射特性与压机刚度的关系。在此基础上,探讨了声发射活动规律与孕震过程及地震活动的内在本质联系。     

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.     

圆钢管约束高强灌浆料局部承压性能研究

以圆钢管约束高强灌浆料局部承压性能试验为基础,基于ABAQUS有限元计算平台,通过建立分段模型(Multi-section Model简称MSM模型)进行数值模拟,探究圆钢管约束高强灌浆料在局部承压作用下的力学性能,并将软件计算得到的荷载-位移曲线与试验结果进行对比分析。结果表明,分段模型适用于圆钢管约束高强灌浆料局部承压作用下的数值模拟,且计算结果与试验结果吻合较好;增大试件高径比、钢管壁厚和高强灌浆料强度均可提高其极限承载力。     

岩石损伤软化的修正本构模型

利用服从Weibull分布的岩石微元强度表示方法,基于Drucker-Prager破坏准则,通过引入损伤变量的修正系数,利用应力-应变关系曲线峰值点处的应力、应变确定Weibull分布参数的关系式,建立不同围压下的损伤软化本构模型。该模型参数较少且易于确定,其参数的确定方法揭示了模型参数的物理意义。与前人的研究成果进行对比分析,结果表明该模型比未修正前有着明显的优越性,且与实测结果吻合较好,分析结果显示出该模型的合理性。     

Identification and mapping of riverbed sediment facies in the Columbia River through integration of field observations and numerical simulations

In the Hanford Reach of the Columbia River, a thin layer of recent alluvium overlies the sedimentary formations that comprise the unconfined groundwater aquifer. Experimental and modelling studies have demonstrated that this alluvial layer exerts significant control on the exchange of groundwater and surface water (hydrologic exchange flux), and is associated with elevated levels of biogeochemical activity. This layer is also observed to be strongly heterogeneous, and quantifying the spatial distribution of properties over the range of scales of interest is challenging. Facies are elements of a sediment classification scheme that groups complex geologic materials into a set of discrete classes according to distinguishing features. Facies classifications have been used as a framework for assigning heterogeneous material properties to grid cells of numerical models of flow and reactive transport in subsurface media. The usefulness of such an approach hinges on being able to relate facies to quantitative properties needed for flow and reactive transport modelling, and on being able to map facies over the domain of interest using readily available information. Although aquifer facies have been used in various modelling contexts, application of this concept to riverbed sediments is relatively new. Here, we describe an approach for categorizing and mapping recent alluvial (riverbed) sediments based on the integration of diverse observations with numerical simulations of river hydrodynamics. The facies have distinct distributions of sediment texture that correspond to variations in hydraulic properties, and therefore provide a useful framework for assigning heterogeneous properties in numerical simulations of hydrologic exchange flows and biogeochemical processes.     

A Damage Model of Acoustic Emission in Rock and Its Initial Application in Seismological Research

The model of the relationship of AE parameter,damage variable and strain is derived by applying the damage theory and micro-element statistical strength theory.The relation between AE characteristics during rock failure and machine stiffness is analyzed under uniaxial compression with the above model.Based on the above analysis,the internal connection among AE activity law and seisraogenic process and earthquake activity is discussed.     

Spatial and temporal characteristics of acoustic emission activity during deformation of rock samples with inhomogeneous fault under biaxial compression

The spatial and temporal distributions of acoustic emission （AE） during the deformation of samples containing an inhomogeneous fault have been studied under biaxial compression. The results show that the fault strength and the duration from loading to failure increase and the failure mode changes from abrupt instability to gradual failure with increase of lateral stress σ2. The pre-setting fault and its heterogeneity play an important role in controlling AE spatial distribution during the deformation. The basic pattern of AE spatial distribution is controlled by the pre-setting fault, especially the parts with inhomogeneous strength and the strong segment, and the localization of fracturing starts from the positions with inhomogeneous strength, With increase of σ5, the dense AE distribution area spreads from the positions with inhomogeneous strength to the whole strong segment gradually. AE temporal sequence is significantly affected by σ2. The fault shows abrupt instability in final failure, which occurs on a background of ＂enhancement-quiescence＂ in microfracturing activity at lower σ2. At higher σ5, the fault shows gradual failure, and AE occurrence rate increases continuously and AER increase exponentially before and after the failure, The effect of σ2 on b-value is also remarkable, b-value shows precursory decrease in the weakening stage when fault behavior is abrupt instability, but it shows balanced change when fault behavior is gradual failure,     

SPATIAL DISTRIBUTION AND FRACTAL STRUCTURE OF AE FOCUSES OF INADA GRANITE UNDER TRIAXIAL COMPRESSION

An experimental study on the spatial distribution of the focuses of acoustic emissions of Inada granite (coarse grained granodiorite) under triaxial compression is reported in this paper. It is pointed out that AE locations is fractal with a dimension of 2.2 in average, lower than the result of Oshima granite (fine grained granodiorite), which is 2.7, under similar deformation condition. For Inada granite a band-limited fractal had been observed. Fractal dimension changed around grai     

Simulation of the Micro-physics of Rocks Using LSMearth

-- The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (M0ra and Place, 1994; Place and Mora, 1999). A new modular and flexible LSM approach has been developed that allows different micro-physics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.     

论地震序列的一种可能的重要成因

本文在Asperity模式和Barrier模式的基础上,提出了另外一种不均匀体模式——障碍体模式,并对该模式进行了数值模拟的研究。结果除证实了我们以前做过的单轴压力下模拟障碍体的实验结果以外,还表明在含有障碍体的断层面上,破裂的扩展与破裂的传播速度、始破裂区域的大小、裂纹的类型以及障碍体所在断层的介质有关。在有围压的情况下,破裂可以先越过或绕过强度较高或尺度较大的障碍体,经一段时间的应力积累过程后,再使障碍体发生破坏。本文通过对震源区中障碍体的性质和影响的研究,解释了海城和邢台两个地震序列形成和发展的可能原因,得到一些较为合理的结果。     

油气田典型岩石三轴压缩变形破坏与声发射活动特征——四川盆地震旦系白云岩及页岩的破坏过程

随着二氧化碳地质封存、深部地热开采、地下储气库建设、页岩气开发、二次驱油/驱气等工业应用的快速发展,与地下流体注入有关的诱发地震活动呈现一定的增加趋势.利用声发射实验观测油气田典型岩石在三轴压缩条件下变形破坏过程与声发射活动特征,对研究注水诱发地震过程有着重要意义.本文利用四川盆地现场采集的震旦系白云岩及页岩,采用实验室声发射技术观测研究岩石三轴压缩变形破坏过程中地震波速度等物性参数及声发射事件时空分布特征.实验结果表明:震旦系白云岩及页岩在变形破坏过程中均有一定的声发射活动.根据声发射定位结果,声发射主要集中在破坏前后的较短时间内,页岩的层理面为结构弱面,控制最终破坏面的形态及声发射特征.根据应力-应变结果,白云岩在压缩的后期阶段有一定的扩容现象,但页岩在整个压缩阶段均没有明显扩容现象.研究结果表明四川盆地较古老的白云岩及页岩具有脆性破坏特征,地下流体注入容易诱发微震活动,形成裂缝,有利于页岩气压裂开采.微震活动有利于监测裂缝的发生发展,但同时在页岩气开发及二氧化碳地质封存时应采取相应的预防控制措施进行安全合理的储盖层管理,避免灾害性诱发地震的发生.     

高孔隙岩石变形的离散单元模型

实验和野外观测表明颗粒破碎显著影响高孔隙岩石的变形特征.为建立高孔隙岩石变形的数值模型,我们以弹性理论为基础并根据问题的特征进行合理简化,给出了一种颗粒破碎机制,并根据此颗粒破碎条件改进传统的离散单元模型.利用改进的离散单元模型,研究了不同压力条件下高孔隙岩石变形的宏观和微观特征.数值模拟结果表明考虑颗粒破碎的模型能重现高孔隙岩石在不同围压下变形的应力-应变关系、声发射特征、应变局部化和剪切增强压缩等特征.我们还发现,与野外观测结果一致,在高孔隙岩石脆性阶段剪切带内部及附近伴有显著的颗粒破碎.     

Shear Band Formation in Numerical Simulations Applying a Continuum Damage Rheology Model

In seismically active regions, faults nucleate, propagate, and form networks that evolve over time. To simulate crustal faulting processes, including the evolution of fault-zone properties, a rheological model must incorporate concepts such as damage rheology that describe the various stages of the seismic cycle (strain localization, subcritical crack growth and macroscopic failure) while accounting for material degradation and healing and off-fault deformation. Here we study the fundamental patterns of strain-localisation within the framework of a continuum damage rheology by performing a shear band analysis (linear instability analysis) and comparing predictions of shear band orientations with numerical results of the nonlinear problem. We find (analytically and numerically) that the angle between the shear band and the less compressive (transverse) direction is between 47° in compression tests with a strain ratio of 0.25 (highly confined compression test), and 60° for a strain ratio higher than 1.4 (axial compression and transverse extension). In addition we find that shear bands exhibit local dilation (I ₁ > 0) in a wide range of strain ratios excluding only simulations with highly confined compression (which yield compacting shear bands or non-localized deformation). Finally, we discuss the applicability of the damage model for simulating deformation in the seismogenic, brittle crust.     

有限元模拟弹性位错的等效体力方法

计算地震位错造成的位移场和应力场,对于估计大地震引起的后续地震活动发展趋势十分重要.本文提出在弹性位错问题的有限元模拟中,用等效体力代替位错源,从而在构建几何模型时不用包含断层,却可以处理包含任意复杂断层的问题,极大降低建模的难度.使用此方法,本文计算并讨论了在球形地球模型下2011年日本Tohoku-Oki特大地震对华北地区断层的影响,结果表明此次地震使华北主要断层趋于稳定.     

Permeability Evolution and Rock Brittle Failure

This paper reports an experimental study of the evolution of permeability during rock brittle failure and a theoretical analysis of rock critical stress level. It is assumed that the rock is a strain-softening medium whose strength can be described by Weibull’s distribution. Based on the two-dimensional renormalization group theory, it is found that the stress level λ_c (the ratio of the stress at the critical point to the peak stress) depends mainly on the homogeneity index or shape parameter m in the Weibull’s distribution for the rock. Experimental results show that the evolution of permeability is closely related to rock deformation stages: the permeability has a rapid increase with the growth of cracks and their surface areas (i.e., onset of fracture coalescence point), and reaches the maximum at rock failure. Both the experimental and analytical results show that this point of rapid increase in permeability on the permeability-pressure curve corresponds to the critical point on the stress-strain curve; for rock compression, the stress at this point is approximately 80% of the peak strength. Thus, monitoring the evolution of permeability may provide a new means of identifying the critical point of rock brittle fracture.     

从断裂力学观点探讨 b 值的物理实质

岩石试验声发射及天然地震的 b 值到底反映了什么样的物理实质,是一个长期以来引起争论的问题.本文运用断裂力学的观点和方法,研究岩石试验试件内微裂纹系中各裂纹的扩展顺序及声发射能量,由此确定整个声发射序列,并进而确定 b 值.如果假设裂纹长度的分布密度函数为 p()=B-.经过力学分析和数学推导,可得到关系式 b=3/2.由此看出:b 值的物理实质是,它反映了介质的断裂构造状态————包括介质中裂纹系的空间分布及其它影响材料中裂纹扩展的物理量(如断裂韧度、摩擦系数等)的空间分布.这一结论与大多数实验及观测结果相符.我们的研究结果与茂木清夫观点有共同之处.茂木清夫认为:介质的不均匀性是决定 b 值的重要因素.介质的不均匀性一词固然也包括了介质的断裂构造状态在内,但是本文之观点比茂木更进一步说明了问题的实质.      

Seismic velocities in fractured rocks: an experimental verification of Hudson's theory1

Flow of fluids in many hydrocarbon reservoirs and aquifers is enhanced by the presence of cracks and fractures. These cracks could be detected by their effects on propagation of compressional and shear waves through the reservoir: several theories, including Hudson's, claim to predict the seismic effects of cracks. Although Hudson's theory has already been used to calculate crack densities from seismic surveys, the predictions of the theory have not yet been tested experimentally on rocks containing a known crack distribution. This paper describes an experimental verification of the theory. The rock used, Carrara marble, was chosen for its uniformity and low porosity, so that the effect of cracks would not be obscured by other influences. Cracks were induced by loading of laboratory specimens. Velocities of compressional and shear waves were measured by ultrasound at 0.85 MHz in dry and water-saturated specimens at high and low effective pressures. The cracks were then counted in polished sections of the specimens. In ‘dry’ specimens with both dry and saturated cracks, Hudson's theory overpredicted observed crack densities by a constant amount that is attributed to the observed value being systematically underestimated. The theory made poor predictions for fully saturated specimens. Shear-wave splitting, caused by anisotropy due to both crystal and crack alignment, was observed. Cracks were seen to follow grain boundaries rather than the direction of maximum compression due to loading. The results demonstrate that Hudson's theory may be used in some cases to determine crack and fracture densities from compressional- and shear-wave velocity data.