首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
天然气在开发过程中,储层有效压力和含气饱和度均会发生变化,研究有效压力和含气饱和度的变化对地震响应特征的影响,在基于时移地震的剩余气分布预测研究中具有重要意义。天然气和石油的声学性质有着明显的差异,油藏时移地震的研究成果不能直接应用于气藏,因此需要开展气藏的时移地震研究。利用Shapiro模型表征干岩石弹性模量随有效压力的变化,借助Batzle-Wang方程描述流体速度随压力的变化关系,联合Gassmann理论进行流体替代,表征饱和流体岩石速度随含气饱和度的变化,建立了饱和流体岩石速度随有效压力和饱和度变化的岩石物理模型。基于该模型,对不同含气饱和度和不同有效压力下的气藏储层模型进行了多波时移地震叠前振幅变化（AVO）模拟。结果表明多波时移地震AVO技术可以有效地区分有效压力变化和含气饱和度变化,为进一步开展气藏多波时移地震流体监测提供了理论参考依据。   相似文献   

2.
Quantitative interpretation of time‐lapse seismic data requires knowledge of the relationship between elastic wave velocities and fluid saturation. This relationship is not unique but depends on the spatial distribution of the fluid in the pore‐space of the rock. In turn, the fluid distribution depends on the injection rate. To study this dependency, forced imbibition experiments with variable injection rates have been performed on an air‐dry limestone sample. Water was injected into a cylindrical sample and was monitored by X‐Ray Computed Tomography and ultrasonic time‐of‐flight measurements across the sample. The measurements show that the P‐wave velocity decreases well before the saturation front approaches the ultrasonic raypath. This decrease is followed by an increase as the saturation front crosses the raypath. The observed patterns of the acoustic response and water saturation as functions of the injection rate are consistent with previous observations on sandstone. The results confirm that the injection rate has significant influence on fluid distribution and the corresponding acoustic response. The complexity of the acoustic response —‐ that is not monotonic with changes in saturation, and which at the same saturation varies between hydrostatic conditions and states of dynamic fluid flow – may have implications for the interpretation of time‐lapse seismic responses.  相似文献   

3.
Quantitative detection of fluid distribution using time-lapse seismic   总被引：1，自引：0，他引：1  
Although previous seismic monitoring studies have revealed several relationships between seismic responses and changes in reservoir rock properties, the quantitative evaluation of time‐lapse seismic data remains a challenge. In most cases of time‐lapse seismic analysis, fluid and/or pressure changes are detected qualitatively by changes in amplitude strength, traveltime and/or Poisson's ratio. We present the steps for time‐lapse seismic analysis, considering the pressure effect and the saturation scale of fluids. We then demonstrate a deterministic workflow for computing the fluid saturation in a reservoir in order to evaluate time‐lapse seismic data. In this approach, we derive the physical properties of the water‐saturated sandstone reservoir, based on the following inputs: VP, VS, ρ and the shale volume from seismic analysis, the average properties of sand grains, and formation‐water properties. Next, by comparing the in‐situ fluid‐saturated properties with the 100% formation‐water‐saturated reservoir properties, we determine the bulk modulus and density of the in‐situ fluid. Solving three simultaneous equations (relating the saturations of water, oil and gas in terms of the bulk modulus, density and the total saturation), we compute the saturation of each fluid. We use a real time‐lapse seismic data set from an oilfield in the North Sea for a case study.  相似文献   

4.
The hydrodynamic characterization of the epikarst, the shallow part of the unsaturated zone in karstic systems, has always been challenging for geophysical methods. This work investigates the feasibility of coupling time‐lapse refraction seismic data with petrophysical and hydrologic models for the quantitative determination of water storage and residence time at shallow depth in carbonate rocks. The Biot–Gassmann fluid substitution model describing the seismic velocity variations with water saturation at low frequencies needs to be modified for this lithology. I propose to include a saturation‐dependent rock‐frame weakening to take into account water–rock interactions. A Bayesian inversion workflow is presented to estimate the water content from seismic velocities measured at variable saturations. The procedure is tested first with already published laboratory measurements on core samples, and the results show that it is possible to estimate the water content and its uncertainty. The validated procedure is then applied to a time‐lapse seismic study to locate and quantify seasonal water storage at shallow depth along a seismic profile. The residence time of the water in the shallow layers is estimated by coupling the time‐lapse seismic measurements with rainfall chronicles, simple flow equations, and the petrophysical model. The daily water input computed from the chronicles is used to constraint the inversion of seismic velocities for the daily saturation state and the hydrodynamic parameters of the flow model. The workflow is applied to a real monitoring case, and the results show that the average residence time of the water in the epikarst is generally around three months, but it is only 18 days near an infiltration pathway. During the winter season, the residence times are three times shorter in response to the increase in the effective rainfall.  相似文献   

