首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.  相似文献   

2.
Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth.  相似文献   

3.
为研究裂缝、裂隙介质中波致流引起的衰减,将裂缝看作背景孔隙岩石中非常薄且孔隙度非常高的层状介质,并等价成White周期层状模型.分别考虑不同类型的裂隙和孔隙之间的挤喷流影响,结合改进的Biot方程,推导得到裂缝裂隙介质的刚度与频率的关系.当缝隙中饱含流体时,介质的衰减和速度频散受裂缝、孔隙之间和裂隙、孔隙之间流体流动的显著影响.在低频极限下,裂缝裂隙介质的性质由各向异性Gassmann理论和挤喷流模型获得;而在非常高的频率时,由于缝隙中的压力来不及达到平衡,波致流的影响可忽略.分析表明,裂隙密度主要影响波的衰减,而裂隙纵横比主要控制优势衰减频率和速度显著变化的频率范围;由于不同裂隙的衰减机制不同,衰减和速度频散大小有所差异,但基本趋势相同.  相似文献   

4.
虑传播效应的多波保幅AVO正演(英文)   总被引:1,自引:1,他引:0  
传统的AVO正演只考虑了单一界面的反射系数对地震波波场振幅的影响,忽略了地震波在介质中传播的各种传播效应。通过引入地震波在介质中传播的几何扩散、吸收衰减以及透射损失等传播效应,提出了基于射线理论的水平层状介质多波保幅AVO正演方法。推导了水平层状介质多波几何扩散校正公式,来描述多波在介质中传播的几何扩散效应。通过直接引入复旅行时,而无需借助复速度,建立了复旅行时与品质因子的关系,来描述粘弹介质的吸收衰减。直接求解Zoeppritz方程计算多波的透射系数,用于描述多波在介质中传播时的透射损失。数值计算表明,几何扩散、吸收衰减以及透射损失对多波振幅的影响是随偏移距变化而变化的,多波保幅AVO正演需要考虑波传播效应对反射波振幅的改造。  相似文献   

5.
Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.  相似文献   

6.
An analysis of amplitudes of refraction records of some shallow refraction profiles shot primarily for detailing the near-surface structure in a granitic terrain has yielded information on refractor properties: reduced amplitudes are plotted on amplitude-distance graphs. The negative power n to which distance should be raised to represent (elastic) amplitude decay with respect to distance due to spreading of the critically refracted wave involved is examined. Computed values of this “spreading index”n are close to n = 2 as predicted by the theory. With this value of n, amplitude data are processed to determine residual attenuation attributable to elastic absorption in the bedrock. A graphical approach for this purpose from comparison of amplitude-distance graphs with the plots of amplitude decay due to spreading which is applicable to flat and horizontal refractor situations is suggested. Assuming residual attenuation to represent absorption in the granite bedrock, the computed coefficients of absorption, which vary from 0.5 to 3.90 km?1 for a frequency of 50 Hz, are obtained. From amplitude graphs of reversed profiles it is shown that the amplitude differences plot bears a relation to lateral velocity changes in the refractor. From comparison of practical amplitude decay graphs with those computed for different subsurface models, it appears possible to detect fractured rock occurrences in the refractor.  相似文献   

7.
An effective medium model for the stress-dependent seismic properties of fractured reservoirs is developed here on the basis of a combination of a general theory of viscoelastic waves in rock-like composites with recently published formulae for deformation of communicating and interacting cavities (interconnected pores/cracks and fractures at finite concentration) under drained loading. The inclusion-based model operates with spheroidal cavities at two different length scales; namely, the microscopic scale of the pores and (grain-boundary) cracks, and the mesoscopic scale of the fractures (controlling the flow of fluid). The different cavity types can in principle have any orientation and aspect ratio, but the microscopic pores/cracks and mesoscopic fractures were here assumed to be randomly and vertically oriented, respectively. By using three different aspect ratios for the relatively round pores (representing the stiff part of the pore space) and a distribution of aspect ratios for the relatively flat cracks (representing the compliant part of the pore space), we obtained a good fit between theoretical predictions and ultrasonic laboratory measurements on an unfractured rock sample under dry conditions. By using a single aspect ratio for the mesoscopic fractures, we arrived at a higher-order microstructural model of fractured porous media which represents a generalization of the first-order model developed by Chapman et al. (2002,2003). The effect of cavity size was here modelled under the assumption that the characteristic time for wave-induced (squirt) flow at the scale of a particular cavity (pore/crack vs. fracture) is proportional with the relevant scale-size. In the modelling, we investigate the effect of a decreasing pore pressure with constant confining pressure (fixed depth), and hence, increasing effective pressure. The analysis shows that the attenuation-peak due to the mesoscopic fractures in the reservoir will move downward in frequency as the effective pressure increases. In the range of seismic frequencies, our modelling indicates that the P-wave velocities may change by more than 20% perpendicular to the fractures and close to 10% parallel to the fractures. In comparison, the vertical S-wave velocities change by only about 5% for both polarization directions (perpendicular and parallel to the fractures) when the effective pressure increases from 0 to 15 MPa. This change is mainly due to the overall change in porosity with pressure. The weak pressure dependence is a consequence of the fact that the S waves will only sense if the fractures are open or not, and since all the fractures have the same aspect ratio, they will close at the same effective pressure (which is outside the analysed interval). Approximate reflection coefficients were computed for a model consisting of the fractured reservoir embedded as a layer in an isotropic shale and analysed with respect to variations in Amplitude Versus Offset and aZimuth (AVOZ) properties at seismic frequencies for increasing effective pressure. For the P-P reflections at the top of the reservoir, it is found that there is a significant dependence on effective pressure, but that the variations with azimuth and offset are small. The lack of azimuthal dependence may be explained from the approximate reflection coefficient formula as a result of cancellation of terms related to the S-wave velocity and the Thomson’s anisotropy parameter δ. For the P-S reflection, the azimuthal dependence is larger, but the pressure dependence is weaker (due to a single aspect ratio for the fractures). Finally, using the effective stiffness tensor for the fractured reservoir model with a visco-elastic finite-difference code, synthetic seismograms and hodograms were computed. From the seismograms, attenuation changes in the P wave reflected at the bottom of the reservoir can be observed as the effective pressure increases. S waves are not much affected by the fractures with respect to attenuation, but azimuthal dependence is stronger than for P waves, and S-wave splitting in the bottom reservoir P-S reflection is clearly seen both in the seismograms and hodograms. From the hodograms, some variation in the P-S reflection with effective pressure can also be observed.  相似文献   

