共查询到20条相似文献，搜索用时 31 毫秒
1.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
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.
Huajian Yao 《地震科学（英文版）》2015,28(1):59-69
Seismic anisotropy provides important constraints on deformation patterns of Earth's material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm(NA) for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depthdependent VSV(or the elastic parameter L) as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polarization azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle. 相似文献
4.
相速度和群速度是研究地震波传播规律和描述介质特性的重要参数,是弹性波传播理论中的核心内容,在理论研究和实际应用中有重要作用.本文根据VTI介质的刚度矩阵,利用Bond变换建立了TTI介质刚度矩阵.再利用TTI介质刚度矩阵,结合弹性动力学的本构方程、牛顿运动微分方程和几何方程,得到了三维TTI介质弹性波波动方程和Christoffel方程.通过本征值方法求解Christoffel方程,推导了三维TTI介质弹性波相速度的解析表达式.利用Berryman和Crampin推导各向异性介质群速度公式,根据三维TTI介质的相速度解析式推导了三维TTI介质群速度解析表达式.数值试例表明,随着各向异性介质参数改变,TI介质弹性波相速度变化较为平缓,群速度变化较为剧烈,qP波和SH波速度变化较为平缓,qSV波速度变化较为剧烈. 相似文献
5.
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. 相似文献
6.
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. 相似文献
7.
8.
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. 相似文献
9.
黏弹各向异性介质中传播不均匀波，其反射、透射模式不仅与介质分界面两侧速度对比有关，还与品质因子Q的对比有关. 用伪谱技术模拟黏弹各向异性介质分界面上波的反射、透射，并与弹性各向异性介质、黏弹各向同性介质和弹性各向同性介质的模拟结果做比较. 计算平面波的反射、透射系数，分析介质的黏弹性和各向异性对反射、透射系数的影响. 数值模拟了一个三层介质模型中的波场，分析两个分界面上产生的反射波的特征. 黏弹各向异性介质中，qS波比qP波衰减程度大. 相似文献
10.
目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础. 相似文献
11.
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好. 相似文献
12.
13.
地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性. 相似文献
14.
The effect of a hydrous phase on P‐wave velocity anisotropy within a detachment shear zone in the slow‐spreading oceanic crust: A case study from the Godzilla Megamullion,Philippine Sea 
下载免费PDF全文
下载免费PDF全文 Katsuyoshi Michibayashi Tohru Watanabe Yumiko Harigane Yasuhiko Ohara 《Island Arc》2016,25(3):209-219
We studied the contributions of plagioclase, clinopyroxene, and amphibole to the P‐wave velocity properties of gabbroic mylonites of the Godzilla Megamullion (site KH07‐02‐D18) in the Parece Vela Rift of the central Parece Vela Basin, Philippine Sea, based on their crystal‐preferred orientations (CPOs), mineral modes, and elastic constants and densities of single crystals. The gabbroic mylonites have been classified into three types based on their microstructures and temperature conditions: HT1, HT2 and medium‐temperature (MT) mylonites. The P‐wave velocity properties of the HT1 mylonite are dominantly influenced by plagioclase CPOs. Secondary amphibole occurred after deformation in the HT1 mylonite, so that its effect on P‐wave velocity anisotropy is minimal due to weak CPOs. Although the HT2 mylonite developed deformation microstructures in the three minerals, the P‐wave velocity properties of the HT2 mylonite are essentially isotropic, resulting from the destructive interference of different P‐wave velocity anisotropy patterns produced by the distinct CPOs of the three constituent minerals (i.e., plagioclase, clinopyroxene, and amphibole). The P‐wave velocity properties of the MT mylonite are influenced mainly by amphibole CPOs, whereas the effect of plagioclase CPOs on P‐wave velocity anisotropy becomes very small with a decrease in the intensity of plagioclase CPOs. As a result, the gabbroic mylonites tend to have weak P‐wave velocity anisotropy in seismic velocity, although their constituent minerals show distinct CPOs. Such weakness in the whole‐rock P‐wave velocity anisotropy could result from the destructive contributions of the different mineral CPOs with respect to the structural framework (foliation and lineation). These results show that amphibole has a high potential for P‐wave velocity anisotropy by aligning both crystallographically and dimensionally during deformation in the hydrous oceanic crust. The results also suggest that the effect of a hydrous phase on P‐wave velocity anisotropy within the detachment shear zone in a slow‐spreading oceanic crust varies depending on the degree of deformation and on the timing of hydrothermal activity. 相似文献
15.
本文对声波与弹性波方程进行有限元法离散,构造有限元法频散关系的一般特征值问题,分析了时间离散格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性. 比较了三种质量矩阵即分布式质量矩阵、集中质量矩阵和混合质量矩阵对有限元法频散的影响;选取四种典型三角网格,分析了混合质量矩阵有限元（MFEM）频散的方向各向异性;数值频散、方向各向异性随插值阶数的增加逐渐减弱,当空间为三阶插值时,频散主要表现为随采样率的变化而几乎无明显方向各向异性, 其频散幅值也较小. 控制其他影响因素不变的情况下,研究了不同波速比介质中弹性波的数值频散. 最后给出了三角网格MFEM的数值耗散性. 相似文献
16.
Characterisation and multifaceted anisotropy assessment of Corvio sandstone for geological CO2 storage studies 下载免费PDF全文
下载免费PDF全文 Ismael Falcon‐Suarez Jacobo Canal‐Vila Jordi Delgado‐Martin Laurence North Angus Best 《Geophysical Prospecting》2017,65(5):1293-1311
We present a comprehensive characterisation of the physical, mineralogical, geomechanical, geophysical, and hydrodynamic properties of Corvio sandstone. This information, together with a detailed assessment of anisotropy, is needed to establish Corvio sandstone as a useful laboratory rock‐testing standard for well‐constrained studies of thermo–hydro–mechanical–chemical coupled phenomena associated with CO2 storage practices and for geological reservoir studies in general. More than 200 core plugs of Corvio sandstone (38.1 and 50 mm diameters, 2:1 length‐to‐diameter ratio) were used in this characterisation study, with a rock porosity of 21.7 ± 1.2%, dry density 2036 ± 32 kg m?3, and unconfined compressive and tensile strengths of 41 ± 3.28 and 2.3 ± 0.14 MPa, respectively. Geomechanical tests show that the rock behaves elastically between ～10 and ～18 MPa under unconfined conditions with associated Young's modulus and Poisson's ratio of 11.8 ± 2.8 GPa and 0.34 ± 0.01 GPa, respectively. Permeability abruptly decreases with confining pressure up to ～10 MPa and then stabilises at ～1 mD. Ultrasonic P‐ and S‐wave velocities vary from about 2.8–3.8 km s?1 and 1.5–2.4 km s?1, respectively, over confining and differential pressures between 0.1 and 35 MPa, allowing derivation of associated dynamic elastic moduli. Anisotropy was investigated using oriented core plugs for electrical resistivity, elastic wave velocity and attenuation, permeability, and tracer injection tests. Corvio sandstone shows weak transverse isotropy (symmetry axis normal to bedding) of <10% for velocity and <20% for attenuation. 相似文献
17.
A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks 总被引：4，自引：0，他引：4
The laboratory ultrasonic pulse‐echo method was used to collect accurate P‐ and S‐wave velocity (±0.3%) and attenuation (±10%) data at differential pressures of 5–50 MPa on water‐saturated core samples of sandstone, limestone and siltstone that were cut parallel and perpendicular to the vertical borehole axis. The results, when expressed in terms of the P‐ and S‐wave velocity and attenuation anisotropy parameters for weakly transversely isotropic media (ɛ, γ, ɛQ, γQ) show complex variations with pressure and lithology. In general, attenuation anisotropy is stronger and more sensitive to pressure changes than velocity anisotropy, regardless of lithology. Anisotropy is greatest (over 20% for velocity, over 70% for attenuation) in rocks with visible clay/organic matter laminations in hand specimens. Pressure sensitivities are attributed to the opening of microcracks with decreasing pressure. Changes in magnitude of velocity and attenuation anisotropy with effective pressure show similar trends, although they can show different signs (positive or negative values of ɛ, ɛQ, γ, γQ). We conclude that attenuation anisotropy in particular could prove useful to seismic monitoring of reservoir pressure changes if frequency‐dependent effects can be quantified and modelled. 相似文献
18.
波动方程有限差分法是波场模拟的一个重要方法,为解决常规有限差分法存在着数值频散的问题,本文从具有垂直对称轴的三维横向各向同性(VTI)介质频率-空间域qP波动方程出发,在常规差分算子的基础上构造了适合三维VTI介质的频率空间域有限差分优化算子,然后利用最优化理论中的Gauss-Newton法求解了优化算子的系数,使差分方程的相速度与波动方程的相速度尽量吻合,从而在理论上使网格数值频散达到极小,精度对比分析及数值测试表明,有限差分优化算子具有较高的波场数值模拟精度,有效地压制了数值频散现象,为三维VTI介质频率一空间域qP波正演模拟研究提供了理论基础. 相似文献
19.
JOS M. CARCIONE 《Geophysical Prospecting》1992,40(7):761-783
When a seismic signal propagates through a finely layered medium, there is anisotropy if the wavelengths are long enough compared to the layer thicknesses. It is well known that in this situation, the medium is equivalent to a transversely isotropic material. In addition to anisotropy, the layers may show intrinsic anelastic behaviour. Under these circumstances, the layered medium exhibits Q anisotropy and anisotropic velocity dispersion. The present work investigates the anelastic effect in the long-wavelength approximation. Backus's theory and the standard linear solid rheology are used as models to obtain the directional properties of anelasticity corresponding to the quasi-compressional mode qP, the quasi-shear mode qSV, and the pure shear mode SH, respectively. The medium is described by a complex and frequency-dependent stiffness matrix. The complex and phase velocities for homogeneous viscoelastic waves are calculated from the Christoffel equation, while the wave-fronts (energy velocities) and quality factor surfaces are obtained from energy considerations by invoking Poynting's theorem. We consider two-constituent stationary layered media, and study the wave characteristics for different material compositions and proportions. Analyses on sequences of sandstone-limestone and shale-limestone with different degrees of anisotropy indicate that the quality factors of the shear modes are more anisotropic than the corresponding phase velocities, cusps of the qSV mode are more pronounced for low frequencies and midrange proportions, and in general, attenuation is higher in the direction perpendicular to layering or close to it, provided that the material with lower velocity is the more dissipative. A numerical simulation experiment verifies the attenuation properties of finely layered media through comparison of elastic and anelastic snapshots. 相似文献