损伤和剪胀效应对裂隙煤体渗透率演化规律的影响研究

煤层瓦斯渗透率是影响瓦斯抽采和动力灾害防治的重要参数。为了研究煤体损伤和剪胀变形对渗透率的影响,首先引入损伤变量反映煤体损伤破坏状态,建立了基于体应变增量的煤体损伤本构模型。并采用Hurst指数表征裂隙表面粗糙度,基于裂隙表面的分形特征,建立了裂隙渗透率在压缩和剪切作用下的演化模型。通过对TOUGH2和FLAC3D软件进行二次开发,建立了基于双重孔隙模型的TOUGH2(CH4)-FLAC气-固耦合数值分析工具。采用本软件对煤样单轴压缩过程进行模拟分析,结果表明：煤体的破坏是损伤单元累积和贯通的结果,最终形成贯通煤体的损伤带是造成煤体失稳破坏的主因;围岩内的渗透率增加区域与损伤区位置基本一致,其中裂隙系统的渗透率增加幅度最大可达2个数量级;剪切破坏区的裂隙发生剪胀变形,引起裂隙渗透率大幅增加。建立的理论模型与数值计算工具为制定瓦斯治理措施提供了理论指导。     

Gas permeability evolution of cataclasite and fault gouge in triaxial compression and implications for changes in fault-zone permeability structure through the earthquake cycle

We report the results of permeability measurements of fault gouge and tonalitic cataclasite from the fault zone of the Median Tectonic Line, Ohshika, central Japan, carried out during triaxial compression tests. The experiments revealed marked effects of deformation on the permeability of the specimens. Permeability of fault gouge decreases rapidly by about two orders of magnitude during initial loading and continues to decrease slowly during further inelastic deformation. The drop in permeability during initial loading is much smaller for cataclasite than for gouge, followed by abrupt increase upon failure, and the overall change in permeability correlates well with change in volumetric strain, i.e., initial, nearly elastic contraction followed by dilatancy upon the initiation of inelastic deformation towards specimen failure. If cemented cataclasite suffers deformation prior to or during an earthquake, a cataclasite zone may change into a conduit for fluid flow. Fault gouge zones, however, are unlikely to switch to very permeable zones upon the initiation of fault slip. Thus, overall permeability structure of a fault may change abruptly prior to or during earthquakes and during the interseismic period. Fault gouge and cataclasite have internal angles of friction of about 36° and 45°, respectively, as is typical for brittle rocks.     

A microstructural study of fault rocks from the SAFOD: Implications for the deformation mechanisms and strength of the creeping segment of the San Andreas Fault

The San Andreas Fault zone in central California accommodates tectonic strain by stable slip and microseismic activity. We study microstructural controls of strength and deformation in the fault using core samples provided by the San Andreas Fault Observatory at Depth (SAFOD) including gouge corresponding to presently active shearing intervals in the main borehole. The methods of study include high-resolution optical and electron microscopy, X-ray fluorescence mapping, X-ray powder diffraction, energy dispersive X-ray spectroscopy, white light interferometry, and image processing.The fault zone at the SAFOD site consists of a strongly deformed and foliated core zone that includes 2–3 m thick active shear zones, surrounded by less deformed rocks. Results suggest deformation and foliation of the core zone outside the active shear zones by alternating cataclasis and pressure solution mechanisms. The active shear zones, considered zones of large-scale shear localization, appear to be associated with an abundance of weak phases including smectite clays, serpentinite alteration products, and amorphous material. We suggest that deformation along the active shear zones is by a granular-type flow mechanism that involves frictional sliding of microlithons along phyllosilicate-rich Riedel shear surfaces as well as stress-driven diffusive mass transfer. The microstructural data may be interpreted to suggest that deformation in the active shear zones is strongly displacement-weakening. The fault creeps because the velocity strengthening weak gouge in the active shear zones is being sheared without strong restrengthening mechanisms such as cementation or fracture sealing. Possible mechanisms for the observed microseismicity in the creeping segment of the SAF include local high fluid pressure build-ups, hard asperity development by fracture-and-seal cycles, and stress build-up due to slip zone undulations.     

