共查询到20条相似文献，搜索用时 31 毫秒
1.
Determination of the ray vector (the unit vector specifying the direction of the group velocity vector) corresponding to a given wave normal (the unit vector parallel to the phase velocity vector or slowness vector) in an arbitrary anisotropic medium can be performed using the exact formula following from the ray tracing equations. The determination of the wave normal from the ray vector is, generally, a more complicated task, which is usually solved iteratively. We present a first-order perturbation formula for the approximate determination of the ray vector from a given wave normal and vice versa. The formula is applicable to qP as well as qS waves in directions, in which the waves can be dealt with separately (i.e. outside singular directions of qS waves). Performance of the approximate formulae is illustrated on models of transversely isotropic and orthorhombic symmetry. We show that the formula for the determination of the ray vector from the wave normal yields rather accurate results even for strong anisotropy. The formula for the determination of the wave normal from the ray vector works reasonably well in directions, in which the considered waves have convex slowness surfaces. Otherwise, it can yield, especially for stronger anisotropy, rather distorted results. 相似文献
2.
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. 相似文献
3.
提出了利用人工爆破P波走时反演地壳介质方位各向异性参数的方法. 在假定介质是弱各向异性介质的情况下, 使用扰动理论得到了线性化的反演公式, 其中待反演的弱各向异性参数是P波走时的线性函数. 如果在反演公式中参考走时取相同震中距接收点的P波平均走时, 那么所获得的弱各向异性参数与参考介质速度的选取无关. 反演得到的弱各向异性参数可以看作是不同震中距和不同深度范围内介质的等效弱各向异性参数. 等效弱各向异性参数在一定程度上反映了不同深度范围内水平方向相速度随方位的变化. 这种变化可能是不同时期构造应力作用的结果. 2007年中国地震局在首都圈怀来地区实施了一次大吨位人工爆破实验, 以爆破点为中心, 布设了高密度的地震观测台网和台阵. 台站相对于爆破点具有360°的全方位覆盖, 所得到的地震记录数据为研究怀来、 延庆地区地壳介质P波方位各向异性提供了必要条件. 我们通过走时反演获得了与水平方位相关的弱各向异性参数, 并对弱各向异性参数进行坐标变换, 得到了能够直观描述岩石弱各向异性的具有水平对称轴的横向各向同性介质, 给出了对应的3个独立弱各向异性参数及其对称轴方位, 讨论了介质各向异性与构造应力场的关系. 结果表明该地区地壳介质存在明显的方位各向异性, 其最大值约为4.6%. 相似文献
4.
5.
— 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. 相似文献
6.
— A P-wave tomographic method for 3-D complex media (3-D distribution of elastic parameters and curved interfaces) with orthorhombic symmetry is presented in this paper. The technique uses an iterative linear approach to the nonlinear travel-time inversion problem. The hypothesis of orthorhombic anisotropy and 3-D inhomogeneity increases the set of parameters describing the model dramatically compared to the isotropic case. Assuming a Factorized Anisotropic Inhomogeneous (FAI) medium and weak anisotropy, we solve the forward problem by a perturbation approach. We use a finite element approach in which the FAI medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, analytical expressions for rays and travel times, valid to first-order, are given for P waves in orthorhombic inhomogeneous media. More complex media can be modeled by introducing interfaces separating FAI media with different elastic properties. Simple formulae are given for the Fréchet derivatives of the travel time with respect to the elastic parameters and the interface parameters. In the weak anisotropy hypothesis the P-wave travel times are sensitive only to a subset of the orthorhombic parameters: the six P-wave elastic parameters and the three Euler angles defining the orientation of the mirror planes of symmetry. The P-wave travel times are inverted by minimizing in terms of least-squares the misfit between the observed and calculated travel times. The solution is approached using a Singular Value Decomposition (SVD). The stability of the inversion is ensured by making use of suitable a priori information and/or by applying regularization. The technique is applied to two synthetic data sets, simulating simple Vertical Seismic Profile (VSP) experiments. The examples demonstrate the necessity of good 3-D ray coverage when considering complex anisotropic symmetry. 相似文献
7.
利用微扰理论推导了弱各向异性(WA)介质参数正反演计算的基本公式;给出了在已知慢度矢量的一个分量和偏振矢量情况下确定WA参数的方法.如果这一方法被用于单一的变井源距垂直地震剖面(walkawayVSP)资料,可以获得9个WA参数.这9个WA参数完全决定了qP波和两个qS波在由剖面和钻井所决定的平面内传播的特性.对单条walkawayVSP观测系统数据的产生和WA参数的反演进行了数值模拟计算,对所能确定的WA参数及其可靠性进行了详细的讨论. 相似文献
8.
—Approximate PP plane wave displacement coefficients of reflection and transmission for weak contrast interfaces separating weakly but arbitrarily anisotropic elastic media are presented. The PP reflection coefficient for such an interface has been derived recently by Vavry?uk and P?en?ík (1997). The PP transmission coefficient presented in this paper was derived by the same approach. The coefficients are given as a sum of the coefficient for the weak contrast interface separating two nearby isotropic media and a term depending linearly on contrasts of the so-called weak anisotropy (WA) parameters (parameters specifying deviation of properties of the medium from isotropy), across the interface. While the reflection coefficient depends only on 8 of the complete set of the WA parameters describing P-wave phase velocity in weakly anisotropic media, the transmission coefficient depends on their complete set. The PP reflection coefficient depends on "shear-wave splitting parameter" γ. Tests of accuracy of the approximate formulae are presented on several models. 相似文献
9.
Wave propagation in a finely layered medium is a very important topic in seismic modelling and inversion. Here we analyse non‐vertical wave propagation in a periodically layered transversely isotropic (VTI) medium and show that the evanescent (attenuation) zones in the frequency‐horizontal slowness domain result in caustics in the group velocity domain. These caustics, which may appear for both the quasi‐compressional (qP) and quasi‐shear (qSV) wave surfaces are frequency dependent but display weak dependence at low frequencies. The caustics computed for a specific frequency differ from those observed at the low‐ and high‐frequency limits. We illustrate these caustics with a few numerical examples and snapshots computed for both qP‐ and qSV‐wave types. 相似文献
10.
11.
Local Determination of Weak Anisotropy Parameters from Walkaway VSP qP-Wave Data in the Java Sea Region 总被引：1，自引：0，他引：1
We apply the inversion scheme of Zheng and Peník (2002) to the walkaway VSP data of Horne and Leaney (2000) collected in the Java Sea region. The goal is a local determination of parameters of the medium surrounding the borehole receiver array. The inversion scheme is based on linearized equations expressing qP-wave slowness and polarization vectors in terms of weak anisotropy (WA) parameters. It thus represents an alternative approach to Horne and Leaney (2000), who based their procedure on inversion of the Christoffel equation using a global optimization method. The presented inversion scheme is independent of structural complexities in the overburden and of the orientation of the borehole. The inverson formula is local, and has therefore potential to separate effects of anisotropy from effects of inhomogeneity. The data used are components of the slowness vector along the receiver array and polarization vectors. The inversion is performed without any assumptions concerning the remaining components of the slowness vector. The inversion is made (a) assuming arbitrary anisotropy, i.e., without any assumptions about symmetry of the medium, (b) assuming transverse isotropy with a vertical axis of symmetry and (c) assuming isotropy of the medium. Inverted are the raw data as well as data, in which weighting is used to reduce the effect of outliers. It is found that the WA parameters
z,
15 and
35 are considerably more stable than the parameters
x and
x. The latter two parameters are also found to be strongly correlated. Weaker correlation is also found between the mentioned two parameters and
z. The results of inversion show clearly that the studied medium is not isotropic. They also seem to indicate that the studied medium does not possess the VTI symmetry. 相似文献
12.
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. 相似文献
13.
A technique allowing inversion of the shale stiffness tensor from standard logging data: sonic velocities, density, porosity and clay content is developed. The inversion is based on the effective medium theory. The testing of the technique on laboratory measurements of the elastic wave velocities in shale samples shows that the inversion makes it possible to predict the elastic wave velocities VP, VS1 and VS2 in any direction within an error of a few per cent. The technique has been applied for the stiffness tensor inversion along a well penetrating a shale formation of the Mississippian age altered by thin layers of limestone. It is demonstrated that the symmetry of a stiffness tensor inverted at the sonic frequency (2 kHz) is slightly orthorhombic and taking into account the experimental errors, can be related to the vertical transverse isotropy symmetry. For the productive interval of the shale formation, the Thomsen parameters ?, γ, and δ average, respectively, 0.32, 0.25 and 0.21, which indicate anelliptic behaviour of the velocities in this shale. The coefficients of anisotropy of this shale interval are around 24% and 20% for the compressional and shear waves, respectively. The values of the inverted velocities in the bedding plane for this interval are in good agreement with the laboratory measurements. The technique also allows inversion of the water saturation of the formation (Sw) and the inverted values are in agreement with the Sw values available for this formation. A Backus‐like upscaling of the inverted stiffness tensors is carried out for the lower and upper bounds of the frequency band used in the crosswell tomography (100 Hz and 500 Hz). These results can serve as an initial velocity model for the microearthquake location during hydrofracking of the shale formation. 相似文献
14.
Preliminary reference Earth model 总被引：29，自引：0，他引：29
A large data set consisting of about 1000 normal mode periods, 500 summary travel time observations, 100 normal mode Q values, mass and moment of inertia have been inverted to obtain the radial distribution of elastic properties, Q values and density in the Earth's interior. The data set was supplemented with a special study of 12 years of ISC phase data which yielded an additional 1.75 × 106 travel time observations for P and S waves. In order to obtain satisfactory agreement with the entire data set we were required to take into account anelastic dispersion. The introduction of transverse isotropy into the outer 220 km of the mantle was required in order to satisfy the shorter period fundamental toroidal and spheroidal modes. This anisotropy also improved the fit of the larger data set. The horizontal and vertical velocities in the upper mantle differ by 2–4%, both for P and S waves. The mantle below 220 km is not required to be anisotropic. Mantle Rayleigh waves are surprisingly sensitive to compressional velocity in the upper mantle. High Sn velocities, low Pn velocities and a pronounced low-velocity zone are features of most global inversion models that are suppressed when anisotropy is allowed for in the inversion.The Preliminary Reference Earth Model, PREM, and auxiliary tables showing fits to the data are presented. 相似文献
15.
We present a parsimonious wave‐equation travel‐time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave‐equation travel‐time inversion. Numerical results show that the parsimonious wave‐equation travel‐time tomogram has about the same accuracy as the tomogram computed by standard wave‐equation travel‐time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location. 相似文献
16.
本文基于弱各向异性（WA）介质的正反演公式和qP波的坐标变换,推导了利用qP波反演任意倾斜对称轴的横向各向同性（TTI）介质的各向异性参数和对称轴方向的公式.理论和数值实验表明,利用2个相互正交的变井源距垂直地震剖面（walkaway VSP）可以完全确定钻井中TTI介质qP波的3个WA参数和对称轴的2个方向参数.我们完成了几个由不同数量剖面组成的walkaway VSP模拟实验,使用TTI模型和一般各向异性模型对模拟数据进行了反演,证明了反演公式的正确性和可靠性.使用这些公式,对来自Java Sea的由3条剖面组成的walkaway VSP观测数据进行了各向异性反演,获得了钻井中接收点处介质的WA参数. 相似文献
17.
In this paper, we have considered the reflection and refraction of a plane wave at an interface between two half-spaces. The
lower half-spaces is composed of highly anisotropic triclinic crystalline material and the upper half-space is homogeneous
and isotropic. It has been assumed that due to incidence of a plane quasi-P (qP) wave, three types of waves, namely, quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated in the lower half space whereas P and S waves will be generated in the upper half space. The phase velocities of all the quasi waves have been calculated. It has
been assumed that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation.
Some specific relations have been established between directions of motion and propagation, respectively. The expressions
for reflection coefficients of qP, qSV, qSH and refracted coefficients of P and SV waves are obtained. Results of reflection and refraction coefficients are presented. 相似文献
18.
提出了一种使用慢度矢量分量和偏振矢量计算变井源距垂直地震剖面(walkaway VSP)钻孔中接收点附近介质弱各向异性(WA)参数的方法. 假定介质是任意弱各向异性介质, 从一般公式中得到了只有一条观测剖面情况下的反演公式. 如果知道了慢度矢量的垂直分量和偏振矢量, 可以通过反演得到与剖面和钻孔所在平面相关的WA参数, 反演过程不用进行射线追踪, 与上覆介质无关. 用合成数据检验了公式和方法的正确性, 并把它们应用于在爪哇海地区得到的一条变井源距垂直地震剖面的弱各向异性参数反演中. 相似文献
19.
20.
本文从一般弱各向异性介质参数反演中得到了使用两条相互正交的变偏VSP测线计算弱各向异性(WA)参数的反演公式. 如果仅仅使用qP波, 则可以确定9个独立的WA参数, 这9个WA参数可以完全地描述井中接收点在两个剖面内介质的各向异性性质. 通过对走时曲线进行最小二乘条件下的三次样条光滑, 可以获得慢度矢量的垂直分量和水平分量. 如果介质是横向非均匀介质, 则水平分量仅仅被用作反演时的约束条件. 为了获得偏振矢量, 本文引入质心计算方法, 该方法计算简单、 稳定, 而且不存在解的奇异问题. 在寻找与一般各向异性介质最接近的高对称性正交各向异性介质和TTI各向异性介质时, 使用qP波各向异性坐标变换方法和最小二乘求解方法, 得到了与一般各向异性介质最接近的正交各向异性和TTI各向异性参数及其对称轴方向参数的计算公式. 使用这些方法, 对瓜哇海地区布设的两条相互正交的变偏VSP测线数据进行各向异性反演, 获得了井中10个接收点处介质的WA参数. 数值计算和实际资料反演表明, 本文所使用的反演方法能够准确地得到VSP井中接收点处介质的WA参数, 这是地震勘探中研究地壳介质各向异性性质最直接和最可靠的方法. 相似文献