Breakthroughs in seismic and borehole characterization of Basalt sequestration targets

AbstractMafic continental flood basalts form a globally important, but under-characterized CO2 sequestration target. The Columbia River Basalt Group (CRBG) in the northwestern U.S. is up to 5 km thick and covers over 168,000 km2. In India, flood basalts are 3 km thick and cover greater than 500,000 km2. Laboratory experiments demonstrate that the CRBG and other basalts react with formation water and super critical (sc) CO2 to precipitate carbonates, thus adding a potential mineral trapping mechanism to the standard trapping mechanisms of most other types of CO2 sequestration reservoirs.Brecciated tops of individual basalt flows in the CRBG form regional aquifers that locally have greater than 30% porosity and three Darcies of permeability. Porous flow tops are potential sites for sequestration of gigatons of sc CO2 in areas where the basalts contain unpotable water and are at depths greater than 800 m. In this paper we report on the U.S. DOE Big Sky Regional Carbon Sequestration Partnership surface seismic and borehole geophysical characterization that supports a field test of capacity, integrity, and geochemical reactivity of CRBG reservoirs in eastern Washington, U.S.A.Traditional surface seismic methods have had little success in imaging basalt features in on-shore areas where the basalt is thinly covered by sediment. Processing of the experimental 6.5 km, 5 line 3C seismic swath included constructing an elastic wavefield model, identifying and separating seismic wave modes, and processing the swath as a single 2D line. Important findings include: (1) a wide variety of shear wave energy modes swamp the P-wave seismic records; (2) except at very short geophone offsets, ground roll overprints P-wave signal; and (3) because of extreme velocity contrasts, P-wave events are refracted at incidence angles greater than 7–15 degrees. Subsequent removal of S-wave and other noise during processing resulted in tremendous improvement in image quality.The application of wireline logging to onshore basalts is underexploited. Full waveform sonic logs and resistivity-based image logs acquired in the 1250 m basalt pilot borehole provide powerful tools for evaluating geomechanics and lithofacies. The azimuth of the fast shear wave is parallel to SH and records the changes through geologic time in basalt flow and tectonic stress tensors. Combined with image log data, azimuthal S-wave data provide a borehole technique for assessing basalt emplacement and cooling history that is related to the development of reservoirs and seals, as well as the orientation of tectonic stresses and fracture systems that could affect CO2 transport or containment. Reservoir and seal properties are controlled by basalt lithofacies, and rescaled P- and S-wave slowness curves, integrated with image logs, provide a tool for improved recognition of subsurface lithofacies

Sub-basalt decoupled walkaway VSP imaging

Walkaway vertical seismic profile (VSP) data acquired in basalt-covered areas can be used to improve knowledge of the sub-basalt structure. A synthetic example and a case study from the North Atlantic (UK) show that elastic two-way downward-continuation migration combined with the stationary-phase principle is well suited to the processing of VSP data. Vector data are processed using decoupled elastic migration algorithms in both isotropic and anisotropic media. To illustrate the value of decoupled imaging equations, conventional PP imaging is carried out on the enhanced VSP data and compared with the decoupled scheme. Decoupled vector migration operates directly on the displacement vector, and uses various wave modes. Downgoing waves are migrated to image basalt lava flows and measure their anisotropy. Upgoing waves are used for high-resolution sub-basalt imaging

Joint inversion of seismic and geoelectric data recorded in an underground coal mine

Until the present time the 'rock-coal-rock' layer sequence and offsets in coal-seams in underground coal mines have been detected with the aid of seismic waves and geoelectric measurements. In order to determine the geometrical and petrophysical parameters of the coal-seam situation, the data recorded using seismic and geoelectric methods have been inverted independently. In consequence, the inversion of partially inaccurate data resulted in a certain degree of ambiguity. This paper presents the first results of a joint inversion scheme to process underground vertical seismic profiling data, geolectric resistivity and resistance data. The joint inversion algorithm makes use of the damped least-squares method and its weighted version to solve the linearized set of equations for the seismic and geolectric unknowns. In order to estimate the accuracy and reliability of the derived geometrical and petrophysical layer parameters, both a model covariance matrix and a correlation matrix are calculated. The weighted least-squares algorithm is based on the method of most frequent values (MFV). The weight factors depend on the difference between measured data and those calculated by an iteration process. The joint inversion algorithm is tested by means of synthetic data. Compared to the damped least-squares algorithm, the MFV inversion leads to smaller estimation errors as well as lower sensitivities due to the choice of the initial model. It is shown that, compared to an independent inversion, the correlation between the model parameters is definitely reduced, while the accuracy of the parameter estimation is appreciably increased by the joint inversion process. Thus the ambiguity is significantly reduced. Finally, the joint inversion algorithm the MFV method is applied to underground field data. The model parameters can be derived with a sufficient degree of accuracy, even in the case of noisy data