首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Carbon capture and storage is a viable greenhouse gas mitigation technology and the Sleipner CO2 sequestration site in the North Sea is an excellent example. Storage of CO2 at the Sleipner site requires monitoring over large areas, which can successfully be accomplished with time lapse seismic imaging. One of the main goals of CO2 storage monitoring is to be able to estimate the volume of the stored CO2 in the reservoir. This requires a parametrization of the subsurface as exact as possible. Here we use elastic 2D time‐domain full waveform inversion in a time lapse manner to obtain a P‐wave velocity constrain directly in the depth domain for a base line survey in 1994 and two post‐injection surveys in 1999 and 2006. By relating velocity change to free CO2 saturation, using a rock physics model, we find that at the considered location the aquifer may have been fully saturated in some places in 1999 and 2006.  相似文献   

2.
Seismic monitoring of reservoir and overburden performance during subsurface CO2 storage plays a key role in ensuring efficiency and safety. Proper interpretation of monitoring data requires knowledge about the rock physical phenomena occurring in the subsurface formations. This work focuses on rock stiffness and elastic velocity changes of a shale overburden formation caused by both reservoir inflation induced stress changes and leakage of CO2 into the overburden. In laboratory experiments, Pierre shale I core plugs were loaded along the stress path representative for the in situ stress changes experienced by caprock during reservoir inflation. Tests were carried out in a triaxial compaction cell combining three measurement techniques and permitting for determination of (i) ultrasonic velocities, (ii) quasistatic rock deformations, and (iii) dynamic elastic stiffness at seismic frequencies within a single test, which allowed to quantify effects of seismic dispersion. In addition, fluid substitution effects connected with possible CO2 leakage into the caprock formation were modelled by the modified anisotropic Gassmann model. Results of this work indicate that (i) stress sensitivity of Pierre shale I is frequency dependent; (ii) reservoir inflation leads to the increase of the overburden Young's modulus and Poisson's ratio; (iii) in situ stress changes mostly affect the P‐wave velocities; (iv) small leakage of the CO2 into the overburden may lead to the velocity changes, which are comparable with one associated with geomechanical influence; (v) non‐elastic effects increase stress sensitivity of an acoustic waves; (iv) and both geomechanical and fluid substitution effects would create significant time shifts, which should be detectable by time‐lapse seismic.  相似文献   

3.
CO2 has been injected into the saline aquifer Utsira Fm at the Sleipner field since 1996. In order to monitor the movement of the CO2 in the sub‐surface, the seventh seismic monitor survey was acquired in 2010, with dual sensor streamers which enabled optimal towing depths compared to previous surveys. We here report both on the time‐lapse observations and on the improved resolution compared to the conventional streamer surveys. This study shows that the CO2 is still contained in the subsurface, with no indications of leakage. The time‐lapse repeatability of the dual sensor streamer data versus conventional data is sufficient for interpreting the time‐lapse effects of the CO2 at Sleipner, and the higher resolution of the 2010 survey has enabled a refinement of the interpretation of nine CO2 saturated layers with improved thickness estimates of the layers. In particular we have estimated the thickness of the uppermost CO2 layer based on an analysis of amplitude strength together with time‐separation of top and base of this layer and found the maximum thickness to be 11 m. This refined interpretation gives a good base line for future time‐lapse surveys at the Sleipner CO2 injection site.  相似文献   

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

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

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

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

8.
4D seismic is widely used to remotely monitor fluid movement in subsurface reservoirs. This technique is especially effective offshore where high survey repeatability can be achieved. It comes as no surprise that the first 4D seismic that successfully monitored the CO2 sequestration process was recorded offshore in the Sleipner field, North Sea. In the case of land projects, poor repeatability of the land seismic data due to low S/N ratio often obscures the time‐lapse seismic signal. Hence for a successful on shore monitoring program improving seismic repeatability is essential. Stage 2 of the CO2CRC Otway project involves an injection of a small amount (around 15,000 tonnes) of CO2/CH4 gas mixture into a saline aquifer at a depth of approximately 1.5 km. Previous studies at this site showed that seismic repeatability is relatively low due to variations in weather conditions, near surface geology and farming activities. In order to improve time‐lapse seismic monitoring capabilities, a permanent receiver array can be utilised to improve signal to noise ratio and hence repeatability. A small‐scale trial of such an array was conducted at the Otway site in June 2012. A set of 25 geophones was installed in 3 m deep boreholes in parallel to the same number of surface geophones. In addition, four geophones were placed into boreholes of 1–12 m depth. In order to assess the gain in the signal‐to‐noise ratio and repeatability, both active and passive seismic surveys were carried out. The surveys were conducted in relatively poor weather conditions, with rain, strong wind and thunderstorms. With such an amplified background noise level, we found that the noise level for buried geophones is on average 20 dB lower compared to the surface geophones. The levels of repeatability for borehole geophones estimated around direct wave, reflected wave and ground roll are twice as high as for the surface geophones. Both borehole and surface geophones produce the best repeatability in the 30–90 Hz frequency range. The influence of burying depth on S/N ratio and repeatability shows that significant improvement in repeatability can be reached at a depth of 3 m. The level of repeatability remains relatively constant between 3 and 12 m depths.  相似文献   

