围压作用下岩石样品中微裂纹的闭合

岩石中普遍存在的微裂纹对其物理性质有十分显著的影响.为了得到不同围压条件下岩样内部微裂纹形态的变化,选取几种具有明显层理的岩石圈岩石样品,在不同围压下沿相互正交的三个方向对纵横波速度和线应变等进行了测量,根据O’Connell模型通过计算得出有效泊松比、裂纹密度、裂纹孔隙度和裂纹的c/a比值（a为硬币形裂纹半长轴的长度,c为其半宽度）.结果表明,随着围压的增加,在垂直于层理面的方向上裂纹闭合最快但c/a比值不减反增.结合显微组构分析可知,岩石中的裂纹表面实际上是凸凹不平的,在围压增加引起裂纹闭合过程中一条裂纹可以变成几个,长度成倍减小但宽度变化不大,所以在垂直于层理面的方向上c/a比值随围压增加而增大.但另两个方向上的变化明显不同,应该与模型本身的适用范围有关.通过同时测量纵、横波速度,计算c/a比值,可以对样品中微裂纹的几何形态及空间分布方式有定性或半定量的了解.     

Elastic wave velocities in a granitic geothermal reservoir

Geothermal resources have potential for providing cost-effective and sustainable energy. Monitoring of production-induced changes in geothermal reservoirs using seismic waves requires understanding of the elastic properties of the rock and how they change due to injection of fluids and opening and closing of natural and hydraulic fractures. P- and S-wave velocities measured in a granitic geothermal reservoir using sonic logging are systematically lower than those predicted using the composition of the rock. Cracks may occur in granitic rocks from tectonic stresses and from the thermal expansion mismatch between differently oriented anisotropic crystals. An isotropic orientation distribution of microcracks causes a significant reduction in both the P- and S-velocities, consistent with the observed sonic P- and S-velocities. Vertical fractures cause a difference in the velocity of vertically propagating shear waves polarized parallel and perpendicular to the fractures. An assumption that the lower measured velocities are caused by the presence of vertical fractures is inconsistent with the sonic data. This is because vertical fractures cause a decrease in slow S-wave velocity that greatly exceeds the decrease in P-wave velocity, in contrast to the observed data. The growth of vertical fractures in the geothermal reservoir may be monitored using the difference in velocity of the fast and slow shear waves, while the change in P-velocity in a crossplot of measured P- and slow S-velocities is useful for estimating the ratio of the normal-to-shear compliance of the fractures.     

两类三轴实验花岗岩的微破裂对比研究

用显微镜对三轴实验后,稻田花岗岩样品内部的微破裂进行了观测.对比了两类三轴实验的微破裂,一种是常规三轴实验(1>2=3),另一种是真三轴实验(1>2>3).观测表明,茂木清夫在真三轴实验中发现的脆性岩石体积膨胀的各向异性效应是由微破裂的各向异性引起的,从而阐明了这一效应的微观机制,并对中等主应力引起的膨胀和微破裂的各向异性效应在地震研究中的意义进行了初步探讨.      

龙马溪页岩弹性各向异性与矿物分布之间的关系探讨

岩石弹性各向异性特征是普遍存在的,但导致岩石弹性各向异性的原因复杂且仍然存在一定争议.本研究以龙马溪页岩为例,试图建立页岩弹性各向异性和矿物分布之间关系.研究使用无损超声波探测获取岩石弹性各向异性参数,并使用背散射技术获取岩石矿物分布特征.研究通过引进变异系数来描述矿物或孔裂隙在不同方向的差异性,并通过2个正交方向的变异系数建立微观异质性指标,用于表征微观尺度上矿物或孔裂隙在不同方向的差异程度.微观背散射图像显示,龙马溪组页岩主要由石英和伊利石构成,且这两种矿物和孔裂隙在X和Y面上都有定向排列特征.相应地,它们的变异系数在X面和Y面上均表现出随角度增大而降低的特征;而在Z面,石英、伊利石和孔裂隙变异系数变化不明显,这与Z面上矿物和孔裂隙无明显方向性的特点一致.无损超声波探测结果显示,波速在X和Y面上随角度增加而减小,这与主要矿物和孔裂隙变异系数变化趋势相同;而在Z面,波速变化不大,与主要矿物和孔裂隙变异系数变化不明显的特征一致.以上观测结果说明,宏观波速与矿物的微观变异系数明显相关,暗示岩石弹性各向异性与矿物分布直接相关.     