8.
含裂隙介质中的视电阻率各向异性变化   总被引:1,自引:0,他引:1       下载免费PDF全文
解滔  卢军 《地球物理学报》2020,63(4):1675-1694
我国50多年的视电阻率连续观测结果表明,大地震前近震中区域的视电阻率呈现出与主压应力方位有关的各向异性变化,即:垂直于主压应力方向观测的变化幅度最大,平行方向最小或不明显,斜交方向介于二者之间.目前我国定点台站视电阻率观测的探测范围主要在浅层沉积层以内,通常含有较多的含水裂隙.本文将地下岩土介质简化为由固体基质和含流体/气体裂隙组成的固液气三相介质,且基质、流体和气体具有标量形式的电阻率,推导出了包含基质和流体电阻率、裂隙率、饱和度和裂隙面积率因子的电阻率张量表达式.以裂隙的扩展/闭合表示应力作用下裂隙的变化,得到了电阻率随裂隙变化的微分形式,电阻率变化对裂隙体积变化放大系数的表达式和裂隙横向变化对纵向电阻率影响的横向权系数的表达式.在此基础上得到了介质电阻率和视电阻率的各向异性变化特征:对于含水裂隙介质,无论裂隙如何变化,均是最小主轴方向电阻率的变化幅度大于其他方向;对于含水孔隙介质,沿孔隙主要变化方向的主轴电阻率变化幅度大于其他方向.对于各向异性变化,视电阻率和介质电阻率存在π/2的方向差异.相较于含水岩石,无水岩石介质电阻率的各向异性变化不显著.本文提出的电阻率表达式可以对实验室和野外实际观测的许多结果做出合理的解释.  相似文献   

9.
根据Chapman理论模型,在各向异性介质(如HTI介质)中,当入射角在0-45。范围内,慢横波会发生较大的衰减和频散,且对流体粘度敏感,而P波和快横波则比较小。对于沿裂隙法向传播的慢横波,其振幅受流体影响很大。因此,在P波响应对流体不敏感的情况下,可利用慢横波来获得裂隙型油气藏的流体信息。本文分析了胜利油田垦71地区三维三分量地震数据,检测出的慢横波振幅和旅行时异常与该区的测井资料十分吻合。分析结果还发现,与含油区相比,含水区会产生更高的横波分裂。在含水区,慢横波振幅会产生明显变化,而在含油区则几乎没有变化。  相似文献   

10.
Recent observations show that the differential amplitudes between the faster and slower split shear-waves in reflection surveys contain information about lateral variations of crack density in cracked reservoirs. However, the variation of amplitude with crack geometry when the crack strike changes with depth has not been reported previously. In this paper, we derive expressions for reflection and transmission coefficients of plane split shear-waves at vertical incidence at an interface separating two cracked (anisotropic) media with different crack strikes. We examine the effects on these coefficients as crack strike and crack density vary. For interfaces with large velocity-contrasts, the reflection coefficients carry little information about crack geometry, and the effects of crack strike varying with depth are negligible. In such cases, the polarization and time-delay of the split shear-waves are the only features which reliably diagnose anisotropy and contain information about the variation of crack strike and density. However, for interfaces with small velocity-contrasts, the effects of any variation of crack strike with depth cannot be neglected. In such cases, in addition to the polarization and time-delays of split shear-waves, both the differential amplitude of faster and slower shear-waves and the amplitude ratio of the two off-diagonal elements in the reflected data matrix after separation of split shear-waves, contain information about the variation of crack strike and crack density. In contrast, effects of crack strike changing with depth on transmitted waves are more sensitive regardless of the velocity-contrast and the degree of anisotropy.  相似文献   

