Upscaled modeling of CO2 injection and migration with coupled thermal processes

A practical modeling approach for CO2 storage over relatively large length and time scales is the vertical-equilibrium model, which solves partially integrated conservation equations for flow in two lateral dimensions. We couple heat transfer within the vertical equilibrium framework for fluid flow, focusing on thermal processes that most impact the CO2 plume. We investigate a simplified representation of heat exchange that also includes transport of heat within the plume. In addition, we explore available CO2 thermodynamic models for reliable prediction of density under different injection pressures and temperatures. The model concept is demonstrated on simplified systems.publishedVersio

Operator splitting methods for systems of convection-diffusion equations: nonlinear error mechanisms and correction strategies

Authors final draft post-refereeing.Many numerical methods for systems of convection-diffusion equations are based upon an operator splitting formulation, where convective and diffusive forces are accounted for in separate substeps. We describe the nonlinear mechanism of the splitting error in such numerical methods in the onedimensional case, a mechanism that is intimately linked to the local linearizations introduced implicitly in the (hyperbolic) convection steps by the use of an entropy condition. For convection-dominated flows, we demonstrate that operator splitting methods typically generate a numerical widening of viscous fronts, unless the splitting step is of the same magnitude as the diffusion scale. To compensate for the potentially damaging splitting error, we propose a corrected operator splitting (COS) method for general systems of convection-diffusion equations with the ability of correctly resolving the nonlinear balance between the convective and diffusive forces. In particular, COS produces viscous shocks with correct structure also when the splitting step is large. A front tracking method for systems of conservation laws, which in turn relies heavily on a Riemann solver, constitutes an important part of our COS strategy. The proposed COS method is successfully applied to a system modelling two-phase, multicomponent flow in porous media and a triangular system modelling three-phase flow

An Ellam Scheme for Advection-Diffusion Equations in Two Dimensions

We develop an Eulerian{Lagrangian localized adjoint method (ELLAM) to solve two-dimensional advection-diusion equations with all combinations of inflow and outflow Dirichlet, Neumann, and flux boundary conditions. The ELLAM formalism provides a systematic framework for implementation of general boundary conditions, leading to mass-conservative numerical schemes. The computational advantages of the ELLAM approximation have been demonstrated for a number of one-dimensional transport systems; practical implementations of ELLAM schemes in multiple spatial dimensions that require careful algorithm development are discussed in detail in this paper. Extensive numerical results are presented to compare the ELLAM scheme with many widely used numerical methods and to demonstrate the strength of the ELLAM scheme

Field-case simulation of CO2-plume migration using vertical-equilibrium models

When injected in deep saline aquifers, CO2 moves radially away from the injection well and progressively higher in the formation because of buoyancy forces. Analyzes have shown that after the injection period, CO2 will potentially migrate over several kilometers in the horizontal direction but only tens of meters in the vertical direction, limited by the aquifer caprock. Because of the large horizontal plume dimensions, three-dimensional numerical simulations of the plume migration over long periods of time are computationally intensive. Thus, to get results within a reasonable time frame, one is typically forced to use coarse meshes and long time steps which result in inaccurate results because of numerical errors in resolving the plume tip. Given the large aspect ratio between the vertical and horizontal plume dimensions, it is reasonable to approximate the CO2 migration using vertically averaged models. Such models can, in many cases, be more accurate than coarse three-dimensional computations. In particular, models based on vertical equilibrium (VE) are attractive to simulate the long-term fate of CO2 sequestered into deep saline aquifers. The reduced spatial dimensionality resulting from the vertical integration ensures that the computational performance of VE models exceeds the performance of standard three-dimensional models. Thus, VE models are suitable to study the long-time and large-scale behavior of plumes in real large-scale CO2-injection projects. We investigate the use of VE models to simulate CO2 migration in a real large-scale field case based on data from the Sleipner site in the North Sea. We discuss the potential and limitations of VE models and show how VE models can be used to give reliable estimates of long-term CO2 migration. In particular, we focus on a VE formulation that incorporates the aquifer geometry and heterogeneity, and that considers the effects of hydrodynamic and residual trapping. We compare the results of VE simulations with standard reservoir simulation tools on test cases and discuss their advantages and limitations and show how, provided that certain conditions are met, they can be used to give reliable estimates of long-term CO2 migration.publishedVersio

The role of entry screening in case finding of tuberculosis among asylum seekers in Norway

<p>Abstract</p> <p>Background</p> <p>Most new cases of active tuberculosis in Norway are presently caused by imported strains and not transmission within the country. Screening for tuberculosis with a Mantoux test of everybody and a chest X-ray of those above 15 years of age is compulsory on arrival for asylum seekers.</p> <p>We aimed to assess the effectiveness of entry screening of a cohort of asylum seekers. Cases detected by screening were compared with cases detected later. Further we have characterized cases with active tuberculosis.</p> <p>Methods</p> <p>All asylum seekers who arrived at the National Reception Centre between January 2005 - June 2006 with an abnormal chest X-ray or a Mantoux test ≥ 6 mm were included in the study and followed through the health care system. They were matched with the National Tuberculosis Register by the end of May 2008.</p> <p>Cases reported within two months after arrival were defined as being detected by screening.</p> <p>Results</p> <p>Of 4643 eligible asylum seekers, 2237 were included in the study. Altogether 2077 persons had a Mantoux ≥ 6 mm and 314 had an abnormal chest X-ray. Of 28 cases with tuberculosis, 15 were detected by screening, and 13 at 4-27 months after arrival. Abnormal X-rays on arrival were more prevalent among those detected by screening. Female gender and Somalian origin increased the risk for active TB.</p> <p>Conclusion</p> <p>In spite of an imperfect follow-up of screening results, a reasonable number of TB cases was identified by the programme, with a predominance of pulmonary TB.</p

A pore network model for calculation of interfacial velocities

Two-phase flow in porous media is characterized by fluid-fluid interfaces that separate the fluid phases at the pore scale. These interfaces support pressure differences between phases, and their dynamics allow for saturation changes within the porous medium. Dynamic porescale network models allow interface dynamics to be modeled explicitly, such that each fluid-fluid interface within a network of pores is tracked explicitly in space and time. Because these models produce a detailed description of both phase and interface dynamics, results from these models can be averaged to provide values for many upscaled variables. These include nontraditional variables such as amounts of interfacial area, and volumeaveraged interfacial velocities. While these upscaled variables provide insights into the underlying dynamics of two-phase flow systems, they also allow new theories involving interfacial area to be tested directly. Results from a test of one such theory reveals that proposed constitutive equations involving interfacial velocities fail to match results from the network model for most cases considered. Therefore the macroscopic equations require additional development before they can be used for macro-scale simulations