The influence of microcrack density on the elastic and fracture mechanical properties of Stripa granite

Microcracks can have a strong influence on the elastic and fracture mechanical properties of rocks if they are numerous, or if they are orientated in unfavourable directions in anisotropic rocks in particular. This paper presents results from a great number of mechanical tests on Stripa granite containing various amounts of microcracks. Variations in the microcrack density were obtained by shock-heating the rock at different temperatures in the range 100–600°C for 3 h.The results presented are obtained from sound velocity measurements, uniaxial compression tests, Brazilian tests and three-point bending tests. The density of microcracks in the heated rock is studied by means of optical microscopy, SEM and differential strain analysis (DSA).Some of the calculated parameters show a maximum value for specimens heat-treated at about 100°C. The tensile strength is, for instance, substantially higher for specimens shock-heated at 100°C than for non-heated ones. Another striking feature is the initial decrease of the diameter observed in specimens heat-treated at 600°C when loaded in uniaxial compression. Both optical microscopy and DSA experiments reveal a large increase in microcracking when the heat-treatment temperature exceeds 300°C.     

岩石正交各向异性的实验观测

岩石速度各向异性由多种因素引起,如果这些因素在岩石中构成正交分布,就会出现正交各向异性.在实验室通过选择了三种不同岩性的岩石,采用超声脉冲测试方法中的横波偏振扫描方式,观测了岩石样品三个坐标方向上的横波速度变化,通过分析三个传播方向上出现横波分裂的差异,以及出现最大或最小快慢横波振幅的偏振方向发现,可从横波偏振方向和速度值判断出纹理砂岩和裂隙大理岩样品存在不同程度的速度正交各向异性和裂隙方向等特征,认为这两种岩石的正交特性的原因是岩石的层理或裂隙与岩石本身的结构(定向颗粒或结晶)构成的.实验室超声波测试可以找出正交各向异性的结构方向.     

Two species of microcracks

We identify two interrelated but independent species of microcracks with differentorigins and different distributions. One species is the classic high-stress microcracksidentified in laboratory stress-cells associated with acoustic emissions as microcracks openwith increasing stress. The other species is the low-stress distributions of closely-spacedstress-aligned fluid-saturated microcracks that observations of shear-wave splitting （SWS）demonstrate pervade almost all in situ rocks in the upper crust, the lower crust, and theuppermost 400 km of the mantle. On some occasions these two sets of microcracks may beinterrelated and similar （hence ＇species＇） but they typically have fundamentally-differentproperties, different distributions, and different implications. The importance for hydrocarbonexploration and recovery is that SWS in hydrocarbon reservoirs monitors crack alignmentsand preferred directions of fluid-flow. The importance for earthquake seismology is that SWSabove small earthquakes monitors the effects of increasing stress on the pervasive low-stressmicrocrack distributions so that stress-accumulation before, possibly distant, earthquakes canbe recognised and impendin~ earthquakes stress-forecast.     

Fracture toughness and subcritical crack growth during high-temperature tensile deformation of Westerly granite and Black gabbro

The double torsion testing method has been used to determine catastrophic and subcritical crack propagation parameters for pre-cracked specimens of Westerly granite and Black gabbro under a number of environmental conditions.The critical stress intensity factor for catastrophic crack propagation (fracture toughness) of granite and gabbro has been measured at temperatures from 20 to 400°C, in a vacuum. At 20°C, the fracture toughness of Westerly granite was $\text{1.79 ± 0.02}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, and for two blocks of Black gabbro it was $\text{3.03 ± 0.08}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ and $\text{2.71 ± 0.15}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, respectively. These values are very close to those reported by other investigators for tests conducted in air of ambient humidity at room temperature. For both rocks, fracture toughness at first increased slightly, and then decreased steadily on raising the temperature above ambient conditions. This behaviour is explained in terms of the density and distribution of thermally induced microcracks, as determined by quantitative optical microscopy.Subcritical crack growth behaviour has been studied at temperatures up to 300°C, and under water vapour at pressures of 0.6 to 15 kPa. Both the load relaxation and incremental constant displacement rate forms of the double torsion testing method were utilised to generate stress intensity factor/crack velocity diagrams. Crack growth was measured over the velocity range 5 × 10^?3 to 10^?7 m · s^?1. Increasing both temperature and water vapour pressure resulted in substantially higher crack growth rates. The overall effect of raising the temperature over the range studied here (20–300°C) was to increase the crack growth rate in granite and gabbro by ～5 and 7 orders of magnitude, respectively, at constant stress intensity factor and vapour pressure of water. For both rocks, the slopes of stress intensity factor/crack velocity curves were sensitive to changes in both temperature and water vapour pressure at low values of the latter parameter. Slopes fell substantially on raising the water vapour pressure, but were relatively insensitive to changes in temperature at these higher pressures. No subcritical crack growth limit was encountered.Estimates of the uncertainty in our experimental data are given. From the results of multiple load relaxation experiments on Westerly granite specimens, we estimate the uncertainty in position of stress intensity factor/crack velocity curves along the stress intensity axis to be c. 10% of the fracture toughness, and the uncertainty in slope of such curves to be c. 12%.Problems associated with the extrapolation of our experimental data to regions of higher effective confining pressure in the Earth's crust are discussed.     

