Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Caracterización de parámetros petrofísicos relacionados con la capacidad de almacenamiento de CO2 en el acuífero salino profundo de Hontomín (Burgos)

La conciencia de la crisis de la modernidad -que comienza ya a finales del siglo XIX- ha cobrado más experiencia debido al conocimiento de los límites del desarrollo económico, ya que como parecía razonable pensar, también los recursos naturales son finitos. En 1972, el Club de Roma analizó las distintas opciones disponibles para conseguir armonizar el desarrollo sostenible y las limitaciones medioambientales. Fue en 1987 cuando la Comisión Mundial para el Medio Ambiente y el Desarrollo de la ONU definía por primera vez el concepto de desarrollo sostenible. Definición que posteriormente fue incorporada en todos los programas de la ONU y sirvió de eje, por ejemplo, a la Cumbre de la Tierra celebrada en Río de Janeiro en 1992. Parece evidente que satisfacer la demanda energética, fundamentalmente desde la Revolución Industrial en el s XIX, trajo consigo un creciente uso de los combustibles fósiles, con la consiguiente emisión de los gases de efecto invernadero (GEI) y el aumento de la temperatura global media terrestre. Esta temperatura se incrementó en los últimos cien años en una media de 0.74ºC. La mayor parte del incremento observado desde la mitad del siglo XX en esta temperatura media se debe, con una probabilidad de al menos el 90%, al aumento observado en los GEI antropogénicos, siendo uno de ellos el CO2 que proviene de la transformación del carbono de los combustibles fósiles durante su combustión. Ante el creciente uso de los combustibles fósiles, los proyectos CAC, proyectos de captura, transporte y almacenamiento, se presentan como una contribución al desarrollo sostenible ya que se trata de una tecnología que permite mitigar el cambio climático. Para valorar si la tecnología CAC es sostenible, habrá que comprobar si existe o no capacidad para almacenar el CO2 en una cantidad mayor a la de producción y durante el tiempo necesario que impone la evolución de la concentración de CO2 en la atmósfera para mantenerla por debajo de las 450ppmv (concentración de CO2 que propone el Panel Intergubernamental para el Cambio Climático). El desarrollo de los proyectos CAC completos pasa por la necesaria selección de adecuados almacenes de CO2 que sean capaces de soportar los efectos de las presiones de inyección, así como asegurar la capacidad de dichos almacenes y la estanqueidad del CO2 en los mismos. La caracterización geológica de un acuífero susceptible de ser almacén de CO2 debe conducir a determinar las propiedades que dicho almacén posee para asegurar un volumen adecuado de almacenamiento, una inyectabilidad del CO2 en el mismo a un ritmo adecuado y la estanqueidad del CO2 en dicho acuífero a largo plazo. El presente trabajo pretende estudiar los parámetros que tienen influencia en el cálculo de la capacidad del almacén, para lo que en primer lugar se ha desarrollado la tecnología necesaria para llevar a cabo la investigación mediante ensayos de laboratorio. Así, se ha desarrollado una patente, "ATAP, equipo para ensayos petrofísicos (P201231913)", con la que se ha llevado a cabo la parte experimental de este trabajo para la caracterización de los parámetros que tienen influencia en el cálculo de la capacidad del almacén. Una vez desarrollada la tecnología, se aborda el estudio de los distintos parámetros que tienen influencia en la capacidad del almacén realizando ensayos con ATAP. Estos ensayos definen el volumen del almacenamiento, llegándose a la conclusión de que en la determinación de este volumen, juegan un papel importante el alcance de los mecanismos trampa, físicos o químicos, del CO2 en el almacén. Ensayos que definen la capacidad del almacén de "aceptar" o "rechazar" el CO2 inyectado, la inyectabilidad, y por último, ensayos encaminados a determinar posibles fugas que se pueden dar a través de los pozos de inyección, definidos estos como caminos preferenciales de fugas en un almacén subterráneo de CO2. Queda de este modo caracterizada la estanqueidad del CO2 en el acuífero a largo plazo y su influencia obvia en la determinación de la capacidad del almacén. Unido al propósito de la estimación de la capacidad del almacén, se encuentra el propósito de asegurar la estanqueidad de dichos almacenes en el tiempo, y adelantarse a la evolución de la pluma de CO2 en el interior de dichos almacenes. Para cumplir este propósito, se ha desarrollado un modelo dinámico a escala de laboratorio, mediante el programa ECLIPSE 300, con el fin de establecer una metodología para el cálculo de la capacidad estimada del almacén, así como el estudio de la evolución de la pluma de CO2 dentro del acuífero a lo largo del tiempo, partiendo de los resultados obtenidos en los ensayos realizados en ATAP y con la modelización de la probeta de roca almacén empleada en dichos ensayos. Presentamos por tanto un trabajo que establece las bases metodológicas para el estudio de la influencia de distintos parámetros petrofísicos en el cálculo de la capacidad del almacén unidos al desarrollo tecnológico de ATAP y su utilización para la determinación de dichos parámetros aplicables a cada acuífero concreto de estudio. ABSTRACT The crisis of modernity –which begins at the end of 19th Century- has been more important due to the knowledge of the limits of economic development, since it appeared to be thought reasonable, the natural resources are finite. In 1972, The Club of Rome analyzed the different options available in order to harmonize the sustainability and the environment development. It was in 1987 when The Global Commission on The Environment and the Development of UN, defined for the first time the concept of Sustainable Development. This definition that was fully incorporated in all the UN programs and it was useful as an axis, for example, in La Cumbre de la Tierra summit in Río de Janeiro in 1992. It seems obvious to satisfy energetic demand, basically after The Industrial Revolution in 19th Century, which represented an increasing use of fossil fuels, therefore greenhouse gases emission and the increasing of global average temperature. This temperature increased in the last 100 years up to 0.74ºC. The major part of the temperature increase is due to the increase observed in Greenhouse gases with human origin, at least with 90% of probability. The most important gas is the CO2 because of its quantity. In the face of the increasing use of fossil fuels, the CCS projects, Carbon Capture and Storage projects, appear as a contribution of sustainable development since it is a technology for avoiding the climate change. In order to evaluate if CCS technology is sustainable, it will be necessary to prove if the capacity for CO2 storage is available or not in a quantity greater than the production one and during the time necessary to keep the CO2 concentration in the atmosphere lower than 450ppmv (concentration imposed by IPCC). The development of full CCS projects goes through the selection of good CO2 storages that are able to support the effects of pressure injection, and assure the capacity of such storages and the watertightness of CO2. The geological characterization of the aquifer that could be potential CO2 storage should lead to determine the properties that such storage has in order to assure the adequate storage volume, the CO2 injectivity in a good rate, and the watertightness of the CO2 in the long term. The present work aims to study the parameters that have influence on the calculation of storage capacity, and for that purpose the appropriate technology has been developed for carrying out the research by mean of laboratory tests. Thus, a patent has been developed, "ATAP, equipo para ensayos petrofísicos (P201231913)", that has been used for developing the experimental part of this work. Once the technology has been developed, the study of different parameters, that have influence on the capacity of the storage, has been addressed developing different tests in ATAP. These tests define the storage volume which is related to the scope of different CO2 trap mechanisms, physical or chemical, in the storage. Tests that define the capacity of the storage to “accept” or “reject” the injected CO2, the injectivity, and tests led to determine possible leakages through injection wells. In this way we could talk about the watertightness in the aquifer in the long term and its influence on the storage capacity estimation. Together with the purpose of the storage capacity estimation, is the purpose of assuring the watertightness of such storages in the long term and anticipating the evolution of CO2 plume inside such aquifers. In order to fulfill this purpose, a dynamic model has been developed with ECLIPSE 300, for stablishing the methodology for the calculation of storage capacity estimation and the evolution of the CO2 plume, starting out with the tests carried out in ATAP. We present this work that establishes the methodology bases for the study of the influence of different petrophysics parameters in the calculation of the capacity of the storage together with the technological development of ATAP and its utilization for the determination of such parameters applicable to each aquifer.