5.
Saturation of porous rocks with a mixture of two fluids has a substantial effect on seismic‐wave propagation. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves due to the mechanism of wave‐induced fluid‐flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a regular way. However, recent experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Based on an approximation for the coherent wavefield in random porous media, we develop a model which assumes a continuous distribution of fluid heterogeneities. As this continuous random media approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental data. We also show how to relate the random functions to experimentally measurable parameters.  相似文献   

6.
Seismic wave propagation through a fluid-saturated poroelastic layer might be strongly affected by media heterogeneities. Via incorporating controlled laboratory simulation experiments, we extend previous studies of time-lapse seismic effects to evaluate the wave scattering influence of the heterogeneous nature of porous permeable media and the associated amplification effects on 4D seismic response characteristics of reservoir fluid substitution. A physical model consisted of stratified thin layers of shale and porous sandstone reservoir with rock heterogeneities was built based on the geological data of a real hydrocarbon-saturated reservoir in Northeast China. Multi-surveys data of good quality were acquired by filling poroelastic reservoir layers with gas, water and oil in sequence. Experimental observations show that reservoir heterogeneity effect causes significantly magnified abnormal responses to the fluid-saturated media. Specifically, reflection signatures of the gas-filled reservoir are dramatically deviated from those of the liquid fluid-filled reservoir, compared with ones of the homogeneous media. By removing the influences unrelated to reservoir property alterations, 4D seismic estimates of travel-time and frequency-dependent characteristic are reasonably consistent with fluid variations. Nevertheless, strong 4D amplitude difference anomalies might not correspond to the regions where fluid variations occur. We also find that 4D seismic difference attributes are evident between oil- and water-filled models, whereas significant between oil- and gas-filled models. Meanwhile, rock physics modelling results reveal the predicted 4D seismic differences are obviously smaller than those calculated from seismic observations. The results in this paper, therefore, implicate that the effect of a reservoir's heterogeneous nature might be beneficial for hydrocarbons detection as well as monitoring small variations in pore fluids.  相似文献   

7.
Common shot ray tracing and finite difference seismic modelling experiments were undertaken to evaluate variations in the seismic response of the Devonian Redwater reef in the Alberta Basin, Canada after replacement of native pore waters in the upper rim of the reef with CO2. This part of the reef is being evaluated for a CO2 storage project. The input geological model was based on well data and the interpretation of depth‐converted, reprocessed 2D seismic data in the area. Pre‐stack depth migration of the ray traced and finite difference synthetic data demonstrate similar seismic attributes for the Mannville, Nisku, Ireton, Cooking Lake, and Beaverhill Lake formations and clear terminations of the Upper Leduc and Middle Leduc events at the reef margin. Higher amplitudes at the base of Upper‐Leduc member are evident near the reef margin due to the higher porosity of the foreslope facies in the reef rim compared to the tidal flat lagoonal facies within the central region of the reef. Time‐lapse seismic analysis exhibits an amplitude difference of about 14% for Leduc reflections before and after CO2 saturation and a travel‐time delay through the reservoir of 1.6 ms. Both the ray tracing and finite difference approaches yielded similar results but, for this particular model, the latter provided more precise imaging of the reef margin. From the numerical study we conclude that time‐lapse surface seismic surveys should be effective in monitoring the location of the CO2 plume in the Upper Leduc Formation of the Redwater reef, although the differences in the results between the two modelling approaches are of similar order to the effects of the CO2 fluid replacement itself.  相似文献   