The effect of cracks on the thermal expansion of rocks

Thermal expansion during the first heating cycle at atmospheric pressure was measured in several directions in seven igneous rocks between 25° and 400°C at slow heating rates. The coefficient of thermal expansion measured under these conditions increases more rapidly as temperature is increased than the average thermal expansion coefficient of the constituent minerals. The “extra” expansion is attributed to the formation of cracks by differential expansion of mineral grains. The presence of such cracks in the rocks during the cooling part of the cycle and during any subsequent heating and cooling cycles will result in a substantial decrease in the coefficient of thermal expansion as compared to that measured during the first heating cycles. The effect of cracks initially present in a rock was studied by measuring the full tensor of the coefficient of thermal expansion on two rocks with anisotropic crack distributions. In these two rocks the coefficient of thermal expansion is least in the direction perpendicular to the plane of greatest crack concentration. The implication of our data is that thermal expansion depends greatly on the fracture state of the rock. Both the fractures in the rock and the boundary conditions on the rock are significant for the interpretation of thermal expansion measurements and for their application to other problems.     

A Model for the Mechanical Behaviour of Bentheim Sandstone in the Brittle Regime

— The mechanical behaviour of Bentheim sandstone, a homogeneous quartz-rich sandstone with porosity of 22.8%, was investigated by triaxial compression tests conducted on dry samples. At confining pressures up to 35 MPa, the failure mode was characterized by a typical brittle deformation regime, as the samples showed dilatancy and failed by strain softening and brittle faulting. Previous studies have shown that the mechanical behaviour and failure mode of brittle porous granular rocks are governed by the time-dependent growth of microcracks. We analyse this process using the “Pore Crack Model” based on fracture mechanics analysis. It is consistent with the microstructure of porous granular rocks since it considers the growth of axial cracks from cylindrical holes in two dimensions. These cracks grow when their stress intensity factors reach the subcritical crack growth limit. Interaction between neighbouring cracks is introduced by calculating the stress intensity factor as the sum of two terms: a component for an isolated crack and an interaction term computed using the method of successive approximations. It depends on crack length, pore radius, pore density, and applied stresses. The simulation of crack growth from cylindrical holes, associated with a failure criterion based on the coalescence of interacting cracks, is used to compare the theoretical stress at the onset of dilatancy and at macroscopic rupture to the experimental determined values. Our approach gives theoretical results in good agreement with experimental data when microstructural parameters consistent with observations are introduced.     

Crustal shear-wave splitting in the epicentral zone of the 2001 Mw 7.7 Bhuj earthquake, Gujarat, India

We report here crustal shear-wave anisotropy, ranging from 1% to 10.76% with an average of 2.4% in the aftershock zone of the 2001 Bhuj earthquake, Gujarat, India, from a study of leading shear-wave polarization directions (LPSDs), which vary on average from NNW–SSE to E–W with a delay of 0.07–0.14 s. The delays in the NNW–SSE to NE–SW directions observed at seven stations, near the seismogenic fault, suggest cracks parallel to the direction of the maximum horizontal regional compressional stress prevailing in the region, suggesting a dilatancy-induced anisotropy resulting from approximately stress-aligned parallel vertical micro-cracks. In contrast, the LPSDs at Ramvav, Rapar and Vondh stations, away from the seismogenic fault, are fault parallel, approximately E–W and almost orthogonal to the stress-aligned polarizations inferred elsewhere. The maximum average time delay of 0.14 s is observed at Lodai, where the fast polarization direction is found to be N338°W. This has been observed from anisotropic poro-elastic (APE) modelling and observations that these are 90° flips in shear-wave polarization, resulting from propagation through micro-cracks containing fluids at critically high pore-fluid pressure surrounding the hypocenter of the 2001 mainshock. The presence of high pore-fluid pressure in the seismogenic fault zone could also explain the observed scatter in shear-wave time delays. Further, the coincidence of the N–S trending intrusive bodies (as inferred from tomographic studies in the area) with the N–S direction of regional maximum horizontal compressional stress supports the interpretation of stress-aligned vertical extensive-dilatant anisotropic (EDA) cracks. The depth distribution of the estimated anisotropy (1–10.76%), b-values and stress drop values suggests an increase at 18–30 km depths, which could be attributed to high pore-fluid pressures resulting from a fluid-filled fractured rock matrix or open micro-cracks (characterized by high crack density and high porosity) coinciding with a low velocity zone (at 18–30 km depths) as delineated from tomographic studies in the area.     