On the role of particle breakage in the shear failure behavior of granular soils by DEM

This article presents a fundamental study on the role of particle breakage on the shear behavior of granular soils using the three‐dimensional (3‐D) discrete element method. The effects of particle breakage on the stress ratio, volumetric strain, plastic deformation, and shear failure behavior of dense crushable specimens undergoing plane strain shearing conditions are thoroughly investigated through a variety of micromechanical analyses and mechanism demonstrations. The simulation of a granular specimen is based on the effective modeling of realistic fracture behavior of single soil particles, which is demonstrated by the qualitative agreement between the results from platen compression simulations and those from physical laboratory tests. The simulation results show that the major effects of particle breakage include the reduction of volumetric dilation and peak stress ratio and more importantly the plastic deformation mechanisms and the shear failure modes vary as a function of soil crushability. Consistent macro‐ and micromechanical evidence demonstrates that shear banding and massive volumetric contraction depict the two end failure modes of a dense specimen, which is dominated by particle rearrangement–induced dilation and particle crushing–induced compression, respectively, with a more general case being the combination and competition of the two failure modes in the medium range of soil crushability and confining stress. However, it is further shown that a highly crushable specimen will eventually develop a shear band at a large strain because of the continuous decay of particle breakage. Copyright © 2011 John Wiley & Sons, Ltd.     

断裂在纯净砂岩中的变形机制及断裂带内部结构

以纯净砂岩为对象,考虑了影响断裂变形的主要因素,包括成岩阶段、孔隙度、温度和围压,系统剖析了不同性质砂岩的断裂变形机制、微构造特征及形成断裂带的内部结构,对研究砂岩内断裂封闭性具有重要的指导意义。研究表明：纯净的砂岩在未固结-半固结成岩阶段发生断裂,变形机制为颗粒边界摩擦滑动,导致颗粒旋转和滚动,即为颗粒流,形成的微构造为解聚带,孔渗性同母岩比没有明显降低,断裂带尽管具有断层核和破碎带二分结构,但渗透率比母岩高,为流体垂向运移的通道;在固结成岩阶段（孔隙度大于15%）发生断裂,变形机制为碎裂作用,颗粒边界摩擦滑动导致颗粒旋转,即为碎裂流,形成的微构造为碎裂带,渗透率同母岩比一般降低1～3个数量级,形成侧向有一定封闭能力、垂向渗透的断裂带;固结阶段（孔隙度小于15%）发生断裂,开始由于破裂作用,形成断层角砾岩,伴随着碎裂流发生,形成碎裂岩,因此早期形成高渗透断裂带,后期断层逐渐封闭。固结成岩的砂岩在抬升过程发生断裂,变形机制为破裂作用,形成无内聚力的角砾岩,为高渗透断裂带。在不同成岩阶段发生变形,形成多类型变形构造的叠加,对于一条晚期形成的断层而言,由于不同深度变形机制及微构造类型不同,导致油气选择性充注,碎裂带和压溶胶结碎裂带阻止油气向高孔隙度砂岩中充注,解聚带会成为油气运移的通道,裂缝有利于油气优先充注。因此,高孔隙性砂岩中孔隙度较低的储集层由于碎裂带不发育常常含油气性最好,而低孔隙性砂岩由于裂缝产生含油气性较好。     

Mechanisms of melt extraction during lower crustal partial melting