9.
Conventional seismic data are band limited and therefore, provide limited geological information. Every method that can push the limits is desirable for seismic data analysis. Recently, time‐frequency decomposition methods are being used to quickly extract geological information from seismic data and, especially, for revealing frequency‐dependent amplitude anomalies. Higher frequency resolution at lower frequencies and higher temporal resolution at higher frequencies are the objectives for different time‐frequency decomposition methods. Continuous wavelet transform techniques, which are the same as narrow‐band spectral analysis methods, provide frequency spectra with high temporal resolution without the windowing process associated with other techniques. Therefore, this technique can be used for analysing geological information associated with low and high frequencies that normally cannot be observed in conventional seismic data. In particular, the continuous wavelet transform is being used to detect thin sand bodies and also as a direct hydrocarbon indicator. This paper presents an application of the continuous wavelet transform method for the mapping of potential channel deposits, as well as remnant natural gas detection by mapping low‐frequency anomalies associated with the gas. The study was carried out at the experimental CO2 storage site at Ketzin, Germany (CO2SINK). Given that reservoir heterogeneity and faulting will have significant impact on the movement and storage of the injected CO2, our results are encouraging for monitoring the migration of CO2 at the site. Our study confirms the efficiency of the continuous wavelet transform decomposition method for the detection of frequency‐dependent anomalies that may be due to gas migration during and after the injection phase and in this way, it can be used for real‐time monitoring of the injected CO2 from both surface and borehole seismics.  相似文献   

10.
In the Norwegian North Sea, the Sleipner field produces gas with a high CO2 content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO2) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO2 plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO2 over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO2 plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO2 injection and the 2006 vintage carried out after a CO2 injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO2 in a reservoir. With the post‐stack inversion workflow applied to CO2 storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.  相似文献   

11.
More than 50 000 tons of CO2 have been injected at Ketzin into the Stuttgart Formation, a saline aquifer, at approximately 620 m depth, as of summer 2011. We present here results from the 1st repeat 3D seismic survey that was performed at the site in autumn 2009, after about 22 000 tons of CO2 had been injected. We show here that rather complex time‐lapse signatures of this CO2 can be clearly observed within a radius of about 300 m from the injection well. The highly irregular amplitude response within this radius is attributed to the heterogeneity of the injection reservoir. Time delays to a reflection below the injection level are also observed. Petrophysical measurements on core samples and geophysical logging of CO2 saturation levels allow an estimate of the total amount of CO2 visible in the seismic data to be made. These estimates are somewhat lower than the actual amount of CO2 injected at the time of the survey and they are dependent upon the choice of a number of parameters. In spite of some uncertainty, the close agreement between the amount injected and the amount observed is encouraging for quantitative monitoring of a CO2 storage site using seismic methods.  相似文献   

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

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

14.
The sequestration of CO2 in subsurface reservoirs constitutes an immediate counter‐measure to reduce anthropogenic emissions of CO2, now recognized by international scientific panels to be the single most critical factor driving the observed global climatic warming. To ensure and verify the safe geological containment of CO2 underground, monitoring of the CO2 site is critical. In the high Arctic, environmental considerations are paramount and human impact through, for instance, active seismic surveys, has to be minimized. Efficient seismic modelling is a powerful tool to test the detectability and imaging capability prior to acquisition and thus improve the characterization of CO2 storage sites, taking both geological setting and seismic acquisition set‐up into account. The unique method presented here avoids the costly generation of large synthetic data sets by employing point spread functions to directly generate pre‐stack depth‐migrated seismic images. We test both a local‐target approach using an analytical filter assuming an average velocity and a full‐field approach accounting for the spatial variability of point spread functions. We assume a hypothetical CO2 plume emplaced in a sloping aquifer inspired by the conditions found at the University of Svalbard CO2 lab close to Longyearbyen, Svalbard, Norway, constituting an unconventional reservoir–cap rock system. Using the local‐target approach, we find that even the low‐to‐moderate values of porosity (5%–18%) measured in the reservoir should be sufficient to induce significant change in seismic response when CO2 is injected. The sensitivity of the seismic response to changes in CO2 saturation, however, is limited once a relatively low saturation threshold of 5% is exceeded. Depending on the illumination angle provided by the seismic survey, the quality of the images of five hypothetical CO2 plumes of varying volume differs depending on the steepness of their flanks. When comparing the resolution of two orthogonal 2D surveys to a 3D survey, we discover that the images of the 2D surveys contain significant artefacts, the CO2‐brine contact is misplaced and an additional reflector is introduced due to the projection of the point spread function of the unresolvable plane onto the imaging plane. All of these could easily lead to a misinterpretation of the behaviour of the injected CO2. Our workflow allows for testing the influence of geological heterogeneities in the target aquifer (igneous intrusions, faults, pervasive fracture networks) by utilizing increasingly complex and more realistic geological models as input as more information on the subsurface becomes available.  相似文献   

