Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

Intra-survey reservoir fluctuations – implications for quantitative 4D seismic analysis

During the time taken for seismic data to be acquired, reservoir pressure may fluctuate as a consequence of field production and operational procedures and fluid fronts may move significantly. These variations prevent accurate quantitative measurement of the reservoir change using 4D seismic data. Modelling studies on the Norne field simulation model using acquisition data from ocean-bottom seismometer and towed streamer systems indicate that the pre-stack intra-survey reservoir fluctuations are important and cannot be neglected. Similarly, the time-lapse seismic image in the post-stack domain does not represent a difference between two states of the reservoir at a unique base and monitor time, but is a mixed version of reality that depends on the sequence and timing of seismic shooting. The outcome is a lack of accuracy in the measurement of reservoir changes using the resulting processed and stacked 4D seismic data. Even for perfect spatial repeatability between surveys, a spatially variant noise floor is still anticipated to remain. For our particular North Sea acquisition data, we find that towed streamer data are more affected than the ocean-bottom seismometer data. We think that this may be typical for towed streamers due to their restricted aperture compared to ocean-bottom seismometer acquisitions, even for a favourable time sequence of shooting and spatial repeatability. Importantly, the pressure signals on the near and far offset stacks commonly used in quantitative 4D seismic inversion are found to be inconsistent due to the acquisition timestamp. Saturation changes at the boundaries of fluid fronts appear to show a similar inconsistency across sub-stacks. We recommend that 4D data are shot in a consistent manner to optimize aerial time coverage, and that additionally, the timestamp of the acquisition should be used to optimize pre-stack quantitative reservoir analysis.     

Three-dimensional geophysical characterization of the La Rambla and Zafra de Záncara anticlines (Loranca Basin,Central Spain)

The Zafra de Záncara anticline (also known as the El Hito anticline), located in the Loranca Cenozoic Basin (part of the Tagus Basin, Central Spain), had been studied by several oil companies during the late 1960s and early 1970s. In 2009, within the ‘Plan for selection and characterization of suitable structures of CO₂ geological storage’, this anticline was selected as a potential CO₂ storage site. A preliminary three-dimensional geological model, based on five geological cross sections that were constrained with the interpretation of the available seismic profiles (that are rather old and do not have very good quality), was created. With the aim of improving the geological knowledge of the Zafra de Záncara anticline and helping to investigate the suitability of a nearby anticline, namely La Rambla, as another structural closure that might make it a possible CO₂ storage site, a local gravity survey (1 station every km²) was carried out in the area, seven new geological cross sections, based on these existing seismic profiles and field geology, were build, and a new three-dimensional geological model that included both anticlines, improved through three-dimensional stochastic gravity inversion, was constructed. The densities needed for the geological formations of the model come from the analysis of rock samples, logging data from El Hito-1 drillhole and petrophysical information from Instituto Geológico y Minero de España database. The inversion has improved the knowledge about the geometry of the anticlines’ traps and seals as well as the geometry of the basement relief and the structural relationship between basement and cover.     

Time‐lapse seismic imaging using regularized full‐waveform inversion with a prior model: which strategy?

Full‐waveform inversion is an appealing technique for time‐lapse imaging, especially when prior model information is included into the inversion workflow. Once the baseline reconstruction is achieved, several strategies can be used to assess the physical parameter changes, such as parallel difference (two separate inversions of baseline and monitor data sets), sequential difference (inversion of the monitor data set starting from the recovered baseline model) and double‐difference (inversion of the difference data starting from the recovered baseline model) strategies. Using synthetic Marmousi data sets, we investigate which strategy should be adopted to obtain more robust and more accurate time‐lapse velocity changes in noise‐free and noisy environments. This synthetic application demonstrates that the double‐difference strategy provides the more robust time‐lapse result. In addition, we propose a target‐oriented time‐lapse imaging using regularized full‐waveform inversion including a prior model and model weighting, if the prior information exists on the location of expected variations. This scheme applies strong prior model constraints outside of the expected areas of time‐lapse changes and relatively less prior constraints in the time‐lapse target zones. In application of this process to the Marmousi model data set, the local resolution analysis performed with spike tests shows that the target‐oriented inversion prevents the occurrence of artefacts outside the target areas, which could contaminate and compromise the reconstruction of the effective time‐lapse changes, especially when using the sequential difference strategy. In a strongly noisy case, the target‐oriented prior model weighting ensures the same behaviour for both time‐lapse strategies, the double‐difference and the sequential difference strategies and leads to a more robust reconstruction of the weak time‐lapse changes. The double‐difference strategy can deliver more accurate time‐lapse variation since it can focus to invert the difference data. However, the double‐difference strategy requires a preprocessing step on data sets such as time‐lapse binning to have a similar source/receiver location between two surveys, while the sequential difference needs less this requirement. If we have prior information about the area of changes, the target‐oriented sequential difference strategy can be an alternative and can provide the same robust result as the double‐difference strategy.