共查询到20条相似文献，搜索用时 31 毫秒
1.
在长波长假设条件下,水平层状地层中发育一组垂直排列的裂缝构成了等效正交各向异性介质.各向异性参数与裂缝弱度参数的估算有助于非均匀各向异性介质的各向异性特征描述,而弹性逆散射理论是非均匀介质参数反演的有效途径.基于地震散射理论,我们首先推导了非均匀正交介质中纵波散射系数方程,并通过引入正交各向异性特征参数,提出了一种新颖的正交各向异性方位弹性阻抗参数化方法.为了提高反演的稳定性与横向连续性,我们发展了贝叶斯框架下的正交各向异性方位弹性阻抗反演方法,同时考虑了柯西稀疏约束正则化和平滑模型约束正则化,最终使用非线性的迭代重加权最小二乘策略实现了各向异性特征参数的稳定估算.模型和实际资料处理表明,反演结果与测井解释数据相吻合,证明了该方法能够稳定可靠地从方位叠前地震资料中获取各向异性特征参数,减小参数估算的不确定性,为非均匀正交介质的各向异性预测提供了一种高可靠性的地震反演方法. 相似文献
2.
长波长假设条件下,各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性（HTI）介质.基于不同观测方位的岩石地震响应特征变化,宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测,同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法,首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数,并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联,并采用马尔科夫链蒙特卡洛（MCMC）抽样方法进行大量样本的随机模拟,使用期望最大化（EM）算法估计后验条件概率分布,最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明,该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数,并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法. 相似文献
3.
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. 相似文献
4.
Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures. 相似文献
5.
Wide-azimuth seismic data can be used to derive anisotropic parameters on the subsurface by observing variation in subsurface seismic response along different azimuths. Layer-based high-resolution estimates of components of the subsurface anisotropic elastic tensor can be reconstructed by using wide-azimuth P-wave data by combining the kinematic information derived from anisotropic velocity analysis with dynamic information obtained from amplitude versus angle and azimuth analysis of wide-azimuth seismic data. Interval P-impedance, S-impedance and anisotropic parameters associated with anisotropic fracture media are being reconstructed using linearized analysis assuming horizontal transverse anisotropy symmetry. In this paper it is shown how additional assumptions, such as the rock model, can be used to reduce the degrees of freedom in the estimation problem and recover all five anisotropic parameters. Because the use of a rock model is needed, the derived elastic parameters are consistent with the rock model and are used to infer fractured rock properties using stochastic rock physics inversion. The inversion is based on stochastic rock physics modelling and maximum a posteriori estimate of both porosity and crack density parameters associated with the observed elastic parameters derived from both velocity and amplitude versus angle and azimuth analysis. While the focus of this study is on the use of P-wave reflection data, we also show how additional information such as shear wave splitting and/or anisotropic well log data can reduce the assumptions needed to derive elastic parameter and rock properties. 相似文献
6.
Based on the long-wavelength approximation, a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study, we propose a reliable method for predicting the Thomsen’s weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen’s weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen’s weak anisotropic parameters and fracture weaknesses. Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth (AVAZ) inversion for Thomsen’s weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares (IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen’s weak anisotropic parameters and fracture weaknesses. 相似文献
7.
地应力是评价页岩气储层是否可压裂成网的重要参数.利用地震资料进行地应力参数预测可以根据水平应力差异比(DHSR,Differential Horizontal Stress Ratio)值大小来判断.现阶段计算水平应力差异比主要基于具有对称轴的横向各向同性介质理论,即HTI介质假设,欠缺对页岩中实际大量存在的VTI特征的考虑.本研究同时考虑页岩气地层VTI和HTI特征的作用,利用正交各向异性介质岩石物理关系,推导出水平应力差异比与弹性参数及各向异性参数之间的关系.首先,在分析应力应变本构方程的基础上,对正交各向异性介质(OA)的本构方程进行推导,通过胡克定律一般形式获得正交各向异性介质的应力和应变关系;其次,利用该应力和应变关系推导出最大水平应力、最小水平应力和水平应力差异比与地层弹性参数和各向异性参数的关系并定义了参数正交各向异性水平应力差异比(ODHSR,Orthorhombic Differential Horizontal Stress Ratio);最后,对正交各向异性水平应力差异比ODHSR与具有水平对称轴的横向各向同性介质的DHSR进行了对比,证明了本文的正交各向异性水平应力差异比(ODHSR)计算公式. 相似文献
8.
Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media 
下载免费PDF全文
下载免费PDF全文 A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account. 相似文献
9.
本文从一般弱各向异性介质参数反演中得到了使用两条相互正交的变偏VSP测线计算弱各向异性(WA)参数的反演公式. 如果仅仅使用qP波, 则可以确定9个独立的WA参数, 这9个WA参数可以完全地描述井中接收点在两个剖面内介质的各向异性性质. 通过对走时曲线进行最小二乘条件下的三次样条光滑, 可以获得慢度矢量的垂直分量和水平分量. 如果介质是横向非均匀介质, 则水平分量仅仅被用作反演时的约束条件. 为了获得偏振矢量, 本文引入质心计算方法, 该方法计算简单、 稳定, 而且不存在解的奇异问题. 在寻找与一般各向异性介质最接近的高对称性正交各向异性介质和TTI各向异性介质时, 使用qP波各向异性坐标变换方法和最小二乘求解方法, 得到了与一般各向异性介质最接近的正交各向异性和TTI各向异性参数及其对称轴方向参数的计算公式. 使用这些方法, 对瓜哇海地区布设的两条相互正交的变偏VSP测线数据进行各向异性反演, 获得了井中10个接收点处介质的WA参数. 数值计算和实际资料反演表明, 本文所使用的反演方法能够准确地得到VSP井中接收点处介质的WA参数, 这是地震勘探中研究地壳介质各向异性性质最直接和最可靠的方法. 相似文献
10.
各向异性介质弹性波方程全波形反演过程中多参数之间的相互耦合,使得弱参数在反演过程中难得到理想的结果.本文以VTI介质为例,在各参数辐射模式分析的基础上,基于改进的散射积分算法实现目标函数梯度的直接求取,进一步构建高斯牛顿方向,实现Hessian矩阵的有效利用,以考虑Hessian矩阵非主对角线元素包含的各参数间的耦合效应,在不使用任何反演策略的情况下实现高精度的VTI介质弹性波方程多参数同步反演.同时,该方法在计算过程中无需存储庞大的核函数矩阵,且无需传统截断牛顿法中额外的正演计算,因此内存占用小,计算效率高.本文数值试验验证了该方法的有效性,为各向异性多参数全波形反演提供了一种新的解决方案. 相似文献
11.
正交各向异性介质P波走时分析及Thomsen参数反演 总被引：10，自引：3，他引：10
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层，正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时，时距曲线表现出强方位各向异性.采用小生境遗传算法，对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明，此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数. 相似文献
12.
— 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. 相似文献
13.
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. 相似文献
14.
Analysing S-wave splitting has become a routine step in processing multicomponent data. Typically, this analysis leads to determining the principal directions of a transversely isotropic medium with a horizontal symmetry axis, which is assumed to be responsible for azimuthal anisotropy, and to the time delays between the fast and slow S-waves. These parameters are commonly estimated layer-by-layer from the top. Errors in layer stripping occurring in shallow layers might propagate to deeper layers. We propose a method for S-wave splitting analysis and compensation that consists of inverting interval values of splitting intensity to obtain a model of anisotropic parameters that vary with time and/or depth. Splitting intensity is a robust attribute with respect to structural variations and is commutative, which means that it can be summed along a ray (or throughout a sensitivity kernel volume) and can be linearly related to anisotropic perturbations at depth. Therefore, it is possible to estimate anisotropic properties within a geological formation (e.g. the reservoir) by analysing the differences of splitting intensity measured at the top and at the bottom of the layer. This allows us to avoid layer stripping, in particular, for shallow layers where anisotropic parameters are difficult to estimate due to poor coverage, and it makes S-wave splitting analysis simpler to apply. We demonstrate this method on synthetic and real data. Because the splitting intensity attribute shows usefulness in S-wave splitting analysis in transversely isotropic media, we extend the splitting intensity theory to lower symmetry classes. It enables the characterization of tilted transversely isotropic and tilted orthorhombic media, opening new opportunities for anisotropic model building. 相似文献
15.
Amir Abbas Babasafari Deva Prasad Ghosh Ahmed Mohammed Ahmed Salim Masoumeh Kordi 《Geophysical Prospecting》2020,68(8):2471-2493
Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification. 相似文献
16.
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. 相似文献
17.
In multi-parameter ray-based anisotropic migration/inversion, it is essential that we have an understanding of the scattering mechanism corresponding to parameter perturbations. Because the complex nonlinearity in the anisotropic inversion problem is intractable, the construction of true-amplitude linearized migration/inversion procedures is needed and important. By using the acoustic medium assumption for transversely isotropic media with a vertical axis of symmetry and representing the anisotropy with P-wave normal moveout velocity, Thomsen parameter δ and anelliptic parameter η, we formalize the linearized inverse scattering problem for three-dimensional pseudo-acoustic equations. Deploying the single-scattering approximation and an elliptically anisotropic background introduces a new linear integral operator that connects the discontinuous perturbation parameters with the multi-shot/multi-offset P-wave scattered data. We further apply the high-frequency asymptotic Green's function and its derivatives to the integral operator, and then the scattering pattern of each perturbation parameter can be explicitly presented. By naturally establishing a connection to generalized Radon transform, the pseudo-inverse of the integral operator can be solved by the generalized Radon transform inversion. In consideration of the structure of this pseudo-inverse operator, the computational implementation is done pointwise by shooting a fan of rays from the target imaging area towards the acquisition system. Results from two-dimensional numerical tests show amplitude-preserving images with high quality. 相似文献
18.
Anisotropic elastic full‐waveform inversion of walkaway vertical seismic profiling data from the Arabian Gulf 下载免费PDF全文
下载免费PDF全文 John C Owusu Marwan Charara Scott Leaney Allan Campbell Shujaat Ali Igor Borodin Les Nutt Henry Menkiti 《Geophysical Prospecting》2016,64(1):38-53
Borehole seismic addresses the need for high‐resolution images and elastic parameters of the subsurface. Full‐waveform inversion of vertical seismic profile data is a promising technology with the potential to recover quantitative information about elastic properties of the medium. Full‐waveform inversion has the capability to process the entire wavefield and to address the wave propagation effects contained in the borehole data—multi‐component measurements; anisotropic effects; compressional and shear waves; and transmitted, converted, and reflected waves and multiples. Full‐waveform inversion, therefore, has the potential to provide a more accurate result compared with conventional processing methods. We present a feasibility study with results of the application of high‐frequency (up to 60 Hz) anisotropic elastic full‐waveform inversion to a walkaway vertical seismic profile data from the Arabian Gulf. Full‐waveform inversion has reproduced the majority of the wave events and recovered a geologically plausible layered model with physically meaningful values of the medium. 相似文献
19.
《中国科学:地球科学(英文版)》2015,(8)
The existing expressions of elastic impedance,as the generalized form of acoustic impedance,represent the resistance of subsurface media to seismic waves of non-normal incidence,and thus include information on the shear-wave velocity.In this sense,conventional elastic impedance is an attribute of the seismic reflection and not an intrinsic physical property of the subsurface media.The derivation of these expressions shares the approximations made for reflectivity,such as weak impedance contrast andisotropic or weakly anisotropic media,which limits the accuracy of reflectivity reconstruction and seismic inversion.In this paper,we derive exact elastic impedance tensors of seismic P-and S-waves for isotropic media based on the stress-velocity law.Each componentof the impedance tensor represents a unique mechanical property of the medium.Approximations of P-wave elastic impedance tensor components are discussed for seismic inversion and interpretation.Application to synthetic data and real data shows the accuracy and robust interpretation capability of the derived elastic impedance in lithology characterizations. 相似文献
20.
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. 相似文献