首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
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.  相似文献   

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

3.
Time-lapse seismic analysis of pressure depletion in the Southern Gas Basin   总被引：1，自引：0，他引：1  
In the Southern Gas Basin (SGB) of the North Sea there are many mature gas fields where time‐lapse monitoring could be very beneficial in extending production life. However, the conditions are not immediately attractive for time‐lapse seismic assessment. This is primarily because the main production effect to be assessed is a pore pressure reduction and frame stiffening because of gas production in tight sandstone reservoirs that also have no real seismic direct hydrocarbon indicators. Modelling, based on laboratory measurements, has shown that such an effect would be small and difficult to detect in seismic data. This paper makes two main contributions. Firstly, this is, to our knowledge, the first time‐lapse study in the SGB and involves a real‐data assessment of the viability for detecting production in such an environment. Secondly, the feasibility of using markedly different legacies of data in such a study is addressed, including an assessment of the factors influencing the crossmatching. From the latter, it is found that significant, spatially varying time shifts need to be, and are successfully, resolved through 3‐D warping. After the warping, the primary factors limiting the crossmatching appear to be residual local phase variations, possibly induced by the differing migration strategies, structure, reverberations and different coherencies of the volumes, caused by differences in acquisition‐structure azimuth and acquisition fold. Despite these differences, a time‐lapse amplitude signature is observed that is attributable to production. The character of the 4‐D amplitude anomalies may also indicate variations in stress sensitivity, e.g. because of zones of fracturing. Additionally, warping‐derived time attributes have been highlighted as a potential additional avenue for detection of pressure depletion in such reservoirs. Although the effects are subtle, they may indicate changes in stress/pressure in and around the reservoir because of production. However, to fully resolve the subtle time‐lapse effects in such a reservoir, the data differences need to be better addressed, which may be possible by full re‐processing and pre‐stack analysis, but more likely dedicated 4‐D acquisition would be required.  相似文献   

4.
We present an approach that creates the possibility of reservoir monitoring on a quasi‐continuous basis using surface seismic data. Current strategies and logistics for seismic data acquisition impose restrictions on the calendar‐time temporal resolution obtainable for a given surface‐seismic time‐lapse monitoring program. One factor that restricts the implementation of a quasi‐continuous monitoring program using conventional strategies is the time it takes to acquire a complete survey. Here quasi‐continuous monitoring describes the process of reservoir monitoring at short‐time intervals. Our approach circumvents the restriction by requiring only a subset of complete survey data each time an image of the reservoir is needed using surface seismic data. Ideally, the time interval between survey subset acquisitions should be short so that changes in the reservoir properties are small. The accumulated data acquired are used to estimate the unavailable data at the monitor survey time and the combined recorded and estimated data are used to produce an image of the subsurface for monitoring. We will illustrate the effectiveness of our approach using 2D and 3D synthetic seismic data and 3D field seismic data. We will explain the benefits and drawbacks of the proposed approach.  相似文献   

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

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

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.
Knowledge about the spatial distribution of the fracture density and the azimuthal fracture orientation can greatly help in optimizing production from fractured reservoirs. Frequency-dependent seismic velocity and attenuation anisotropy data contain information about the fractures present in the reservoir. In this study, we use the measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuths to infer information about the multiple fracture sets present in the reservoir. We consider a reservoir model with two sets of vertical fractures characterized by unknown azimuthal fracture orientations and fracture densities. Frequency-dependent seismic velocity and attenuation anisotropy data is computed using the effective viscoelastic stiffness tensor and solving the Christoffel equation. A Bayesian inversion method is then applied to measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuth to estimate the azimuthal fracture orientations and the fracture densities, as well as their uncertainties. Our numerical examples suggest that velocity anisotropy data alone cannot recover the unknown fracture parameters. However, an improved estimation of the unknown fracture parameters can be obtained by joint inversion of velocity and attenuation anisotropy data.  相似文献   

9.
Most sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis.  相似文献   

