共查询到20条相似文献，搜索用时 265 毫秒
1.
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. 相似文献
2.
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. 相似文献
3.
把遗传算法引入到了TTI介质AVO信息反演各向异性参数的过程中,依据TTI介质PP波、PS波反射系数公式,建立Thomsen参数和TTI介质对称轴倾角、方位角的目标函数,分别通过PP波和PS波的反射系数反演出了各向异性参数和对称轴倾角、方位角等信息.文中对反演结果的精确度和稳定性进行了分析,发现PS波的反演结果优于PP波反演结果;对称轴倾角的反演准确性明显优于对称轴方位角.本文通过模型正演合理解释了这一现象的原因.最后,本文通过对PS波AVO梯度的研究,提出了利用PS波振幅定性分析TTI介质对称轴倾角的方法. 相似文献
4.
5.
When a porous layer is permeated by mesoscale fractures, wave-induced fluid flow between pores and fractures can cause significant attenuation and dispersion of velocities and anisotropy parameters in the seismic frequency band. This intrinsic dispersion due to fracturing can create frequency-dependent reflection coefficients in the layered medium. In this study, we derive the frequency-dependent PP and PS reflection coefficients versus incidence angle in the fractured medium. We consider a two-layer vertical transverse isotropy model constituted by an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce the intrinsic dispersion due to fracturing in the sand layer. Based on the series coefficients that control the behaviour of velocity and anisotropy parameters in the fractured medium at low frequencies, we extend the conventional amplitude-versus-offset equations into frequency domain and derive frequency-dependent amplitude-versus-offset equations at the elastic–anelastic surface. Increase in fracture length or fracture density can enlarge the frequency dependence of amplitude-versus-offset attributes of PP and PS waves. Also, the frequency dependence of magnitude and phase angle of PP and PS reflection coefficients increases as fracture length or fracture density increases. Amplitude-versus-offset type of PP and PS reflection varies with fracture parameters and frequency. What is more, fracture length shows little impact on the frequency-dependent critical phase angle, while the frequency dependence of the critical phase angle increases with fracture density. 相似文献
6.
Joint PP‐ and PSV‐wave amplitudes versus offset and azimuth inversion for fracture compliances in horizontal transversely isotropic media 
下载免费PDF全文
下载免费PDF全文 We propose a robust approach for the joint inversion of PP‐ and PSV‐wave angle gathers along different azimuths for the elastic properties of the homogeneous isotropic host rock and excess compliances due to the presence of fractures. Motivated by the expression of fluid content indicator in fractured reservoirs and the sensitivity of Lamé impedances to fluid type, we derive PP‐ and PSV‐wave reflection coefficients in terms of Lamé impedances, density, and fracture compliances for an interface separating two horizontal transversely isotropic media. Following a Bayesian framework, we construct an objective function that includes initial models. We employ the iteratively reweighted least‐squares algorithm to solve the inversion problem to estimate unknown parameters (i.e., Lamé impedances, density, and fracture compliances) from PP‐ and PSV‐wave angle gathers along different azimuths. Synthetic tests reveal that the unknown parameters estimated using the joint inversion approach match true values better than those estimated using a PP‐wave amplitude inversion only. A real data test indicates that reasonable results for subsurface fracture detection are obtained from the joint inversion approach. 相似文献
7.
Based on knowledge of a commutative group calculation of the rock stiffness and on some geophysical assumptions, the simplest fractured medium may be regarded as a fracture embedded in an isotropic background medium, and the fracture interface can be simulated as a linear slip interface that satisfies non‐welded contact boundary conditions: the kinematic displacements are discontinuous across the interface, whereas the dynamic stresses are continuous across the interface. The finite‐difference method with boundary conditions explicitly imposed is advantageous for modelling wave propagation in fractured discontinuous media that are described by the elastic equation of motion and non‐welded contact boundary conditions. In this paper, finite‐difference schemes for horizontally, vertically, and orthogonally fractured media are derived when the fracture interfaces are aligned with the boundaries of the finite‐difference grid. The new finite‐difference schemes explicitly have an additional part that is different from the conventional second‐order finite‐difference scheme and that directly describes the contributions of the fracture to the wave equation of motion in the fractured medium. The numerical seismograms presented, to first order, show that the new finite‐difference scheme is accurate and stable and agrees well with the results of previously published finite‐difference schemes (the Coates and Schoenberg method). The results of the new finite‐difference schemes show how the amplitude of the reflection produced by the fracture varies with the fracture compliances. Later, comparisons with the reflection coefficients indicate that the reflection coefficients of the fracture are frequency dependent, whereas the reflection coefficients of the impedance contrast interface are frequency independent. In addition, the numerical seismograms show that the reflections of the fractured medium are equal to the reflections of the background medium plus the reflections of the fracture in the elastic fractured medium. 相似文献
8.
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. 相似文献
9.
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. 相似文献
10.
长波长假设条件下,各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性（HTI）介质.基于不同观测方位的岩石地震响应特征变化,宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测,同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法,首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数,并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联,并采用马尔科夫链蒙特卡洛（MCMC）抽样方法进行大量样本的随机模拟,使用期望最大化（EM）算法估计后验条件概率分布,最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明,该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数,并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法. 相似文献
11.
Seismic imaging for mineral exploration in hardrock environment has typically relied on P-waves recorded using vertical component geophones. Here we demonstrate the potential of using P-to-S converted waves for mineral exploration by applying converted wave processing to high-resolution 3-component data acquired in the Flin Flon mining camp, Canada. Inferences regarding the relative value of PS data are made through comparison with corresponding PP images, interpretation with geological constraints, and by comparison with modeling results. Modeling results show that the horizontal component data provide a robust means of separating the PS conversions from the PP reflections in converted wave processing whereas contamination by PP energy is more severe if the vertical component data are used. A customized processing workflow has been developed in which horizontal component data are sorted into common-reflection point gathers by the application of the converted wave dip-moveout operator. The PS images are characterized by lower signal-to-noise ratios than the PP images due to the smaller S-wave impedance contrasts associated with the main lithologic units and the angular dependence of P-to-S conversion. An exception to this occurs for PS reflections from the ore zones which are accentuated on the PS images relative to the background host-rock reflectivity. There are discrepancies in the reflection dips and positions observed on corresponding PP and PS sections that are not fully explained. They can be partly attributed to the non-coincidence of subsurface spatial sampling associated with common-depth points versus common-conversion points, with an additional complication associated with the known 3D geology. Vp/Vs ratios determined by correlation of reflections observed on PP and PS images can provide a limited means of lithological discrimination. 相似文献
12.
Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t2(x2) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (VP0 and VS0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (VP0, VS0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves. 相似文献
13.
Seismic amplitude variation with offset and azimuth (AVOaz) inversion is well known as a popular and pragmatic tool utilized to estimate fracture parameters. A single set of vertical fractures aligned along a preferred horizontal direction embedded in a horizontally layered medium can be considered as an effective long-wavelength orthorhombic medium. Estimation of Thomsen’s weak-anisotropy (WA) parameters and fracture weaknesses plays an important role in characterizing the orthorhombic anisotropy in a weakly anisotropic medium. Our goal is to demonstrate an orthorhombic anisotropic AVOaz inversion approach to describe the orthorhombic anisotropy utilizing the observable wide-azimuth seismic reflection data in a fractured reservoir with the assumption of orthorhombic symmetry. Combining Thomsen’s WA theory and linear-slip model, we first derive a perturbation in stiffness matrix of a weakly anisotropic medium with orthorhombic symmetry under the assumption of small WA parameters and fracture weaknesses. Using the perturbation matrix and scattering function, we then derive an expression for linearized PP-wave reflection coefficient in terms of P- and S-wave moduli, density, Thomsen’s WA parameters, and fracture weaknesses in such an orthorhombic medium, which avoids the complicated nonlinear relationship between the orthorhombic anisotropy and azimuthal seismic reflection data. Incorporating azimuthal seismic data and Bayesian inversion theory, the maximum a posteriori solutions of Thomsen’s WA parameters and fracture weaknesses in a weakly anisotropic medium with orthorhombic symmetry are reasonably estimated with the constraints of Cauchy a priori probability distribution and smooth initial models of model parameters to enhance the inversion resolution and the nonlinear iteratively reweighted least squares strategy. The synthetic examples containing a moderate noise demonstrate the feasibility of the derived orthorhombic anisotropic AVOaz inversion method, and the real data illustrate the inversion stabilities of orthorhombic anisotropy in a fractured reservoir. 相似文献
14.
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. 相似文献
15.
Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing. 相似文献
16.
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. 相似文献
17.
Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies. 相似文献
18.
19.
Ilya Tsvankin 《Geophysical Prospecting》1997,45(3):479-512
The transversely isotropic (TI) model with a tilted axis of symmetry may be typical, for instance, for sediments near the flanks of salt domes. This work is devoted to an analysis of reflection moveout from horizontal and dipping reflectors in the symmetry plane of TI media that contains the symmetry axis. While for vertical and horizontal transverse isotropy zero-offset reflections exist for the full range of dips up to 90°, this is no longer the case for intermediate axis orientations. For typical homogeneous models with a symmetry axis tilted towards the reflector, wavefront distortions make it impossible to generate specular zero-offset reflected rays from steep interfaces. The ‘missing’ dipping planes can be imaged only in vertically inhomogeneous media by using turning waves. These unusual phenomena may have serious implications in salt imaging. In non-elliptical TI media, the tilt of the symmetry axis may have a drastic influence on normal-moveout (NMO) velocity from horizontal reflectors, as well as on the dependence of NMO velocity on the ray parameter p (the ‘dip-moveout (DMO) signature’). The DMO signature retains the same character as for vertical transverse isotropy only for near-vertical and near-horizontal orientation of the symmetry axis. The behaviour of NMO velocity rapidly changes if the symmetry axis is tilted away from the vertical, with a tilt of ±20° being almost sufficient to eliminate the influence of the anisotropy on the DMO signature. For larger tilt angles and typical positive values of the difference between the anisotropic parameters ε and δ, the NMO velocity increases with p more slowly than in homogeneous isotropic media; a dependence usually caused by a vertical velocity gradient. Dip-moveout processing for a wide range of tilt angles requires application of anisotropic DMO algorithms. The strong influence of the tilt angle on P-wave moveout can be used to constrain the tilt using P-wave NMO velocity in the plane that includes the symmetry axis. However, if the azimuth of the axis is unknown, the inversion for the axis orientation cannot be performed without a 3D analysis of reflection traveltimes on lines with different azimuthal directions. 相似文献
20.
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. 相似文献