共查询到20条相似文献，搜索用时 487 毫秒
1.
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane. 相似文献
2.
Angular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous Media 总被引：1，自引：0，他引：1
José M. Carcione Boris Gurevich Juan E. Santos Stefano Picotti 《Pure and Applied Geophysics》2013,170(11):1673-1683
Wave-induced fluid flow generates a dominant attenuation mechanism in porous media. It consists of energy loss due to P-wave conversion to Biot (diffusive) modes at mesoscopic-scale inhomogeneities. Fractured poroelastic media show significant attenuation and velocity dispersion due to this mechanism. The theory has first been developed for the symmetry axis of the equivalent transversely isotropic (TI) medium corresponding to a poroelastic medium containing planar fractures. In this work, we consider the theory for all propagation angles by obtaining the five complex and frequency-dependent stiffnesses of the equivalent TI medium as a function of frequency. We assume that the flow direction is perpendicular to the layering plane and is independent of the loading direction. As a consequence, the behaviour of the medium can be described by a single relaxation function. We first consider the limiting case of an open (highly permeable) fracture of negligible thickness. We then compute the associated wave velocities and quality factors as a function of the propagation direction (phase and ray angles) and frequency. The location of the relaxation peak depends on the distance between fractures (the mesoscopic distance), viscosity, permeability and fractures compliances. The flow induced by wave propagation affects the quasi-shear (qS) wave with levels of attenuation similar to those of the quasi-compressional (qP) wave. On the other hand, a general fracture can be modeled as a sequence of poroelastic layers, where one of the layers is very thin. Modeling fractures of different thickness filled with CO2 embedded in a background medium saturated with a stiffer fluid also shows considerable attenuation and velocity dispersion. If the fracture and background frames are the same, the equivalent medium is isotropic, but strong wave anisotropy occurs in the case of a frameless and highly permeable fracture material, for instance a suspension of solid particles in the fluid. 相似文献
3.
Michael A. Schoenberg 《Geophysical Prospecting》2009,57(2):169-185
A vertically fractured transversely isotropic (VFTI) elastic medium is one in which any number of sets of vertical aligned fractures (each set has its normal lying in the horizontal x1, x2‐plane) pervade the medium and the sets of aligned fractures are the only features of the medium disturbing the axi‐symmetry about the x3‐axis implying that in the absence of fractures, the background medium is transversely isotropic (TI). Under the assumptions of long wavelength equivalent medium theory, the compliance matrix of a fractured medium is the sum of the background medium's compliance matrix and a fracture compliance matrix. For sets of parallel rotationally symmetric fractures (on average), the fracture compliance matrix is dependent on 3 parameters − its normal and tangential compliance and its strike direction. When one fracture set is present, the medium is orthorhombic and the analysis is straightforward. When two (non‐orthogonal) or more sets are present, the overall medium is in general elastically monoclinic; its compliance tensor components are subject to two equalities yielding an 11 parameter monoclinic medium. Constructing a monoclinic VFTI medium with n embedded vertical fracture sets, requires 5 TI parameters plus 3×n fracture set parameters. A deconstruction of such an 11 parameter monoclinic medium involves using its compliance tensor to find a background transversely isotropic medium and several sets of vertical fractures which, in the long wavelength limit, will behave exactly as the original 11 parameter monoclinic medium. A minimal deconstruction, would be to determine, from the 11 independent components, the transversely isotropic background (5 parameters) and two fracture sets (specified by 2 × 3 = 6 parameters). Two of the background TI medium's compliance matrix components are known immediately by inspection, leaving nine monoclinic components to be used in the minimal deconstruction of the VFTI medium. The use of the properties of a TI medium, which are linear relations on its compliance components, allows the deconstruction to be reduced to solving a pair of non‐linear equations on the orientations of two fracture sets. A single root yielding a physically meaningful minimum deconstruction yields a unique minimal representation of the monoclinic medium as a VFTI medium. When no such root exists, deconstruction requires an additional fracture set and uniqueness is lost. The boundary between those monoclinic media that have a unique minimal representation and those that do not is yet to be determined. 相似文献
4.
Michael Schoenberg 《Geophysical Prospecting》1994,42(8):885-915
Any set of isotropic layers is equivalent, in the long wavelength limit, to a unique transversely isotropic (TI) layer; to find the elastic moduli of that layer is a solved problem. The converse problem is to find a set of isotropic layers equivalent to a given TI media. Here, explicit necessary and sufficient conditions on the TI stiffness moduli for the existence of an equivalent set of isotropic layers are found by construction of a minimal decomposition consisting of either two or three isotropic constituent layers. When only two constituents are required, their elastic properties are uniquely determined. When three constituents are required, two have the same Poisson's ratio and the same thickness fraction, and even then there is a one-parameter family of satisfactory minimal decompositions. The linear slip model for fractured rock (aligned fractures in an isotropic background) yields a restricted range of transverse isotropy dependent on only four independent parameters. If the ratio of the normal to tangential fracture compliance is small enough, the medium is equivalent to thin isotropic layering and in general its minimal decomposition consists of three constituents. 相似文献
5.
A model of parallel slip interfaces simulates the behaviour of a fracture system composed of large, closely spaced, aligned joints. The model admits any fracture system anisotropy: triclinic (the most general), monoclinic, orthorhombic or transversely isotropic, and this is specified by the form of the 3 × 3 fracture system compliance matrix. The fracture system may be embedded in an anisotropic elastic background with no restrictions on the type of anisotropy. To compute the long wavelength equivalent moduli of the fractured medium requires at most the inversion of two 3 × 3 matrices. When the fractures are assumed on average to have rotational symmetry (transversely isotropic fracture system behaviour) and the background is assumed isotropic, the resulting equivalent medium is transversely isotropic and the effect of the additional compliance of the fracture system may be specified by two parameters (in addition to the two isotropic parameters of the isotropic background). Dilute systems of flat aligned microcracks in an isotropic background yield an equivalent medium of the same form as that of the isotropic medium with large joints, i.e. there are two additional parameters due to the presence of the microcracks which play roles in the stress-strain relations of the equivalent medium identical to those played by the parameters due to the presence of large joints. Thus, knowledge of the total of four parameters describing the anisotropy of such a fractured medium tells nothing of the size or concentration of the aligned fractures but does contain information as to the overall excess compliance due to the fracture system and its orientation. As the aligned microcracks, which were assumed to be ellipsoidal, with very small aspect ratio are allowed to become non-fiat, i.e. have a growing aspect ratio, the moduli of the equivalent medium begin to diverge from the standard form of the moduli for flat cracks. The divergence is faster for higher crack densities but only becomes significant for microcracks of aspect ratios approaching 0.3. 相似文献
6.
裂缝诱导HTI双孔隙介质模型是将一组垂直排列的裂缝系统嵌入到统计各向同性的孔隙岩石基质系统中而建立的.为了研究裂缝参数对地震波在该模型中传播规律的影响,本文分别对裂缝弱度、裂缝孔隙度和裂缝渗透率这三个主要的裂缝参数进行了分析研究.数值结果表明,裂缝诱导HTI双孔隙介质中,裂缝弱度越大,介质的各向异性强度越强;与基质孔隙系统相比,裂缝系统孔隙度对介质等效孔隙度的影响很小,而裂缝系统渗透率的增大则将显著提高介质在裂缝发育方向上的等效渗透率,这符合对裂缝系统"低孔"、"高渗"特性的认识.此外,裂缝系统渗透率的增大也使慢纵波的振幅显著增强. 相似文献
7.
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. 相似文献
8.
Colin M. Sayers 《Geophysical Prospecting》2002,50(4):393-401
Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. A shear wave propagating in an azimuthally anisotropic medium splits into two components with different polarizations if the source polarization is not aligned with the principal axes of the medium. For vertical propagation of shear waves in a horizontally layered medium containing vertical fractures, the shear‐wave splitting depends on the shear compliance of the fractures, but is independent of their normal compliance. If the fractures are not perfectly vertical, the shear‐wave splitting also depends on the normal compliance of the fractures. The normal compliance of a fluid‐filled fracture decreases with increasing fluid bulk modulus. For dipping fractures, this results in a decrease in shear‐wave splitting and an increase in shear‐wave velocity with increasing fluid bulk modulus. The sensitivity of the shear‐wave splitting to fluid bulk modulus depends on the interconnectivity of the fracture network, the permeability of the background medium and on whether the fracture is fully or partially saturated. 相似文献
9.
Although clay is composed of disconnected anisotropic clay platelets, many rock physics models treat the clay platelets in shale as interconnected. However, the clay matrix in shales can be modelled as anisotropic clay platelets embedded within a soft isotropic interplatelet region, allowing the influence of disconnected clay platelets on the elastic properties of the clay matrix to be analysed. In this model, properties of the interplatelet region are governed by its effective bulk and shear moduli, whereas the effective properties of the clay platelets are governed by their volume fraction, aspect ratio and elastic stiffness tensor. Together, these parameters implicitly account for variations in clay and fluid properties, as well as fluid saturation. Elastic stiffnesses of clay platelets are obtained from the literature, including both experimental measurements and first-principles calculations of the full anisotropic (monoclinic or triclinic) elastic stiffness tensors of layered silicates. These published elastic stiffness tensors are used to compile a database of equivalent transverse isotropic elastic stiffness tensors, and other physical properties, for eight common varieties of layered silicates. Clay matrix anisotropy is then investigated by examining the influence of these different elastic stiffnesses, and of varying model parameters, upon the effective transverse isotropic elastic stiffness tensor of the clay matrix. The relationship between the different clay minerals and their associated anisotropy parameters is studied, and their impact on the resulting anisotropy of the clay matrix is analysed. 相似文献
10.
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. 相似文献
11.
Understanding the effects of in situ fluid content and fracture parameters on seismic characteristics is important for the subsurface exploration and production of fractured porous rocks. The ratio of normal-to-shear fracture compliance is typically utilized as a fluid indicator to evaluate anisotropy and identify fluids filling the fractures, but it represents an underdetermined problem because this fluid indicator varies as a function of both fracture geometry and fluid content. On the bases of anisotropic Gassmann's equation and linear-slip model, we suggest an anisotropic poroelasticity model for fractured porous reservoirs. By combining a perturbed stiffness matrix and asymptotic ray theory, we then construct a direct relationship between the PP-wave reflection coefficients and characteristic parameters of fluids(P-and S-wave moduli) and fractures(fracture quasi-weaknesses), thereby decoupling the effects of fluid and fracture properties on seismic reflection characterization.By incorporating fracture quasi-weakness parameters, we propose a novel parameterization method for elastic impedance variation with offset and azimuth(EIVOA). By incorporating wide-azimuth observable seismic reflection data with regularization constraints, we utilize Bayesian seismic inversion to estimate the fluid content and fracture parameters of fractured porous rocks. Tests on synthetic and real data demonstrate that fluid and fracture properties can be reasonably estimated directly from azimuthal seismic data and the proposed approach provides a reliable method for fluid identification and fracture characterization in a gas-saturated fractured porous reservoir. 相似文献
12.
结合有限差分方法和等效介质理论，模拟了离散分布裂缝介质中地震波的传播. 基于等效介质理论，利用二维有限差分实现封闭裂缝的离散分布；裂缝可以处理成固体岩石中的高度柔性界面，并可以用线性滑动或者位移间断模型进行裂缝的物理描述. 对于含有多组裂隙的破裂固体，其有效柔度可以认为是固体骨架背景柔度和裂缝附加柔度之和. 在一阶近似条件下，固体骨架和裂缝参数可以通过有效各向异性系数联系起来，有效各向异性系数决定了各向异性（裂缝效应）对于地震波传播的影响. 通过与射线理论方法的对比检验，说明本文提出的模拟方法的有效性，并通过几个数值算例说明本方法可有效模拟不同的裂缝分布效应. 结果表明，即使在裂缝密度很小的情况下，具有相同裂缝密度的不同的空间分布可以产生不同的波场特征. 同时，也验证了不同裂缝尺度对波长的不同影响，以及裂缝尺度具有幂率分布（分形）时，尺度对波场的影响. 最后得出结论：在运用建立在等效介质理论基础上的地震各向异性概念来描述裂缝固体的特征时，要倍加小心，等效介质理论中尚未合理处理的裂缝尺度和空间分布对波的传播特征具有重要的影响. 相似文献
13.
The effect of anisotropy on the onset of convection in a saturated porous medium is discussed. In particular, the case of time-dependent density-driven convection is examined. The applicability of the value of an equivalent Rayleigh number as the criterion for the onset of convection is discussed. 相似文献
14.
Quantifying Damage, Saturation and Anisotropy in Cracked Rocks by Inverting Elastic Wave Velocities 总被引：1，自引：0，他引：1
Alexandre Schubnel Philip M. Benson Ben D. Thompson Jim F. Hazzard R. Paul Young 《Pure and Applied Geophysics》2006,163(5-6):947-973
Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive
crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms
of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic
properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived.
Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching
up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range
of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering
symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in
order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and
give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms
of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More
importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture)
with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect
ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic
and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and
mean crack fabric orientations we recovered are consistent with in situ geophysical imaging.
Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor
and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate
the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using
such an inversion tool opens the door in linking elastic properties, variations to permeability. 相似文献
15.
Computing effective medium properties is very important when upscaling data measured at small scale. In the presence of stratigraphic layering, seismic velocities and anisotropy parameters are scale and frequency dependent. For a porous layer permeated by aligned fractures, wave-induced fluid flow between pores and fractures can also cause significant dispersion in velocities and anisotropy parameters. In this study, we compare the dispersion of anisotropy parameters due to fracturing and layering at low frequencies. We consider a two-layer model consisting of an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce anisotropy parameters dispersion due to fractures (meso-scale) in the sand layer. This intrinsic dispersion is added to anisotropy parameters dispersion induced by layering (macro-scale) at low frequencies. We derive the series coefficients that control the behaviour of anisotropy parameters at low frequencies. We investigate the influences of fracture length and fracture density on fracturing effect, layering effect and combined effect versus frequency and volume fraction of sand layer. Numerical modelling results indicate that the frequency dependence due to layering is not always the dominant effect of the effective properties of the medium. The intrinsic dispersion is not negligible compared with the layering effect while evaluating the frequency-dependent properties of the layered medium. 相似文献
16.
假设横向各向同性(TI)介质的对称轴是垂直的(VTI)或者水平的(HTI)能给实际资料处理带来便利,然而实际TI介质的对称轴往往是倾斜的(TTI),忽略对称轴倾角可能给各向异性参数提取和成像带来偏差,因此需要研究是否能、以及什么条件下能忽略TTI介质对称轴倾角.本文通过理论研究和数值分析研究了与TTI介质弹性性质最接近的VTI介质(OAVTI)的弹性常数和各向异性参数与原TTI介质的弹性常数和各向异性参数之间的联系与差别.结果表明:OAVTI介质各向异性参数与原TTI介质各向异性参数之间的差别可统一表示成F(α0/β0,ε,δ,γ)ξ2的形式,其中F(α0/β0,ε,δ,γ)是无量纲各向异性参数(ε, δ, γ)的线性函数,ξ是对称轴倾角;ξ的大小对各参数的误差起主导作用,一般不建议忽略20°~25°以上的对称轴倾角;当ξ较小时,即使是对强各向异性的TTI介质作VTI近似,引起的P波各向异性参数误差也很小,因此在纵波资料处理中忽略TTI介质对称轴倾角通常是可行的;即使在小ξ条件下,倾斜对称轴对SV波也有显著影响,因此在转换波资料处理中,不建议忽略TTI介质的对称轴倾角.本文的研究为分析忽略TTI介质对称轴倾角的可行性提供了理论依据和简便的判据. 相似文献
17.
Fractures in elastic media add compliance to a rock in the direction normal to the fracture strike. Therefore, elastic wave velocities in a fractured rock will vary as a function of the energy propagation direction relative to the orientation of the aligned fracture set. Anisotropic Thomson–Haskell matrix Rayleigh-wave equations for a vertically transverse isotropic media can be used to model surface-wave dispersion along the principal axes of a vertically fractured and transversely isotropic medium. Furthermore, a workflow combining first-break analysis and azimuthal anisotropic Rayleigh-wave inversion can be used to estimate P-wave and S-wave velocities, Thomsen's ε, and Thomsen's δ along the principal axes of the orthorhombic symmetry. In this work, linear slip theory is used to map our inversion results to the equivalent vertically fractured and transversely isotropic medium coefficients. We carried out this inversion on a synthetic example and a field example. The synthetic data example results show that joint estimation of S-wave velocities with Thomsen's parameters ε and δ along normal and parallel to the vertical fracture set is reliable and, when mapped to the corresponding vertically fractured and transversely isotropic medium, provides insight into the fracture compliances. When the inversion was carried out on the field data, results indicated that the fractured rock is more compliant in the azimuth normal to the visible fracture set orientation and that the in situ normal fracture compliance to tangential fracture compliance ratio is less than half, which implies some cementation may have occurred along the fractures. Such an observation has significant implications when modelling the transport properties of the rock and its strength. Both synthetic and field examples show the potential of azimuthal anisotropic Rayleigh-wave inversion as the method can be further expanded to a more general case where the vertical fracture set orientation is not known a priori. 相似文献
18.
Transmitting boundary for water-saturated transversely isotropic strata based on the u–U formulation
A new transmitting boundary in a cylindrical coordinate system has been developed for modeling the elastic waves radiating out to an infinite boundary in water-saturated transversely isotropic soil strata over a rigid bedrock. The saturated soil strata are assumed to consist of a porous material and modeled as a transversely isotropic two-phase medium, based on the u−U formulation. The newly developed transmitting boundary is combined with the finite-elements model of the near-field region, using the same u−U formulation, and applied to the study of the dynamics of a rigid circular foundation in porous isotropic or transversely isotropic layered strata, either fully or partly saturated with water. The verification and application examples give valuable insights into new and interesting aspects of the dynamic behavior of rigid circular foundations in fully or partly saturated two-phase ground in terms of permeability, transverse anisotropy, and ground-water table level. 相似文献
19.
Effects of fracture intersections on seismic dispersion: theoretical predictions versus numerical simulations 
下载免费PDF全文
下载免费PDF全文 Junxin Guo J. Germán Rubino Stanislav Glubokovskikh Boris Gurevich 《Geophysical Prospecting》2017,65(5):1264-1276
The detection and characterisation of domains of intersecting fractures are important goals in several disciplines of current interest, including exploration and production of unconventional reservoirs, nuclear waste storage, CO2 sequestration, and groundwater hydrology, among others. The objective of this study is to propose a theoretical framework for quantifying the effects of fracture intersections on the frequency‐dependent elastic properties of fluid‐saturated porous and fractured rocks. Three characteristic frequency regimes for fluid pressure communication are identified. In the low‐frequency limit, fractures are in full pressure communication with the embedding porous matrix and with other fractures. Conversely, in the high‐frequency limit, fractures are hydraulically isolated from the matrix and from other fractures. At intermediate frequencies, fractures are hydraulically isolated from the matrix porosity but can be in hydraulic communication with each other, depending on whether fracture sets are intersecting. For each frequency regime, the effective stiffness coefficients are derived using the linear‐slip theory and anisotropic Gassmann equations. Explicit mathematical expressions for the two characteristic frequencies that separate the three frequency regimes are also determined. Theoretical predictions are then applied to two synthetic 2D samples, each containing two orthogonal fracture sets: one with and another without intersections. The resulting stiffness coefficients, Thomsen‐style anisotropy parameters, and the transition frequencies show good agreement with corresponding numerical simulations. The theoretical results are applicable not only to 2D but also to 3D fracture systems and are amenable to being employed in inversion schemes designed to characterise fracture systems. 相似文献
20.
We have implemented a 3D finite‐difference scheme to simulate wave propagation in arbitrary anisotropic media. The anisotropic media up to orthorhombic symmetry were modelled using a standard staggered grid scheme and beyond (monoclinic and triclinic) using a rotated staggered grid scheme. The rationale of not using rotated staggered grid for all types of anisotropic media is that the rotated staggered grid schemes are more expensive than standard staggered grid schemes. For a 1D azimuthally anistropic medium, we show a comparison between the seismic data generated by our finite‐difference code and by the reflectivity algorithm; they are in excellent agreement. We conducted a study on zero‐offset shear‐wave splitting using the finite‐difference modelling algorithm using the rotated staggered grid scheme. Our S‐wave splitting study is mainly focused on fractured media. On the scale of seismic wavelenghts, small aligned fractures behave as an equivalent anisotropic medium. We computed the equivalent elastic properties of the fractures and the background in which the fractures were embedded, using low‐frequency equivalent media theories. Wave propagation was simulated for both rotationally invariant and corrugated fractures embedded in an isotropic background for one, or more than one, set of fluid‐filled and dry fractures. S‐wave splitting was studied for dipping fractures, two vertical non‐orthogonal fractures and corrugated fractures. Our modelling results confirm that S‐wave splitting can reveal the fracture infill in the case of dipping fractures. S‐wave splitting has the potential to reveal the angle between the two vertical fractures. We also notice that in the case of vertical corrugated fractures, S‐wave splitting is sensitive to the fracture infill. 相似文献