Table 2 & Appendix A. Properties of operations injecting CO2 into saline aquifers

The experience from CO2 injection at pilot projects (Frio, Ketzin, Nagaoka, US Regional Partnerships) and existing commercial operations (Sleipner, Snøhvit, In Salah, acid-gas injection) demonstrates that CO2 geological storage in saline aquifers is technologically feasible. Monitoring and verification technologies have been tested and demonstrated to detect and track the CO2 plume in different subsurface geological environments. By the end of 2008, approximately 20 Mt of CO2 had been successfully injected into saline aquifers by existing operations. Currently, the highest injection rate and total storage volume for a single storage operation are approximately 1 Mt CO2/year and 25 Mt, respectively. If carbon capture and storage (CCS) is to be an effective option for decreasing greenhouse gas emissions, commercial-scale storage operations will require orders of magnitude larger storage capacity than accessed by the existing sites. As a result, new demonstration projects will need to develop and test injection strategies that consider multiple injection wells and the optimisation of the usage of storage space. To accelerate large-scale CCS deployment, demonstration projects should be selected that can be readily employed for commercial use; i.e. projects that fully integrate the capture, transport and storage processes at an industrial emissions source.

Optical assessment of impact and recovery of sedimentary pH profiles in ocean acidification and carbon capture and storage research

