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.     

含流体致密砂岩的纵波频散及衰减:基于双重双重孔隙结构模型描述的特征分析

致密砂岩普遍具有低孔、低渗及微裂缝发育的地质特征,并且呈现出很强的非均匀性.致密砂岩储层与常规砂岩储层比较,具有明显的岩石物理性质、渗流力学性质方面的差异.致密砂岩内部的非均匀性对弹性波频散、衰减有显著影响,其中包括孔隙结构的非均匀性,即岩石内部孔隙参数的不均一性,以及孔隙内部不相混溶流体的非均匀分布;此外,非均匀性的尺度也决定了波出现显著频散与衰减的频段.综合考虑致密砂岩孔隙结构非均匀性及流体斑块状饱和的非均匀性,本文采用双双重孔隙介质结构模拟了致密砂岩的弹性波响应,分析了同时具备两类非均质性岩石中的波传播特征.调查分析了两组分别来自中国鄂尔多斯盆地苏里格气田及四川盆地广安气田的不同类型致密砂岩储层的岩芯超声波实验数据,给出了岩石样本的弹性波速度频散与衰减曲线.结果显示理论模型预测结果与完全饱和、部分饱和岩石的实验数据吻合良好.对两个地区致密砂岩岩芯数据进行对比分析,苏里格致密砂岩样本总体上比广安致密砂岩渗透率高,在各孔隙度范围内,特征模拟显示苏里格样本的裂隙尺寸明显大于广安样本.广安致密砂岩在低孔隙度范围内发育了更多、更小的颗粒裂隙/接触.致密砂岩的速度频散与衰减结果受流体黏度、晶体破裂及流体斑块状饱和的共同影响.此外,孔隙度越大,部分饱和岩石中斑块状饱和机制对总衰减的贡献越低,与之相对,结构非均质性所占的比重则有所增强.     

Simulation of full-waveform log in saturated cracked formations using Hudson's approach

We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.     

有效应力变化对致密砂岩孔隙结构及弹性波响应的影响规律

致密砂岩气藏具有裂缝发育和有效应力高的特征,研究不同有效压力下孔、裂隙介质地震波传播特征,有利于地震解释与地下储层的识别.但是前人的研究较少考虑岩石内部微观孔隙结构特征与孔隙、裂隙间流体流动的关系.本文首先通过选取四川盆地典型致密砂岩岩样,在不同有效压力下对岩石样本进行超声波实验测量.然后基于实验测得的纵、横波速度进行裂隙参数反演,得到不同有效压力下致密砂岩样本的裂隙孔隙度.再将裂隙孔隙度和样本岩石物理参数代入双重孔隙介质模型,模拟得到不同有效压力下饱水致密砂岩样本纵横波速度频散和衰减的变化规律.结果表明模型预测的速度频散曲线与纵波速度实验测量结果能够较好的吻合.最后统计分析了致密砂岩裂隙参数,得到了致密砂岩储层裂隙参数随有效压力及孔隙度变化特征.依据实际岩石物理参数建立模型,其裂隙参数三维拟合结果能够较好描述致密砂岩裂隙结构与孔隙度、应力的关联,可为实际地震勘探中预测储层裂缝性质提供基础依据.     

Comparison of stress-dependent geophysical,hydraulic and mechanical properties of synthetic and natural sandstones for reservoir characterization and monitoring studies

Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi-physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P- and S-wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny-shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies.     

INDUCED POLARIZATION,A METHOD TO STUDY WATER-COLLECTING PROPERTIES OF ROCKS*