裂缝型储层流体识别方法

裂缝型储层的描述包括预测裂缝分布特征和识别裂隙充填物.依据等效介质理论计算的纵波速度随裂缝密度的增大而减小.正演地震记录显示,裂缝介质含气时反射振幅最大,且变化程度比含油或含水时大.叠前方位AVO反演所得的各向异性梯度Bani与裂缝密度成正比,可用于描述有效裂缝发育强度.对于不同的裂缝密度,各向异性梯度Bani与各向同性梯度Biso的比值I(1)fluid近似为常数,且对流体敏感.经裂缝纵横比和背景介质拉梅常数修正后,流体因子Ifluid既不随纵横比变化,又不受背景介质的影响,是裂缝型储层敏感的流体识别因子.在塔里木盆地塔北哈拉哈塘地区热瓦普区块碳酸盐岩储层裂缝发育区域,运用该参数在井点处的流体识别效果与钻井结果一致.     

小江断裂带北段地壳浅层P波速度各向 异性观测的原理与数值试验

以二维情形下观测速度场为各向同性场和各向异性场的叠加为前提, 提出了一种利用走时残差估算地震波速度各向异性的方法, 即剩余慢度矢量法. 利用小江断裂带北段巧家流动地震台阵24个台站记录的3181次地震事件的P波走时残差, 采用剩余慢度矢量法计算了各观测台站周围水平方向上尺度为0.5°×0.5°, 震源深度为0—5 km的剩余慢度矢量, 由此得到了P波快波和慢波方向. 计算结果表明, 大部分观测台站周围的P波速度方向性较为一致, 快波方向为ESE向, 慢波方向为NNE向. 快波方向与小江断裂带北段应力场P轴方向较为一致, 而慢波方向与应力场T轴方向一致, 表明应力的长期作用可能是导致P波速度各向异性的重要原因.      

低孔隙度泥页岩应力依赖的各向异性裂纹演化特性研究

本文在实验室对所获取的东营地区层理发育的低孔隙度页岩和泥岩的各向异性裂纹演化特性进行了研究,获得了各向同性条件下泥页岩的力学与超声波响应特性,分析了应力幅度对于页岩声波速度和各向异性的影响.主要结论包括:(1)泥页岩在循环载荷下存在滞后效应,表明其经历了去压实或油气产生导致的超压;(2)泥岩和页岩具有不同程度的各向异性,随着各向同性压力的增高微裂隙逐渐闭合,样品的各向异性程度减弱;(3)分析了岩石韧度和裂纹损伤参数随压力的变化特征,相比泥岩,页岩各向异性程度更高,随压力变化更明显,其裂纹导致的附加各向异性更强;(4)分析了各向异性岩石的动态弹性模量特征,由于软裂隙空间的闭合,动态弹性模量在低压条件下都随着围压的增加有硬化趋势.     

TTI介质弹性波相速度与偏振特征分析

相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好.     

工程结构三维疲劳裂纹最大应力强度因子计算

介绍了裂纹的类型、裂纹尖端应力场的奇异性。以一维问题为例,推导论证了奇异单元能够很好的反映裂纹尖端应力场的奇异性。应力强度因子一般表达式表明应力强度因子与载荷呈线性关系,并依赖于物体和裂纹的几何形状和尺寸。本文借助大型通用有限元软件ANSYS,采用位移外插法计算了三维表面裂纹前沿不同位置处的应力强度因子,并与《应力强度因子手册》基于实验的理论公式计算结果相比较。结果表明:有限元结果与理论解误差较小,裂纹最深处应力强度因子最大。     

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.     