Available methods for measuring the impact of ocean acidification (OA) and leakage from carbon capture and storage (CCS) on marine sedimentary pH profiles are unsuitable for replicated experimental setups. To overcome this issue, a novel optical sensor application is presented, using off-the-shelf optode technology (MOPP). The application is validated using microprofiling, during a CCS leakage experiment, where the impact and recovery from a high CO2 plume was investigated in two types of natural marine sediment. MOPP offered user-friendliness, speed of data acquisition, robustness to sediment type, and large sediment depth range. This ensemble of characteristics overcomes many of the challenges found with other pH measuring methods, in OA and CCS research. The impact varied greatly between sediment types, depending on baseline pH variability and sediment permeability. Sedimentary pH profile recovery was quick, with profiles close to control conditions 24 h after the cessation of the leak. However, variability of pH within the finer sediment was still apparent 4 days into the recovery phase. Habitat characteristics need therefore to be considered, to truly disentangle high CO2 perturbation impacts on benthic systems. Impacts on natural communities depend not only on the pH gradient caused by perturbation, but also on other processes that outlive the perturbation, adding complexity to recovery.

Geochemistry and radiogenic isotopes of ODP Hole 157-953C lavas

The Canary Island primitive basaltic magmas are thought to be derived from an HIMU-type upwelling mantle containing isotopically depleted (NMORB)-type component having interacted with an enriched (EM)-type component, the origin of which is still a subject of debate. We studied the relationships between Ni, Mn and Ca concentrations in olivine phenocrysts (85.6-90.0 mol.% Fo, 1,722-3,915 ppm Ni, 1,085-1,552 ppm Mn, 1,222-3,002 ppm Ca) from the most primitive subaerial and ODP Leg 157 high-silica (picritic to olivine basaltic) lavas with their bulk rock Sr-Nd-Pb isotope compositions (87Sr/86Sr = 0.70315-0.70331, 143Nd/144Nd = 0.51288-0.51292, 206Pb/204Pb = 19.55-19.93, 207Pb/204Pb = 15.60-15.63, 208Pb/204Pb = 39.31-39.69). Our data point toward the presence of both a peridotitic and a pyroxenitic component in the magma source. Using the model (Sobolev et al., 2007, Science Vol 316) in which the reaction of Si-rich melts originated during partial melting of eclogite (a high pressure product of subducted oceanic crust) with ambient peridotitic mantle forms olivine-free reaction pyroxenite, we obtain an end member composition for peridotite with 87Sr/86Sr = 0.70337, 143Nd/144Nd = 0.51291, 206Pb/204Pb = 19.36, 207Pb/204Pb = 15.61 and 208Pb/204Pb = 39.07 (EM-type end member), and pyroxenite with 87Sr/86Sr = 0.70309, 143Nd/144Nd = 0.51289, 206Pb/204Pb = 20.03, 207Pb/204Pb = 15.62 and 208Pb/204Pb = 39.84 (HIMU-type end member). Mixing of melts from these end members in proportions ranging from 70% peridotite and 30% pyroxenite to 28% peridotite and 72% pyroxenite derived melt fractions can generate the compositions of the most primitive Gran Canaria shield stage lavas. Combining our results with those from the low-silica rocks from the western Canary Islands (Gurenko et al., 2009, doi:10.1016/j.epsl.2008.11.013), at least four distinct components are required. We propose that they are (1) HIMU-type pyroxenitic component (representing recycled ocean crust of intermediate age) from the plume center, (2) HIMU-type peridotitic component (ancient recycled ocean crust stirred into the ambient mantle) from the plume margin, (3) depleted, MORB-type pyroxenitic component (young recycled oceanic crust) in the upper mantle entrained by the plume, and (4) EM-type peridotitic component from the asthenosphere or lithosphere above the plume center.