During the last decade many hypotheses were suggested to explain the phenomenon of induced electrical polarization in ionic conductive media. The most reliable of these is Fredricksberg's. Fredricksberg (1962) supposed that the pore spaces of a rock is composed of successively narrow (active zones) and wide (inactive zones). He simulated these pore spaces by a synthetic material that has an extremely high resistance. The pore spaces were generally in tube forms which exhibited some constrictions. He saturated these tubes with an electrolyte of a given concentration. An electric current was passed through this model. He observed an induced polarization voltage after current interruption. He attributed the formation of this voltage to a concentration gradient which took place due to the presence of excess charges in the active zones. Fredricksberg introduced a parameter (9) which described the relation among the lengths and cross sectional areas of the wide zones, the number of ions within each zone after current interruption with the recorded polarizability. The aim of this work is to correlate Fredricksberg's parameter with a parameter determined for natural rocks and to show experimentally the validity of this hypothesis when applying for some varieties of sandstones and volcanic rocks. The new parameter will help to evaluate a relationship between the polarizability and the water-collecting properties of rocks. Herein, we used the tortuosity T of sandstone samples instead of the parameter φ which was used by Fredricksberg to represent the pore geometry within his model (tortuosity of the passes within the model). It was shown that both φ and T have the same relationship with the polarizability ν of the rock samples and if φ or T have very low or very high values the polarizability ν tends to its minimum value, i.e. the curve representing the relation between ν and T has a maximum point corresponding to an intermediate value of T. This result supports Fredricksberg's hypothesis and confirm his results on synthetic models. For volcanic rocks the formation factor F was used since it was difficult to determine the porosity of the samples and consequently to calculate the tortuosity T as for sandstone samples. Experimental results confirm those obtained from sandstone. The grain constituents of sandstone samples were represented on equilateral triangle and the magnitude of induced polarization ν of each sample was deduced and represented on this triangle. Equipolarizability values ν drawn on this triangle showed that TJ will increase as the silty fractions of the rock increase, where the center of this triangle (represents minimum porosity) has polarizability less than 0.25%. An attempt was made to determine the coefficient of anisotropy of volcanic rock samples using the induced polarization method. For this reason the polarizability was deduced by measuring the induced polarization voltage for two perpendicular directions in a fractured cubes of andesitic basalt samples the coefficient of anisotropy was found to be equal 1.18.     

多重孔岩石微分等效介质模型及其干燥情形下的解析近似式

经典的微分等效介质(DEM)理论可用于确定多孔介质的弹性性质,但由于缺乏多重孔DEM方程,其估计的多重孔岩石的等效弹性模量依赖于包裹体(即不同孔隙纵横比的孔或缝)的添加顺序.本文首先从Kuster-Toksöz理论出发建立了Zimmermann和Norris两种形式的多重孔DEM方程.Norris形式的多重孔DEM方程预测的等效弹性模量总是位于Hashin-Shtrikman上下限内,而Zimmermann形式的多重孔DEM方程有时会越界.然后,通过使用干燥岩石模量比的解析近似式,对两个相互耦合的Norris形式DEM方程进行解耦得到干燥多重孔岩石的体积和剪切模量解析式.用全DEM方程的数值解对解析近似式的有效性进行了测试,解析公式的计算结果在整个孔隙度分布区间与数值解吻合良好.对实验室测量数据在假设岩石含有双重孔隙的情形下用双重孔DEM解析公式对岩石的弹性模量进行了预测,结果表明,解析式准确地预测了弹性模量随孔隙度的变化.双重孔(即软、硬孔)DEM解析模型可用来反演各孔隙类型的孔隙体积比,它可以通过实验室测量与理论预测之间的平方误差最小反演得到.砂岩样品的反演结果揭示,软孔的孔隙体积百分比与粘土含量没有明显的相关性.     

Velocity–porosity relationships: Predictive velocity model for cemented sands composed of multiple mineral phases

Computer simulations are used to calculate the elastic properties of model cemented sandstones composed of two or more mineral phases. Two idealized models are considered – a grain‐overlap clay/quartz mix and a pore‐lining clay/quartz mix. Unlike experimental data, the numerical data exhibit little noise yet cover a wide range of quartz/cement ratios and porosities. The results of the computations are in good agreement with experimental data for clay‐bearing consolidated sandstones. The effective modulus of solid mineral mixtures is found to be relatively insensitive to microstructural detail. It is shown that the Hashin–Shtrikman average is a good estimate for the modulus of the solid mineral mixtures. The distribution of the cement phase is found to have little effect on the computed modulus–porosity relationships. Numerical data for dry and saturated states confirm that Gassmann's equations remain valid for porous materials composed of multiple solid constituents. As noted previously, the Krief relationship successfully describes the porosity dependence of the dry shear modulus, and a recent empirical relationship provides a good estimate for the dry‐rock Poisson's ratio. From the numerical computations, a new empirical model, which requires only a knowledge of system mineralogy, is proposed for the modulus–porosity relationship of isotropic dry or fluid‐saturated porous materials composed of multiple solid constituents. Comparisons with experimental data for clean and shaly sandstones and computations for more complex, three‐mineral (quartz/dolomite/clay) systems show good agreement with the proposed model over a very wide range of porosities.     

基于数字岩心的含裂隙储层砂岩介电性质研究