15.
Natural depletion of subsurface petroleum liquids releases energy in the form of heat. The rate of natural source zone depletion (NSZD) can be derived from subsurface temperature data. An energy balance is performed to resolve NSZD‐generated energy in terms of W/m2. Biodegradation rates are resolved by dividing the NSZD energy by the heat of reaction in joules/mol. Required temperature data are collected using data loggers, wireless connections, and automated data storage and analysis. Continuous thermal resolution of monthly NSZD rates at a field site indicates that apparent monthly NSZD rates vary through time, ranging from 10,000 to 77,000 L/ha/year. Temporal variations in observed apparent NSZD rates are attributed to processes governing the conversion of CH4 to CO2, as opposed to the actual rates of NSZD. Given a year or more of continuous NSZD rate data, it is anticipated that positive and negative biases in apparent NSZD rates will average out, and averaged apparent NSZD rates will converge to true NSZD rates. An 8.4% difference between average apparent NSZD rates over a 31‐month period using the thermal monitoring method and seven rounds of CO2 efflux measurements using CO2 traps supports the validity of both CO2 trap and thermal monitoring methods. A promising aspect of thermal monitoring methods is that continuous data provide a rigorous approach to resolving the true mean NSZD rates as compared to temporally sparse CO2 trap NSZD rate measurements. Overall, a vision is advanced of real‐time sensor‐based groundwater monitoring that can provide better data at lower costs and with greater safety, security, and sustainability.  相似文献   

16.
A series of time‐lapse seismic cross‐well and single‐well experiments were conducted in a diatomite reservoir to monitor the injection of CO2 into a hydrofracture zone, based on P‐ and S‐wave data. A high‐frequency piezo‐electric P‐wave source and an orbital‐vibrator S‐wave source were used to generate waves that were recorded by hydrophones as well as 3‐component geophones. During the first phase the set of seismic experiments was conducted after the injection of water into the hydrofractured zone. The set of seismic experiments was repeated after a time period of seven months during which CO2 was injected into the hydrofractured zone. The questions to be answered ranged from the detectability of the geological structure in the diatomic reservoir to the detectability of CO2 within the hydrofracture. Furthermore, it was intended to determine which experiment (cross‐well or single‐well) is best suited to resolve these features. During the pre‐injection experiment, the P‐wave velocities exhibited relatively low values between 1700 and 1900 m/s, which decreased to 1600–1800 m/s during the post‐injection phase (?5%). The analysis of the pre‐injection S‐wave data revealed slow S‐wave velocities between 600 and 800 m/s, while the post‐injection data revealed velocities between 500 and 700 m/s (?6%). These velocity estimates produced high Poisson's ratios between 0.36 and 0.46 for this highly porous (～50%) material. Differencing post‐ and pre‐injection data revealed an increase in Poisson's ratio of up to 5%. Both velocity and Poisson's ratio estimates indicate the dissolution of CO2 in the liquid phase of the reservoir accompanied by an increase in pore pressure. The single‐well data supported the findings of the cross‐well experiments. P‐ and S‐wave velocities as well as Poisson's ratios were comparable to the estimates of the cross‐well data. The cross‐well experiment did not detect the presence of the hydrofracture but appeared to be sensitive to overall changes in the reservoir and possibly the presence of a fault. In contrast, the single‐well reflection data revealed an arrival that could indicate the presence of the hydrofracture between the source and receiver wells, while it did not detect the presence of the fault, possibly due to out‐of‐plane reflections.  相似文献   

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

18.
The injection of CO2 at the Ketzin pilot CO2 storage site started in June 2008 and ended in August 2013. During the 62 months of injection, a total amount of about 67 kt of CO2 was injected into a saline aquifer. A third repeat three‐dimensional seismic survey, serving as the first post‐injection survey, was acquired in 2015, aiming to investigate the recent movement of the injected CO2. Consistent with the previous two time‐lapse surveys, a predominantly west–northwest migration of the gaseous CO2 plume in the up‐dip direction within the reservoir is inferred in this first post‐injection survey. No systematic anomalies are detected through the reservoir overburden. The extent of the CO2 plume west of the injection site is almost identical to that found in the 2012 second repeat survey (after injection of 61 kt); however, there is a significant decrease in its size east of the injection site. Assessment of the CO2 plume distribution suggests that the decrease in the size of the anomaly may be due to multiple factors, such as limited vertical resolution, CO2 dissolution, and CO2 migration into thin layers, in addition to the effects of ambient noise. Four‐dimensional seismic modelling based on dynamic flow simulations indicates that a dynamic balance between the newly injected CO2 after the second repeat survey and the CO2 migrating into thin layers and being dissolved was reached by the time of the first post‐injection survey. In view of the significant uncertainties in CO2 mass estimation, both patchy and non‐patchy saturation models for the Ketzin site were taken into consideration.  相似文献   

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

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

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

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