ALBERT

Description: Reliable quantification of carbon dioxide ([Formula: see text]) properties and saturation is crucial in the monitoring of [Formula: see text] underground storage projects. We have focused on quantitative seismic characterization of [Formula: see text] at the Sleipner storage pilot site. We evaluate a methodology combining high-resolution seismic waveform tomography, with uncertainty quantification and rock physics inversion. We use full-waveform inversion (FWI) to provide high-resolution estimates of P-wave velocity [Formula: see text] and perform an evaluation of the reliability of the derived model based on posterior covariance matrix analysis. To get realistic estimates of [Formula: see text] saturation, we implement advanced rock physics models taking into account effective fluid phase theory and patchy saturation. We determine through sensitivity tests that the estimation of [Formula: see text] saturation is possible even when using only the P-wave velocity as input. After a characterization of rock frame properties based on log data prior to the [Formula: see text] injection at Sleipner, we apply our two-step methodology. The FWI result provides clear indications of the injected [Formula: see text] plume being observed as low-velocity zones corresponding to thin [Formula: see text] filled layers. Several tests, varying the rock physics model and [Formula: see text] properties, are then performed to estimate [Formula: see text] saturation. The results suggest saturations reaching 30%–35% in the thin sand layers and up to 75% when patchy mixing is considered. We have carried out a joint estimation of saturation with distribution type and, even if the inversion is not well-constrained due to limited input data, we conclude that the [Formula: see text] has an intermediate pattern between uniform and patchy mixing, which leads to saturation levels of approximately [Formula: see text]. It is worth noting that the 2D section used in this work is located 533 m east of the injection point. We also conclude that the joint estimation of [Formula: see text] properties with saturation is not crucial and consequently that knowing the pressure and temperature state of the reservoir does not prevent reliable estimation of [Formula: see text] saturation.

Description: Risk assessment of CO2 storage requires the use of geophysical monitoring techniques to quantify changes in selected reservoir properties such as CO2 saturation, pore pressure and porosity. Conformance monitoring and associated decision-making rest upon the quantified properties derived from geophysical data, with uncertainty assessment. A general framework combining seismic and controlled source electromagnetic inversions with rock physics inversion is proposed with fully Bayesian formulations for proper quantification of uncertainty. The Bayesian rock physics inversion rests upon two stages. First, a search stage consists in exploring the model space and deriving models with associated probability density function (PDF). Second, an appraisal or importance sampling stage is used as a "correction" step to ensure that the full model space is explored and that the estimated posterior PDF can be used to derive quantities like marginal probability densities. Both steps are based on the neighbourhood algorithm. The approach does not require any linearization of the rock physics model or assumption about the model parameters distribution. After describing the CO2 storage context, the available data at the Sleipner field before and after CO2 injection (baseline and monitor), and the rock physics models, we perform an extended sensitivity study. We show that prior information is crucial, especially in the monitor case. We demonstrate that joint inversion of seismic and CSEM data is also key to quantify CO2 saturations properly. We finally apply the full inversion strategy to real data from Sleipner. We obtain rock frame moduli, porosity, saturation and patchiness exponent distributions and associated uncertainties along a 1D profile before and after injection. The results are consistent with geology knowledge and reservoir simulations, i.e., that the CO2 saturations are larger under the caprock confirming the CO2 upward migration by buoyancy effect. The estimates of patchiness exponent have a larger uncertainty, suggesting semi-patchy mixing behaviour.