首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
We extend the frequency‐ and angle‐dependent poroelastic reflectivity to systematically analyse the characteristic of seismic waveforms for highly attenuating reservoir rocks. It is found that the mesoscopic fluid pressure diffusion can significantly affect the root‐mean‐square amplitude, frequency content, and phase signatures of seismic waveforms. We loosely group the seismic amplitude‐versus‐angle and ‐frequency characteristics into three classes under different geological circumstances: (i) for Class‐I amplitude‐versus‐angle and ‐frequency, which corresponds to well‐compacted reservoirs having Class‐I amplitude‐versus‐offset characteristic, the root‐mean‐square amplitude at near offset is boosted at high frequency, whereas seismic energy at far offset is concentrated at low frequency; (ii) for Class‐II amplitude‐versus‐angle and ‐frequency, which corresponds to moderately compacted reservoirs having Class‐II amplitude‐versus‐offset characteristic, the weak seismic amplitude might exhibit a phase‐reversal trend, hence distorting both the seismic waveform and energy distribution; (iii) for Class‐III amplitude‐versus‐angle and ‐frequency, which corresponds to unconsolidated reservoir having Class‐III amplitude‐versus‐offset characteristic, the mesoscopic fluid flow does not exercise an appreciable effect on the seismic waveforms, but there exists a non‐negligible amplitude decay compared with the elastic seismic responses based on the Zoeppritz equation.  相似文献   

2.
Various models have been proposed to link partial gas saturation to seismic attenuation and dispersion, suggesting that the reflection coefficient should be frequency‐dependent in many cases of practical importance. Previous approaches to studying this phenomenon typically have been limited to single‐interface models. Here, we propose a modelling technique that allows us to incorporate frequency‐dependent reflectivity into convolutional modelling. With this modelling framework, seismic data can be synthesised from well logs of velocity, density, porosity, and water saturation. This forward modelling could act as a basis for inversion schemes aimed at recovering gas saturation variations with depth. We present a Bayesian inversion scheme for a simple thin‐layer case and a particular rock physics model and show that, although the method is very sensitive to prior information and constraints, both gas saturation and layer thickness theoretically can be estimated in the case of interfering reflections.  相似文献   

3.
Hydrocarbon prediction from seismic amplitude and amplitude‐versus‐offset is a daunting task. Amplitude interpretation is ambiguous due to the effects of lithology and pore fluid. In this paper, we propose a new attribute “J” based on a Gassmann–Biot fluid substitution to reduce ambiguity. Constrained by seismic and rock physics, the J attribute has good ability to detect hydrocarbons from seismic data. There are currently many attributes for hydrocarbon prediction. Among the existing attributes, far‐minus‐near times far and fluid factor are commonly used. In this paper, the effectiveness of these two existing attributes was compared with the new attribute. Numerical modelling was used to test the new attribute “J” and to compare “J” with the two existing attributes. The results showed that the J attribute can predict the existence of hydrocarbon in different porosity scenarios with less ambiguity than the other two attributes. Tests conducted with real seismic data demonstrated the effectiveness of the J attribute. The J attribute has performed well in scenarios in which the other two attributes gave inaccurate predictions. The proposed attribute “J” is fast and simple, and it could be used as a first step in hydrocarbon analysis for exploration.  相似文献   

4.
优质气层与低含气饱和度的差气层、气水同层均呈现"亮点"强振幅的地震响应特征,利用振幅类属性或者常规弹性参数难以进行有效判识。考虑到密度参数随含气饱和度的高敏感性,本文在贝叶斯叠前反演框架下,以三变量柯西分布为先验正则约束进行叠前地震反演,减小了纵波速度、横波速度与密度三参数统计相关性造成的参数估算误差,进一步提高了弹性参数的可靠性,模型测试表明新方法反演的密度精度及稳定性更高。将方法应用于马达班湾盆地X区块,在流体替代及岩石物理统计分析的指导下,利用基于三变量柯西分布的叠前地震反演方法提取密度参数,预测优质气层富集区,实钻结果与预测结果吻合度高,表明该方法具有较高的实用性。   相似文献   