The Effect of Varying Damage History in Crystalline Rocks on the P- and S-Wave Velocity under Hydrostatic Confining Pressure

Cracks play a very important role in many geotechnical issues and in a number of processes in the Earth’s crust. Elastic waves can be used as a remote sensing tool for determining crack density. The effect of varying crack density in crystalline rock on the P- and S-wave velocity and dynamic elastic properties under confining pressure has been quantified. The evolution of P- and S-wave velocity were monitored as a suite of dry Westerly granite samples were taken to 60, 70, 80 and 90 % of the unconfined uniaxial strength of the sample. The damaged samples were then subjected to hydrostatic confining pressure from 2 MPa to 200 MPa to quantify the effect of varying crack density on the P- and S-wave velocity and elastic properties under confining pressure. The opening and propagation of microcracks predominantly parallel to the loading direction during uniaxial loading caused a 0.5 and 6.3 % decrease in the P- and S-wave velocity, respectively. During hydrostatic loading, microcracks are closed at 130 MPa confining pressure. At lower pressures the amount of crack damage in the samples has a small but measureable effect. We observed a systematic 6 and 4 % reduction in P- and S-wave velocity, respectively, due to an increase in the fracture density at 2 MPa confining pressure. The overall reduction in the P- and S-wave velocity decreased to 2 and 1 %, respectively, at 50 MPa. The elastic wave velocities of samples that have a greater amount of microcrack damage are more sensitive to pressure. Effective medium modelling was used to invert elastic wave velocities and infer crack density evolution. Comparing the crack density results with experimental data on Westerly granite samples shows that the effective medium modelling used gave interpretable and reasonable results. Changes in crack density can be interpreted as closure or opening of cracks and crack growth.     

含裂隙介质中的视电阻率各向异性变化

我国50多年的视电阻率连续观测结果表明,大地震前近震中区域的视电阻率呈现出与主压应力方位有关的各向异性变化,即:垂直于主压应力方向观测的变化幅度最大,平行方向最小或不明显,斜交方向介于二者之间.目前我国定点台站视电阻率观测的探测范围主要在浅层沉积层以内,通常含有较多的含水裂隙.本文将地下岩土介质简化为由固体基质和含流体/气体裂隙组成的固液气三相介质,且基质、流体和气体具有标量形式的电阻率,推导出了包含基质和流体电阻率、裂隙率、饱和度和裂隙面积率因子的电阻率张量表达式.以裂隙的扩展/闭合表示应力作用下裂隙的变化,得到了电阻率随裂隙变化的微分形式,电阻率变化对裂隙体积变化放大系数的表达式和裂隙横向变化对纵向电阻率影响的横向权系数的表达式.在此基础上得到了介质电阻率和视电阻率的各向异性变化特征:对于含水裂隙介质,无论裂隙如何变化,均是最小主轴方向电阻率的变化幅度大于其他方向;对于含水孔隙介质,沿孔隙主要变化方向的主轴电阻率变化幅度大于其他方向.对于各向异性变化,视电阻率和介质电阻率存在π/2的方向差异.相较于含水岩石,无水岩石介质电阻率的各向异性变化不显著.本文提出的电阻率表达式可以对实验室和野外实际观测的许多结果做出合理的解释.     

Determination of directions of horizontal principal stress and identification of active faults in Kachchh (India) by electromagnetic radiation (EMR)

The Kachchh basin in the western India is known for its recent high seismicity. This study presents an application of the geogenic Electromagnetic Radiation (EMR) technique for deciphering the directions of principal horizontal stress in the eastern Kachchh. The principal direction of horizontal stress obtained from EMR differs from those obtained from earthquake focal plane solutions. The major horizontal principal stress based on the EMR study shows an azimuth of N60°E ± 10°. The principal directions of EMR emissions are parallel to the acute bisector of conjugate microcracks. The azimuthal distribution of EMR signal and dimension of microcracks suggest that the EM emissions are transversely polarized.The study also deals with the first application of electromagnetic radiation emissions to identify active fracture planes in sandstones that could become potential active faults later, which might be seismogenic or nonseimogenic. This study is based on linear profiling at six different places across two major faults, the Kachchh Mainland Fault (KMF) and the South Wagad Fault (SWF) in the eastern Kachchh. Anomalously, high EMR emissions are observed in the eastern part of the KMF, indicating active surface deformation.