CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion     

Hong Yan  Bastien Dupuy  Anouar Romdhane  Børge Arntsen 《Geophysical Prospecting》2019,67(4):1055-1071

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.  相似文献   

Landslide Susceptibility Mapping Using GIS-based Fuzzy Logic and the Analytical Hierarchical Processes Approach: A Case Study in Constantine (North-East Algeria)     

Abdı  Amina  Bouamrane  Ali  Karech  Toufik  Dahri  Noura  Kaouachi  Anouar 《Geotechnical and Geological Engineering》2021,39(8):5675-5691

The main purpose of this study was to compare and evaluate the performance of two multicriteria models for landslide susceptibility assessment in Constantine, north-east of Algeria. The landslide susceptibility maps were produced using the analytic hierarchy process (AHP) and Fuzzy AHP (FAHP) via twelve landslides conditioning factors, including the slope gradient, lithology, land cover, distance from drainage network, distance from the roads, distance from faults, topographic wetness index, stream power index, slope curvature, Normalized Difference Vegetation Index, slope aspect and elevation. In this study, the mentioned models were used to derive the weighting value of the conditioning factors. For the validation process of these models, the receiver operating characteristic analysis, and the area under the curve (AUC) were applied by comparing the obtained results to The landslide inventory map which prepared using the archives of scientific publications, reports of local authorities, and field survey as well as analyzing satellite imagery. According to the AUC values, the FAHP model had the highest value (0.908) followed by the AHP model (0.777). As a result, the FAHP model is more consistent and accurate than the AHP in this case study. The outcome of this paper may be useful for landslide susceptibility assessment and land use management.

  相似文献   

Paleoenvironmental conditions and strontium isotope stratigraphy in the Paleogene Gafsa Basin (Tunisia) deduced from geochemical analyses of phosphatic fossils     

László Kocsis  Anouar Ounis  Fredj Chaabani  Neili Mohamed Salah 《International Journal of Earth Sciences》2013,102(4):1111-1129

Fossil shark teeth and coprolites from three major phosphorite occurrences in the Gafsa Basin (southwestern Tunisia) were investigated for their geochemical compositions to improve local stratigraphy and to better assess paleoenvironmental conditions. ⁸⁷Sr/⁸⁶Sr isotope ratios of shark teeth from the Early Maastrichtian El Haria Formation and from the Early Eocene Métlaoui s.s. Formation yielded Sr isotope ages of 68 ± 1 and 47.9 ± 1.3 Ma, respectively, which accord with the expected stratigraphic positions of these sediments. Conversely, shark teeth from the Paleocene–Eocene Chouabine Formation have large variation in Sr isotope ratios even within individual layers. After statistical treatment and then elimination of certain outlier samples, three age-models are proposed and discussed. The most reasonable solution includes three subsequent Sr ages of 61.8 ± 2.2 Ma, 57.2 ± 1.8 and 54.6 ± 1.6 for layer IX, layers VIII–V and layers IV–0, respectively. Three scenarios are discussed for explanation of the presence of the outliers: (1) diagenesis, (2) re-working and (3) locally controlled seawater Sr isotope ratio. The most plausible account for the higher ⁸⁷Sr/⁸⁶Sr ratios relative to the global ocean in some fossils is enhanced intrabasinal re-working due to low sea level. Conversely, the sample with lower ⁸⁷Sr/⁸⁶Sr than the global seawater may link to diagenesis or to seawater influenced by weathering of Late Cretaceous marine carbonates, which latter is supported by model calculation as well. The ε_Nd values of these fossils are very similar to those reported for Paleogene and Late Cretaceous Tethyan seawater and are compatible with the above interpretations. The relatively low oxygen isotope values in shark teeth from the topmost phosphate bed of the Chouabine Formation, together with the Sr isotope results, point toward recovering better connections with the open sea. These δ¹⁸O data reflect elevated ambient temperature, which may link to the Early Eocene Climatic Optimum.  相似文献   

Offset dependence of overburden time-shifts from ultrasonic data     

Audun Bakk  Rune M. Holt  Andreas Bauer  Bastien Dupuy  Anouar Romdhane 《Geophysical Prospecting》2020,68(6):1847-1863

Depletion or injection into a reservoir implies stress changes and strains in the reservoir and its surroundings. This may lead to measurable time-shifts for seismic waves propagating in the subsurface. To better understand the offset dependence of time-shifts in the overburden, we have systematically quantified the time-shifts of three different overburden shales in controlled laboratory tests. These experiments may be viewed as an analogue to the time-shifts recorded from seismic field surveys. For a range of different stress paths, defined as the ratio between the horizontal and the vertical stress changes, the changes of the P-wave velocities in different directions were measured such that the offset dependence of time-shifts for different stress paths could be studied. The time-shifts are stress path dependent, which is particularly pronounced at large offsets. For all stress paths, the time-shifts exhibit a linearly decreasing trend with increasing offset, that is, a negative offset gradient. At zero offset, for which the ray path is normal to the bedding, the time-shifts are similar for all investigated stress paths. The isotropic stress path is associated with the smallest offset gradient of the time-shifts. In contrast, the constant-mean-stress path shows the largest gradient with a flip in the polarity of the time-shifts for the largest offsets. The separate contributions from the strain and velocity changes to the time-shifts were also quantified. The time-shifts for the isotropic stress path are dominated by the contribution from velocity changes at all offsets. In contrast, the strain contributes significantly to the time-shifts at small offsets for the constant-mean-stress path. This shows that the offset dependence in pre-stack seismic data may be a key to understand the changes of subsurface stresses, pore pressure and strain upon depletion or injection. To utilize this knowledge from laboratory experiments, calibrated rock physics models and correlations are needed to constrain the seismic time-shifts and to obtain an adequately updated geological model reflecting the true anisotropic nature of the subsurface. This may have important implications for improved recovery and safety, particularly in mature fields.  相似文献