8.
Knowledge about saturation and pressure distributions in a reservoir can help in determining an optimal drainage pattern, and in deciding on optimal well designs to reduce risks of blow‐outs and damage to production equipment. By analyzing time‐lapse PP AVO or time‐lapse multicomponent seismic data, it is possible to separate the effects of production related saturation and pressure changes on seismic data. To be able to utilize information about saturation and pressure distributions in reservoir model building and simulation, information about uncertainty in the estimates is useful. In this paper we present a method to estimate changes in saturation and pressure from time‐lapse multicomponent seismic data using a Bayesian estimation technique. Results of the estimations will be probability density functions (pdfs), giving immediate information about both parameter values and uncertainties. Linearized rock physical models are linked to the changes in saturation and pressure in the prior probability distribution. The relationship between the elastic parameters and the measured seismic data is described in the likelihood model. By assuming Gaussian distributed prior uncertainties the posterior distribution of the saturation and pressure changes can be calculated analytically. Results from tests on synthetic seismic data show that this method produces more precise estimates of changes in effective pressure than a similar methodology based on only PP AVO time‐lapse seismic data. This indicates that additional information about S‐waves obtained from converted‐wave seismic data is useful for obtaining reliable information about the pressure change distribution.  相似文献   

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

10.
P‐wave data from a time‐lapse 3D OBC survey have been analysed to estimate and interpret azimuthal seismic anisotropy. This is achieved by careful processing to preserve the azimuthal signature. The survey images a major reservoir body in a channelized turbidite field in the Gulf of Mexico. Three distinct and significant anisotropy anomalies are discovered on or around this particular ‘4500‐ft sand’, all of which change intensity but not orientation with hydrocarbon production. These anomalies are distributed along the highest concentration of cumulative sand thickness, with their symmetry axes aligned with the main channel axis. We suspect that this time‐lapse anisotropy could be caused by the alignment of the depositional grain fabric. Theoretical calculation predicts that this mechanism, when combined with fluid‐saturation changes, can generate the observed pattern of behaviour. If further supported by other researchers, this result would indicate that appropriately designed seismic surveys could be a useful tool for palaeo‐direction studies in clastic reservoirs and also a useful constraint for directional permeability in the reservoir flow simulation model.  相似文献   

11.
Ghawar, the largest oilfield in the world, produces oil from the Upper Jurassic Arab‐D carbonate reservoir. The high rigidity of the limestone–dolomite reservoir rock matrix and the small contrast between the elastic properties of the pore fluids, i.e. oil and water, are responsible for the weak 4D seismic effect due to oil production. A feasibility study was recently completed to quantify the 4D seismic response of reservoir saturation changes as brine replaced oil. The study consisted of analysing reservoir rock physics, petro‐acoustic data and seismic modelling. A seismic model of flow simulation using fluid substitution concluded that time‐lapse surface seismic or conventional 4D seismic is unlikely to detect the floodfront within the repeatability of surface seismic measurements. Thus, an alternative approach to 4D seismic for reservoir fluid monitoring is proposed. Permanent seismic sensors could be installed in a borehole and on the surface for passive monitoring of microseismic activity from reservoir pore‐pressure perturbations. Reservoir production and injection operations create these pressure or stress perturbations. Reservoir heterogeneities affecting the fluid flow could be mapped by recording the distribution of epicentre locations of these microseisms or small earthquakes. The permanent borehole sensors could also record repeated offset vertical seismic profiling surveys using a surface source at a fixed location to ensure repeatability. The repeated vertical seismic profiling could image the change in reservoir properties with production.  相似文献   

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

13.
We address the issue of linearity and scale dependence in forward modelling of seismic data from well logs, for large ray parameters, wide angles or large offsets. We present a forward model, within the context of seismic‐to‐well matching, that is linearized in the elastic properties of the earth. This model preserves linearity at large ray parameters and can handle fine‐layering effects such as induced anisotropy. Starting from a low‐contrast small‐ray‐parameter model, we extend it to a large‐ray‐parameter model by fully linearizing the elastic‐property contrasts. Overall linearity of the forward model is extended by partitioning the compressional‐wave and shear‐wave velocity fields into two fundamental scales: a kinematic scale that governs wavefield propagation effects and a dynamic scale that governs wavefield scattering effects. This analysis reveals that the standard practice in forward modelling of strongly filtering the ratio of compressional‐wave velocity to shear‐wave velocity is well founded in the underlying physics. The partitioning of the velocity fields also leads naturally to forward modelling that accounts fully for stretch effects, to resolution of the angle‐of‐incidence versus ray‐parameter dichotomy in seismic‐amplitude analysis, and to full accounting for induced anisotropy and dispersion effects due to fine‐layering of isotropic media. With the onset of routine long‐offset acquisition and the compelling need to optimize asset management in order to maximize reserve recovery, this forward model recognizes the physics of seismic wave propagation and enables a more complete exploitation of amplitude information in pre‐critical seismic data.  相似文献   

