首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Legacy streamer data and newer 3D ocean‐bottom‐cable data are cross‐matched and analysed for time‐lapse analysis of geomechanical changes due to production in the Valhall Field. The issues relating to time‐lapse analysis using two such distinctly different data sets are addressed to provide an optimal cross‐matching workflow that includes 3D warping. Additionally an assessment of the differences between the imaging using single‐azimuth streamer and multi‐azimuth ocean‐bottom‐cable data is provided. The 3D warping utilized in the cross‐matching procedure is sensitive to acquisition and processing differences but is also found to provide valuable insight into the geometrical changes that occur in the subsurface due to production. As such, this work also provides a demonstration of the use of high‐resolution 3D interpreted warping to resolve the 3D heterogeneity of the compaction and subsidence. This is an important tool for Valhall, and possibly other fields, where compaction and subsidence (and monitoring thereof) are key factors in the reservoir management since the predominant observed production‐induced changes are compaction of the soft, high‐porosity chalk reservoir, due to pore‐pressure reduction, and the resultant overburden subsidence. Such reservoir compaction could have significant implications for production by changing permeabilities and production rates. Furthermore the subsidence effects could impact upon subsea installations and well‐bore stability. Geomechanical studies that have previously been used to model such subsidence and compaction are only constrained by observed surface displacements and measured reservoir pressure changes, with the geological overburden being largely neglected. The approaches suggested herein provide the potential for monitoring and assessment in three dimensions, including the probable heterogeneity and shearing, that is needed for full understanding of reservoir compaction and the resultant effects on the overburden to, for example, mitigate well‐bore failures.  相似文献   

2.
This article addresses the question whether time‐lapse seismic reflection techniques can be used to follow and quantify the effects of solution salt mining. Specifically, the production of magnesium salts as mined in the north of the Netherlands is considered. The use of seismic time‐lapse techniques to follow such a production has not previously been investigated. For hydrocarbon production and CO2 storage, time‐lapse seismics are used to look at reservoir changes mainly caused by pressure and saturation changes in large reservoirs, while for solution mining salt is produced from caverns with a limited lateral extent, with much smaller production volumes and a fluid (brine) replacing a solid (magnesium salt). In our approach we start from the present situation of the mine and then study three different production scenarios, representing salt production both in vertical and lateral directions of the mine. The present situation and future scenarios have been transformed into subsurface models that were input to an elastic finite‐difference scheme to create synthetic seismic data. These data have been analysed and processed up to migrated seismic images, such that time‐lapse analyses of intermediate and final results could be done. From the analyses, it is found that both vertical and lateral production is visible well above the detection threshold in difference data, both at pre‐imaging and post‐imaging stages. In quantitative terms, an additional production of the mine of 6 m causes time‐shifts in the order of 2 ms (pre‐imaging) and 4 ms (post‐imaging) and amplitude changes of above 20% in the imaged sections. A laterally oriented production causes even larger amplitude changes at the edge of the cavern due to replacement of solid magnesium salt with brine introducing a large seismic contrast. Overall, our pre‐imaging and post‐imaging time‐lapse analysis indicates that the effects of solution salt mining can be observed and quantified on seismic data. The effects seem large enough to be observable in real seismic data containing noise.  相似文献   

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

4.
Two formulae are developed for estimating horizontal permeability directly from maps of 4D seismic signatures. The choice of the formula used depends on whether the seismic is dominated by changes of pressure or saturation. However, pressure derived from time‐lapse seismic, or seismic amplitudes controlled predominantly by pressure are to be preferred for optimal ‘illumination’ of the reservoir. The permeability is predicted to be dependent on porosity but weighted by a 4D term related to the magnitude and spatial gradient of the 4D signature. Tests performed on model‐based synthetic seismic data affirm the validity and accuracy of this approach. Application to field data from the UK continental shelf reveals a large‐scale permeability variation similar to the existing simulation model, but with additional fine‐scale detail. The technique thus has the potential of providing extra information with which to update the simulation model. The resultant permeability estimates have been successfully ground‐truthed against the results of two well tests. As non‐repeatable noise in the time‐lapse seismic data diminishes with improved 4D‐related acquisition, it will become increasingly possible to make robust permeability estimates using this approach.  相似文献   