Progressive vapour‐absent partial melting of a closed rock system increases melt pressure due to an expansion in the volume of the mineral plus melt assemblage. For a locally closed system, we quantify the melt pressure increase per increment of partial melting of a metapelite using phase equilibria modelling and combine it with Mohr–Coulomb theory to examine the interplay between melt pressure and fracture behaviour. It is shown that very small increments of vapour‐absent partial melting (<1%) increase melt pore pressure by tens of MPa leading to inevitable brittle failure of locally closed systems. Fracturing will affect these systems, even if initially limited to the scale of a few grains, and a connected microfracture network will enhance permeability as partial melting progresses. This will lead to a conditionally open system, potentially limiting accumulation of melt in the source. Repeated and cyclic fracture as temperature progressively increases will drive migration of the melt into sites of low fluid pressure at all scales. Crystal‐plastic creep processes create deformation‐induced dilatancy gradients that dominate over buoyancy forces at all scales in the melt source. Brittle and ductile deformation therefore cooperate in the extraction of melt. Enhanced porosity and permeability in ductile shear zones result in lower fluid pressure, providing a potentially important driving force for melt migration and drainage ‘up’ shear zones and along larger scale fluid pressure gradients in the crust.     

Groundwater targeting in a hard-rock terrain using fracture-pattern modeling, Niva River basin, Andhra Pradesh, India

In hard-rock terrain, due to the lack of primary porosity in the bedrock, joints, fault zones, and weathered zones are the sources for groundwater occurrence and movement. To study the groundwater potential in the hard-rock terrain and drought-prone area in the Niva River basin, southern Andhra Pradesh state, India, Landsat 5 photographic data were used to prepare an integrated hydrogeomorphology map. Larsson's integrated deformation model was applied to identify the various fracture systems, to pinpoint those younger tensile fracture sets that are the main groundwater reservoirs, and to understand the importance of fracture density in groundwater prospecting. N35°–55°E fractures were identified as tensile and N35°–55°W fractures as both tensile and shear in the study area. Apparently, these fractures are the youngest open fractures. Wherever N35°–55°E and N35°–55°W fracture densities are high, weathered-zone thickness is greater, water-table fluctuations are small, and well yields are high. Groundwater-potential zones were delineated and classified as very good, good to very good, moderate to good, and poor. Electronic Publication     

郯庐断裂带南段桴槎山韧性剪切带糜棱岩的变形条件和组分迁移系

本文对郯庐断裂带南段主干断裂典型的韧性剪切带进行了系统的剖析。从糜棱岩塑性变形的亚颗粒化、动态重结晶和矿物成分特征及岩石组分迁移变化等入手进行系统的研究工作，计算了岩石形成的温度、压力和流动应力和流变速率参数；模拟计算了岩石在剪切变形作用下的体积亏损及组分迁移的量值，探讨了变形－变质及流体的相关关系     

Cataclastic deformation of quartzite in the moine thrust zone

Deformed Cambrian quartzites from the Moine thrust zone on Skye have been examined under cathodoluminescence in the Scanning Electron Microscope (SEM). A sequence of progressive deformation from intact quartzite through protobreccia to breccia and ultrabreccia can be established for samples from a large fold with a wavelength of 2 km, and a comparable cataclasite sequence exists for samples approaching the Ord thrust plane. The microstructural evolution involves the development of extension microcracks by impingement at grain contacts, followed by the formation of small shear faults by linking of the extension microcracks. Further strain is localized on large breccia zones; the proportion of fine-grained matrix increases and fluid flow is concentrated into the dilatant zones, depositing cement and iron oxides. The dominant deformation mechanism for both folding and faulting in this part of the Moine thrust zone is cataclasis. Extension microcracks and shear faults have close analogues in experimental deformation, but the formation of the protobreccia in the hinge of the fold is an example of cataclastic flow on a much larger scale than a laboratory specimen. The conditions for the initiation of the shear faults and breccia zones are suggested by the Rudnicki and Rice model of a dilatant material. Work-hardening must occur in both cases; possible mechanisms include changes in the dilatancy factor and bulk modulus, and syntectonic cementation.     

甘孜嘎拉脆—韧性剪切带的构造变形特征及对金矿的控制

本文阐述了甘孜嘎拉脆—韧性剪切带的构造变形特征、运动学标志及对金矿的控制,提出甘孜嘎拉金矿床主要的控矿构造为与甘孜—理塘断裂带平行展布的、具多期活动复合叠加的次级脆—韧性剪切变质带。     

