共查询到18条相似文献，搜索用时 109 毫秒
1.
P波通过裂缝地层传播时,其响应会随着地面观测方位的变化而变化,表现出明显的方位各向异性.这种随方位变化的裂缝响应信息(如振幅等)广泛存在于叠前方位一偏移距二维空间域中,检测此域中的裂缝响应,通过椭圆拟合的方法识别各向异性的存在,进而计算裂缝的方向——椭圆长轴的方向,和裂缝密度——椭圆长轴与短轴的比.本文研究出了一种基于常规叠前地震反射振幅来检测裂缝的三维叠前裂缝储层预测方法.它采用面元组合的方法扩大方位角和偏移距范围;采用多级数据控制技术获取所需的信息;采用稳妥的技术步骤来预测.通过大港油田实际工区的应用,结果表明本文研究出的方法技术应用效果好、精度高、经济实用. 相似文献
2.
《地球物理学进展》2017,(3)
叠前地震裂缝检测技术是目前裂缝型储层定量描述的主要手段之一,利用地质、测井资料和岩石物理参数,根据岩石物理理论建立裂缝储层的理论地质模型,模拟地震波在已知裂缝模型中传播时各个方位的地震响应,分析某区段由裂缝引起的地震振幅属性的方位AVO(AVAZ)特征,利用该技术流程,编写了一套基于HTI介质的正演模拟软件,在软件开发过程中,研发了基于Qt Concurrent编程模型组件的并行处理技术.通过该项技术的研究与开发可以:1)利用地震属性数据实时分析地下裂缝型储层的地震响应规律;2)解决地震振幅数据方位各向异性特征的不确定性导致的各向异性椭圆解释成果错误问题,确定椭圆的正确方位指示;3)为叠前裂缝检测软件提供正演分析工具,完善并扩充检测流程.在实际生产应用的效果表明,正演的解释结论能够为叠前裂缝检测提供有效的理论依据,增强了反演裂缝方位数据的准确性与可信度. 相似文献
3.
穿透含裂隙、裂缝地壳8s视周期的SV波的理论地震图研究表明,当地壳平均裂隙密度高于0.01即横波各向异性高于1%时,非对称面内不同方位的SKS波均发生分裂;地震图中直接的记录显示是切向T分量上出现SKS波的振动,其振幅随地壳平均裂隙密度的增大而增强,甚至能与径向R分量上的振幅相当.局限于上地壳的强裂缝各向异性同样能引起SKS分裂.长周期SKS波分裂对地壳内裂隙、裂缝的分布缺乏分辨率.直立平行排列裂隙、裂缝使得SKS分裂T分量记录特征具有方位对称性,这来自于HTI介质中快、慢波偏振和到时差随方位变化的对称性;而倾斜裂隙、裂缝使得该方位对称性丧失.对实际观测SKS分裂的偏振解释需要考虑地壳裂隙各向异性,特别是断裂附近的强裂缝各向异性. 相似文献
4.
为了实现三维倾斜界面真正的共反射点叠加,获得更加精确的P波速度模型,进行可靠的AVO/AVA分析,形成与三维多分量PS转换波共转换点道集配套的技术,在常速介质共反射点精确走时的基础上,推导并给出三维倾斜界面P波共反射点道集近似双曲时距关系及NMO速度;认识到近似时距在偏移距与深度比小于2时具有较高的拟合精度;指出P波共反射点道集NMO速度是P波速度、界面倾角、界面倾向和测线方位的函数,且随测线方位变化而呈现椭圆特征.与P波共中心点道集NMO速度不同.P波共反射点道集NMO速度小于介质速度.且随界面倾角增大而减小;椭圆长轴方向为倾斜界面走向方位,短轴方向为倾斜界面倾向方位. 相似文献
5.
利用地震数据较为准确地预测地下介质的裂缝发育信息,是裂缝型储层预测的关键手段,也是非常规页岩油气储层压裂改造的关键参数.方位AVO(AVAZ)在描述HTI介质的裂缝分布及发育方向等方面有重要应用,本文提供一种基于全方位地震数据的方位AVO(AVAZ)反演方法,为裂缝预测提供参考资料.首先利用Ruger近似方程正演计算不同方位角及入射角下的反射系数,继而与方位角度道集中地震反射振幅建立目标函数,然后基于正演方程及目标函数采用最小二乘法进行AVAZ反演.实现了从全方位地震资料中同步反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向.模型数据和实际地震资料应用表明该方法能够有效的预测HTI介质的裂缝分布及方向. 相似文献
6.
地震波场的方位特征对于裂缝性油气藏的地震预测方法研究有着重要的意义.为使裂缝模型更具一般性,本文基于线性滑动理论和Bond变换构造了两组任意夹角竖直裂缝模型;考虑到波场二维模拟中突出方位特征,依据弹性动力学的基本方程和Bond变换,推导出了含方位角的弹性波传播速度应力方程;使用高精度交错网格有限差分法对几种裂缝介质进行了数值模拟,并分析其波场特征的变化.模拟结果表明,在两组参数相同的任意夹角竖直裂缝模型中,裂缝内夹角各个方位的波场特征变化不明显,各向异性强度较低,裂缝外夹角各个方位的波场特征变化剧烈,裂缝正交时,各处波场特征相近;针对单组裂缝模型,分析了裂缝参数变化对波场特征的影响,总结出了三种波场形态,得出了裂缝垂直面的波场特征与裂缝法向弱度和切向弱度的相对大小有关,从裂缝垂直面到裂缝平行面的波场朝着相同的形态变化,与裂缝参数无关等认识.这些分析结果有助于进一步认识和应用裂缝介质的波场方位特征. 相似文献
7.
针对裂隙型储集层中更具代表性的各向异性介质模型,即在各向同性背景介质中含有两组斜交的垂直裂隙所构成的单斜各向异性介质模型,利用时间和空间上可达任意阶的高阶交错网格有限差分技术,对具有不同裂隙填充物性质的单斜介质中波的传播快照进行了模拟.结果证实各向异性介质中波的传播速度随传播方向的不同而产生明显的差异;裂隙填充物的性质对于速度各向异性具有很大的影响.另外,利用坐标旋转法,对水平层状各向异性介质中多方位地面三分量记录进行了模拟,结果表明了方位各向异性介质中,波的传播速度不仅随入射角的变化而变化,同时也随观测方位的不同而产生差异.数值模拟结果为进一步利用地面多方位地震属性进行各向异性参数的反演及裂隙参数的描述提供理论基础. 相似文献
8.
在含有定向排列裂缝组的多孔岩石地层中,特别是碳酸盐岩储集层,弹性波的传播和衰减,因受诸如裂缝闭合、尺寸、及其填充的流体和走向等因素影响,具有各向异性特征。在本次研究中,我们采用了基于Chapman等效介质模型,利用数值实验方法研究了这类储层中P波传播特性诸如速度、衰减和横波分裂等各向异性随频率和方位的变化规律。模拟结果表明,当一组裂缝闭合与张开时随方位变化的速度、衰减和各向异性是有所不同的。随方位变化衰减最小值或P波速度最大值与张开裂缝的走向趋于一致,闭合裂缝的P波速度大于张开裂缝,而衰减和各向异性则张开裂缝的大于闭合裂缝;不同尺度裂缝,速度的最大值和衰减最小值与裂缝组平均方位对应,小尺度裂缝对波的传播有小的影响,方位依赖的各向异性相比其他裂缝属性有更小的影响;裂缝密度对P波速度、衰减和各向异性有更大的影响,而衰减比速度和各向异性更敏感裂缝尺寸;在地震勘探频段油、气饱和的衰减与盐水饱和不同,充填油和气随方位变化值比较低。两组裂缝有相同的密度,快横波偏振方位线性的决定于一组裂缝的方位。 相似文献
9.
裂缝性储层是乍得Bongor盆地潜山油藏主要储集类型，研究区地震资料较差，基底成像不清，利用地震资料进行裂缝预测符合率不高.近几年，在乍得多个潜山实施了"宽频、宽方位、高密度"（以下称"两宽一高"）三维地震勘探，基底成像效果有了明显改善.由于"两宽一高"地震勘探近几年刚刚兴起，目前利用"两宽一高"OVT处理资料预测裂缝的方法较少，主要是基于振幅随偏移距及方位角的变化的原理（AVAZ），目前主要的方法为基于椭圆拟合的各向异性强度预测法，本文首先应用椭圆拟合法预测了Bongor盆地P潜山的各向异性强度，对结果进行了分析，指出了椭圆拟合法的局限性，进而首次在Bongor盆地研发并应用了基于统计法的各向异性强度预测技术，即首先通过道集规则化将不规则的道集处理成螺旋道集，然后在道集中统计振幅或双程走时的方差进行裂缝预测，在乍得P潜山取得了较好的效果，预测符合率达到80%以上，平面规律性明显. 相似文献
10.
《应用地球物理》2015,(3)
基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体,其地震反射波形序列包含如下动力学信息,即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减,以及来自顶底界面的反射波的调谐与干涉等。计算表明,速度频散和衰减增顶界面反射波大入射角反射时的振幅,而减弱底界面反射振幅。对于固定入射角的方位地震方位地震响应,PP波反射特征表现为随方位角的增加反射波形序列延续时间变长,而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小,表现为PSV波反射振幅随方位角增加而增加,PSH波在0°和90°方位无反射能量,在45°方位反射振幅最强。 相似文献
11.
舒梦珵 《CT理论与应用研究》2017,26(1):45-52
对于裂缝型油气藏，基于叠前方位角地震数据体，利用AVAZ方法进行裂缝属性预测是目前较为有效的裂缝描述手段，然而，由于地震数据通常为窄方位采集，使得不同方位角数据覆盖次数差异较大，导致各方位振幅能量分布不均，最终造成裂缝属性预测中产生趋势性误差。振幅归一化处理技术应用到现有的AVAZ方法中，该方法能够有效地降低因数据采集所造成的趋势性误差，预测结果会更为真实、可靠。在实际油田的应用中，振幅归一化技术显著改善了AVAZ方法对裂缝属性预测的准确性，取得了较好的效果。 相似文献
12.
13.
Weak‐anisotropy approximation for P‐wave reflection coefficient at the boundary between two tilted transversely isotropic media 
下载免费PDF全文
下载免费PDF全文 Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy. 相似文献
14.
Understanding fracture orientations is important for optimal field development of fractured reservoirs because fractures can act as conduits for fluid flow. This is especially true for unconventional reservoirs (e.g., tight gas sands and shale gas). Using walkaround Vertical Seismic Profiling (VSP) technology presents a unique opportunity to identify seismic azimuthal anisotropy for use in mapping potential fracture zones and their orientation around a borehole. Saudi Aramco recently completed the acquisition, processing and analysis of a walkaround VSP survey through an unconventional tight gas sand reservoir to help characterize fractures. In this paper, we present the results of the seismic azimuthal anisotropy analysis using seismic traveltime, shear‐wave splitting and amplitude attenuation. The azimuthal anisotropy results are compared to the fracture orientations derived from dipole sonic and image logs. The image log interpretation suggests that an orthorhombic fracture system is present. VSP data show that the P‐wave traveltime anisotropy direction is NE to SW. This is consistent with the cemented fractures from the image log interpretation. The seismic amplitude attenuation anisotropy direction is NW to SE. This is consistent with one of the two orientations obtained using transverse to radial amplitude ratio analysis, with the dipole sonic and with open fracture directions interpreted from image log data. 相似文献
15.
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP- and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities. 相似文献
16.
Colin M. Sayers 《Geophysical Prospecting》2009,57(2):187-192
Natural fractures in reservoirs play an important role in determining fluid flow during production and knowledge of the orientation and density of fractures is required to optimize production. Variations in reflection amplitude versus offset (AVO) are sensitive to the presence of fractures but current models used to invert the seismic response often make simplified assumptions that prevent fractured reservoirs from being characterized correctly. For example, many models assume a single set of perfectly aligned fractures, whereas most reservoirs contain several fracture sets with variable orientation within a given fracture set. In addition, many authors only consider the azimuthal variation in the small offset amplitude versus offset and azimuth response (the variation in AVO gradient with azimuth), while the effect of fractures on amplitude versus offset and azimuth increases with increasing offset. In this paper, the variation in the reflection coefficient of seismic P -waves as a function of azimuth and offset due to the presence of multiple sets of fractures with variable orientation within any fracture set is used to determine the components of a second-rank fracture compliance tensor α ij . The variation in the trace of this tensor as a function of position in the reservoir can be used to estimate the variation in fracture density with position in the reservoir and may be used to choose the location of infill wells in the field. The principal axes of α ij reveal the most compliant direction within the reservoir and may be used to optimize the trajectory of deviated wells. The determination of the principal axes of α ij requires wide azimuth acquisition and the use of the small-offset amplitude versus offset and azimuth (the azimuthal variation of the AVO gradient) may give misleading results. 相似文献
17.
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. 相似文献
18.
Rayleigh wave dispersion can be induced in an anisotropic medium or a layered isotropic medium. For a layered azimuthally anisotropic structure, traditional wave equation of layered structure can be modified to describe the dispersion behavior of Rayleigh waves. Numerical stimulation results show that for layered azimuthal anisotropy both the dispersion velocities and anisotropic parameters depend principally on anisotropic S-wave velocities. The splitting S-wave velocities may produce dispersion splitting of Rayleigh waves. Such dispersion splitting appears noticeable at azimuthal angle 45°. This feature was confirmed by the measured results of a field test. The fundamental mode splits into two branches at azimuthal angle 45° to the symmetry axis for some frequencies, and along the same direction the difference of splitting-phase velocities of the fundamental model reaches the maximum. Dispersion splitting of Rayleigh waves was firstly displayed for anisotropy study in dispersion image by means of multichannel analysis of surface waves, the image of which provides a new window for studying the anisotropic property of media. 相似文献