Numerical analysis of thermal effects during carbon dioxide injection with enhanced gas recovery: a theoretical case study for the Altmark gas field

Prediction about reservoir temperature change during carbon dioxide injection requires consideration of all, often subtle, thermal effects. In particular, Joule?CThomson cooling (JTC) and the viscous heat dissipation (VHD) effect are factors that cause flowing fluid temperature to differ from the static formation temperature. In this work, warm-back behavior (thermal recovery after injection completed), as well as JTC and VHD effects, at a multi-layered depleted gas reservoir are demonstrated numerically. OpenGeoSys (OGS) is able to solve coupled partial differential equations for pressure, temperature and mole-fraction of each component of the mixture with a combination of monolithic and staggered approaches. The Galerkin finite element approach is adapted for space discretization of governing equations, whereas for temporal discretization, a generalized implicit single-step scheme is used. For numerical modeling of warm-back behavior, we chose a simplified test case of carbon dioxide injection. This test case is numerically solved by using OGS and FeFlow simulators independently. OGS differs from FeFlow in the capability of representing multi-componential effects on warm-back behavior. We verify both code results by showing the close comparison of shut-in temperature profiles along the injection well. As the JTC cooling rate is inversely proportional to the volumetric heat capacity of the solid matrix, the injection layers are cooled faster as compared to the non-injection layers. The shut-in temperature profiles are showing a significant change in reservoir temperature; hence it is important to account for thermal effects in injection monitoring.     

Europe’s longest-operating on-shore CO2 storage site at Ketzin, Germany: a progress report after three years of injection

The Ketzin pilot site, led by the GFZ German Research Centre for Geosciences, is Europe??s longest-operating on-shore CO₂ storage site with the aim of increasing the understanding of geological storage of CO₂ in saline aquifers. Located near Berlin, the Ketzin pilot site is an in situ laboratory for CO₂ storage in an anticlinal structure in the Northeast German Basin. Starting research within the framework of the EU project CO₂SINK in 2004, Ketzin is Germany??s first CO₂ storage site and fully in use since the injection began in June 2008. After 39?months of operation, about 53,000 tonnes of CO₂ have been stored in 630?C650?m deep sandstone units of the Upper Triassic Stuttgart Formation. An extensive monitoring program integrates geological, geophysical and geochemical investigations at Ketzin for a comprehensive characterization of the reservoir and the CO₂ migration at various scales. Integrating a unique field and laboratory data set, both static geological modeling and dynamic simulations are regularly updated. The Ketzin project successfully demonstrates CO₂ storage in a saline aquifer on a research scale. The results of monitoring and modeling can be summarized as follows: (1) Since the start of the CO₂ injection in June 2008, the operation has been running reliably and safely. (2) Downhole pressure data prove correlation between the injection rate and the reservoir pressure and indicates the presence of an overall dynamic equilibrium within the reservoir. (3) The extensive geochemical and geophysical monitoring program is capable of detecting CO₂ on different scales and gives no indication for any leakage. (4) Numerical simulations (history matching) are in good agreement with the monitoring results.     

Modeling, parameterization and evaluation of monitoring methods for CO2 storage in deep saline formations: the CO2-MoPa project

Capture and geological sequestration of CO₂ from large industrial sources is considered a measure for reducing anthropogenic emissions of CO₂ and thus mitigating climate change. One of the main storage options proposed are deep saline formations, as they provide the largest potential storage capacities among the geologic options. A thorough assessment of this type of storage site therefore is required. The CO₂-MoPa project aims at contributing to the dimensioning of CO₂ storage projects and to evaluating monitoring methods for CO₂ injection by an integrated approach. For this, virtual, but realistic test sites are designed geometrically and fully parameterized. Numerical process models are developed and then used to simulate the effects of a CO₂ injection into the virtual test sites. Because the parameterization of the virtual sites is known completely, investigation as well as monitoring methods can be closely examined and evaluated by comparing the virtual monitoring result with the simulation. To this end, the monitoring or investigation method is also simulated, and the (virtual) measurements are recorded and evaluated like real data. Application to a synthetic site typical for the north German basin showed that pressure response has to be evaluated taking into account the layered structure of the storage system. Microgravimetric measurements are found to be promising for detecting the CO₂ phase distribution. A combination of seismic and geoelectric measurements can be used to constrain the CO₂ phase distribution for the anticline system used in the synthetic site.     