11.
The bottom simulating reflector (BSR) in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional (fast) and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional (slow) waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.  相似文献   

12.
弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础.  相似文献   

13.
地球物理勘探中,地下应力的估测是非常重要的.本文利用应力诱导的弱各向异性的各向异性参数,借助已有的具有水平对称轴横向各向同性介质顾及高阶小量的横波反射系数近似公式,针对下介质为各向同性受单轴应力作用情况随应力的增加求取应力精度降低的问题,从横波反射系数近似公式的推导出发,联合SH波与SV波反射系数给出新的改进的单轴应力与横波反射系数的直观表达式;计算并对比了受单轴应力作用的各向同性介质的横波近似和严格的反射系数,并利用理论得到的横波反射系数对应力进行了估值.结果表明,利用改进后横波反射系数与应力的直观表达式,当应力增加时,应力估值精度有了明显的改善,这为利用横波反射数据直观估测应力提供理论依据.  相似文献   

14.
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.  相似文献   

15.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大.  相似文献   

16.
高角度缝隙充填的碳酸盐岩储层可以等效为具有水平对称轴的横向各向同性介质.本文提出了适用于裂缝型碳酸盐岩的岩石物理模型构建流程,重点介绍了在碳酸盐岩各向同性背景中,综合利用微小裂隙模型和线性滑动模型添加缝隙系统,并分析了当缝隙充填不同流体时,各向异性参数随纵横波速比的变化特征.同时本文讨论了裂缝密度和缝隙充填流体对地震反射系数的影响,推导了不同类型流体充填时储层反射系数与裂缝密度的近似关系式,阐述了各向异性流体替换理论,最终实现饱含流体碳酸盐岩裂缝储层的纵横波速度和各向异性参数的估测.选取某碳酸盐岩工区A井对该方法进行试算,结果表明基于碳酸盐岩裂缝岩石物理模型估算的纵横波速度值与测井值吻合较好,而且估测所得的各向异性参数值也能够较好地反映出裂缝储层位置.  相似文献   

17.
裂缝型储层流体识别方法   总被引:2,自引:0,他引:2       下载免费PDF全文
裂缝型储层的描述包括预测裂缝分布特征和识别裂隙充填物.依据等效介质理论计算的纵波速度随裂缝密度的增大而减小.正演地震记录显示,裂缝介质含气时反射振幅最大,且变化程度比含油或含水时大.叠前方位AVO反演所得的各向异性梯度Bani与裂缝密度成正比,可用于描述有效裂缝发育强度.对于不同的裂缝密度,各向异性梯度Bani与各向同性梯度Biso的比值I(1)fluid近似为常数,且对流体敏感.经裂缝纵横比和背景介质拉梅常数修正后,流体因子Ifluid既不随纵横比变化,又不受背景介质的影响,是裂缝型储层敏感的流体识别因子.在塔里木盆地塔北哈拉哈塘地区热瓦普区块碳酸盐岩储层裂缝发育区域,运用该参数在井点处的流体识别效果与钻井结果一致.  相似文献   

18.
正交各向异性介质P波走时分析及Thomsen参数反演   总被引:10,自引:3,他引:10  
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层,正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时,时距曲线表现出强方位各向异性.采用小生境遗传算法,对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明,此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数.  相似文献   

19.
Frequency attenuation occurs when seismic waves propagate through the porous reservoirs containing hydrocarbons. Current researches on the seismic frequency attenuation mainly focus on the post-stack domain instead of the prestack domain. Here we propose the frequency attenuation gradient vs. offset (FAGVO) based on the amplitude variation with offset and frequency attenuation integral equations. We derive the FAGVO equation that equals to zero in a full-elastic medium and is negative in a viscoelastic medium. FAGVO is affected by the viscosity of the medium, the coefficients of reflection, the frequency variation, and high-frequency attenuation. FAGVO uses the differences of partially stacked data to decrease the interference caused by subsurface strata affecting the frequency attenuation, highlights the frequency attenuation gradient anomalies in hydrocarbon-bearing reservoir pores, and finally realizes the hydrocarbon fluid identification. The method was verified using a two-dimensional wave equation forward model and was found to be cost effective. Furthermore, the method does not require well information, which can be applied in the stage of seismic exploration, especially, in the exploration of a none-well project.  相似文献   

20.
黏弹各向异性介质中波的反射与透射问题分析   总被引:6,自引:1,他引:5       下载免费PDF全文
黏弹各向异性介质中传播不均匀波,其反射、透射模式不仅与介质分界面两侧速度对比有关,还与品质因子Q的对比有关. 用伪谱技术模拟黏弹各向异性介质分界面上波的反射、透射,并与弹性各向异性介质、黏弹各向同性介质和弹性各向同性介质的模拟结果做比较. 计算平面波的反射、透射系数,分析介质的黏弹性和各向异性对反射、透射系数的影响. 数值模拟了一个三层介质模型中的波场,分析两个分界面上产生的反射波的特征. 黏弹各向异性介质中,qS波比qP波衰减程度大.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号