Fractures,hydrothermal alterations and permeability in the Soultz Enhanced Geothermal System

Borehole studies of the Soultz-sous-Forêts granite are dedicated to deep geothermics. The hydraulic properties of the reservoir are mainly controlled by the occurrence of some altered cataclastic shear zones showing a low natural permeability characterized by the occurrence of brines. Those zones show a fracture cluster organisation with sealed fractures of various types (post-filled joints, sheared fractures, veins). The main hydrothermal deposits observed within the permeable zones are geodic quartz, carbonates, illite and more locally sulphides. The fracture wall–rocks are intensely transformed: dissolution of igneous minerals, crystallization of new minerals, porosity and permeability increase. It is important to characterize the newly-formed minerals in order to choose the reagents used to improve the permeability of the exchanger by chemical stimulations. This article represents a synthesis of the studies completed by the authors between 1990 and 2008 on the fracture networks, hydrothermal alterations and mineral crystallizations they induced and data about the flow pathways in the exchanger.     

Strain and Crystallographic Fabric in Mesoscopic Ductile Shear Zones of Garhwal Himalaya

The MCT Zone of Bhagirathi valley of Garhwal Himalaya is characterized by numerous mesoscopic ductile shear zones. These shear zones are developed in response to nearly NNE-SSW maximum horizontal compression and provide an opportunity to study the variation in strain and crystallographic fabrics within the ductile shear zones.The grain shape and orientation of quartz under microscope reflect that strain is higher in the center and it progressively decreases towards the shear zone boundary. The preferred orientation of quartz c-axes across the shear zone suggests that the single girdle of the quartz c-axes are probably first developed at the shear zone boundary and become prominent in the center of shear zone with increase in the intensity of deformation. The strong crystallographic preferred orientation normal to foliation suggests that the internal deformation of the quartz might have taken place by dislocation creep mechanism exhibiting a non-coaxial deformation history.     

Deformation path in high-strain zones,with reference to slip partitioning in transpressional plate-boundary regions

The current status of the kinematics and strain geometry of high-strain zone studies is briefly summarized. A general high-strain zone has a triclinic deformation path, and monoclinic shear zones are special end member cases. Fabrics observed in natural shear zones and theoretical considerations based on continuum mechanics are compatible with this conclusion. Non-steady deformation paths remain difficult to deal with, and may ultimately rely on a realistic mechanical treatment of high-strain zones which may be possible when our knowledge of the mechanical behavior of rocks under natural deformation conditions is improved. An examination of the phenomenon of slip partitioning in transpressional plate boundary regions shows that the bulk deformation path in the forearc area (trench-parallel high-strain zone) is generally triclinic. The Alpine Fault in the South Island of New Zealand provides an example of a currently active triclinic shear zone. The Southern Knee Lake shear zone of Manitoba, Canada, provides an Archean example of a triclinic shear zone.     

Alteration-weakening leading to localized deformation in a damage aureole adjacent to a dormant shear zone

Deformation adjacent to faults and shear zones is traditionally thought to correlate with slip. Inherited structures may control damage geometry, localizing fluid flow and deformation in a damage aureole around structures, even after displacement has ceased. In this paper we document a post-shearing anastomosing foliation and fracture network that developed to one side of the Mesoarchean Marmion Shear Zone. This fracture network hosts the low-grade, disseminated Hammond Reef gold deposit. The shear zone juxtaposed a greenstone belt against tonalite gneiss and was locked by an intrusion that was emplaced during the final stages of suturing. After cessation of activity, fluids channeled along fault- and intrusion-related fractures led to the pervasive sericitization of feldspars. Foliated zones resulted from flattening in the weaker sericite-rich tonalite during progressive alteration without any change in the regional NW-SE shortening direction. The anastomosing pattern may have been inherited from an earlier ductile fabric, but sericite alteration and flattening fabrics all formed post-shearing. Thus, the apparent foliated fracture network adjacent to the Marmion Shear Zone is a second-order effect of shear-related damage, distinct in time from shear activity, adjacent to an effectively dormant shear zone. This phenomenon has implications for understanding the relative timing of fault zone activity, alteration and (in this case) gold mineralization related to long-term fault zone permeability.     

