Search

Abstract: Induced seismicity is one of the main risks for gigaton-scale geological storage of carbon dioxide (CO2). Thus, passive seismic monitoring is often recommended as a necessary component of the monitoring systems for CO2 storage projects, with a particular forcus on risk mitigation. We present the first field study, CO2CRC Otway Project Stage 3 (Victoria, Australia), where distributed acoustic sensing (DAS) enabled high-precision tracking of the induced seismicity triggered by a small CO2 injection and also informed the reservoir models. In 610 days of passive seismic monitoring of the Stage 3 injection, we detected 17 microseismic events (maximum moment magnitude Mw 0.1) using five deep boreholes equipped with enhanced-sensitivity optical fiber. The DAS array has sensitivity sufficient for detection and location of induced events with Mw∼−2 in a monitoring borehole located up to 1500 m away. Thanks to the dense spatial sampling by the DAS, we were able to estimate the focal mechanisms for events with Mw>−1.5; although the monitoring boreholes provided very limited angular coverage. The main cluster of the events has the same location and source mechanism as the one triggered by the previous CO2 injection at the Otway Project site, Stage 2C. Surprizingly, the Stage 2C and Stage 3 events closely followed the actual movement of the CO2 saturation plume front (not the pressure front), as observed using controlled-source reflection seismic images. The nature of the plume-fault interaction remains unclear, but some alteration of the fault gouge by CO2 might be responsible for the faults’ reactivation by the pressure perturbation. Importantly, the seismogenic fault could not be identified in the seismic images and was only revealed by DAS observations, which also demonstrated the signature of fluid–rock interaction, that may control the CO2 flow.

Seismic velocities and elastic moduli of rocks are known to vary significantly with applied stress, which indicates that these materials exhibit nonlinear elasticity. Monochromatic waves in nonlinear elastic media are known to generate higher harmonics and combinational frequencies. Such effects have the potential to be used for broadening the frequency band of seismic sources, characterization of the subsurface, and safety monitoring of civil engineering infrastructure. However, knowledge on nonlinear seismic effects is still scarce, which impedes the development of their practical applications. To explore the potential of nonlinear seismology, we performed three experiments: two in the field and one in the laboratory. The first field experiment used two vibroseis sources generating signals with two different monochromatic frequencies. The second field experiment used a surface orbital vibrator with two eccentric motors working at different frequencies. In both experiments, the generated wavefield was recorded in a borehole using a fiber-optic distributed acoustic sensing cable. Both experiments showed combinational frequencies, harmonics, and other intermodulation products of the fundamental frequencies both on the surface and at depth. Laboratory experiments replicated the setup of the field test with vibroseis sources and showed similar nonlinear combinations of fundamental frequencies. Amplitudes of the nonlinear signals observed in the laboratory showed variation with the saturating fluid. These results confirm that nonlinear components of the wavefield propagate as body waves, are likely to generate within rock formations, and can be potentially used for reservoir fluid characterization.

Active time-lapse seismic is widely employed for monitoring CO2 geosequestration due to its ability to track the distribution of fluids in space and time. However, standard 4D seismic monitoring suffers from several challenges, including high cost, disruption to other land uses, and, consequently, relatively large intervals between monitor surveys. Some of these challenges can be mitigated using permanently installed sources and receivers. Such an approach was tested at the CO2CRC Otway site by continuous offset VSP monitoring of 15,000 t of supercritical CO2 injected into an aquifer 1,500 m deep with nine permanent seismic sources (surface orbital vibrators or SOVs) and five downhole fibre-optic receivers. This continuous monitoring is complemented by multi-well 4D VSP using a mobile vibroseis source and the same DAS receivers, which included one baseline and two monitor surveys after injection of 4,000 and 12,000 t of CO2. The continuous DAS-SOV monitoring detected an abrupt increase of travel times below the injection interval on the second day of injection (after injection of 300 t of CO2) and tracked the growth of the areal CO2 plume by mapping changes of reflection amplitudes. The plume is also detected by time-lapse changes of reflection amplitudes in multi-well 4D VSPs. The plume images obtained from continuous offset VSP and 4D VSP are broadly consistent with each other but with some differences due to differences in illumination, lateral variations of velocities and seismic anisotropy. These differences also serve as a measure of uncertainty of 4D VSP images.

A permanent automated continuous seismic CO2 geosequestration monitoring system for was installed at CO2CRC Otway Project site (Victoria, Australia) in early 2020. The system is composed of five deviated ∼1600 m deep wells equipped with distributed acoustic sensing (DAS) acting as seismic receivers and nine seismic orbital vibrators (SOV) as seismic sources. DAS recording is performed continuously by three iDASv3 units. Each SOV operates for 2.5 h at a time, and hence all SOVs operating sequentially (during daytime only) produce in a single vintage every two days. Each vintage consists of 45 offset VSP transects covering predicted CO2 plume migration paths over ∼0.7 km2 area. An automated data processing implemented on-site reduces data size from ∼1.3 TB/day to ∼500 MB/day with the results transmitted to the office daily. The repeatability analysis based on pre-injection data (acquired from May to October 2020 before the injection start in December 2020) shows that variability of SOV performance is the main source of non-repeatability while borehole measurements are stable. An SOV waveform could reach NRMS value from 20 to 100 % within a few days. However, deconvolution of the seismograms with the waveform of the direct wave reduces the repeatability to within 10–15 % NRMS.

CO2 sequestration projects benefit from quantitative assessment of saturation distribution and plume extent for field development and leakage prevention. In this work, we carry out quantitative analysis of time-lapse seismic by using rock physics and seismic modelling tools. We investigate the suitability of Gassmann's equation for a CO2 sequestration project with 1600 tons of CO2 injected into high-porosity, brine-saturated sandstone. We analyze the observed time delays and amplitude changes in a time-lapse vertical seismic profile dataset. Both reflected and transmitted waves are analyzed qualitatively and quantitatively. To interpret the changes obtained from the vertical seismic profile, we perform a 2.5D elastic, finite-difference modelling study. The results show a P-wave velocity reduction of 750 m/s in the proximity of the injection well evident by the first arrivals (travel-time delays and amplitude change) and reflected wave amplitude changes. These results do not match with our rock physics model using Gassmann's equation predictions even when taking uncertainty in CO2 saturation and grain properties into account. We find that time-lapse vertical seismic profile data integrated with other information (e.g., core and well log) can be used to constrain the velocity–saturation relation and verify the applicability of theoretical models such as Gassmann's equation with considerable certainty. The study shows that possible nonelastic factors are in play after CO2 injection (e.g., CO2–brine–rock interaction and pressure effect) as Gassmann's equation underestimated the velocity reduction in comparison with field data for all three sets of time-lapse vertical seismic profile attributes. Our work shows the importance of data integration to validate the applicability of theoretical models such as Gassmann's equation for quantitative analysis of time-lapse seismic data.