5.
含气饱和度预测是天然气储层地震解释工作的重要目标.本文将岩石物理分析与地震物理模拟技术相结合,构建了部分;饱和砂岩储层物理模型并进行含气饱和度预测分析.物理模型中设置了高孔渗常规砂岩和低孑孔渗致密砂岩两种模拟储层,每种储层都是由具有不同含水饱和度的气-水双相饱和砂体组成.岩石物理分析结果显示在低孔渗致密砂岩中气-水混合流体更加倾向于非均匀的斑块分布,而结合了Brie等效流体公式的Gassmann流体替换理论可以更准确地描述纵波速度随含水饱和度的变化趋势.对物理模型进行地震资料采集处理后,对比了AVO特征和叠前同步反演结果对两种砂岩储层含气饱和度预测能力的差异.AVO特征结果显示,对于混合流体均匀分布的高孔渗砂岩储层,AVO响应曲线和属性变化很难对含气饱和度进行估算;对于混合流体斑块分布的致密砂岩储层,AVO特征可以定性地分辨出储层是否为高、中、低含气情况.反演结果显示,密度及纵横波速度比分别对高孔渗及致密砂岩储层的含气饱和度有着较好的指示能力.  相似文献   

6.
For 3‐D shallow‐water seismic surveys offshore Abu Dhabi, imaging the target reflectors requires high resolution. Characterization and monitoring of hydrocarbon reservoirs by seismic amplitude‐versus‐offset techniques demands high pre‐stack amplitude fidelity. In this region, however, it still was not clear how the survey parameters should be chosen to satisfy the required data quality. To answer this question, we applied the focal‐beam method to survey evaluation and design. This subsurface‐ and target‐oriented approach enables quantitative analysis of attributes such as the best achievable resolution and pre‐stack amplitude fidelity at a fixed grid point in the subsurface for a given acquisition geometry at the surface. This method offers an efficient way to optimize the acquisition geometry for maximum resolution and minimum amplitude‐versus‐offset imprint. We applied it to several acquisition geometries in order to understand the effects of survey parameters such as the four spatial sampling intervals and apertures of the template geometry. The results led to a good understanding of the relationship between the survey parameters and the resulting data quality and identification of the survey parameters for reflection imaging and amplitude‐versus‐offset applications.  相似文献   

7.
Fluid depletion within a compacting reservoir can lead to significant stress and strain changes and potentially severe geomechanical issues, both inside and outside the reservoir. We extend previous research of time‐lapse seismic interpretation by incorporating synthetic near‐offset and full‐offset common‐midpoint reflection data using anisotropic ray tracing to investigate uncertainties in time‐lapse seismic observations. The time‐lapse seismic simulations use dynamic elasticity models built from hydro‐geomechanical simulation output and a stress‐dependent rock physics model. The reservoir model is a conceptual two‐fault graben reservoir, where we allow the fault fluid‐flow transmissibility to vary from high to low to simulate non‐compartmentalized and compartmentalized reservoirs, respectively. The results indicate time‐lapse seismic amplitude changes and travel‐time shifts can be used to qualitatively identify reservoir compartmentalization. Due to the high repeatability and good quality of the time‐lapse synthetic dataset, the estimated travel‐time shifts and amplitude changes for near‐offset data match the true model subsurface changes with minimal errors. A 1D velocity–strain relation was used to estimate the vertical velocity change for the reservoir bottom interface by applying zero‐offset time shifts from both the near‐offset and full‐offset measurements. For near‐offset data, the estimated P‐wave velocity changes were within 10% of the true value. However, for full‐offset data, time‐lapse attributes are quantitatively reliable using standard time‐lapse seismic methods when an updated velocity model is used rather than the baseline model.  相似文献   

8.
Seismic conditioning of static reservoir model properties such as porosity and lithology has traditionally been faced as a solution of an inverse problem. Dynamic reservoir model properties have been constrained by time‐lapse seismic data. Here, we propose a methodology to jointly estimate rock properties (such as porosity) and dynamic property changes (such as pressure and saturation changes) from time‐lapse seismic data. The methodology is based on a full Bayesian approach to seismic inversion and can be divided into two steps. First we estimate the conditional probability of elastic properties and their relative changes; then we estimate the posterior probability of rock properties and dynamic property changes. We apply the proposed methodology to a synthetic reservoir study where we have created a synthetic seismic survey for a real dynamic reservoir model including pre‐production and production scenarios. The final result is a set of point‐wise probability distributions that allow us to predict the most probable reservoir models at each time step and to evaluate the associated uncertainty. Finally we also show an application to real field data from the Norwegian Sea, where we estimate changes in gas saturation and pressure from time‐lapse seismic amplitude differences. The inverted results show the hydrocarbon displacement at the times of two repeated seismic surveys.  相似文献   

