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.     

P- and S-wave anisotropy of a synthetic sandstone with controlled crack geometryl1

A method is reported for the production of synthetic porous sandstones containing cracks of known dimensions and geometry with respect to the matrix. A synthetic sandstone was manufactured from Sand cemented with an epoxy glue. The cracks of known geometry were introduced into the material in the manufacturing stage, by emplacing thin metallic discs in the Sand-epoxy matrix. These discs were chemically leached out of the consolidated porous sandstone. Acoustic anisotropy. and shear-wave splitting were observed in the synthetic sandstones. For the dry sample the observed angular dependence of the P- and S- wave velocities (at 100 kHz) compares well, qualitatively, with the theoretical models of Hudson and of Thomsen. Quantitatively, however, the experimental data fits Hudson's model better. For the case of a saturated sample the experimental results are in excellent agreement with Thomsen's model. Hudson's model, on the other hand, predicts a different angular dependence for P-waves. This demonstrates that the concept of fluid transfer between cracks and the ambient porosity can be a significant process. The results reported here are from the first successful experiment in which the theoretical models were tested on a porous material containing a known crack geometry.     

The porosity of saturated shallow sediments from seismic compressional and shear wave velocities

This study devises a new analytical relationship to determine the porosity of water-saturated soils at shallow depth using seismic compressional and shear wave velocities. Seismic refraction surveys together with soil sample collection were performed in selected areas containing water-saturated clay–silt, sand and gravely soils. Classification of clay–silt, sand and gravel dense soils provided the coefficient of experimental equation between the data sets, namely, Poisson's ratio, shear modulus and porosity values. This study presents a new analytical relationship between Poisson's ratio and shear modulus values, which are obtained from seismic velocities and porosity values of water-saturated material computed from water content and grain densities, which are determined by laboratory analysis of disturbed samples. The analytical relationship between data sets indicates that when the shear modulus of water-saturated loose soil increases, porosity decreases logarithmically. If shear modulus increases in dense or solid saturated soils, porosity decreases linearly.     

A PHYSICAL MODEL STUDY OF SCATTERING OF WAVES BY ALIGNED CRACKS: COMPARISON BETWEEN EXPERIMENT AND THEORY1

An approximation to plane-wave propagation through a composite material is examined using a physical model with oriented but randomly distributed penny-shaped rubber inclusions within an isotropic epoxy resin matrix. A pulse transmission method is used to determine velocities of shear and compressional waves as a function of angle of incidence and crack density. The experimental and theoretical results of Hudson were compared and limitations within the crack parameters used in this study have been determined. Results from both polarized shear waves (S1, S2) compare favourably with the theory for a composite with up to 7% crack density, but theory and experiment diverge at higher crack densities. On the other hand, compressional-wave velocities at low crack densities (1% and 3%) compare favourably with the theory. It is also shown that the velocity ratio V_p/V_s for two extreme cases, i.e. propagation normal and parallel to the cracks, as a function of crack density and porosity, has a strong directional dependence.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

水饱和裂纹对地壳岩样中地震波速及各向异性的影响

