共查询到20条相似文献，搜索用时 62 毫秒
1.
Wavefront construction (WFC) methods provide robust tools for computing ray theoretical traveltimes and amplitudes for multivalued wavefields. They simulate a wavefront propagating through a model using a mesh that is refined adaptively to ensure accuracy as rays diverge during propagation. However, an implementation for quasi-shear (qS) waves in anisotropic media can be very difficult, since the two qS slowness surfaces and wavefronts often intersect at shear-wave singularities. This complicates the task of creating the initial wavefront meshes, as a particular wavefront will be the faster qS-wave in some directions, but slower in others. Analogous problems arise during interpolation as the wavefront propagates, when an existing mesh cell that crosses a singularity on the wavefront is subdivided. Particle motion vectors provide the key information for correctly generating and interpolating wavefront meshes, as they will normally change slowly along a wavefront. Our implementation tests particle motion vectors to ensure correct initialization and propagation of the mesh for the chosen wave type and to confirm that the vectors change gradually along the wavefront. With this approach, the method provides a robust and efficient algorithm for modeling shear-wave propagation in a 3-D, anisotropic medium. We have successfully tested the qS-wave WFC in transversely isotropic models that include line singularities and kiss singularities. Results from a VTI model with a strong vertical gradient in velocity also show the accuracy of the implementation. In addition, we demonstrate that the WFC method can model a wavefront with a triplication caused by intrinsic anisotropy and that its multivalued traveltimes are mapped accurately. Finally, qS-wave synthetic seismograms are validated against an independent, full-waveform solution. 相似文献
2.
Inversion of Travel Times in Weakly Anisotropic Rock Samples 总被引：1，自引：0，他引：1
Based on the perturbation theory, inversion formulae for travel time of qP and qS waves in arbitrary weak anisotropy media are presented. The inversion formulae are linear expressions of elastic parameters expressed in terms of weak anisotropy (WA) parameters. The formulae of qS1 and qS2 waves have the same form and they can be used without identifying which wave is considered. A synthetic experiment similar to the measurement of rock sample in the laboratory is carried out to illustrate the efficiency of the presented inversion formulae. Two data sets for qP wave travel time from rock samples in the laboratory are inverted and 15 WA parameters are obtained. 相似文献
3.
Analytic expressions for the velocity sensitivity to the elastic moduli for the most general anisotropic media 总被引：1，自引：0，他引：1
For non‐linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density‐normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS1, qS2). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP‐, qS1‐ and qS2‐wave computations and mitigates the singularity problem for the two quasi‐shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes. 相似文献
4.
— We propose an algorithm for local evaluation of weak anisotropy (WA) parameters from measurements of slowness vector components and/or of particle motions of q P waves at individual receivers in a borehole in a multi-azimuthal multiple-source offset VSP experiment. As a byproduct the algorithm yields approximate angular variation of q P-wave phase velocity. The formulae are derived under assumption of weak but arbitrary anisotropy and lateral inhomogeneity of the medium. The algorithm is thus independent of structural complexities between the source and the receiver. If complete slowness vector is determinable from observed data, then the information about polarization can be used as an independent additional constraint. If only the component of the slowness along the borehole can be determined from observations (which is mostly the case), the inversion without information about polarization is impossible. We present several systems of equations which can be used when different numbers of components of the slowness vector are available. The SVD algorithm is used to solve an overdetermined system of linear equations for WA parameters for two test examples of synthetic multi-azimuthal multiple-source offset VSP data. The system of equations results from approximate first-order perturbation equations for the slowness and polarization vectors of the q P wave. Analysis of singular values and of variances of WA parameters is used for the estimation of chances to recover the sought parameters. Effects of varying number of profiles with sources and of noise added to “observed” data are illustrated. An important observation is that although, due to insufficient data, we often cannot recover all individual WA parameters with sufficient accuracy, angular phase velocity variation can be recovered rather well. 相似文献
5.
6.
以二维情形下观测速度场为各向同性场和各向异性场的叠加为前提, 提出了一种利用走时残差估算地震波速度各向异性的方法, 即剩余慢度矢量法. 利用小江断裂带北段巧家流动地震台阵24个台站记录的3181次地震事件的P波走时残差, 采用剩余慢度矢量法计算了各观测台站周围水平方向上尺度为0.5°×0.5°, 震源深度为0—5 km的剩余慢度矢量, 由此得到了P波快波和慢波方向. 计算结果表明, 大部分观测台站周围的P波速度方向性较为一致, 快波方向为ESE向, 慢波方向为NNE向. 快波方向与小江断裂带北段应力场P轴方向较为一致, 而慢波方向与应力场T轴方向一致, 表明应力的长期作用可能是导致P波速度各向异性的重要原因. 相似文献
7.
V. Červený 《Studia Geophysica et Geodaetica》2007,51(3):391-410
The reflection/transmission laws (R/T laws) of plane waves at a plane interface between two homogeneous anisotropic viscoelastic
(dissipative) halfspaces are discussed. Algorithms for determining the slowness vectors of reflected/transmitted plane waves
from the known slowness vector of the incident wave are proposed. In viscoelastic media, the slowness vectors of plane waves
are complex-valued, p = P + iA, where P is the propagation vector, and A the attenuation vector. The proposed algorithms may be applied to bulk plane waves (A = 0), homogeneous plane waves (A ≠ 0, P and A parallel), and inhomogeneous plane waves (A ≠ 0, P and A non-parallel). The manner, in which the slowness vector is specified, plays an important role in the algorithms. For unrestricted
anisotropy and viscoelasticity, the algorithms require an algebraic equation of the sixth degree to be solved in each halfspace.
The degree of the algebraic equation decreases to four or two for simpler cases (isotropic media, plane waves in symmetry
planes of anisotropic media). The physical consequences of the proposed algorithms are discussed in detail.
vcerveny@seis.karlov.mff.cuni.cz 相似文献
8.
Acoustic full waveforms recorded in wells are the simplest way to get the velocity of P, S, and Stoneley waves in situ. Processing and interpretation of acoustic full waveforms in hard formations does not generate problems with identification
packets of waves and calculation of their slowness and arrivals, and determination of the elastic parameter of rocks. But
in shallow intervals of wells, in soft formations, some difficulties arise with proper evaluation of the S-wave velocity due to the lack of refracted S wave in case when its velocity is lower than the velocity of mud. Dynamic approach to selection of a proper value of semblance
to determine the proper slowness and arrival is presented. Correlation between the results obtained from the proposed approach
and the theoretical modeling is a measure of the correctness of the method. 相似文献
9.
Refracted arrivals are analysed to estimate the near‐surface anisotropy of marine sediments using a vertical‐cable (VC) configuration. In the presence of dip, the horizontal and vertical ray‐slownesses are obtained from the observed apparent slownesses in the up‐ and downdip directions using a sum or difference at each azimuth. The multiple azimuths generated by a VC geometry permit the ray‐slowness distribution of the marine sediments to be determined. An inversion procedure is developed to provide dip and anisotropy parameters for refractive layers from the measured refraction traveltimes in multilayered azimuthally isotropic and anisotropic media. Two sets of transversely isotropic models are used to analyse the azimuthal variations of apparent and ray slownesses. In the first set, we fix the anisotropic parameters of the models but vary the dip (0°, 5° and 10°) to test the effects of the presence of dip. In the second set, we vary the P‐wave anisotropy strength (5.2%, 10.3%, 15.8% and 22.0%) to examine the sensitivity and accuracy of ray‐slowness approximations which are independent of dip. We test this inversion procedure on synthetic P‐wave VC data calculated for six different models by a finite‐difference method. The results of applications to real VC data acquired from the North Sea are also presented. 相似文献
10.
In viscoelastic media, the slowness vector p of plane waves is complex-valued, p = P + iA. The real-valued vectors P and A are usually called the propagation and the attenuation vector, repectively. For P and A nonparallel, the plane wave is called inhomogeneousThree basic approaches to the determination of the slowness vector of an inhomogeneous plane wave propagating in a homogeneous viscoelastic anisotropic medium are discussed. They differ in the specification of the mathematical form of the slowness vector p. We speak of directional specification, componental specification and mixed specification of the slowness vector. Individual specifications lead to the eigenvalue problems for 3 × 3 or 6 × 6 complex-valued matrices.In the directional specification of the slowness vector, the real-valued unit vectors N and M in the direction of P and A are assumed to be known. This has been the most common specification of the slowness vector used in the seismological literature. In the componental specification, the real-valued unit vectors N and M are not known in advance. Instead, the complex-valued vactorial component p
of slowness vector p into an arbitrary plane with unit normal n is assumed to be known. Finally, the mixed specification is a special case of the componental specification with p
purely imaginary. In the mixed specification, plane represents the plane of constant phase, so that N = ±n. Consequently, unit vector N is known, similarly as in the directional specification. Instead of unit vector M, however, the vectorial component d of the attenuation vector in the plane of constant phase is known.The simplest, most straightforward and transparent algorithms to determine the phase velocities and slowness vectors of inhomogeneous plane waves propagating in viscoelastic anisotropic media are obtained, if the mixed specification of the slowness vector is used. These algorithms are based on the solution of a conventional eigenvalue problem for 6 × 6 complex-valued matrices. The derived equations are quite general and universal. They can be used both for homogeneous and inhomogeneous plane waves, propagating in elastic or viscoelastic, isotropic or anisotropic media. Contrary to the mixed specififcation, the directional specification can hardly be used to determine the slowness vector of inhomogeneous plane waves propagating in viscoelastic anisotropic media. Although the procedure is based on 3 × 3 complex-valued matrices, it yields a cumbersome system of two coupled equations. 相似文献
11.
Despite the complexity of wave propagation in anisotropic media, reflection moveout on conventional common-midpoint (CMP) spreads is usually well described by the normal-moveout (NMO) velocity defined in the zero-offset limit. In their recent work, Grechka and Tsvankin showed that the azimuthal variation of NMO velocity around a fixed CMP location generally has an elliptical form (i.e. plotting the NMO velocity in each azimuthal direction produces an ellipse) and is determined by the spatial derivatives of the slowness vector evaluated at the CMP location. This formalism is used here to develop exact solutions for the NMO velocity in anisotropic media of arbitrary symmetry. For the model of a single homogeneous layer above a dipping reflector, we obtain an explicit NMO expression valid for all pure modes and any orientation of the CMP line with respect to the reflector strike. The contribution of anisotropy to NMO velocity is contained in the slowness components of the zero-offset ray (along with the derivatives of the vertical slowness with respect to the horizontal slownesses) — quantities that can be found in a straightforward way from the Christoffel equation. If the medium above a dipping reflector is horizontally stratified, the effective NMO velocity is determined through a Dix-type average of the matrices responsible for the ‘interval’ NMO ellipses in the individual layers. This generalized Dix equation provides an analytic basis for moveout inversion in vertically inhomogeneous, arbitrarily anisotropic media. For models with a throughgoing vertical symmetry plane (i.e. if the dip plane of the reflector coincides with a symmetry plane of the overburden), the semi-axes of the NMO ellipse are found by the more conventional rms averaging of the interval NMO velocities in the dip and strike directions. Modelling of normal moveout in general heterogeneous anisotropic media requires dynamic ray tracing of only one (zero-offset) ray. Remarkably, the expressions for geometrical spreading along the zero-offset ray contain all the components necessary to build the NMO ellipse. This method is orders of magnitude faster than multi-azimuth, multi-offset ray tracing and, therefore, can be used efficiently in traveltime inversion and in devising fast dip-moveout (DMO) processing algorithms for anisotropic media. This technique becomes especially efficient if the model consists of homogeneous layers or blocks separated by smooth interfaces. The high accuracy of our NMO expressions is illustrated by comparison with ray-traced reflection traveltimes in piecewise-homogeneous, azimuthally anisotropic models. We also apply the generalized Dix equation to field data collected over a fractured reservoir and show that P-wave moveout can be used to find the depth-dependent fracture orientation and to evaluate the magnitude of azimuthal anisotropy. 相似文献
12.
地下介质中普遍存在着各向异性,当前基于各向异性的地震波射线追踪多是在弱各向异性介质中进行且采用群速度近似表示方法,这些近似方法在强各项异性介质中会导致很大误差而无法真正模拟地震波的传播规律。根据地下普遍存在各向异性的事实和地震波基本传播规律,提出利用牛顿迭代法高效求解群速度,基于Paraview平台自动化构建三维地质模型,采用最短路径法进行地震波射线追踪模拟及可视化,实现对复杂三维地质的速度不均匀性和各向异性的表达,为三维地质模型的构建和地震波射线追踪模拟及可视化提供一种新思路,并以华北克拉通山西断陷带北部局部区域为例进行研究。结果表明,该方法能够减少由各向异性对地震波传播模拟造成的影响,清晰表达了研究区地质结构和各向异性特点,在对复杂三维地质结构的解读中能够较好应用。 相似文献
13.
由于构造运动等作用,TI介质对称轴往往沿空间任意方向分布,具有任意空间取向对称轴的TI（ATI）介质更符合实际地质情况.VTI介质与ATI介质的相速度在形式上具有一致性,VTI介质中地震波的相角对应ATI介质对称轴与地震波传播方向的夹角.本文基于Tsvankin的VTI介质精确相速度公式,利用TI介质对称轴和地震波传播方向上单位向量的数量积和向量积来计算ATI介质的精确相速度.根据弱各向异性假设,导出qP波和qSV波的近似相速度,分析了近似公式的误差,讨论总结了ATI介质qP波和qSV波的相速度特征.本文中的单位向量采用观测坐标系表示,通过相角关系,可以较为方便地由ATI介质近似相速度导出频散关系,然后借助傅里叶逆变换推导出时间-波数域qP波和qSV波解耦的波动方程.数值算例表明本文的波动方程是qP波和qSV波解耦的,波场计算结果稳定,未出现明显的数值频散,验证了本文方法的有效性. 相似文献
14.
采用三维有限差分方法模拟了正交偶极子声源在含偏心钻铤的充液井孔中激发的声场,研究了钻铤偏心对模式波的种类、激发幅度、以及频散特征的影响.研究结果表明,钻铤偏心导致偶极子声源激发的声场的模式不唯一,观察到了除偶极模式外的单极模式波和四极模式波;钻铤偏心导致偶极模式波出现分裂现象,尤其是快速地层F2模式和慢速地层弯曲波的稍高频率的部分,且快、慢波所对应的两个方位为偏心的方位和与偏心方位垂直的方向;井孔折射横波以及快速地层F2模式的低频部分的慢度基本未受到钻铤偏心的影响,仍然能够正确反映地层的横波慢度及各向异性;对于本文研究的慢速地层井孔模型,当偏心距离l小于等于0.01 m时,弯曲波的慢度和激发幅度受钻铤偏心的影响很小,从快、慢弯曲波中提取的快、慢横波慢度基本能够反映地层的各向异性特征. 相似文献
15.
S-wave splitting from records of local micro-earthquakes in West Bohemia/Vogtland: An indicator of complex crustal anisotropy 总被引：1，自引：1，他引：0
Most S-wave particle motions of local micro-earthquakes in the West Bohemia/Vogtland region display S-wave splitting. The
split S waves are usually well defined, being separated in time and polarized in roughly perpendicular directions in the horizontal
projection. In most cases, the polarization of the fast S wave is aligned NW-SE (referred to as “normal splitting”), which
is close to the direction of the maximum horizontal compression in the region. However, for some ray directions, the polarization
of the fast S wave is aligned NE-SW (referred to as “reverse splitting”). The pattern of normal/reverse splitting on a focal
sphere is station-dependent, indicating the presence of inhomogeneities in anisotropy. For some stations, the normal/reverse
splitting pattern is asymmetric with respect to the vertical axis, indicating the symmetry axes of anisotropy are probably
inclined. The presence of inclined anisotropy is confirmed by observations of directionally dependent delay times between
split S waves. A complex and station-dependent anisotropy pattern is probably the result of a complicated anisotropic crust
characterized by diverse geological structures. The spatial variation of anisotropy probably reflects the presence of a variety
of different types of anisotropic rocks in the region. 相似文献
16.
地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质. 相似文献
17.
几种弱向异性介质中qSV波速度和偏振异常 总被引：2，自引：2，他引：0
研究了在具有生趣对称轴的横向各向同性介质（ＴＩＶ介质）中ｑＳＶ波群速度及偏振矢量,给出了相应的精确和近似公式,进一步讨论了用较简单的近似公式来代替复杂精确公式的可靠性;最后展示了地球内部几种常见各向异性岩石矿物中地震ｓＳＶ波速度各向异性因子、偏振信息及其与传播方向之间的偏差。 相似文献
18.
19.
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. 相似文献
20.
Yuriy Ivanov 《Geophysical Prospecting》2019,67(9):2287-2297
Degeneracies of the slowness surfaces of shear (and compressional) waves in low-symmetry anisotropic media (such as orthorhombic), known as point singularities, pose difficulties during modelling and inversion, but can be potentially used in the latter as model parameter constraints. I analyse the quantity and spatial arrangement of point singularities in orthorhombic media, as well as their relation to the overall strength of velocity anisotropy. A classification scheme based on the number and spatial distribution of singularity directions is proposed. In normal orthorhombic models (where the principal shear moduli are smaller than the principal compressional moduli), point singularities can only be arranged in three distinct patterns, and media with the theoretical minimum (0) and maximum (16) number of singularities are not possible. In orthorhombic models resulting from embedding vertical fractures in transversely isotropic background, only two singularity distributions are possible, in contrast to what was previously thought. Although the total number of singularities is independent of the overall anisotropy strength, for general (non-normal) orthorhombic models, different spatial distributions of singularities become more probable with increasing magnitude of anisotropy. 相似文献