9.
Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.  相似文献   

10.
Elastic behaviour of North Sea chalk: A well-log study   总被引：1，自引：1，他引：0  
We present two different elastic models for, respectively, cemented and uncemented North Sea chalk well‐log data. We find that low Biot coefficients correlate with anomalously low cementation factors from resistivity measurements at low porosity and we interpret this as an indication of cementation. In contrast, higher Biot coefficients and correspondingly higher cementation factors characterize uncemented chalk for the same (low) porosity. Accordingly, the Poisson's ratio–porosity relationship for cemented chalk is different from that of uncemented chalk. We have tested the application of the self‐consistent approximation, which here represents the unrelaxed scenario where the pore spaces of the rock are assumed to be isolated, and the Gassmann theory, which assumes that pore spaces are connected, as tools for predicting the effect of hydrocarbons from the elastic properties of brine‐saturated North Sea reservoir chalk. In the acoustic impedance–Poisson's ratio plane, we forecast variations in porosity and hydrocarbon saturation from their influence on the elastic behaviour of the chalk. The Gassmann model and the self‐consistent approximation give roughly similar predictions of the effect of fluid on acoustic impedance and Poisson's ratio, but we find that the high‐frequency self‐consistent approach gives a somewhat smaller predicted fluid‐saturation effect on Poisson's ratio than the low‐frequency Gassmann model. The Gassmann prediction for the near and potentially invaded zone corresponds more closely to logging data than the Gassmann prediction for the far, virgin zone. We thus conclude that the Gassmann approach predicts hydrocarbons accurately in chalk in the sonic‐frequency domain, but the fluid effects as recorded by the acoustic tool are significantly affected by invasion of mud filtrate. The amplitude‐versus‐angle (AVA) response for the general North Sea sequence of shale overlying chalk is predicted as a function of porosity and pore‐fill. The AVA response of both cemented and uncemented chalk generally shows a declining reflectivity coefficient versus offset and a decreasing normal‐incidence reflectivity with increasing porosity. However, for the uncemented model, a phase reversal will appear at a relatively lower porosity compared to the cemented model.  相似文献   

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

12.
In heterogeneous natural gas reservoirs, gas is generally present as small patch-like pockets embedded in the water-saturated host matrix. This type of heterogeneity, alsocalled ＂patchy saturation＂, causes significant seismic velocity dispersion and attenuation. Toestablish the relation between seismic response and type of fluids, we designed a rock physicsmodel for carbonates. First, we performed CT scanning and analysis of the fluid distributionin the partially saturated rocks. Then, we predicted the quantitative relation between the waveresponse at different frequency ranges and the basic lithological properties and pore fluids.A rock physics template was constructed based on thin section analysis of pore structuresand seismic inversion. This approach was applied to the limestone gas reservoirs of the rightbank block of the Amu Darya River. Based on poststack wave impedance and prestack elasticparameter inversions, the seismic data were used to estimate rock porosity and gas saturation.The model results were in ~ood a~reement with the production regime of the wells.  相似文献   

13.
We study the appraisal problem for the joint inversion of seismic and controlled source electro‐magnetic (CSEM) data and utilize rock‐physics models to integrate these two disparate data sets. The appraisal problem is solved by adopting a Bayesian model and we incorporate four representative sources of uncertainty. These are uncertainties in 1) seismic wave velocity, 2) electric conductivity, 3) seismic data and 4) CSEM data. The uncertainties in porosity and water saturation are quantified by a posterior random sampling in the model space of porosity and water saturation in a marine one‐dimensional structure. We study the relative contributions from the four individual sources of uncertainty by performing several statistical experiments. The uncertainties in the seismic wave velocity and electric conductivity play a more significant role on the variation of posterior uncertainty than do the seismic and CSEM data noise. The numerical simulations also show that the uncertainty in porosity is most affected by the uncertainty in the seismic wave velocity and that the uncertainty in water saturation is most influenced by the uncertainty in electric conductivity. The framework of the uncertainty analysis presented in this study can be utilized to effectively reduce the uncertainty of the porosity and water saturation derived from the integration of seismic and CSEM data.  相似文献   