Geotechnical investigations at the Kalyani dam site, Chitoor District, Andhra Pradesh, India

Details are presented of site investigations carried out for the Kalyani Dam, Chittoor District, Andhra Pradesh, India. Geologically, the dam-site region forms a triple junction of the Dharwar, Easternghat and Cuddapah orogenies and it is tectonically disturbed; the course of the Kalyani river is guided by a major strike-slip fault. The geophysical studies reveal the presence of several fault and shear zones in the area. Laboratory studies on the physical and engineering properties of the rocks disclose distinct differences in the sheared (Zone I) rocks as against those that are relatively undisturbed (Zone III). Zone I samples, in general, show high porosity and elastic an isotropy, low values of density, elastic wave velocity, elastic moduli, fracture strength and high absorption, while the reverse is noticed in the undisturbed samples of Zone III. Samples from the shear zone further show large variations in compressional velocity and amplitude, as a function of compressive stress. The laboratory velocities are in good agreement with the field seismic refraction results. These studies helped in estimating the soundness of the rock formations, delineating the structurally weak zones and suggesting remedial measures.     

Small-scale ductile shear zones as transposed rheologically weak domains: A numerical modeling investigation and practical application

We test the hypothesis that small ductile shear zones are developed from initial rheologically weak domains. We regard weak domains as ellipsoidal inhomogeneities and apply Eshelby's formalism extended for power-law viscous materials to investigate the kinematics and finite strain evolution of the partitioned flow fields in weak domains. We show that, under an imposed bulk flow field, weak domains, regardless of their initial shapes and orientations, deform into zone-like features at relatively small bulk strains and the deformation paths inside weak domains have all characteristics expected in ductile shear zones. We apply our model to the Cap de Creus area Spain, where abundant small ductile shear zones exist. To relate the fabrics inside these shear zones with the regional deformation, we take a multi-scale approach. We assume that the area contained many weak domains which were randomly shaped and oriented initially and were deformed into shear zones eventually in response to the regional flow. We constrain the regional flow field by the fabric patterns across the area and compute numerically the partitioned flow fields in individual weak domains. The latter are related to fabrics inside shear zones. Our model reproduces first-order features of field-observed fabrics. Although the deformation path of each shear zone was close to simple shearing, the deformation of the whole belt was not. Our approach also resolves the strain compatibility problem for a finite-sized shear zone embedded in a far less deformed country rock.     

Structural controls of lode-gold mineralization by mafic dykes in late-Paleozoic granitoids of the Kochkar district, southern Urals, Russia

The Kochkar gold district in the East Uralian Zone of the southern Urals is located in late-Paleozoic granite gneisses of the Plast massif. Gold mineralization is associated with tabular quartz lodes that are preferentially developed along the margins of easterly trending mafic dykes. Fabric development indicates that dykes had a profound influence on the development of shear zones in granitoids. ENE- and SE-trending dykes have been reactivated as dextral and sinistral oblique strike-slip shear zones, respectively, forming a set of approximately conjugate shear zones related to the Permian, regional-scale E-W directed shortening. Dyke-shear zone relationships in the Plast massif are the result of strain refraction due to the presence of biotite-rich, incompetent dykes in more competent granite-gneisses. Deformation and the formation of associated gold-quartz lodes occurred close to peak-metamorphic, upper-greenschist to lower-amphibolite facies conditions. Strain refraction has resulted in partitioning of the bulk strain into a component of non-coaxial mainly ductile shear in mafic dykes, and a component of layer-normal pure shear in surrounding granitoids where deformation was brittle-ductile. Brittle fracturing in granitoids has resulted in the formation of fracture permeabilities adjacent to sheared dykes, that together with the layer-normal dilational component, promoted the access of mineralizing fluids. Both ore-controlling dykes and gold-quartz lodes were subsequently overprinted by lower greenschist-facies, mainly brittle fault zones and associated hydrothermal alteration that post-date gold mineralization. Received: 15 October 1998 / Accepted: 18 August 1999     