10.
Time‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (?30 m/s, or ?1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep).  相似文献   

11.
A method to provide an improved time‐lapse seismic attribute for dynamic interpretation is presented. This is based on the causal link between the time‐lapse seismic response and well production activity taken over time. The resultant image is obtained by computing correlation coefficients between sequences of time‐lapse seismic changes extracted over different time intervals from multiply repeated seismic and identical time sequences of cumulative fluid volumes produced or injected from the wells. Maps of these cross‐correlations show localized, spatially contiguous signals surrounding individual wells or a specific well group. These may be associated with connected regions around the selected well or well group. Application firstly to a synthetic data set reveals that hydraulic compartments may be delineated using this method. A second application to a field data set provides empirical evidence that a connected well‐centric fault block and active geobody can be detected. It is concluded that uniting well data and time‐lapse seismic using our proposed method delivers a new attribute for dynamic interpretation and potential updating of the model for the producing reservoir.  相似文献   

12.
The possibility of using 4D seismic data for monitoring pressure depletion in the low‐porosity, tight gas‐bearing Rotliegende sandstones of the UK Southern Gas Basin is investigated. The focus here is on whether fractures in the upper part of the reservoir, known to enhance productivity, can also enhance the time‐lapse seismic response. The study uses laboratory data to evaluate core‐plug stress sensitivity, published data for the stress behaviour of the fractures, followed by petro‐elastic and 4D seismic modelling of both the fractured and unfractured formation. The magnitude of the resultant 4D signatures suggests that production‐induced changes in the unfractured sands are unlikely to be observed except perhaps with highly repeatable time‐lapse surveys. On the other hand, the presence of fractures could render production effects visible in dedicated 4D acquisition or prestack parallel processed data. If present however, the signature will be sporadic, as fractures in the area are known to exist in clusters. The 4D signature may be enhanced further by certain classes of vertical geological variability and also areas of high reservoir pressure. The strongest evidence of depletion is expected to be time‐shifts seen at the base of the Rotliegende reservoir.  相似文献   

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

14.
The frequent time‐lapse observations from the life of field seismic system across the Valhall field provide a wealth of information. The responses from the production and injection wells can be observed through time‐shift and amplitude changes. These observations can be compared to modelled synthetic seismic responses from a reservoir simulation model of the Valhall Field. The observed differences between the observations and the modelling are used to update and improve the history match of the reservoir model. The uncertainty of the resulting model is reduced and a more confident prediction of future reservoir performance is provided. A workflow is presented to convert the reservoir model to a synthetic seismic response and compare the results to the observed time‐lapse responses for any time range and area of interest. Correlation based match quality factors are calculated to quantify the visual differences. This match quality factor allows us to quantitatively compare alternative reservoir models to help identify the parameters that best match the seismic observations. Three different case studies are shown where this workflow has helped to reduce the uncertainty range associated with specific reservoir parameters. By updating various reservoir model parameters we have been able to improve the match to the observations and thereby improve the overall reservoir model predictability. The examples show positive results in a range of different reservoir modelling issues, which indicates the flexibility of this workflow and the ability to have an impact in most reservoir modelling challenges.  相似文献   

15.
We present the analysis of a multi-azimuth vertical seismic profiling data set that has been acquired in a tight gas field with the objective of characterizing fracture distributions using seismic anisotropy. We investigate different measurements of anisotropy, which are shear-wave splitting, P-wave traveltime anisotropy and azimuthal amplitude variation with offset. We find that for our field case shear-wave splitting is the most robust measure of azimuthal anisotropy, which is clearly observed over two distinct intervals in the target. We compare the results of the vertical seismic profiling analysis with other borehole data from the same well. Cross-dipole sonic and Formation MicroImager data from the reservoir section suggest that no open fractures intersect the well or are present within half a metre of the borehole wall. Furthermore, a detailed dispersion analysis of the sonic scanner data provides no indication of stress-induced seismic anisotropy along the logged borehole section. We therefore explain the azimuthal anisotropy measured in the vertical seismic profiling data with a model that contains discrete fracture corridors, which do not intersect the well itself but lie within the vertical seismic profiling investigation radius. We show that such a model can reproduce some basic characteristics of azimuthal anisotropy observed in the vertical seismic profiling data. The model is also consistent with well test data that suggest the presence of a fracture corridor away from the well. With this study we demonstrate the necessity of integrating different data types that investigate different scales of rock volume and can provide complementary information for understanding the characteristics of fracture networks in the subsurface.  相似文献   

