Site characterization for CO<Subscript>2</Subscript> geologic storage and vice versa: the Frio brine pilot,Texas, USA as a case study

Careful site characterization is critical for successful geologic storage of carbon dioxide (CO₂) because of the many physical and chemical processes impacting CO₂ movement and containment under field conditions. Traditional site characterization techniques such as geological mapping, geophysical imaging, well logging, core analyses, and hydraulic well testing provide the basis for judging whether or not a site is suitable for CO₂ storage. However, only through the injection and monitoring of CO₂ itself can the coupling between buoyancy flow, geologic heterogeneity, and history-dependent multi-phase flow effects be observed and quantified. CO₂ injection and monitoring can therefore provide a valuable addition to the site-characterization process. Additionally, careful monitoring and verification of CO₂ plume development during the early stages of commercial operation should be performed to assess storage potential and demonstrate permanence. The Frio brine pilot, a research project located in Dayton, Texas (USA) is used as a case study to illustrate the concept of an iterative sequence in which traditional site characterization is used to prepare for CO₂ injection and then CO₂ injection itself is used to further site-characterization efforts, constrain geologic storage potential, and validate understanding of geochemical and hydrological processes. At the Frio brine pilot, in addition to traditional site-characterization techniques, CO₂ movement in the subsurface is monitored by sampling fluid at an observation well, running CO₂-saturation-sensitive well logs periodically in both injection and observation wells, imaging with crosswell seismic in the plane between the injection and observation wells, and obtaining vertical seismic profiles to monitor the CO₂ plume as it migrates beyond the immediate vicinity of the wells. Numerical modeling plays a central role in integrating geological, geophysical, and hydrological field observations.     

Seal integrity and feasibility of CO<Subscript>2</Subscript> sequestration in the Teapot Dome EOR pilot: geomechanical site characterization

This paper reports a preliminary investigation of CO₂ sequestration and seal integrity at Teapot Dome oil field, Wyoming, USA, with the objective of predicting the potential risk of CO₂ leakage along reservoir-bounding faults. CO₂ injection into reservoirs creates anomalously high pore pressure at the top of the reservoir that could potentially hydraulically fracture the caprock or trigger slip on reservoir-bounding faults. The Tensleep Formation, a Pennsylvanian age eolian sandstone is evaluated as the target horizon for a pilot CO₂ EOR-carbon storage experiment, in a three-way closure trap against a bounding fault, termed the S1 fault. A preliminary geomechanical model of the Tensleep Formation has been developed to evaluate the potential for CO₂ injection inducing slip on the S1 fault and thus threatening seal integrity. Uncertainties in the stress tensor and fault geometry have been incorporated into the analysis using Monte Carlo simulation. The authors find that even the most pessimistic risk scenario would require ∼10 MPa of excess pressure to cause the S1 fault to reactivate and provide a potential leakage pathway. This would correspond to a CO₂ column height of ∼1,500 m, whereas the structural closure of the Tensleep Formation in the pilot injection area does not exceed 100 m. It is therefore apparent that CO₂ injection is not likely to compromise the S1 fault stability. Better constraint of the least principal stress is needed to establish a more reliable estimate of the maximum reservoir pressure required to hydrofracture the caprock.     

Dynamics of CO<Subscript>2</Subscript> fluxes and concentrations during a shallow subsurface CO<Subscript>2</Subscript> release

A field facility located in Bozeman, Montana provides the opportunity to test methods to detect, locate, and quantify potential CO₂ leakage from geologic storage sites. From 9 July to 7 August 2008, 0.3 t CO₂ day⁻¹ were injected from a 100-m long, ~2.5-m deep horizontal well. Repeated measurements of soil CO₂ fluxes on a grid characterized the spatio-temporal evolution of the surface leakage signal and quantified the surface leakage rate. Infrared CO₂ concentration sensors installed in the soil at 30 cm depth at 0–10 m from the well and at 4 cm above the ground at 0 and 5 m from the well recorded surface breakthrough of CO₂ leakage and migration of CO₂ leakage through the soil. Temporal variations in CO₂ concentrations were correlated with atmospheric and soil temperature, wind speed, atmospheric pressure, rainfall, and CO₂ injection rate.     

A stochastic inversion workflow for monitoring the distribution of CO<Subscript>2</Subscript> injected into deep saline aquifers