选择４种地壳岩石样品,经干燥或水饱和处理后在不同围压条件下测量了在其中传播的纵、横波的速度及其各向异性．在大气压条件下低孔隙度(＜１％岩样中,水饱和样品中的纵波速度明显地比干燥样品中的高,但横波速度的差别不大．因为在低孔隙度岩样中纵波速度对孔隙流体的反应比横波速度敏感,可以用泊松比的变化来反映随着围压的增加晶粒间流体对弹性波传播特性的影响．根据实验数据,按Ｏ’Ｃｏｎｎｅｌｌ模型分别计算了干燥和水饱和岩样中的裂纹密度,与通过实测体应变曲线得到的裂纹孔隙度十分吻合．利用横波的速度和偏振特性可以推断岩样中定向排列微裂纹的空间取向情况．研究表明,同时测量在岩样中传播的纵、横波的速度,通过Ｖｐ／Ｖｓ比值可以给出有关颗粒边界流体的证据,也可以估计岩样中的裂纹密度．     

含流体砂岩地震波频散实验研究

为了研究孔隙流体对不同渗透率岩石地震波速度的影响,在实验室利用跨频带岩石弹性参数测试系统得到了应变幅值10-6的2~2000Hz频段下的地震波速度和1 MHz频率下的超声波速度,利用差分共振声谱法得到了频率600Hz岩石干燥和完全饱水情况下岩石声学参数.实验表明,在低饱和度下,致密砂岩在地震和超声频段下没有明显的频散;在高饱和度下纵波速度的频散变得明显.从干燥到完全水饱和条件,不同频率测量的致密砂岩的体积模量随岩石孔隙度增高而降低,且体积模量的变化量受岩石微观孔隙结构的影响较大.高孔、高渗砂岩无论在低含水度下还是在高含水饱和度下频散微弱,并且在地震频段下围压对于岩石纵横波速度的影响要大于频率的影响.高孔、高渗砂岩和致密砂岩不同含水饱和度下的频散差异可应用于储层预测,油气检测等方面,同时该研究可以更好地帮助理解岩石的黏弹性行为,促进岩石物理频散理论的发展,提高地震解释的精度.     

Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory

Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given.     

Effects of pore fluids on quasi-static shear modulus caused by pore-scale interfacial phenomena

It is evident from the laboratory experiments that shear moduli of different porous isotropic rocks may show softening behaviour upon saturation. The shear softening means that the shear modulus of dry samples is higher than of saturated samples. Shear softening was observed both at low (seismic) frequencies and high (ultrasonic) frequencies. Shear softening is stronger at seismic frequencies than at ultrasonic frequencies, where the softening is compensated by hardening due to unrelaxed squirt flow. It contradicts to Gassmann's theory suggesting that the relaxed shear modulus of isotropic rock should not depend upon fluid saturation, provided that no chemical reaction between the solid frame and the pore fluid. Several researchers demonstrated that the shear softening effect is reversible during re-saturation of rock samples, suggesting no permanent chemical reaction between the solid frame and the pore fluid. Therefore, it is extremely difficult to explain this fluid–rock interaction mechanism theoretically, because it does not contradict to the assumptions of Gassmann's theory, but contradicts to its conclusions. We argue that the observed shear softening of partially saturated rocks by different pore fluids is related to pore-scale interfacial phenomena effects, typically neglected by the rock physics models. These interface phenomena effects are dependent on surface tension between immiscible fluids, rock wettability, aperture distribution of microcracks, compressibility of microcracks, porosity of microcracks, elastic properties of rock mineral, fluid saturation, effective stress and wave amplitude. Derived equations allow to estimate effects of pore fluids and saturation on the shear modulus and mechanical strength of rocks.     

Seismic-frequency Laboratory Measurements of Shear Mode Viscoelasticity in Crustal Rocks II: Thermally Stressed Quartzite and Granite

—Forced torsional oscillation techniques have been used to explore the seismic-frequency shear mode viscoelasticity of specimens of two crustal rocks (Cape Sorell quartzite and Delegate aplite), cycled between room temperature and 700°C under conditions of moderate confining pressure. The anisotropy and intergranular inhomogeneity of thermal expansivity in these materials give rise to large deviatoric stresses, resulting in thermal cracking at temperatures above a pressure-dependent threshold temperature, associated with the onset of very pronounced temperature sensitivity of the shear modulus, in general accord with the predictions of fracture mechanics models. For Delegate aplite in particular, the shear modulus behaves reproducibly during multiple thermal cycles at different confining pressures, consistent with the notion that the thermal cracks are of low aspect ratio (minimum/maximum dimension), and are therefore readily closed by the prevailing confining pressure once the thermal stresses are removed. Marked frequency-dependent dissipation of shear strain energy is observed on heating each rock to temperatures ≥ 500°C, although the attenuation varies significantly with prior thermal history, probably as a result of progressive dehydration and relaxation of deviatoric stresses. Temperature and pressure dependent crack densities for Delegate aplite have been estimated by comparison of the observed shear moduli with those expected for a crack-free aggregate. In parallel with the forced oscillation tests, measurements have been made of the rate at which (argon) pore pressure equilibrium is re-established following a perturbation. Combination of these results, which provide a proxy for permeability, with the inferred crack densities indicates that the variation of permeability with crack density is well described by a percolation model with a threshold crack density of ～0.2.     

地震作用下夯土城墙动力响应分析——以山丹明长城为例

2022年1月8日门源M6.9地震造成山丹明长城局部破坏。为研究此次地震作用下夯土城墙的动力响应与破坏特征,基于地震现场考察结果,采用振幅等效处理后的记录地震波为输入地震动,开展双向地震荷载作用下夯土城墙的动力响应数值分析,研究不同位置测点的最大位移、峰值加速度与墙体应力分布特征,探讨地震导致夯土城墙破坏的主要内因。研究结果表明：双向地震荷载作用下,墙体位移和峰值加速度(PGA)随着高度的增加逐渐增加,但距墙体底部0.5 m高度范围内PGA放大效应不明显,最大位移、加速度均出现在墙体顶部裂缝位置处;水平地震荷载作用下墙体的地震动响应更为显著;墙体的最大主应力、最大剪应力均出现在有裂缝处的底端掏蚀悬空部位,墙体裂缝、夯筑搭接、掏蚀悬空处应力集中明显;裂缝对夯土城墙的地震动放大效应在一定高度范围内表现为弱化作用,但随高度增加逐渐过渡为强化作用;裂缝可显著增强墙体顶部地震动响应,可能是本次地震诱发城墙破坏的主要内因。研究成果可为古城墙遗址的加固修缮提供科学指导。     

考虑劣化屈服老化裂纹后混凝土水库堤坝排水边坡抗震性研究

针对在水库堤坝排水边坡混凝土裂纹的抗震性研究中,未考虑岩土体抗剪强度参数的劣化屈服效应以及混凝土裂纹的老化,存在抗震性判断结果准确率较差等问题,提出水库堤坝排水边坡混凝土考虑老化后产生裂纹的抗震性能研究方法。模拟强震下边坡混凝土的开裂破坏过程,根据D-P屈服准则,实现对闸墩混凝土材料的屈服判断。采用薄层整体单元模拟和分离式裂纹单元,实现混凝土裂纹的数值模拟,加载地震波后,获取混凝土裂纹的强震响应规律与破坏特征。实验结果可知,本文方法对坝体位移变化的研究精度高,得到的混凝土裂纹扩展范围更为准。运用本文方法对水库堤坝排水边坡混凝土的抗震性研究准确率以及可信度较高,说明本文方法具有一定的可取性。     

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.     

Distinguishing oil and water layers by interpreting acoustic logging data with changing well diameters

During surveys, water layers may interfere with the detection of oil layers. In order to distinguish between oil and water layers, research on the properties of well diameters and oil and water layers and their relation to acoustic logging rules is essential. Using Hudson's crack theory, we simulated oil and water layers with different well diameters or crack parameters (angle and number density). We found that when the well radius increases from 0.03 m to 0.05 m, the variation ratio of compressional wave amplitude for the oil layer is less than that for the water layer. The difference of Stoneley wave amplitude between the crack parameters (angle and number density) is greater in the case of the water layer than in the case of the oil layer. The response sensitivity of wave energy is greater for the water layer than that for the oil layer. When the well radius increases from 0.05 m to 0.14 m, the maximum excitation intensity for oil layer is greater than that for the water layer. We conclude that the propagation of an elastic wave is affected by medium composition and well diameter, and the influence has certain regularity. These results can guide further reservoir logging field exploration work.     

非穿透裂纹平板试件三维破裂的实验研究

在研究地震孕育及震源破裂的问题时,用三维破裂来模拟地震破裂过程,比二维破裂与实际更加接近。我们设计了一组实验,采用几种不同的材料（玻璃、有机玻璃和大理石等）做成平板试件,在其中心预制非穿透裂纹,并在单轴压下进行了实验研究。 实验结果表明,非穿透裂纹与穿透裂纹的扩展情况有质的差别。非穿透裂纹的扩展并不沿原裂纹面延伸,而是沿裂纹前缘扩展出许多裂纹面,每一个破裂面都是复杂的扭曲、拐折的曲面。这可能意味着实际地震断层的扩展情况也是极为复杂的。     

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.     

Penny-shaped fractures revisited

All methods of seismic characterization of fractured reservoirs are based on effective media theories that relate geometrical and material properties of fractures and surrounding rock to the effective stiffnesses. In exploration seismology, the first-order theory of Hudson is the most popular. It describes the effective model caused by the presence of a single set of thin, aligned vertical fractures in otherwise isotropic rock. This model is known to be transversely isotropic with a horizontal symmetry axis (HTI). Following the theory, one can invert the effective anisotropy for the crack density and type of fluid infill of fractures, the quantities of great importance for reservoir appraisal and management.Here I compute effective media numerically using the finite element method. I deliberately construct models that contain a single set of vertical, ellipsoidal, non-intersecting and non-interconnected fractures to check validity of the first-order Hudson’s theory and establish the limits of its applicability. Contrary to conventional wisdom that Hudson’s results are accurate up to crack density e ≈ 0.1, I show that they consistently overestimate the magnitudes of all effective anisotropic coefficients ε^(V), δ^(V), and γ^(V). Accuracy of theoretically derived anisotropy depends on the type of fluid infill and typically deteriorates as e grows. While the theory gives | ε^(V)|, |δ^(V)|, |γ^(V)| and close to the upper bound of the corresponding numerically obtained values for randomly distributed liquid-filled fractures, theoretical predictions of ε^(V), δ^(V) are not supported by numerical computations when the cracks are dry. This happens primarily because the first-order Hudson’s theory makes no attempt to account for fracture interaction which contributes to the final result much stronger for gas- than for liquid-filled cracks. I find that Mori-Tanaka’s theory is superior to Hudson’s for all examined crack densities and both types of fluid infill.The paper was presented at the 11th International Workshop on Seismic Anisotropy (11IWSA) held in St. John’s, Canada in 2004.     

Experimental study on the seismic performance of new sandwich masonry walls

Sandwich masonry walls are widely used as energy-saving panels since the interlayer between the outer leaves can act as an insulation layer.New types of sandwich walls are continually being introduced in research and applications,and due to their unique bond patterns,experimental studies have been performed to investigate their mechanical properties,especially with regard to their seismic performance.In this study,three new types of sandwich masonry wall have been designed,and cyclic lateral loading tests were carried out on five specimens.The results showed that the specimens failed mainly due to slippage along the bottom cracks or the development of diagonal cracks,and the failure patterns were considerably influenced by the aspect ratio.Analysis was undertaken on the seismic response of the new walls,which included ductility,stiffness degradation and energy dissipation capacity,and no obvious difference was observed between the seismic performance of the new walls and traditional walls.Comparisons were made between the experimental results and the calculated results of the shear capacity.It is concluded that the formulas in the two Chinese codes(GB 50011 and GB 50003) are suitable for the calculation of the shear capacity for the new types of walls,and the formula in GB 50011 tends to be more conservative.     

基于二维三分量伪谱法模拟数据的EDA介质中横波分裂研究

为了检测定向裂隙介质中横波分裂的方位属性特征,分析地震属性随裂隙密度和方位变化,采用人工吸收边界和反周期扩展边界,用伪谱法获得不同裂隙密度和不同方位地质模型三分量地面记录;应用时频分析和剪切波偏振分析研究由于裂隙方位和密度引起的横波分裂.结果显示,裂隙密度和方位决定着横波分裂的时差和偏振.快慢横波的延迟时间随裂隙密度增大而增加,不同方位相同裂隙密度的横波分裂时差有微小的变化.在45°方位检测时间延迟时间最大.通过时频分析,可以看到不同方位的瞬时主频有显著的变化,在横波分裂处瞬时主频有明显变化.因此,瞬时主频和快横波的偏振以及延迟时间可以作为裂隙方位和密度的指示.     

Biot dispersion for P‐ and S‐wave velocities in partially and fully saturated sandstones

Ultrasonic compressional‐ and shear‐wave velocities have been measured on 34 samples of sandstones from hydrocarbon reservoirs. The sandstones are all of low clay content, high porosity, and cover a wide range of permeabilities. They were measured dry and brine‐saturated under hydrostatic effective stresses of 10, 20 and 40 MPa. For eight of the sandstones, ultrasonic velocity measurements were made at different partial water saturations in the range from dry to fully saturated. The Gassmann–Biot theory is found to account for most of the changes in velocities at high effective stress levels when the dry sandstones are fully saturated with brine, provided the lower velocities resulting when the dry sandstone initially adsorbs small amounts of moisture are used to determine the elastic properties of the ‘dry’ sandstone. At lower effective stress levels, local flow phenomena due to the presence of open microcracks are assumed to be responsible for measured velocities higher than those predicted by the theory. The partial saturation results are modelled fairly closely by the Gassmann–Biot theory, assuming heterogeneous saturation for P‐waves.