共查询到20条相似文献，搜索用时 609 毫秒
1.
Fluid flow in many hydrocarbon reservoirs is controlled by aligned fractures which make the medium anisotropic on the scale of seismic wavelength. Applying the linear‐slip theory, we investigate seismic signatures of the effective medium produced by a single set of ‘general’ vertical fractures embedded in a purely isotropic host rock. The generality of our fracture model means the allowance for coupling between the normal (to the fracture plane) stress and the tangential jump in displacement (and vice versa). Despite its low (triclinic) symmetry, the medium is described by just nine independent effective parameters and possesses several distinct features which help to identify the physical model and estimate the fracture compliances and background velocities. For example, the polarization vector of the vertically propagating fast shear wave S1 and the semi‐major axis of the S1‐wave normal‐moveout (NMO) ellipse from a horizontal reflector always point in the direction of the fracture strike. Moreover, for the S1‐wave both the vertical velocity and the NMO velocity along the fractures are equal to the shear‐wave velocity in the host rock. Analysis of seismic signatures in the limit of small fracture weaknesses allows us to select the input data needed for unambiguous fracture characterization. The fracture and background parameters can be estimated using the NMO ellipses from horizontal reflectors and vertical velocities of P‐waves and two split S‐waves, combined with a portion of the P‐wave slowness surface reconstructed from multi‐azimuth walkaway vertical seismic profiling (VSP) data. The stability of the parameter‐estimation procedure is verified by performing non‐linear inversion based on the exact equations. 相似文献
2.
Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities. 相似文献
3.
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. 相似文献
4.
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. 相似文献
5.
Weak‐anisotropy approximation for P‐wave reflection coefficient at the boundary between two tilted transversely isotropic media 
下载免费PDF全文
下载免费PDF全文 Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy. 相似文献
6.
长波长假设条件下,各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性（HTI）介质.基于不同观测方位的岩石地震响应特征变化,宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测,同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法,首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数,并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联,并采用马尔科夫链蒙特卡洛（MCMC）抽样方法进行大量样本的随机模拟,使用期望最大化（EM）算法估计后验条件概率分布,最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明,该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数,并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法. 相似文献
7.
By analogy with P- and S-wave impedances, the product of Young's modulus and density can be termed as Young's impedance, which indicates the rock lithology and brittleness of unconventional hydrocarbon reservoirs. Poisson's ratio is also an effective indicator of rock brittleness and fluid property of unconventional reservoirs, and fracture weaknesses indicate the fracture properties (fracturing intensity and fracture fillings) in fracture-induced unconventional reservoirs. We aim to simultaneously estimate the Young's impedance, Poisson's ratio and fracture weaknesses from wide-azimuth surface seismic data in a fracture-induced shale gas reservoir, and use the horizontal transversely isotropic model to characterize the fractures. First, the linearized PP-wave reflection coefficient in terms of Young's impedance, Poisson's ratio, density and fracture weaknesses is derived for the case of a weak-contrast interface separating two weakly horizontal transversely isotropic media. In addition, an orthorhombic anisotropic case is also discussed in this paper. Then a Bayesian amplitude variation with incident angle and azimuth scheme with a model constraint is used to stably estimate Young's impedance, Poisson's ratio and fracture weaknesses with only PP-wave azimuthal seismic data. The proposed approach is finally demonstrated on both synthetic and real data sets with reasonable results. 相似文献
8.
Anisotropy is often observed due to the thin layering or aligned micro‐structures, like small fractures. At the scale of cross‐well tomography, the anisotropic effects cannot be neglected. In this paper, we propose a method of full‐wave inversion for transversely isotropic media and we test its robustness against structured noisy data. Optimization inversion techniques based on a least‐square formalism are used. In this framework, analytical expressions of the misfit function gradient, based on the adjoint technique in the time domain, allow one to solve the inverse problem with a high number of parameters and for a completely heterogeneous medium. The wave propagation equation for transversely isotropic media with vertical symmetry axis is solved using the finite difference method on the cylindrical system of coordinates. This system allows one to model the 3D propagation in a 2D medium with a revolution symmetry. In case of approximately horizontal layering, this approximation is sufficient. The full‐wave inversion method is applied to a crosswell synthetic 2‐component (radial and vertical) dataset generated using a 2D model with three different anisotropic regions. Complex noise has been added to these synthetic observed data. This noise is Gaussian and has the same amplitude f?k spectrum as the data. Part of the noise is localized as a coda of arrivals, the other part is not localized. Five parameter fields are estimated, (vertical) P‐wave velocity, (vertical) S‐wave velocity, volumetric mass and the Thomsen anisotropic parameters epsilon and delta. Horizontal exponential correlations have been used. The results show that the full‐wave inversion of cross‐well data is relatively robust for high‐level noise even for second‐order parameters such as Thomsen epsilon and delta anisotropic parameters. 相似文献
9.
Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation. 相似文献
10.
11.
The moveout approximations play an important role in seismic data processing. The standard hyperbolic moveout approximation is based on an elliptical background model with two velocities: vertical and normal moveout. We propose a new set of moveout approximations based on a perturbation series in terms of anellipticity parameters using the alternative elliptical background model defined by vertical and horizontal velocities. We start with a transversely isotropic medium with a vertical symmetry axis. Then, we extend this approach to a homogeneous orthorhombic medium. To define the perturbation coefficients for a new background, we solve the eikonal equation with horizontal velocities in transversely isotropic medium with a vertical symmetry axis and orthorhombic media. To stabilise the perturbation series and improve the accuracy, the Shanks transform is applied for all the cases. We select different parameterisations for both velocities and anellipticity parameters for an orthorhombic model. From the comparison in traveltime error, the new moveout approximations result in better accuracy comparing with the standard perturbation‐based methods and other approximations. 相似文献
12.
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. The symmetry of a rock with horizontal bedding that contains two or more non-orthogonal sets of vertical fractures may be approximated as monoclinic with a horizontal plane of mirror symmetry. For offsets that are small compared with the depth of the reflector, the azimuthal variation in P-wave AVO gradient for such a medium varies with azimuth as where φ is the azimuth measured with respect to the fast polarization direction for a vertically polarized shear wave. φ 2 depends on both the normal compliance B N and the shear compliance B T of the fractures and may differ from zero if B N B T varies significantly between fracture sets. If B N B T is the same for all fractures, and the principal axes of the azimuthal variation in P-wave AVO for fixed offset are determined by the polarization directions of a vertically propagating shear wave. At larger offsets, terms in and are required to describe the azimuthal variation in AVO accurately. φ 4 and φ 6 also depend on B N B T . For gas-filled open fractures but a lower value of B N B T may result from the presence of a fluid with non-zero bulk modulus. 相似文献
13.
In seismic data processing, serious problems could be caused by the existence of triplication and need to be treated properly for tomography and other inversion methods. The triplication in transversely isotropic medium with a vertical symmetry axis has been well studied and concluded that the triplicated traveltime only occurs for S wave and there is no triplication for P and converted PS waves since the P wave convexity slowness always compensates the S wave slowness concavity. Compared with the vertical symmetry axis model, the research of the triplication in transversely isotropic medium with a tilted symmetry axis is still keeping blank. In order to analyse the triplication for the converted wave in the tilted symmetry axis model, we examine the traveltime of the triplication from the curvature of averaged P and S wave slowness. Three models are defined and tested in the numerical examples to illustrate the behaviour of the tilted symmetry axis model for the triplicated traveltime with the change of the rotation angle. Since the orientation of an interface is related to the orientation of the symmetry axis, the triplicated traveltime is encountered for the converted wave in the tilted symmetry axis model assuming interfaces to be planar and horizontal. The triplicated region is influenced by the place and level of the concave curvature of the P and S wave slowness. 相似文献
14.
Tariq Alkhalifah 《Geophysical Prospecting》2016,64(2):505-513
While velocity contrasts are responsible for most of the events recorded in our data, the long wavelength behavior of the velocity model is responsible for the geometrical shape of these events. For isotropic acoustic materials, the wave dependency on the long (wave propagation) and short (scattering) wavelength velocity components is stationary with the propagation angle. On the other hand, in representing a transversely isotropic with a vertical symmetry axis medium with the normal moveout velocity, the anellepticity parameter η, the vertical scaling parameter δ, and the sensitivity of waves vary with the polar angle for both the long and short wavelength features of the anisotropic dimensionless medium parameters (δ and η). For horizontal reflectors at reasonable depths, the long wavelength features of the η model is reasonably constrained by the long offsets, whereas the short wavelength features produce very week reflections at even reasonable offsets. Thus, for surface acquired seismic data, we could mainly invert for smooth η responsible for the geometrical shape of reflections. On the other hand, while the δ long wavelength components mildly affects the recorded data, its short wavelength variations can produce reflections at even zero offset, with a behavior pattern synonymous to density. The lack of the long wavelength δ information will mildly effect focusing but will cause misplacement of events in depth. With low enough frequencies (very low), we may be able to recover the long wavelength δ using full waveform inversion. However, unlike velocity, the frequencies needed for that should be ultra‐low to produce long‐wavelength scattering‐based model information as δ perturbations do not exert scattering at large offsets. For a combination given by the horizontal velocity, η, and ε, the diving wave influence of η is absorbed by the horizontal velocity, severely limiting the η influence on the data and full waveform inversion. As a result, with a good smooth η estimation, for example, from tomography, we can focus the full waveform inversion to invert for only the horizontal velocity and maybe ε as a parameter to fit the amplitude. This is possibly the most practical parametrization for inversion of surface seismic data in transversely isotropic with vertical symmetry axis media. 相似文献
15.
Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media 下载免费PDF全文
下载免费PDF全文 A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account. 相似文献
16.
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. 相似文献
17.
Seismic wavefield modeling in media with fluid-filled fractures and surface topography 总被引：3，自引：0，他引：3
We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures. 相似文献
18.
O.H. Al‐Harrasi A. Al‐Anboori A. Wüstefeld J.‐M. Kendall 《Geophysical Prospecting》2011,59(2):227-243
The study of seismic anisotropy in exploration seismology is gaining interest as it provides valuable information about reservoir properties and stress directions. In this study we estimate anisotropy in a petroleum field in Oman using observations of shear‐wave splitting from microseismic data. The data set was recorded by arrays of borehole geophones deployed in five wells. We analyse nearly 3400 microearthquakes, yielding around 8500 shear‐wave splitting measurements. Stringent quality control reduces the number of reliable measurements to 325. Shear‐wave splitting modelling in a range of rock models is then used to guide the interpretation. The difference between the fast and slow shear‐wave velocities along the raypath in the field ranges between 0–10% and it is controlled both by lithology and proximity to the NE‐SW trending graben fault system that cuts the field formations. The anisotropy is interpreted in terms of aligned fractures or cracks superimposed on an intrinsic vertical transversely isotropic (VTI) rock fabric. The highest magnitudes of anisotropy are within the highly fractured uppermost unit of the Natih carbonate reservoir. Anisotropy decreases with depth, with the lowest magnitudes found in the deep part of the Natih carbonate formation. Moderate amounts of anisotropy are found in the shale cap rock. Anisotropy also varies laterally with the highest anisotropy occurring either side of the south‐eastern graben fault. The predominant fracture strikes, inferred from the fast shear‐wave polarizations, are consistent with the trends of the main faults (NE‐SW and NW‐SE). The majority of observations indicate subvertical fracture dip (>70° ). Cumulatively, these observations show how studies of shear‐wave splitting using microseismic data can be used to characterize fractures, important information for the exploitation of many reservoirs. 相似文献
19.
V. Grechka 《Studia Geophysica et Geodaetica》2005,49(3):365-381
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. 相似文献
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. 相似文献