了解储层岩石的介电特性在石油工业的各个方面都有重要的应用.小尺度裂隙是影响岩石介电性质的地质因素之一,获得裂隙对含裂隙岩石介电性质影响的定量关系具有重要的理论和实践意义.以含裂隙人造砂岩的三维微观数字结构为基础,通过基于三维有限差分算法计算的岩石介电性质与实验数据的对比验证数值计算方法的有效性.在此基础上,通过理论模型获得不同孔隙度基质的介电性质,并在不含裂隙人造砂岩的三维微观数字结构中人为添加以裂隙密度和纵横比为定量表征参数的裂隙,应用验证后的数值算法模拟随频率变化的含裂隙砂岩的介电性质,分析和研究不同孔隙度基质中定向排列裂隙对砂岩介电性质的影响.结果表明,当裂隙孔隙度随裂隙纵横比或裂隙密度发生改变时,含裂隙砂岩的介电性质与裂隙密度以及裂隙纵横比呈正相关关系,而当裂隙孔隙度保持不变时,含裂隙砂岩的介电性质随裂隙纵横比的减小而增大;裂隙参数的改变对不同基质孔隙度的含裂隙砂岩的介电性质的影响趋势较为一致,但随着基质孔隙度的减小,裂隙对砂岩介电性质的影响逐渐增大.裂隙参数和基质孔隙度对含裂隙砂岩介电性质影响的研究结果为基于介电特性的裂缝性油气储层的定量表征提供了依据,在油气勘探开发中具有重要的应用前景.     

The influence of local fluid flow and the microstructure on elastic and anelastic rock properties

The frequency dependent mechanism of local fluid flow was found to be the decisive absorption and dispersion mechanism in fluid containing sandstones. In the ultrasonic frequency range local fluid flow and grain surface effects control the behaviour of highly porous and highly permeable rock if a pore fluid is present. Both mechanisms depend less on macroscopic rock parameters like porosity and permeability than essentially on microscopic parameters like crack size, crack density and grain contact properties. To demonstrate directly the important influence of the microstructure on the rock elastic and anelastic properties the microstructure of a sandstone was artificially changed by thermal cracking. The cracked rock exhibits a clearly changed behaviour at low uniaxial as well as at high hydrostatic pressure despite small changes of porosity and permeability. Fluid effects increase due to cracking. The experimental results are explained by means of a rock, model and local fluid flow. These results emphasize that it is the microstructure which controls the elastic and anelastic rock behaviour, even at high hydrostatic pressure.     

Elastic properties of double-porosity rocks using the differential effective medium model

An approach to determining the effective elastic moduli of rocks with double porosity is presented. The double‐porosity medium is considered to be a heterogeneous material composed of a homogeneous matrix with primary pores and inclusions that represent secondary pores. Fluid flows in the primary‐pore system and between primary and secondary pores are neglected because of the low permeability of the primary porosity. The prediction of the effective elastic moduli consists of two steps. Firstly, we calculate the effective elastic properties of the matrix with the primary small‐scale pores (matrix homogenization). The porous matrix is then treated as a homogeneous isotropic host in which the large‐scale secondary pores are embedded. To calculate the effective elastic moduli at each step, we use the differential effective medium (DEM) approach. The constituents of this composite medium – primary pores and secondary pores – are approximated by ellipsoidal or spheroidal inclusions with corresponding aspect ratios. We have applied this technique in order to compute the effective elastic properties for a model with randomly orientated inclusions (an isotropic medium) and aligned inclusions (a transversely isotropic medium). Using the special tensor basis, the solution of the one‐particle problem with transversely isotropic host was obtained in explicit form. The direct application of the DEM method for fluid‐saturated pores does not account for fluid displacement in pore systems, and corresponds to a model with isolated pores or the high‐frequency range of acoustic waves. For the interconnected secondary pores, we have calculated the elastic moduli for the dry inclusions and then applied Gassmann's tensor relationships. The simulation of the effective elastic characteristic demonstrated that the fluid flow between the connected secondary pores has a significant influence only in porous rocks containing cracks (flattened ellipsoids). For pore shapes that are close to spherical, the relative difference between the elastic velocities determined by the DEM method and by the DEM method with Gassmann's corrections does not exceed 2%. Examples of the calculation of elastic moduli for water‐saturated dolomite with both isolated and interconnected secondary pores are presented. The simulations were verified by comparison with published experimental data.     

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.     

Predicting horizontal velocities from well data

