• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.10-18
• /
• 2006

We measured the changes in P-wave velocity that occur when injecting $CO_2$ in gaseous, liquid, and supercritical phases into water-saturated anisotropic sandstones. P-wave velocities were measured in two cylindrical samples of Tako Sandstone, drilled along directions normal and parallel to the bedding plane, using a piezo-electric transducer array system. The velocity changes caused by $CO_2$ injection are typically -6% on average, with maximum values about -16% for the case of supercritical $CO_2$ injection. P-wave velocity tomograms obtained by the differential arrival-time method clearly show that $CO_2$ migration behaviour is more complex when $CO_2$ flows normal to the bedding plane than when it flows parallel to bedding. We also found that the differences in P-wave velocity images were associated both with the $CO_2$ phases and with heterogeneity of pore distribution in the rocks. Seismic images showed that the highest velocity reduction occurred for supercritical $CO_2$ injection, compared with gaseous or liquid $CO_$ injection. This result may justify the use of the seismic method for $CO_2$ monitoring in geological sequestration.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.58-68
• /
• 2011

As part of basic studies of monitoring carbon dioxide ($CO_2$) storage using electrical and seismic surveys, laboratory experiments have been conducted to measure resistivity and P-wave velocity changes due to the injection of $CO_2$ into water-saturated sandstone. The rock sample used is a cylinder of Berea sandstone. $CO_2$ was injected under supercritical conditions (10 MPa, $40^{\circ}C$). The experimental results show that resistivity increases monotonously throughout the injection period, while P-wave velocity and amplitude decrease drastically due to the supercritical $CO_2$ injection. A reconstructed P-wave velocity tomogram clearly images $CO_2$ migration in the sandstone sample. Both resistivity and seismic velocity are useful for monitoring $CO_2$ behaviour. P-wave velocity, however, is less sensitive than resistivity when the $CO_2$ saturation is greater than ~20%. The result indicates that the saturation estimation from resistivity can effectively complement the difficulty of $CO_2$ saturation estimations from seismic velocity variations. By combining resistivity and seismic velocity we were able to estimate $CO_2$ saturation distribution and the injected $CO_2$ behaviour in our sample.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.85-98
• /
• 2009

The tau-p inversion algorithm is widely employed to generate starting models with most computer programs, which implement refraction tomography. This algorithm emphasises the vertical resolution of many layers, and as a result, it frequently fails to detect even large lateral variations in seismic velocities, such as the decreases which are indicative of shear zones. This study demonstrates the failure of the tau-p inversion algorithm to detect or define a major shear zone which is 50m or 10 stations wide. Furthermore, the majority of refraction tomography programs parameterise the seismic velocities within each layer with vertical velocity gradients. By contrast, the Generalized Reciprocal Method (GRM) inversion algorithms emphasise the lateral resolution of individual layers. This study demonstrates the successful detection and definition of the 50m wide shear zone with the GRM inversion algorithms. The existence of the shear zone is confirmed by a 2D analysis of the head wave amplitudes and by numerous closely spaced orthogonal seismic profiles carried out as part of a later 3D refraction investigation. Furthermore, an analysis of the shot record amplitudes indicates that a reversal in the seismic velocities, rather than vertical velocity gradients, occurs in the weathered layers. The major conclusion reached in this study is that while all seismic refraction operations should aim to provide as accurate depth estimates as is practical, those which emphasise the lateral resolution of individual layers generate more useful results for geotechnical and environmental applications. The advantages of the improved lateral resolution are obtained with 2D traverses in which the structural features can be recognised from the magnitudes of the variations in the seismic velocities. Furthermore, the spatial patterns obtained with 3D investigations facilitate the recognition of structural features such as faults which do not display any intrinsic variation or 'signature' in seismic velocities.