We present a two-step stochastic inversion approach for monitoring the distribution of CO₂ injected into deep saline aquifers for the typical scenario of one single injection well and a database comprising a common suite of well logs as well as time-lapse vertical seismic profiling (VSP) data. In the first step, we compute several sets of stochastic models of the elastic properties using conventional sequential Gaussian co-simulations (SGCS) representing the considered reservoir before CO₂ injection. All realizations within a set of models are then iteratively combined using a modified gradual deformation algorithm aiming at reducing the mismatch between the observed and simulated VSP data. In the second step, these optimal static models then serve as input for a history matching approach using the same modified gradual deformation algorithm for minimizing the mismatch between the observed and simulated VSP data following the injection of CO₂. At each gradual deformation step, the injection and migration of CO₂ is simulated and the corresponding seismic traces are computed and compared with the observed ones. The proposed stochastic inversion approach has been tested for a realistic, and arguably particularly challenging, synthetic case study mimicking the geological environment of a potential CO₂ injection site in the Cambrian-Ordivician sedimentary sequence of the St. Lawrence platform in Southern Québec. The results demonstrate that the proposed two-step reservoir characterization approach is capable of adequately resolving and monitoring the distribution of the injected CO₂. This finds its expression in optimized models of P- and S-wave velocities, density, and porosity, which, compared to conventional stochastic reservoir models, exhibit a significantly improved structural similarity with regard to the corresponding reference models. The proposed approach is therefore expected to allow for an optimal injection forecast by using a quantitative assimilation of all available data from the appraisal stage of a CO₂ injection site.     

Using hyperspectral plant signatures for CO<Subscript>2</Subscript> leak detection during the 2008 ZERT CO<Subscript>2</Subscript> sequestration field experiment in Bozeman,Montana

