Time-lapse elastic full waveform inversion is used to monitor the spatio-temporal evolution of the CO2 plume during and after supercritical CO2 injection based on a series of time-lapse (repeated) cross-well seismic monitoring datasets obtained at the Nagaoka Carbon Capture and Storage (CCS) site in Japan. The full waveform inversion method successfully estimates the time-lapse velocity decrease of up to 30% within a thin 12 m layer, which is consistent with the magnitude and thickness of the well-log measurements. After the second monitoring survey, the velocity decrease becomes stable and gradually extends down dip along pre-existing geological structures. The full waveform inversion results starkly contrast with the previous estimates based on traveltime tomography. The previous traveltime tomography applications only used the traveltime-delays and resulted in low resolution with few percentage change which was not adequate to correctly resolve CO2 injection changes. The datasets pose significant challenges due to background noise, tube waves, apparent non-isotropic source radiation patterns, apparent reservoir velocity anisotropy and missing key acquisition parameters such as the number of stacks per shot point. To overcome these obstacles, we meticulously perform careful data preprocessing integrating both the body waves and tube waves. We develop waveform-based source mechanism estimation to represent non-isotropic source excitation, and then conduct forward modeling studies to constrain the anisotropy model.

Our limited understanding of mineral reactive surface area contributes to significant uncertainties in quantitative simulations of reactive chemical transport in subsurface processes. Continuum formulations for reactive transport typically use a number of different approximations for reactive surface area, including geometric, specific, and effective surface area. In this study, reactive surface area estimates are developed and evaluated for their ability to predict dissolution rates in a well-stirred flow-through reactor experiment using disaggregated samples from the Nagaoka pilot CO2 injection site (Japan). The disaggregated samples are reacted with CO2 acidified synthetic brine under conditions approximating the field conditions and the evolution of solute concentrations in the reactor effluent is tracked over time. The experiments, carried out in fluid-dominated conditions at a pH of 3.2 for 650 h, resulted in substantial dissolution of the sample and release of a disproportionately large fraction of the divalent cations. Traditional reactive surface area estimation methods, including an adjusted geometric surface area and a BET-based surface area, are compared to a newly developed image-based method. Continuum reactive transport modeling is used to determine which of the reactive surface area models provides the best match with the effluent chemistry from the well-stirred reactor. The modeling incorporates laboratory derived mineral dissolution rates reported in the literature and the initial modal mineralogy of the Nagaoka sediment was determined from scanning electron microscopy (SEM) characterization. The closest match with the observed steady-state effluent concentrations was obtained using specific surface area estimates from the image-based approach supplemented by literature-derived BET measurements. To capture the evolving effluent chemistry, particularly over the first 300 h of the experiment, it was also necessary to account for the grain size distribution in the sediment and the presence of a highly reactive volcanic glass phase that shows preferential cation leaching.

The rates of mineral dissolution reactions in porous media are difficult to predict, in part because of a lack of understanding of mineral reactive surface area in natural porous media. Common estimates of mineral reactive surface area used in reactive transport models for porous media are typically ad hoc and often based on average grain size, increased to account for surface roughness or decreased by several orders of magnitude to account for reduced surface reactivity of field as opposed to laboratory samples. In this study, accessible mineral surface areas are determined for a sample from the reservoir formation at the Nagaoka pilot CO2 injection site (Japan) using a multi-scale image analysis based on synchrotron X-ray microCT, SEM QEMSCAN, XRD, SANS, and FIB-SEM. This analysis not only accounts for accessibility of mineral surfaces to macro-pores, but also accessibility through connected micro-pores in smectite, the most abundant clay mineral in this sample. While the imaging analysis reveals that most of the micro- and macro-pores are well connected, some pore regions are unconnected and thus inaccessible to fluid flow and diffusion. To evaluate whether mineral accessible surface area accurately reflects reactive surface area a flow-through core experiment is performed and modeled at the continuum scale. The core experiment is performed under conditions replicating the pilot site and the evolution of effluent solutes in the aqueous phase is tracked. Various reactive surface area models are evaluated for their ability to capture the observed effluent chemistry, beginning with parameter values determined as a best fit to a disaggregated sediment experiment (Beckingham et al., 2016) described previously. Simulations that assume that all mineral surfaces are accessible (as in the disaggregated sediment experiment) over-predict the observed mineral reaction rates, suggesting that a reduction of RSA by a factor of 10–20 is required to match the core flood experimental data. While the fit of the effluent chemistry (and inferred mineral dissolution rates) greatly improve when the pore-accessible mineral surface areas are used, it was also necessary to include highly reactive glass phases to match the experimental observations, in agreement with conclusions from the disaggregated sediment experiment. It is hypothesized here that the 10–20 reduction in reactive surface areas based on the limited pore accessibility of reactive phases in core flood experiment may be reasonable for poorly sorted and cemented sediments like those at the Nagaoka site, although this reflects pore rather than larger scale heterogeneity.