共查询到20条相似文献，搜索用时 15 毫秒
1.
We present a method for inversion of fracture compliance matrix components from wide‐azimuth noisy synthetic PS reflection data and quantitatively show that reflection amplitude variations with offset and azimuth for converted PS‐waves are more informative than P‐waves for fracture characterization. We consider monoclinic symmetry for fractured reservoir (parameters chosen from Woodford Shale), which can be formed by two or more sets of vertical fractures embedded in a vertically transverse isotropic background. Components of effective fracture compliance matrices for a medium with monoclinic symmetry are related to the characteristics of the fractured medium. Monte Carlo simulation results show that inversion of PS reflection data is more robust than that of PP reflection data to uncertainties in our a priori knowledge (vertically transverse isotropic parameters of unfractured rock) than PP reflection data. We also show that, while inversion of PP reflections is sensitive to contrasts in elastic properties of upper and lower media, inversion of PS reflections is robust with respect to such contrasts. 相似文献
2.
We use linear slip theory to evaluate seismic reflections at non‐welded interfaces, such as faults or fractures, sandwiched between general anisotropic media and show that at low frequencies the real parts of the reflection coefficients can be approximated by the responses of equivalent welded interfaces, whereas the imaginary parts can be related directly to the interface compliances. The imaginary parts of low frequency seismic reflection coefficients at fault zones can be used to estimate the interface compliances, which can be related to fault properties upon using a fault model. At normal incidence the expressions uncouple and the complex‐valued P‐wave reflection coefficient can be related linearly to the normal compliance. As the normal compliance is highly sensitive to the infill of the interface, it can be used for gas/fluid identification in the fault plane. Alternatively, the tangential compliance of a fault can be estimated from the complex‐valued S‐wave reflection coefficient. The tangential compliance can provide information on the crack density in a fault zone. Coupling compliances can be identified and quantified by the observation of PS conversion at normal incidence, with a comparable linear relationship. 相似文献
3.
When a porous layer is permeated by mesoscale fractures, wave-induced fluid flow between pores and fractures can cause significant attenuation and dispersion of velocities and anisotropy parameters in the seismic frequency band. This intrinsic dispersion due to fracturing can create frequency-dependent reflection coefficients in the layered medium. In this study, we derive the frequency-dependent PP and PS reflection coefficients versus incidence angle in the fractured medium. We consider a two-layer vertical transverse isotropy model constituted by an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce the intrinsic dispersion due to fracturing in the sand layer. Based on the series coefficients that control the behaviour of velocity and anisotropy parameters in the fractured medium at low frequencies, we extend the conventional amplitude-versus-offset equations into frequency domain and derive frequency-dependent amplitude-versus-offset equations at the elastic–anelastic surface. Increase in fracture length or fracture density can enlarge the frequency dependence of amplitude-versus-offset attributes of PP and PS waves. Also, the frequency dependence of magnitude and phase angle of PP and PS reflection coefficients increases as fracture length or fracture density increases. Amplitude-versus-offset type of PP and PS reflection varies with fracture parameters and frequency. What is more, fracture length shows little impact on the frequency-dependent critical phase angle, while the frequency dependence of the critical phase angle increases with fracture density. 相似文献
4.
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. 相似文献
5.
Extended reflectivity method for modelling the propagation of diffusive–viscous wave in dip‐layered media 
下载免费PDF全文
下载免费PDF全文 The reflectivity method plays an important role in seismic modelling. It has been used to model different types of waves propagating in elastic and anelastic media. The diffusive–viscous wave equation was proposed to investigate the relationship between frequency dependence of reflections and fluid saturation. It is also used to describe the attenuation property of seismic wave in a fluid‐saturated medium. The attenuation of diffusive–viscous wave is mainly characterised by the effective attenuation parameters in the equation. Thus, it is essential to obtain those parameters and further characterise the features of the diffusive–viscous wave. In this work, we use inversion method to obtain the effective attenuation parameters through quality factor to investigate the characteristics of diffusive–viscous wave by comparing with those of the viscoacoustic wave. Then, the reflection/transmission coefficients in a dip plane‐layered medium are studied through coordinate transform and plane‐wave theory. Consequently, the reflectivity method is extended to compute seismograms of diffusive–viscous wave in a dip plane multi‐layered medium. Finally, we present two models to simulate the propagation of diffusive–viscous wave in a dip plane multi‐layered medium by comparing the results with those in a viscoacoustic medium. The numerical results demonstrate the validity of our extension of reflectivity method to the diffusive–viscous medium. The numerical examples in both time domain and time–frequency domain show that the reflections from a dip plane interface have significant phase shift and amplitude change compared with the results of horizontal plane interface due to the differences in reflection/transmission coefficients. Moreover, the modelling results show strong attenuation and phase shift in the diffusive–viscous wave compared to those of the viscoacoustic wave. 相似文献
6.
We recently proposed an efficient hybrid scheme to absorb boundary reflections for acoustic wave modelling that could attain nearly perfect absorptions. This scheme uses weighted averaging of wavefields in a transition area, between the inner area and the model boundaries. In this paper we report on the extension of this scheme to 2D elastic wave modelling with displacement‐stress formulations on staggered grids using explicit finite‐difference, pseudo‐implicit finite‐difference and pseudo‐spectral methods. Numerical modelling results of elastic wave equations with hybrid absorbing boundary conditions show great improvement for modelling stability and significant absorption for boundary reflections, compared with the conventional Higdon absorbing boundary conditions, demonstrating the effectiveness of this scheme for elastic wave modelling. The modelling results also show that the hybrid scheme works well in 2D rotated staggered‐grid modelling for isotropic medium, 2D staggered‐grid modelling for vertically transversely isotropic medium and 2D rotated staggered‐grid modelling for tilted transversely isotropic medium. 相似文献
7.
This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results. 相似文献
8.
R.-G. FERBER 《Geophysical Prospecting》1985,33(2):212-223
A main problem in computing reflection coefficients from seismograms is the instability of the inversion procedure due to noise. This problem is attacked for two well-known inversion schemes for normal-incidence reflection seismograms. The crustal model consists of a stack of elastic, laterally homogeneous layers between two elastic half-spaces. The first method, which directly computes the reflection coefficients from the seismogram is called “Dynamic Deconvolution”. The second method, here called “Inversion Filtering”, is a two-stage procedure. The first stage is the construction of a causal filter by factorization of the spectral function via Levinson-recursion. Filtering the seismogram is the second stage. The filtered seismogram is a good approximation for the reflection coefficients sequence (unless the coefficients are too large). In the non-linear terms of dynamic deconvolution and Levinson-recursion the noise could play havoc with the computation. In order to stabilize the algorithms, the bias of these terms is estimated and removed. Additionally incorporated is a statistical test for the reflection coefficients in dynamic deconvolution and the partial correlation coefficients in Levinson-recursion, which are set to zero if they are not significantly different from noise. The result of stabilization is demonstrated on synthetic seismograms. For unit spike source pulse and white noise, dynamic deconvolution outperforms inversion filtering due to its exact nature and lesser computational burden. On the other hand, especially in the more realistic bandlimited case, inversion filtering has the great advantage that the second stage acts linearly on the seismogram, which allows the calculation of the effect of the inversion procedure on the wavelet shape and the noise spectrum. 相似文献
9.
传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的. 相似文献
10.
11.
The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case. 相似文献
12.
利用横波响应进行裂缝性各向异性介质的检测在实际应用中取得了很好的效果,但技术复杂、成本较高使该方法的广泛使用受到限制,而纵波资料采集和处理技术的精细有效保持了纵波的各种属性,这为直接利用纵波资料进行裂缝检测创造了条件.地下垂直定向裂缝通常用HTI介质模型来描述,为此,本文利用射线追踪和反射率法计算了层状各向同性介质背景下的HTI介质顶、底界面的反射纵波旅行时和反射系数,并分析了这些属性随观测方位的变化规律.研究表明,HTI介质底界面反射纵波旅行时和HTI介质顶界面反射系数表现出了明显的方位各向异性;旅行时、振幅和AVO梯度属性均在0°观测方位和90°观测方位上存在最大差异,可以用多种属性联合来精确判定裂缝的发育方向. 相似文献
13.
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. 相似文献
14.
In this paper we propose a method for the characterization of naturally fractured reservoirs by quantitative integration of seismic and production data. The method is based on a consistent theoretical frame work to model both effective hydraulic and elastic properties of fractured porous media and a (non‐linear) Bayesian method of inversion that provides information about uncertainties as well as mean (or maximum likelihood) values. We model a fractured reservoir as a porous medium containing a single set of vertical fractures characterized by an unknown fracture density, azimuthal orientation and aperture. We then look at the problem of fracture parameter estimation as a non‐linear inverse problem and try to estimate the unknown fracture parameters by joint inversion of seismic amplitude versus angle and azimuth data and dynamic production data. Once the fracture parameters have been estimated the corresponding effective stiffness and permeability tensors can be estimated using consistent models. A synthetic example is provided to clearly explain and test the workflow. It shows that seismic and production data complement each other, in the sense that the seismic data resolve a non‐uniqueness in the fracture orientation and the production data help to recover the true fracture aperture and permeability, because production data are more sensitive to the fracture aperture than the seismic data. 相似文献
15.
弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础. 相似文献
16.
本文基于Biot理论,推导出Zoeppritz形式的双相介质分界面上弹性波的反射与透射公式,对单相Zoeppritz公式与双相反射系数公式进行了比较,对双相介质含油、水、气不同流体时的反射规律以及孔隙度、渗透率、饱和度等储层参数对纵波反射的影响进行了研究.数值模拟分析表明,双相与单相反射公式的主要差异在于双相介质反射公式中考虑了液相、固液耦合相弹性模量的影响;油砂、水砂、气砂岩的慢P波反射差异明显;快P波反射对孔隙度的变化敏感,饱和度次之,对渗透率和频率的变化不敏感. 相似文献
17.
油气勘探实践表明, 裂缝是油气的储存空间或运移通道, 裂缝介质地震波场的研究越来越受到关注. 实际地层中裂缝形成受多种因素控制, 物理属性比较复杂, 表现出强烈的各向异性. 由于上覆地层的压力使得水平或低陡倾角裂缝存在较少, 大多数为高陡倾角裂缝, 利用线性滑动裂缝介质的等效理论将高陡倾角裂缝介质等效为横向各向同性介质, 便于实际应用. 本文采用各向异性弹性波旋转交错网格模拟方法对含裂缝介质单炮记录进行模拟与分析. 结果表明: 裂缝的存在相当于人为增加反射界面; 裂缝密度越大, 裂缝纵横比越小, 裂缝充填物与背景介质弹性性质差别越大, 引起的反射波能量变化越大. 本文模拟结果为利用地震数据进行裂缝介质参数反演与储层识别及油气预测提供了依据. 相似文献
18.
Average elastic properties of a fluid‐saturated fractured rock are discussed in association with the extremely slow and dispersive Krauklis wave propagation within individual fractures. The presence of the Krauklis wave increases P‐wave velocity dispersion and attenuation with decreasing frequency. Different laws (exponential, power, fractal, and gamma laws) of distribution of the fracture length within the rock show more velocity dispersion and attenuation of the P‐wave for greater fracture density, particularly at low seismic frequencies. The results exhibit a remarkable difference in the P‐wave reflection coefficient for frequency and angular dependency from the fractured layer in comparison with the homogeneous layer. The biggest variation in behaviour of the reflection coefficient versus incident angle is observed at low seismic frequencies. The proposed approach and results of calculations allow an interpretation of abnormal velocity dispersion, high attenuation, and special behaviour of reflection coefficients versus frequency and angle of incidence as the indicators of fractures. 相似文献
19.
We propose a novel computational method for the efficient simulation of two-phase flow in fractured porous media. Instead of refining the grid to capture the flow along the faults or fractures, we represent the latter as immersed interfaces, using a reduced model for the flow and suitable coupling conditions. We allow for non matching grids between the porous matrix and the fractures to increase the flexibility of the method in realistic cases. We employ the extended finite element method for the Darcy problem and a finite volume method that is able to handle cut cells and matrix-fracture interactions for the saturation equation. Moreover, we address through numerical experiments the problem of the choice of a suitable numerical flux in the case of a discontinuous flux function at the interface between the fracture and the porous matrix. A wrong approximate solution of the Riemann problem can yield unphysical solutions even in simple cases. 相似文献
20.
本文在各向同性介质假设下,计算了薄层顶、底反射P波和PS波的理论地震图.理论模拟发现薄层的反射P波与PS波是一复合波,包括层内的多次透射和反射转换波型,且具有类似单一界面反射的、脉冲式的波形特征.两类薄层反射整体振幅随着薄层厚度的降低而缩小;在地震子波主频40 Hz条件下,2 m左右厚度的薄层反射相比单一界面反射具有等同的振幅水平.1 m以下极薄层仅有弱反射甚至无反射;薄层反射复合波振幅与炮检距的关系依然成立,但不惟一;单界面反射AVO原理、方法不适用于薄层反射解释与反演. 相似文献