The buoyant plume above a heat source

The boundary-layer type conservation equations of mass, momentum and energy for the steady free turbulent flow in gravitational convection over heat sources are set up for both two-dimensional and axisymmetric cases. These are reduced to ordinary differential equations in a similarity parameter by suitable transformations. The three classical hypotheses of turbulent diffusion-the Constant Exchange Coefficient hypothesis, Prandtl's Momentum Transfer theory and Taylor's Vorticity Transfer theory-are then incorporated into these equations in succession. The resulting equations are solved numerically and the results compared with some experimental results on gravitational convection over heat sources reported by Rouse et al.

Free atmospheric CO2 enrichment increased above ground biomass but did not affect symbiotic N2-fixation and soil carbon dynamics in a mixed deciduous stand in Wales

Through increases in net primary production (NPP), elevated CO2 is hypothesizes to increase the amount of plant litter entering the soil. The fate of this extra carbon on the forest floor or in mineral soil is currently not clear. Moreover, increased rates of NPP can be maintained only if forests can escape nitrogen limitation. In a Free atmospheric CO2 Enrichment (FACE) experiment near Bangor, Wales, 4 ambient CO2 and 4 FACE plots were planted with patches of Betula pendula, Alnus glutinosa and Fagus sylvatica on a former arable field. Four years after establishment, only a shallow L forest floor litter layer had formed due to intensive bioturbation. Total soil C and N contents increased irrespective of treatment and species as a result of afforestation. We could not detect an additional C sink in the soil, nor were soil C stabilization processes affected by FACE. We observed a decrease of leaf N content in Betula and Alnus under FACE, while the soil C/N ratio decreased regardless of CO2 treatment. The ratio of N taken up from the soil and by N2-fixation in Alnus was not affected by FACE. We infer that increased nitrogen use efficiency is the mechanism by which increased NPP is sustained under elevated CO2 at this site.

Coppicing shifts CO2 stimulation of poplar productivity to above-ground pools: a synthesis of leaf to stand level results from the POP/EUROFACE experiment

A poplar short rotation coppice (SRC) grown for the production of bioenergy can combine carbon (C) storage with fossil fuel substitution. Here, we summarize the responses of a poplar (Populus) plantation to 6 yr of free air CO2 enrichment (POP/EUROFACE consisting of two rotation cycles). We show that a poplar plantation growing in nonlimiting light, nutrient and water conditions will significantly increase its productivity in elevated CO2 concentrations ([CO2]). Increased biomass yield resulted from an early growth enhancement and photosynthesis did not acclimate to elevated [CO2]. Sufficient nutrient availability, increased nitrogen use efficiency (NUE) and the large sink capacity of poplars contributed to the sustained increase in C uptake over 6 yr. Additional C taken up in high [CO2] was mainly invested into woody biomass pools. Coppicing increased yield by 66% and partly shifted the extra C uptake in elevated [CO2] to above-ground pools, as fine root biomass declined and its [CO2] stimulation disappeared. Mineral soil C increased equally in ambient and elevated [CO2] during the 6 yr experiment. However, elevated [CO2] increased the stabilization of C in the mineral soil. Increased productivity of a poplar SRC in elevated [CO2] may allow shorter rotation cycles, enhancing the viability of SRC for biofuel production.