14.
频变AVO含气性识别技术研究与应用   总被引：2，自引：6，他引：2       下载免费PDF全文
在常规AVO理论的基础上,频变AVO属性计算方法利用多尺度裂缝介质模型中物性参数具有频变特征的优势,基于Zoeppritz方程建立反射系数与频率之间的数学关系,推导出截距、梯度、碳氢检测因子、流体检测因子、拟泊松比等AVO属性与频率之间的数学关系.应用地震反演方法,综合地质、地震、测井等数据,反演出高精度的频变AVO属性,在天然气敏感属性分析的基础上建立起频变AVO含气性识别技术.将该技术应用到川西新场陆相深层须家河组碎屑岩储层的含气性识别中,在孔隙度通常在1%~4%,渗透率普遍低于0.06×10-3µm2的致密背景下,较准确地预测了孔隙度大于4％、渗透性偏高的富气优质储层分布带,为该区含气性识别难题的解决和钻井成功率的提高,提供了重要的技术支撑.  相似文献   

15.
An approach is developed to estimate pore‐pressure changes in a compacting chalk reservoir directly from time‐lapse seismic attributes. It is applied to data from the south‐east flank of the Valhall field. The time‐lapse seismic signal of the reservoir in this area is complex, despite the fact that saturation changes do not have an influence. This complexity reflects a combination of pressure depletion, compaction and stress re‐distribution throughout the reservoir and into the surrounding rocks. A simple relation is found to link the time‐lapse amplitude and time‐shift attributes to variations in the key controlling parameter of initial porosity. This relation is sufficient for an accurate estimation of pore‐pressure change in the inter‐well space. Although the time‐lapse seismic estimates mostly agree with reservoir simulation, unexplained mismatches are apparent at a small number of locations with lower porosities (less than 38%). The areas of difference between the observations and predictions suggest possibilities for simulation model updating or a better understanding of the physics of the reservoir.  相似文献   

16.
The injection of CO2 at the Ketzin pilot site commenced in June 2008 and was terminated in August 2013 after 67 kT had been injected into a saline formation at a depth of 630–650 m. As part of the site monitoring program, four 3D surface seismic surveys have been acquired to date, one baseline and three repeats, of which two were conducted during the injection period, and one during the post‐injection phase. The surveys have provided the most comprehensive images of the spreading CO2 plume within the reservoir layer. Both petrophysical experiments on core samples from the Ketzin reservoir and spectral decomposition of the 3D time‐lapse seismic data show that the reservoir pore pressure change due to CO2 injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO2 saturation. In this study, amplitude versus offset analysis has been applied to investigate the amplitude versus offset response from the top of the sandstone reservoir during the injection and post‐injection phases, and utilize it to obtain a more quantitative assessment of the CO2 gaseous saturation changes. Based on the amplitude versus offset modelling, a prominent decrease in the intercept values imaged at the top of the reservoir around the injection well is indeed associated solely with the CO2 saturation increase. Any change in the gradient values, which would, in case it was positive, be the only signature induced by the reservoir pressure variations, has not been observed. The amplitude versus offset intercept change is, therefore, entirely ascribed to CO2 saturation and used for its quantitative assessment. The estimated CO2 saturation values around the injection area in the range of 40%–60% are similar to those obtained earlier from pulsed neutron‐gamma logging. The highest values of 80% are found in the second seismic repeat in close vicinity to the injection and observation wells.  相似文献   