Hyperspectral plant signatures can be used as a short-term, as well as long-term (100-year timescale) monitoring technique to verify that CO₂ sequestration fields have not been compromised. An influx of CO₂ gas into the soil can stress vegetation, which causes changes in the visible to near-infrared reflectance spectral signature of the vegetation. For 29 days, beginning on July 9, 2008, pure carbon dioxide gas was released through a 100-m long horizontal injection well, at a flow rate of 300 kg day⁻¹. Spectral signatures were recorded almost daily from an unmown patch of plants over the injection with a “FieldSpec Pro” spectrometer by Analytical Spectral Devices, Inc. Measurements were taken both inside and outside of the CO₂ leak zone to normalize observations for other environmental factors affecting the plants. Four to five days after the injection began, stress was observed in the spectral signatures of plants within 1 m of the well. After approximately 10 days, moderate to high amounts of stress were measured out to 2.5 m from the well. This spatial distribution corresponded to areas of high CO₂ flux from the injection. Airborne hyperspectral imagery, acquired by Resonon, Inc. of Bozeman, MT using their hyperspectral camera, also showed the same pattern of plant stress. Spectral signatures of the plants were also compared to the CO₂ concentrations in the soil, which indicated that the lower limit of soil CO₂ needed to stress vegetation is between 4 and 8% by volume.     

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (P_n and S_n). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (P_cc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P_n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.     

中国大陆科学钻探孔区反射地震剖面的数值模拟与分析

本初步分析了苏鲁超高压变质带大陆科学钻探孔区的地震波速度结构，在此基础上对在东海孔区三维反射地震勘探过井孔的NNW-SSE主剖面和该周围地区的DH99二维反射地震剖面进行了数值模拟分析。结果表明，数值理论模拟结果与三维反射主剖面地震波到时吻合良好，模拟法能实现有效的时深转换，也表明数值模拟中采用的速度结构是合理的。东海钻井孔区实施反射地震勘测，所得地震记录剖面图是单一P波记录，在超高压变质岩区反射地震记录图中S波的影响甚微。DH99二维反射地震剖面的模拟结果表明东海孔区2000米深度以上的密度，波速结构不随深度加大而变化；而且苏鲁超高压变质岩区地表浅层的地震波速和密度较高，是高速高密地区。岩心结果反映的反射地震阻抗差大的地段，数值模拟的理论剖面上强反射面与孔区反射地震剖面图中的强反射体基本吻合．它意味着在非沉积层的结晶岩区，反射地震勘测手段是研究地下复杂结构的有效方法。数值模拟方法可以进一步提取地震反射记录信息，是利用地震反射方法研究地球构造的一种有效的方法。     

Crust–upper mantle seismic velocity structure across Southeastern China

One in-line wide-angle seismic profile was conducted in 1990 in the course of the Southeastern China Continental Dynamics project aimed at the study of the contact between the Cathaysia block and the Yangtze block. This 380-km-long profile extended in NW–SE direction from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. Five in-line shots were fired and recorded at seismic stations with spacing of about 3 km along the recording line. We have used two-dimensional ray tracing to model P- and S-wave arrivals and provide constraints on the velocity structure of the upper crust, middle crust, lower crust, Moho discontinuity, and the top part of the lithospheric mantle. P-wave velocity, S-wave velocity and V_P/V_S ratio are mapped. The crust is 36-km thick on average, albeit it gradually thins from the northwest end to the southeast end (offshore) of the profile. The average crustal velocity is 6.26 km/s for P-waves but 3.6 km/s for S-waves. A relatively narrow low-velocity layer of about 4 km of thickness, with P- and S-wave velocities of 6.2 km/s and 3.5 km/s, respectively, marks the bottom of the middle crust at a depth of 23-km northwest and 17-km southeast. At the crust–mantle transition, the P- and S-wave velocity change quickly from 7.4 to 7.8 km/s (northwest) and 8.0 to 8.2 km/s (southeast) and from 3.9 to 4.2 km/s (northwest) and 3.9 to 4.5 km/s (southeast), respectively. This result implies a lateral contrast in the upper mantle velocity along the 140 km sampled by the profile approximately. The average V_P/V_S ratio ranges from 1.68–1.8 for the upper crust to 1.75 for the middle and 1.75–1.85 for lower crust. With the interpretation of the wide-angle seismic data, Jiangshan–Shaoxin fault is considered as the boundary between the Yangtze and the Cathaysia block.     

CO<Subscript>2</Subscript> absorption by alkaline soils and its implication to the global carbon cycle

Motivated by the rapid increase in atmospheric CO₂ due to human activities since the Industrial Revolution, and the climate changes it produced, the world’s concerned scientific community has made a huge effort to investigate the global carbon cycle. However, the results reveal that the global CO₂ budget cannot be balanced, unless a “missing sink” is invoked. Although numerous studies claimed to find the “missing sink”, none of those claims has been widely accepted. This current study showed that alkaline soil on land are absorbing CO₂ at a rate of 0.3–3.0 μmol m⁻² s⁻¹ with an inorganic, non-biological process. The intensity of this CO₂ absorption is determined by the salinity, alkalinity, temperature and water content of the saline/alkaline soils, which are widely distributed on land. Further studies revealed that high salinity or alkalinity positively affected the CO₂ absorbing intensity, while high temperature and water content had a negative effect on the CO₂ absorbing intensity of these soils. This inorganic, non-biological process of CO₂ absorption by alkaline soils might have significant implications to the global carbon budget accounting.     

Lower crustal fluid distribution in the northeastern Japan arc revealed by high-resolution 3D seismic tomography

The Ou Backbone Range strikes northwards through the central northeastern Japan arc and is bounded on both sides by the active reverse Uwandaira and Sen'ya faults. We have applied a traveltime inversion method (seismic tomography) with spatial velocity correlation to active and passive seismic data in order to investigate a three-dimensional (3-D) velocity structure. The data set contains 33,993 P- and 18,483 S-wave arrivals from 706 natural sources and 40 blasts, as well as 2803 P-wave traveltime data from 10 explosions detonated during the 1997 controlled source experiment. The traveltime inversion reveals a zone beneath the Ou Backbone Range in which P-wave velocities (V_P) are approximately 6–8% lower than the average velocity at equivalent depths. The low V_P and a low V_P to S-wave velocity (V_S) ratio (V_P/V_S) of about 1.65 suggest the presence of aqueous fluids in the middle crust.     

Fault throw determination using 4 component VSP: Aigion fault (Greece) case study

The present paper presents the results of a technique based on Vertical Seismic Profiling (VSP) to quantify the offset of a fault in a context where seismic reflection profiles do not image any reflectors. The case study is located in Greece, in the Aigion area, on the south border of the Corinth Gulf. The Aigion fault, oriented East–West, and dipping at 60°N, has been intersected by the scientific well Aig-10 at 760 m in depth, but some uncertainty remained concerning its exact offset since the depth of the pre-rift sequence is open to debate due to the lack of subsurface data. The pre-rift consists of a Mesozoic low-porosity series (carbonates and radiolarite), while the syn-rift consists of poorly compacted conglomerates and turbidites.

Seismic diffractions on fault edges are expected to be present along the fault surface. We find that a few diffracted events are effectively recorded in P and S wave mode by the 4 components of the VSP survey, and critically refracted arrivals along the Aigion fault plane. These singular events have been used in order to refine the fault geometry and to determine its throw. Additionally, results from a 2D finite difference elastic seismic model of a single fault step geometry have been closely examined in order to illustrate the generation of the singular seismic events observed on the field data, such as P-wave and S-wave refracted and diffracted events, and support their interpretation with higher confidence. As a result, the seismic arrival patterns from full waveform seismic modelling confirm the characterisation of the main fault geometry (mainly its throw: about 200 m) in the well vicinity derived from the actual analysis of the 4 Component Aigion VSP data.     

First-principles simulation of high-pressure polymorphs in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript>

We have used density functional theory to investigate the stability of MgAl₂O₄ polymorphs under pressure. Our results can reasonably explain the transition sequence of MgAl₂O₄ polymorphs observed in previous experiments. The spinel phase (stable at ambient conditions) dissociates into periclase and corundum at 14 GPa. With increasing pressure, a phase change from the two oxides to a calcium-ferrite phase occurs, and finally transforms to a calcium-titanate phase at 68 GPa. The calcium-titanate phase is stable up to at least 150 GPa, and we did not observe a stability field for a hexagonal phase or periclase + Rh₂O₃(II)-type Al₂O₃. The bulk moduli of the phases calculated in this study are in good agreement with those measured in high-pressure experiments. Our results differ from those of a previous study using similar methods. We attribute this inconsistency to an incomplete optimization of a cell shape and ionic positions at high pressures in the previous calculations.     

P,S wave velocity model of the crust and upper most mantle of Albania region

This paper describes the one-dimensional (1D) velocity model computed by VELEST in the SEISAN seismic analysis system, inverting re-picked P-wave and S-wave arrival times recorded during 2002–2006 by the Albanian, Montenegro, Thessalonica and Macedonia seismic networks. The re-picked data yield P-wave and S-wave velocities proved to be more suitable compared to bulletin data for this detailed inversion study. Seismic phases recorded by the Albania seismic network and integrated with data from the Montenegro, Thessalonica and Macedonia networks are used to prepare the Albanian seismic bulletin. Earthquake hypocenters from the Albanian bulletins have also location errors that are negligible for civil protection purposes, large scale seismotectonic analyses and more accurate hypocentral determinations which are necessary for detailed seismotectonic and geodynamic studies.It was noted that the smoothness of the velocity variation increased with depth. A velocity of 5.5 km/s was calculated for the upper crust, 6.1 km/s was calculated for the middle crust and 6.9 km/s was computed for the lower crust. P wave velocity was 7.85 km/s at depth of 50 km and for the upper mantle it is 8.28 km/s. Using the improved velocity model, the earthquakes which occurred in Albania in the past 5 years were able to be relocated, achieving constrained hypocentral determinations for events in Albania. The interpretation of the 1 D velocity models infers interesting features of the deep structure of Albania. These results represent an important step towards more detailed seismotectonic analyses.     

Applications of Engineering Seismology for Site Characterization

We determined the seismic model of the soil column within a residential project site in Istanbul, Turkey. Specifically, we conducted a refraction seismic survey at 20 locations using a receiver spread with 484.5-Hz vertical geophones at 2-m intervals. We applied nonlinear tomography to first-arrival times to estimate the P-wave velocity-depth profiles and performed Rayleigh-wave inversion to estimate the S-wave velocity-depth profiles down to a depth of 30 m at each of the locations. We then combined the seismic velocities with the geotechnical borehole information regarding the lithology of the soil column and determined the site-specific geotechnical earthquake engineering parameters for the site. Specifically, we computed the maximum soil amplification ratio, maximum surface-bedrock acceleration ratio, depth interval of significant acceleration, maximum soil-rock response ratio, and design spectrum periods TA-TB. We conducted reflection seismic surveys along five line traverses with lengths between 150 and 300 m and delineated landslide failure surfaces within the site. We recorded shot gathers at 2-m intervals along each of the seismic line traverses using a receiver spread with 4 840-Hz vertical geophones at 2-m intervals. We applied nonlinear tomograpby to first-arrival times to estimate a P-wave velocity-depth model and analyzed the reflected waves to obtain a seismic image of the deep near-surface along each of the line traverses.     

Gas–water–rock interactions in sedimentary basins: CO2 sequestration in the Frio Formation, Texas, USA

To investigate the potential for the geologic storage of CO₂ in saline sedimentary aquifers, 1600 ton of CO₂ were injected at 1500 m depth into a 24-m sandstone section of the Frio Formation — a regional reservoir in the US Gulf Coast. Fluid samples obtained from the injection and observation wells before, during and after CO₂ injection show a Na–Ca–Cl type brine with 93,000 mg/L TDS and near saturation of CH₄ at reservoir conditions. As injected CO₂ gas reached the observation well, results showed sharp drops in pH (6.5 to 5.7), pronounced increases in alkalinity (100 to 3000 mg/L as HCO₃) and Fe (30 to 1100 mg/L), and significant shifts in the isotopic compositions of H₂O and DIC. Geochemical modeling indicates that brine pH would have dropped lower, but for buffering by dissolution of calcite and Fe oxyhydroxides. Post-injection results show the brine gradually returning to its pre-injection composition.     

Comparison of P- and S-wave velocity profiles obtained from surface seismic refraction/reflection and downhole data

High-resolution seismic-reflection/refraction data were acquired on the ground surface at six locations to compare with near-surface seismic-velocity downhole measurements. Measurement sites were in Seattle, WA, the San Francisco Bay Area, CA, and the San Fernando Valley, CA. We quantitatively compared the data in terms of the average shear-wave velocity to 30-m depth (V_s30), and by the ratio of the relative site amplification produced by the velocity profiles of each data type over a specified set of quarter-wavelength frequencies. In terms of V_s30, similar values were determined from the two methods. There is <15% difference at four of the six sites. The V_s30 values at the other two sites differ by 21% and 48%. The relative site amplification factors differ generally by less than 10% for both P- and S-wave velocities. We also found that S-wave reflections and first-arrival phase delays are essential for identifying velocity inversions. The results suggest that seismic reflection/refraction data are a fast, non-invasive, and less expensive alternative to downhole data for determining V_s30. In addition, we emphasize that some P- and S-wave reflection travel times can directly indicate the frequencies of potentially damaging earthquake site resonances. A strong correlation between the simple S-wave first-arrival travel time/apparent velocity on the ground surface at 100 m offset from the seismic source and the V_s30 value for that site is an additional unique feature of the reflection/refraction data that could greatly simplify V_s30 determinations.     

CO<Subscript>2</Subscript> source-sink matching in the lower 48 United States,with examples from the Texas Gulf Coast and Permian Basin

Documenting geographic distribution and spatial linkages between CO₂ sources and potential sinks in areas with significant levels of CO₂ emissions is important when considering carbon-management strategies such as geologic sequestration or enhanced oil recovery (EOR). For example, the US Gulf Coast overlies a thick succession (>6,000 m [>20,000 ft]) of highly porous and permeable sandstone formations separated by thick, regionally extensive shale aquitards. The Gulf Coast and Permian Basin also have a large potential for EOR, in which CO₂ injected into suitable oil reservoirs could be followed by long-term storage of CO₂ in nonproductive formations below reservoir intervals. For example, >6 billion barrels (Bbbl) of oil from 182 large reservoirs is technically recoverable in the Permian Basin as a result of miscible-CO₂ flooding_. The Gulf Coast also contains an additional 4.5 Bbbl of oil that could be produced by using miscible CO₂. Although the CO₂ pipeline infrastructure is well-developed in the Permian Basin, east Texas and the Texas Gulf Coast may have a greater long-term potential for deep, permanent storage of CO₂ because of thick brine-bearing formations near both major subsurface and point sources of CO₂.     

Relationship between soil CO<Subscript>2</Subscript> flux and volcanic tremor at Mt. Etna: Implications for magma dynamics

Large variations of the CO₂ flux through the soil were observed between November 2002 and January 2006 at Mt. Etna volcano. In many cases, the CO₂ flux was strongly influenced by changes in air temperature and atmospheric pressure. A new filtering method was then developed to remove the atmospheric influences on soil CO₂ flux and, at the same time, to highlight the variations strictly related to volcanic activity. Successively, the CO₂ corrected data were quantitatively compared with the spectral amplitude of the volcanic tremor by cross correlation function, cross-wavelet spectrum and wavelet coherence. These analyses suggested that the soil CO₂ flux variations preceded those of volcanic tremor by about 50 days. Given that volcanic tremor is linked to the shallow (a few kilometer) magma dynamics and soil CO₂ flux related to the deeper (~12 km b.s.l.) magma dynamics, the “delayed similarity” between the CO₂ flux and the volcanic tremor amplitude was used to assess the average speed in the magma uprising into the crust, as about 170–260 m per day. Finally, the large amount of CO₂ released before the onset of the 2004–2005 eruption indicated a deep ingression of new magma, which might have triggered such an eruption.