CLEAN: project overview on CO2 large-scale enhanced gas recovery in the Altmark natural gas field (Germany)

The joint research project CLEAN was conducted in the years 2008?C2011 by a German research and development (R&D) alliance of 16 partners from science and industry. The project was set-up as pilot project to investigate the processes relevant to enhanced gas recovery (EGR) by the injection of CO₂ into a subfield of the almost depleted Altmark natural gas field. Despite the setback that permission for active injection was not issued by the mining authority during the period of the project, important results fostering the understanding of processes linked with EGR were achieved. Work carried out led to a comprehensive evaluation of the EGR potential of the Altmark field and the Altensalzwedel subfield in particular. The calculated safety margins emphasize that technical well integrity of the 12 examined boreholes is given for EGR without a need for any further intervention. The laboratory and field tests confirm that the Altensalzwedel subfield is suitable for the injection of 100,000?t of CO₂. Numerical simulations provide sound predictions for the efficiency and safety of the EGR technology based on the CO₂ injection. The development and testing of different monitoring techniques facilitate an improved surveying of CO₂ storage sites in general. The CLEAN results provide the technological, logistic and conceptual prerequisites for implementing a CO₂-based EGR project in the Altmark and provide a benchmark for similar projects in the world.     

Impact of shore-parallel breakwaters on the microtidal coast of Hammam-Lif (North-East of Tunisia)

The Hammam Lif coast (North-East of Tunisia) suffered from erosion mainly after the storm that damaged the zone in 1981. To protect it, eight shore-parallel breakwaters were built. The present study is interested in presenting the effect of the protection works on the coastline behavior by using aerial photographs geo-referenced by the digital photogrammetric method and empirical relationships. The successive breakwaters made a general accretion with formation of seven tombolos and salient. The response of the coast to these structures was related to the characteristics of these latter (length, distance to coast, and spacing). The effectiveness of such protection works is shown by the evaluation of the accretion areas formed behind them. The breakwaters which are the longest, less distant from the shoreline, and less spaced formed well-developed tombolos. Whereas, the structures which are short, very far from the coastline, and far from each other caused very small tombolos and even a salient.     

Using Two-Point Set Statistics to Estimate the Diameter Distribution in Boolean Models with Circular Grains

This paper deals with the counting and measuring of overlapping circular objects in binary images, a problem that arises in the mineral processing industry when estimating the distribution of bubble diameters in flotation cells in order to predict the geometallurgical performance of the flotation process. By viewing the images as realizations of a stationary planar Boolean model with circular grains and by using two-point set statistics, it is possible to jointly estimate the average number of objects per unit area and their diameter distribution. When applied to real and simulated images, this method is shown to give accurate estimates, to be robust to the presence of noise contaminating the images and of moderate drifts in the number of objects per unit area, and to speed up the processing time with respect to currently used methods. Combined, these properties serve to make the method suitable for an on-line monitoring of the flotation process.     

Feasibility of geoelectrical monitoring and multiphase modeling for process understanding of gaseous CO2 injection into a shallow aquifer

