共查询到20条相似文献，搜索用时 359 毫秒
1.
本文基于任意空间取向TI介质中坐标变换的方法,扩展研究了任意强弱、具有任意空间取向对称轴的TI介质中体波相速度和群速度的方位变化;通过模型数值计算,获得了相速度与群速度方位变化图案,并比较了二者偏差.研究发现,各向异性越强群速度与相速度的偏差越大;并且,TI对称轴的空间取向和测线方位影响速度方位变化.相对于TI对称轴,速度方位变化图案不变;但是,随着TI对称轴空间取向和测线方位的改变,其速度方位变化呈现一定的规律性.研究结果可以精确地预测任意空间取向TI介质中速度方位变化,有利于地震各向异性数据处理与解释. 相似文献
2.
同类反射波（非转换波）走时偏离双曲,称为非双曲或四次时差,在长排列各向异性地震资料处理中需要校正．本文基于我们导出的水平界面任意空间取向TI (ATI)介质中同类反射波四次时差系数(A4)的精确解析解,数值计算研究P波四次时差系数特征．正演结果表明,ATI条件下A4系数随CMP测线方位变化的特征不仅与TI介质的各向异性参数有关,而且与TI对称轴的空间取向密切相关; TI介质的各向异性参数和TI对称轴的倾角决定了A4变化特征,而且TI对称轴的方位决定了A4随测线方位变化的对称性．此研究结果将对各向异性解释及参数反演有参考意义． 相似文献
3.
The horizontal transversely isotropic model, with arbitrary symmetry axis orientation, is the simplest effective representative that explains the azimuthal behaviour of seismic data. Estimating the anisotropy parameters of this model is important in reservoir characterisation, specifically in terms of fracture delineation. We propose a travel‐time‐based approach to estimate the anellipticity parameter η and the symmetry axis azimuth ? of a horizontal transversely isotropic medium, given an inhomogeneous elliptic background model (which might be obtained from velocity analysis and well velocities). This is accomplished through a Taylor's series expansion of the travel‐time solution (of the eikonal equation) as a function of parameter η and azimuth angle ?. The accuracy of the travel time expansion is enhanced by the use of Shanks transform. This results in an accurate approximation of the solution of the non‐linear eikonal equation and provides a mechanism to scan simultaneously for the best fitting effective parameters η and ?, without the need for repetitive modelling of travel times. The analysis of the travel time sensitivity to parameters η and ? reveals that travel times are more sensitive to η than to the symmetry axis azimuth ?. Thus, η is better constrained from travel times than the azimuth. Moreover, the two‐parameter scan in the homogeneous case shows that errors in the background model affect the estimation of η and ? differently. While a gradual increase in errors in the background model leads to increasing errors in η, inaccuracies in ?, on the other hand, depend on the background model errors. We also propose a layer‐stripping method valid for a stack of arbitrary oriented symmetry axis horizontal transversely isotropic layers to convert the effective parameters to the interval layer values. 相似文献
4.
5.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
6.
7.
针对粘弹性裂隙型单斜介质,本文应用高阶交错网格差分法,对其传播快照以及多方位地面三分量记录进行模拟,结果证明了裂隙填充物的性质引起的各向异性与粘弹性性质对单斜介质波场有明显的影响,总结了不同观测方位地震记录呈现出的规律性,这为进一步了解波场在介质中的传播情况,提出更合理的理论模型,提供一种尝试. 相似文献
8.
An important cause of seismic anisotropic attenuation is the interbedding of thin viscoelastic layers. However, much less attention has been devoted to layer‐induced anisotropic attenuation. Here, we derive a group of unified weighted average forms for effective attenuation from a binary isotropic, transversely isotropic‐ and orthorhombic‐layered medium in the zero‐frequency limit by using the Backus averaging/upscaling method and analyse the influence of interval parameters on effective attenuation. Besides the corresponding interval attenuation and the real part of stiffness, the contrast in the real part of the complex stiffness is also a key factor influencing effective attenuation. A simple linear approximation can be obtained to calculate effective attenuation if the contrast in the real part of stiffness is very small. In a viscoelastic medium, attenuation anisotropy and velocity anisotropy may have different orientations of symmetry planes, and the symmetry class of the former is not lower than that of the latter. We define a group of more general attenuation‐anisotropy parameters to characterize not only the anisotropic attenuation with different symmetry classes from the anisotropic velocity but also the elastic case. Numerical tests reveal the influence of interval attenuation anisotropy, interval velocity anisotropy and the contrast in the real part of stiffness on effective attenuation anisotropy. Types of effective attenuation anisotropy for interval orthorhombic attenuation and interval transversely isotropic attenuation with a vertical symmetry (vertical transversely isotropic attenuation) are controlled only by the interval attenuation anisotropy. A type of effective attenuation anisotropy for interval TI attenuation with a horizontal symmetry (horizontal transversely isotropic attenuation) is controlled by the interval attenuation anisotropy and the contrast in the real part of stiffness. The type of effective attenuation anisotropy for interval isotropic attenuation is controlled by all three factors. The magnitude of effective attenuation anisotropy is positively correlated with the contrast in the real part of the stiffness. Effective attenuation even in isotropic layers with identical isotropic attenuation is anisotropic if the contrast in the real part of stiffness is non‐zero. In addition, if the contrast in the real part of stiffness is very small, a simple linear approximation also can be performed to calculate effective attenuation‐anisotropy parameters for interval anisotropic attenuation. 相似文献
9.
This study used SKS waveforms from the International Deep Profiling of Tibet and the Himalayas (INDEPTH) III dataset and a new 2D method for modeling seismic waves in anisotropic media to construct an image of anisotropic structures beneath central Tibet. A preferred model revealed three-segment anisotropic structures in the upper mantle beneath the study region. Waveform modeling demonstrated that the anisotropy was mainly generated by the lithosphere but not the asthenosphere, and that an anisotropic model with a flatter axis of symmetry provides a more consistent interpretation of the observations than models having steeply dipping symmetry axes. A relatively low velocity zone may underlie or intermingle with the anisotropic structures in the northern portion of the region. Synthetic tests also indicate that variations in the elastic constants and depth extent of the anisotropy assumed by the calculations do not affect the general conclusions, although trade-offs exist among certain model parameters. The modeling results suggest that the complex seismic structures in central Tibet were associated with underthrusting of the Indian lithosphere beneath the Asian lithosphere; the inferred flat symmetry axis of the anisotropy was likely generated during this collision process. If this were not the case, the inherited anisotropy would exhibit a steeply dipping axis of symmetry, parallel to the direction of underthrusting. 相似文献
10.
11.
非双曲（远偏移距）时差为各向异性介质中正、反演研究,特别是各向异性参数估计提供了重要信息.本文在任意空间取向TI（ATI）介质水平界面同类反射波四次时差系数(A4)精确解的基础上,进一步讨论我们推导得出的ATI介质中四次时差系数解析近似解,比较随CMP测线方位变化的近似解与精确解之间的差别,为利用近似解来解析研究ATI介质中非双曲时距以及参数反演提供有价值的信息．结合实际岩性资料的数值研究表明,ATI条件下四次时差系数近似解与精确解之间存在差别,不仅表现在A4数的大小及符号特征上,更突出地表现在A4系数随方位的变化特征上;在强各向异性条件下,近似解相比精确解存在较大误差．但在各向异性参数满足0<ε-δ<0.15、｜δ｜<0.20的情况下,对于TI对称轴的特殊倾角范围(75°～80°),近似解与精确解的差别很小,可用近似解进行各向异性观测解释及参数反演． 相似文献
12.
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. 相似文献
13.
The conventional intersection method for earthquake location in isotropic media is developed in the case of transversely isotropic
media with a tilted symmetry axis (TTI media). The hypocenter is determined using its loci, which are calculated through a
minimum travel time tree algorithm for ray tracing in TTI media. There are no restrictions on the structural complexity of
the model or on the anisotropy strength of the medium. The location method is validated by its application to determine the
hypocenter and origin time of an event in a complex TTI structure, in accordance with four hypotheses or study cases: (a)
accurate model and arrival times, (b) perturbed model with randomly variable elastic parameter, (c) noisy arrival time data,
and (d) incomplete set of observations from the seismic stations. Furthermore, several numerical tests demonstrate that the
orientation of the symmetry axis has a significant effect on the hypocenter location when the seismic anisotropy is not very
weak. Moreover, if the hypocentral determination is based on an isotropic reference model while the real medium is anisotropic,
the resultant location errors can be considerable even though the anisotropy strength does not exceed 6.10%. 相似文献
14.
Patrick N.J. Rasolofosaon 《Geophysical Prospecting》2010,58(4):637-655
Seismic anisotropy in geological media is now widely accepted. Parametrizations and explicit approximations for the velocities in such media, considered as purely elastic and moderately anisotropic, are now standards and have even been extended to arbitrary types of anisotropy. In the case of attenuating media, some authors have also recently published different parametrizations and velocity and attenuation approximations in viscoelastic anisotropic media of particular symmetry type (e.g., transversely isotropic or orthorhombic). This paper extends such work to media of arbitrary anisotropy type, that is to say to triclinic media. In the case of homogeneous waves and using the so‐called ‘correspondence principle’, it is shown that the viscoelastic equations (for the phase velocities, phase slownesses, moduli, wavenumbers, etc.) are formally identical to the corresponding purely elastic equations available in the literature provided that all the corresponding quantities are complex (except the unit vector in the propagation direction that remains real). In contrast to previous work, the new parametrization uses complex anisotropy parameters and constitutes a simple extension to viscoelastic media of previous work dealing with non‐attenuating elastic media of arbitrary anisotropy type. We make the link between these new complex anisotropy parameters and measurable parameters, as well as with previously published anisotropy parameters, demonstrating the usefulness of the new parametrization. We compute the explicit complete directional dependence of the exact and of the approximate (first and higher‐order perturbation) complex phase velocities of the three body waves (qP, qS1 and qS2). The exact equations are successfully compared with the ultrasonic phase velocities and phase attenuations of the three body waves measured in a strongly attenuating water‐saturated sample of Vosges sandstone exhibiting moderate velocity anisotropy but very strong attenuation anisotropy. The approximate formulas are checked on experimental data. Compared to the exact solutions, the errors observed on the first‐order approximate velocities are small (<1%) for qP‐waves and moderate (<10%) for qS‐waves. The corresponding errors on the quality factor Q are moderate (<6%) for qP‐waves but critically large (up to 160%) for the qS‐waves. The use of higher‐order approximations substantially improves the accuracy, for instance typical maximum relative errors do not exceed 0.06% on all the velocities and 0.6% on all the quality factors Q, for third‐order approximations. All the results obtained on other rock samples confirm the results obtained on this rock. The simplicity of the derivations and the generality of the results are striking and particularly convenient for practical applications. 相似文献
15.
在裂缝诱导各向异性理论研究中通常使用等效HTI介质来近似多组裂缝所引起的综合效应.由于构造运动的复杂性,多组裂缝普遍存在于地壳与油气储层中.为了研究多组裂缝的地震属性特征,分析常用的等效HTI模型对于多组裂缝近似精度及附加裂缝对介质属性特征的影响,本文利用线性滑移模型进行了多组垂直裂缝的单斜各向异性等效介质理论计算,并利用空间搜索方法求取与其最为接近的HTI介质各向异性弹性参数.重点研究了在两种各向异性介质中纵波速度、快慢横波速度和极化特征及其差异,量化分析附加裂缝对于地震属性如速度、极化方向和走时等的影响,研究对附加裂缝敏感的地震属性.此研究结果和方法为进一步研究多组裂缝的反演及识别方法提供基础,同时对于将高阶对称性各向异性介质中已存在的计算方法应用于低阶对称性时的适用程度、精度分析及相关方法研究具有重要作用. 相似文献
16.
对于倾斜叠层和非垂直裂隙岩层,用具有任意空间取向对称轴的TI (ATI)模型来描述更符合实际观测. 本文基于坐标变换的方法,研究任意强弱ATI介质中体波速度的角散和方位变化特征. 研究结果表明,ATI介质中体波速度随传播方向变化的速度图案相对TI对称轴确定,此确定的速度图案与TI的Thomsen参数相关;速度特征只依赖于传播矢量与对称轴的夹角. 因此,随着TI对称轴取向的空间变化和测线方位的变化,体波速度图案呈现多样性变化,并具有一定的对称性、渐变和重复性. 研究结果有助于进一步的理论研究和各向异性资料处理解释. 相似文献
17.
针对裂隙型储集层中更具代表性的各向异性介质模型,即在各向同性背景介质中含有两组斜交的垂直裂隙所构成的单斜各向异性介质模型,利用时间和空间上可达任意阶的高阶交错网格有限差分技术,对具有不同裂隙填充物性质的单斜介质中波的传播快照进行了模拟.结果证实各向异性介质中波的传播速度随传播方向的不同而产生明显的差异;裂隙填充物的性质对于速度各向异性具有很大的影响.另外,利用坐标旋转法,对水平层状各向异性介质中多方位地面三分量记录进行了模拟,结果表明了方位各向异性介质中,波的传播速度不仅随入射角的变化而变化,同时也随观测方位的不同而产生差异.数值模拟结果为进一步利用地面多方位地震属性进行各向异性参数的反演及裂隙参数的描述提供理论基础. 相似文献
18.
具有垂直对称轴的横向各向同性介质模型(VTI)是目前各向异性理论研究和多波多分量地震资料叠前成像处理中最常用的一种各向异性模型.VTI介质中反射 P波时距曲线一般不再是双曲线.基于不同的相速度近似公式会得到不同的时距关系式.文中对几种典型的非双曲时距曲线与射线追踪得到的准确时距曲线在不同各向异性强度下进行了对比，结果表明Muir等和Stovas等提出的非双曲时距公式由于过高地考虑了横波垂直速度的影响与精确的时距曲线有很大偏差；Tsvankin等提出的弱各向异性非双曲时距公式在ε-δ<0时误差增大；Alkhalifah等提出的非双曲时距公式在大炮检距任意各向异性强度下都具有较高的精度，适于在实际资料处理中应用. 相似文献
19.
Fritz Gassmann 《Pure and Applied Geophysics》1964,58(1):63-112
Summary Section 1 (and 11) develops the concepts of the front velocity, the front gradient, the travel time in space and on seismometric profiles, the profile velocity and the profile gradient in connection with the propagation of the fronts of elastic waves in solid isotropic and anisotropic media. The sectional velocity and the sectional gradient are defined in terms of the motion of the curve of intersection of a front with a fixed surface. Section 2 (and 12) relates the coefficients of elasticity of the medium, the front types, and their respective rays. In section 12, the theory of fronts of arbitrary shape and of the corresponding rays for any anisotropic, homogeneous or inhomogeneous solid medium is summarized. In section 3 (and 13), the law of reflection and refraction of fronts on surfaces of discontinuity of arbitrary shape is presented. Sections 4 to 6 (and 14 to 16) treat some elementary applications of seismic travel time methods to homogeneous, uniaxially anisotropic media (=transverse isotropy) in greater detail. In section 4 (and 14), the travel time of a direct front generated by a point source is considered and it is shown how the coefficients of elasticity of the medium can be found based on travel time measurements. The seismic prospection of a plane reflector and of a reflecting boundary of arbitrary shape and position are discussed in section 5 (and 15). In section 6 (and 16), the seismic refraction method is used to locate a plane boundary between a homogeneous, uniaxially anisotropic and a homogeneous isotropic medium, where the boundary is perpendicular or at an arbitrary angle to the direction of anisotropy. 相似文献
20.
Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations 总被引：2，自引：0，他引：2
In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the
elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying
the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover,
the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or
crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the
non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times
and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave
velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out
a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results
derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first
procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion
method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern
China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy
signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia
blocks. 相似文献