5.
Hydrocarbon production and fluid injection affect the level of subsurface stress and physical properties of the subsurface, and can cause reservoir‐related issues, such as compaction and subsidence. Monitoring of oil and gas reservoirs is therefore crucial. Time‐lapse seismic is used to monitor reservoirs and provide evidence of saturation and pressure changes within the reservoir. However, relative to background velocities and reflector depths, the time‐lapse changes in velocity and geomechanical properties are typically small between consecutive surveys. These changes can be measured by using apparent displacement between migrated images obtained from recorded data of multiple time‐lapse surveys. Apparent displacement measurements by using the classical cross‐correlation method are poorly resolved. Here, we propose the use of a phase‐correlation method, which has been developed in satellite imaging for sub‐pixel registration of the images, to overcome the limitations of cross‐correlation. Phase correlation provides both vertical and horizontal displacements with a much better resolution. After testing the method on synthetic data, we apply it to a real dataset from the Norne oil field and show that the phase‐correlation method can indeed provide better resolution.  相似文献   

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

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

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

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

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

14.
Pre‐stack seismic data are indicative of subsurface elastic properties within the amplitude versus offset characteristic and can be used to detect elastic rock property changes caused by injection. We perform time‐lapse pre‐stack 3‐D seismic data analysis for monitoring sequestration at Cranfield. The time‐lapse amplitude differences of Cranfield datasets are found entangled with time‐shifts. To disentangle these two characters, we apply a local‐correlation‐based warping method to register the time‐lapse pre‐stack datasets, which can effectively separate the time‐shift from the time‐lapse seismic amplitude difference without changing the original amplitudes. We demonstrate the effectiveness of our registration method by evaluating the inverted elastic properties. These inverted time‐lapse elastic properties can be reliably used for monitoring plumes.  相似文献   

15.
  相似文献   

16.
With the increasing use of permanently installed seismic installations, many of the issues in time‐lapse seismic caused by the lack of repeatability can be reduced. However, a number of parameters still influence the degree of reliability of 4D seismic data. In this paper, the specific impact of seawater velocity variations on time‐lapse repeatability is investigated in a synthetic study. A zero‐lag time‐lapse seabed experiment with no change in the subsurface but with velocity changes in the water column is simulated. The velocity model in the water column is constant for the baseline survey while the model for the repeat survey is heterogeneous, designed from sea salinity and temperature measurements in the West of Shetlands. The difference section shows up to 80% of residual amplitude, which highlights the poor repeatability. A new dynamic correction which removes the effect of seawater velocity variations specifically for permanent installations is developed. When applied to the synthetic data, it reduces the difference residual amplitude to about 3%. This technique shows substantial improvement in repeatability beyond conventional time‐lapse cross‐equalization.  相似文献   

17.
The rough sea surface causes perturbations in the seismic data that can be significant for time‐lapse studies. The perturbations arise because the reflection response of the non‐flat sea perturbs the seismic wavelet. In order to remove these perturbations from the received seismic data, special deconvolution methods can be used, but these methods require, as input, the time varying wave elevation above each hydrophone in the streamer. In addition, the vertical displacement of the streamer itself must also be known at the position of each hydrophone and at all times. This information is not available in conventional seismic acquisition. However, it can be obtained from the hydrophone measurements provided that the hydrophones are recorded individually (not grouped), that the recording bandwidth is extended down to 0.05 Hz and that data are recorded without gaps between the shot records. The sea surface elevation, and also the wave‐induced vertical displacement of the streamer, can be determined from the time‐varying pressure that the sea waves cause in the hydrophone measurements. When this was done experimentally, using a single sensor seismic streamer without a conventional low cut filter, the wave induced pressure variations were easily detected. The inversion of these experimental data gives results for the sea surface elevation that are consistent with the weather and sea state at the time of acquisition. A high tension approximation allows a simplified solution of the equations that does not demand a knowledge of the streamer tension. However, best results at the tail end of the streamer are obtained using the general equation.  相似文献   

18.
Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO2 sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO2 sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO2 sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO2 sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.  相似文献   

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

20.
A calendar time interpolation method for 2D seismic amplitude maps, done in two steps, is presented. The contour interpolation part is formulated as a quadratic programming problem, whereas the amplitude value interpolation is based on a conditional probability formulation. The method is applied on field data from the Sleipner CO2 storage project. The output is a continuous image (movie) of the CO2 plume. Besides visualization, the output can be used to better couple 4D seismic to other types of data acquired. The interpolation uncertainty increases with the time gap between consecutive seismic surveys and is estimated by leaving a survey out (blind test). Errors from such tests can be used to identify problems in understanding the flow and possibly improve the interpolation scheme for a given case. Field‐life cost of various acquisition systems and repeat frequencies are linked to the time‐lapse interpolation errors. The error in interpolated amplitudes increased by 3%‐4% per year of interpolation gap for the Sleipner case. Interpolation can never fully replace measurements.  相似文献   

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

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