Fluxes of CO2 above a sugarcane plantation in Brazil

Fluxes of CO2 were measured above a sugarcane plantation using the eddy-covariance method covering two growth cycles, representing the second and third re-growth (ratoons) harvested with stubble burning. The total net ecosystem exchange (NEE) in the first cycle (second ratoon, 393 days long) was −1964 ± 44 g C m−2; the gross ecosystem productivity (GEP) was 3612 ± 46 g C m−2 and the ecosystem respiration (RE) was 1648 ± 14 g C m−2. The NEE and GEP totals in the second cycle (third ratoon, 374 days long) decreased 51% and 25%, respectively and RE increased 7%. Accounting for the carbon emitted during biomass burning and the removal of stalks at harvest, net ecosystem carbon balance (NECB) totals were 102 ± 130 g C m−2 and 403 ± 84 g C m−2 in each cycle respectively. Thus the sugarcane agrosystem was approximately carbon neutral in the second ratoon. Yield in stalks fresh weight (SFW) attained the regional average (8.3 kg SFW m−2). Although it was a carbon source to the atmosphere, observed productivity (6.2 kg SFW m−2) of the third ratoon was 19% lower than the regional average due to the lower water availability observed during the initial 120 days of re-growth. However, the overall water use efficiency (WUE) achieved in the first cycle (4.3 g C kg−1 H2O) decreased only 5% in the second cycle. © 2013 Elsevier B.V. All rights reserved

(Table 1) Turbidity, temperature, and salinity in the hydrothermal plume and background waters above the northern part of the Mohn's Ridge at Station AMK36-3573-5

In July 1995 geological and biological studies in the axial zone of the northern part of the Mohn's Ridge (72°20'N) were carried out during Cruise 36 of R/V Akademik Mstislav Keldysh. Slopes of the neovolcanic zone, as well as a caldera on its crest were investigated with use of deep-sea manned submersibles Mir, geological and biological samples were also collected. Use of the Rosette sounding complex provided recognition of several major hydrothermal plumes. Bottom sediments of the marginal depression are enriched in metals characteristic for hydrothermal metalliferous sediments. Thus, a new unknown hydrothermal field was found.

(Table 1) H2O and CO2 contents in submarine basaltic glasses from Ontong Java Plateau

The Ontong Java Plateau in the western Pacific is anomalous compared to other oceanic large igneous provinces in that it appears to have never formed a large subaerial plateau. Paleoeruption depths (at 122 Ma) estimated from dissolved H2O and CO2 in submarine basaltic glass pillow rims vary from ~1100 m below sea level (mbsl) on the central part of the plateau to 2200-3000 mbsl on the northeastern edge. Our results suggest maximum initial uplift for the plateau of 2500-3600 m above the surrounding seafloor and 1500+/-400 m of postemplacement subsidence since 122 Ma. Our estimates of uplift and subsidence for the plateau are significantly less than predictions from thermal models of oceanic lithosphere, and thus our results are inconsistent with formation of the plateau by a high-temperature mantle plume. Two controversial possibilities to explain the anomalous uplift and subsidence are that the plateau (1) formed as a result of a giant bolide impact, or (2) formed from a mantle plume but has a lower crust of dense garnet granulite and/or eclogite; neither of these possibilities is fully consistent with all available geological, geophysical, and geochemical data. The origin of the largest magmatic event on Earth in the past 200 m.y. thus remains an enigma.