Microstructural study of a natural slip zone: quantification and deformation history

Microstructure of a natural slip zone was comprehensively examined using a combination of images captured systematically by optical microscopy (OPM) and backscattered electron microscopy (BEM) techniques. Microstructural features identified on these images were processed and evaluated using an advanced image analysis system, which proved that quantitative analyses could considerably enhance the understanding of shear behavior of slip zones. It was found that variations of porosity, abundance of platy clay particles and alignments of particles are significant indicators revealing nature of deformation processes. These indicators show that global mechanical behavior of the investigated slip zone can be conceptualized as that of normally consolidated clayey soils under drained conditions.

The geometric patterns of the microstructure of the slip zone are similar to the S–C fabrics seen in tectonic shear zones. It is suggested that combined progressive bulk simple shear and pure shear modes enable to realistically reconstruct the kinematic history of the slip zone, through which particle movements and microstructural evolution were accomplished via various types of particulate flows. The results of this study show that clay mineralogy plays a more important role in the development of the slip zone than abundance of clay-size particles, while both clay mineralogy and relative proportions of each particle size fraction control the response of particles to shear deformation. Among the fractions present in the slip zone, fine silts are the strongest indicator of global shear stress characterized by their highest degree of alignment, whereas clay particles are the weakest. Highest degree of shape preferred orientation is also found within fine silt domains.     

八宝山断裂带浅部岩体水力耦合模型

八宝山断裂带是北京西部的一条活动断裂,观测表明其形变、地下水位和降水之间存在强相关性。从大灰厂观测站的断裂带条件和观测系统布置来判断,这种强相关性是浅部破碎岩体渗流-应力耦合机制的结果。断裂带浅部岩体可以分为无压区和承压区。降水能够直接入渗补给无压区,也能够渗入到承压区上部的上盘岩体风化带,从而同时影响断裂带和风化带的地下水位。一个7参数的线性水力耦合模型给出了与实测动态基本一致的垂向形变、地下水位模拟结果。模型分析表明观测到的垂向变形主要由无压区地下水位升降引起的膨胀和收缩变形构成,而上盘湿润或疏干引起的重量变化对垂向变形贡献很小。由于断裂带的地下水排泄基准面和水准基台值可能存在某种外在因素引起的变化,模型参数可能不够稳定,需要进一步调查研究。     

Effects of non‐uniform traction and specimen height in the direct shear test on stress and deformation in a rock fracture

In the direct shear test (DST), an internal moment is distributed within the rock specimen by non‐coaxial shear loads applied to the specimen, which cause non‐uniform distributions of both the traction on the loading planes and the stress and deformation in the specimen. To examine the validity of the DST for a rock fracture and to clarify the effect of specimen height, both the stress and deformation in a fracture in the DST were analyzed for specimens with three different heights using a three‐dimensional finite element method with quadratic joint elements for a fracture model. The constitutive law of the fracture considers the dependence of the non‐linear behavior of closure on shear displacement and that of shear stiffness on normal stress and was implemented in simulation code to give a conceptional fracture with uniform mechanical properties to extract only the effect of non‐uniform traction on the stress and deformation in the fracture. The results showed that both normal and shear stresses are concentrated near the end edges of the fracture, and these stress concentrations decrease with a decrease in the specimen height according to the magnitude of the moment produced by the non‐coaxial shear loads. Furthermore, although closure is greater near the end edges of the fracture, where normal stress is concentrated, this concentration of closure is not so significant within the range of this study because of the non‐linear behavior of closure, that is, closure does not significantly increase with an increase in normal stress at large normal stresses. Copyright © 2012 John Wiley & Sons, Ltd.