14.
Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated.  相似文献   

15.
Combined time‐lapse reservoir simulation and seismic modelling has been performed on both 1D and 3D models of a channelized turbidite reservoir. The models have been built using core, log, laboratory and seismic data from the Nelson Field (central North Sea) as a template. Oil and water movement in the main channels, channel margins and interchannel regions is investigated, with a particular focus being the effect of poor net‐to‐gross. The analysis confirms that saturation effects dominate the response whilst stress‐sensitivity effects play a minor role. The trough–peak signature in the seismic difference volumes formed by the sweep of the water can be continued and mapped slightly further than the channel margins. This characteristic 4D signature remains roughly intact, despite the complicated depositional architecture, and accurately delineates the area of moved fluid, but it needs additional calibration to combat the detrimental influence of the low net‐to‐gross. Signal strength is largely dependent on reservoir quality, but is also compromised by the net‐to‐gross, fluid distribution and, more critically, by the exact timing of the seismic survey. For example, a region of bypassed oil zone remains undetected as it forms early during the production. This work demonstrates that to understand fully the 4D signature at a quantitative level requires adequate knowledge of the fluid properties, but also, more critically, the geology.  相似文献   

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

17.
Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C33 (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C33. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.  相似文献   

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

19.
基于传播矩阵理论开发砂泥岩薄互层地震合成记录算法,与褶积算法、基于界面模型的Zoeppritz方法以及波动方程等方法相比,该方法更适用于具有复杂结构的薄互层模型,能够在充分考虑地震反射波动力学因素的同时不受网格间距的限制.基于正演算法开发了基于波形对比的砂泥岩薄互层地震反演技术,由地震反射波形特征的变化反演薄地层单元中砂体含量与空间位置等参数,进而确定薄互层段砂泥岩组合结构以及砂体的空间展布.理论模型验证了反演方法的有效性.通过测井分析建立薄互层地震地质模型,并将该技术应用于研究区实际地震数据,反演的砂体空间分布与测井资料进行对比分析,验证了反演方法的实用性.  相似文献   

20.
介观尺度孔隙流体流动是地震频段岩石表现出较强速度频散与衰减的主要作用.利用周期性层状孔隙介质模型,基于准静态孔弹性理论给出了模型中孔隙压力、孔隙流体相对运动速度以及固体骨架位移等物理量的数学解析表达式,同时利用Biot理论将其扩展至全频段条件下,克服了传统White模型中介质分界面处流体压力不连续的假设. 在此基础上对准静态与全频段下模型介质中孔隙压力、孔隙流体相对运动速度变化形式及其对弹性波传播特征的影响进行了讨论,为更有效理解介观尺度下流体流动耗散和频散机制提供物理依据.研究结果表明,低频条件下快纵波孔压在介质层内近于定值,慢纵波通过流体扩散改变总孔隙压力, 随频率的增加慢波所形成的流体扩散作用逐渐减弱致使介质中总孔压逐渐接近于快纵波孔压,在较高频率下孔压与应力的二次耦合作用使总孔压超过快纵波孔压.介质中孔隙流体相对运动速度与慢纵波形成的流体相对运动速度变化形式一致;随频率的增加孔隙流体逐渐从排水的弛豫状态过渡到非弛豫状态,其纵波速度-含水饱和度变化形式也从符合孔隙流体均匀分布模式过渡到斑块分布模式,同时介质在不同含水饱和度下的衰减峰值与慢纵波所形成的孔隙流体相对流动速度具有明显的相关性.  相似文献   

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

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