16.
Shear‐wave polarization and time delay are attributes commonly used for fracture detection and characterization. In time‐lapse analysis these parameters can be used as indicators of changes in the fracture orientation and density. Indeed, changes in fracture characteristics provide key information for increased reservoir characterization and exploitation. However, relative to the data uncertainty, is the comparison of these parameters over time statistically meaningful? We present the uncertainty in shear‐wave polarization and time delay as a function of acquisition uncertainties, such as receiver and source misorientation, miscoupling and band‐limited random noise. This study is applied to a time‐lapse borehole seismic survey, recorded in Vacuum Field, New Mexico. From the estimated uncertainties for each survey, the uncertainty in the difference between the two surveys is 31° for the shear‐wave polarization angle and 4 ms for the shear‐wave time delay. Any changes in these parameters greater than these error estimates can be interpreted with confidence. This analysis can be applied to any time‐lapse measurement to provide an interval of confidence in the interpretation of shear‐wave polarization angles and time splitting.  相似文献   

17.
Numerous examples of reservoir fields from continental and marine environments involve thin‐bedded geology, yet, the inter‐relationship between thin‐bedded geology, fluid flow and seismic wave propagation is poorly understood. In this paper, we explore the 4D seismic signature due to saturation changes of gas within thin layers, and address the challenge of identifying the relevant scales and properties, which correctly define the geology, fluid flow and seismic wave propagation in the field. Based on the study of an outcrop analogue for a thin‐bedded turbidite, we model the time‐lapse seismic response to fluid saturation changes for different levels of model scale, and explore discrepancies in quantitative seismic attributes caused by upscaling. Our model reflects the geological complexity associated with thin‐bedded turbidites, and its coupling to fluid flow, which in turn affects the gas saturation distribution in space, and its time‐lapse seismic imprint. Rock matrix and fluid properties are modelled after selected fields to reproduce representative field models with realistic impedance contrasts. In addition, seismic modelling includes multiples, in order to assess their contribution in seismic propagation through thin gas layers. Our results show that multiples could contribute significantly to the measured amplitudes in the case of thin‐bedded geology. This suggests that forward/inverse modelling involving the flow simulation and seismic domains used in time‐lapse seismic interpretation should account for thin layers, when these are present in the geological setting.  相似文献   

18.
Time‐lapse seismic surveying has become an accepted tool for reservoir monitoring applications, thus placing a high premium on data repeatability. One factor affecting data repeatability is the influence of the rough sea‐surface on the ghost reflection and the resulting seismic wavelets of the sources and receivers. During data analysis, the sea‐surface is normally assumed to be stationary and, indeed, to be flat. The non‐flatness of the sea‐surface introduces amplitude and phase perturbations to the source and receiver responses and these can affect the time‐lapse image. We simulated the influence of rough sea‐surfaces on seismic data acquisition. For a typical seismic line with a 48‐fold stack, a 2‐m significant‐wave‐height sea introduces RMS errors of about 5–10% into the stacked data. This level of error is probably not important for structural imaging but could be significant for time‐lapse surveying when the expected difference anomaly is small. The errors are distributed differently for sources and receivers because of the different ways they are towed. Furthermore, the source wavelet is determined by the sea shape at the moment the shot is fired, whereas the receiver wavelet is time‐varying because the sea moves significantly during the seismic record.  相似文献   

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

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

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

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