Potential pathways in the subsurface may allow upwardly migrating gaseous CO₂ from deep geological storage formations to be released into near surface aquifers. Consequently, the availability of adequate methods for monitoring potential CO₂ releases in both deep geological formations and the shallow subsurface is a prerequisite for the deployment of Carbon Capture and Storage technology. Geoelectrical surveys are carried out for monitoring a small-scale and temporally limited CO₂ injection experiment in a pristine shallow aquifer system. Additionally, the feasibility of multiphase modeling was tested in order to describe both complex non-linear multiphase flow processes and the electrical behavior of partially saturated heterogeneous porous media. The suitability of geoelectrical methods for monitoring injected CO₂ and geochemically altered groundwater was proven. At the test site, geoelectrical measurements reveal significant variations in electrical conductivity in the order of 15?C30?%. However, site-specific conditions (e.g., geological settings, groundwater composition) significantly influence variations in subsurface electrical conductivity and consequently, the feasibility of geoelectrical monitoring. The monitoring results provided initial information concerning gaseous CO₂ migration and accumulation processes. Geoelectrical monitoring, in combination with multiphase modeling, was identified as a useful tool for understanding gas phase migration and mass transfer processes that occur due to CO₂ intrusions in shallow aquifer systems.     

Combined L1/L2 Kalman filter-based tracking scheme for weak signal environments

A Kalman filter-based method combining the energy of both L1 C/A and L2C GPS signals in a combined tracking loop method to enhance performance under adverse conditions is developed. Standard tracking methods and the ionospheric effect on GPS signals are reviewed and compared to a new Kalman filter that simultaneously estimates delay, phase and total electron content by combining L1 C/A and L2C code and phase discriminator outputs. The new filter is tested and compared to standard methods for tracking L1 C/A and L2C using both simulated and real data. The new method is found to have improved sensitivity of 3 dB compared to standard L1 tracking and 4.5 dB compared to standard L2C tracking while at the same time providing an accurate estimate of the total electron content along the signal path.     

VLBI-derived troposphere parameters during CONT08

Time-series of zenith wet and total troposphere delays as well as north and east gradients are compared, and zenith total delays (ZTD) are combined on the level of parameter estimates. Input data sets are provided by ten Analysis Centers (ACs) of the International VLBI Service for Geodesy and Astrometry (IVS) for the CONT08 campaign (12?C26 August 2008). The inconsistent usage of meteorological data and models, such as mapping functions, causes systematics among the ACs, and differing parameterizations and constraints add noise to the troposphere parameter estimates. The empirical standard deviation of ZTD among the ACs with regard to an unweighted mean is 4.6?mm. The ratio of the analysis noise to the observation noise assessed by the operator/software impact (OSI) model is about 2.5. These and other effects have to be accounted for to improve the intra-technique combination of VLBI-derived troposphere parameters. While the largest systematics caused by inconsistent usage of meteorological data can be avoided and the application of different mapping functions can be considered by applying empirical corrections, the noise has to be modeled in the stochastic model of intra-technique combination. The application of different stochastic models shows no significant effects on the combined parameters but results in different mean formal errors: the mean formal errors of the combined ZTD are 2.3?mm (unweighted), 4.4?mm (diagonal), 8.6?mm [variance component (VC) estimation], and 8.6?mm (operator/software impact, OSI). On the one hand, the OSI model, i.e. the inclusion of off-diagonal elements in the cofactor-matrix, considers the reapplication of observations yielding a factor of about two for mean formal errors as compared to the diagonal approach. On the other hand, the combination based on VC estimation shows large differences among the VCs and exhibits a comparable scaling of formal errors. Thus, for the combination of troposphere parameters a combination of the two extensions of the stochastic model is recommended.     

The economics of climate change impacts and policy benefits at city scale: a conceptual framework

Cities are particularly vulnerable to climate change and climate extremes in part because they concentrate many activities, people and wealth in limited areas. As a result they represent an important scale for assessment and understanding of climate change impacts. This paper provides a conceptual and methodological framework for urban economic impact assessment of climate change. The focus of the paper is on model-based analysis of future scenarios, including a framing of uncertainty for these projections, as one valuable input into the decision-making process. The paper highlights the main assessment difficulties, methods and tools, and selected examples across these areas. A number of challenges are unique to climate change impact assessment and others are unique to the problem of working at local scales. The paper also identifies the need for additional research, including the need for more integrated and systemic approaches to address climate change as a part of the urban development challenge as well as the need to assess the economic impacts of climate change and response policy at local scale.