The Imperial College borehole test site consists of four boreholes with depths lying between 260 and 280 m. The boreholes intersect several cyclical sequences of sandstones, mudstones and limestones. The formations are highly laminated and ultrasonic measurements on preserved core have shown that the mudstones are intrinsically anisotropic. Little or no anisotropy is associated with the sandstones and limestones. A scheme is proposed to predict synthetic vertical and horizontal P- and S-wave logs. Combining (an)isotropic effective medium theories, the Gassmann equation and Backus averaging, the scheme extends previous sand-shale models to transversely isotropic rock formations. The model assumes that the anisotropy is due to layering and due to the preferred horizontal orientation of the clay minerals, pores and cracks within the mudstones. The pores and cracks within the sandstones and limestones are randomly orientated. After fitting the model to the ultrasonic data to obtain the unknown parameters, the model successfully predicts the sonic log and the direct arrival times from a cross-hole survey.     

基于等效介质理论的孔隙纵横比分布反演

孔隙纵横比是描述多孔岩石微观孔隙结构特征的重要参数,目前用于获取岩石完整孔隙纵横比分布的经典模型为David-Zimmerman(D-Z)孔隙结构模型,该模型假设岩石由固体矿物基质、一组纵横比相等的硬孔隙以及多组纵横比不等的微裂隙构成,并认为固体矿物基质和硬孔隙均不受压力影响,在此基础上,利用超声纵横波速度的压力依赖性...     

Origin of gray-green sandstone in ore bed of sandstone type uranium deposit in north Ordos Basin

Dongsheng sandstone-type uranium deposit is located in the northern part of Ordos Basin, occurring in the transitional zones between gray-green and gray sandstones of Jurassic Zhiluo Formation. Sandstones in oxidized zone of the ore bed look gray-green, being of unique signature and different from one of ordinary inter-layered oxidation zone of sandstone-type uranium deposits. The character and origin of gray-green sandstones are systematically studied through their petrology, mineralogy and geochemistry. It is pointed out that this color of sandstones is originated from secondary oil-gas reduction processes after paleo-oxidation, being due to acicular-leaf chlorite covering surfaces of the sandstone grains. To find out the origin of gray-green sandstone and recognize paleo-oxidation zones in the ore bed are of not only theoretical significance for understanding metallogenesis of this kind of sandstone-type uranium deposit, but also very importantly practical significance for prospecting for similar kind of sandstone-type uranium deposit.     

Origin of gray-green sandstone in ore bed of sandstone type uranium deposit in north Ordos Basin

Dongsheng sandstone-type uranium deposit is located in the northern part of Ordos Basin, occurring in the transitional zones between gray-green and gray sandstones of Jurassic Zhiluo Formation. Sandstones in oxidized zone of the ore bed look gray-green, being of unique signature and different from one of ordinary inter-layered oxidation zone of sandstone-type uranium deposits. The character and origin of gray-green sandstones are systematically studied through their petrology, mineralogy and geochemistry. It is pointed out that this color of sandstones is originated from secondary oil-gas reduction processes after paleo-oxidation, being due to acicular-leaf chlorite covering surfaces of the sandstone grains. To find out the origin of gray-green sandstone and recognize paleo-oxidation zones in the ore bed are of not only theoretical significance for understanding metallogenesis of this kind of sandstone-type uranium deposit, but also very importantly practical significance for prospecting for similar kind of sandstone-type uranium deposit.     

Seismoelectric measurements in a porous quartz‐sand sample with anisotropic permeability

Seismoelectric coupling coefficients are difficult to predict theoretically because they depend on a large numbers of rock properties, including porosity, permeability, tortuosity, etc. The dependence of the coupling coefficient on rock properties such as permeability requires experimental data. In this study, we carry out a set of laboratory measurements to determine the dependence of seismoelectric coupling coefficient on permeability. We use both an artificial porous “sandstone” sample, with cracks, built using quartz‐sand and Berea sandstone samples. The artificial sample is a cube with 39% porosity. Its permeability levels are anisotropic: 14.7 D, 13.8 D, and 8.3 D in the x‐, y‐, and z‐directions, respectively. Seismoelectric measurements are performed in a water tank in the frequency range of 20 kHz–90 kHz. A piezoelectric P‐wave source is used to generate an acoustic wave that propagates through the sample from the three different (x, y, and z) directions. The amplitudes of the seismoelectric signal induced by the acoustic waves vary with the direction. The highest signal is in the direction of the highest permeability, and the lowest signal is in the direction of the lowest permeability. Since the porosity of the sample is constant, the results directly show the dependence of seismoelectric coefficients on permeability. Seismoelectric measurements with natural rocks are performed using Berea sandstone 500 and 100 samples. Because the Berea samples are nearly isotropic in permeability, the amplitudes of the seismoelectric signals induced in the different directions are the same within the measurement error. Because the permeability of Berea 500 is higher than that of Berea 100, the amplitude of the seismoelectric signals induced in Berea 500 is higher than those in Berea 100. To determine the relative contributions of porosity and permeability on seismoelectric conversion, we carried out an analysis, using Pride (1994) formulation and Kozeny–Carman relationship; the normalized amplitudes of seismoelectric coupling coefficients in three directions are calculated and compared with the experimental results. The results show that the seismoelectric conversion is related to permeability in the frequency range of measurements. This is an encouraging result since it opens the possibility of determining the permeability of a formation from seismoelectric measurements.     