17.
In fractured reservoirs, seismic wave velocity and amplitude depend on frequency and incidence angle. Frequency dependence is believed to be principally caused by the wave‐induced flow of pore fluid at the mesoscopic scale. In recent years, two particular phenomena, i.e., patchy saturation and flow between fractures and pores, have been identified as significant mechanisms of wave‐induced flow. However, these two phenomena are studied separately. Recently, a unified model has been proposed for a porous rock with a set of aligned fractures, with pores and fractures filled with two different fluids. Existing models treat waves propagating perpendicular to the fractures. In this paper, we extend the model to all propagation angles by assuming that the flow direction is perpendicular to the layering plane and is independent of the loading direction. We first consider the limiting cases through poroelastic Backus averaging, and then we obtain the five complex and frequency‐dependent stiffness values of the equivalent transversely isotropic medium as a function of the frequency. The numerical results show that, when the bulk modulus of the fracture‐filling fluid is relatively large, the dispersion and attenuation of P‐waves are mainly caused by fractures, and the values decrease as angles increase, almost vanishing when the incidence angle is 90° (propagation parallel to the fracture plane). While the bulk modulus of fluid in fractures is much smaller than that of matrix pores, the attenuation due to the “partial saturation” mechanism makes the fluid flow from pores into fractures, which is almost independent of the incidence angle.  相似文献   

18.
Time‐lapse 3D seismic reflection data, covering the CO2 storage operation at the Snøhvit gas field in the Barents Sea, show clear amplitude and time‐delay differences following injection. The nature and extent of these changes suggest that increased pore fluid pressure contributes to the observed seismic response, in addition to a saturation effect. Spectral decomposition using the smoothed pseudo‐Wigner–Ville distribution has been used to derive discrete‐frequency reflection amplitudes from around the base of the CO2 storage reservoir. These are utilized to determine the lateral variation in peak tuning frequency across the seismic anomaly as this provides a direct proxy for the thickness of the causative feature. Under the assumption that the lateral and vertical extents of the respective saturation and pressure changes following CO2 injection will be significantly different, discrete spectral amplitudes are used to distinguish between the two effects. A clear spatial separation is observed in the distribution of low‐ and high‐frequency tuning. This is used to discriminate between direct fluid substitution of CO2, as a thin layer, and pressure changes that are distributed across a greater thickness of the storage reservoir. The results reveal a striking correlation with findings derived from pressure and saturation discrimination algorithms based on amplitude versus offset analysis.  相似文献   

19.
The cross‐calibration of different vintage data is an important prerequisite in attempting to determine the time‐lapse seismic effects induced by hydrocarbon production in a reservoir. This paper reports the preprocessing and cross‐calibration procedures adopted to modify the data of four seismic vintages (1982, 1989, 1992 and 1999) from the Oseberg field in the North Sea, for optimal conditions for a time‐lapse seismic amplitude analysis. The final results, in terms of time‐lapse variations, of acoustic impedance and of amplitude‐versus‐offset, are illustrated for selected data sets. The application of preprocessing to each individual vintage data set reduces the effects of the different acquisition and noise conditions, and leads to consistency in the amplitude response of the four vintages. This consistency facilitates the final amplitude cross‐calibration that is carried out using, as reference, the Cretaceous horizon reflections above the Brent reservoir. Such cross‐calibration can be considered as vintage‐consistent residual amplitude correction. Acoustic impedance sections, intercept and gradient amplitude‐versus‐offset attributes and coherent amplitude‐versus‐offset estimates are computed on the final cross‐calibrated data. The results, shown for three spatially coincident 2D lines selected from the 1982, 1989 and 1999 data sets, clearly indicate gas‐cap expansion resulting from oil production. Such expansion is manifested as a decrease in acoustic impedance and a modification of the amplitude‐versus‐offset trends in the apical part of the reservoir.  相似文献   

20.
The propagation of seismic waves through a saturated reservoir compresses the fluid in the pore spaces. During this transition, parts of seismic energy would be attenuated because of intrinsic absorption. Rock physics models make the bridge between the seismic properties and petrophysical reality in the earth. Attenuation is one of the significant seismic attributes used to describe the fluid behaviour in the reservoirs. We examined the core samples using ultrasonic experiments at the reservoir conditions. Given the rock properties of the carbonate reservoir and experiment results, the patchy saturation mechanism was solved for substituted fluid using the theory of modulus frequency. The extracted relationship between the seismic attenuation and water saturation was used in time–frequency analysis. We performed the peak frequency method to estimate the Q factor in the Gabor domain and determined the water saturation based on the computed rock physics model. The results showed how the probable fault in the reservoir has stopped the fluid movement in the reservoir and caused touching the water‐bearing zone through drilling.  相似文献   

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

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