Failure Mode and Spatial Distribution of Damage in Rothbach Sandstone in the Brittle-ductile Transition

— To elucidate the spatial complexity of damage and evolution of localized failure in the transitional regime from brittle faulting to cataclastic ductile flow in a porous sandstone, we performed a series of triaxial compression experiments on Rothbach sandstone (20% porosity). Quantitative microstructural analysis and X-ray computed tomography (CT) imaging were conducted on deformed samples. Localized failure was observed in samples at effective pressures ranging from 5 MPa to 130 MPa. In the brittle faulting regime, dilating shear bands were observed. The CT images and stereological measurements reveal the geometric complexity and spatial heterogeneity of damage in the failed samples. In the transitional regime (at effective pressures between 45 MPa and 130 MPa), compacting shear bands at high angles and compaction bands perpendicular to the maximum compression direction were observed. The laboratory results suggest that these complex localized features can be pervasive in sandstone formations, not just limited to the very porous aeolian sandstone in which they were first documented. The microstructural observations are in qualitative agreement with theoretical predictions of bifurcation analyses, except for the occurrence of compaction bands in the sample deformed at effective pressure of 130 MPa. The bifurcation analysis with the constitutive model used in this paper is nonadequate to predict compaction band formation, may be due to the neglect of bedding anisotropy of the rock and multiple yield mechanisms in the constitutive model.     

Experimental verification of effective anisotropic crack theories in variable crack aspect ratio medium

Physical modelling of cracked/fractured media using downscaled laboratory experiments has been used with great success as a useful alternative for understanding the effect of anisotropy in the hydrocarbon reservoir characterization and in the crustal and mantle seismology. The main goal of this work was to experimentally verify the predictions of effective elastic parameters in anisotropic cracked media by Hudson and Eshelby–Cheng's effective medium models. For this purpose, we carried out ultrasonic measurements on synthetic anisotropic samples with low crack densities and different aspect ratios. Twelve samples were prepared with two different crack densities, 5% and 8%. Three samples for each crack density presented cracks with only one crack aspect ratio, whereas other three samples for each crack density presented cracks with three different aspect ratios in their composition. It results in samples with aspect ratio values varying from 0.13 to 0.26. All the cracked samples were simulated by penny‐shaped rubber inclusions in a homogeneous isotropic matrix made with epoxy resin. Moreover, an isotropic sample for reference was constructed with epoxy resin only. Regarding velocity predictions performed by the theoretical models, Eshelby–Cheng shows a better fit when compared with the experimental results for samples with single and mix crack aspect ratio (for both crack densities). From velocity values, our comparisons were also performed in terms of the ε, γ, and δ parameters (Thomsen parameters). The results show that Eshelby–Cheng effective medium model fits better with the measurements of ε and γ parameters for crack samples with only one type of crack aspect ratio.     

Small faults formed as deformation bands in sandstone

Small faults with displacements of a few millimeters or centimeters are abundant in the Entrada and Navajo Sandstones, in the San Rafael Desert, Utah, where they are important primary structures, preceding the development of large faults with displacements of several meters or tens of meters. The small faults contain no surfaces of discontinuity, rather they occur asdeformation bands about one millimeter and tens or hundreds of meters long, and across which the displacements are distributed. Two zones with different modes of deformation can be distinguished within a deformation band: an outer zone where the matrix, including pores and matrix material, deforms; and an inner zone, about 0.5 mm thick, where the sand grains fracture and further consolidation takes place. Fracturing of the grains is controlled by contact geometry; the grains tend to split into subgrains along lines connecting contact points between the grains. Then the angular subgrains, which are readily fractured, are further granulated and mixed with the matrix. The final product is the deformation band, with much smaller grain size, poorer sorting, and lower porosity than the original parent sandstone. The sandstone on either side of a deformation band is almost undisturbed-fractures are rare there — so that deformation is highly localized within the band. The material within a deformation band apparently strain hardens as a result of the deformation; perhaps this is why the shear displacement across a deformation band is at most a few centimeters.